Massive changes in genome architecture accompany the transition to self-fertility in the filamentous fungus Neurospora tetrasperma

Stepwise recombination suppression around the mating-type locus in an ascomycete fungus with self-fertile spores

Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.

Nuclear-specific gene expression in heterokaryons of the filamentous ascomycete Neurospora tetrasperma

Heterokaryosis is a system in which genetically distinct nuclei coexist within the same cytoplasm. While heterokaryosis dominates the life cycle of many fungal species, the transcriptomic changes associated with the transition from homokaryosis to heterokaryosis is not well understood. Here, we analyse gene expression profiles of homokaryons and heterokaryons from three phylogenetically and reproductively isolated lineages of the filamentous ascomycete Neurospora tetrasperma. We show that heterokaryons are transcriptionally distinct from homokaryons in the sexual stage of development, but not in the vegetative stage, suggesting that the phenotypic switch to fertility in heterokaryons is associated with major changes in gene expression. Heterokaryon expression is predominantly defined by additive effects of its two nuclear components. Furthermore, allele-specific expression analysis of heterokaryons with varying nuclear ratios show patterns of expression ratios strongly dependent on nuclear ratios in the vegetative stage. By contrast, in the sexual stage, strong deviations of expression ratios indicate a co-regulation of nuclear gene expression in all three lineages. Taken together, our results show two levels of expression control: additive effects suggest a nuclear level of expression, whereas co-regulation of gene expression indicate a heterokaryon level of control.

FunOrder: A robust and semi-automated method for the identification of essential biosynthetic genes through computational molecular co-evolution

Secondary metabolites (SMs) are a vast group of compounds with different structures and properties that have been utilized as drugs, food additives, dyes, and as monomers for novel plastics. In many cases, the biosynthesis of SMs is catalysed by enzymes whose corresponding genes are co-localized in the genome in biosynthetic gene clusters (BGCs). Notably, BGCs may contain so-called gap genes, that are not involved in the biosynthesis of the SM. Current genome mining tools can identify BGCs, but they have problems with distinguishing essential genes from gap genes. This can and must be done by expensive, laborious, and time-consuming comparative genomic approaches or transcriptome analyses. In this study, we developed a method that allows semi-automated identification of essential genes in a BGC based on co-evolution analysis. To this end, the protein sequences of a BGC are blasted against a suitable proteome database. For each protein, a phylogenetic tree is created. The trees are compared by treeKO to detect co-evolution. The results of this comparison are visualized in different output formats, which are compared visually. Our results suggest that co-evolution is commonly occurring within BGCs, albeit not all, and that especially those genes that encode for enzymes of the biosynthetic pathway are co-evolutionary linked and can be identified with FunOrder. In light of the growing number of genomic data available, this will contribute to the studies of BGCs in native hosts and facilitate heterologous expression in other organisms with the aim of the discovery of novel SMs.

The discovery and description of novel fungal secondary metabolites promises novel antibiotics, pharmaceuticals, and other useful compounds. A way to identify novel secondary metabolites is to express the corresponding genes in a suitable expression host. Consequently, a detailed knowledge or an accurate prediction of these genes is necessary. In fungi, the genes are co-localized in so-called biosynthetic gene clusters. Notably, the clusters may also contain genes that are not necessary for the biosynthesis of the secondary metabolites, so-called gap genes. We developed a method to detect co-evolved genes within the clusters and demonstrated that essential genes are co-evolving and can thus be differentiated from the gap genes. This adds an additional layer of information, which can support researchers with their decisions on which genes to study and express for the discovery of novel secondary metabolites.

Genomic diversity, chromosomal rearrangements, and interspecies hybridization in the Ogataea polymorpha species complex

The methylotrophic yeast Ogataea polymorpha has long been a useful system for recombinant protein production, as well as a model system for methanol metabolism, peroxisome biogenesis, thermotolerance, and nitrate assimilation. It has more recently become an important model for the evolution of mating-type switching. Here, we present a population genomics analysis of 47 isolates within the O. polymorpha species complex, including representatives of the species O. polymorpha, Ogataea parapolymorpha, Ogataea haglerorum, and Ogataea angusta. We found low levels of nucleotide sequence diversity within the O. polymorpha species complex and identified chromosomal rearrangements both within and between species. In addition, we found that one isolate is an interspecies hybrid between O. polymorpha and O. parapolymorpha and present evidence for loss of heterozygosity following hybridization.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Recombination suppression and evolutionary strata around mating-type loci in fungi: documenting patterns and understanding evolutionary and mechanistic causes

Genomic regions determining sexual compatibility often display recombination suppression, as occurs in sex chromosomes, plant self-incompatibility loci and fungal mating-type loci. Regions lacking recombination can extend beyond the genes determining sexes or mating types, by several successive steps of recombination suppression. Here we review the evidence for recombination suppression around mating-type loci in fungi, sometimes encompassing vast regions of the mating-type chromosomes. The suppression of recombination at mating-type loci in fungi has long been recognized and maintains the multiallelic combinations required for correct compatibility determination. We review more recent evidence for expansions of recombination suppression beyond mating-type genes in fungi ('evolutionary strata'), which have been little studied and may be more pervasive than commonly thought. We discuss testable hypotheses for the ultimate (evolutionary) and proximate (mechanistic) causes for such expansions of recombination suppression, including (1) antagonistic selection, (2) association of additional functions to mating-type, such as uniparental mitochondria inheritance, (3) accumulation in the margin of nonrecombining regions of various factors, including deleterious mutations or transposable elements resulting from relaxed selection, or neutral rearrangements resulting from genetic drift. The study of recombination suppression in fungi could thus contribute to our understanding of recombination suppression expansion across a broader range of organisms.

Genome-Wide Analyses of Repeat-Induced Point Mutations in the Ascomycota

The Repeat-Induced Point (RIP) mutation pathway is a fungus-specific genome defense mechanism that mitigates the deleterious consequences of repeated genomic regions and transposable elements (TEs). RIP mutates targeted sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence and extent of RIP with a sliding-window approach. Using genome-wide RIP data and two sets of control groups, the association between RIP, TEs, and GC content were contrasted in organisms capable and incapable of RIP. Based on these data, we then set out to determine the extent and occurrence of RIP in 58 representatives of the Ascomycota. The findings were summarized by placing each of the fungi investigated in one of six categories based on the extent of genome-wide RIP. In silico RIP analyses, using a sliding-window approach with stringent RIP parameters, implemented simultaneously within the same genetic context, on high quality genome assemblies, yielded superior results in determining the genome-wide RIP among the Ascomycota. Most Ascomycota had RIP and these mutations were particularly widespread among classes of the Pezizomycotina, including the early diverging Orbiliomycetes and the Pezizomycetes. The most extreme cases of RIP were limited to representatives of the Dothideomycetes and Sordariomycetes. By contrast, the genomes of the Taphrinomycotina and Saccharomycotina contained no detectable evidence of RIP. Also, recent losses in RIP combined with controlled TE proliferation in the Pezizomycotina subphyla may promote substantial genome enlargement as well as the formation of sub-genomic compartments. These findings have broadened our understanding of the taxonomic range and extent of RIP in Ascomycota and how this pathway affects the genomes of fungi harboring it.

Occurrence and possible roles of polysaccharides in fungi and their influence on the development of new technologies

The article summarizes the roles of polysaccharides in the biology of fungi and their relationship in the development of new technologies. The comparative approach between the evolution of fungi and the chemistry of glycobiology elucidated relevant aspects about the role of polysaccharides in fungi. Also, based on the knowledge of fungal glycobiology, it was possible to address the development of new technologies, such as the production of new anti-tumor drugs, vaccines, biomaterials, and applications in the field of robotics. We conclude that polysaccharides activate pathways of apoptosis, secretion of pro-inflammatory substances, and macrophage, inducing anticancer activity. Also, the activation of the immune system, which opens the way for the production of vaccines. The development of biomaterials and parts for robotics is a promising and little-explored field. Finally, the article is multidisciplinary, with a different and integrated approach to the role of nature in the sustainable development of new technologies.

Comparative analysis of genome-wide DNA methylation in Neurospora

DNA methylation is an epigenetic mark that plays an important role in genetic regulation in eukaryotes. Major progress has been made in dissecting the molecular pathways that regulate DNA methylation. Yet, little is known about DNA methylation variation over evolutionary time. Here we present an investigation of the variation of DNA methylation and transposable element (TE) content in species of the filamentous ascomycetes Neurospora. We generated genome-wide DNA methylation data at single-base resolution, together with genomic TE content and gene expression data, of 10 individuals representing five closely related Neurospora species. We found that the methylation levels were low (ranging from 1.3% to 2.5%) and varied among the genomes in a species-specific way. Furthermore, we found that the TEs over 400 bp long were targeted by DNA methylation, and in all genomes, high methylation correlated with low GC, confirming a conserved link between DNA methylation and Repeat Induced Point (RIP) mutations in this group of fungi. Both TE content and DNA methylation pattern showed phylogenetic signal, and the species with the highest TE load (N. crassa) also exhibited the highest methylation level per TE. Our results suggest that DNA methylation is an evolvable trait and indicate that the genomes of Neurospora are shaped by an evolutionary arms race between TEs and host defence.

Genomic Features of Cladobotryum dendroides, Which Causes Cobweb Disease in Edible Mushrooms, and Identification of Genes Related to Pathogenicity and Mycoparasitism

Cladobotryum dendroides, which causes cobweb disease in edible mushrooms, is one of the major fungal pathogens. Our previous studies focused on the genetic and morphological characterization of this fungus, as well as its pathogenicity and the identification of appropriate fungicides. However, little is known about the genome characters, pathogenic genes, and molecular pathogenic mechanisms of C. dendroides. Herein, we reported a high-quality de novo genomic sequence of C. dendroides and compared it with closely-related fungi. The assembled C. dendroides genome was 36.69 Mb, consisting of eight contigs, with an N50 of 4.76 Mb. This genome was similar in size to that of C. protrusum, and shared highly conserved syntenic blocks and a few inversions with C. protrusum. Phylogenetic analysis revealed that, within the Hypocreaceae, Cladobotryum was closer to Mycogone than to Trichoderma, which is consistent with phenotypic evidence. A significant number of the predicted expanded gene families were strongly associated with pathogenicity, virulence, and adaptation. Our findings will be instrumental for the understanding of fungi-fungi interactions, and for exploring efficient management strategies to control cobweb disease.

Repeat-Induced Point Mutations Drive Divergence between Fusarium circinatum and Its Close Relatives

The Repeat-Induced Point (RIP) mutation pathway is a fungal-specific genome defense mechanism that counteracts the deleterious effects of transposable elements. This pathway permanently mutates its target sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence of RIP in the pitch canker pathogen, Fusarium circinatum, and its close relatives in the Fusarium fujikuroi species complex (FFSC). Our results showed that the examined fungi all exhibited hallmarks of RIP, but that they differed in terms of the extent to which their genomes were affected by this pathway. RIP mutations constituted a large proportion of all the FFSC genomes, including both core and dispensable chromosomes, although the latter were generally more extensively affected by RIP. Large RIP-affected genomic regions were also much more gene sparse than the rest of the genome. Our data further showed that RIP-directed sequence diversification increased the variability between homologous regions of related species, and that RIP-affected regions can interfere with homologous recombination during meiosis, thereby contributing to post-mating segregation distortion. Taken together, these findings suggest that RIP can drive the independent divergence of chromosomes, alter chromosome architecture, and contribute to the divergence among F. circinatum and other members of this economically important group of fungi.

Programmed Cell Death in Neurospora crassa Is Controlled by the Allorecognition Determinant rcd-1

Nonself recognition following cell fusion between genetically distinct individuals of the same species in filamentous fungi often results in a programmed cell death (PCD) reaction, where the heterokaryotic fusion cell is compartmentalized and rapidly killed. The allorecognition process plays a key role as a defense mechanism that restricts genome exploitation, resource plundering, and the spread of deleterious senescence plasmids and mycoviruses. Although a number of incompatibility systems have been described that function in mature hyphae, less is known about the PCD pathways in asexual spores, which represent the main infectious unit in various human and plant fungal pathogens. Here, we report the identification of regulator of cell death-1 (rcd-1), a novel allorecognition gene, controlling PCD in germinating asexual spores of Neurospora crassa; rcd-1 is one of the most polymorphic genes in the genomes of wild N. crassa isolates. The coexpression of two antagonistic rcd-1-1 and rcd-1-2 alleles was necessary and sufficient to trigger cell death in fused germlings and in hyphae. Based on analysis of wild populations of N. crassa and N. discreta, rcd-1 alleles appeared to be under balancing selection and associated with trans-species polymorphisms. We shed light on genomic rearrangements that could have led to the emergence of the incompatibility system in Neurospora and show that rcd-1 belongs to a much larger gene family in fungi. Overall, our work contributes toward a better understanding of allorecognition and PCD in an underexplored developmental stage of filamentous fungi.

Evolution and functional characterization of pectate lyase PEL12, a member of a highly expanded Clonostachys rosea polysaccharide lyase 1 family

Background

Pectin is one of the major and most complex plant cell wall components that needs to be overcome by microorganisms as part of their strategies for plant invasion or nutrition. Microbial pectinolytic enzymes therefore play a significant role for plant-associated microorganisms and for the decomposition and recycling of plant organic matter. Recently, comparative studies revealed significant gene copy number expansion of the polysaccharide lyase 1 (PL1) pectin/pectate lyase gene family in the Clonostachys rosea genome, while only low numbers were found in Trichoderma species. Both of these fungal genera are widely known for their ability to parasitize and kill other fungi (mycoparasitism) and certain species are thus used for biocontrol of plant pathogenic fungi.

Results

In order to understand the role of the high number of pectin degrading enzymes in Clonostachys, we studied diversity and evolution of the PL1 gene family in C. rosea compared with other Sordariomycetes with varying nutritional life styles. Out of 17 members of C. rosea PL1, we could only detect two to be secreted at acidic pH. One of them, the pectate lyase pel12 gene was found to be strongly induced by pectin and, to a lower degree, by polygalacturonic acid. Heterologous expression of the PEL12 in a PL1-free background of T. reesei revealed direct enzymatic involvement of this protein in utilization of pectin at pH 5 without a requirement for Ca²⁺. The mutants showed increased utilization of pectin compounds, but did not increase biocontrol ability in detached leaf assay against the plant pathogen Botrytis cinerea compared to the wild type.

Conclusions

In this study, we aimed to gain insight into diversity and evolution of the PL1 gene family in C. rosea and other Sordariomycete species in relation to their nutritional modes. We show that C. rosea PL1 expansion does not correlate with its mycoparasitic nutritional mode and resembles those of strong plant pathogenic fungi. We further investigated regulation, specificity and function of the C. rosea PEL12 and show that this enzyme is directly involved in degradation of pectin and pectin-related compounds, but not in C. rosea biocontrol.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1310-9) contains supplementary material, which is available to authorized users.

Convergent evolution of complex genomic rearrangements in two fungal meiotic drive elements

Meiotic drive is widespread in nature. The conflict it generates is expected to be an important motor for evolutionary change and innovation. In this study, we investigated the genomic consequences of two large multi-gene meiotic drive elements, Sk-2 and Sk-3, found in the filamentous ascomycete Neurospora intermedia. Using long-read sequencing, we generated the first complete and well-annotated genome assemblies of large, highly diverged, non-recombining regions associated with meiotic drive elements. Phylogenetic analysis shows that, even though Sk-2 and Sk-3 are located in the same chromosomal region, they do not form sister clades, suggesting independent origins or at least a long evolutionary separation. We conclude that they have in a convergent manner accumulated similar patterns of tandem inversions and dense repeat clusters, presumably in response to similar needs to create linkage between genes causing drive and resistance.

Ceratocystidaceae exhibit high levels of recombination at the mating-type (MAT) locus

Mating is central to many fungal life cycles and is controlled by genes at the mating-type (MAT) locus. These genes determine whether the fungus will be self-sterile (heterothallic) or self-fertile (homothallic). Species in the ascomycete family Ceratocystidaceae have different mating strategies, making them interesting to consider with regards to their MAT loci. The aim of this study was to compare the composition of the MAT locus flanking regions in 11 species of Ceratocystidaceae representing four genera. Genome assemblies for each species were examined to identify the MAT locus and determine the structure of the flanking regions. Large contigs containing the MAT locus were then functionally annotated and analysed for the presence of transposable elements. Genes typically flanking the MAT locus in sordariomycetes were found to be highly conserved in the Ceratocystidaceae. The different genera in the Ceratocystidaceae displayed little synteny outside of the immediate MAT locus flanking genes. Even though species ofCeratocystis did not show much synteny outside of the immediate MAT locus flanking genes, species of Huntiella and Endoconidiophora were comparatively syntenic. Due to the high number of transposable elements present in Ceratocystis MAT flanking regions, we hypothesise that Ceratocystis species may have undergone recombination in this region.

Accessories Make the Outfit: Accessory Chromosomes and Other Dispensable DNA Regions in Plant-Pathogenic Fungi

Fungal pathogen genomes can often be divided into core and accessory regions. Accessory regions ARs) may be comprised of either ARs (within core chromosomes (CCs) or wholly dispensable (accessory) chromosomes (ACs). Fungal ACs and ARs typically accumulate mutations and structural rearrangements more rapidly over time than CCs and many harbor genes relevant to host-pathogen interactions. These regions are of particular interest in plant pathology and include host-specific virulence factors and secondary metabolite synthesis gene clusters. This review outlines known ACs and ARs in fungal genomes, methods used for their detection, their common properties that differentiate them from the core genome, and what is currently known of their various roles in pathogenicity. Reports on the evolutionary processes generating and shaping AC and AR compartments are discussed, including repeat induced point mutation and breakage fusion bridge cycles. Previously ACs have been studied extensively within key genera, including Fusarium, Zymoseptoria, and Alternaria, but are growing in frequency of observation and perceived importance across a wider range of fungal species. Recent advances in sequencing technologies permit affordable genome assembly and resequencing of populations that will facilitate further discovery and routine screening of ACs.

Evolution of recombination rates between sex chromosomes

In species with genetic sex-determination, the chromosomes carrying the sex-determining genes have often evolved non-recombining regions and subsequently evolved the full set of characteristics denoted by the term 'sex chromosomes'. These include size differences, creating chromosomal heteromorphism, and loss of gene functions from one member of the chromosome pair. Such characteristics and changes have been widely reviewed, and underlie molecular genetic approaches that can detect sex chromosome regions. This review deals mainly with the evolution of new non-recombining regions, focusing on how certain evolutionary situations select for suppressed recombination (rather than the proximate mechanisms causing suppressed recombination between sex chromosomes). Particularly important is the likely involvement of sexually antagonistic polymorphisms in genome regions closely linked to sex-determining loci. These may be responsible for the evolutionary strata of sex chromosomes that have repeatedly formed by recombination suppression evolving across large genome regions. More studies of recently evolved non-recombining sex-determining regions should help to test this hypothesis empirically, and may provide evidence about whether other situations can sometimes lead to sex-linked regions evolving. Similarities with other non-recombining genome regions are discussed briefly, to illustrate common features of the different cases, though no general properties apply to all of them.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

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.

Large-scale suppression of recombination predates genomic rearrangements in Neurospora tetrasperma

A common feature of eukaryote genomes is large chromosomal regions where recombination is absent or strongly reduced, but the factors that cause this reduction are not well understood. Genomic rearrangements have often been implicated, but they may also be a consequence of recombination suppression rather than a cause. In this study, we generate eight high-quality genomic data sets of the filamentous ascomycete Neurospora tetrasperma, a fungus that lacks recombination over most of its largest chromosome. The genomes surprisingly reveal collinearity of the non-recombining regions and although large inversions are enriched in these regions, we conclude these inversions to be derived and not the cause of the suppression. To our knowledge, this is the first time that non-recombining, genic regions as large as 86% of a full chromosome (or 8 Mbp), are shown to be collinear. These findings are of significant interest for our understanding of the evolution of sex chromosomes and other supergene complexes.

Genomic structure and evolution of the mating type locus in the green seaweed Ulva partita

The evolution of sex chromosomes and mating loci in organisms with UV systems of sex/mating type determination in haploid phases via genes on UV chromosomes is not well understood. We report the structure of the mating type (MT) locus and its evolutionary history in the green seaweed Ulva partita, which is a multicellular organism with an isomorphic haploid-diploid life cycle and mating type determination in the haploid phase. Comprehensive comparison of a total of 12.0 and 16.6 Gb of genomic next-generation sequencing data for mt- and mt+ strains identified highly rearranged MT loci of 1.0 and 1.5 Mb in size and containing 46 and 67 genes, respectively, including 23 gametologs. Molecular evolutionary analyses suggested that the MT loci diverged over a prolonged period in the individual mating types after their establishment in an ancestor. A gene encoding an RWP-RK domain-containing protein was found in the mt- MT locus but was not an ortholog of the chlorophycean mating type determination gene MID. Taken together, our results suggest that the genomic structure and its evolutionary history in the U. partita MT locus are similar to those on other UV chromosomes and that the MT locus genes are quite different from those of Chlorophyceae.

A-to-I RNA editing is developmentally regulated and generally adaptive for sexual reproduction in Neurospora crassa

Although fungi lack adenosine deaminase acting on RNA (ADAR) enzymes, adenosine to inosine (A-to-I) RNA editing was reported recently in Fusarium graminearum during sexual reproduction. In this study, we profiled the A-to-I editing landscape and characterized its functional and adaptive properties in the model filamentous fungus Neurospora crassa A total of 40,677 A-to-I editing sites were identified, and approximately half of them displayed stage-specific editing or editing levels at different sexual stages. RNA-sequencing analysis with the Δstc-1 and Δsad-1 mutants confirmed A-to-I editing occurred before ascus development but became more prevalent during ascosporogenesis. Besides fungal-specific sequence and secondary structure preference, 63.5% of A-to-I editing sites were in the coding regions and 81.3% of them resulted in nonsynonymous recoding, resulting in a significant increase in the proteome complexity. Many genes involved in RNA silencing, DNA methylation, and histone modifications had extensive recoding, including sad-1, sms-3, qde-1, and dim-2. Fifty pseudogenes harbor premature stop codons that require A-to-I editing to encode full-length proteins. Unlike in humans, nonsynonymous editing events in N. crassa are generally beneficial and favored by positive selection. Almost half of the nonsynonymous editing sites in N. crassa are conserved and edited in Neurospora tetrasperma Furthermore, hundreds of them are conserved in F. graminearum and had higher editing levels. Two unknown genes with editing sites conserved between Neurospora and Fusarium were experimentally shown to be important for ascosporogenesis. This study comprehensively analyzed A-to-I editing in N. crassa and showed that RNA editing is stage-specific and generally adaptive, and may be functionally related to repeat induced point mutation and meiotic silencing by unpaired DNA.

Fungal Sex: The Ascomycota

This article provides an overview of sexual reproduction in the ascomycetes, a phylum of fungi that is named after the specialized sacs or "asci" that hold the sexual spores. They have therefore also been referred to as the Sac Fungi due to these characteristic structures that typically contain four to eight ascospores. Ascomycetes are morphologically diverse and include single-celled yeasts, filamentous fungi, and more complex cup fungi. The sexual cycles of many species, including those of the model yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe and the filamentous saprobes Neurospora crassa, Aspergillus nidulans, and Podospora anserina, have been examined in depth. In addition, sexual or parasexual cycles have been uncovered in important human pathogens such as Candida albicans and Aspergillus fumigatus, as well as in plant pathogens such as Fusarium graminearum and Cochliobolus heterostrophus. We summarize what is known about sexual fecundity in ascomycetes, examine how structural changes at the mating-type locus dictate sexual behavior, and discuss recent studies that reveal that pheromone signaling pathways can be repurposed to serve cellular roles unrelated to sex.

Diverse data supports the transition of filamentous fungal model organisms into the post-genomics era

Filamentous fungi have been important as model organisms since the beginning of modern biological inquiry and have benefitted from open data since the earliest genetic maps were shared. From early origins in simple Mendelian genetics of mating types, parasexual genetics of colony colour, and the foundational demonstration of the segregation of a nutritional requirement, the contribution of research systems utilising filamentous fungi has spanned the biochemical genetics era, through the molecular genetics era, and now are at the very foundation of diverse omics approaches to research and development. Fungal model organisms have come from most major taxonomic groups although Ascomycete filamentous fungi have seen the most major sustained effort. In addition to the published material about filamentous fungi, shared molecular tools have found application in every area of fungal biology. Similarly, shared data has contributed to the success of model systems. The scale of data supporting research with filamentous fungi has grown by 10 to 12 orders of magnitude. From genetic to molecular maps, expression databases, and finally genome resources, the open and collaborative nature of the research communities has assured that the rising tide of data has lifted all of the research systems together.

OcculterCut: A Comprehensive Survey of AT-Rich Regions in Fungal Genomes

We present a novel method to measure the local GC-content bias in genomes and a survey of published fungal species. The method, enacted as "OcculterCut" (https://sourceforge.net/projects/occultercut, last accessed April 30, 2016), identified species containing distinct AT-rich regions. In most fungal taxa, AT-rich regions are a signature of repeat-induced point mutation (RIP), which targets repetitive DNA and decreases GC-content though the conversion of cytosine to thymine bases. RIP has in turn been identified as a driver of fungal genome evolution, as RIP mutations can also occur in single-copy genes neighboring repeat-rich regions. Over time RIP perpetuates "two speeds" of gene evolution in the GC-equilibrated and AT-rich regions of fungal genomes. In this study, genomes showing evidence of this process are found to be common, particularly among the Pezizomycotina. Further analysis highlighted differences in amino acid composition and putative functions of genes from these regions, supporting the hypothesis that these regions play an important role in fungal evolution. OcculterCut can also be used to identify genes undergoing RIP-assisted diversifying selection, such as small, secreted effector proteins that mediate host-microbe disease interactions.

The Fast-Evolving phy-2 Gene Modulates Sexual Development in Response to Light in the Model Fungus Neurospora crassa

Rapid responses to changes in incident light are critical to the guidance of behavior and development in most species. Phytochrome light receptors in particular play key roles in bacterial physiology and plant development, but their functions and regulation are less well understood in fungi. Nevertheless, genome-wide expression measurements provide key information that can guide experiments that reveal how genes respond to environmental signals and clarify their role in development. We performed functional genomic and phenotypic analyses of the two phytochromes in Neurospora crassa, a fungal model adapted to a postfire environment that experiences dramatically variable light conditions. Expression of phy-1 and phy-2 was low in early sexual development and in the case of phy-2 increased in late sexual development. Under light stimulation, strains with the phytochromes deleted exhibited increased expression of sexual development-related genes. Moreover, under red light, the phy-2 knockout strain commenced sexual development early. In the evolution of phytochromes within ascomycetes, at least two duplications have occurred, and the faster-evolving phy-2 gene has frequently been lost. Additionally, the three key cysteine sites that are critical for bacterial and plant phytochrome function are not conserved within fungal phy-2 homologs. Through the action of phytochromes, transitions between asexual and sexual reproduction are modulated by light level and light quality, presumably as an adaptation for fast asexual growth and initiation of sexual reproduction of N. crassa in exposed postfire ecosystems.

Environmental signals, including light, play critical roles in regulating fungal growth and pathogenicity, and balance of asexual and sexual reproduction is critical in fungal pathogens’ incidence, virulence, and distribution. Red light sensing by phytochromes is well known to play critical roles in bacterial physiology and plant development. Homologs of phytochromes were first discovered in the fungal model Neurospora crassa and then subsequently in diverse other fungi, including many plant pathogens. Our study investigated the evolution of red light sensors in ascomycetes and confirmed—using the model fungus Neurospora crassa—their roles in modulating the asexual-sexual reproduction balance in fungi. Our findings also provide a key insight into one of the most poorly understood aspects of fungal biology, suggesting that further study of the function of phytochromes in fungi is critical to reveal the genetic basis of the asexual-sexual switch responsible for fungal growth and distribution, including diverse and destructive plant pathogens.

Introgression maintains the genetic integrity of the mating-type determining chromosome of the fungus Neurospora tetrasperma

Genome evolution is driven by a complex interplay of factors, including selection, recombination, and introgression. The regions determining sexual identity are particularly dynamic parts of eukaryotic genomes that are prone to molecular degeneration associated with suppressed recombination. In the fungus Neurospora tetrasperma, it has been proposed that this molecular degeneration is counteracted by the introgression of nondegenerated DNA from closely related species. In this study, we used comparative and population genomic analyses of 92 genomes from eight phylogenetically and reproductively isolated lineages of N. tetrasperma, and its three closest relatives, to investigate the factors shaping the evolutionary history of the genomes.We found that suppressed recombination extends across at least 6 Mbp (∼ 63%) of the mating-type (mat) chromosome in N. tetrasperma and is associated with decreased genetic diversity, which is likely the result primarily of selection at linked sites. Furthermore, analyses of molecular evolution revealed an increased mutational load in this region, relative to recombining regions. However, comparative genomic and phylogenetic analyses indicate that the mat chromosomes are temporarily regenerated via introgression from sister species; six of eight lineages show introgression into one of their mat chromosomes, with multiple Neurospora species acting as donors. The introgressed tracts have been fixed within lineages, suggesting that they confer an adaptive advantage in natural populations, and our analyses support the presence of selective sweeps in at least one lineage. Thus, these data strongly support the previously hypothesized role of introgression as a mechanism for the maintenance of mating-type determining chromosomal regions.

Evolution of sexual reproduction: a view from the Fungal Kingdom supports an evolutionary epoch with sex before sexes

Sexual reproduction is conserved throughout each supergroup within the eukaryotic tree of life, and therefore thought to have evolved once and to have been present in the last eukaryotic common ancestor (LECA). Given the antiquity of sex, there are features of sexual reproduction that are ancient and ancestral, and thus shared in diverse extant organisms. On the other hand, the vast evolutionary distance that separates any given extant species from the LECA necessarily implies that other features of sex will be derived. While most types of sex we are familiar with involve two opposite sexes or mating types, recent studies in the fungal kingdom have revealed novel and unusual patterns of sexual reproduction, including unisexual reproduction. In this mode of reproduction a single mating type can on its own undergo self-fertile/homothallic reproduction, either with itself or with other members of the population of the same mating type. Unisexual reproduction has arisen independently as a derived feature in several different lineages. That a myriad of different types of sex determination and sex determinants abound in animals, plants, protists, and fungi suggests that sex specification itself may not be ancestral and instead may be a derived trait. If so, then the original form of sexual reproduction may have been unisexual, onto which sexes were superimposed as a later feature. In this model, unisexual reproduction is both an ancestral and a derived trait. In this review, we consider what is new and what is old about sexual reproduction from the unique vantage point of the fungal kingdom.

Contrasted patterns in mating-type chromosomes in fungi: hotspots versus coldspots of recombination

It is striking that, while central to sexual reproduction, the genomic regions determining sex or mating-types are often characterized by suppressed recombination that leads to a decrease in the efficiency of selection, shelters genetic load, and inevitably contributes to their genic degeneration. Research on model and lesser-explored fungi has revealed similarities in recombination suppression of the genomic regions involved in mating compatibility across eukaryotes, but fungi also provide opposite examples of enhanced recombination in the genomic regions that determine their mating types. These contrasted patterns of genetic recombination (sensu lato, including gene conversion and ectopic recombination) in regions of the genome involved in mating compatibility point to important yet complex processes occurring in their evolution. A number of pieces in this puzzle remain to be solved, in particular on the unclear selective forces that may cause the patterns of recombination, prompting theoretical developments and experimental studies. This review thus points to fungi as a fascinating group for studying the various evolutionary forces at play in the genomic regions involved in mating compatibility.

Clonal reproduction in fungi

Research over the past two decades shows that both recombination and clonality are likely to contribute to the reproduction of all fungi. This view of fungi is different from the historical and still commonly held view that a large fraction of fungi are exclusively clonal and that some fungi have been exclusively clonal for hundreds of millions of years. Here, we first will consider how these two historical views have changed. Then we will examine the impact on fungal research of the concept of restrained recombination [Tibayrenc M, Ayala FJ (2012) Proc Natl Acad Sci USA 109 (48):E3305-E3313]. Using animal and human pathogenic fungi, we examine extrinsic restraints on recombination associated with bottlenecks in genetic variation caused by geographic dispersal and extrinsic restraints caused by shifts in reproductive mode associated with either disease transmission or hybridization. Using species of the model yeast Saccharomyces and the model filamentous fungus Neurospora, we examine intrinsic restraints on recombination associated with mating systems that range from strictly clonal at one extreme to fully outbreeding at the other and those that lie between, including selfing and inbreeding. We also consider the effect of nomenclature on perception of reproductive mode and a means of comparing the relative impact of clonality and recombination on fungal populations. Last, we consider a recent hypothesis suggesting that fungi thought to have the most severe intrinsic constraints on recombination actually may have the fewest.

Maintaining heterokaryosis in pseudo-homothallic fungi

Among all the strategies displayed by fungi to reproduce and propagate, some species have adopted a peculiar behavior called pseudo-homothallism. Pseudo-homothallic fungi are true heterothallics, i.e., they need 2 genetically-compatible partners to mate, but they produce self-fertile mycelium in which the 2 different nuclei carrying the compatible mating types are present. This lifestyle not only enables the fungus to reproduce without finding a compatible partner, but also to cross with any mate it may encounter. However, to be fully functional, pseudo-homothallism requires maintaining heterokaryosis at every stage of the life cycle. We recently showed that neither the structure of the mating-type locus nor hybrid-enhancing effect due to the presence of the 2 mating types accounts for the maintenance of heterokaryosis in the pseudo-homothallic fungus P. anserina. In this addendum, we summarize the mechanisms creating heterokaryosis in P. anserina and 2 other well-known pseudo-homothallic fungi, Neurospora tetrasperma and Agaricus bisporus. We also discuss mechanisms potentially involved in maintaining heterokaryosis in these 3 species.

The status of supergenes in the 21st century: recombination suppression in Batesian mimicry and sex chromosomes and other complex adaptations

I review theoretical models for the evolution of supergenes in the cases of Batesian mimicry in butterflies, distylous plants and sex chromosomes. For each of these systems, I outline the genetic evidence that led to the proposal that they involve multiple genes that interact during ‘complex adaptations’, and at which the mutations involved are not unconditionally advantageous, but show advantages that trade‐off against some disadvantages. I describe recent molecular genetic studies of these systems and questions they raise about the evolution of suppressed recombination. Nonrecombining regions of sex chromosomes have long been known, but it is not yet fully understood why recombination suppression repeatedly evolved in systems in distantly related taxa, but does not always evolve. Recent studies of distylous plants are tending to support the existence of recombination‐suppressed genome regions, which may include modest numbers of genes and resemble recently evolved sex‐linked regions. For Batesian mimicry, however, molecular genetic work in two butterfly species suggests a new supergene scenario, with a single gene mutating to produce initial adaptive phenotypes, perhaps followed by modifiers specifically refining and perfecting the new phenotype.

Endogenous Small RNA Mediates Meiotic Silencing of a Novel DNA Transposon

Genome defense likely evolved to curtail the spread of transposable elements and invading viruses. A combination of effective defense mechanisms has been shown to limit colonization of the Neurospora crassa genome by transposable elements. A novel DNA transposon named Sly1-1 was discovered in the genome of the most widely used laboratory "wild-type" strain FGSC 2489 (OR74A). Meiotic silencing by unpaired DNA, also simply called meiotic silencing, prevents the expression of regions of the genome that are unpaired during karyogamy. This mechanism is posttranscriptional and is proposed to involve the production of small RNA, so-called masiRNAs, by proteins homologous to those involved in RNA interference-silencing pathways in animals, fungi, and plants. Here, we demonstrate production of small RNAs when Sly1-1 was unpaired in a cross between two wild-type strains. These small RNAs are dependent on SAD-1, an RNA-dependent RNA polymerase necessary for meiotic silencing. We present the first case of endogenously produced masiRNA from a novel N. crassa DNA transposable element.

Chaos of Rearrangements in the Mating-Type Chromosomes of the Anther-Smut Fungus Microbotryum lychnidis-dioicae

Sex chromosomes in plants and animals and fungal mating-type chromosomes often show exceptional genome features, with extensive suppression of homologous recombination and cytological differentiation between members of the diploid chromosome pair. Despite strong interest in the genetics of these chromosomes, their large regions of suppressed recombination often are enriched in transposable elements and therefore can be challenging to assemble. Here we show that the latest improvements of the PacBio sequencing yield assembly of the whole genome of the anther-smut fungus, Microbotryum lychnidis-dioicae (the pathogenic fungus causing anther-smut disease of Silene latifolia), into finished chromosomes or chromosome arms, even for the repeat-rich mating-type chromosomes and centromeres. Suppressed recombination of the mating-type chromosomes is revealed to span nearly 90% of their lengths, with extreme levels of rearrangements, transposable element accumulation, and differentiation between the two mating types. We observed no correlation between allelic divergence and physical position in the nonrecombining regions of the mating-type chromosomes. This may result from gene conversion or from rearrangements of ancient evolutionary strata, i.e., successive steps of suppressed recombination. Centromeres were found to be composed mainly of copia-like transposable elements and to possess specific minisatellite repeats identical between the different chromosomes. We also identified subtelomeric motifs. In addition, extensive signs of degeneration were detected in the nonrecombining regions in the form of transposable element accumulation and of hundreds of gene losses on each mating-type chromosome. Furthermore, our study highlights the potential of the latest breakthrough PacBio chemistry to resolve complex genome architectures.

Neurospora Heterokaryons with Complementary Duplications and Deficiencies in Their Constituent Nuclei Provide an Approach to Identify Nucleus-Limited Genes

Introgression is the transfer of genes or genomic regions from one species into another via hybridization and back-crosses. We have introgressed four translocations (EB4, IBj5, UK14-1, and B362i) from Neurospora crassa into N. tetrasperma. This enabled us to construct two general types of heterokaryons with mat-A and mat-a nuclei of different genotypes: one type is [T + N] (with one translocation nucleus and one normal sequence nucleus), and the other is [Dp + Df] (with one nucleus carrying a duplication of the translocation region and the other being deleted for the translocation region). Self-crossing these heterokaryons again produced [T + N] and [Dp + Df] progeny. From conidia (vegetative spores) produced by the heterokaryotic mycelia, we obtained self-fertile (heterokaryotic) and self-sterile (homokaryotic) derivative strains. [T + N] heterokaryons produced homokaryotic conidial derivatives of both mating types, but [Dp + Df] heterokaryons produced viable conidial homokaryons of only the mating type of the Dp nucleus. All four [T + N] heterokaryons and three [Dp + Df] heterokaryons produced both self-sterile and self-fertile conidial derivatives, but the [Dp(B362i) + Df(B362i)] heterokaryons produced only self-sterile ones. Conceivably, the Df(B362i) nuclei may be deleted for a nucleus-limited gene required for efficient mitosis or nuclear division, and whose deficit is not complemented by the neighboring Dp(B362i) nuclei. A cross involving Dp(EB4) showed repeat-induced point mutation (RIP). Because RIP can occur in self-crosses of [Dp + Df] but not [T + N] heterokaryons, RIP alteration of a translocated segment would depend on the relative numbers of [Dp + Df] vs. [T + N] ancestors.

Degeneration of the nonrecombining regions in the mating-type chromosomes of the anther-smut fungi

Dimorphic mating-type chromosomes in fungi are excellent models for understanding the genomic consequences of recombination suppression. Their suppressed recombination and reduced effective population size are expected to limit the efficacy of natural selection, leading to genomic degeneration. Our aim was to identify the sequences of the mating-type chromosomes (a1 and a2) of the anther-smut fungi and to investigate degeneration in their nonrecombining regions. We used the haploid a1 Microbotryum lychnidis-dioicae reference genome sequence. The a1 and a2 mating-type chromosomes were both isolated electrophoretically and sequenced. Integration with restriction-digest optical maps identified regions of recombination and nonrecombination in the mating-type chromosomes. Genome sequence data were also obtained for 12 other Microbotryum species. We found strong evidence of degeneration across the genus in the nonrecombining regions of the mating-type chromosomes, with significantly higher rates of nonsynonymous substitution (dN/dS) than in nonmating-type chromosomes or in recombining regions of the mating-type chromosomes. The nonrecombining regions of the mating-type chromosomes also showed high transposable element content, weak gene expression, and gene losses. The levels of degeneration did not differ between the a1 and a2 mating-type chromosomes, consistent with the lack of homogametic/heterogametic asymmetry between them, and contrasting with X/Y or Z/W sex chromosomes.

Recent and massive expansion of the mating-type-specific region in the smut fungus Microbotryum

The presence of large genomic regions with suppressed recombination (SR) is a key shared property of some sex- and mating-type determining (mat) chromosomes identified to date in animals, plants, and fungi. Why such regions form and how they evolve remain central questions in evolutionary genetics. The smut fungus Microbotryum lychnis-dioicae is a basidiomycete fungus in which dimorphic mat chromosomes have been reported, but the size, age, and evolutionary dynamics of the SR region remains unresolved. To identify the SR region in M. lychnis-dioicae and to study its evolution, we sequenced 12 genomes (6 per mating type) of this species and identified the genomic contigs that show fixed sequence differences between the mating types. We report that the SR region spans more than half of the mat chromosome (>2.3 Mbp) and that it is of very recent origin (∼2 × 10(6) years) as the average sequence divergence between mating types was only 2% in the SR region. This contrasts with a much higher divergence in and around the mating-type determining pheromone receptor locus in the SR, suggesting a recent and massive expansion of the SR region. Our results comprise the first reported case of recent massive SR expansion documented in a basidiomycete fungus.

Evolutionary Perspectives on Human Fungal Pathogens

Owing to their small size and paucity of phenotypic characters, progress in the evolutionary biology of microbes in general, and human pathogenic fungi in particular, has been linked to a series of advances in DNA sequencing over the past quarter century. Phylogenetics was the first area to benefit, with the achievement of a basic understanding of fungal phylogeny. Population genetics was the next advance, finding cryptic species everywhere, and recombination in species previously thought to be asexual. Comparative genomics saw the next advance, in which variation in gene content and changes in gene family size were found to be important sources of variation. Fungal population genomics is showing that gene flow among closely related populations and species provides yet another source of adaptive, genetic variation. Now, two means to associate genetic variation with phenotypic variation, "reverse ecology" for adaptive phenotypes, and genome-wide association of any phenotype, are letting evolutionary biology make a profound contribution to molecular developmental biology of pathogenic fungi.

Intron evolution in Neurospora: the role of mutational bias and selection

We used comparative and population genomics to study intron evolutionary dynamics in the fungal model genus Neurospora. For our investigation, we used well-annotated genomes of N. crassa, N. discreta, and N. tetrasperma, and 92 resequenced genomes of N. tetrasperma from natural populations. By analyzing the four well-annotated genomes, we identified 9495 intron sites in 7619 orthologous genes. Our data supports nonhomologous end joining (NHEJ) and tandem duplication as mechanisms for intron gains in the genus and the RT-mRNA process as a mechanism for intron loss. We found a moderate intron gain rate (5.78–6.89 × 10⁻¹³ intron gains per nucleotide site per year) and a high intron loss rate (7.53–13.76 × 10⁻¹⁰ intron losses per intron sites per year) as compared to other eukaryotes. The derived intron gains and losses are skewed to high frequencies, relative to neutral SNPs, in natural populations of N. tetrasperma, suggesting that selection is involved in maintaining a high intron turnover. Furthermore, our analyses of the association between intron population-level frequency and genomic features suggest that selection is involved in shaping a 5′ intron position bias and a low intron GC content. However, intron sequence analyses suggest that the gained introns were not exposed to recent selective sweeps. Taken together, this work contributes to our understanding of the importance of mutational bias and selection in shaping the intron distribution in eukaryotic genomes.

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.

Prevalence of transcription factors in ascomycete and basidiomycete fungi

Background

Gene regulation underlies fungal physiology and therefore is a major factor in fungal biodiversity. Analysis of genome sequences has revealed a large number of putative transcription factors in most fungal genomes. The presence of fungal orthologs for individual regulators has been analysed and appears to be highly variable with some regulators widely conserved and others showing narrow distribution. Although genome-scale transcription factor surveys have been performed before, no global study into the prevalence of specific regulators across the fungal kingdom has been presented.

Results

In this study we have analysed the number of members for 37 regulator classes in 77 ascomycete and 31 basidiomycete fungal genomes and revealed significant differences between ascomycetes and basidiomycetes. In addition, we determined the presence of 64 regulators characterised in ascomycetes across these 108 genomes. This demonstrated that overall the highest presence of orthologs is in the filamentous ascomycetes. A significant number of regulators lacked orthologs in the ascomycete yeasts and the basidiomycetes. Conversely, of seven basidiomycete regulators included in the study, only one had orthologs in ascomycetes.

Conclusions

This study demonstrates a significant difference in the regulatory repertoire of ascomycete and basidiomycete fungi, at the level of both regulator class and individual regulator. This suggests that the current regulatory systems of these fungi have been mainly developed after the two phyla diverged. Most regulators detected in both phyla are involved in central functions of fungal physiology and therefore were likely already present in the ancestor of the two phyla.

Maintaining two mating types: structure of the mating type locus and its role in heterokaryosis in Podospora anserina

Pseudo-homothallism is a reproductive strategy elected by some fungi producing heterokaryotic sexual spores containing genetically different but sexually compatible nuclei. This lifestyle appears as a compromise between true homothallism (self-fertility with predominant inbreeding) and complete heterothallism (with exclusive outcrossing). However, pseudohomothallic species face the problem of maintaining heterokaryotic mycelia to fully benefit from this lifestyle, as homokaryons are self-sterile. Here, we report on the structure of chromosome 1 in mat+ and mat- isolates of strain S of the pseudohomothallic fungus Podospora anserina. Chromosome 1 contains either one of the mat+ and mat- mating types of P. anserina, which is mostly found in nature as a mat+/mat- heterokaryotic mycelium harboring sexually compatible nuclei. We identified a "mat" region ∼0.8 Mb long, devoid of meiotic recombination and containing the mating-type idiomorphs, which is a candidate to be involved in the maintenance of the heterokaryotic state, since the S mat+ and S mat- strains have different physiology that may enable hybrid-vigor-like phenomena in the heterokaryons. The mat region contains 229 coding sequences. A total of 687 polymorphisms were detected between the S mat+ and S mat- chromosomes. Importantly, the mat region is colinear between both chromosomes, which calls for an original mechanism of recombination inhibition. Microarray analyses revealed that 10% of the P. anserina genes have different transcriptional profiles in S mat+ and S mat-, in line with their different phenotypes. Finally, we show that the heterokaryotic state is faithfully maintained during mycelium growth of P. anserina, yet mat+/mat+ and mat-/mat- heterokaryons are as stable as mat+/mat- ones, evidencing a maintenance of heterokaryosis that does not rely on fitness-enhancing complementation between the S mat+ and S mat- strains.

Species and population level molecular profiling reveals cryptic recombination and emergent asymmetry in the dimorphic mating locus of C. reinhardtii

Heteromorphic sex-determining regions or mating-type loci can contain large regions of non-recombining sequence where selection operates under different constraints than in freely recombining autosomal regions. Detailed studies of these non-recombining regions can provide insights into how genes are gained and lost, and how genetic isolation is maintained between mating haplotypes or sex chromosomes. The Chlamydomonas reinhardtii mating-type locus (MT) is a complex polygenic region characterized by sequence rearrangements and suppressed recombination between its two haplotypes, MT+ and MT−. We used new sequence information to redefine the genetic contents of MT and found repeated translocations from autosomes as well as sexually controlled expression patterns for several newly identified genes. We examined sequence diversity of MT genes from wild isolates of C. reinhardtii to investigate the impacts of recombination suppression. Our population data revealed two previously unreported types of genetic exchange in Chlamydomonas MT—gene conversion in the rearranged domains, and crossover exchanges in flanking domains—both of which contribute to maintenance of genetic homogeneity between haplotypes. To investigate the cause of blocked recombination in MT we assessed recombination rates in crosses where the parents were homozygous at MT. While normal recombination was restored in MT+×MT+ crosses, it was still suppressed in MT−×MT− crosses. These data revealed an underlying asymmetry in the two MT haplotypes and suggest that sequence rearrangements are insufficient to fully account for recombination suppression. Together our findings reveal new evolutionary dynamics for mating loci and have implications for the evolution of heteromorphic sex chromosomes and other non-recombining genomic regions.

Sex chromosomes and mating-type loci are often atypical in their structure and evolutionary dynamics. One distinguishing feature is the absence of recombination that results in genetic isolation and promotes rapid evolution and sometimes degeneration. We investigated gene content, sex-regulated expression, and recombination of mating locus (MT) genes in the unicellular alga Chlamydomonas reinhardtii. Despite the lack of observable recombination in and around Chlamydomonas MT, genes from its two mating types are far more similar to each other than expected for a non-recombining region. This discrepancy is explained by our finding evidence of genetic exchange between the two mating types within wild populations. In addition, we observed an unexpected asymmetry in the recombination behavior of the two mating types that may have contributed to the preferential expansion of one MT haplotype over the other through insertion of new genes. Our data suggest a mechanism to explain the emergence of heteromorphic sex chromosomes in haploid organisms by asymmetric expansion rather than by loss or degeneration as occurs in some Y or W chromosomes from diploid organisms. Our observations support a revised view of recombination in sex-determining regions as a quantitative phenomenon that can significantly affect rates of evolution and sex-linked genetic diversification.

Neurospora and the dead-end hypothesis: genomic consequences of selfing in the model genus

It is becoming increasingly evident that adoption of different reproductive strategies, such as sexual selfing and asexuality, greatly impacts genome evolution. In this study, we test theoretical predictions on genomic maladaptation of selfing lineages using empirical data from the model fungus Neurospora. We sequenced the genomes of four species representing distinct transitions to selfing within the history of the genus, as well as the transcriptome of one of these, and compared with available data from three outcrossing species. Our results provide evidence for a relaxation of purifying selection in protein-coding genes and for a reduced efficiency of transposable element silencing by Repeat Induced Point mutation. A reduction in adaptive evolution was also identified in the form of reduced codon usage bias in highly expressed genes of selfing Neurospora, but this result may be confounded by mutational bias. Potentially counteracting these negative effects, the nucleotide substitution rate and the spread of transposons is reduced in selfing species. We suggest that differences in substitution rate relate to the absence, in selfing Neurospora, of the asexual pathway producing conidia. Our results support the dead-end theory and show that Neurospora genomes bear signatures of both sexual and asexual reproductive mode.

Sex-linked transcriptional divergence in the hermaphrodite fungus Neurospora tetrasperma

In the filamentous ascomycete Neurospora tetrasperma, a large (approx. 7 Mbp) region of suppressed recombination surrounds the mating-type (mat) locus. While the remainder of the genome is largely homoallelic, this region of recombinational suppression, extending over 1500 genes, is associated with sequence divergence. Here, we used microarrays to examine how the molecular phenotype of gene expression level is linked to this divergent region, and thus to the mating type. Culturing N. tetrasperma on agar media that induce sexual/female or vegetative/male tissue, we found 196 genes significantly differentially expressed between mat A and mat a mating types. Our data show that the genes exhibiting mat-linked expression are enriched in the region genetically linked to mating type, and sequence and expression divergence are positively correlated. Our results indicate that the phenotype of mat A strains is optimized for traits promoting sexual/female development and the phenotype of mat a strains for vegetative/male development. This discovery of differentially expressed genes associated with mating type provides a link between genotypic and phenotypic divergence in this taxon and illustrates a fungal analogue to sexual dimorphism found among animals and plants.

Regional control of histone H3 lysine 27 methylation in Neurospora

Trimethylated lysine 27 on histone H3 (H3K27me3) is present in Drosophila, Arabidopsis, worms, and mammals, but is absent from yeasts that have been examined. We identified and analyzed H3K27me3 in the filamentous fungus Neurospora crassa and in other Neurospora species. H3K27me3 covers 6.8% of the N. crassa genome, encompassing 223 domains, including 774 genes, all of which are transcriptionally silent. N. crassa H3K27me3-marked genes are less conserved than unmarked genes and only ∼35% of genes marked by H3K27me3 in N. crassa are also H3K27me3-marked in Neurospora discreta and Neurospora tetrasperma. We found that three components of the Neurospora Polycomb repressive complex 2 (PRC2)--[Su-(var)3-9; E(z); Trithorax] (SET)-7, embryonic ectoderm development (EED), and SU(Z)12 (suppressor of zeste12)--are required for H3K27me3, whereas the fourth component, Neurospora protein 55 (an N. crassa homolog of p55/RbAp48), is critical for H3K27me3 only at subtelomeric domains. Loss of H3K27me3, caused by deletion of the gene encoding the catalytic PRC2 subunit, set-7, resulted in up-regulation of 130 genes, including genes in both H3K27me3-marked and unmarked regions.

Plant sex chromosome evolution

It is now well established that plants have an important place in studies of sex chromosome evolution because of the repeated independent evolution of separate sexes and sex chromosomes. There has been considerable recent progress in studying plant sex chromosomes. In this review, I focus on how these recent studies have helped clarify or answer several important questions about sex chromosome evolution, and I shall also try to clarify some common misconceptions. I also outline future work that will be needed to make further progress, including testing some important ideas by genetic, molecular, and developmental approaches. Systems with different ages can clearly help show the time course of events during changes from an ancestral co-sexual state (hermaphroditism or monoecy), and I will also explain how different questions can be studied in lineages whose dioecy or sex chromosomes evolved at different times in the past.