Rediscovery by Whole Genome Sequencing: Classical Mutations and Genome Polymorphisms in Neurospora crassa

Unveiling a classical mutant in the context of the GH3 β-glucosidase family in Neurospora crassa

Classical fungal mutant strains obtained by mutagenesis have helped to elucidate fundamental metabolic pathways in the past. In the filamentous fungus Neurospora crassa, the gluc-1 strain was isolated long ago and characterized by its low level of β-glucosidase activity, which is essential for the degradation of cellulose, the most abundant biopolymer on Earth and the main polymeric component of the plant cell wall. Based on genomic resequencing, we hypothesized that the causative mutation resides in the β-glucosidase gene gh3-3 (bgl6, NCU08755). In this work, growth patterns, enzymatic activities and sugar utilization rates were analyzed in several mutant and overexpression strains related to gluc-1 and gh3-3. In addition, different mutants affected in the degradation and transport of cellobiose were analyzed. While overexpression of gh3-3 led to the recovery of β-glucosidase activity in the gluc-1 mutant, as well as normal utilization of cellobiose, the full gene deletion strain Δgh3-3 was found to behave differently than gluc-1 with lower secreted β-glucosidase activity, indicating a dominant role of the amino acid substitution in the point mutated gh3-3 gene of gluc-1. Our results furthermore confirm that GH3-3 is the major extracellular β-glucosidase in N. crassa and demonstrate that the two cellodextrin transporters CDT-1 and CDT-2 are essential for growth on cellobiose when the three main N. crassa β-glucosidases are absent. Overall, these findings provide valuable insight into the mechanisms of cellulose utilization in filamentous fungi, being an essential step in the efficient production of biorefinable sugars from agricultural and forestry plant biomass.

Identification of a common secondary mutation in the Neurospora crassa knockout collection conferring a cell fusion-defective phenotype

Gene-deletion mutants represent a powerful tool to study gene function. The filamentous fungus Neurospora crassa is a well-established model organism, and features a comprehensive gene knockout strain collection. While these mutant strains have been used in numerous studies, resulting in the functional annotation of many Neurospora genes, direct confirmation of gene-phenotype relationships is often lacking, which is particularly relevant given the possibility of background mutations, sample contamination, and/or strain mislabeling. Indeed, spontaneous mutations resulting in phenotypes resembling many cell fusion mutants have long been known to occur at relatively high frequency in N. crassa, and these secondary mutations are common in the Neurospora deletion collection. The identity of these mutations, however, is largely unknown. Here, we report that the Δada-3 strain from the N. crassa knockout collection, which exhibits a cell fusion defect, harbors a secondary mutation responsible for this phenotype. Through whole-genome sequencing and genetic analyses, we found a ~30-Kb deletion in this strain affecting a known cell fusion-related gene, so/ham-1, and show that it is the absence of this gene-and not of ada-3-that underlies its cell fusion defect. We additionally found three other knockout strains harboring the same deletion, suggesting that this mutation may be common in the collection and could have impacted previous studies. Our findings provide a cautionary note and highlight the importance of proper functional validation of strains from mutant collections. We discuss our results in the context of the spread of cell fusion-defective cheater variants in N. crassa cultures. IMPORTANCE This study emphasizes the need for careful and detailed characterization of strains from mutant collections. Specifically, we found a common deletion in various strains from the Neurospora crassa gene knockout collection that results in a cell fusion-defective phenotype. This is noteworthy because this collection is known to contain background mutations-of a largely unclear nature-that produce cell fusion-defective phenotypes. Our results describe an example of such mutations, and highlight how this common genetic defect could have impacted previous studies that have used the affected strains. Furthermore, they provide a cautionary note about the use of Neurospora strains with similar phenotypes. Lastly, these findings offer additional details relevant to our understanding of the origin and spread of cell fusion-defective cheater variants in N. crassa cultures.

The Role of Chromatin and Transcriptional Control in the Formation of Sexual Fruiting Bodies in Fungi

Fungal fruiting bodies are complex, three-dimensional structures that arise from a less complex vegetative mycelium. Their formation requires the coordinated action of many genes and their gene products, and fruiting body formation is accompanied by major changes in the transcriptome. In recent years, numerous transcription factor genes as well as chromatin modifier genes that play a role in fruiting body morphogenesis were identified, and through research on several model organisms, the underlying regulatory networks that integrate chromatin structure, gene expression, and cell differentiation are becoming clearer. This review gives a summary of the current state of research on the role of transcriptional control and chromatin structure in fruiting body development. In the first part, insights from transcriptomics analyses are described, with a focus on comparative transcriptomics. In the second part, examples of more detailed functional characterizations of the role of chromatin modifiers and/or transcription factors in several model organisms (Neurospora crassa, Aspergillus nidulans, Sordaria macrospora, Coprinopsis cinerea, and Schizophyllum commune) that have led to a better understanding of regulatory networks at the level of chromatin structure and transcription are discussed.

Comparative Genomic Analysis Uncovers the Chloroplast Genome Variation and Phylogenetic Relationships of Camellia Species

Camellia is the largest genus in the family Theaceae. Due to phenotypic diversity, frequent hybridization, and polyploidization, an understanding of the phylogenetic relationships between Camellia species remains challenging. Comparative chloroplast (cp) genomics provides an informative resource for phylogenetic analyses of Camellia. In this study, 12 chloroplast genome sequences from nine Camellia species were determined using Illumina sequencing technology via de novo assembly. The cp genome sizes ranged from 156,545 to 157,021 bp and were organized into quadripartite regions with the typical angiosperm cp genomes. Each genome harbored 87 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes in the same order and orientation. Differences in long and short sequence repeats, SNPs, and InDels were detected across the 12 cp genomes. Combining with the complete cp sequences of seven other species in the genus Camellia, a total of nine intergenic sequence divergent hotspots and 14 protein-coding genes with high sequence polymorphism were identified. These hotspots, especially the InDel (~400 bp) located in atpH-atpI region, had sufficient potential to be used as barcode markers for further phylogenetic analysis and species identification. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These cp genomes could facilitate the development of new molecular markers, accurate species identification, and investigations of the phylogenomic relationships of the genus Camellia.

Complete Genome Sequences and Genome-Wide Characterization of Trichoderma Biocontrol Agents Provide New Insights into their Evolution and Variation in Genome Organization, Sexual Development, and Fungal-Plant Interactions

Trichoderma spp. represent one of the most important fungal genera to mankind and in natural environments. The genus harbors prolific producers of wood-decaying enzymes, biocontrol agents against plant pathogens, plant-growth-promoting biofertilizers, as well as model organisms for studying fungal-plant-plant pathogen interactions. Pursuing highly accurate, contiguous, and chromosome-level reference genomes has become a primary goal of fungal research communities. Here, we report the chromosome-level genomic sequences and whole-genome annotation data sets of four strains used as biocontrol agents or biofertilizers (Trichoderma virens Gv29-8, Trichoderma virens FT-333, Trichoderma asperellum FT-101, and Trichoderma atroviride P1). Our results provide comprehensive categorization, correct positioning, and evolutionary detail of both nuclear and mitochondrial genomes, including telomeres, AT-rich blocks, centromeres, transposons, mating-type loci, nuclear-encoded mitochondrial sequences, as well as many new secondary metabolic and carbohydrate-active enzyme gene clusters. We have also identified evolutionarily conserved core genes contributing to plant-fungal interactions, as well as variations potentially linked to key behavioral traits such as sex, genome defense, secondary metabolism, and mycoparasitism. The genomic resources we provide herein significantly extend our knowledge not only of this economically important fungal genus, but also fungal evolution and basic biology in general. IMPORTANCE Telomere-to-telomere and gapless reference genome assemblies are necessary to ensure that all genomic variants are studied and discovered, including centromeres, telomeres, AT-rich blocks, mating type loci, biosynthetic, and metabolic gene clusters. Here, we applied long-range sequencing technologies to determine the near-completed genome sequences of four widely used biocontrol agents or biofertilizers: Trichoderma virens Gv29-8 and FT-333, Trichoderma asperellum FT-101, and Trichoderma atroviride P1. Like those of three Trichoderma reesei wild isolates [QM6a, CBS999.97(MAT1-1) and CBS999.97(MAT1-2)] we reported previously, these four biocontrol agent genomes each contain seven nuclear chromosomes and a circular mitochondrial genome. Substantial intraspecies and intragenus diversities are also discovered, including single nucleotide polymorphisms, chromosome shuffling, as well as genomic relics derived from historical transposition events and repeat-induced point (RIP) mutations.

Important role of melanin for fertility in the fungus Podospora anserina

Melanins are pigments used by fungi to withstand various stresses and to strengthen vegetative and reproductive structures. In Sordariales fungi, their biosynthesis starts with a condensation step catalyzed by an evolutionary-conserved polyketide synthase. Here we show that complete inactivation of this enzyme in the model ascomycete Podospora anserina through targeted deletion of the PaPks1 gene results in reduced female fertility, in contrast to a previously analyzed nonsense mutation in the same gene that retains full fertility. We also show the utility of PaPks1 mutants for detecting rare genetic events in P. anserina, such as parasexuality and possible fertilization and/or apomixis of nuclei devoid of mating-type gene.

The F-box protein gene exo-1 is a target for reverse engineering enzyme hypersecretion in filamentous fungi

Carbohydrate active enzymes (CAZymes) are vital for the lignocellulose-based biorefinery. The development of hypersecreting fungal protein production hosts is therefore a major aim for both academia and industry. However, despite advances in our understanding of their regulation, the number of promising candidate genes for targeted strain engineering remains limited. Here, we resequenced the genome of the classical hypersecreting Neurospora crassa mutant exo-1 and identified the causative point of mutation to reside in the F-box protein-encoding gene, NCU09899. The corresponding deletion strain displayed amylase and invertase activities exceeding those of the carbon catabolite derepressed strain Δcre-1, while glucose repression was still mostly functional in Δexo-1 Surprisingly, RNA sequencing revealed that while plant cell wall degradation genes are broadly misexpressed in Δexo-1, only a small fraction of CAZyme genes and sugar transporters are up-regulated, indicating that EXO-1 affects specific regulatory factors. Aiming to elucidate the underlying mechanism of enzyme hypersecretion, we found the high secretion of amylases and invertase in Δexo-1 to be completely dependent on the transcriptional regulator COL-26. Furthermore, misregulation of COL-26, CRE-1, and cellular carbon and nitrogen metabolism was confirmed by proteomics. Finally, we successfully transferred the hypersecretion trait of the exo-1 disruption by reverse engineering into the industrially deployed fungus Myceliophthora thermophila using CRISPR-Cas9. Our identification of an important F-box protein demonstrates the strength of classical mutants combined with next-generation sequencing to uncover unanticipated candidates for engineering. These data contribute to a more complete understanding of CAZyme regulation and will facilitate targeted engineering of hypersecretion in further organisms of interest.

A global search for novel transcription factors impacting the Neurospora crassa circadian clock

Eukaryotic circadian oscillators share a common circuit architecture, a negative feedback loop in which a positive element activates the transcription of a negative one that then represses the action of the former, inhibiting its own expression. While studies in mammals and insects have revealed additional transcriptional inputs modulating the expression of core clock components, this has been less characterized in the model Neurospora crassa, where the participation of other transcriptional components impacting circadian clock dynamics remains rather unexplored. Thus, we sought to identify additional transcriptional regulators modulating the N. crassa clock, following a reverse genetic screen based on luminescent circadian reporters and a collection of transcription factors (TFs) knockouts, successfully covering close to 60% of them. Besides the canonical core clock components WC-1 and -2, none of the tested transcriptional regulators proved to be essential for rhythmicity. Nevertheless, we identified a set of 23 TFs that when absent lead to discrete, but significant, changes in circadian period. While the current level of analysis does not provide mechanistic information about how these new players modulate circadian parameters, the results of this screen reveal that an important number of light and clock-regulated TFs, involved in a plethora of processes, are capable of modulating the clockworks. This partial reverse genetic clock screen also exemplifies how the N. crassa knockout collection continues to serve as an expedite platform to address broad biological questions.

Repeat-induced point mutation in Neurospora crassa causes the highest known mutation rate and mutational burden of any cellular life

BACKGROUND

Repeat-induced point (RIP) mutation in Neurospora crassa degrades transposable elements by targeting repeats with C→T mutations. Whether RIP affects core genomic sequence in important ways is unknown.

RESULTS

By parent-offspring whole genome sequencing, we estimate a mutation rate (3.38 ×  10-6 per bp per generation) that is two orders of magnitude higher than reported for any non-viral organism, with 93-98% of mutations being RIP-associated. RIP mutations are, however, relatively rare in coding sequence, in part because RIP preferentially attacks GC-poor long duplicates that interact in three dimensional space, while coding sequence duplicates are rare, GC-rich, short, and tend not to interact. Despite this, with over 5 coding sequence mutations per genome per generation, the mutational burden is an order of magnitude higher than the previously highest observed. Unexpectedly, the majority of these coding sequence mutations appear not to be the direct product of RIP being mostly in non-duplicate sequence and predominantly not C→T mutations. Nonetheless, RIP-deficient strains have over an order of magnitude fewer coding sequence mutations outside of duplicated domains than RIP-proficient strains.

CONCLUSIONS

Neurospora crassa has the highest mutation rate and mutational burden of any non-viral life. While the high rate is largely due to the action of RIP, the mutational burden appears to be RIP-associated but not directly caused by RIP.

Whole Genome Sequencing Reveals Major Deletions in the Genome of M7, a Gamma Ray-Induced Mutant of Trichoderma virens That Is Repressed in Conidiation, Secondary Metabolism, and Mycoparasitism

Trichoderma virens is a commercial biofungicide used in agriculture. We have earlier isolated a mutant of T. virens using gamma ray-induced mutagenesis. This mutant, designated as M7, is defective in morphogenesis, secondary metabolism, and mycoparasitism. The mutant does not produce conidia, and the colony is hydrophilic. M7 cannot utilize cellulose and chitin as a sole carbon source and is unable to parasitize the plant pathogens Rhizoctonia solani and Pythium aphanidermatum in confrontation assay. Several volatile (germacrenes, beta-caryophyllene, alloaromadendrene, gamma-muurolene) and non-volatile (viridin, viridiol, gliovirin, heptelidic acid) metabolites are not detected in M7. In transcriptome analysis, many genes related to secondary metabolism, carbohydrate metabolism, hydrophobicity, and transportation, among others, were found to be downregulated in the mutant. Using whole genome sequencing, we identified five deletions in the mutant genome, totaling about 250 kb (encompassing 71 predicted ORFs), which was confirmed by PCR. This study provides novel insight into genetics of morphogenesis, secondary metabolism, and mycoparasitism and eventually could lead to the identification of novel regulators of beneficial traits in plant beneficial fungi Trichoderma spp. We also suggest that this mutant can be developed as a microbial cell factory for the production of secondary metabolites and proteins.

The Fungal Cell Death Regulator czt-1 Is Allelic to acr-3

Fungal infections have far-reaching implications that range from severe human disease to a panoply of disruptive agricultural and ecological effects, making it imperative to identify and understand the molecular pathways governing the response to antifungal compounds. In this context, CZT-1 (cell death-activated zinc cluster transcription factor) functions as a master regulator of cell death and drug susceptibility in Neurospora crassa. Here we provide evidence indicating that czt-1 is allelic to acr-3, a previously described locus that we now found to harbor a point mutation in its coding sequence. This nonsynonymous amino acid substitution in a low complexity region of CZT-1/ACR-3 caused a robust gain-of-function that led to reduced sensitivity to acriflavine and staurosporine, and increased expression of the drug efflux pump abc-3. Thus, accumulating evidence shows that CZT-1 is an important broad regulator of the cellular response to various antifungal compounds that appear to share common molecular targets.

The Genetics and Biochemistry of Cell Wall Structure and Synthesis in Neurospora crassa, a Model Filamentous Fungus

This review discusses the wealth of information available for the N. crassa cell wall. The basic organization and structure of the cell wall is presented and how the wall changes during the N. crassa life cycle is discussed. Over forty cell wall glycoproteins have been identified by proteomic analyses. Genetic and biochemical studies have identified many of the key enzymes needed for cell wall biogenesis, and the roles these enzymes play in cell wall biogenesis are discussed. The review includes a discussion of how the major cell wall components (chitin, β-1,3-glucan, mixed β-1,3-/ β-1,4- glucans, glycoproteins, and melanin) are synthesized and incorporated into the cell wall. We present a four-step model for how cell wall glycoproteins are covalently incorporated into the cell wall. In N. crassa, the covalent incorporation of cell wall glycoproteins into the wall occurs through a glycosidic linkage between lichenin (a mixed β-1,3-/β-1,4- glucan) and a "processed" galactomannan that has been attached to the glycoprotein N-linked oligosaccharides. The first step is the addition of the galactomannan to the N-linked oligosaccharide. Mutants affected in galactomannan formation are unable to incorporate glycoproteins into their cell walls. The second step is carried out by the enzymes from the GH76 family of α-1,6-mannanases, which cleave the galactomannan to generate a processed galactomannan. The model suggests that the third and fourth steps are carried out by members of the GH72 family of glucanosyltransferases. In the third step the glucanosyltransferases cleave lichenin and generate enzyme/substrate intermediates in which the lichenin is covalently attached to the active site of the glucanosyltransferases. In the final step, the glucanosyltransferases attach the lichenin onto the processed galactomannans, which creates new glycosidic bonds and effectively incorporates the glycoproteins into the cross-linked cell wall glucan/chitin matrix.

A RID-like putative cytosine methyltransferase homologue controls sexual development in the fungus Podospora anserina

DNA methyltransferases are ubiquitous enzymes conserved in bacteria, plants and opisthokonta. These enzymes, which methylate cytosines, are involved in numerous biological processes, notably development. In mammals and higher plants, methylation patterns established and maintained by the cytosine DNA methyltransferases (DMTs) are essential to zygotic development. In fungi, some members of an extensively conserved fungal-specific DNA methyltransferase class are both mediators of the Repeat Induced Point mutation (RIP) genome defense system and key players of sexual reproduction. Yet, no DNA methyltransferase activity of these purified RID (RIP deficient) proteins could be detected in vitro. These observations led us to explore how RID-like DNA methyltransferase encoding genes would play a role during sexual development of fungi showing very little genomic DNA methylation, if any. To do so, we used the model ascomycete fungus Podospora anserina. We identified the PaRid gene, encoding a RID-like DNA methyltransferase and constructed knocked-out ΔPaRid defective mutants. Crosses involving P. anserina ΔPaRid mutants are sterile. Our results show that, although gametes are readily formed and fertilization occurs in a ΔPaRid background, sexual development is blocked just before the individualization of the dikaryotic cells leading to meiocytes. Complementation of ΔPaRid mutants with ectopic alleles of PaRid, including GFP-tagged, point-mutated and chimeric alleles, demonstrated that the catalytic motif of the putative PaRid methyltransferase is essential to ensure proper sexual development and that the expression of PaRid is spatially and temporally restricted. A transcriptomic analysis performed on mutant crosses revealed an overlap of the PaRid-controlled genetic network with the well-known mating-types gene developmental pathway common to an important group of fungi, the Pezizomycotina.

Characterization of the Neurospora crassa DHN melanin biosynthetic pathway in developing ascospores and peridium cells

Neurospora crassa contains all four enzymes for the synthesis of DHN (dihydroxynaphthalene), the substrate for melanin formation. We show that the DHN melanin pathway functions during N. crassa female development to generate melanized peridium and ascospore cell walls. N. crassa contains one polyketide synthase (PER-1), two polyketide hydrolases (PKH-1 and PKH-2), two THN (tetrahydroxynaphthalene) reductases (PKR-1 and PKR-2), and one scytalone dehydratase (SCY-1). We show that the PER-1, PKH-1, PKR-1 and SCY-1 are required for ascospoer melanization. We also identified the laccase that functions in the conversion of DHN into melanin via a free radical oxidative polymerization reaction, and have named the gene lacm-1 (laccase for melanin formation-1). In maturing perithecia, we show that LACM-1 is localized to the peridium cell wall space while the DHN pathway enzymes are localized to intracellular vesicles. We present a model for melanin formation in which melanin is formed within the cell wall space and the cell wall structure is similar to "reinforced concrete" with the cell wall glucan, chitin, and glycoproteins encased within the melanin polymer. This arrangement provides for a very strong and resilient cell wall and protects the glucan/chitin/glycoprotein matrix from digestion from enzymes and damage from free radicals.

Parasexual Crossings for Bulk Segregant Analysis in Aspergillus niger to Facilitate Mutant Identification Via Whole Genome Sequencing

The industrially important fungus Aspergillus niger is known to reproduce only asexually. The parasexual cycle of fungi can be used for crossing two different strains to produce segregants or progeny with combined mutations even in fungi without a known sexual cycle. In A. niger, the parasexual cycle has been extensively used to establish linkage groups and to generate genetic maps. With the advent of whole genome sequencing, the parasexual cycle has received renewed attention as a method to create segregants for bulk segregant analysis. Bulk segregant analysis is a genetic technique used to link and ultimately identify the mutation associated with a particular phenotype. In this chapter we describe the procedure for setting up parasexual crossings in A. niger. The segregants obtained with this method can be used in combination with next-generation sequencing to map mutations in the organism.

Functional Profiling of Transcription Factor Genes in Neurospora crassa

Regulation of gene expression by DNA-binding transcription factors is essential for proper control of growth and development in all organisms. In this study, we annotate and characterize growth and developmental phenotypes for transcription factor genes in the model filamentous fungus Neurospora crassa. We identified 312 transcription factor genes, corresponding to 3.2% of the protein coding genes in the genome. The largest class was the fungal-specific Zn₂Cys₆ (C6) binuclear cluster, with 135 members, followed by the highly conserved C2H2 zinc finger group, with 61 genes. Viable knockout mutants were produced for 273 genes, and complete growth and developmental phenotypic data are available for 242 strains, with 64% possessing at least one defect. The most prominent defect observed was in growth of basal hyphae (43% of mutants analyzed), followed by asexual sporulation (38%), and the various stages of sexual development (19%). Two growth or developmental defects were observed for 21% of the mutants, while 8% were defective in all three major phenotypes tested. Analysis of available mRNA expression data for a time course of sexual development revealed mutants with sexual phenotypes that correlate with transcription factor transcript abundance in wild type. Inspection of this data also implicated cryptic roles in sexual development for several cotranscribed transcription factor genes that do not produce a phenotype when mutated.

Revitalization of a Forward Genetic Screen Identifies Three New Regulators of Fungal Secondary Metabolism in the Genus Aspergillus

The study of aflatoxin in Aspergillus spp. has garnered the attention of many researchers due to aflatoxin's carcinogenic properties and frequency as a food and feed contaminant. Significant progress has been made by utilizing the model organism Aspergillus nidulans to characterize the regulation of sterigmatocystin (ST), the penultimate precursor of aflatoxin. A previous forward genetic screen identified 23 A. nidulans mutants involved in regulating ST production. Six mutants were characterized from this screen using classical mapping (five mutations in mcsA) and complementation with a cosmid library (one mutation in laeA). The remaining mutants were backcrossed and sequenced using Illumina and Ion Torrent sequencing platforms. All but one mutant contained one or more sequence variants in predicted open reading frames. Deletion of these genes resulted in identification of mutant alleles responsible for the loss of ST production in 12 of the 17 remaining mutants. Eight of these mutations were in genes already known to affect ST synthesis (laeA, mcsA, fluG, and stcA), while the remaining four mutations (in laeB, sntB, and hamI) were in previously uncharacterized genes not known to be involved in ST production. Deletion of laeB, sntB, and hamI in A. flavus results in loss of aflatoxin production, confirming that these regulators are conserved in the aflatoxigenic aspergilli. This report highlights the multifaceted regulatory mechanisms governing secondary metabolism in Aspergillus Additionally, these data contribute to the increasing number of studies showing that forward genetic screens of fungi coupled with whole-genome resequencing is a robust and cost-effective technique.IMPORTANCE In a postgenomic world, reverse genetic approaches have displaced their forward genetic counterparts. The techniques used in forward genetics to identify loci of interest were typically very cumbersome and time-consuming, relying on Mendelian traits in model organisms. The current work was pursued not only to identify alleles involved in regulation of secondary metabolism but also to demonstrate a return to forward genetics to track phenotypes and to discover genetic pathways that could not be predicted through a reverse genetics approach. While identification of mutant alleles from whole-genome sequencing has been done before, here we illustrate the possibility of coupling this strategy with a genetic screen to identify multiple alleles of interest. Sequencing of classically derived mutants revealed several uncharacterized genes, which represent novel pathways to regulate and control the biosynthesis of sterigmatocystin and of aflatoxin, a societally and medically important mycotoxin.

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.

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.

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.

Current challenges of research on filamentous fungi in relation to human welfare and a sustainable bio-economy: a white paper

The EUROFUNG network is a virtual centre of multidisciplinary expertise in the field of fungal biotechnology. The first academic-industry Think Tank was hosted by EUROFUNG to summarise the state of the art and future challenges in fungal biology and biotechnology in the coming decade. Currently, fungal cell factories are important for bulk manufacturing of organic acids, proteins, enzymes, secondary metabolites and active pharmaceutical ingredients in white and red biotechnology. In contrast, fungal pathogens of humans kill more people than malaria or tuberculosis. Fungi are significantly impacting on global food security, damaging global crop production, causing disease in domesticated animals, and spoiling an estimated 10 % of harvested crops. A number of challenges now need to be addressed to improve our strategies to control fungal pathogenicity and to optimise the use of fungi as sources for novel compounds and as cell factories for large scale manufacture of bio-based products. This white paper reports on the discussions of the Think Tank meeting and the suggestions made for moving fungal bio(techno)logy forward.

period-1 encodes an ATP-dependent RNA helicase that influences nutritional compensation of the Neurospora circadian clock

Mutants in the period-1 (prd-1) gene, characterized by a recessive allele, display a reduced growth rate and period lengthening of the developmental cycle controlled by the circadian clock. We refined the genetic location of prd-1 and used whole genome sequencing to find the mutation defining it, confirming the identity of prd-1 by rescuing the mutant circadian phenotype via transformation. PRD-1 is an RNA helicase whose orthologs, DDX5 [DEAD (Asp-Glu-Ala-Asp) Box Helicase 5] and DDX17 in humans and DBP2 (Dead Box Protein 2) in yeast, are implicated in various processes, including transcriptional regulation, elongation, and termination, ribosome biogenesis, and mRNA decay. Although prd-1 mutants display a long period (∼25 h) circadian developmental cycle, they interestingly display a WT period when the core circadian oscillator is tracked using a frq-luciferase transcriptional fusion under conditions of limiting nutritional carbon; the core oscillator in the prd-1 mutant strain runs with a long period under glucose-sufficient conditions. Thus, PRD-1 clearly impacts the circadian oscillator and is not only part of a metabolic oscillator ancillary to the core clock. PRD-1 is an essential protein, and its expression is neither light-regulated nor clock-regulated. However, it is transiently induced by glucose; in the presence of sufficient glucose, PRD-1 is in the nucleus until glucose runs out, which elicits its disappearance from the nucleus. Because circadian period length is carbon concentration-dependent, prd-1 may be formally viewed as a clock mutant with defective nutritional compensation of circadian period length.

The United States Culture Collection Network (USCCN): Enhancing Microbial Genomics Research through Living Microbe Culture Collections

The mission of the United States Culture Collection Network (USCCN; http://usccn.org) is "to facilitate the safe and responsible utilization of microbial resources for research, education, industry, medicine, and agriculture for the betterment of human kind." Microbial culture collections are a key component of life science research, biotechnology, and emerging global biobased economies. Representatives and users of several microbial culture collections from the United States and Europe gathered at the University of California, Davis, to discuss how collections of microorganisms can better serve users and stakeholders and to showcase existing resources available in public culture collections.

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.

Genome sequencing of the Trichoderma reesei QM9136 mutant identifies a truncation of the transcriptional regulator XYR1 as the cause for its cellulase-negative phenotype

Background

Trichoderma reesei is the main industrial source of cellulases and hemicellulases required for the hydrolysis of biomass to simple sugars, which can then be used in the production of biofuels and biorefineries. The highly productive strains in use today were generated by classical mutagenesis. As byproducts of this procedure, mutants were generated that turned out to be unable to produce cellulases. In order to identify the mutations responsible for this inability, we sequenced the genome of one of these strains, QM9136, and compared it to that of its progenitor T. reesei QM6a.

Results

In QM9136, we detected a surprisingly low number of mutagenic events in the promoter and coding regions of genes, i.e. only eight indels and six single nucleotide variants. One of these indels led to a frame-shift in the Zn₂Cys₆ transcription factor XYR1, the general regulator of cellulase and xylanase expression, and resulted in its C-terminal truncation by 140 amino acids. Retransformation of strain QM9136 with the wild-type xyr1 allele fully recovered the ability to produce cellulases, and is thus the reason for the cellulase-negative phenotype. Introduction of an engineered xyr1 allele containing the truncating point mutation into the moderate producer T. reesei QM9414 rendered this strain also cellulase-negative. The correspondingly truncated XYR1 protein was still able to enter the nucleus, but failed to be expressed over the basal constitutive level.

Conclusion

The missing 140 C-terminal amino acids of XYR1 are therefore responsible for its previously observed auto-regulation which is essential for cellulases to be expressed. Our data present a working example of the use of genome sequencing leading to a functional explanation of the QM9136 cellulase-negative phenotype.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1526-0) contains supplementary material, which is available to authorized users.

Whole-genome analysis of Fusarium graminearum insertional mutants identifies virulence associated genes and unmasks untagged chromosomal deletions

Background

Identifying pathogen virulence genes required to cause disease is crucial to understand the mechanisms underlying the pathogenic process. Plasmid insertion mutagenesis of fungal protoplasts is frequently used for this purpose in filamentous ascomycetes. Post transformation, the mutant population is screened for loss of virulence to a specific plant or animal host. Identifying the insertion event has previously met with varying degrees of success, from a cleanly disrupted gene with minimal deletion of nucleotides at the insertion point to multiple-copy insertion events and large deletions of chromosomal regions. Currently, extensive mutant collections exist in laboratories globally where it was hitherto impossible to identify all the affected genes.

Results

We used a whole-genome sequencing (WGS) approach using Illumina HiSeq 2000 technology to investigate DNA tag insertion points and chromosomal deletion events in mutagenised, reduced virulence F. graminearum isolates identified in disease tests on wheat (Triticum aestivum). We developed the FindInsertSeq workflow to localise the DNA tag insertions to the nucleotide level. The workflow was tested using four mutants showing evidence of single and multi-copy insertions in DNA blot analysis. FindInsertSeq was able to identify both single and multi-copy concatenation insertion sites. By comparing sequencing coverage, unexpected molecular recombination events such as large tagged and untagged chromosomal deletions, and DNA amplification were observed in three of the analysed mutants. A random data sampling approach revealed the minimum genome coverage required to survey the F. graminearum genome for alterations.

Conclusions

This study demonstrates that whole-genome re-sequencing to 22x fold genome coverage is an efficient tool to characterise single and multi-copy insertion mutants in the filamentous ascomycete Fusarium graminearum. In some cases insertion events are accompanied with large untagged chromosomal deletions while in other cases a straight-forward insertion event could be confirmed. The FindInsertSeq analysis workflow presented in this study enables researchers to efficiently characterise insertion and deletion mutants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1412-9) contains supplementary material, which is available to authorized users.

Draft Genome Sequence of Neurospora crassa Strain FGSC 73

We report the elucidation of the complete genome of the Neurospora crassa (Shear and Dodge) strain FGSC 73, a mat-a, trp-3 mutant strain. The genome sequence around the idiotypic mating type locus represents the only publicly available sequence for a mat-a strain. 40.42 Megabases are assembled into 358 scaffolds carrying 11,978 gene models.

Characterization of the mutagenic spectrum of 4-nitroquinoline 1-oxide (4-NQO) in Aspergillus nidulans by whole genome sequencing

4-Nitroquinoline 1-oxide (4-NQO) is a highly carcinogenic chemical that induces mutations in bacteria, fungi, and animals through the formation of bulky purine adducts. 4-NQO has been used as a mutagen for genetic screens and in both the study of DNA damage and DNA repair. In the model eukaryote Aspergillus nidulans, 4-NQO-based genetic screens have been used to study diverse processes, including gene regulation, mitosis, metabolism, organelle transport, and septation. Early work during the 1970s using bacterial and yeast mutation tester strains concluded that 4-NQO was a guanine-specific mutagen. However, these strains were limited in their ability to determine full mutagenic potential, as they could not identify mutations at multiple sites, unlinked suppressor mutations, or G:C to C:G transversions. We have now used a whole genome resequencing approach with mutant strains generated from two independent genetic screens to determine the full mutagenic spectrum of 4-NQO in A. nidulans. Analysis of 3994 mutations from 38 mutant strains reveals that 4-NQO induces substitutions in both guanine and adenine residues, although with a 19-fold preference for guanine. We found no association between mutation load and mutagen dose and observed no sequence bias in the residues flanking the mutated purine base. The mutations were distributed randomly throughout most of the genome. Our data provide new evidence that 4-NQO can potentially target all base pairs. Furthermore, we predict that current practices for 4-NQO-induced mutagenesis are sufficient to reach gene saturation for genetic screens with feasible identification of causative mutations via whole genome resequencing.

Neurospora crassa female development requires the PACC and other signal transduction pathways, transcription factors, chromatin remodeling, cell-to-cell fusion, and autophagy

Using a screening protocol we have identified 68 genes that are required for female development in the filamentous fungus Neurospora crassa. We find that we can divide these genes into five general groups: 1) Genes encoding components of the PACC signal transduction pathway, 2) Other signal transduction pathway genes, including genes from the three N. crassa MAP kinase pathways, 3) Transcriptional factor genes, 4) Autophagy genes, and 5) Other miscellaneous genes. Complementation and RIP studies verified that these genes are needed for the formation of the female mating structure, the protoperithecium, and for the maturation of a fertilized protoperithecium into a perithecium. Perithecia grafting experiments demonstrate that the autophagy genes and the cell-to-cell fusion genes (the MAK-1 and MAK-2 pathway genes) are needed for the mobilization and movement of nutrients from an established vegetative hyphal network into the developing protoperithecium. Deletion mutants for the PACC pathway genes palA, palB, palC, palF, palH, and pacC were found to be defective in two aspects of female development. First, they were unable to initiate female development on synthetic crossing medium. However, they could form protoperithecia when grown on cellophane, on corn meal agar, or in response to the presence of nearby perithecia. Second, fertilized perithecia from PACC pathway mutants were unable to produce asci and complete female development. Protein localization experiments with a GFP-tagged PALA construct showed that PALA was localized in a peripheral punctate pattern, consistent with a signaling center associated with the ESCRT complex. The N. crassa PACC signal transduction pathway appears to be similar to the PacC/Rim101 pathway previously characterized in Aspergillus nidulans and Saccharomyces cerevisiae. In N. crassa the pathway plays a key role in regulating female development.

Using next-generation sequencing to isolate mutant genes from forward genetic screens

The long-lasting success of forward genetic screens relies on the simple molecular basis of the characterized phenotypes, which are typically caused by mutations in single genes. Mapping the location of causal mutations using genetic crosses has traditionally been a complex, multistep procedure, but next-generation sequencing now allows the rapid identification of causal mutations at single-nucleotide resolution even in complex genetic backgrounds. Recent advances of this mapping-by-sequencing approach include methods that are independent of reference genome sequences, genetic crosses and any kind of linkage information, which make forward genetics amenable for species that have not been considered for forward genetic screens so far.

Thirteen Camellia chloroplast genome sequences determined by high-throughput sequencing: genome structure and phylogenetic relationships

BACKGROUND

Camellia is an economically and phylogenetically important genus in the family Theaceae. Owing to numerous hybridization and polyploidization, it is taxonomically and phylogenetically ranked as one of the most challengingly difficult taxa in plants. Sequence comparisons of chloroplast (cp) genomes are of great interest to provide a robust evidence for taxonomic studies, species identification and understanding mechanisms that underlie the evolution of the Camellia species.

RESULTS

The eight complete cp genomes and five draft cp genome sequences of Camellia species were determined using Illumina sequencing technology via a combined strategy of de novo and reference-guided assembly. The Camellia cp genomes exhibited typical circular structure that was rather conserved in genomic structure and the synteny of gene order. Differences of repeat sequences, simple sequence repeats, indels and substitutions were further examined among five complete cp genomes, representing a wide phylogenetic diversity in the genus. A total of fifteen molecular markers were identified with more than 1.5% sequence divergence that may be useful for further phylogenetic analysis and species identification of Camellia. Our results showed that, rather than functional constrains, it is the regional constraints that strongly affect sequence evolution of the cp genomes. In a substantial improvement over prior studies, evolutionary relationships of the section Thea were determined on basis of phylogenomic analyses of cp genome sequences.

CONCLUSIONS

Despite a high degree of conservation between the Camellia cp genomes, sequence variation among species could still be detected, representing a wide phylogenetic diversity in the genus. Furthermore, phylogenomic analysis was conducted using 18 complete cp genomes and 5 draft cp genome sequences of Camellia species. Our results support Chang's taxonomical treatment that C. pubicosta may be classified into sect. Thea, and indicate that taxonomical value of the number of ovaries should be reconsidered when classifying the Camellia species. The availability of these cp genomes provides valuable genetic information for accurately identifying species, clarifying taxonomy and reconstructing the phylogeny of the genus Camellia.

Global analysis of serine/threonine and tyrosine protein phosphatase catalytic subunit genes in Neurospora crassa reveals interplay between phosphatases and the p38 mitogen-activated protein kinase

Protein phosphatases are integral components of the cellular signaling machinery in eukaryotes, regulating diverse aspects of growth and development. The genome of the filamentous fungus and model organism Neurospora crassa encodes catalytic subunits for 30 protein phosphatase genes. In this study, we have characterized 24 viable N. crassa phosphatase catalytic subunit knockout mutants for phenotypes during growth, asexual development, and sexual development. We found that 91% of the mutants had defects in at least one of these traits, whereas 29% possessed phenotypes in all three. Chemical sensitivity screens were conducted to reveal additional phenotypes for the mutants. This resulted in the identification of at least one chemical sensitivity phenotype for 17 phosphatase knockout mutants, including novel chemical sensitivities for two phosphatase mutants lacking a growth or developmental phenotype. Hence, chemical sensitivity or growth/developmental phenotype was observed for all 24 viable mutants. We investigated p38 mitogen-activated protein kinase (MAPK) phosphorylation profiles in the phosphatase mutants and identified nine potential candidates for regulators of the p38 MAPK. We demonstrated that the PP2C class phosphatase pph-8 (NCU04600) is an important regulator of female sexual development in N. crassa. In addition, we showed that the Δcsp-6 (ΔNCU08380) mutant exhibits a phenotype similar to the previously identified conidial separation mutants, Δcsp-1 and Δcsp-2, that lack transcription factors important for regulation of conidiation and the circadian clock.

Investment into the future of microbial resources: culture collection funding models and BRC business plans for biological resource centres

Through their long history of public service, diverse microbial Biological Resource Centres (mBRCs) have made myriad contributions to society and science. They have enabled the maintenance of specimens isolated before antibiotics, made available strains showing the development and change of pathogenicity toward animals, humans and plants, and have maintained and provided reference strains to ensure quality and reproducibility of science. However, this has not been achieved without considerable financial commitment. Different collections have unique histories and their support is often tied to their origins. However many collections have grown to serve large constituencies and need to develop novel funding mechanisms. Moreover, several international initiatives have described mBRCs as a factor in economic development and have led to the increased professionalism among mBRCs.

The filamentous fungus Sordaria macrospora as a genetic model to study fruiting body development

Filamentous fungi are excellent experimental systems due to their short life cycles as well as easy and safe manipulation in the laboratory. They form three-dimensional structures with numerous different cell types and have a long tradition as genetic model organisms used to unravel basic mechanisms underlying eukaryotic cell differentiation. The filamentous ascomycete Sordaria macrospora is a model system for sexual fruiting body (perithecia) formation. S. macrospora is homothallic, i.e., self-fertile, easily genetically tractable, and well suited for large-scale genomics, transcriptomics, and proteomics studies. Specific features of its life cycle and the availability of a developmental mutant library make it an excellent system for studying cellular differentiation at the molecular level. In this review, we focus on recent developments in identifying gene and protein regulatory networks governing perithecia formation. A number of tools have been developed to genetically analyze developmental mutants and dissect transcriptional profiles at different developmental stages. Protein interaction studies allowed us to identify a highly conserved eukaryotic multisubunit protein complex, the striatin-interacting phosphatase and kinase complex and its role in sexual development. We have further identified a number of proteins involved in chromatin remodeling and transcriptional regulation of fruiting body development. Furthermore, we review the involvement of metabolic processes from both primary and secondary metabolism, and the role of nutrient recycling by autophagy in perithecia formation. Our research has uncovered numerous players regulating multicellular development in S. macrospora. Future research will focus on mechanistically understanding how these players are orchestrated in this fungal model system.

Molecular analysis of intragenic recombination at the tryptophan synthetase locus in Neurospora crassa

Fifteen different classically generated and mapped mutations at the tryptophan synthetase locus in Neurospora crassa have been characterized to the level of the primary sequence of the gene. This sequence analysis has demonstrated that intragenic recombination is accurate to order mutations within one open reading frame. While classic genetic analysis correctly ordered the mutations, the position of mutations characterized by gene sequence analysis was more accurate. A leaky mutation was found to have a wild-type primary sequence. The presence of unique polymorphisms in the primary sequence of the trp-3 gene from strain 861 confirms that it has a unique history relative to the other strains studied. Most strains that were previously shown to be immunologically nonreactive with antibody preparations raised against tryptophan synthetase protein were shown to have nonsense mutations. This work defines 14 alleles of the N. crassa trp-3 gene.

MycoCosm portal: gearing up for 1000 fungal genomes

MycoCosm is a fungal genomics portal (http://jgi.doe.gov/fungi), developed by the US Department of Energy Joint Genome Institute to support integration, analysis and dissemination of fungal genome sequences and other ‘omics’ data by providing interactive web-based tools. MycoCosm also promotes and facilitates user community participation through the nomination of new species of fungi for sequencing, and the annotation and analysis of resulting data. By efficiently filling gaps in the Fungal Tree of Life, MycoCosm will help address important problems associated with energy and the environment, taking advantage of growing fungal genomics resources.

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.

Comparative genome structure, secondary metabolite, and effector coding capacity across Cochliobolus pathogens

The genomes of five Cochliobolus heterostrophus strains, two Cochliobolus sativus strains, three additional Cochliobolus species (Cochliobolus victoriae, Cochliobolus carbonum, Cochliobolus miyabeanus), and closely related Setosphaeria turcica were sequenced at the Joint Genome Institute (JGI). The datasets were used to identify SNPs between strains and species, unique genomic regions, core secondary metabolism genes, and small secreted protein (SSP) candidate effector encoding genes with a view towards pinpointing structural elements and gene content associated with specificity of these closely related fungi to different cereal hosts. Whole-genome alignment shows that three to five percent of each genome differs between strains of the same species, while a quarter of each genome differs between species. On average, SNP counts among field isolates of the same C. heterostrophus species are more than 25× higher than those between inbred lines and 50× lower than SNPs between Cochliobolus species. The suites of nonribosomal peptide synthetase (NRPS), polyketide synthase (PKS), and SSP-encoding genes are astoundingly diverse among species but remarkably conserved among isolates of the same species, whether inbred or field strains, except for defining examples that map to unique genomic regions. Functional analysis of several strain-unique PKSs and NRPSs reveal a strong correlation with a role in virulence.

Centromeres of filamentous fungi

How centromeres are assembled and maintained remains one of the fundamental questions in cell biology. Over the past 20 years, the idea of centromeres as precise genetic loci has been replaced by the realization that it is predominantly the protein complement that defines centromere localization and function. Thus, placement and maintenance of centromeres are excellent examples of epigenetic phenomena in the strict sense. In contrast, the highly derived "point centromeres" of the budding yeast Saccharomyces cerevisiae and its close relatives are counter-examples for this general principle of centromere maintenance. While we have learned much in the past decade, it remains unclear if mechanisms for epigenetic centromere placement and maintenance are shared among various groups of organisms. For that reason, it seems prudent to examine species from many different phylogenetic groups with the aim to extract comparative information that will yield a more complete picture of cell division in all eukaryotes. This review addresses what has been learned by studying the centromeres of filamentous fungi, a large, heterogeneous group of organisms that includes important plant, animal and human pathogens, saprobes, and symbionts that fulfill essential roles in the biosphere, as well as a growing number of taxa that have become indispensable for industrial use.

Molecular analysis of mutants of the Neurospora adenylosuccinate synthetase locus

The ad-8 gene of Neurospora crassa, in addition to being used for the study of purine biology, has been extensively studied as a model for gene structure, mutagenesis and intralocus recombination. Because of this there is an extensive collection of well-characterized N. crassa ad-8 mutants in the Fungal Genetics Stock Center collection. Among these are spontaneous mutants and mutants induced with X-ray, UV or chemical mutagens. The specific lesions in these mutants have been genetically mapped at high resolution. We have sequenced the ad-8 locus from 13 of these mutants and identified the molecular nature of the mutation in each strain. We compare the historical fine-structure map to the DNA and amino acid sequence of each allele. The placement of the individual lesions in the fine-structure map was more accurate at the 5' end of the gene and no mutants were identified in the 3' untranslated region of this gene. We additionally analysed ad-8(+) alleles in 18 N. crassa strains subjected to whole-genome sequence analysis and describe the variability among Neurospora strains and among fungi and other organisms.

Analysis of complete nucleotide sequences of 12 Gossypium chloroplast genomes: origin and evolution of allotetraploids

BACKGROUND

Cotton (Gossypium spp.) is a model system for the analysis of polyploidization. Although ascertaining the donor species of allotetraploid cotton has been intensively studied, sequence comparison of Gossypium chloroplast genomes is still of interest to understand the mechanisms underlining the evolution of Gossypium allotetraploids, while it is generally accepted that the parents were A- and D-genome containing species. Here we performed a comparative analysis of 13 Gossypium chloroplast genomes, twelve of which are presented here for the first time.

METHODOLOGY/PRINCIPAL FINDINGS

The size of 12 chloroplast genomes under study varied from 159,959 bp to 160,433 bp. The chromosomes were highly similar having >98% sequence identity. They encoded the same set of 112 unique genes which occurred in a uniform order with only slightly different boundary junctions. Divergence due to indels as well as substitutions was examined separately for genome, coding and noncoding sequences. The genome divergence was estimated as 0.374% to 0.583% between allotetraploid species and A-genome, and 0.159% to 0.454% within allotetraploids. Forty protein-coding genes were completely identical at the protein level, and 20 intergenic sequences were completely conserved. The 9 allotetraploids shared 5 insertions and 9 deletions in whole genome, and 7-bp substitutions in protein-coding genes. The phylogenetic tree confirmed a close relationship between allotetraploids and the ancestor of A-genome, and the allotetraploids were divided into four separate groups. Progenitor allotetraploid cotton originated 0.43-0.68 million years ago (MYA).

CONCLUSION

Despite high degree of conservation between the Gossypium chloroplast genomes, sequence variations among species could still be detected. Gossypium chloroplast genomes preferred for 5-bp indels and 1-3-bp indels are mainly attributed to the SSR polymorphisms. This study supports that the common ancestor of diploid A-genome species in Gossypium is the maternal source of extant allotetraploid species and allotetraploids have a monophyletic origin. G. hirsutum AD1 lineages have experienced more sequence variations than other allotetraploids in intergenic regions. The available complete nucleotide sequences of 12 Gossypium chloroplast genomes should facilitate studies to uncover the molecular mechanisms of compartmental co-evolution and speciation of Gossypium allotetraploids.

Variation in mitochondrial genome primary sequence among whole-genome-sequenced strains of Neurospora crassa

Eighteen classical mutant strains of Neurospora crassa were subject to whole genome sequence analysis and the mitochondrial genome is analyzed. Overall, the mitochondrial genomes of the classical mutant strains are 99.45 to 99.98 % identical to the reference genome. Two-thirds of the SNPs and three-fourths of indels identified in this analysis are shared among more than one strain. Most of the limited variability in mitochondrial genome sequence is neutral with regard to protein structure. Despite the fact that the mitochondrial genome is present in multiple copies per cell, many of the polymorphisms were homozygous within each strain. Conversely, some polymorphisms, especially those associated with large scale rearrangements are only present in a fraction of the reads covering each region. The impact of this variation is unknown and further studies will be necessary to ascertain if this level of polymorphism is common among fungi and whether it reflects the impact of ageing cultures.

Whole-Genome Sequencing of Sordaria macrospora Mutants Identifies Developmental Genes

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.

A versatile toolkit for high throughput functional genomics with Trichoderma reesei

BACKGROUND

The ascomycete fungus, Trichoderma reesei (anamorph of Hypocrea jecorina), represents a biotechnological workhorse and is currently one of the most proficient cellulase producers. While strain improvement was traditionally accomplished by random mutagenesis, a detailed understanding of cellulase regulation can only be gained using recombinant technologies.

RESULTS

Aiming at high efficiency and high throughput methods, we present here a construction kit for gene knock out in T. reesei. We provide a primer database for gene deletion using the pyr4, amdS and hph selection markers. For high throughput generation of gene knock outs, we constructed vectors using yeast mediated recombination and then transformed a T. reesei strain deficient in non-homologous end joining (NHEJ) by spore electroporation. This NHEJ-defect was subsequently removed by crossing of mutants with a sexually competent strain derived from the parental strain, QM9414.

CONCLUSIONS

Using this strategy and the materials provided, high throughput gene deletion in T. reesei becomes feasible. Moreover, with the application of sexual development, the NHEJ-defect can be removed efficiently and without the need for additional selection markers. The same advantages apply for the construction of multiple mutants by crossing of strains with different gene deletions, which is now possible with considerably less hands-on time and minimal screening effort compared to a transformation approach. Consequently this toolkit can considerably boost research towards efficient exploitation of the resources of T. reesei for cellulase expression and hence second generation biofuel production.