Filamentous fungi are of great importance in ecology, agriculture, medicine, and biotechnology. Thus, it is not surprising that genomes for more than 100 filamentous fungi have been sequenced, most of them by Sanger sequencing. While next-generation sequencing techniques have revolutionized genome resequencing, e.g. for strain comparisons, genetic mapping, or transcriptome and ChIP analyses, de novo assembly of eukaryotic genomes still presents significant hurdles, because of their large size and stretches of repetitive sequences. Filamentous fungi contain few repetitive regions in their 30–90 Mb genomes and thus are suitable candidates to test de novo genome assembly from short sequence reads. Here, we present a high-quality draft sequence of the Sordaria macrospora genome that was obtained by a combination of Illumina/Solexa and Roche/454 sequencing. Paired-end Solexa sequencing of genomic DNA to 85-fold coverage and an additional 10-fold coverage by single-end 454 sequencing resulted in ∼4 Gb of DNA sequence. Reads were assembled to a 40 Mb draft version (N50 of 117 kb) with the Velvet assembler. Comparative analysis with Neurospora genomes increased the N50 to 498 kb. The S. macrospora genome contains even fewer repeat regions than its closest sequenced relative, Neurospora crassa. Comparison with genomes of other fungi showed that S. macrospora, a model organism for morphogenesis and meiosis, harbors duplications of several genes involved in self/nonself-recognition. Furthermore, S. macrospora contains more polyketide biosynthesis genes than N. crassa. Phylogenetic analyses suggest that some of these genes may have been acquired by horizontal gene transfer from a distantly related ascomycete group. Our study shows that, for typical filamentous fungi, de novo assembly of genomes from short sequence reads alone is feasible, that a mixture of Solexa and 454 sequencing substantially improves the assembly, and that the resulting data can be used for comparative studies to address basic questions of fungal biology.

Species within the human pathogenic Cryptococcus species complex are major threats to public health, causing approximately 1 million annual infections globally. Cryptococcus amylolentus is the most closely known related species of the pathogenic Cryptococcus species complex, and it is non-pathogenic. Additionally, while pathogenic Cryptococcus species have bipolar mating systems with a single large mating type (MAT) locus that represents a derived state in Basidiomycetes, C. amylolentus has a tetrapolar mating system with 2 MAT loci (P/R and HD) located on different chromosomes. Thus, studying C. amylolentus will shed light on the transition from tetrapolar to bipolar mating systems in the pathogenic Cryptococcus species, as well as its possible link with the origin and evolution of pathogenesis. In this study, we sequenced, assembled, and annotated the genomes of 2 C. amylolentus isolates, CBS6039 and CBS6273, which are sexual and interfertile. Genome comparison between the 2 C. amylolentus isolates identified the boundaries and the complete gene contents of the P/R and HD MAT loci. Bioinformatic and chromatin immunoprecipitation sequencing (ChIP-seq) analyses revealed that, similar to those of the pathogenic Cryptococcus species, C. amylolentus has regional centromeres (CENs) that are enriched with species-specific transposable and repetitive DNA elements. Additionally, we found that while neither the P/R nor the HD locus is physically closely linked to its centromere in C. amylolentus, and the regions between the MAT loci and their respective centromeres show overall synteny between the 2 genomes, both MAT loci exhibit genetic linkage to their respective centromere during meiosis, suggesting the presence of recombinational suppressors and/or epistatic gene interactions in the MAT-CEN intervening regions. Furthermore, genomic comparisons between C. amylolentus and related pathogenic Cryptococcus species provide evidence that multiple chromosomal rearrangements mediated by intercentromeric recombination have occurred during descent of the 2 lineages from their common ancestor. Taken together, our findings support a model in which the evolution of the bipolar mating system was initiated by an ectopic recombination event mediated by similar repetitive centromeric DNA elements shared between chromosomes. This translocation brought the P/R and HD loci onto the same chromosome, and further chromosomal rearrangements then resulted in the 2 MAT loci becoming physically linked and eventually fusing to form the single contiguous MAT locus that is now extant in the pathogenic Cryptococcus species.

During sexual development, mycelial cells from most filamentous fungi differentiate into typical fruiting bodies. Here, we describe the isolation and characterization of the Sordaria macrospora developmental mutant per5, which exhibits a sterile phenotype with defects in fruiting body maturation. Cytological investigations revealed that the mutant strain forms only ascus precursors without any mature spores. Using an indexed cosmid library, we were able to complement the mutant to fertility by DNA-mediated transformation. A single cosmid clone, carrying a 3.5-kb region able to complement the mutant phenotype, has been identified. Sequencing of the 3.5-kb region revealed an open reading frame of 2.1 kb interrupted by a 66-bp intron. The predicted polypeptide (674 amino acids) shows significant homology to eukaryotic ATP citrate lyases (ACLs), with 62 to 65% amino acid identity, and the gene was named acl1. The molecular mass of the S. macrospora ACL1 polypeptide is 73 kDa, as was verified by Western blot analysis with a hemagglutinin (HA) epitope-tagged ACL1 polypeptide. Immunological in situ detection of the HA-tagged polypeptide demonstrated that ACL is located within the cytosol. Sequencing of the mutant acl1 gene revealed a 1-nucleotide transition within the coding region, resulting in an amino acid substitution within the predicted polypeptide. Further evidence that ACL1 is essential for fruiting body maturation comes from experiments in which truncated and mutated versions of the acl1 gene were used for transformation. None of these copies was able to reconstitute the fertile phenotype in transformed per5 recipient strains. ACLs are usually involved in the formation of cytosolic acetyl coenzyme A (acetyl-CoA), which is used for the biosynthesis of fatty acids and sterols. Protein extracts from the mutant strain showed a drastic reduction in enzymatic activity compared to values obtained from the wild-type strain. Investigation of the time course of ACL expression suggests that ACL is specifically induced at the beginning of the sexual cycle and produces acetyl-CoA, which most probably is a prerequisite for fruiting body formation during later stages of sexual development. We discuss the contribution of ACL activity to the life cycle of S. macrospora.The fruiting body maturation of filamentous ascomycetes is an attractive model system to study multicellular development in eukaryotes. It involves the formation of the outer structures of the fruiting body but also development of mature ascospores within the fruiting body itself (for reviews, see references 31 and 47). Ascus development starts with the formation of female gametangia called ascogonia. The ascogenous cells are enveloped by sterile hyphae to form fruiting body precursors. Subsequent tissue differentiation gives rise to an outer pigmented peridial tissue, and following caryogamy, inner ascus initials embedded in sterile paraphyses are formed. Mature fruiting bodies from most ascomycetes harbor 200 to 400 asci, which after meiosis and postmeiotic di...

BackgroundDuring sexual development, filamentous ascomycetes form complex, three-dimensional fruiting bodies for the protection and dispersal of sexual spores. Fruiting bodies contain a number of cell types not found in vegetative mycelium, and these morphological differences are thought to be mediated by changes in gene expression. However, little is known about the spatial distribution of gene expression in fungal development. Here, we used laser microdissection (LM) and RNA-seq to determine gene expression patterns in young fruiting bodies (protoperithecia) and non-reproductive mycelia of the ascomycete Sordaria macrospora.ResultsQuantitative analysis showed major differences in the gene expression patterns between protoperithecia and total mycelium. Among the genes strongly up-regulated in protoperithecia were the pheromone precursor genes ppg1 and ppg2. The up-regulation was confirmed by fluorescence microscopy of egfp expression under the control of ppg1 regulatory sequences. RNA-seq analysis of protoperithecia from the sterile mutant pro1 showed that many genes that are differentially regulated in these structures are under the genetic control of transcription factor PRO1.ConclusionsWe have generated transcriptional profiles of young fungal sexual structures using a combination of LM and RNA-seq. This allowed a high spatial resolution and sensitivity, and yielded a detailed picture of gene expression during development. Our data revealed significant differences in gene expression between protoperithecia and non-reproductive mycelia, and showed that the transcription factor PRO1 is involved in the regulation of many genes expressed specifically in sexual structures. The LM/RNA-seq approach will also be relevant to other eukaryotic systems in which multicellular development is investigated.

The advent of complex multicellularity (CM) was a pivotal event in the evolution of animals, plants and fungi. In the fungal Ascomycota, CM is based on hyphal filaments and arose in the Pezizomycotina. The genus Neolecta defines an enigma: phylogenetically placed in a related group containing mostly yeasts, Neolecta nevertheless possesses Pezizomycotina-like CM. Here we sequence the Neolecta irregularis genome and identify CM-associated functions by searching for genes conserved in Neolecta and the Pezizomycotina, which are absent or divergent in budding or fission yeasts. This group of 1,050 genes is enriched for functions related to diverse endomembrane systems and their organization. Remarkably, most show evidence for divergence in both yeasts. Using functional genomics, we identify new genes involved in fungal complexification. Together, these data show that rudimentary multicellularity is deeply rooted in the Ascomycota. Extensive parallel gene divergence during simplification and constraint leading to CM suggest a deterministic process where shared modes of cellular organization select for similarly configured organelle- and transport-related machineries.

The study of mutants to elucidate gene functions has a long and successful history; however, to discover causative mutations in mutants that were generated by random mutagenesis often takes years of laboratory work and requires previously generated genetic and/or physical markers, or resources like DNA libraries for complementation. Here, we present an alternative method to identify defective genes in developmental mutants of the filamentous fungus Sordaria macrospora through Illumina/Solexa whole-genome sequencing. We sequenced pooled DNA from progeny of crosses of three mutants and the wild type and were able to pinpoint the causative mutations in the mutant strains through bioinformatics analysis. One mutant is a spore color mutant, and the mutated gene encodes a melanin biosynthesis enzyme. The causative mutation is a G to A change in the first base of an intron, leading to a splice defect. The second mutant carries an allelic mutation in the pro41 gene encoding a protein essential for sexual development. In the mutant, we detected a complex pattern of deletion/rearrangements at the pro41 locus. In the third mutant, a point mutation in the stop codon of a transcription factor-encoding gene leads to the production of immature fruiting bodies. For all mutants, transformation with a wild type-copy of the affected gene restored the wild-type phenotype. Our data demonstrate that whole-genome sequencing of mutant strains is a rapid method to identify developmental genes in an organism that can be genetically crossed and where a reference genome sequence is available, even without prior mapping information.

The filamentous fungus Sordaria macrospora forms complex three-dimensional fruiting bodies that protect the developing ascospores and ensure their proper discharge. Several regulatory genes essential for fruiting body development were previously isolated by complementation of the sterile mutants pro1, pro11 and pro22. To establish the genetic relationships between these genes and to identify downstream targets, we have conducted cross-species microarray hybridizations using cDNA arrays derived from the closely related fungus Neurospora crassa and RNA probes prepared from wild-type S. macrospora and the three developmental mutants. Of the 1,420 genes which gave a signal with the probes from all the strains used, 172 (12%) were regulated differently in at least one of the three mutants compared to the wild type, and 17 (1.2%) were regulated differently in all three mutant strains. Microarray data were verified by Northern analysis or quantitative real time PCR. Among the genes that are up- or down-regulated in the mutant strains are genes encoding the pheromone precursors, enzymes involved in melanin biosynthesis and a lectin-like protein. Analysis of gene expression in double mutants revealed a complex network of interaction between the pro gene products.

NADPH oxidase (NOX)-derived reactive oxygen species (ROS) act as signaling determinants that induce different cellular processes. To characterize NOX function during fungal development, we utilized the genetically tractable ascomycete Sordaria macrospora. Genome sequencing of a sterile mutant led us to identify the NADPH oxidase encoding nox1 as a gene required for fruiting body formation, regular hyphal growth, and hyphal fusion. These phenotypes are shared by Δnor1, lacking the NOX regulator NOR1. Further phenotypic analyses revealed a high correlation between increased ROS production and hyphal fusion deficiencies in Δnox1 and other sterile mutants. A genome-wide transcriptional profiling analysis of mycelia and isolated protoperithecia from wild type and Δnox1 revealed that nox1 inactivation affects the expression of genes related to cytoskeleton remodeling, hyphal fusion, metabolism, and mitochondrial respiration. Genetic analysis of Δnox2, lacking the NADPH oxidase 2 gene, Δnor1, and transcription factor deletion mutant Δste12, revealed a strict melanin-dependent ascospore germination defect, indicating a common genetic pathway for these three genes. We report that gsa3, encoding a G-protein a-subunit, and sac1, encoding cAMP-generating adenylate cyclase, act in a separate pathway during the germination process. The finding that cAMP inhibits ascospore germination in a melanin-dependent manner supports a model in which cAMP inhibits NOX2 activity, thus suggesting a link between both pathways. Our results expand the current knowledge on the role of NOX enzymes in fungal development and provide a frame to define upstream and downstream components of the NOX signaling pathways in fungi. D URING sexual reproduction, filamentous fungi generate complex fruiting bodies that contain and protect meiosporangia. We used the ascomycetous model fungus Sordaria macrospora to identify genes directly involved in fruiting body development (Kück et al. 2009;Engh et al. 2010;Kück et al. 2009). Due to its homothallic life style, S. macrospora is able to complete the sexual life cycle without the mating of strains with opposite sex, and therefore, fruiting body-deficient mutants can be recognized directly without the need for crossing experiments. In earlier work, we generated sterile mutants showing a developmental block after formation of young fruiting bodies (protoperithecia), but being unable to generate mature perithecia, and referred to these mutants as pro. Recently, we have applied next-generation genome re-sequencing to identify the genes affected in some of these mutants . Based on this approach, we now have characterized mutant pro32 and show that it carries a mutation in the nox1 gene encoding NAPDH oxidase 1 (NOX1).NADPH oxidase (NOX) enzymes are transmembrane proteins that are highly conserved among eukaryotes and produce reactive oxygen species (ROS) through the oxidation of NADPH (Lambeth 2004;Kawahara and Lambeth 2007). ROS have long been recognized as damaging agents due to uncontrolled oxidizing re...