Clustering analysis of large-scale phenotypic data in the model filamentous fungus Neurospora crassa.

BACKGROUND: With 9730 protein-coding genes and a nearly complete gene knockout strain collection, Neurospora crassa is a major model organism for filamentous fungi. Despite this abundance of information, the phenotypes of these gene knockout mutants have not been categorized to determine whether there are broad correlations between phenotype and any genetic features. RESULTS: Here, we analyze data for 10 different growth or developmental phenotypes that have been obtained for 1168 N. crassa knockout mutants. Of these mutants, 265 (23%) are in the normal range, while 903 (77%) possess at least one mutant phenotype. With the exception of unclassified functions, the distribution of functional categories for genes in the mutant dataset mirrors that of the N. crassa genome. In contrast, most genes do not possess a yeast ortholog, suggesting that our analysis will reveal functions that are not conserved in Saccharomyces cerevisiae. To leverage the phenotypic data to identify pathways, we used weighted Partitioning Around Medoids (PAM) approach with 40 clusters. We found that genes encoding metabolic, transmembrane and protein phosphorylation-related genes are concentrated in subsets of clusters. Results from K-Means clustering of transcriptomic datasets showed that most phenotypic clusters contain multiple expression profiles, suggesting that co-expression is not generally observed for genes with shared phenotypes. Analysis of yeast orthologs of genes that co-clustered in MAPK signaling cascades revealed potential networks of interacting proteins in N. crassa. CONCLUSIONS: Our results demonstrate that clustering analysis of phenotypes is a promising tool for generating new hypotheses regarding involvement of genes in cellular pathways in N. crassa. Furthermore, information about gene clusters identified in N. crassa should be applicable to other filamentous fungi, including saprobes and pathogens.

Regulator of G Protein Signaling Proteins Control Growth, Development and Cellulase Production in Neurospora crassa.

Heterotrimeric (αßγ) G protein signaling pathways are critical environmental sensing systems found in eukaryotic cells. Exchange of GDP for GTP on the Gα subunit leads to its activation. In contrast, GTP hydrolysis on the Gα is accelerated by Regulator of G protein Signaling (RGS) proteins, resulting in a return to the GDP-bound, inactive state. Here, we analyzed growth, development and extracellular cellulase production in strains with knockout mutations in the seven identified RGS genes (rgs-1 to rgs-7) in the filamentous fungus, Neurospora crassa. We compared phenotypes to those of strains with either knockout mutations or expressing predicted constitutively activated, GTPase-deficient alleles for each of the three Gα subunit genes (gna-1Q204L, gna-2Q205L or gna-3Q208L). Our data revealed that six RGS mutants have taller aerial hyphae than wild type and all seven mutants exhibit reduced asexual sporulation, phenotypes shared with strains expressing the gna-1Q204L or gna-3Q208L allele. In contrast, Δrgs-1 and Δrgs-3 were the only RGS mutants with a slower growth rate phenotype, a defect in common with gna-1Q204L strains. With respect to female sexual development, Δrgs-1 possessed defects most similar to gna-3Q208L strains, while those of Δrgs-2 mutants resembled strains expressing the gna-1Q204L allele. Finally, we observed that four of the seven RGS mutants had significantly different extracellular cellulase levels relative to wild type. Of interest, the Δrgs-2 mutant had no detectable activity, similar to the gna-3Q208L strain. In contrast, the Δrgs-1 and Δrgs-4 mutants and gna-1Q204L and gna-2Q205L strains exhibited significantly higher cellulase activity than wild type. With the exception of sexual development, our results demonstrate the greatest number of genetic interactions between rgs-1 and gna-1 and rgs-2 and gna-3 in N. crassa.

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 Zn2Cys6 (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.

Global Analysis of Predicted G Protein-Coupled Receptor Genes in the Filamentous Fungus, Neurospora crassa.

G protein-coupled receptors (GPCRs) regulate facets of growth, development, and environmental sensing in eukaryotes, including filamentous fungi. The largest predicted GPCR class in these organisms is the Pth11-related, with members similar to a protein required for disease in the plant pathogen Magnaporthe oryzae. However, the Pth11-related class has not been functionally studied in any filamentous fungal species. Here, we analyze phenotypes in available mutants for 36 GPCR genes, including 20 Pth11-related, in the model filamentous fungus Neurospora crassa. We also investigate patterns of gene expression for all 43 predicted GPCR genes in available datasets. A total of 17 mutants (47%) possessed at least one growth or developmental phenotype. We identified 18 mutants (56%) with chemical sensitivity or nutritional phenotypes (11 uniquely), bringing the total number of mutants with at least one defect to 28 (78%), including 15 mutants (75%) in the Pth11-related class. Gene expression trends for GPCR genes correlated with the phenotypes observed for many mutants and also suggested overlapping functions for several groups of co-transcribed genes. Several members of the Pth11-related class have phenotypes and/or are differentially expressed on cellulose, suggesting a possible role for this gene family in plant cell wall sensing or utilization.

The guanine nucleotide exchange factor RIC8 regulates conidial germination through Gα proteins in Neurospora crassa.

Heterotrimeric G protein signaling is essential for normal hyphal growth in the filamentous fungus Neurospora crassa. We have previously demonstrated that the non-receptor guanine nucleotide exchange factor RIC8 acts upstream of the Gα proteins GNA-1 and GNA-3 to regulate hyphal extension. Here we demonstrate that regulation of hyphal extension results at least in part, from an important role in control of asexual spore (conidia) germination. Loss of GNA-3 leads to a drastic reduction in conidial germination, which is exacerbated in the absence of GNA-1. Mutation of RIC8 leads to a reduction in germination similar to that in the Δgna-1, Δgna-3 double mutant, suggesting that RIC8 regulates conidial germination through both GNA-1 and GNA-3. Support for a more significant role for GNA-3 is indicated by the observation that expression of a GTPase-deficient, constitutively active gna-3 allele in the Δric8 mutant leads to a significant increase in conidial germination. Localization of the three Gα proteins during conidial germination was probed through analysis of cells expressing fluorescently tagged proteins. Functional TagRFP fusions of each of the three Gα subunits were constructed through insertion of TagRFP in a conserved loop region of the Gα subunits. The results demonstrated that GNA-1 localizes to the plasma membrane and vacuoles, and also to septa throughout conidial germination. GNA-2 and GNA-3 localize to both the plasma membrane and vacuoles during early germination, but are then found in intracellular vacuoles later during hyphal outgrowth.