The filamentous fungal gene expression database (FFGED)

Rust expression browser: an open source database for simultaneous analysis of host and pathogen gene expression profiles with expVIP

Background

Transcriptomics is being increasingly applied to generate new insight into the interactions between plants and their pathogens. For the wheat yellow (stripe) rust pathogen (Puccinia striiformis f. sp. tritici, Pst) RNA-based sequencing (RNA-Seq) has proved particularly valuable, overcoming the barriers associated with its obligate biotrophic nature. This includes the application of RNA-Seq approaches to study Pst and wheat gene expression dynamics over time and the Pst population composition through the use of a novel RNA-Seq based surveillance approach called “field pathogenomics”. As a dual RNA-Seq approach, the field pathogenomics technique also provides gene expression data from the host, giving new insight into host responses. However, this has created a wealth of data for interrogation.

Results

Here, we used the field pathogenomics approach to generate 538 new RNA-Seq datasets from Pst-infected field wheat samples, doubling the amount of transcriptomics data available for this important pathosystem. We then analysed these datasets alongside 66 RNA-Seq datasets from four Pst infection time-courses and 420 Pst-infected plant field and laboratory samples that were publicly available. A database of gene expression values for Pst and wheat was generated for each of these 1024 RNA-Seq datasets and incorporated into the development of the rust expression browser (http://www.rust-expression.com). This enables for the first time simultaneous ‘point-and-click’ access to gene expression profiles for Pst and its wheat host and represents the largest database of processed RNA-Seq datasets available for any of the three Puccinia wheat rust pathogens. We also demonstrated the utility of the browser through investigation of expression of putative Pst virulence genes over time and examined the host plants response to Pst infection.

Conclusions

The rust expression browser offers immense value to the wider community, facilitating data sharing and transparency and the underlying database can be continually expanded as more datasets become publicly available.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07488-3.

Fungi.guru: Comparative genomic and transcriptomic resource for the fungi kingdom

The fungi kingdom is composed of eukaryotic heterotrophs, which are responsible for balancing the ecosystem and play a major role as decomposers. They also produce a vast diversity of secondary metabolites, which have antibiotic or pharmacological properties. However, our lack of knowledge of gene function in fungi precludes us from tailoring them to our needs and tapping into their metabolic diversity. To help remedy this, we gathered genomic and gene expression data of 19 most widely-researched fungi to build an online tool, fungi.guru, which contains tools for cross-species identification of conserved pathways, functional gene modules, and gene families. We exemplify how our tool can elucidate the molecular function, biological process and cellular component of genes involved in various biological processes, by identifying a secondary metabolite pathway producing gliotoxin in Aspergillus fumigatus, the catabolic pathway of cellulose in Coprinopsis cinerea and the conserved DNA replication pathway in Fusarium graminearum and Pyricularia oryzae. The tool is available at www.fungi.guru.

Compartmentalized gene regulatory network of the pathogenic fungus Fusarium graminearum

Head blight caused by Fusarium graminearum threatens world-wide wheat production, resulting in both yield loss and mycotoxin contamination. We reconstructed the global F. graminearum gene regulatory network (GRN) from a large collection of transcriptomic data using Bayesian network inference, a machine-learning algorithm. This GRN reveals connectivity between key regulators and their target genes. Focusing on key regulators, this network contains eight distinct but interwoven modules. Enriched for unique functions, such as cell cycle, DNA replication, transcription, translation and stress responses, each module exhibits distinct expression profiles. Evolutionarily, the F. graminearum genome can be divided into core regions shared with closely related species and variable regions harboring genes that are unique to F. graminearum and perform species-specific functions. Interestingly, the inferred top regulators regulate genes that are significantly enriched from the same genomic regions (P < 0.05), revealing a compartmentalized network structure that may reflect network rewiring related to specific adaptation of this plant pathogen. This first-ever reconstructed filamentous fungal GRN primes our understanding of pathogenicity at the systems biology level and provides enticing prospects for novel disease control strategies involving the targeting of master regulators in pathogens. The program can be used to construct GRNs of other plant pathogens.

Conservation and divergence of the cyclic adenosine monophosphate-protein kinase A (cAMP-PKA) pathway in two plant-pathogenic fungi: Fusarium graminearum and F. verticillioides

The cyclic adenosine monophosphate-protein kinase A (cAMP-PKA) pathway is a central signalling cascade that transmits extracellular stimuli and governs cell responses through the second messenger cAMP. The importance of cAMP signalling in fungal biology has been well documented and the key conserved components, adenylate cyclase (AC) and the catalytic subunit of PKA (CPKA), have been functionally characterized. However, other genes involved in this signalling pathway and their regulation are not well understood in filamentous fungi. Here, we performed a comparative transcriptomics analysis of AC and CPKA mutants in two closely related fungi: Fusarium graminearum (Fg) and F. verticillioides (Fv). Combining available Fg transcriptomics and phenomics data, we reconstructed the Fg cAMP signalling pathway. We developed a computational program that combines sequence conservation and patterns of orthologous gene expression to facilitate global transcriptomics comparisons between different organisms. We observed highly correlated expression patterns for most orthologues (80%) between Fg and Fv. We also identified a subset of 482 (6%) diverged orthologues, whose expression under all conditions was at least 50% higher in one genome than in the other. This enabled us to dissect the conserved and unique portions of the cAMP-PKA pathway. Although the conserved portions controlled essential functions, such as metabolism, the cell cycle, chromatin remodelling and the oxidative stress response, the diverged portions had species-specific roles, such as the production and detoxification of secondary metabolites unique to each species. The evolution of the cAMP-PKA signalling pathway seems to have contributed directly to fungal divergence and niche adaptation.

Gene expression differences among three Neurospora species reveal genes required for sexual reproduction in Neurospora crassa

Many fungi form complex three-dimensional fruiting bodies, within which the meiotic machinery for sexual spore production has been considered to be largely conserved over evolutionary time. Indeed, much of what we know about meiosis in plant and animal taxa has been deeply informed by studies of meiosis in Saccharomyces and Neurospora. Nevertheless, the genetic basis of fruiting body development and its regulation in relation to meiosis in fungi is barely known, even within the best studied multicellular fungal model Neurospora crassa. We characterized morphological development and genome-wide transcriptomics in the closely related species Neurospora crassa, Neurospora tetrasperma, and Neurospora discreta, across eight stages of sexual development. Despite diverse life histories within the genus, all three species produce vase-shaped perithecia. Transcriptome sequencing provided gene expression levels of orthologous genes among all three species. Expression of key meiosis genes and sporulation genes corresponded to known phenotypic and developmental differences among these Neurospora species during sexual development. We assembled a list of genes putatively relevant to the recent evolution of fruiting body development by sorting genes whose relative expression across developmental stages increased more in N. crassa relative to the other species. Then, in N. crassa, we characterized the phenotypes of fruiting bodies arising from crosses of homozygous knockout strains of the top genes. Eight N. crassa genes were found to be critical for the successful formation of perithecia. The absence of these genes in these crosses resulted in either no perithecium formation or in arrested development at an early stage. Our results provide insight into the genetic basis of Neurospora sexual reproduction, which is also of great importance with regard to other multicellular ascomycetes, including perithecium-forming pathogens, such as Claviceps purpurea, Ophiostoma ulmi, and Glomerella graminicola.

The aging biological clock in Neurospora crassa

The biological clock affects aging through ras-1 (bd) and lag-1, and these two longevity genes together affect a clock phenotype and the clock oscillator in Neurospora crassa. Using an automated cell-counting technique for measuring conidial longevity, we show that the clock-associated genes lag-1 and ras-1 (bd) are true chronological longevity genes. For example, wild type (WT) has an estimated median life span of 24 days, while the double mutant lag-1, ras-1 (bd) has an estimated median life span of 120 days for macroconidia. We establish the biochemical function of lag-1 by complementing LAG1 and LAC1 in Saccharomyces cerevisiae with lag-1 in N. crassa. Longevity genes can affect the clock as well in that, the double mutant lag-1, ras-1 (bd) can stop the circadian rhythm in asexual reproduction (i.e., banding in race tubes) and lengthen the period of the frequency oscillator to 41 h. In contrast to the ras-1 (bd), lag-1 effects on chronological longevity, we find that this double mutant undergoes replicative senescence (i.e., the loss of replication function with time), unlike WT or the single mutants, lag-1 and ras-1 (bd). These results support the hypothesis that sphingolipid metabolism links aging and the biological clock through a common stress response

Global gene expression and focused knockout analysis reveals genes associated with fungal fruiting body development in Neurospora crassa

Fungi can serve as highly tractable models for understanding genetic basis of sexual development in multicellular organisms. Applying a reverse-genetic approach to advance such a model, we used random and multitargeted primers to assay gene expression across perithecial development in Neurospora crassa. We found that functionally unclassified proteins accounted for most upregulated genes, whereas downregulated genes were enriched for diverse functions. Moreover, genes associated with developmental traits exhibited stage-specific peaks of expression. Expression increased significantly across sexual development for mating type gene mat a-1 and for mat A-1 specific pheromone precursor ccg-4. In addition, expression of a gene encoding a protein similar to zinc finger, stc1, was highly upregulated early in perithecial development, and a strain with a knockout of this gene exhibited arrest at the same developmental stage. A similar expression pattern was observed for genes in RNA silencing and signaling pathways, and strains with knockouts of these genes were also arrested at stages of perithecial development that paralleled their peak in expression. The observed stage specificity allowed us to correlate expression upregulation and developmental progression and to identify regulators of sexual development. Bayesian networks inferred from our expression data revealed previously known and new putative interactions between RNA silencing genes and pathways. Overall, our analysis provides a fine-scale transcriptomic landscape and novel inferences regarding the control of the multistage development process of sexual crossing and fruiting body development in N. crassa.

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.

Abundant gene-by-environment interactions in gene expression reaction norms to copper within Saccharomyces cerevisiae

Genetic variation for plastic phenotypes potentially contributes phenotypic variation to populations that can be selected during adaptation to novel ecological contexts. However, the basis and extent of plastic variation that manifests in diverse environments remains elusive. Here, we characterize copper reaction norms for mRNA abundance among five Saccharomyces cerevisiae strains to 1) describe population variation across the full range of ecologically relevant copper concentrations, from starvation to toxicity, and 2) to test the hypothesis that plastic networks exhibit increased population variation for gene expression. We find that although the vast majority of the variation is small in magnitude (considerably <2-fold), not just some, but most genes demonstrate variable expression across environments, across genetic backgrounds, or both. Plastically expressed genes included both genes regulated directly by copper-binding transcription factors Mac1 and Ace1 and genes indirectly responding to the downstream metabolic consequences of the copper gradient, particularly genes involved in copper, iron, and sulfur homeostasis. Copper-regulated gene networks exhibited more similar behavior within the population in environments where those networks have a large impact on fitness. Nevertheless, expression variation in genes like Cup1, important to surviving copper stress, was linked with variation in mitotic fitness and in the breadth of differential expression across the genome. By revealing a broader and deeper range of population variation, our results provide further evidence for the interconnectedness of genome-wide mRNA levels, their dependence on environmental context and genetic background, and the abundance of variation in gene expression that can contribute to future evolution.

Analyses of expressed sequence tags in Neurospora reveal rapid evolution of genes associated with the early stages of sexual reproduction in fungi

Background

The broadly accepted pattern of rapid evolution of reproductive genes is primarily based on studies of animal systems, although several examples of rapidly evolving genes involved in reproduction are found in diverse additional taxa. In fungi, genes involved in mate recognition have been found to evolve rapidly. However, the examples are too few to draw conclusions on a genome scale.

Results

In this study, we performed microarray hybridizations between RNA from sexual and vegetative tissues of two strains of the heterothallic (self-sterile) filamentous ascomycete Neurospora intermedia, to identify a set of sex-associated genes in this species. We aligned Expressed Sequence Tags (ESTs) from sexual and vegetative tissue of N. intermedia to orthologs from three closely related species: N. crassa, N. discreta and N. tetrasperma. The resulting four-species alignments provided a dataset for molecular evolutionary analyses. Our results confirm a general pattern of rapid evolution of fungal sex-associated genes, compared to control genes with constitutive expression or a high relative expression during vegetative growth. Among the rapidly evolving sex-associated genes, we identified candidates that could be of importance for mating or fruiting-body development. Analyses of five of these candidate genes from additional species of heterothallic Neurospora revealed that three of them evolve under positive selection.

Conclusions

Taken together, our study represents a novel finding of a genome-wide pattern of rapid evolution of sex-associated genes in the fungal kingdom, and provides a list of candidate genes important for reproductive isolation in Neurospora.

Diversification of a protein kinase cascade: IME-2 is involved in nonself recognition and programmed cell death in Neurospora crassa

Kinase cascades and the modification of proteins by phosphorylation are major mechanisms for cell signaling and communication, and evolution of these signaling pathways can contribute to new developmental or environmental response pathways. The Saccharomyces cerevisiae kinase Ime2 has been well characterized for its role in meiosis. However, recent studies have revealed alternative functions for Ime2 in both S. cerevisiae and other fungi. In the filamentous fungus Neurospora crassa, the IME2 homolog (ime-2) is not required for meiosis. Here we determine that ime-2 interacts genetically with a transcription factor vib-1 during nonself recognition and programmed cell death (PCD). Mutations in vib-1 (Δvib-1) suppress PCD due to nonself recognition events; however, a Δvib-1 Δime-2 mutant restored wild-type levels of cell death. A role for ime-2 in the post-translational processing and localization of a mitochondrial matrix protein was identified, which may implicate mitochondria in N. crassa nonself recognition and PCD. Further, Δvib-1 strains do not produce extracellular proteases, but protease secretion reverted to near wild-type levels in a Δvib-1 Δime-2 strain. Mass spectrometry analysis revealed that the VIB-1 protein is phosphorylated at several sites, including a site that matches the IME-2 consensus. The genetic and biochemical data for ime-2 and vib-1 indicate that IME-2 is a negative regulator of VIB-1 and suggest parallel negative regulation by IME-2 of a cell death pathway in N. crassa that functions in concert with the VIB-1 cell death pathway. Thus, IME2 kinase function has evolved following the divergence of S. cerevisiae and N. crassa and provides insight into the evolution of kinases and their regulatory targets.

Transcriptome analyses during fruiting body formation in Fusarium graminearum and Fusarium verticillioides reflect species life history and ecology

Fusarium graminearum and F. verticillioides are devastating cereal pathogens with very different life history and ecological characteristics. F. graminearum is homothallic, and sexual spores are an important component of its life cycle, responsible for disease initiation. F. verticilloides is heterothallic, and produces only modest numbers of fruiting bodies, which are not a significant source of inoculum. To identify corresponding differences in the transcriptional program underlying fruiting body development in the two species, comparative expression was performed, analyzing six developmental stages. To accompany the transcriptional analysis, detailed morphological characterization of F. verticillioides development was performed and compared to a previous morphological analysis of F. graminearum. Morphological development was similar between the two species, except for the observation of possible trichogynes in F. verticillioides ascogonia, which have not been previously reported for any Fusarium species. Expression of over 9000 orthologous genes were measured for the two species. Functional assignments of highly expressed orthologous genes at each time-point revealed the majority of highly expressed genes fell into the "unclassified proteins" category, reflecting the lack of characterization of genes for sexual development in both species. Simultaneous examination of morphological development and stage-specific gene expression suggests that degeneration of the paraphyses during sexual development is an apoptotic process. Expression of mating type genes in the two species differed, possibly reflecting the divergent roles they play in sexual development. Overall, the differences in gene expression reflect the greater role of fruiting bodies in the life cycle and ecology of F. graminearum versus F. verticillioides.

Sex-specific gene expression during asexual development of Neurospora crassa

The impact of loci that determine sexual identity upon the asexual, dominant stage of fungal life history has been well studied. To investigate their impact, expression differences between strains of different mating type during asexual development were assayed, with RNA sampled from otherwise largely isogenic mat A and mat a strains of Neurospora crassa at early, middle, and late clonal stages of development. We observed significant differences in overall gene expression between mating types across clonal development, especially at late development stages. The expression levels of mating-type genes and pheromone genes were assayed by reverse transcription and quantitative PCR, revealing expression of pheromone and receptor genes in strains of both mating types in all development stages, and revealing that mating type (mat) genes were increasingly expressed over the course of asexual development. Interestingly, among differentially expressed genes, the mat A genotype more frequently exhibited a higher expression level than mat a, and demonstrated greater transcriptional regulatory dynamism. Significant up-regulation of expression was observed for many late light-responsive genes at late asexual development stages. Further investigation of the impact of light and the roles of light response genes in asexual development of both mating types are warranted.

Phenotypic diversification is associated with host-induced transposon derepression in the sudden oak death pathogen Phytophthora ramorum

The oomycete pathogen Phytophthora ramorum is responsible for sudden oak death (SOD) in California coastal forests. P. ramorum is a generalist pathogen with over 100 known host species. Three or four closely related genotypes of P. ramorum (from a single lineage) were originally introduced in California forests and the pathogen reproduces clonally. Because of this the genetic diversity of P. ramorum is extremely low in Californian forests. However, P. ramorum shows diverse phenotypic variation in colony morphology, colony senescence, and virulence. In this study, we show that phenotypic variation among isolates is associated with the host species from which the microbe was originally cultured. Microarray global mRNA profiling detected derepression of transposable elements (TEs) and down-regulation of crinkler effector homologs (CRNs) in the majority of isolates originating from coast live oak (Quercus agrifolia), but this expression pattern was not observed in isolates from California bay laurel (Umbellularia californica). In some instances, oak and bay laurel isolates originating from the same geographic location had identical genotypes based on multilocus simples sequence repeat (SSR) marker analysis but had different phenotypes. Expression levels of the two marker genes analyzed by quantitative reverse transcription PCR were correlated with originating host species, but not with multilocus genotypes. Because oak is a nontransmissive dead-end host for P. ramorum, our observations are congruent with an epi-transposon hypothesis; that is, physiological stress is triggered on P. ramorum while colonizing oak stems and disrupts epigenetic silencing of TEs. This then results in TE reactivation and possibly genome diversification without significant epidemiological consequences. We propose the P. ramorum-oak host system in California forests as an ad hoc model for epi-transposon mediated diversification.

Differential impact of nutrition on developmental and metabolic gene expression during fruiting body development in Neurospora crassa

Fungal fruiting body size and form are influenced by the ecology of the species, including diverse environmental stimuli. Accordingly, nutritional resources available to the fungus during development can be vital to successful production of fruiting bodies. To investigate the effect of nutrition, perithecial development of Neurospora crassa was induced on two different media, a chemically sparsely nutritive Synthetic Crossing Medium (SCM) and a natural Carrot Agar (CA). Protoperithecia were collected before crossing, and perithecia were collected at 2, 24, 48, 72, 96, 120, and at full maturity 144 h after crossing. No differences in fruiting body morphology were observed between the two media at any time point. A circuit of microarray hybridizations comparing cDNA from all neighboring stages was performed. For a majority of differentially expressed genes, expression was higher in SCM than in CA, and expression of core metabolic genes was particularly affected. Effects of nutrition were highest in magnitude before crossing, lowering in magnitude during early perithecial development. Interestingly, metabolic effects of the media were also large in magnitude during late perithecial development, at which stage the lower expression in CA presumably reflected the continued intake of diverse complex initial compounds, diminishing the need for expression of anabolic pathways. However, for genes with key regulatory roles in sexual development, including pheromone precursor ccg-4 and poi2, expression patterns were similar between treatments. When possible, a common nutritional environment is ideal for comparing transcriptional profiles between different fungi. Nevertheless, the observed consistency of the developmental program across media, despite considerable metabolic differentiation is reassuring. This result facilitates comparative studies that will require different nutritional resources for sexual development in different fungi.

Systems biology of the qa gene cluster in Neurospora crassa

An ensemble of genetic networks that describe how the model fungal system, Neurospora crassa, utilizes quinic acid (QA) as a sole carbon source has been identified previously. A genetic network for QA metabolism involves the genes, qa-1F and qa-1S, that encode a transcriptional activator and repressor, respectively and structural genes, qa-2, qa-3, qa-4, qa-x, and qa-y. By a series of 4 separate and independent, model-guided, microarray experiments a total of 50 genes are identified as QA-responsive and hypothesized to be under QA-1F control and/or the control of a second QA-responsive transcription factor (NCU03643) both in the fungal binuclear Zn(II)2Cys6 cluster family. QA-1F regulation is not sufficient to explain the quantitative variation in expression profiles of the 50 QA-responsive genes. QA-responsive genes include genes with products in 8 mutually connected metabolic pathways with 7 of them one step removed from the tricarboxylic (TCA) Cycle and with 7 of them one step removed from glycolysis: (1) starch and sucrose metabolism; (2) glycolysis/glucanogenesis; (3) TCA Cycle; (4) butanoate metabolism; (5) pyruvate metabolism; (6) aromatic amino acid and QA metabolism; (7) valine, leucine, and isoleucine degradation; and (8) transport of sugars and amino acids. Gene products both in aromatic amino acid and QA metabolism and transport show an immediate response to shift to QA, while genes with products in the remaining 7 metabolic modules generally show a delayed response to shift to QA. The additional QA-responsive cutinase transcription factor-1β (NCU03643) is found to have a delayed response to shift to QA. The series of microarray experiments are used to expand the previously identified genetic network describing the qa gene cluster to include all 50 QA-responsive genes including the second transcription factor (NCU03643). These studies illustrate new methodologies from systems biology to guide model-driven discoveries about a core metabolic network involving carbon and amino acid metabolism in N. crassa.

Fungal pathogenicity genes in the age of 'omics'

The identification of the fungal genes essential for disease underpins the development of disease control strategies. Improved technologies for gene identification and functional analyses, as well as a plethora of sequenced fungal genomes, have led to the characterization of hundreds of genes, denoted as pathogenicity genes, which are required by fungi to cause disease. We describe recent technologies applied to characterize the fungal genes involved in disease and focus on some genes that are likely to attract continuing research activity.

Getting the most out of your fungal microarray data: two cost- and time-effective methods

Advances in genome sequencing technologies have facilitated production of a wealth of fungal data; within the last 5 years, experimental costs and labor have diminished, shifting the production bottleneck from genomic data generation to data analysis. Genome sequences and microarrays now exist for many fungi, and transcriptional profiling has been shown to be an efficient way to examine how the entire genome changes in response to many different environments or treatments. Multiple platforms, programs, and protocols exist for analyzing such data, making this task daunting for the bench-based scientist. Furthermore, many existing programs are expensive and require license renewals on a yearly basis for each user in the laboratory. Costs may be prohibitively high for bench-based scientists in academia. Our combined experiences with this kind of analysis have favored two programs, depending upon whether the scientist is working with single- or dual-channel hybridization data. Our protocols are aimed toward helping the bench-based PI get the most possible information from their data, without the need for expensive software or an experienced bioinformaticist.

Multi-targeted priming for genome-wide gene expression assays

BACKGROUND

Complementary approaches to assaying global gene expression are needed to assess gene expression in regions that are poorly assayed by current methodologies. A key component of nearly all gene expression assays is the reverse transcription of transcribed sequences that has traditionally been performed by priming the poly-A tails on many of the transcribed genes in eukaryotes with oligo-dT, or by priming RNA indiscriminately with random hexamers. We designed an algorithm to find common sequence motifs that were present within most protein-coding genes of Saccharomyces cerevisiae and of Neurospora crassa, but that were not present within their ribosomal RNA or transfer RNA genes. We then experimentally tested whether degenerately priming these motifs with multi-targeted primers improved the accuracy and completeness of transcriptomic assays.

RESULTS

We discovered two multi-targeted primers that would prime a preponderance of genes in the genomes of Saccharomyces cerevisiae and Neurospora crassa while avoiding priming ribosomal RNA or transfer RNA. Examining the response of Saccharomyces cerevisiae to nitrogen deficiency and profiling Neurospora crassa early sexual development, we demonstrated that using multi-targeted primers in reverse transcription led to superior performance of microarray profiling and next-generation RNA tag sequencing. Priming with multi-targeted primers in addition to oligo-dT resulted in higher sensitivity, a larger number of well-measured genes and greater power to detect differences in gene expression.

CONCLUSIONS

Our results provide the most complete and detailed expression profiles of the yeast nitrogen starvation response and N. crassa early sexual development to date. Furthermore, our multi-targeting priming methodology for genome-wide gene expression assays provides selective targeting of multiple sequences and counter-selection against undesirable sequences, facilitating a more complete and precise assay of the transcribed sequences within the genome.

Meiotic regulators Ndt80 and ime2 have different roles in Saccharomyces and Neurospora

Meiosis is a highly regulated process in eukaryotic species. The filamentous fungus Neurospora crassa has been shown to be missing homologs of a number of meiotic initiation genes conserved in Saccharomyces cerevisiae, but has three homologs of the well-characterized middle meiotic transcriptional regulator NDT80. In this study, we evaluated the role of all three NDT80 homologs in the formation of female reproductive structures, sexual development, and meiosis. We found that none of the NDT80 homologs were required for meiosis and that even the triple mutant was unaffected. However, strains containing mutations in NCU09915 (fsd-1) were defective in female sexual development and ascospore maturation. vib-1 was a major regulator of protoperithecial development in N. crassa, and double mutants carrying deletions of both vib-1 (NCU03725) and fsd-1 exhibited a synergistic effect on the timing of female reproductive structure (protoperithecia) formation. We further evaluated the role of the N. crassa homolog of IME2, a kinase involved in initiation of meiosis in S. cerevisiae. Strains containing mutations in ime-2 showed unregulated development of protoperithecia. Genetic analysis indicated that mutations in vib-1 were epistatic to ime-2, suggesting that IME-2 may negatively regulate VIB-1 activity. Our data indicate that the IME2/NDT80 pathway is not involved in meiosis in N. crassa, but rather regulates the formation of female reproductive structures.

Characterization of the temperature-sensitive mutations un-7 and png-1 in Neurospora crassa

The model filamentous fungus Neurospora crassa has been studied for over fifty years and many temperature-sensitive mutants have been generated. While most of these have been mapped genetically, many remain anonymous. The mutation in the N. crassa temperature-sensitive lethal mutant un-7 was identified by a complementation based approach as being in the open reading frame designated NCU00651 on linkage group I. Other mutations in this gene have been identified that lead to a temperature-sensitive morphological phenotype called png-1. The mutations underlying un-7 result in a serine to phenylalanine change at position 273 and an isoleucine to valine change at position 390, while the mutation in png-1 was found to result in a serine to leucine change at position 279 although there were other conservative changes in this allele. The overall morphology of the strain carrying the un-7 mutation is compared to strains carrying the png-1 mutation and these mutations are evaluated in the context of other temperature-sensitive mutants in Neurospora.