Relationship between phylogenetic distribution and genomic features in Neurospora crassa

Lineage-specific genes are clustered with HET-domain genes and respond to environmental and genetic manipulations regulating reproduction in Neurospora

Lineage-specific genes (LSGs) have long been postulated to play roles in the establishment of genetic barriers to intercrossing and speciation. In the genome of Neurospora crassa, most of the 670 Neurospora LSGs that are aggregated adjacent to the telomeres are clustered with 61% of the HET-domain genes, some of which regulate self-recognition and define vegetative incompatibility groups. In contrast, the LSG-encoding proteins possess few to no domains that would help to identify potential functional roles. Possible functional roles of LSGs were further assessed by performing transcriptomic profiling in genetic mutants and in response to environmental alterations, as well as examining gene knockouts for phenotypes. Among the 342 LSGs that are dynamically expressed during both asexual and sexual phases, 64% were detectable on unusual carbon sources such as furfural, a wildfire-produced chemical that is a strong inducer of sexual development, and the structurally-related furan 5-hydroxymethyl furfural (HMF). Expression of a significant portion of the LSGs was sensitive to light and temperature, factors that also regulate the switch from asexual to sexual reproduction. Furthermore, expression of the LSGs was significantly affected in the knockouts of adv-1 and pp-1 that regulate hyphal communication, and expression of more than one quarter of the LSGs was affected by perturbation of the mating locus. These observations encouraged further investigation of the roles of clustered lineage-specific and HET-domain genes in ecology and reproduction regulation in Neurospora, especially the regulation of the switch from the asexual growth to sexual reproduction, in response to dramatic environmental conditions changes.

The Lost and Found: Unraveling the Functions of Orphan Genes

Orphan Genes (OGs) are a mysterious class of genes that have recently gained significant attention. Despite lacking a clear evolutionary history, they are found in nearly all living organisms, from bacteria to humans, and they play important roles in diverse biological processes. The discovery of OGs was first made through comparative genomics followed by the identification of unique genes across different species. OGs tend to be more prevalent in species with larger genomes, such as plants and animals, and their evolutionary origins remain unclear but potentially arise from gene duplication, horizontal gene transfer (HGT), or de novo origination. Although their precise function is not well understood, OGs have been implicated in crucial biological processes such as development, metabolism, and stress responses. To better understand their significance, researchers are using a variety of approaches, including transcriptomics, functional genomics, and molecular biology. This review offers a comprehensive overview of the current knowledge of OGs in all domains of life, highlighting the possible role of dark transcriptomics in their evolution. More research is needed to fully comprehend the role of OGs in biology and their impact on various biological processes.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

Secondary Metabolism Gene Clusters Exhibit Increasingly Dynamic and Differential Expression during Asexual Growth, Conidiation, and Sexual Development in Neurospora crassa

Secondary metabolite clusters (SMCs) encode the machinery for fungal toxin production. However, understanding their function and analyzing their products requires investigation of the developmental and environmental conditions in which they are expressed. Gene expression is often restricted to specific and unexamined stages of the life cycle. Therefore, we applied comparative genomics analyses to identify SMCs in Neurospora crassa and analyzed extensive transcriptomic data spanning nine independent experiments from diverse developmental and environmental conditions to reveal their life cycle-specific gene expression patterns. We reported 20 SMCs comprising 177 genes-a manageable set for investigation of the roles of SMCs across the life cycle of the fungal model N. crassa-as well as gene sets coordinately expressed in 18 predicted SMCs during asexual and sexual growth under three nutritional and two temperature conditions. Divergent activity of SMCs between asexual and sexual development was reported. Of 126 SMC genes that we examined for knockout phenotypes, al-2 and al-3 exhibited phenotypes in asexual growth and conidiation, whereas os-5, poi-2, and pmd-1 exhibited phenotypes in sexual development. SMCs with annotated function in mating and crossing were actively regulated during the switch between asexual and sexual growth. Our discoveries call for attention to roles that SMCs may play in the regulatory switches controlling mode of development, as well as the ecological associations of those developmental stages that may influence expression of SMCs. IMPORTANCE Secondary metabolites (SMs) are low-molecular-weight compounds that often mediate interactions between fungi and their environments. Fungi enriched with SMs are of significant research interest to agriculture and medicine, especially from the aspects of pathogen ecology and environmental epidemiology. However, SM clusters (SMCs) that have been predicted by comparative genomics alone have typically been poorly defined and insufficiently functionally annotated. Therefore, we have investigated coordinate expression in SMCs in the model system N. crassa, and our results suggest that SMCs respond to environmental signals and to stress that are associated with development. This study examined SMC regulation at the level of RNA to integrate observations and knowledge of these genes in various growth and development conditions, supporting combining comparative genomics and inclusive transcriptomics to improve computational annotation of SMCs. Our findings call for detailed study of the function of SMCs during the asexual-sexual switch, a key, often-overlooked developmental stage.

Transcriptional Divergence Underpinning Sexual Development in the Fungal Class Sordariomycetes

Gene expression divergence through evolutionary processes is thought to be important for achieving programmed development in multicellular organisms. To test this premise in filamentous fungi, we investigated transcriptional profiles of 3,942 single-copy orthologous genes (SCOGs) in five related sordariomycete species that have morphologically diverged in the formation of their flask-shaped perithecia. We compared expression of the SCOGs to inferred gene expression levels of the most recent common ancestor of the five species, ranking genes from their largest increases to smallest increases in expression during perithecial development in each of the five species. We found that a large proportion of the genes that exhibited evolved increases in gene expression were important for normal perithecial development in Fusarium graminearum. Many of these genes were previously uncharacterized, encoding hypothetical proteins without any known functional protein domains. Interestingly, the developmental stages during which aberrant knockout phenotypes appeared largely coincided with the elevated expression of the deleted genes. In addition, we identified novel genes that affected normal perithecial development in Magnaporthe oryzae and Neurospora crassa, which were functionally and transcriptionally diverged from the orthologous counterparts in F. graminearum. Furthermore, comparative analysis of developmental transcriptomes and phylostratigraphic analysis suggested that genes encoding hypothetical proteins are generally young and transcriptionally divergent between related species. This study provides tangible evidence of shifts in gene expression that led to acquisition of novel function of orthologous genes in each lineage and demonstrates that several genes with hypothetical function are crucial for shaping multicellular fruiting bodies. IMPORTANCE The fungal class Sordariomycetes includes numerous important plant and animal pathogens. It also provides model systems for studying fungal fruiting body development, as its members develop fruiting bodies with a few well-characterized tissue types on common growth media and have rich genomic resources that enable comparative and functional analyses. To understand transcriptional divergence of key developmental genes between five related sordariomycete fungi, we performed targeted knockouts of genes inferred to have evolved significant upward shifts in expression. We found that many previously uncharacterized genes play indispensable roles at different stages of fruiting body development, which have undergone transcriptional activation in specific lineages. These novel genes are predicted to be phylogenetically young and tend to be involved in lineage- or species-specific function. Transcriptional activation of genes with unknown function seems to be more frequent than ever thought, which may be crucial for rapid adaption to changing environments for successful sexual reproduction.

Codon usage and protein length-dependent feedback from translation elongation regulates translation initiation and elongation speed

Essential cellular functions require efficient production of many large proteins but synthesis of large proteins encounters many obstacles in cells. Translational control is mostly known to be regulated at the initiation step. Whether translation elongation process can feedback to regulate initiation efficiency is unclear. Codon usage bias, a universal feature of all genomes, plays an important role in determining gene expression levels. Here, we discovered that there is a conserved but codon usage-dependent genome-wide negative correlation between protein abundance and CDS length. The codon usage effects on protein expression and ribosome flux on mRNAs are influenced by CDS length; optimal codon usage preferentially promotes production of large proteins. Translation of mRNAs with long CDS and non-optimal codon usage preferentially induces phosphorylation of initiation factor eIF2α, which inhibits translation initiation efficiency. Deletion of the eIF2α kinase CPC-3 (GCN2 homolog) in Neurospora preferentially up-regulates large proteins encoded by non-optimal codons. Surprisingly, CPC-3 also inhibits translation elongation rate in a codon usage and CDS length-dependent manner, resulting in slow elongation rates for long CDS mRNAs. Together, these results revealed a codon usage and CDS length-dependent feedback mechanism from translation elongation to regulate both translation initiation and elongation kinetics.

Transposon-mediated telomere destabilization: a driver of genome evolution in the blast fungus

The fungus Magnaporthe oryzae causes devastating diseases of crops, including rice and wheat, and in various grasses. Strains from ryegrasses have highly unstable chromosome ends that undergo frequent rearrangements, and this has been associated with the presence of retrotransposons (Magnaporthe oryzae Telomeric Retrotransposons-MoTeRs) inserted in the telomeres. The objective of the present study was to determine the mechanisms by which MoTeRs promote telomere instability. Targeted cloning, mapping, and sequencing of parental and novel telomeric restriction fragments (TRFs), along with MinION sequencing of genomic DNA allowed us to document the precise molecular alterations underlying 109 newly-formed TRFs. These included truncations of subterminal rDNA sequences; acquisition of MoTeR insertions by 'plain' telomeres; insertion of the MAGGY retrotransposons into MoTeR arrays; MoTeR-independent expansion and contraction of subtelomeric tandem repeats; and a variety of rearrangements initiated through breaks in interstitial telomere tracts that are generated during MoTeR integration. Overall, we estimate that alterations occurred in approximately sixty percent of chromosomes (one in three telomeres) analyzed. Most importantly, we describe an entirely new mechanism by which transposons can promote genomic alterations at exceptionally high frequencies, and in a manner that can promote genome evolution while minimizing collateral damage to overall chromosome architecture and function.

Genetic differentiation and phylogenetic potential of Ty3/Gypsy LTR retrotransposon markers in soil and plant pathogenic fungi

Genetic diversity studies are crucial for understanding the genetic structure and evolutionary dynamics of fungal species and communities. Fungal genomes are often reshaped by their repetitive components such as transposable elements. These elements are key players in genomic rearrangements and are ideal targets for genetic diversity and evolutionary studies. Herein, we used three Ty3/Gypsy long terminal repeat retrotransposons, Grasshopper, Maggy, and Pyret, for genetic differentiation and diversity in soil and plant pathogenic fungi, representing diverse species, order, and phyla. Pyret DNA markers showed the highest gene diversity and Shannon's information indices, followed by Maggy and Grasshopper. The observed high levels of multilocus polymorphism indicate the continuous mobility of these elements after their transfer in the new host. In conclusion, this study presents novel markers for genetic differentiation and evolutionary studies of fungi, and sheds light on the prevalence of gene acquisition phenomenon in field fungi.

Evolution and comparative genomics of the most common Trichoderma species

BACKGROUND

The growing importance of the ubiquitous fungal genus Trichoderma (Hypocreales, Ascomycota) requires understanding of its biology and evolution. Many Trichoderma species are used as biofertilizers and biofungicides and T. reesei is the model organism for industrial production of cellulolytic enzymes. In addition, some highly opportunistic species devastate mushroom farms and can become pathogens of humans. A comparative analysis of the first three whole genomes revealed mycoparasitism as the innate feature of Trichoderma. However, the evolution of these traits is not yet understood.

RESULTS

We selected 12 most commonly occurring Trichoderma species and studied the evolution of their genome sequences. Trichoderma evolved in the time of the Cretaceous-Palaeogene extinction event 66 (±15) mya, but the formation of extant sections (Longibrachiatum, Trichoderma) or clades (Harzianum/Virens) happened in Oligocene. The evolution of the Harzianum clade and section Trichoderma was accompanied by significant gene gain, but the ancestor of section Longibrachiatum experienced rapid gene loss. The highest number of genes gained encoded ankyrins, HET domain proteins and transcription factors. We also identified the Trichoderma core genome, completely curated its annotation, investigated several gene families in detail and compared the results to those of other fungi. Eighty percent of those genes for which a function could be predicted were also found in other fungi, but only 67% of those without a predictable function.

CONCLUSIONS

Our study presents a time scaled pattern of genome evolution in 12 Trichoderma species from three phylogenetically distant clades/sections and a comprehensive analysis of their genes. The data offer insights in the evolution of a mycoparasite towards a generalist.

Functional diversification accompanies gene family expansion of MED2 homologs in Candida albicans

Gene duplication facilitates functional diversification and provides greater phenotypic flexibility to an organism. Expanded gene families arise through repeated gene duplication but the extent of functional divergence that accompanies each paralogous gene is generally unexplored because of the difficulty in isolating the effects of single family members. The telomere-associated (TLO) gene family is a remarkable example of gene family expansion, with 14 members in the more pathogenic Candida albicans relative to two TLO genes in the closely-related species C. dubliniensis. TLO genes encode interchangeable Med2 subunits of the major transcriptional regulatory complex Mediator. To identify biological functions associated with each C. albicans TLO, expression of individual family members was regulated using a Tet-ON system and the strains were assessed across a range of phenotypes involved in growth and virulence traits. All TLOs affected multiple phenotypes and a single phenotype was often affected by multiple TLOs, including simple phenotypes such as cell aggregation and complex phenotypes such as virulence in a Galleria mellonella model of infection. No phenotype was regulated by all TLOs, suggesting neofunctionalization or subfunctionalization of ancestral properties among different family members. Importantly, regulation of three phenotypes could be mapped to individual polymorphic sites among the TLO genes, including an indel correlated with two phenotypes, growth in sucrose and macrophage killing. Different selective pressures have operated on the TLO sequence, with the 5’ conserved Med2 domain experiencing purifying selection and the gene/clade-specific 3’ end undergoing extensive positive selection that may contribute to the impact of individual TLOs on phenotypic variability. Therefore, expansion of the TLO gene family has conferred unique regulatory properties to each paralog such that it influences a range of phenotypes. We posit that the genetic diversity associated with this expansion contributed to C. albicans success as a commensal and opportunistic pathogen.

Gene duplication is a rapid mechanism to generate additional sequences for natural selection to act upon and confer greater organismal fitness. If additional copies of the gene are beneficial, this process may be repeated to produce an expanded gene family containing many copies of related sequences. Following duplication, individual gene family members may retain functions of the ancestral gene or acquire new functions through mutation. How functional diversification accompanies expansion into large gene families remains largely unexplored due to the difficulty in assessing individual genes in the presence of the remaining family members. Here, we addressed this question using an inducible promoter to regulate expression of individual genes of the TLO gene family in the commensal yeast and opportunistic pathogen Candida albicans, which encode components of a major transcriptional regulator. Induced expression of individual TLOs affected a wide range of phenotypes such that significant functional overlap occurred among TLO genes and most phenotypes were affected by more than one TLO. Induced expression of individual TLOs did not produce massive phenotypic effects in most cases, suggesting that functional overlap among TLO genes may buffer new mutations that arise. Specific sequence variants among the TLO genes correlated with certain phenotypes and these sequence variants did not necessarily correlate with sequence similarity across the entire gene. Therefore, individual TLO family members evolved specific functional roles following duplication that likely reflect a combination of inherited function and new mutation.

Codon usage is an important determinant of gene expression levels largely through its effects on transcription

Codon usage biases are found in all eukaryotic and prokaryotic genomes, and preferred codons are more frequently used in highly expressed genes. The effects of codon usage on gene expression were previously thought to be mainly mediated by its impacts on translation. Here, we show that codon usage strongly correlates with both protein and mRNA levels genome-wide in the filamentous fungus Neurospora Gene codon optimization also results in strong up-regulation of protein and RNA levels, suggesting that codon usage is an important determinant of gene expression. Surprisingly, we found that the impact of codon usage on gene expression results mainly from effects on transcription and is largely independent of mRNA translation and mRNA stability. Furthermore, we show that histone H3 lysine 9 trimethylation is one of the mechanisms responsible for the codon usage-mediated transcriptional silencing of some genes with nonoptimal codons. Together, these results uncovered an unexpected important role of codon usage in ORF sequences in determining transcription levels and suggest that codon biases are an adaptation of protein coding sequences to both transcription and translation machineries. Therefore, synonymous codons not only specify protein sequences and translation dynamics, but also help determine gene expression levels.

Analysis of Magnaporthe oryzae Genome Reveals a Fungal Effector, Which Is Able to Induce Resistance Response in Transgenic Rice Line Containing Resistance Gene, Pi54

Rice blast caused by Magnaporthe oryzae is one of the most important diseases of rice. Pi54, a rice gene that imparts resistance to M. oryzae isolates prevalent in India, was already cloned but its avirulent counterpart in the pathogen was not known. After decoding the whole genome of an avirulent isolate of M. oryzae, we predicted 11440 protein coding genes and then identified four candidate effector proteins which are exclusively expressed in the infectious structure, appresoria. In silico protein modeling followed by interaction analysis between Pi54 protein model and selected four candidate effector proteins models revealed that Mo-01947_9 protein model encoded by a gene located at chromosome 4 of M. oryzae, interacted best at the Leucine Rich Repeat domain of Pi54 protein model. Yeast-two-hybrid analysis showed that Mo-01947_9 protein physically interacts with Pi54 protein. Nicotiana benthamiana leaf infiltration assay confirmed induction of hypersensitive response in the presence of Pi54 gene in a heterologous system. Genetic complementation test also proved that Mo-01947_9 protein induces avirulence response in the pathogen in presence of Pi54 gene. Here, we report identification and cloning of a new fungal effector gene which interacts with blast resistance gene Pi54 in rice.

Nonoptimal codon usage influences protein structure in intrinsically disordered regions

Synonymous codons are not used with equal frequencies in most genomes. Codon usage has been proposed to play a role in regulating translation kinetics and co-translational protein folding. The relationship between codon usage and protein structures and the in vivo role of codon usage in eukaryotic protein folding is not clear. Here, we show that there is a strong codon usage bias in the filamentous fungus Neurospora. Importantly, we found genome-wide correlations between codon choices and predicted protein secondary structures: Nonoptimal codons are preferentially used in intrinsically disordered regions, and more optimal codons are used in structured domains. The functional importance of such correlations in vivo was confirmed by structure-based codon manipulation of codons in the Neurospora circadian clock gene frequency (frq). The codon optimization of the predicted disordered, but not well-structured regions of FRQ impairs clock function and altered FRQ structures. Furthermore, the correlations between codon usage and protein disorder tendency are conserved in other eukaryotes. Together, these results suggest that codon choices and protein structures co-evolve to ensure proper protein folding in eukaryotic organisms.

Real-Time Evolution of a Subtelomeric Gene Family in Candida albicans

Subtelomeric regions of the genome are notable for high rates of sequence evolution and rapid gene turnover. Evidence of subtelomeric evolution has relied heavily on comparisons of historical evolutionary patterns to infer trends and frequencies of these events. Here, we describe evolution of the subtelomeric TLO gene family in Candida albicans during laboratory passaging for over 4000 generations. C. albicans is a commensal and opportunistic pathogen of humans and the TLO gene family encodes a subunit of the Mediator complex that regulates transcription and affects a range of virulence factors. We identified 16 distinct subtelomeric recombination events that altered the TLO repertoire. Ectopic recombination between subtelomeres on different chromosome ends occurred approximately once per 5000 generations and was often followed by loss of heterozygosity, resulting in the complete loss of one TLO gene sequence with expansion of another. In one case, recombination within TLO genes produced a novel TLO gene sequence. TLO copy number changes were biased, with some TLOs preferentially being copied to novel chromosome arms and other TLO genes being frequently lost. The majority of these nonreciprocal recombination events occurred either within the 3′ end of the TLO coding sequence or within a conserved 50-bp sequence element centromere-proximal to TLO coding sequence. Thus, subtelomeric recombination is a rapid mechanism of generating genotypic diversity through alterations in the number and sequence of related gene family members.

Filamentous pathogen effector functions: of pathogens, hosts and microbiomes

Microorganisms play essential roles in almost every environment on earth. For instance, microbes decompose organic material, or establish symbiotic relationships that range from pathogenic to mutualistic. Symbiotic relationships have been particularly well studied for microbial plant pathogens and have emphasized the role of effectors; secreted molecules that support host colonization. Most effectors characterized thus far play roles in deregulation of host immunity. Arguably, however, pathogens not only deal with immune responses during host colonization, but also encounter other microbes including competitors, (myco)parasites and even potential co-operators. Thus, part of the effector catalog may target microbiome co-inhabitants rather than host physiology.

Discovering functions of unannotated genes from a transcriptome survey of wild fungal isolates

Most fungal genomes are poorly annotated, and many fungal traits of industrial and biomedical relevance are not well suited to classical genetic screens. Assigning genes to phenotypes on a genomic scale thus remains an urgent need in the field. We developed an approach to infer gene function from expression profiles of wild fungal isolates, and we applied our strategy to the filamentous fungus Neurospora crassa. Using transcriptome measurements in 70 strains from two well-defined clades of this microbe, we first identified 2,247 cases in which the expression of an unannotated gene rose and fell across N. crassa strains in parallel with the expression of well-characterized genes. We then used image analysis of hyphal morphologies, quantitative growth assays, and expression profiling to test the functions of four genes predicted from our population analyses. The results revealed two factors that influenced regulation of metabolism of nonpreferred carbon and nitrogen sources, a gene that governed hyphal architecture, and a gene that mediated amino acid starvation resistance. These findings validate the power of our population-transcriptomic approach for inference of novel gene function, and we suggest that this strategy will be of broad utility for genome-scale annotation in many fungal systems. IMPORTANCE Some fungal species cause deadly infections in humans or crop plants, and other fungi are workhorses of industrial chemistry, including the production of biofuels. Advances in medical and industrial mycology require an understanding of the genes that control fungal traits. We developed a method to infer functions of uncharacterized genes by observing correlated expression of their mRNAs with those of known genes across wild fungal isolates. We applied this strategy to a filamentous fungus and predicted functions for thousands of unknown genes. In four cases, we experimentally validated the predictions from our method, discovering novel genes involved in the metabolism of nutrient sources relevant for biofuel production, as well as colony morphology and starvation resistance. Our strategy is straightforward, inexpensive, and applicable for predicting gene function in many fungal species.

Genomics in Coccidioides: insights into evolution, ecology, and pathogenesis

Coccidioides immitis and C. posadasii, the causative agents of the mammalian disease coccidioidomycosis, are dimorphic fungal pathogens distributed throughout desert-like environments in North and South America. Coccidioides spp. are members of the Onygenales, a diverse group of pathogenic and nonpathogenic fungi. Recently, full genomes have been published for Coccidioides and a number of other Onygenales species. Phylogenomic comparisons and additional studies in Coccidioides population genomics and gene expression have shed light on the ecology and pathogenesis of Coccidioides and the other medically important species in this clade. Observed patterns of gene family expansion/contraction and evidence of gene flow have provided insight to the evolution of Coccidioides and greatly broadened our understanding of the diversity and sources of genetic variation found in fungi. In the future, expansion of the number of sequenced isolates from all populations will allow deeper insight into the evolutionary processes that have shaped this unique human pathogen. In addition, deep sequencing of isolates from a single Coccidioides population and pairing of those data with phenotype information on growth and pathogenicity for genome-wide association analysis will allow researchers to find genes responsible for any phenotype, virulence included, that shows variation in the population.

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.

Comparative transcriptomics reveals different strategies of Trichoderma mycoparasitism

BACKGROUND

Trichoderma is a genus of mycotrophic filamentous fungi (teleomorph Hypocrea) which possess a bright variety of biotrophic and saprotrophic lifestyles. The ability to parasitize and/or kill other fungi (mycoparasitism) is used in plant protection against soil-borne fungal diseases (biological control, or biocontrol). To investigate mechanisms of mycoparasitism, we compared the transcriptional responses of cosmopolitan opportunistic species and powerful biocontrol agents Trichoderma atroviride and T. virens with tropical ecologically restricted species T. reesei during confrontations with a plant pathogenic fungus Rhizoctonia solani.

RESULTS

The three Trichoderma spp. exhibited a strikingly different transcriptomic response already before physical contact with alien hyphae. T. atroviride expressed an array of genes involved in production of secondary metabolites, GH16 ß-glucanases, various proteases and small secreted cysteine rich proteins. T. virens, on the other hand, expressed mainly the genes for biosynthesis of gliotoxin, respective precursors and also glutathione, which is necessary for gliotoxin biosynthesis. In contrast, T. reesei increased the expression of genes encoding cellulases and hemicellulases, and of the genes involved in solute transport. The majority of differentially regulated genes were orthologues present in all three species or both in T. atroviride and T. virens, indicating that the regulation of expression of these genes is different in the three Trichoderma spp. The genes expressed in all three fungi exhibited a nonrandom genomic distribution, indicating a possibility for their regulation via chromatin modification.

CONCLUSION

This genome-wide expression study demonstrates that the initial Trichoderma mycotrophy has differentiated into several alternative ecological strategies ranging from parasitism to predation and saprotrophy. It provides first insights into the mechanisms of interactions between Trichoderma and other fungi that may be exploited for further development of biofungicides.

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.

The HaloTag: Improving Soluble Expression and Applications in Protein Functional Analysis

Technological and methodological advances have been critical for the rapidly evolving field of proteomics. The development of fusion tag systems is essential for purification and analysis of recombinant proteins. The HaloTag is a 34 KDa monomeric protein derived from a bacterial haloalkane dehalogenase. The majority of fusion tags in use today utilize a reversible binding interaction with a specific ligand. The HaloTag system is unique in that it forms a covalent linkage to its chloroalkane ligand. This linkage permits attachment of the HaloTag to a variety of functional reporters, which can be used to label and immobilize recombinant proteins. The success rate for HaloTag expression of soluble proteins is very high and comparable to maltose binding protein (MBP) tag. Furthermore, cleavage of the HaloTag does not result in protein insolubility that often is observed with the MBP tag. In the present report, we describe applications of the HaloTag system in our ongoing investigation of protein-protein interactions of the Y. pestis Type 3 secretion system on a custom protein microarray. We also describe the utilization of affinity purification/mass spectroscopy (AP/MS) to evaluate the utility of the Halo Tag system to characterize DNA binding activity and protein specificity.

Intrinsically disordered proteins aggregate at fungal cell-to-cell channels and regulate intercellular connectivity

Like animals and plants, multicellular fungi possess cell-to-cell channels (septal pores) that allow intercellular communication and transport. Here, using a combination of MS of Woronin body-associated proteins and a bioinformatics approach that identifies related proteins based on composition and character, we identify 17 septal pore-associated (SPA) proteins that localize to the septal pore in rings and pore-centered foci. SPA proteins are not homologous at the primary sequence level but share overall physical properties with intrinsically disordered proteins. Some SPA proteins form aggregates at the septal pore, and in vitro assembly assays suggest aggregation through a nonamyloidal mechanism involving mainly α-helical and disordered structures. SPA loss-of-function phenotypes include excessive septation, septal pore degeneration, and uncontrolled Woronin body activation. Together, our data identify the septal pore as a complex subcellular compartment and focal point for the assembly of unstructured proteins controlling diverse aspects of intercellular connectivity.

Telomere-targeted retrotransposons in the rice blast fungus Magnaporthe oryzae: agents of telomere instability

The fungus Magnaporthe oryzae is a serious pathogen of rice and other grasses. Telomeric restriction fragments in Magnaporthe isolates that infect perennial ryegrass (prg) are hotspots for genomic rearrangement and undergo frequent, spontaneous alterations during fungal culture. The telomeres of rice-infecting isolates are very stable by comparison. Sequencing of chromosome ends from a number of prg-infecting isolates revealed two related non-LTR retrotransposons (M. oryzae Telomeric Retrotransposons or MoTeRs) inserted in the telomere repeats. This contrasts with rice pathogen telomeres that are uninterrupted by other sequences. Genetic evidence indicates that the MoTeR elements are responsible for the observed instability. MoTeRs represent a new family of telomere-targeted transposons whose members are found exclusively in fungi.

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.

Correlation of gene expression and protein production rate - a system wide study

BACKGROUND

Growth rate is a major determinant of intracellular function. However its effects can only be properly dissected with technically demanding chemostat cultivations in which it can be controlled. Recent work on Saccharomyces cerevisiae chemostat cultivations provided the first analysis on genome wide effects of growth rate. In this work we study the filamentous fungus Trichoderma reesei (Hypocrea jecorina) that is an industrial protein production host known for its exceptional protein secretion capability. Interestingly, it exhibits a low growth rate protein production phenotype.

RESULTS

We have used transcriptomics and proteomics to study the effect of growth rate and cell density on protein production in chemostat cultivations of T. reesei. Use of chemostat allowed control of growth rate and exact estimation of the extracellular specific protein production rate (SPPR). We find that major biosynthetic activities are all negatively correlated with SPPR. We also find that expression of many genes of secreted proteins and secondary metabolism, as well as various lineage specific, mostly unknown genes are positively correlated with SPPR. Finally, we enumerate possible regulators and regulatory mechanisms, arising from the data, for this response.

CONCLUSIONS

Based on these results it appears that in low growth rate protein production energy is very efficiently used primarly for protein production. Also, we propose that flux through early glycolysis or the TCA cycle is a more fundamental determining factor than growth rate for low growth rate protein production and we propose a novel eukaryotic response to this i.e. the lineage specific response (LSR).

Causes and consequences of genome expansion in fungi

Fungi display a large diversity in genome size and complexity, variation that is often considered to be adaptive. But because nonadaptive processes can also have important consequences on the features of genomes, we investigated the relationship of genetic drift and genome size in the phylum Ascomycota using multiple indicators of genetic drift. We detected a complex relationship between genetic drift and genome size in fungi: genetic drift is associated with genome expansion on broad evolutionary timescales, as hypothesized for other eukaryotes; but within subphyla over smaller timescales, the opposite trend is observed. Moreover, fungi and bacteria display similar patterns of genome degradation that are associated with initial effects of genetic drift. We conclude that changes in genome size within Ascomycota have occurred using two different routes: large-scale genome expansions are catalyzed by increasing drift as predicted by the mutation-hazard model of genome evolution and small-scale modifications in genome size are independent of drift.

Expression of biomass-degrading enzymes is a major event during conidium development in Trichoderma reesei

The conidium plays a critical role in the life cycle of many filamentous fungi, being the primary means for survival under unfavorable conditions. To investigate the transcriptional changes taking place during the transition from growing hyphae to conidia in Trichoderma reesei, microarray experiments were performed. A total of 900 distinct genes were classified as differentially expressed, relative to their expression at time zero of conidiation, at least at one of the time points analyzed. The main functional categories (FunCat) overrepresented among the upregulated genes were those involving solute transport, metabolism, transcriptional regulation, secondary metabolite synthesis, lipases, proteases, and, particularly, cellulases and hemicellulases. Categories overrepresented among the downregulated genes were especially those associated with ribosomal and mitochondrial functions. The upregulation of cellulase and hemicellulase genes was dependent on the function of the positive transcriptional regulator XYR1, but XYR1 exerted no influence on conidiation itself. At least 20% of the significantly regulated genes were nonrandomly distributed within the T. reesei genome, suggesting an epigenetic component in the regulation of conidiation. The significant upregulation of cellulases and hemicellulases during this process, and thus cellulase and hemicellulase content in the spores of T. reesei, contributes to the hypothesis that the ability to hydrolyze plant biomass is a major trait of this fungus enabling it to break dormancy and reinitiate vegetative growth after a period of facing unfavorable conditions.

Targeted cloning of fungal telomeres

Telomeres are the sequences that form the ends of eukaryotic chromosomes and are essential structures that confer genome stability and guide chromosome behavior. In addition, the terminal regions of the chromosomes tend to house genes with predicted roles in ecological adaptation. Unfortunately, however, most fungal genome assemblies contain very few telomeres and, therefore, the identities of genes residing near the chromosome ends are often unknown. In an effort to develop a complete understanding of the organization and gene content of chromosome ends in a number of fungi, we developed efficient methods for the identification and targeted cloning of telomeres. This chapter describes the basic steps and shows exemplary results from the targeted cloning of Epichloë festucae telomeres.

Novel genes exhibit distinct patterns of function acquisition and network integration

BACKGROUND

Genes are created by a variety of evolutionary processes, some of which generate duplicate copies of an entire gene, while others rearrange pre-existing genetic elements or co-opt previously non-coding sequence to create genes with 'novel' sequences. These novel genes are thought to contribute to distinct phenotypes that distinguish organisms. The creation, evolution, and function of duplicated genes are well-studied; however, the genesis and early evolution of novel genes are not well-characterized. We developed a computational approach to investigate these issues by integrating genome-wide comparative phylogenetic analysis with functional and interaction data derived from small-scale and high-throughput experiments.

RESULTS

We examine the function and evolution of new genes in the yeast Saccharomyces cerevisiae. We observed significant differences in the functional attributes and interactions of genes created at different times and by different mechanisms. Novel genes are initially less integrated into cellular networks than duplicate genes, but they appear to gain functions and interactions more quickly than duplicates. Recently created duplicated genes show evidence of adapting existing functions to environmental changes, while young novel genes do not exhibit enrichment for any particular functions. Finally, we found a significant preference for genes to interact with other genes of similar age and origin.

CONCLUSIONS

Our results suggest a strong relationship between how and when genes are created and the roles they play in the cell. Overall, genes tend to become more integrated into the functional networks of the cell with time, but the dynamics of this process differ significantly between duplicate and novel genes.

Temporal and spatial regulation of gene expression during asexual development of Neurospora crassa

In this study we profiled spatial and temporal transcriptional changes during asexual sporulation in the filamentous fungus Neurospora crassa. Aerial tissue was separated from the mycelium to allow detection of genes specific to each tissue. We identified 2641 genes that were differentially expressed during development, which represents ∼25% of the predicted genes in the genome of this model fungus. On the basis of the distribution of functional annotations of 1102 of these genes, we identified gene expression patterns that define key physiological events during conidial development. Not surprisingly, genes encoding transcription factors, cell wall remodeling proteins, and proteins involved in signal transduction were differentially regulated during asexual development. Among the genes differentially expressed in aerial tissues the majority were unclassified and tended to be unique to ascomycete genomes. This finding is consistent with the view that these genes evolved for asexual development in the Pezizomycotina. Strains containing deletions of several differentially expressed genes encoding transcription factors exhibited asexual development-associated phenotypes. Gene expression patterns during asexual development suggested that cAMP signaling plays a critical role in the transition from aerial growth to proconidial chain formation. This observation prompted us to characterize a deletion of the gene encoding a high-affinity cAMP phosphodiesterase (NCU00478). NCU00478 was determined to be allelic to aconidiate-2, a previously identified genetic locus controlling conidiation.

Fungal evo-devo: organelles and multicellular complexity

Peroxisome-derived Woronin bodies of the Ascomycota phyla, and the endoplasmic reticulum (ER)-derived septal pore cap (SPC) of the Basidiomycota, are both fungal organelles that prevent cytoplasmic bleeding when multicellular hyphal filaments are wounded. Analysis of Woronin body constituent proteins suggests that these organelles evolved in part through gene duplication and co-opting of non-essential genes for new functions, indicating that new organelles can arise through typical evolutionary mechanisms. Interestingly, clades possessing the Woronin body and SPC also produce the largest and most complex multicellular fungal reproductive structures. Certain Woronin body and SPC mutants have defects in growth and development, suggesting functions beyond cellular wound healing. I argue that studying these specialized systems will help to reveal the basis for fungal diversity and provide general principles for co-evolution of organelles and multicellular complexity.

Conserved components, but distinct mechanisms for the placement and assembly of the cell division machinery in unicellular and filamentous ascomycetes

Cytokinesis is essential for cell proliferation, yet its molecular description is challenging, because >100 conserved proteins must be spatially and temporally co-ordinated. Despite the high importance of a tight co-ordination of cytokinesis with chromosome and organelle segregation, the mechanism for determining the cell division plane is one of the least conserved aspects of cytokinesis in eukaryotic cells. Budding and fission yeast have developed fundamentally distinct mechanisms to ensure proper nuclear segregation. The extent to which these pathways are conserved in multicellular fungi remains unknown. Recent progress indicates common components, but different mechanisms that are required for proper selection of the septation site in the different groups of Ascomycota. Cortical cues are used in yeast- and filament-forming species of the Saccharomycotina clade that are established at the incipient bud site or the hyphal tip respectively. In contrast, septum formation in the filament-forming Pezizomycotina species Aspergillus nidulans and Neurospora crassa seems more closely related to the fission yeast programme in that they may combine mitotic signals with a cell end-based marker system and Rho GTPase signalling. Thus, significant differences in the use and connection of conserved signalling modules become apparent that reflect the phylogenetic relationship of the analysed models.

Relaxed purifying selection and possibly high rate of adaptation in primate lineage-specific genes

Genes in the same organism vary in the time since their evolutionary origin. Without horizontal gene transfer, young genes are necessarily restricted to a few closely related species, whereas old genes can be broadly distributed across the phylogeny. It has been shown that young genes evolve faster than old genes; however, the evolutionary forces responsible for this pattern remain obscure. Here, we classify human–chimp protein-coding genes into different age classes, according to the breath of their phylogenetic distribution. We estimate the strength of purifying selection and the rate of adaptive selection for genes in different age classes. We find that older genes carry fewer and less frequent nonsynonymous single-nucleotide polymorphisms than younger genes suggesting that older genes experience a stronger purifying selection at the protein-coding level. We infer the distribution of fitness effects of new deleterious mutations and find that older genes have proportionally more slightly deleterious mutations and fewer nearly neutral mutations than younger genes. To investigate the role of adaptive selection of genes in different age classes, we determine the selection coefficient (γ = 2N_es) of genes using the MKPRF approach and estimate the ratio of the rate of adaptive nonsynonymous substitution to synonymous substitution (ω_A) using the DoFE method. Although the proportion of positively selected genes (γ > 0) is significantly higher in younger genes, we find no correlation between ω_A and gene age. Collectively, these results provide strong evidence that younger genes are subject to weaker purifying selection and more tenuous evidence that they also undergo adaptive evolution more frequently.

De novo assembly of a 40 Mb eukaryotic genome from short sequence reads: Sordaria macrospora, a model organism for fungal morphogenesis

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

Systematic overrepresentation of DNA termini and underrepresentation of subterminal regions among sequencing templates prepared from hydrodynamically sheared linear DNA molecules

Background

Analysis of fungal genome sequence assemblies reveals that telomeres are poorly represented even though telomeric reads tend to be superabundant. We surmised that the problem might lie in the DNA shearing conditions used to create clone libraries for genome sequencing.

Results

A shotgun strategy was used to sequence and assemble circular and linear cosmid DNAs sheared using conditions typical for a genome project. The DNA sheared in circular form assembled into a single sequence contig. However, the linearized cosmid produced an incomplete assembly because the two DNA termini, though greatly overrepresented in the clone library used for sequencing, were separated from neighboring sequences by gaps of ~1.4 and 1.8 kb. These gap sizes were reduced, but not eliminated, by shearing the linear cosmid into smaller fragments. Mapping of shearing breakpoints revealed a paucity of breaks in the subterminal regions of the linearized cosmid and also near chromosome ends of the fungus Neurospora crassa.

Conclusion

Together, our data indicate that the ends of linear DNA molecules are recalcitrant to hydrodynamic shearing. We propose that this causes DNA termini to be overrepresented in the resulting fragment population but ultimately prevents their incorporation into sequence assemblies.