A-to-I RNA editing is developmentally regulated and generally adaptive for sexual reproduction in Neurospora crassa

Comparative genomic analyses reveal evidence for adaptive A-to-I RNA editing in insect Adar gene

Although A-to-I RNA editing leads to similar effects to A-to-G DNA mutation, nonsynonymous RNA editing (recoding) is believed to confer its adaptiveness by 'epigenetically' regulating proteomic diversity in a temporospatial manner, avoiding the pleiotropic effect of genomic mutations. Recent discoveries on the evolutionary trajectory of Ser>Gly auto-editing site in insect Adar gene demonstrated a selective advantage to having an editable codon compared to uneditable ones. However, apart from pure observations, quantitative approaches for justifying the adaptiveness of individual RNA editing sites are still lacking. We performed a comparative genomic analysis on 113 Diptera species, focusing on the Adar Ser>Gly auto-recoding site in Drosophila. We only found one species having a derived Gly at the corresponding site, and this occurrence was significantly lower than genome-wide random expectation. This suggests that the Adar Ser>Gly site is unlikely to be genomically replaced with G during evolution, and thus indicating the advantage of editable status over hardwired genomic alleles. Similar trends were observed for the conserved Ile>Met recoding in gene Syt1. In the light of evolution, we established a comparative genomic approach for quantitatively justifying the adaptiveness of individual editing sites. Priority should be given to such adaptive editing sites in future functional studies.

Narrowing down the candidates of beneficial A-to-I RNA editing by comparing the recoding sites with uneditable counterparts

Adar-mediated adenosine-to-inosine (A-to-I) RNA editing mainly occurs in nucleus and diversifies the transcriptome in a flexible manner. It has been a challenging task to identify beneficial editing sites from the sea of total editing events. The functional Ser>Gly auto-recoding site in insect Adar gene has uneditable Ser codons in ancestral nodes, indicating the selective advantage to having an editable status. Here, we extended this case study to more metazoan species, and also looked for all Drosophila recoding events with potential uneditable synonymous codons. Interestingly, in D. melanogaster, the abundant nonsynonymous editing is enriched in the codons that have uneditable counterparts, but the Adar Ser>Gly case suggests that the editable orthologous codons in other species are not necessarily edited. The use of editable versus ancestral uneditable codon is a smart way to infer the selective advantage of RNA editing, and priority might be given to these editing sites for functional studies due to the feasibility to construct an uneditable allele. Our study proposes an idea to narrow down the candidates of beneficial recoding sites. Meanwhile, we stress that the matched transcriptomes are needed to verify the conservation of editing events during evolution.

DNA cytosine methyltransferases differentially regulate genome-wide hypermutation and interhomolog recombination in Trichoderma reesei meiosis

Trichoderma reesei is an economically important enzyme producer with several unique meiotic features. spo11, the initiator of meiotic double-strand breaks (DSBs) in most sexual eukaryotes, is dispensable for T. reesei meiosis. T. reesei lacks the meiosis-specific recombinase Dmc1. Rad51 and Sae2, the activator of the Mre11 endonuclease complex, promote DSB repair and chromosome synapsis in wild-type and spo11Δ meiosis. DNA methyltransferases (DNMTs) perform multiple tasks in meiosis. Three DNMT genes (rid1, dim2 and dimX) differentially regulate genome-wide cytosine methylation and C:G-to-T:A hypermutations in different chromosomal regions. We have identified two types of DSBs: type I DSBs require spo11 or rid1 for initiation, whereas type II DSBs do not rely on spo11 and rid1 for initiation. rid1 (but not dim2) is essential for Rad51-mediated DSB repair and normal meiosis. rid1 and rad51 exhibit a locus heterogeneity (LH) relationship, in which LH-associated proteins often regulate interconnectivity in protein interaction networks. This LH relationship can be suppressed by deleting dim2 in a haploid rid1Δ (but not rad51Δ) parental strain, indicating that dim2 and rid1 share a redundant function that acts earlier than rad51 during early meiosis. In conclusion, our studies provide the first evidence of the involvement of DNMTs during meiotic initiation and recombination.

Genetic and selective constraints on the optimization of gene product diversity

RNA and protein expressed from the same gene can have diverse isoforms due to various post-transcriptional and post-translational modifications. For the vast majority of alternative isoforms, It is unknown whether they are adaptive or simply biological noise. As we cannot experimentally probe the function of each isoform, we can ask whether the distribution of isoforms across genes and across species is consistent with expectations from different evolutionary processes. However, there is currently no theoretical framework that can generate such predictions. To address this, we developed a mathematical model where isoform abundances are determined collectively by cis-acting loci, trans-acting factors, gene expression levels, and isoform decay rates to predict isoform abundance distributions across species and genes in the face of mutation, genetic drift, and selection. We found that factors beyond selection, such as effective population size and the number of cis-acting loci, significantly influence evolutionary outcomes. Notably, suboptimal phenotypes are more likely to evolve when the population is small and/or when the number of cis-loci is large. We also explored scenarios where modification processes have both beneficial and detrimental effects, revealing a non-monotonic relationship between effective population size and optimization, demonstrating how opposing selection pressures on cis- and trans-acting loci can constrain the optimization of gene product diversity. As a demonstration of the power of our theory, we compared the expected distribution of A-to-I RNA editing levels in coleoids and found this to be largely consistent with non-adaptive explanations.

Hyphal editing of the conserved premature stop codon in CHE1 is stimulated by oxidative stress in Fusarium graminearum

Although genome-wide A-to-I editing mediated by adenosine-deaminase-acting-on-tRNA (ADAT) occurs during sexual reproduction in the presence of stage-specific cofactors, RNA editing is not known to occur during vegetative growth in filamentous fungi. Here we identified 33 A-to-I RNA editing events in vegetative hyphae of Fusarium graminearum and functionally characterized one conserved hyphal-editing site. Similar to ADAT-mediated editing during sexual reproduction, majority of hyphal-editing sites are in coding sequences and nonsynonymous, and have strong preference for U at -1 position and hairpin loops. Editing at TA437G, one of the hyphal-specific editing sites, is a premature stop codon correction (PSC) event that enables CHE1 gene to encode a full-length zinc fingertranscription factor. Manual annotations showed that this PSC site is conserved in CHE1 orthologs from closely-related Fusarium species. Whereas the che1 deletion and CHE1TAA (G438 to A) mutants had no detectable phenotype, the CHE1TGG (A437 to G) mutant was defective in hyphal growth, conidiation, sexual reproduction, and plant infection. However, the CHE1TGG mutant was increased in tolerance against oxidative stress and editing of TA437G in CHE1 was stimulated by H2O2 treatment in F. graminearum. These results indicate that fixation of the premature stop codon in CHE1 has a fitness cost on normal hyphal growth and reproduction but provides a benefit to tolerance against oxidative stress. Taken together, A-to-I editing events, although rare (not genome-wide), occur during vegetative growth and editing in CHE1 plays a role in response to oxidative stress in F. graminearum and likely in other fungal pathogens.

Functions and mechanisms of A-to-I RNA editing in filamentous ascomycetes

Although lack of ADAR (adenosine deaminase acting on RNA) orthologs, genome-wide A-to-I editing occurs specifically during sexual reproduction in a number of filamentous ascomycetes, including Fusarium graminearum and Neurospora crassa. Unlike ADAR-mediated editing in animals, fungal A-to-I editing has a strong preference for hairpin loops and U at -1 position, which leads to frequent editing of UAG and UAA stop codons. Majority of RNA editing events in fungi are in the coding region and cause amino acid changes. Some of these editing events have been experimentally characterized for providing heterozygote and adaptive advantages in F. graminearum. Recent studies showed that FgTad2 and FgTad3, 2 ADAT (adenosine deaminase acting on tRNA) enzymes that normally catalyze the editing of A34 in the anticodon of tRNA during vegetative growth mediate A-to-I mRNA editing during sexual reproduction. Stage specificity of RNA editing is conferred by stage-specific expression of short transcript isoforms of FgTAD2 and FgTAD3 as well as cofactors such as AME1 and FIP5 that facilitate the editing of mRNA in perithecia. Taken together, fungal A-to-I RNA editing during sexual reproduction is catalyzed by ADATs and it has the same sequence and structural preferences with editing of A34 in tRNA.

Unveiling the A-to-I mRNA editing machinery and its regulation and evolution in fungi

A-to-I mRNA editing in animals is mediated by ADARs, but the mechanism underlying sexual stage-specific A-to-I mRNA editing in fungi remains unknown. Here, we show that the eukaryotic tRNA-specific heterodimeric deaminase FgTad2-FgTad3 is responsible for A-to-I mRNA editing in Fusarium graminearum. This editing capacity relies on the interaction between FgTad3 and a sexual stage-specific protein called Ame1. Although Ame1 orthologs are widely distributed in fungi, the interaction originates in Sordariomycetes. We have identified key residues responsible for the FgTad3-Ame1 interaction. The expression and activity of FgTad2-FgTad3 are regulated through alternative promoters, alternative translation initiation, and post-translational modifications. Our study demonstrates that the FgTad2-FgTad3-Ame1 complex can efficiently edit mRNA in yeasts, bacteria, and human cells, with important implications for the development of base editors in therapy and agriculture. Overall, this study uncovers mechanisms, regulation, and evolution of RNA editing in fungi, highlighting the role of protein-protein interactions in modulating deaminase function.

Sexual stage-specific A-to-I mRNA editing is mediated by tRNA-editing enzymes in fungi

A-to-I RNA editing catalyzed by adenosine-deaminase-acting-on-RNA (ADARs) was assumed to be unique to metazoans because fungi and plants lack ADAR homologs. However, genome-wide messenger RNA (mRNA) editing was found to occur specifically during sexual reproduction in filamentous ascomycetes. Because systematic characterization of adenosine/cytosine deaminase genes has implicated the involvement of TAD2 and TAD3 orthologs in A-to-I editing, in this study, we used genetic and biochemical approaches to characterize the role of FgTAD2, an essential adenosine-deaminase-acting-on-tRNA (ADAT) gene, in mRNA editing in Fusarium graminearum. FgTAD2 had a sexual-stage-specific isoform and formed heterodimers with enzymatically inactive FgTAD3. Using a repeat-induced point (RIP) mutation approach, we identified 17 mutations in FgTAD2 that affected mRNA editing during sexual reproduction but had no effect on transfer RNA (tRNA) editing and vegetative growth. The functional importance of the H352Y and Q375*(nonsense) mutations in sexual reproduction and mRNA editing were confirmed by introducing specific point mutations into the endogenous FgTAD2 allele in the wild type. An in vitro assay was developed to show that FgTad2-His proteins purified from perithecia, but not from vegetative hyphae, had mRNA editing activities. Moreover, the H352Y mutation affected the enzymatic activity of FgTad2 to edit mRNA but had no effect on its ADAT activity. We also identified proteins co-purified with FgTad2-His by mass spectrometry analysis and found that two of them have the RNA recognition motif. Taken together, genetic and biochemical data from this study demonstrated that FgTad2, an ADAT, catalyzes A-to-I mRNA editing with the stage-specific isoform and cofactors during sexual reproduction in fungi.

The first A-to-I RNA editome of hemipteran species Coridius chinensis reveals overrepresented recoding and prevalent intron editing in early-diverging insects

BACKGROUND

Metazoan adenosine-to-inosine (A-to-I) RNA editing resembles A-to-G mutation and increases proteomic diversity in a temporal-spatial manner, allowing organisms adapting to changeable environment. The RNA editomes in many major animal clades remain unexplored, hampering the understanding on the evolution and adaptation of this essential post-transcriptional modification.

METHODS

We assembled the chromosome-level genome of Coridius chinensis belonging to Hemiptera, the fifth largest insect order where RNA editing has not been studied yet. We generated ten head RNA-Seq libraries with DNA-Seq from the matched individuals.

RESULTS

We identified thousands of high-confidence RNA editing sites in C. chinensis. Overrepresentation of nonsynonymous editing was observed, but conserved recoding across different orders was very rare. Under cold stress, the global editing efficiency was down-regulated and the general transcriptional processes were shut down. Nevertheless, we found an interesting site with "conserved editing but non-conserved recoding" in potassium channel Shab which was significantly up-regulated in cold, serving as a candidate functional site in response to temperature stress.

CONCLUSIONS

RNA editing in C. chinensis largely recodes the proteome. The first RNA editome in Hemiptera indicates independent origin of beneficial recoding during insect evolution, which advances our understanding on the evolution, conservation, and adaptation of RNA editing.

Phylogenetic and functional analyses of N6-methyladenosine RNA methylation factors in the wheat scab fungus Fusarium graminearum

In eukaryotes, N6-methyladenosine (m6A) RNA modification plays a crucial role in governing the fate of RNA molecules and has been linked to various developmental processes. However, the phyletic distribution and functions of genetic factors responsible for m6A modification remain largely unexplored in fungi. To get insights into the evolution of m6A machineries, we reconstructed global phylogenies of potential m6A writers, readers, and erasers in fungi. Substantial copy number variations were observed, ranging from up to five m6A writers in early-diverging fungi to a single copy in the subphylum Pezizomycotina, which primarily comprises filamentous fungi. To characterize m6A factors in a phytopathogenic fungus Fusarium graminearum, we generated knockout mutants lacking potential m6A factors including the sole m6A writer MTA1. However, the resulting knockouts did not exhibit any noticeable phenotypic changes during vegetative and sexual growth stages. As obtaining a homozygous knockout lacking MTA1 was likely hindered by its essential role, we generated MTA1-overexpressing strains (MTA1-OE). The MTA1-OE5 strain showed delayed conidial germination and reduced hyphal branching, suggesting its involvement during vegetative growth. Consistent with these findings, the expression levels of MTA1 and a potential m6A reader YTH1 were dramatically induced in germinating conidia, followed by the expression of potential m6A erasers at later vegetative stages. Several genes including transcription factors, transporters, and various enzymes were found to be significantly upregulated and downregulated in the MTA1-OE5 strain. Overall, our study highlights the functional importance of the m6A methylation during conidial germination in F. graminearum and provides a foundation for future investigations into m6A modification sites in filamentous fungi.IMPORTANCEN6-methyladenosine (m6A) RNA methylation is a reversible posttranscriptional modification that regulates RNA function and plays a crucial role in diverse developmental processes. This study addresses the knowledge gap regarding phyletic distribution and functions of m6A factors in fungi. The identification of copy number variations among fungal groups enriches our knowledge regarding the evolution of m6A machinery in fungi. Functional characterization of m6A factors in a phytopathogenic filamentous fungus Fusarium graminearum provides insights into the essential role of the m6A writer MTA1 in conidial germination and hyphal branching. The observed effects of overexpressing MTA1 on fungal growth and gene expression patterns of m6A factors throughout the life cycle of F. graminearum further underscore the importance of m6A modification in conidial germination. Overall, this study significantly advances our understanding of m6A modification in fungi, paving the way for future research into its roles in filamentous growth and potential applications in disease control.

Adaptive advantages of restorative RNA editing in fungi for resolving survival-reproduction trade-offs

RNA editing in various organisms commonly restores RNA sequences to their ancestral state, but its adaptive advantages are debated. In fungi, restorative editing corrects premature stop codons in pseudogenes specifically during sexual reproduction. We characterized 71 pseudogenes and their restorative editing in Fusarium graminearum, demonstrating that restorative editing of 16 pseudogenes is crucial for germ tissue development in fruiting bodies. Our results also revealed that the emergence of premature stop codons is facilitated by restorative editing and that premature stop codons corrected by restorative editing are selectively favored over ancestral amino acid codons. Furthermore, we found that ancestral versions of pseudogenes have antagonistic effects on reproduction and survival. Restorative editing eliminates the survival costs of reproduction caused by antagonistic pleiotropy and provides a selective advantage in fungi. Our findings highlight the importance of restorative editing in the evolution of fungal complex multicellularity and provide empirical evidence that restorative editing serves as an adaptive mechanism enabling the resolution of genetic trade-offs.

A non-pheromone GPCR is essential for meiosis and ascosporogenesis in the wheat scab fungus

Meiosis is essential for generating genetic diversity and sexual spores, but the regulation of meiosis and ascosporogenesis is not clear in filamentous fungi, in which dikaryotic and diploid cells formed inside fruiting bodies are not free living and independent of pheromones or pheromone receptors. In this study, Gia1, a non-pheromone GPCR (G protein-coupled receptor) with sexual-specific expression in Fusarium graminearum, is found to be essential for ascosporogenesis. The gia1 mutant was normal in perithecium development, crozier formation, and karyogamy but failed to undergo meiosis, which could be partially rescued by a dominant active mutation in GPA1 and activation of the Gpmk1 pathway. GIA1 orthologs have conserved functions in regulating meiosis and ascosporogenesis in Sordariomycetes. GIA1 has a paralog, GIP1, in F. graminearum and other Hypocreales species which is essential for perithecium formation. GIP1 differed from GIA1 in expression profiles and downstream signaling during sexual reproduction. Whereas the C-terminal tail and IR3 were important for intracellular signaling, the N-terminal region and EL3 of Gia1 were responsible for recognizing its ligand, which is likely a protein enriched in developing perithecia, particularly in the gia1 mutant. Taken together, these results showed that GIA1 encodes a non-pheromone GPCR that regulates the entry into meiosis and ascosporogenesis via the downstream Gpmk1 MAP kinase pathway in F. graminearum and other filamentous ascomycetes.

The Genome-Wide Characterization of Alternative Splicing and RNA Editing in the Development of Coprinopsis cinerea

Coprinopsis cinerea is one of the model species used in fungal developmental studies. This mushroom-forming Basidiomycetes fungus has several developmental destinies in response to changing environments, with dynamic developmental regulations of the organism. Although the gene expression in C. cinerea development has already been profiled broadly, previous studies have only focused on a specific stage or process of fungal development. A comprehensive perspective across different developmental paths is lacking, and a global view on the dynamic transcriptional regulations in the life cycle and the developmental paths is far from complete. In addition, knowledge on co- and post-transcriptional modifications in this fungus remains rare. In this study, we investigated the transcriptional changes and modifications in C. cinerea during the processes of spore germination, vegetative growth, oidiation, sclerotia formation, and fruiting body formation by inducing different developmental paths of the organism and profiling the transcriptomes using the high-throughput sequencing method. Transition in the identity and abundance of expressed genes drive the physiological and morphological alterations of the organism, including metabolism and multicellularity construction. Moreover, stage- and tissue-specific alternative splicing and RNA editing took place and functioned in C. cinerea. These modifications were negatively correlated to the conservation features of genes and could provide extra plasticity to the transcriptome during fungal development. We suggest that C. cinerea applies different molecular strategies in its developmental regulation, including shifts in expressed gene sets, diversifications of genetic information, and reversible diversifications of RNA molecules. Such features would increase the fungal adaptability in the rapidly changing environment, especially in the transition of developmental programs and the maintenance and balance of genetic and transcriptomic divergence. The multi-layer regulatory network of gene expression serves as the molecular basis of the functioning of developmental regulation.

A-to-I RNA editing shows dramatic up-regulation in osteosarcoma and broadly regulates tumor-related genes by altering microRNA target regions

A-to-I RNA editing is a prevalent type of RNA modification in animals. The dysregulation of RNA editing has led to multiple human cancers. However, the role of RNA editing has never been studied in osteosarcoma, a complex bone cancer with unknown molecular basis. We retrieved the RNA-sequencing data from 24 primary osteosarcoma patients and 3 healthy controls. We systematically profiled the RNA editomes in these samples and quantitatively identified reliable differential editing sites (DES) between osteosarcoma and normal samples. RNA editing efficiency is dramatically increased in osteosarcoma, presumably due to the significant up-regulation of editing enzymes ADAR1 and ADAR2. Up-regulated DES in osteosarcoma are enriched in 3'UTRs. Strikingly, such 3'UTR sites are further enriched in microRNA binding regions of gene EMP2 and other oncogenes, abolishing the microRNA suppression on target genes. Accordingly, the expression of these tumor-promoting genes is elevated in osteosarcoma. There might be an RNA editing-dependent pathway leading to osteosarcoma. We expanded our knowledge on the potential roles of RNA editing in oncogenesis. Based on these molecular features, our work is valuable for future prognosis and diagnosis of osteosarcoma.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

The Sordariomycetes: an expanding resource with Big Data for mining in evolutionary genomics and transcriptomics

Advances in genomics and transcriptomics accompanying the rapid accumulation of omics data have provided new tools that have transformed and expanded the traditional concepts of model fungi. Evolutionary genomics and transcriptomics have flourished with the use of classical and newer fungal models that facilitate the study of diverse topics encompassing fungal biology and development. Technological advances have also created the opportunity to obtain and mine large datasets. One such continuously growing dataset is that of the Sordariomycetes, which exhibit a richness of species, ecological diversity, economic importance, and a profound research history on amenable models. Currently, 3,574 species of this class have been sequenced, comprising nearly one-third of the available ascomycete genomes. Among these genomes, multiple representatives of the model genera Fusarium, Neurospora, and Trichoderma are present. In this review, we examine recently published studies and data on the Sordariomycetes that have contributed novel insights to the field of fungal evolution via integrative analyses of the genetic, pathogenic, and other biological characteristics of the fungi. Some of these studies applied ancestral state analysis of gene expression among divergent lineages to infer regulatory network models, identify key genetic elements in fungal sexual development, and investigate the regulation of conidial germination and secondary metabolism. Such multispecies investigations address challenges in the study of fungal evolutionary genomics derived from studies that are often based on limited model genomes and that primarily focus on the aspects of biology driven by knowledge drawn from a few model species. Rapidly accumulating information and expanding capabilities for systems biological analysis of Big Data are setting the stage for the expansion of the concept of model systems from unitary taxonomic species/genera to inclusive clusters of well-studied models that can facilitate both the in-depth study of specific lineages and also investigation of trait diversity across lineages. The Sordariomycetes class, in particular, offers abundant omics data and a large and active global research community. As such, the Sordariomycetes can form a core omics clade, providing a blueprint for the expansion of our knowledge of evolution at the genomic scale in the exciting era of Big Data and artificial intelligence, and serving as a reference for the future analysis of different taxonomic levels within the fungal kingdom.

A crucial role for dynamic expression of components encoding the negative arm of the circadian clock

In the Neurospora circadian system, the White Collar Complex (WCC) drives expression of the principal circadian negative arm component frequency (frq). FRQ interacts with FRH (FRQ-interacting RNA helicase) and CKI, forming a stable complex that represses its own expression by inhibiting WCC. In this study, a genetic screen identified a gene, designated as brd-8, that encodes a conserved auxiliary subunit of the NuA4 histone acetylation complex. Loss of brd-8 reduces H4 acetylation and RNA polymerase (Pol) II occupancy at frq and other known circadian genes, and leads to a long circadian period, delayed phase, and defective overt circadian output at some temperatures. In addition to strongly associating with the NuA4 histone acetyltransferase complex, BRD-8 is also found complexed with the transcription elongation regulator BYE-1. Expression of brd-8, bye-1, histone h2a.z, and several NuA4 subunits is controlled by the circadian clock, indicating that the molecular clock both regulates the basic chromatin status and is regulated by changes in chromatin. Taken together, our data identify auxiliary elements of the fungal NuA4 complex having homology to mammalian components, which along with conventional NuA4 subunits, are required for timely and dynamic frq expression and thereby a normal and persistent circadian rhythm.

A crucial role for dynamic expression of components encoding the negative arm of the circadian clock

In the Neurospora circadian system, the White Collar Complex (WCC) drives expression of the principal circadian negative arm component frequency ( frq ). FRQ interacts with FRH (FRQ-interacting helicase) and CK-1 forming a stable complex that represses its own expression by inhibiting WCC. In this study, a genetic screen identified a gene, designated as brd-8 , that encodes a conserved auxiliary subunit of the NuA4 histone acetylation complex. Loss of brd-8 reduces H4 acetylation and RNA polymerase (Pol) II occupancy at frq and other known circadian genes, and leads to a long circadian period, delayed phase, and defective overt circadian output at some temperatures. In addition to strongly associating with the NuA4 histone acetyltransferase complex, BRD-8 is also found complexed with the transcription elongation regulator BYE-1. Expression of brd-8, bye-1, histone hH2Az , and several NuA4 subunits is controlled by the circadian clock, indicating that the molecular clock both regulates the basic chromatin status and is regulated by changes in chromatin. Taken together, our data identify new auxiliary elements of the fungal NuA4 complex having homology to mammalian components, which along with conventional NuA4 subunits, are required for timely and dynamic frq expression and thereby a normal and persistent circadian rhythm.

Experimental evidence for the functional importance and adaptive advantage of A-to-I RNA editing in fungi

Adenosine-to-inosine (A-to-I) editing is the most prevalent type of RNA editing in animals, and it occurs in fungi specifically during sexual reproduction. However, it is debatable whether A-to-I RNA editing is adaptive. Deciphering the functional importance of individual editing sites is essential for the mechanistic understanding of the adaptive advantages of RNA editing. Here, by performing gene deletion for 17 genes with conserved missense editing (CME) sites and engineering underedited (ue) and overedited (oe) mutants for 10 CME sites using site-specific mutagenesis at the native locus in Fusarium graminearum, we demonstrated that two CME sites in CME5 and CME11 genes are functionally important for sexual reproduction. Although the overedited mutant was normal in sexual reproduction, the underedited mutant of CME5 had severe defects in ascus and ascospore formation like the deletion mutant, suggesting that the CME site of CME5 is co-opted for sexual development. The preediting residue of Cme5 is evolutionarily conserved across diverse classes of Ascomycota, while the postediting one is rarely hardwired into the genome, implying that editing at this site leads to higher fitness than a genomic A-to-G mutation. More importantly, mutants expressing only the underedited or the overedited allele of CME11 are defective in ascosporogenesis, while those expressing both alleles displayed normal phenotypes, indicating that concurrently expressing edited and unedited versions of Cme11 is more advantageous than either. Our study provides convincing experimental evidence for the long-suspected adaptive advantages of RNA editing in fungi and likely in animals.

The occurrence, characteristics, and adaptation of A-to-I RNA editing in bacteria: A review

A-to-I RNA editing is a very important post-transcriptional modification or co-transcriptional modification that creates isoforms and increases the diversity of proteins. In this process, adenosine (A) in RNA molecules is hydrolyzed and deaminated into inosine (I). It is well known that ADAR (adenosine deaminase acting on RNA)-dependent A-to-I mRNA editing is widespread in animals. Next, the discovery of A-to-I mRNA editing was mediated by TadA (tRNA-specific adenosine deaminase) in Escherichia coli which is ADAR-independent event. Previously, the editing event S128P on the flagellar structural protein FliC enhanced the bacterial tolerance to oxidative stress in Xoc. In addition, the editing events T408A on the enterobactin iron receptor protein XfeA act as switches by controlling the uptake of Fe³⁺ in response to the concentration of iron in the environment. Even though bacteria have fewer editing events, the great majority of those that are currently preserved have adaptive benefits. Interestingly, it was found that a TadA-independent A-to-I RNA editing event T408A occurred on xfeA, indicating that there may be other new enzymes that perform a function like TadA. Here, we review recent advances in the characteristics, functions, and adaptations of editing in bacteria.

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

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

Characterization of Neurospora crassa GH16, GH17, and GH72 gene families of cell wall crosslinking enzymes

GH16 chitin transferases, GH17 β-1,3-glucan transferases, and GH72 β-1,3-glucan/lichenin transferases are important fungal cell wall crosslinking enzymes. The Neurospora crassa genome encodes three genes from the GH17 gene family and five members in the GH16 subfamily 18 and 19 fungal chitin transferases. We created deletion mutants lacking all three GH17 genes and determined that they had wild type morphology and are more sensitive to cell wall perturbation reagents than the wild type. We also created deletion mutants lacking all five GH16 subfamily 18 and 19 genes and found that they had wild type morphology and are more sensitive to cell wall perturbation reagents than the wild type. We conclude that GH16 and GH17 enzymes play roles in cell wall biogenesis. In N. crassa, GH72 enzymes have been reported to be lichenin transferases, while in other fungi they have been shown to be the β-1,3-glucan transferases. Neurospora triple GH72 deletions give rise to a tight colonial morphology, sensitivity to cell wall perturbation reagents, and release of cell wall proteins into the medium. To ask if GH72 and GH17 enzymes might be redundant in N. crassa, we created sextuple mutants lacking the three GH72 genes and the three GH17 genes and found that they were indistinguishable from the GH72 triple mutant. We also found that a recombinant GH72 enzyme is able to form a lichenin-enzyme intermediate demonstrating that GH72 enzymes are lichenin transferases. The N. crassa GH72 enzymes are lichenin transferases and are not redundant with the GH17 β-1,3-glucan transferases.

Gene product diversity: adaptive or not?

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

Uncovering Cis-Regulatory Elements Important for A-to-I RNA Editing in Fusarium graminearum

Adenosine-to-inosine (A-to-I) RNA editing independent of adenosine deaminase acting on RNA (ADAR) enzymes was discovered in fungi recently, and shown to be crucial for sexual reproduction. However, the underlying mechanism for editing is unknown. Here, we combine genome-wide comparisons, proof-of-concept experiments, and machine learning to decipher cis-regulatory elements of A-to-I editing in Fusarium graminearum. We identified plenty of RNA primary sequences and secondary structural features that affect editing specificity and efficiency. Although hairpin loop structures contribute importantly to editing, unlike in animals, the primary sequences have more profound influences on editing than secondary structures. Nucleotide preferences at adjacent positions of editing sites are the most important features, especially preferences at the -1 position. Unexpectedly, besides the number of positions with preferred nucleotides, the combination of preferred nucleotides with depleted ones at different positions are also important for editing. Some cis-sequence features have distinct importance for editing specificity and efficiency. Machine learning models built from diverse sequence and secondary structural features can accurately predict genome-wide editing sites but not editing levels, indicating that the cis-regulatory principle of editing efficiency is more complex than that of editing specificity. Nevertheless, our model interpretation provides insights into the quantitative contribution of each feature to the prediction of both editing sites and levels. We found that efficient editing of FG3G34330 transcripts depended on the full-length RNA molecule, suggesting that additional RNA structural elements may also contribute to editing efficiency. Our work uncovers multidimensional cis-regulatory elements important for A-to-I RNA editing in F. graminearum, helping to elucidate the fungal editing mechanism. IMPORTANCE A-to-I RNA editing is a new epigenetic phenomenon that is crucial for sexual reproduction in fungi. Deciphering cis-regulatory elements of A-to-I RNA editing can help us elucidate the editing mechanism and develop a model that accurately predicts RNA editing. In this study, we discovered multiple RNA sequence and secondary structure features important for A-to-I editing in Fusarium graminearum. We also identified the cis-sequence features with distinct importance for editing specificity and efficiency. The potential importance of full-length RNA molecules for editing efficiency is also revealed. This study represents the first comprehensive investigation of the cis-regulatory principles of A-to-I RNA editing in fungi.

Sordaria macrospora Sterile Mutant pro34 Is Impaired in Respiratory Complex I Assembly

The formation of fruiting bodies is a highly regulated process that requires the coordinated formation of different cell types. By analyzing developmental mutants, many developmental factors have already been identified. Yet, a complete understanding of fruiting body formation is still lacking. In this study, we analyzed developmental mutant pro34 of the filamentous ascomycete Sordaria macrospora. Genome sequencing revealed a deletion in the pro34 gene encoding a putative mitochondrial complex I assembly factor homologous to Neurospora crassa CIA84. We show that PRO34 is required for fast vegetative growth, fruiting body and ascospore formation. The pro34 transcript undergoes adenosine to inosine editing, a process correlated with sexual development in fruiting body-forming ascomycetes. Fluorescence microscopy and western blot analysis showed that PRO34 is a mitochondrial protein, and blue-native PAGE revealed that the pro34 mutant lacks mitochondrial complex I. Inhibitor experiments revealed that pro34 respires via complexes III and IV, but also shows induction of alternative oxidase, a shunt pathway to bypass complexes III and IV. We discuss the hypothesis that alternative oxidase is induced to prevent retrograde electron transport to complex I intermediates, thereby protecting from oxidative stress.

FgCsn12 Is Involved in the Regulation of Ascosporogenesis in the Wheat Scab Fungus Fusarium graminearum

Fusarium head blight (FHB), caused by the fungal pathogen Fusarium graminearum, is a destructive disease worldwide. Ascospores are the primary inoculum of F. graminearum, and sexual reproduction is a critical step in its infection cycle. In this study, we characterized the functions of FgCsn12. Although the ortholog of FgCsn12 in budding yeast was reported to have a direct interaction with Csn5, which served as the core subunit of the COP9 signalosome, the interaction between FgCsn12 and FgCsn5 was not detected through the yeast two-hybrid assay. The deletion of FgCSN12 resulted in slight defects in the growth rate, conidial morphology, and pathogenicity. Instead of forming four-celled, uninucleate ascospores, the Fgcsn12 deletion mutant produced oval ascospores with only one or two cells and was significantly defective in ascospore discharge. The 3'UTR of FgCsn12 was dispensable for vegetative growth but essential for sexual reproductive functions. Compared with those of the wild type, 1204 genes and 2240 genes were up- and downregulated over twofold, respectively, in the Fgcsn12 mutant. Taken together, FgCsn12 demonstrated an important function in the regulation of ascosporogenesis in F. graminearum.

On the Mechanistic Basis of Killer Meiotic Drive in Fungi

Spore killers are specific genetic elements in fungi that kill sexual spores that do not contain them. A range of studies in the last few years have provided the long-awaited first insights into the molecular mechanistic aspects of spore killing in different fungal models, including both yeast-forming and filamentous Ascomycota. Here we describe these recent advances, focusing on the wtf system in the fission yeast Schizosaccharomyces pombe; the Sk spore killers of Neurospora species; and two spore-killer systems in Podospora anserina, Spok and [Het-s]. The spore killers appear thus far mechanistically unrelated. They can involve large genomic rearrangements but most often rely on the action of just a single gene. Data gathered so far show that the protein domains involved in the killing and resistance processes differ among the systems and are not homologous. The emerging picture sketched by these studies is thus one of great mechanistic and evolutionary diversity of elements that cheat during meiosis and are thereby preferentially inherited over sexual generations.

Identification of Candidate Genes Associated with Trichothecene Biosynthesis in Fusarium graminearum Species Complex Combined with Transcriptomic and Proteomic Analysis

The Fusarium graminearum species complex is the main causal agent of wheat head blight worldwide. Trichothecenes produced by the pathogen in infected grains have important food safety implications. Previously reported studies on trichothecene production have all focused on the conditions conducive to mycotoxin production, while the molecular mechanisms of trichothecene biosynthesis in Fusarium strains under normal or non-inducing conditions are still unclear. Here, a global analysis of the fungal gene expression of three strains using the Affymetrix Fusarium GeneChip under non-inducing conditions is reported. Differentially expressed genes were identified among strains with different trichothecene-production ability, and some novel genes associated with trichothecene biosynthesis were found by bioinformatics analysis. To verify the transcriptome results, proteomic analyses of the three strains were conducted under the same culture conditions. In total, 69 unique fungal proteins were identified in 77 protein spots. Combined with transcriptome and proteome analysis, 27 novel genes were predicted to be associated with trichothecene mycotoxin production. A protein, encoded by FGSG_01403, was found to be associated with trichothecene production via proteome analysis. Gene knock-out mutations of FGSG_01403 resulted in mutants with increased production of trichothecenes. Future functional analysis of the candidate genes identified in this study may reveal new insights into the negative regulation of trichothecene production in the Fusarium graminearum species complex.

Genome-Wide Identification and Characterization of RNA/DNA Differences Associated with Fusarium graminearum Infection in Wheat

RNA/DNA difference (RDD) is a post-transcriptional modification playing a crucial role in regulating diverse biological processes in eukaryotes. Although it has been extensively studied in plant chloroplast and mitochondria genomes, RDDs in plant nuclear genomes are not well studied at present. Here, we investigated the RDDs associated with fusarium head blight (FHB) through a novel method by comparing the RNA-seq data between Fusarium-infected and control samples of four wheat genotypes. A total of 187 high-confidence unique RDDs in 36 genes were identified, representing the first landscape of the FHB-responsive RDD in wheat. The majority (26) of these 36 RDD genes were correlated either positively or negatively with FHB levels. Effects of these RDDs on RNA and protein sequences have been identified, their editing frequency and the expression level of the corresponding genes provided, and the prediction of the effect on the minimum folding free energy of mRNA, miRNA binding, and colocation of RDDs with conserved domains presented. RDDs were predicted to induce modifications in the mRNA and protein structures of the corresponding genes. In two genes, TraesCS1B02G294300 and TraesCS3A02G263900, editing was predicted to enhance their affinity with tae-miR9661-5p and tae-miR9664-3p, respectively. To our knowledge, this study is the first report of the association between RDD and FHB in wheat; this will contribute to a better understanding of the molecular basis underlying FHB resistance, and potentially lead to novel strategies to improve wheat FHB resistance through epigenetic methods.

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.

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Genetic Networks That Govern Sexual Reproduction in the Pezizomycotina

Sexual development in filamentous fungi is a complex process that relies on the precise control of and interaction between a variety of genetic networks and pathways. The mating-type (MAT) genes are the master regulators of this process and typically act as transcription factors, which control the expression of genes involved at all stages of the sexual cycle. In many fungi, the sexual cycle typically begins when the mating pheromones of one mating type are recognized by a compatible partner, followed by physical interaction and fertilization. Subsequently, highly specialized sexual structures are formed, within which the sexual spores develop after rounds of meiosis and mitosis. These spores are then released and germinate, forming new individuals that initiate new cycles of growth. This review provides an overview of the known genetic networks and pathways that are involved in each major stage of the sexual cycle in filamentous ascomycete fungi.

RNA Modification and Its Implication in Plant Pathogenic Fungi

Interaction of a pathogen with its host plant requires both flexibility and rapid shift in gene expression programs in response to environmental cues associated with host cells. Recently, a growing volume of data on the diversity and ubiquity of internal RNA modifications has led to the realization that such modifications are highly dynamic and yet evolutionarily conserved system. This hints at these RNA modifications being an additional regulatory layer for genetic information, culminating in epitranscriptome concept. In plant pathogenic fungi, however, the presence and the biological roles of RNA modifications are largely unknown. Here we delineate types of RNA modifications, and provide examples demonstrating roles of such modifications in biology of filamentous fungi including fungal pathogens. We also discuss the possibility that RNA modification systems in fungal pathogens could be a prospective target for new agrochemicals.

Quantitative single molecule RNA-FISH and RNase-free cell wall digestion in Neurospora crassa

Single molecule RNA-FISH (smFISH) is a valuable tool for analysis of mRNA spatial patterning in fixed cells that is underutilized in filamentous fungi. A primary complication for fixed-cell imaging in filamentous fungi is the need for enzymatic cell wall permeabilization, which is compounded by considerable variability in cell wall composition between species. smFISH adds another layer of complexity due to a requirement for RNase free conditions. Here, we describe the cloning, expression, and purification of a chitinase suitable for supplementation of a commercially available RNase-free enzyme preparation for efficient permeabilization of the Neurospora cell wall. We further provide a method for smFISH in Neurospora which includes a tool for generating numerical data from images that can be used in downstream customized analysis protocols.

An introgressed gene causes meiotic drive in Neurospora sitophila

In order to survive, most organisms must deal with parasites. Such parasites can be other organisms or, sometimes, selfish genes found within the host genome itself. While much is known about parasitic organisms, the interaction with their hosts, and their ability to spread within and between species, much less is known about selfish genes. We here identify a selfish “spore killer” gene in the fungus Neurospora sitophila. The gene appears to have evolved within the genus but has entered the species through hybridization and introgression. We also show that the host can counteract the gene through RNA interference. These results shed light on the diversity of selfish genes in terms of origin, evolution, and host interactions.

Meiotic drive elements cause their own preferential transmission following meiosis. In fungi, this phenomenon takes the shape of spore killing, and in the filamentous ascomycete Neurospora sitophila, the Sk-1 spore killer element is found in many natural populations. In this study, we identify the gene responsible for spore killing in Sk-1 by generating both long- and short-read genomic data and by using these data to perform a genome-wide association test. We name this gene Spk-1. Through molecular dissection, we show that a single 405-nt-long open reading frame generates a product that both acts as a poison capable of killing sibling spores and as an antidote that rescues spores that produce it. By phylogenetic analysis, we demonstrate that the gene has likely been introgressed from the closely related species Neurospora hispaniola, and we identify three subclades of N. sitophila, one where Sk-1 is fixed, another where Sk-1 is absent, and a third where both killer and sensitive strain are found. Finally, we show that spore killing can be suppressed through an RNA interference-based genome defense pathway known as meiotic silencing by unpaired DNA. Spk-1 is not related to other known meiotic drive genes, and similar sequences are only found within Neurospora. These results shed light on the diversity of genes capable of causing meiotic drive, their origin and evolution, and their interaction with the host genome.

Stage-specific regulation of purine metabolism during infectious growth and sexual reproduction in Fusarium graminearum

Ascospores generated during sexual reproduction are the primary inoculum for the wheat scab fungus Fusarium graminearum. Purine metabolism is known to play important roles in fungal pathogens but its lifecycle stage-specific regulation is unclear. By characterizing the genes involved in purine de novo and salvage biosynthesis pathways, we showed that de novo syntheses of inosine, adenosine and guanosine monophosphates (IMP, AMP and GMP) are important for vegetative growth, sexual/asexual reproduction, and infectious growth, whereas purine salvage synthesis is dispensable for these stages in F. graminearum. Addition of GMP rescued the defects of the Fgimd1 mutant in vegetative growth and conidiation but not sexual reproduction, whereas addition of AMP rescued all of these defects of the Fgade12 mutant, suggesting that the function of de novo synthesis of GMP rather than AMP is distinct in sexual stages. Moreover, Acd1, an ortholog of AMP deaminase, is dispensable for growth but essential for ascosporogenesis and pathogenesis, suggesting that AMP catabolism has stage-specific functions during sexual reproduction and infectious growth. The expression of almost all the genes involved in de novo purine synthesis is downregulated during sexual reproduction and infectious growth relative to vegetative growth. This study revealed that F. graminearum has stage-specific regulation of purine metabolism during infectious growth and sexual reproduction.

Discovering RNA Editing Events in Fungi

RNA editing is an important posttranscriptional process that alters the genetic information of RNA encoded by genomic DNA. Adenosine-to-inosine (A-to-I) editing is the most prevalent type of RNA editing in animal kingdom, catalyzed by adenosine deaminases acting on RNA (ADARs). Recently, genome-wide A-to-I RNA editing is discovered in fungi, involving adenosine deamination mechanisms distinct from animals. Aiming to draw more attention to RNA editing in fungi, here we discuss the considerations for deep sequencing data preparation and the available various methods for detecting RNA editing, with a special emphasis on their usability for fungal RNA editing detection. We describe computational protocols for the identification of candidate RNA editing sites in fungi by using two software packages REDItools and RES-Scanner with RNA sequencing (RNA-Seq) and genomic DNA sequencing (DNA-Seq) data.

FgBUD14 is important for ascosporogenesis and involves both stage-specific alternative splicing and RNA editing during sexual reproduction

In wheat head blight fungus Fusarium graminearum, A-to-I RNA editing occurs specifically during sexual reproduction. Among the genes with premature stop codons (PSCs) that require RNA editing to encode full-length proteins, FgBUD14 also had alternative splicing events in perithecia. In this study, we characterized the functions of FgBUD14 and its post-transcriptional modifications during sexual reproduction. The Fgbud14 deletion mutant was slightly reduced in growth, conidiation and virulence. Although deletion of FgBUD14 had no effect on perithecium morphology, the Fgbud14 mutant was defective in crozier formation and ascus development. The FgBud14-GFP localized to the apex of ascogenous hyphae and croziers, which may be related to its functions during early sexual development. During vegetative growth and asexual reproduction, FgBud14-GFP localized to hyphal tips and both ends of conidia. Furthermore, mutations blocking the splicing of intron 2 that has the PSC site had no effect on the function of FgBUD14 during sexual reproduction but caused a similar defect in growth with Fgbud14 mutant. Expression of the non-editable FgBUD14^Intron2-TAA mutant allele also failed to complement the Fgbud14 mutant. Taken together, FgBUD14 plays important roles in ascus development, and both alternative splicing and RNA editing occur specifically to its transcripts during sexual reproduction in F. graminearum.

Structures, Biosynthesis, and Physiological Functions of (1,3;1,4)-β-D-Glucans

(1,3;1,4)-β-d-Glucans, also named as mixed-linkage glucans, are unbranched non-cellulosic polysaccharides containing both (1,3)- and (1,4)-β-linkages. The linkage ratio varies depending upon species origin and has a significant impact on the physicochemical properties of the (1,3;1,4)-β-d-glucans. (1,3;1,4)-β-d-Glucans were thought to be unique in the grasses family (Poaceae); however, evidence has shown that (1,3;1,4)-β-d-glucans are also synthesized in other taxa, including horsetail fern Equisetum, algae, lichens, and fungi, and more recently, bacteria. The enzyme involved in (1,3;1,4)-β-d-glucan biosynthesis has been well studied in grasses and cereal. However, how this enzyme is able to assemble the two different linkages remains a matter of debate. Additionally, the presence of (1,3;1,4)-β-d-glucan across the species evolutionarily distant from Poaceae but absence in some evolutionarily closely related species suggest that the synthesis is either highly conserved or has arisen twice as a result of convergent evolution. Here, we compare the structure of (1,3;1,4)-β-d-glucans present across various taxonomic groups and provide up-to-date information on how (1,3;1,4)-β-d-glucans are synthesized and their functions.

A-to-I RNA editing in bacteria increases pathogenicity and tolerance to oxidative stress

Adenosine-to-inosine (A-to-I) RNA editing is an important posttranscriptional event in eukaryotes; however, many features remain largely unexplored in prokaryotes. This study focuses on a serine-to-proline recoding event (S128P) that originated in the mRNA of fliC, which encodes a flagellar filament protein; the editing event was observed in RNA-seq samples exposed to oxidative stress. Using Sanger sequencing, we show that the S128P editing event is induced by H₂O₂. To investigate the in vivo interaction between RNAs and TadA, which is the principal enzyme for A-to-I editing, genome-wide RNA immunoprecipitation–coupled high-throughput sequencing (iRIP-Seq) analysis was performed using HA-tagged TadA from Xanthomonas oryzae pv. oryzicola. We found that TadA can bind to the mRNA of fliC and the binding motif is identical to that previously reported by Bar-Yaacov and colleagues. This editing event increased motility and enhanced tolerance to oxidative stress due to changes in flagellar filament structure, which was modelled in 3D and measured by TEM. The change in filament structure due to the S128P mutant increased biofilm formation, which was measured by the 3D laser scanning confocal microscopy. RNA-seq revealed that a gene cluster that contributes to siderophore biosynthesis and Fe³⁺ uptake was upregulated in S128P compared with WT. Based on intracellular levels of reactive oxygen species and an oxidative stress survival assay, we found that this gene cluster can contribute to the reduction of the Fenton reaction and increases biofilm formation and bacterial virulence. This oxidative stress response was also confirmed in Pseudomonas putida. Overall, our work demonstrates that A-to-I RNA editing plays a role in bacterial pathogenicity and adaptation to oxidative stress.

Adenosine-to-inosine (A-to-I) RNA editing is an important posttranscriptional event in eukaryotes that has only been recently documented in bacteria. In this study, we use multiple ‘omic’ approaches to show that A-to-I RNA editing can occur in fliC, a flagellar filament protein. We show that TadA, which encodes adenosine deaminase, can directly bind to mRNA of target genes through recognition of a GACG motif. This editing event changes a single amino acid residue from serine to proline in FliC, resulting in a structural change in the flagellar filament. This posttranscriptional editing event contributes to virulence and increases tolerance to oxidative stress by enhancing biofilm formation. Our results provide insight into a new mechanism that bacterial pathogens use to adapt to oxidative stress, which can also increase virulence.

Pig-specific RNA editing during early embryo development revealed by genome-wide comparisons

Posttranscriptional modification of mRNA sequences through RNA editing can increase transcriptome and proteome diversity in eukaryotes. Studies of fetal and adult tissues showed that adenosine‐to‐inosine RNA editing plays a crucial role in early human development, but there is a lack of global understanding of dynamic RNA editing during mammalian early embryonic development. Therefore, here we used RNA sequencing data from human, pig and mouse during early embryonic development to detect edited genes that may regulate stem cell pluripotency. We observed that although most of the RNA editing sites are located in intergenic, intron and UTR, a few editing sites are in coding regions and may result in nonsynonymous amino acid changes. Some editing sites are predicted to change the structure of a protein. We also report that HNF1A, TBX3, ACLY, ECI1 and ERDR1 are related to embryonic development and cell division.

Fungal RNA editing: who, when, and why?

Abstract

RNA editing occurs in all kingdoms of life and in various RNA species. The editing of nuclear protein-coding transcripts has long been known in metazoans, but was only recently detected in fungi. In contrast to many metazoan species, fungal editing sites occur mostly in coding regions, and therefore, fungal editing can change protein sequences and lead to modified or new functions of proteins. Indeed, mRNA editing is thought to be generally adaptive on fungi. Although RNA editing has been detected in both, Ascomycota and Basidiomycota, there seem to be considerable differences between these two classes of fungi concerning the types, the timing, and the purpose of editing. This review summarizes the characteristics of RNA editing in fungi and compares them to metazoan species and bacteria. In particular, it will review cellular processes affected by editing and speculate on the purpose of editing for fungal biology with a focus on the filamentous ascomycetes.

Key Points

• Fungi show various types of mRNA editing in nuclear transcripts.

• Fungal editing leads to proteome diversification.

• Filamentous ascomycetes may require editing for sexual sporulation.

• Wood-degrading basidiomycetes may use editing for adaptation to different substrates.

Sordaria macrospora: 25 years as a model organism for studying the molecular mechanisms of fruiting body development

Abstract

Fruiting bodies are among the most complex multicellular structures formed by fungi, and the molecular mechanisms that regulate their development are far from understood. However, studies with a number of fungal model organisms have started to shed light on this developmental process. One of these model organisms is Sordaria macrospora, a filamentous ascomycete from the order Sordariales. This fungus has been a genetic model organism since the 1950s, but its career as a model organism for molecular genetics really took off in the 1990s, when the establishment of a transformation protocol, a mutant collection, and an indexed cosmid library provided the methods and resources to start revealing the molecular mechanisms of fruiting body development. In the 2000s, “omics” methods were added to the S. macrospora tool box, and by 2020, 58 developmental genes have been identified in this fungus. This review gives a brief overview of major method developments for S. macrospora, and then focuses on recent results characterizing different processes involved in regulating development including several regulatory protein complexes, autophagy, transcriptional and chromatin regulation, and RNA editing.

Key points

•Sordaria macrospora is a model system for analyzing fungal fruiting body development.

•More than 100 developmental mutants are available for S. macrospora.

•More than 50 developmental genes have been characterized in S. macrospora.

The SR-protein FgSrp2 regulates vegetative growth, sexual reproduction and pre-mRNA processing by interacting with FgSrp1 in Fusarium graminearum

Serine/arginine (SR) proteins play significant roles in pre-mRNA splicing in eukaryotes. To investigate how gene expression influences fungal development and pathogenicity in Fusarium graminearum, a causal agent of Fusarium head blight (FHB) of wheat and barley, our previous study identified a SR protein FgSrp1 in F. graminearum, and showed that it is important for conidiation, plant infection and pre-mRNA processing. In this study, we identified another SR protein FgSrp2 in F. graminearum, which is orthologous to Schizosaccharomyces pombe Srp2. Our data showed that, whereas yeast Srp2 is essential for growth, deletion of FgSRP2 resulted in only slight defects in vegetative growth and perithecia melanization. FgSrp2 localized to the nucleus and both its N- and C-terminal regions were important for the localization to the nucleus. FgSrp2 interacted with FgSrp1 to form a complex in vivo. Double deletion of FgSRP1 and FgSRP2 revealed that they had overlapping functions in vegetative growth and sexual reproduction. RNA-seq analysis revealed that, although deletion of FgSRP2 alone had minimal effects, deletion of both FgSRP1 and FgSRP2 caused significant changes in gene transcription and RNA splicing. Overall, our results indicated that FgSrp2 regulates vegetative growth, sexual reproduction and pre-mRNA processing by interacting with FgSrp1.

Investigating RNA editing in deep transcriptome datasets with REDItools and REDIportal

RNA editing is a widespread post-transcriptional mechanism able to modify transcripts through insertions/deletions or base substitutions. It is prominent in mammals, in which millions of adenosines are deaminated to inosines by members of the ADAR family of enzymes. A-to-I RNA editing has a plethora of biological functions, but its detection in large-scale transcriptome datasets is still an unsolved computational task. To this aim, we developed REDItools, the first software package devoted to the RNA editing profiling in RNA-sequencing (RNAseq) data. It has been successfully used in human transcriptomes, proving the tissue and cell type specificity of RNA editing as well as its pervasive nature. Outcomes from large-scale REDItools analyses on human RNAseq data have been collected in our specialized REDIportal database, containing more than 4.5 million events. Here we describe in detail two bioinformatic procedures based on our computational resources, REDItools and REDIportal. In the first procedure, we outline a workflow to detect RNA editing in the human cell line NA12878, for which transcriptome and whole genome data are available. In the second procedure, we show how to identify dysregulated editing at specific recoding sites in post-mortem brain samples of Huntington disease donors. On a 64-bit computer running Linux with ≥32 GB of random-access memory (RAM), both procedures should take ~76 h, using 4 to 24 cores. Our protocols have been designed to investigate RNA editing in different organisms with available transcriptomic and/or genomic reads. Scripts to complete both procedures and a docker image are available at https://github.com/BioinfoUNIBA/REDItools.

Advanced strategy for metabolite exploration in filamentous fungi

Filamentous fungi comprise an abundance of gene clusters that encode high-value metabolites, whereas affluent gene clusters remain silent during laboratory conditions. Complex cellular metabolism further limits these metabolite yields. Therefore, diverse strategies such as genetic engineering and chemical mutagenesis have been developed to activate these cryptic pathways and improve metabolite productivity. However, lower efficiencies of gene modifications and screen tools delayed the above processes. To address the above issues, this review describes an alternative design-construction evaluation optimization (DCEO) approach. The DCEO tool provides theoretical and practical principles to identify potential pathways, modify endogenous pathways, integrate exogenous pathways, and exploit novel pathways for their diverse metabolites and desirable productivities. This DCEO method also offers different tactics to balance the cellular metabolisms, facilitate the genetic engineering, and exploit the scalable metabolites in filamentous fungi.

Integrative Activity of Mating Loci, Environmentally Responsive Genes, and Secondary Metabolism Pathways during Sexual Development of Chaetomium globosum

Fungal diversity has amazed evolutionary biologists for decades. One societally important aspect of this diversity manifests in traits that enable pathogenicity. The opportunistic pathogen Chaetomium globosum is well adapted to a high-humidity environment and produces numerous secondary metabolites that defend it from predation. Many of these chemicals can threaten human health. Understanding the phases of the C. globosum life cycle in which these products are made enables better control and even utilization of this fungus. Among its intriguing traits is that it both is self-fertile and lacks any means of propagule-based asexual reproduction. By profiling genome-wide gene expression across the process of sexual reproduction in C. globosum and comparing it to genome-wide gene expression in the model filamentous fungus N. crassa and other closely related fungi, we revealed associations among mating-type genes, sexual developmental genes, sexual incompatibility regulators, environmentally responsive genes, and secondary metabolic pathways.

The origins and maintenance of the rich fungal diversity have been longstanding issues in evolutionary biology. To investigate how differences in expression regulation contribute to divergences in development and ecology among closely related species, transcriptomes were compared between Chaetomium globosum, a homothallic pathogenic fungus thriving in highly humid ecologies, and Neurospora crassa, a heterothallic postfire saprotroph. Gene expression was quantified in perithecia at nine distinct morphological stages during nearly synchronous sexual development. Unlike N. crassa, expression of all mating loci in C. globosum was highly correlated. Key regulators of the initiation of sexual development in response to light stimuli—including orthologs of N. crassasub-1, sub-1-dependent gene NCU00309, and asl-1—showed regulatory dynamics matching between C. globosum and N. crassa. Among 24 secondary metabolism gene clusters in C. globosum, 11—including the cochliodones biosynthesis cluster—exhibited highly coordinated expression across perithecial development. C. globosum exhibited coordinately upregulated expression of histidine kinases in hyperosmotic response pathways—consistent with gene expression responses to high humidity we identified in fellow pathogen Fusarium graminearum. Bayesian networks indicated that gene interactions during sexual development have diverged in concert with the capacities both to reproduce asexually and to live a self-compatible versus self-incompatible life cycle, shifting the hierarchical roles of genes associated with conidiation and heterokaryon incompatibility in N. crassa and C. globosum. This divergence supports an evolutionary history of loss of conidiation due to unfavorable combinations of heterokaryon incompatibility in homothallic species.