Neurospora from Natural Populations: A Global Study

The First Observation of the Filamentous Fungus Neurospora crassa Growing in the Roots of the Grass Brachypodium distachyon

The model organism Neurospora crassa has been cultivated in laboratories since the 1920s and its saprotrophic lifestyle has been established for decades. However, beyond their role as saprotrophs, fungi engage in intricate relationships with plants, showcasing diverse connections ranging from mutualistic to pathogenic. Although N. crassa has been extensively investigated under laboratory conditions, its ecological characteristics remain largely unknown. In contrast, Brachypodium distachyon, a sweet grass closely related to significant crops, demonstrates remarkable ecological flexibility and participates in a variety of fungal interactions, encompassing both mutualistic and harmful associations. Through a comprehensive microscopic analysis using electron, fluorescence, and confocal laser scanning microscopy, we discovered a novel endophytic interaction between N. crassa and B. distachyon roots, where fungal hyphae not only thrive in the apoplastic space and vascular bundle but also may colonize plant root cells. This new and so far hidden trait of one of the most important fungal model organisms greatly enhances our view of N. crassa, opening new perspectives concerning the fungus' ecological role. In addition, we present a new tool for studying plant-fungus interspecies communication, combining two well-established model systems, which improves our possibilities of experimental design on the molecular level.

Longevity of Fungal Mycelia and Nuclear Quality Checks: a New Hypothesis for the Role of Clamp Connections in Dikaryons

This paper addresses the stability of mycelial growth in fungi and differences between ascomycetes and basidiomycetes. Starting with general evolutionary theories of multicellularity and the role of sex, we then discuss individuality in fungi. Recent research has demonstrated the deleterious consequences of nucleus-level selection in fungal mycelia, favoring cheaters with a nucleus-level benefit during spore formation but a negative effect on mycelium-level fitness. Cheaters appear to generally be loss-of-fusion (LOF) mutants, with a higher propensity to form aerial hyphae developing into asexual spores. Since LOF mutants rely on heterokaryosis with wild-type nuclei, we argue that regular single-spore bottlenecks can efficiently select against such cheater mutants. We then zoom in on ecological differences between ascomycetes being typically fast-growing but short-lived with frequent asexual-spore bottlenecks and basidiomycetes being generally slow-growing but long-lived and usually without asexual-spore bottlenecks. We argue that these life history differences have coevolved with stricter nuclear quality checks in basidiomycetes. Specifically, we propose a new function for clamp connections, structures formed during the sexual stage in ascomycetes and basidiomycetes but during somatic growth only in basidiomycete dikaryons. During dikaryon cell division, the two haploid nuclei temporarily enter a monokaryotic phase, by alternatingly entering a retrograde-growing clamp cell, which subsequently fuses with the subapical cell to recover the dikaryotic cell. We hypothesize that clamp connections act as screening devices for nuclear quality, with both nuclei continuously testing each other for fusion ability, a test that LOF mutants will fail. By linking differences in longevity of the mycelial phase to ecology and stringency of nuclear quality checks, we propose that mycelia have a constant and low lifetime cheating risk, irrespective of their size and longevity.

Influence of scat ageing on the gut microbiome: how old is too old?

BACKGROUND

The study of the host-microbiome by the collection of non-invasive samples has the potential to become a powerful tool for conservation monitoring and surveillance of wildlife. However, multiple factors can bias the quality of data recovered from scats, particularly when field-collected samples are used given that the time of defecation is unknown. Previous studies using scats have shown that the impact of aerobic exposure on the microbial composition is species-specific, leading to different rates of change in microbial communities. However, the impact that this aging process has on the relationship between the bacterial and fungal composition has yet to be explored. In this study, we measured the effects of time post-defecation on bacterial and fungal compositions in a controlled experiment using scat samples from the endangered koala (Phascolarctos cinereus).

RESULTS

We found that the bacterial composition remained stable through the scat aging process, while the fungal composition did not. The absence of an increase in facultative anaerobes and the stable population of obligate anaerobic bacteria were likely due to our sampling from the inner portion of the scat. We report a cluster of fungal taxa that colonises scats after defecation which can dilute the genetic material from the autochthonous mycoflora and inhibit recovery.

CONCLUSION

We emphasize the need to preserve the integrity of scat samples collected in the wild and combat the effects of time and provide strategies for doing so.

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.

Watershed and fire severity are stronger determinants of soil chemistry and microbiomes than within-watershed woody encroachment in a tallgrass prairie system

Fire can impact terrestrial ecosystems by changing abiotic and biotic conditions. Short fire intervals maintain grasslands and communities adapted to frequent, low-severity fires. Shrub encroachment that follows longer fire intervals accumulates fuel and can increase fire severity. This patchily distributed biomass creates mosaics of burn severities in the landscape-pyrodiversity. Afforded by a scheduled burn of a watershed protected from fires for 27 years, we investigated effects of woody encroachment and burn severity on soil chemistry and soil-inhabiting bacteria and fungi. We compared soils before and after fire within the fire-protected, shrub-encroached watershed and soils in an adjacent, annually burned and non-encroached watershed. Organic matter and nutrients accumulated in the fire-protected watershed but responded less to woody encroachment within the encroached watershed. Bioavailable nitrogen and phosphorus and fungal and bacterial communities responded to high-severity burn regardless of encroachment. Low-severity fire effects on soil nutrients differed, increased bacterial but decreased fungal diversity and effects of woody encroachment within the encroached watershed were minimal. High-severity burns in the fire-protected watershed led to a novel soil system state distinct from non-encroached and encroached soil systems. We conclude that severe fires may open grassland restoration opportunities to manipulate soil chemistry and microbial communities in shrub-encroached habitats.

Soil microbial trait-based strategies drive metabolic efficiency along an altitude gradient

Trait-based approaches provide a candidate framework for linking soil microbial community to ecosystem processes, yet how the trade-offs in different microbial traits regulate the community-level metabolic efficiency remains unknown. Herein we assessed the roles of the microbial taxa with particular trait strategies in mediating soil microbial metabolic efficiency along an altitude gradient on the Tibetan Plateau. Results showed that soil microbial metabolic efficiency declined with increasing altitude, as indicated by the increasing metabolic quotient (microbial respiration per unit biomass, qCO2) and decreasing carbon use efficiency (CUE). Both qCO2 and CUE were predominantly predicted by microbial physiological and taxonomic attributes after considering key environmental factors including soil pH, substrate quantity and quality. Specifically, the reduced metabolic efficiency was associated with higher investment into nutrient (particularly for phosphorus) acquisitions via enzymes. Furthermore, we identified key microbial assemblies selected by harsh environments (low substrate quality and temperature) as important predictors of metabolic efficiency. These results suggest that particular microbial assemblies adapted to nutrient limited and cold habitats, but at the expense of lower metabolic efficient at higher altitude. Our findings provide a candidate mechanism underlying community-level metabolic efficiency, which has important implications for microbial-mediated processes such as carbon dynamics under global climate changes.

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.

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

The ascomycete filamentous fungus Neurospora crassa played a historic role in experimental biology and became a model system for genetic research. Stimulated by a systematic effort to collect wild strains initiated by Stanford geneticist David Perkins, the genus Neurospora has also become a basic model for the study of evolutionary processes, speciation, and population biology. In this chapter, we will first trace the history that brought Neurospora into the era of population genomics. We will then cover the major contributions of population genomic investigations using Neurospora to our understanding of microbial biogeography and speciation, and review recent work using population genomics and genome-wide association mapping that illustrates the unique potential of Neurospora as a model for identifying the genetic basis of (potentially adaptive) phenotypes in filamentous fungi. The advent of population genomics has contributed to firmly establish Neurospora as a complete model system and we hope our review will entice biologists to include Neurospora in their research.

Wild Isolates of Neurospora crassa Reveal Three Conidiophore Architectural Phenotypes

The vegetative life cycle in the model filamentous fungus, Neurospora crassa, relies on the development of conidiophores to produce new spores. Environmental, temporal, and genetic components of conidiophore development have been well characterized; however, little is known about their morphological variation. We explored conidiophore architectural variation in a natural population using a wild population collection of 21 strains from Louisiana, United States of America (USA). Our work reveals three novel architectural phenotypes, Wild Type, Bulky, and Wrap, and shows their maintenance throughout the duration of conidiophore development. Furthermore, we present a novel image-classifier using a convolutional neural network specifically developed to assign conidiophore architectural phenotypes in a high-throughput manner. To estimate an inheritance model for this discrete complex trait, crosses between strains of each phenotype were conducted, and conidiophores of subsequent progeny were characterized using the trained classifier. Our model suggests that conidiophore architecture is controlled by at least two genes and has a heritability of 0.23. Additionally, we quantified the number of conidia produced by each conidiophore type and their dispersion distance, suggesting that conidiophore architectural phenotype may impact N. crassa colonization capacity.

Quantitative genetics of temperature performance curves of Neurospora crassa

Earth's temperature is increasing due to anthropogenic CO 2 emissions; and organisms need either to adapt to higher temperatures, migrate into colder areas, or face extinction. Temperature affects nearly all aspects of an organism's physiology via its influence on metabolic rate and protein structure, therefore genetic adaptation to increased temperature may be much harder to achieve compared to other abiotic stresses. There is still much to be learned about the evolutionary potential for adaptation to higher temperatures, therefore we studied the quantitative genetics of growth rates in different temperatures that make up the thermal performance curve of the fungal model system Neurospora crassa. We studied the amount of genetic variation for thermal performance curves and examined possible genetic constraints by estimating the G-matrix. We observed a substantial amount of genetic variation for growth in different temperatures, and most genetic variation was for performance curve elevation. Contrary to common theoretical assumptions, we did not find strong evidence for genetic trade-offs for growth between hotter and colder temperatures. We also simulated short-term evolution of thermal performance curves of N. crassa, and suggest that they can have versatile responses to selection.

Low Amplitude Boom-and-Bust Cycles Define the Septoria Nodorum Blotch Interaction

INTRODUCTION

Septoria nodorum blotch (SNB) is a complex fungal disease of wheat caused by the Dothideomycete fungal pathogen Parastagonospora nodorum. The fungus infects through the use of necrotrophic effectors (NEs) that cause necrosis on hosts carrying matching dominant susceptibility genes. The Western Australia (WA) wheatbelt is a SNB "hot spot" and experiences significant under favorable conditions. Consequently, SNB has been a major target for breeders in WA for many years.

MATERIALS AND METHODS

In this study, we assembled a panel of 155 WA P. nodorum isolates collected over a 44-year period and compared them to 23 isolates from France and the USA using 28 SSR loci.

RESULTS

The WA P. nodorum population was clustered into five groups with contrasting properties. 80% of the studied isolates were assigned to two core groups found throughout the collection location and time. The other three non-core groups that encompassed transient and emergent populations were found in restricted locations and time. Changes in group genotypes occurred during periods that coincided with the mass adoption of a single or a small group of widely planted wheat cultivars. When introduced, these cultivars had high scores for SNB resistance. However, the field resistance of these new cultivars often declined over subsequent seasons prompting their replacement with new, more resistant varieties. Pathogenicity assays showed that newly emerged isolates non-core are more pathogenic than old isolates. It is likely that the non-core groups were repeatedly selected for increased virulence on the contemporary popular cultivars.

DISCUSSION

The low level of genetic diversity within the non-core groups, difference in virulence, low abundance, and restriction to limited locations suggest that these populations more vulnerable to a population crash when the cultivar was replaced by one that was genetically different and more resistant. We characterize the observed pattern as a low-amplitude boom-and-bust cycle in contrast with the classical high amplitude boom-and-bust cycles seen for biotrophic pathogens where the contrast between resistance and susceptibility is typically much greater. Implications of the results are discussed relating to breeding strategies for more sustainable SNB resistance and more generally for pathogens with NEs.

Preserving US microbe collections sparks future discoveries

Collections of micro‐organisms are a crucial element of life science research infrastructure but are vulnerable to loss and damage caused by natural or man‐made disasters, the untimely death or retirement of personnel, or the loss of research funding. Preservation of biological collections has risen in priority due to a new appreciation for discoveries linked to preserved specimens, emerging hurdles to international collecting and decreased funding for new collecting. While many historic collections have been lost, several have been preserved, some with dramatic rescue stories. Rescued microbes have been used for discoveries in areas of health, biotechnology and basic life science. Suggestions for long‐term planning for microbial stocks are listed, as well as inducements for long‐term preservation.

Identification and manipulation of Neurospora crassa genes involved in sensitivity to furfural

BACKGROUND

Biofuels derived from lignocellulosic biomass are a viable alternative to fossil fuels required for transportation. Following plant biomass pretreatment, the furan derivative furfural is present at concentrations which are inhibitory to yeasts. Detoxification of furfural is thus important for efficient fermentation. Here, we searched for new genetic attributes in the fungus Neurospora crassa that may be linked to furfural tolerance. The fact that furfural is involved in the natural process of sexual spore germination of N. crassa and that this fungus is highly amenable to genetic manipulations makes it a rational candidate for this study.

RESULTS

Both hypothesis-based and unbiased (random promotor mutagenesis) approaches were performed to identify N. crassa genes associated with the response to furfural. Changes in the transcriptional profile following exposure to furfural revealed that the affected processes were, overall, similar to those observed in Saccharomyces cerevisiae. N. crassa was more tolerant (by ~ 30%) to furfural when carboxymethyl cellulose was the main carbon source as opposed to sucrose, indicative of a link between carbohydrate metabolism and furfural tolerance. We also observed increased tolerance in a Δcre-1 mutant (CRE-1 is a key transcription factor that regulates the ability of fungi to utilize non-preferred carbon sources). In addition, analysis of aldehyde dehydrogenase mutants showed that ahd-2 (NCU00378) was involved in tolerance to furfural as well as the predicted membrane transporter NCU05580 (flr-1), a homolog of FLR1 in S. cerevisiae. Further to the rational screening, an unbiased approach revealed additional genes whose inactivation conferred increased tolerance to furfural: (i) NCU02488, which affected the abundance of the non-anchored cell wall protein NCW-1 (NCU05137), and (ii) the zinc finger protein NCU01407.

CONCLUSIONS

We identified attributes in N. crassa associated with tolerance or degradation of furfural, using complementary research approaches. The manipulation of the genes involved in furan sensitivity can provide a means for improving the production of biofuel producing strains. Similar research approaches can be utilized in N. crassa and other filamentous fungi to identify additional attributes relevant to other furans or toxic chemicals.

Mycotoxin concentrations in rice from three climatic locations in Africa as affected by grain quality, production site, and storage duration

Information on the mycotoxin contamination of rice in Africa is limited although the risk of contamination is high. In this study, domestic milled rice processed by actors using suboptimal methods was purchased and total fumonisin (FUM), zearalenone, and aflatoxin concentrations determined at 0, 90, and 180 days after storage. Three different climatic locations, Cotonou (Benin) in the Guinea savanna, Yaoundé (Cameroon) in the Tropical forest, and N'diaye (Senegal) in the Sahel, were selected as storage sites. Subsets of the samples collected from Glazoue (Benin), Ndop (Cameroon), and Dagana (Senegal) were stored in plastic woven bags under room conditions in the respective sites with or without calcium oxide (burnt scallop shell-BSS, 0.1% w/w) treatment. Multivariance analysis showed that FUM concentration was positively influenced by the duration of storage only while zearalenone concentration was negatively influenced by relative humidity and head rice but positively by impurities. Zearalenone concentration was also influenced by sample collection/storage location, processing type, and duration of storage. Aflatoxin concentration was influenced negatively by storage room temperature and head rice but positively by impurities and chalky grains. In addition, aflatoxin concentration was influenced by collection/storage location and processing type. BSS treatment followed by storage for 6 months had no effect on the concentration of the three assessed mycotoxins. Strategies to reduce the risk of mycotoxin contamination in study sites will include the improvement of physical rice quality through better pre- and postharvest practices and proper packaging of both treated rice and untreated rice in hermetic systems before marketing and storage.

Effects of acclimation time and epigenetic mechanisms on growth of Neurospora in fluctuating environments

Reaction norms or tolerance curves have often been used to predict how organisms deal with fluctuating environments. A potential drawback is that reaction norms measured in different constant environments may not capture all aspects of organismal responses to fluctuating environments. We examined growth of the filamentous fungus Neurospora crassa in fluctuating temperatures and tested if growth in fluctuating temperatures can be explained simply by the growth in different constant temperatures or if more complex models are needed. In addition, as previous studies on fluctuating environments have revealed that past temperatures that organisms have experienced can affect their response to current temperature, we tested the roles of different epigenetic mechanisms in response to fluctuating environments using different mutants. We found that growth of Neurospora can be predicted in fluctuating temperatures to some extent if acclimation times are taken into account in the model. Interestingly, while fluctuating environments have been linked with epigenetic responses, we found only some evidence of involvement of epigenetic mechanisms on tolerating fluctuating temperatures. Mutants which lacked H3K4 or H3K36 methylation had slightly impaired response to temperature fluctuations, in addition the H3K4 methylation mutant and a mutant in the RNA interference pathway had altered acclimation times.

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

The first eukaryotic genome to be sequenced was fungal, and there continue to be more sequenced genomes in the kingdom Fungi than in any other eukaryotic kingdom. Comparison of these genomes reveals many sources of genetic variation, from single nucleotide polymorphisms to horizontal gene transfer and on to changes in the arrangement and number of chromosomes, not to mention endofungal bacteria and viruses. Population genomics shows that all sources generate variation all the time and implicate natural selection as the force maintaining genome stability. Variation in wild populations is a rich resource for associating genetic variation with phenotypic variation, whether through quantitative trait locus mapping, genome-wide association studies, or reverse ecology. Subjects of studies associating genetic and phenotypic variation include model fungi, e.g., Saccharomyces and Neurospora, but pioneering studies have also been made with fungi pathogenic to plants, e.g., Pyricularia (= Magnaporthe), Zymoseptoria, and Fusarium, and to humans, e.g., Coccidioides, Cryptococcus, and Candida.

Metagenomic analysis of orange colored protrusions from the muscle of Queen Conch Lobatus gigas (Linnaeus, 1758)

The endangered marine gastropod, Lobatus gigas, is an important fishery resource in the Caribbean region. Microbiological and parasitological research of this species have been poorly addressed despite its role in ecological fitness, conservation status and prevention of potential pathogenic infections. This study identified taxonomic groups associated with orange colored protrusions in the muscle of queen conchs using histological analysis, 454 pyrosequencing, and a combination of PCR amplification and automated Sanger sequencing. The molecular approaches indicate that the etiological agent of the muscle protrusions is a parasite belonging to the subclass Digenea. Additionally, the scope of the molecular technique allowed the detection of bacterial and fungi clades in the assignment analysis. This is the first evidence of a digenean infection in the muscle of this valuable Caribbean resource.

Natural Variation of the Circadian Clock in Neurospora

Most living organisms on earth experience daily and expected changes from the rotation of the earth. For an organism, the ability to predict and prepare for incoming stresses or resources is a very important skill for survival. This cellular process of measuring daily time of the day is collectively called the circadian clock. Because of its fundamental role in survival in nature, there is a great interest in studying the natural variation of the circadian clock. However, characterizing the genetic and molecular mechanisms underlying natural variation of circadian clocks remains a challenging task. In this chapter, we will summarize the progress in studying natural variation of the circadian clock in the successful eukaryotic model Neurospora, which led to discovering many design principles of the molecular mechanisms of the eukaryotic circadian clock. Despite the success of the system in revealing the molecular mechanisms of the circadian clock, Neurospora has not been utilized to extensively study natural variation. We will review the challenges that hindered the natural variation studies in Neurospora, and how they were overcome. We will also review the advantages of Neurospora for natural variation studies. Since Neurospora is the model fungal species for circadian study, it represents over 5 million species of fungi on earth. These fungi play important roles in ecosystems on earth, and as such Neurospora could serve as an important model for understanding the ecological role of natural variation in fungal circadian clocks.

The ancestral levels of transcription and the evolution of sexual phenotypes in filamentous fungi

Changes in gene expression have been hypothesized to play an important role in the evolution of divergent morphologies. To test this hypothesis in a model system, we examined differences in fruiting body morphology of five filamentous fungi in the Sordariomycetes, culturing them in a common garden environment and profiling genome-wide gene expression at five developmental stages. We reconstructed ancestral gene expression phenotypes, identifying genes with the largest evolved increases in gene expression across development. Conducting knockouts and performing phenotypic analysis in two divergent species typically demonstrated altered fruiting body development in the species that had evolved increased expression. Our evolutionary approach to finding relevant genes proved far more efficient than other gene deletion studies targeting whole genomes or gene families. Combining gene expression measurements with knockout phenotypes facilitated the refinement of Bayesian networks of the genes underlying fruiting body development, regulation of which is one of the least understood processes of multicellular development.

A Review of Living Collections with Special Emphasis on Sustainability and Its Impact on Research Across Multiple Disciplines

Formal living collections have unique characteristics that distinguish them from other types of biorepositories. Comprising diverse resources, microbe culture collections, crop and biodiversity plant germplasm collections, and animal germplasm repositories are commonly allied with specific research communities or stakeholder groups. Among living collections, microbial culture collections have very long and unique life histories, with some being older than 100 years. Regulatory, financial, and technical developments have impacted living collections in many ways. International treaty obligations and restrictions on release of genetically modified organisms complicate the activities of living collections. Funding for living collections is a continuing challenge and threatens to create a two-tier system where medically relevant collections are well funded and all other collections are underfunded and hence understaffed. Molecular, genetic, and whole genome sequence analysis of contents of microbes and other living resource collections bring additional value to living collections.

Functional Evolution in Orthologous Cell-encoded RNA-dependent RNA Polymerases

Many eukaryotic organisms encode more than one RNA-dependent RNA polymerase (RdRP) that probably emerged as a result of gene duplication. Such RdRP paralogs often participate in distinct RNA silencing pathways and show characteristic repertoires of enzymatic activities in vitro However, to what extent members of individual paralogous groups can undergo functional changes during speciation remains an open question. We show that orthologs of QDE-1, an RdRP component of the quelling pathway in Neurospora crassa, have rapidly diverged in evolution at the amino acid sequence level. Analyses of purified QDE-1 polymerases from N. crassa (QDE-1(Ncr)) and related fungi, Thielavia terrestris (QDE-1(Tte)) and Myceliophthora thermophila (QDE-1(Mth)), show that all three enzymes can synthesize RNA, but the precise modes of their action differ considerably. Unlike their QDE-1(Ncr) counterpart favoring processive RNA synthesis, QDE-1(Tte) and QDE-1(Mth) produce predominantly short RNA copies via primer-independent initiation. Surprisingly, a 3.19 Å resolution crystal structure of QDE-1(Tte) reveals a quasisymmetric dimer similar to QDE-1(Ncr) Further electron microscopy analyses confirm that QDE-1(Tte) occurs as a dimer in solution and retains this status upon interaction with a template. We conclude that divergence of orthologous RdRPs can result in functional innovation while retaining overall protein fold and quaternary structure.

Genomic sequencing reveals historical, demographic and selective factors associated with the diversification of the fire-associated fungus Neurospora discreta

Delineating microbial populations, discovering ecologically relevant phenotypes and identifying migrants, hybrids or admixed individuals have long proved notoriously difficult, thereby limiting our understanding of the evolutionary forces at play during the diversification of microbial species. However, recent advances in sequencing and computational methods have enabled an unbiased approach whereby incipient species and the genetic correlates of speciation can be identified by examining patterns of genomic variation within and between lineages. We present here a population genomic study of a phylogenetic species in the Neurospora discreta species complex, based on the resequencing of full genomes (~37 Mb) for 52 fungal isolates from nine sites in three continents. Population structure analyses revealed two distinct lineages in South-East Asia, and three lineages in North America/Europe with a broad longitudinal and latitudinal range and limited admixture between lineages. Genome scans for selective sweeps and comparisons of the genomic landscapes of diversity and recombination provided no support for a role of selection at linked sites on genomic heterogeneity in levels of divergence between lineages. However, demographic inference indicated that the observed genomic heterogeneity in divergence was generated by varying rates of gene flow between lineages following a period of isolation. Many putative cases of exchange of genetic material between phylogenetically divergent fungal lineages have been discovered, and our work highlights the quantitative importance of genetic exchanges between more closely related taxa to the evolution of fungal genomes. Our study also supports the role of allopatric isolation as a driver of diversification in saprobic microbes.

Direct target network of the Neurospora crassa plant cell wall deconstruction regulators CLR-1, CLR-2, and XLR-1

Fungal deconstruction of the plant cell requires a complex orchestration of a wide array of intracellular and extracellular enzymes. In Neurospora crassa, CLR-1, CLR-2, and XLR-1 have been identified as key transcription factors regulating plant cell wall degradation in response to soluble sugars. The XLR-1 regulon was defined using a constitutively active mutant allele, resulting in hemicellulase gene expression and secretion under noninducing conditions. To define genes directly regulated by CLR-1, CLR-2, and XLR-1, we performed chromatin immunoprecipitation and next-generation sequencing (ChIPseq) on epitope-tagged constructs of these three transcription factors. When N. crassa is exposed to plant cell wall material, CLR-1, CLR-2, and XLR-1 individually bind to the promoters of the most strongly induced genes in their respective regulons. These include promoters of genes encoding cellulases for CLR-1 and CLR-2 (CLR-1/CLR-2) and promoters of genes encoding hemicellulases for XLR-1. CLR-1 bound to its regulon under noninducing conditions; however, this binding alone did not translate into gene expression and enzyme secretion. Motif analysis of the bound genes revealed conserved DNA binding motifs, with the CLR-2 motif matching that of its closest paralog in Saccharomyces cerevisiae, Gal4p. Coimmunoprecipitation studies showed that CLR-1 and CLR-2 act in a homocomplex but not as a CLR-1/CLR-2 heterocomplex.

Understanding fungal regulation of complex plant cell wall deconstruction pathways in response to multiple environmental signals via interconnected transcriptional circuits provides insight into fungus/plant interactions and eukaryotic nutrient sensing. Coordinated optimization of these regulatory networks is likely required for optimal microbial enzyme production.

Impact of clay mineral, wood sawdust or root organic matter on the bacterial and fungal community structures in two aged PAH-contaminated soils

The high organic pollutant concentration of aged polycyclic aromatic hydrocarbon (PAH)-contaminated wasteland soils is highly recalcitrant to biodegradation due to its very low bioavailability. In such soils, the microbial community is well adapted to the pollution, but the microbial activity is limited by nutrient availability. Management strategies could be applied to modify the soil microbial functioning as well as the PAH contamination through various amendment types. The impact of amendment with clay minerals (montmorillonite), wood sawdust and organic matter plant roots on microbial community structure was investigated on two aged PAH-contaminated soils both in laboratory and 1-year on-site pot experiments. Total PAH content (sum of 16 PAHs of the US-EPA list) and polar polycyclic aromatic compounds (pPAC) were monitored as well as the available PAH fraction using the Tenax method. The bacterial and fungal community structures were monitored using fingerprinting thermal gradient gel electrophoresis (TTGE) method. The abundance of bacteria (16S rRNA genes), fungi (18S rRNA genes) and PAH degraders (PAH-ring hydroxylating dioxygenase and catechol dioxygenase genes) was followed through qPCR assays. Although the treatments did not modify the total and available PAH content, the microbial community density, structure and the PAH degradation potential changed when fresh organic matter was provided as sawdust and under rhizosphere influence, while the clay mineral only increased the percentage of catechol-1,2-dioxygenase genes. The abundance of bacteria and fungi and the percentage of fungi relative to bacteria were enhanced in soil samples supplemented with wood sawdust and in the plant rhizospheric soils. Two distinct fungal populations developed in the two soils supplemented with sawdust, i.e. fungi related to Chaetomium and Neurospora genera and Brachyconidiellopsis and Pseudallescheria genera, in H and NM soils respectively. Wood sawdust amendment favoured the development of PAH-degrading bacteria holding Gram-negative PAH-ring hydroxylating dioxygenase, catechol-1,2-dioxygenase and catechol-2,3-dioxygenase genes. Regarding the total community structure, bacteria closely related to Thiobacillus (β-Proteobacteria) and Steroidobacter (γ-Proteobacteria) genera were favoured by wood sawdust amendment. In both soils, plant rhizospheres induced the development of fungi belonging to Ascomycota and related to Alternaria and Fusarium genera. Bacteria closely related to Luteolibacter (Verrucomicrobia) and Microbacterium (Actinobacteria) were favoured in alfalfa and ryegrass rhizosphere.

VIB1, a link between glucose signaling and carbon catabolite repression, is essential for plant cell wall degradation by Neurospora crassa

Filamentous fungi that thrive on plant biomass are the major producers of hydrolytic enzymes used to decompose lignocellulose for biofuel production. Although induction of cellulases is regulated at the transcriptional level, how filamentous fungi sense and signal carbon-limited conditions to coordinate cell metabolism and regulate cellulolytic enzyme production is not well characterized. By screening a transcription factor deletion set in the filamentous fungus Neurospora crassa for mutants unable to grow on cellulosic materials, we identified a role for the transcription factor, VIB1, as essential for cellulose utilization. VIB1 does not directly regulate hydrolytic enzyme gene expression or function in cellulosic inducer signaling/processing, but affects the expression level of an essential regulator of hydrolytic enzyme genes, CLR2. Transcriptional profiling of a Δvib-1 mutant suggests that it has an improper expression of genes functioning in metabolism and energy and a deregulation of carbon catabolite repression (CCR). By characterizing new genes, we demonstrate that the transcription factor, COL26, is critical for intracellular glucose sensing/metabolism and plays a role in CCR by negatively regulating cre-1 expression. Deletion of the major player in CCR, cre-1, or a deletion of col-26, did not rescue the growth of Δvib-1 on cellulose. However, the synergistic effect of the Δcre-1; Δcol-26 mutations circumvented the requirement of VIB1 for cellulase gene expression, enzyme secretion and cellulose deconstruction. Our findings support a function of VIB1 in repressing both glucose signaling and CCR under carbon-limited conditions, thus enabling a proper cellular response for plant biomass deconstruction and utilization.

Secret lifestyles of Neurospora crassa

Neurospora crassa has a long history as an excellent model for genetic, cellular, and biochemical research. Although this fungus is known as a saprotroph, it normally appears on burned vegetations or trees after forest fires. However, due to a lack of experimental evidence, the nature of its association with living plants remains enigmatic. Here we report that Scots pine (Pinus sylvestris) is a host plant for N. crassa. The endophytic lifestyle of N. crassa was found in its interaction with Scots pine. Moreover, the fungus can switch to a pathogenic state when its balanced interaction with the host is disrupted. Our data reveal previously unknown lifestyles of N. crassa, which are likely controlled by both environmental and host factors. Switching among the endophytic, pathogenic, and saprotrophic lifestyles confers upon fungi phenotypic plasticity in adapting to changing environments and drives the evolution of fungi and associated plants.

A global multilocus analysis of the model fungus Neurospora reveals a single recent origin of a novel genetic system

The large diversity of mating systems observed in the fungal kingdom underlines the importance of mating system change in fungal evolution. The selfing species Neurospora tetrasperma has evolved a novel method of achieving self-fertility by a mating system referred to as pseudohomothallism. However, little is known about the origin of N. tetrasperma and its relationship to the self-sterile, heterothallic, Neurospora species. In this study, we used a combination of phylogenetic and population genetic analyses to reconstruct the evolutionary history of N. tetrasperma and its heterothallic relatives. We sequenced 9 unlinked nuclear loci from 106 strains of N. tetrasperma sampled from across the globe, and a sample of 28 heterothallic strains of Neurospora. Our analyses provide strong support for monophyly of N. tetrasperma, but reject the monophyly of N. crassa. We estimate that N. tetrasperma is of a recent origin and that it diverged from the heterothallic species ∼1 million years ago. We also extend previous findings on the diversification within the N. tetrasperma clade, with 10 lineages identified. Taken together, these findings indicate that N. tetrasperma is younger than has been previously reported and that a rapid diversification of lineages has occurred within the N. tetrasperma clade.

Discovering functions of unannotated genes from a transcriptome survey of wild fungal isolates

Most fungal genomes are poorly annotated, and many fungal traits of industrial and biomedical relevance are not well suited to classical genetic screens. Assigning genes to phenotypes on a genomic scale thus remains an urgent need in the field. We developed an approach to infer gene function from expression profiles of wild fungal isolates, and we applied our strategy to the filamentous fungus Neurospora crassa. Using transcriptome measurements in 70 strains from two well-defined clades of this microbe, we first identified 2,247 cases in which the expression of an unannotated gene rose and fell across N. crassa strains in parallel with the expression of well-characterized genes. We then used image analysis of hyphal morphologies, quantitative growth assays, and expression profiling to test the functions of four genes predicted from our population analyses. The results revealed two factors that influenced regulation of metabolism of nonpreferred carbon and nitrogen sources, a gene that governed hyphal architecture, and a gene that mediated amino acid starvation resistance. These findings validate the power of our population-transcriptomic approach for inference of novel gene function, and we suggest that this strategy will be of broad utility for genome-scale annotation in many fungal systems. IMPORTANCE Some fungal species cause deadly infections in humans or crop plants, and other fungi are workhorses of industrial chemistry, including the production of biofuels. Advances in medical and industrial mycology require an understanding of the genes that control fungal traits. We developed a method to infer functions of uncharacterized genes by observing correlated expression of their mRNAs with those of known genes across wild fungal isolates. We applied this strategy to a filamentous fungus and predicted functions for thousands of unknown genes. In four cases, we experimentally validated the predictions from our method, discovering novel genes involved in the metabolism of nutrient sources relevant for biofuel production, as well as colony morphology and starvation resistance. Our strategy is straightforward, inexpensive, and applicable for predicting gene function in many fungal species.

Evidence for transceptor function of cellodextrin transporters in Neurospora crassa

Neurospora crassa colonizes burnt grasslands and metabolizes both cellulose and hemicellulose from plant cell walls. When switched from a favored carbon source to cellulose, N. crassa dramatically up-regulates expression and secretion of genes encoding lignocellulolytic enzymes. However, the means by which N. crassa and other filamentous fungi sense the presence of cellulose in the environment remains unclear. Previously, we have shown that a N. crassa mutant carrying deletions of three β-glucosidase enzymes (Δ3βG) lacks β-glucosidase activity, but efficiently induces cellulase gene expression and cellulolytic activity in the presence of cellobiose as the sole carbon source. These observations indicate that cellobiose, or a modified version of cellobiose, functions as an inducer of lignocellulolytic gene expression and activity in N. crassa. Here, we show that in N. crassa, two cellodextrin transporters, CDT-1 and CDT-2, contribute to cellulose sensing. A N. crassa mutant carrying deletions for both transporters is unable to induce cellulase gene expression in response to crystalline cellulose. Furthermore, a mutant lacking genes encoding both the β-glucosidase enzymes and cellodextrin transporters (Δ3βGΔ2T) does not induce cellulase gene expression in response to cellobiose. Point mutations that severely reduce cellobiose transport by either CDT-1 or CDT-2 when expressed individually do not greatly impact cellobiose induction of cellulase gene expression. These data suggest that the N. crassa cellodextrin transporters act as "transceptors" with dual functions - cellodextrin transport and receptor signaling that results in downstream activation of cellulolytic gene expression. Similar mechanisms of transceptor activity likely occur in related ascomycetes used for industrial cellulase production.

A comparative systems analysis of polysaccharide-elicited responses in Neurospora crassa reveals carbon source-specific cellular adaptations

Filamentous fungi are powerful producers of hydrolytic enzymes for the deconstruction of plant cell wall polysaccharides. However, the central question of how these sugars are perceived in the context of the complex cell wall matrix remains largely elusive. To address this question in a systematic fashion we performed an extensive comparative systems analysis of how the model filamentous fungus Neurospora crassa responds to the three main cell wall polysaccharides: pectin, hemicellulose and cellulose. We found the pectic response to be largely independent of the cellulolytic one with some overlap to hemicellulose, and in its extent surprisingly high, suggesting advantages for the fungus beyond being a mere carbon source. Our approach furthermore allowed us to identify carbon source-specific adaptations, such as the induction of the unfolded protein response on cellulose, and a commonly induced set of 29 genes likely involved in carbon scouting. Moreover, by hierarchical clustering we generated a coexpression matrix useful for the discovery of new components involved in polysaccharide utilization. This is exemplified by the identification of lat-1, which we demonstrate to encode for the physiologically relevant arabinose transporter in Neurospora. The analyses presented here are an important step towards understanding fungal degradation processes of complex biomass.

Large-scale introgression shapes the evolution of the mating-type chromosomes of the filamentous ascomycete Neurospora tetrasperma

The significance of introgression as an evolutionary force shaping natural populations is well established, especially in animal and plant systems. However, the abundance and size of introgression tracts, and to what degree interspecific gene flow is the result of adaptive processes, are largely unknown. In this study, we present medium coverage genomic data from species of the filamentous ascomycete Neurospora, and we use comparative genomics to investigate the introgression landscape at the genomic level in this model genus. We revealed one large introgression tract in each of the three investigated phylogenetic lineages of Neurospora tetrasperma (sizes of 5.6 Mbp, 5.2 Mbp, and 4.1 Mbp, respectively). The tract is located on the chromosome containing the locus conferring sexual identity, the mating-type (mat) chromosome. The region of introgression is confined to the region of suppressed recombination and is found on one of the two mat chromosomes (mat a). We used Bayesian concordance analyses to exclude incomplete lineage sorting as the cause for the observed pattern, and multilocus genealogies from additional species of Neurospora show that the introgression likely originates from two closely related, freely recombining, heterothallic species (N. hispaniola and N. crassa/N. perkinsii). Finally, we investigated patterns of molecular evolution of the mat chromosome in Neurospora, and we show that introgression is correlated with reduced level of molecular degeneration, consistent with a shorter time of recombination suppression. The chromosome specific (mat) and allele specific (mat a) introgression reported herein comprise the largest introgression tracts reported to date from natural populations. Furthermore, our data contradicts theoretical predictions that introgression should be less likely on sex-determining chromosomes. Taken together, the data presented herein advance our general understanding of introgression as a force shaping eukaryotic genomes.

Introgression is a process by which genetic material from one species becomes infiltrated into another, genetically distinct species. Introgression usually occurs via sexual reproduction: individuals of two species mate and produce a hybrid offspring, then the offspring repeatedly backcross with one of the parental species. Introgression has long been recognized as a key process in evolution, as it may contribute to speciation, diversification, and adaptation to new environments. The importance and prevalence of introgression has been well established in plant and animal systems, and in this study we use a fungal model system, Neurospora, to study the introgression at the genomic level. We gathered genomic data from six genomes, and by comparative genomics we revealed genetic transfer of DNA regions of unprecedentedly large sizes, covering over 50% of the mating-type chromosomes, and used phylogenetic analyses to reveal the origin and direction of the transfer. Introgression was found solely on the mating-type chromosomes, which contradicts theoretical predictions for sex-determining chromosomes. We argue that this unexpected pattern is due to the fact that fungi do not have differentiated sexes (female/male) and thereby are free from sex-biased evolutionary forces. Instead, we suggest that introgression between fungal species may result in reinvigoration of genomic regions exposed to suppressed recombination.

Quantifying functional heterothallism in the pseudohomothallic ascomycete Neurospora tetrasperma

Neurospora tetrasperma is a pseudohomothallic filamentous ascomycete that has evolved from heterothallic ancestors. Throughout its life cycle, it is predominantly heterokaryotic for mating type, and thereby self-fertile. However, studies of N. tetrasperma have revealed the occasional production of self-sterile asexual and sexual spores of a single-mating type, indicating that it can be functionally heterothallic. Here, we report the extensive sampling and isolation of natural, heterokaryotic, strains of N. tetrasperma from the United Kingdom (UK): 99 strains were collected from Surrey, England, and four from Edinburgh, Scotland. We verified by phylogenetic analyses that these strains belong to N. tetrasperma. We isolated cultures from single germinated asexual spores (conidia) from 17 of these newly sampled UK strains from Surrey, and 16 previously sampled strains of N. tetrasperma from New Zealand (NZ). Our results show that the N. tetrasperma strains from the UK population produced a significantly greater proportion of self-sterile, homokaryotic conidia than the NZ population: the proportion of homokaryotic conidia was 42.6 % (133/312 spores) and 15.3 % (59/386) from the UK and the NZ populations, respectively. Although homokaryons recovered from several strains show a bias for one of the mating types, the total ratio of mat A to mat a mating type in homokaryons (UK: 72/61, NZ 28/31) did not deviate significantly from the expected 1:1 ratio for either of these populations. These results indicate that different populations exhibit differences in their life cycle characteristics, and that a higher degree of outcrossing might be expected from the UK population. This study points to the importance of studying multiple strains and populations when investigating life history traits of an organism with a complex life cycle, as previously undetected differences between populations may be revealed.

Conserved and essential transcription factors for cellulase gene expression in ascomycete fungi

Rational engineering of filamentous fungi for improved cellulase production is hampered by our incomplete knowledge of transcriptional regulatory networks. We therefore used the model filamentous fungus Neurospora crassa to search for uncharacterized transcription factors associated with cellulose deconstruction. A screen of a N. crassa transcription factor deletion collection identified two uncharacterized zinc binuclear cluster transcription factors (clr-1 and clr-2) that were required for growth and enzymatic activity on cellulose, but were not required for growth or hemicellulase activity on xylan. Transcriptional profiling with next-generation sequencing methods refined our understanding of the N. crassa transcriptional response to cellulose and demonstrated that clr-1 and clr-2 were required for the bulk of that response, including induction of all major cellulase and some major hemicellulase genes. Functional CLR-1 was necessary for expression of clr-2 and efficient cellobiose utilization. Phylogenetic analyses showed that CLR-1 and CLR-2 are conserved in the genomes of most filamentous ascomycete fungi capable of degrading cellulose. In Aspergillus nidulans, a strain carrying a deletion of the clr-2 homolog (clrB) failed to induce cellulase gene expression and lacked cellulolytic activity on Avicel. Further manipulation of this control system in industrial production strains may significantly improve yields of cellulases for cellulosic biofuel production.

Unravelling the molecular basis for light modulated cellulase gene expression - the role of photoreceptors in Neurospora crassa

BACKGROUND

Light represents an important environmental cue, which exerts considerable influence on the metabolism of fungi. Studies with the biotechnological fungal workhorse Trichoderma reesei (Hypocrea jecorina) have revealed an interconnection between transcriptional regulation of cellulolytic enzymes and the light response. Neurospora crassa has been used as a model organism to study light and circadian rhythm biology. We therefore investigated whether light also regulates transcriptional regulation of cellulolytic enzymes in N. crassa.

RESULTS

We show that the N. crassa photoreceptor genes wc-1, wc-2 and vvd are involved in regulation of cellulase gene expression, indicating that this phenomenon is conserved among filamentous fungi. The negative effect of VVD on production of cellulolytic enzymes is thereby accomplished by its role in photoadaptation and hence its function in White collar complex (WCC) formation. In contrast, the induction of vvd expression by the WCC does not seem to be crucial in this process. Additionally, we found that WC-1 and WC-2 not only act as a complex, but also have individual functions upon growth on cellulose.

CONCLUSIONS

Genome wide transcriptome analysis of photoreceptor mutants and evaluation of results by analysis of mutant strains identified several candidate genes likely to play a role in light modulated cellulase gene expression. Genes with functions in amino acid metabolism, glycogen metabolism, energy supply and protein folding are enriched among genes with decreased expression levels in the wc-1 and wc-2 mutants. The ability to properly respond to amino acid starvation, i. e. up-regulation of the cross pathway control protein cpc-1, was found to be beneficial for cellulase gene expression. Our results further suggest a contribution of oxidative depolymerization of cellulose to plant cell wall degradation in N. crassa.

A relationship between carotenoid accumulation and the distribution of species of the fungus Neurospora in Spain

The ascomycete fungus Neurospora is present in many parts of the world, in particular in tropical and subtropical areas, where it is found growing on recently burned vegetation. We have sampled the Neurospora population across Spain. The sampling sites were located in the region of Galicia (northwestern corner of the Iberian peninsula), the province of Cáceres, the city of Seville, and the two major islands of the Canary Islands archipelago (Tenerife and Gran Canaria, west coast of Africa). The sites covered a latitude interval between 27.88° and 42.74°. We have identified wild-type strains of N. discreta, N. tetrasperma, N. crassa, and N. sitophila and the frequency of each species varied from site to site. It has been shown that after exposure to light Neurospora accumulates the orange carotenoid neurosporaxanthin, presumably for protection from UV radiation. We have found that each Neurospora species accumulates a different amount of carotenoids after exposure to light, but these differences did not correlate with the expression of the carotenogenic genes al-1 or al-2. The accumulation of carotenoids in Neurospora shows a correlation with latitude, as Neurospora strains isolated from lower latitudes accumulate more carotenoids than strains isolated from higher latitudes. Since regions of low latitude receive high UV irradiation we propose that the increased carotenoid accumulation may protect Neurospora from high UV exposure. In support of this hypothesis, we have found that N. crassa, the species that accumulates more carotenoids, is more resistant to UV radiation than N. discreta or N. tetrasperma. The photoprotection provided by carotenoids and the capability to accumulate different amounts of carotenoids may be responsible, at least in part, for the distribution of Neurospora species that we have observed across a range of latitudes.

Deciphering transcriptional regulatory mechanisms associated with hemicellulose degradation in Neurospora crassa

Hemicellulose, the second most abundant plant biomass fraction after cellulose, is widely viewed as a potential substrate for the production of liquid fuels and other value-added materials. Degradation of hemicellulose by filamentous fungi requires production of many different enzymes, which are induced by biopolymers or its derivatives and regulated mainly at the transcriptional level through transcription factors (TFs). Neurospora crassa, a model filamentous fungus, expresses and secretes enzymes required for plant cell wall deconstruction. To better understand genes specifically associated with degradation of hemicellulose, we applied secretome and transcriptome analysis to N. crassa grown on beechwood xylan. We identified 34 secreted proteins and 353 genes with elevated transcription on xylan. The xylanolytic phenotype of strains with deletions in genes identified from the secretome and transcriptome analysis of the wild type was assessed, revealing functions for known and unknown proteins associated with hemicellulose degradation. By evaluating phenotypes of strains containing deletions of predicted TF genes in N. crassa, we identified a TF (XLR-1; xylan degradation regulator 1) essential for hemicellulose degradation that is an ortholog to XlnR/XYR1 in Aspergillus and Trichoderma species, respectively, a major transcriptional regulator of genes encoding both cellulases and hemicellulases. Deletion of xlr-1 in N. crassa abolished growth on xylan and xylose, but growth on cellulose and cellulolytic activity were only slightly affected. To determine the regulatory mechanisms for hemicellulose degradation, we explored the transcriptional regulon of XLR-1 under xylose, xylanolytic, and cellulolytic conditions. XLR-1 regulated only some predicted hemicellulase genes in N. crassa and was required for a full induction of several cellulase genes. Hemicellulase gene expression was induced by a combination of release from carbon catabolite repression (CCR) and induction. This systematic analysis illustrates the similarities and differences in regulation of hemicellulose degradation among filamentous fungi.

Molecular Diversity of Fungi from Marine Oxygen-Deficient Environments (ODEs)

Molecular diversity surveys of marine fungi have demonstrated that the species richness known to date is just the tip of the iceberg and that there is a large extent of unknown fungal diversity in marine habitats. Reports of novel fungal lineages at higher taxonomic levels are documented from a large number of marine habitats, including the various marine oxygen-deficient environments (ODEs). In the past few years, a strong focus of eukaryote diversity research has been on a variety of ODEs, as these environments are considered to harbor a large number of organisms, which are highly divergent to known diversity and could provide insights into the early eukaryotic evolution. ODEs that have been targeted so far include shallow water sediments, hydrothermal vent systems, deep-sea basins, intertidal habitats, and fjords. Most, if not all, molecular diversity studies in marine ODEs have shown, that contrary to previous assumptions, fungi contribute significantly to the micro-eukaryotic community in such habitats. In this chapter, we have reanalyzed the environmental fungal sequences obtained from the molecular diversity survey in 14 different sites to obtain a comprehensive picture of fungal diversity in these marine habitats. The phylogenetic analysis of the fungal environmental sequences from various ODEs have grouped these sequences into seven distinct clades (Clade 1-7) clustering with well-known fungal taxa. Apart from this, four environmental clades (EnvClade A, B, C, and D) with exclusive environmental sequences were also identified. This has provided information on the positioning of the environmental sequences at different taxonomic levels within the major fungal phylums. The taxonomic distribution of these environmental fungal sequences into clusters and clades has also shown that they are not restricted by geographical boundaries. The distribution pattern together with the reports on the respiratory abilities of fungi under reduced oxygen conditions shows that they are highly adaptive and may have a huge ecological role in these oxygen deficient habitats.