The presence of GC-AG introns in Neurospora crassa and other euascomycetes determined from analyses of complete genomes: implications for automated gene prediction

Comment on Dimou et al. Profile of Membrane Cargo Trafficking Proteins and Transporters Expressed under N Source Derepressing Conditions in Aspergillus nidulans. J. Fungi 2021, 7, 560

Contrary to the opinion recently offered by Dimou et al., our previously published biochemical, subcellular and genetic data supported our contention that AN11127 corresponds to the A. nidulans gene encoding Sec12, which is the guanine nucleotide exchange factor (GEF) specific for SAR1. We add here additional bioinformatics evidence that fully disprove the otherwise negative evidence reported by Dimou et al., highlighting the dangers associated with the lax interpretation of genomic data. On the positive side, we establish guidelines for the identification of this key secretory gene in other species of Ascomycota and Basidiomycota, including species of medical and applied interest.

Rab GDP-dissociation inhibitor gdiA is an essential gene required for cell wall chitin deposition in Aspergillus niger

The cell wall is a distinctive feature of filamentous fungi, providing them with structural integrity and protection from both biotic and abiotic factors. Unlike plant cell walls, fungi rely on structurally strong hydrophobic chitin core for mechanical strength together with alpha- and beta-glucans, galactomannans and glycoproteins. Cell wall stress conditions are known to alter the cell wall through the signaling cascade of the cell wall integrity (CWI) pathway and can result in increased cell wall chitin deposition. A previously isolated set of Aspergillus niger cell wall mutants was screened for increased cell wall chitin deposition. UV-mutant RD15.8#16 was found to contain approximately 60% more cell wall chitin than the wild type. In addition to the chitin phenotype, RD15.8#16 exhibits a compact colony morphology and increased sensitivity towards SDS. RD15.8#16 was subjected to classical genetic approach for identification of the underlying causative mutation, using co-segregation analysis and SNP genotyping. Genome sequencing of RD15.8#16 revealed eight SNPs in open reading frames (ORF) which were individually checked for co-segregation with the associated phenotypes, and showed the potential relevance of two genes located on chromosome IV. In situ re-creation of these ORF-located SNPs in a wild type background, using CRISPR/Cas9 genome editing, showed the importance Rab GTPase dissociation inhibitor A (gdiA) for the phenotypes of RD15.8#16. An alteration in the 5' donor splice site of gdiA reduced pre-mRNA splicing efficiency, causing aberrant cell wall assembly and increased chitin levels, whereas gene disruption attempts showed that a full gene deletion of gdiA is lethal.

Complementation of CTB7 in the Maize Pathogen Cercospora zeina Overcomes the Lack of In Vitro Cercosporin Production

Gray leaf spot (GLS), caused by the sibling species Cercospora zeina or Cercospora zeae-maydis, is cited as one of the most important diseases threatening global maize production. C. zeina fails to produce cercosporin in vitro and, in most cases, causes large coalescing lesions during maize infection, a symptom generally absent from cercosporin-deficient mutants in other Cercospora spp. Here, we describe the C. zeina cercosporin toxin biosynthetic (CTB) gene cluster. The oxidoreductase gene CTB7 contained several insertions and deletions as compared with the C. zeae-maydis ortholog. We set out to determine whether complementing the defective CTB7 gene with the full-length gene from C. zeae-maydis could confer in vitro cercosporin production. C. zeina transformants containing C. zeae-maydis CTB7 were generated by Agrobacterium tumefaciens-mediated transformation and were evaluated for in vitro cercosporin production. When grown on nitrogen-limited medium in the light-conditions conducive to cercosporin production in other Cercospora spp.-one transformant accumulated a red pigment that was confirmed to be cercosporin by the KOH assay, thin-layer chromatography, and ultra performance liquid chromatography-quadrupole-time-of-flight mass spectrometry. Our results indicated that C. zeina has a defective CTB7, but all other necessary machinery required for synthesizing cercosporin-like molecules and, thus, C. zeina may produce a structural variant of cercosporin during maize infection.

Evolution of Fungal U3 snoRNAs: Structural Variation and Introns

The U3 small nucleolar RNA (snoRNA) is an essential player in the initial steps of ribosomal RNA biogenesis which is ubiquitously present in Eukarya. It is exceptional among the small nucleolar RNAs in its size, the presence of multiple conserved sequence boxes, a highly conserved secondary structure core, its biogenesis as an independent gene transcribed by polymerase III, and its involvement in pre-rRNA cleavage rather than chemical modification. Fungal U3 snoRNAs share many features with their sisters from other eukaryotic kingdoms but differ from them in particular in their 5’ regions, which in fungi has a distinctive consensus structure and often harbours introns. Here we report on a comprehensive homology search and detailed analysis of the evolution of sequence and secondary structure features covering the entire kingdom Fungi.

U6 snRNA intron insertion occurred multiple times during fungi evolution

U6 small nuclear RNAs are part of the splicing machinery. They exhibit several unique features setting them appart from other snRNAs. Reports of introns in structured non-coding RNAs have been very rare. U6 genes, however, were found to be interrupted by an intron in several Schizosaccharomyces species and in 2 Basidiomycota. We conducted a homology search across 147 currently available fungal genome and identified the U6 genes in all but 2 of them. A detailed comparison of their sequences and predicted secondary structures showed that intron insertion events in the U6 snRNA were much more common in the fungal lineage than previously thought. Their positional distribution across the entire mature snRNA strongly suggests a large number of independent events. All the intron sequences reported here show canonical splice site and branch site motifs indicating that they require the splicesomal pathway for their removal.

Basidioascus undulatus: genome, origins, and sexuality

Basidioascus undulatus is a soil basidiomycete belonging to the order Geminibasidiales. The taxonomic status of the order was unclear as originally it was only tentatively classified in the class Wallemiomycetes. The fungi in Geminibasidiales have an ambiguously defined sexual cycle. In this study, we sequenced the genome of B. undulatus to gain insights into its sexuality and evolutionary origins. The assembled genome draft was approximately 32 Mb in size, had a median nucleotide coverage of 24X, and contained 6123 predicted genes. Previous morphological descriptions of B. undulatus relied on interpretation of putative sexual structures. In this study, nuclear staining and confocal microscopy showed meiosis occurring in basidia and genome analysis confirmed the existence of genes involved in meiosis and mating. Using 35 protein-coding genes extracted from genomic information, phylogenomic and molecular dating analyses confirmed that B. undulatus indeed belongs to a lineage distantly related to Wallemia while retaining a basal position in Agaricomycotina. These results, combined with differences in septal pore morphology, led us to move the order Geminibasidiales out of the Wallemiomycetes and into the new class Geminibasidiomycetes cl. nov. Finally, the concept of Agaricomycotina is emended to include both Wallemiomycetes and Geminibasidiomycetes.

Improving the gene structure annotation of the apicomplexan parasite Neospora caninum fulfils a vital requirement towards an in silico-derived vaccine

Neospora caninum is an apicomplexan parasite which can cause abortion in cattle, instigating major economic burden. Vaccination has been proposed as the most cost-effective control measure to alleviate this burden. Consequently the overriding aspiration for N. caninum research is the identification and subsequent evaluation of vaccine candidates in animal models. To save time, cost and effort, it is now feasible to use an in silico approach for vaccine candidate prediction. Precise protein sequences, derived from the correct open reading frame, are paramount and arguably the most important factor determining the success or failure of this approach. The challenge is that publicly available N. caninum sequences are mostly derived from gene predictions. Annotated inaccuracies can lead to erroneously predicted vaccine candidates by bioinformatics programs. This study evaluates the current N. caninum annotation for potential inaccuracies. Comparisons with annotation from a closely related pathogen, Toxoplasma gondii, are also made to distinguish patterns of inconsistency. More importantly, a mRNA sequencing (RNA-Seq) experiment is used to validate the annotation. Potential discrepancies originating from a questionable start codon context and exon boundaries were identified in 1943 protein coding sequences. We conclude, where experimental data were available, that the majority of N. caninum gene sequences were reliably predicted. Nevertheless, almost 28% of genes were identified as questionable. Given the limitations of RNA-Seq, the intention of this study was not to replace the existing annotation but to support or oppose particular aspects of it. Ideally, many studies aimed at improving the annotation are required to build a consensus. We believe this study, in providing a new resource on gene structure and annotation, is a worthy contributor to this endeavour.

Fusarium sambucinum astA gene expressed during potato infection is a functional orthologue of Aspergillus nidulans astA

Sulfate assimilation plays a vital role in prototrophic organisms. Orthologues of the alternative sulfate transporter (AstA) gene from Aspergillus nidulans were identified in the fungal plant pathogens Fusarium sambucinum and Fusarium graminearum. By physiological and biochemical analyses, the AstA orthologues were determined to be able to uptake sulfate from the environment. Similarly to astA in A. nidulans, the FsastA gene was found to be regulated by sulfur metabolite repression (SMR) in a sulfur-dependent manner. In contrast, the FgastA transcript was undetectable, however, when the FgastA gene was expressed heterologously in A. nidulans, the translated FgAstA protein acted as a sulfate transporter. Interestingly, F. sambucinum astA expression was remarkably augmented in infected potato tubers, despite the presence abundant sulfate and was found not to be correlated with plant resistance.

BcAtf1, a global regulator, controls various differentiation processes and phytotoxin production in Botrytis cinerea

Atf1-homologous basic region leucine zipper (bZIP) transcription factors are known to act downstream of the stress-activated mitogen-activated protein kinase (SAPK) cascade in mammals, as well as in several fungi; they regulate the transcription of genes involved in the general stress response. Functional analyses of BcAtf1 in Botrytis cinerea show that it is also connected to the SAPK BcSak1, as it shares several stress response target genes. However, Δbcatf1 mutants are not hypersensitive to osmotic or oxidative stress, as are Δbcsak1 mutants. Both BcSak1 and BcAtf1 are regulators of differentiation, but their roles in these processes are almost inverse as, in contrast with Δbcsak1, Δbcatf1 mutants are significantly impaired in conidia production and do not differentiate any sclerotia. They show extremely vigorous growth in axenic culture, with a thick layer of aerial hyphae and a marked increase in colonization efficiency on different host plants and tissues. In addition, the sensitivity to cell wall-interfering agents is increased strongly. Microarray analyses demonstrate that the loss of BcAtf1 leads to extensive transcriptional changes: apart from stress response genes, the expression of a broad set of genes, probably involved in primary metabolism, cell wall synthesis and development, is affected by BcAtf1. Unexpectedly, BcAtf1 also controls secondary metabolism: the mutant contains significantly elevated levels of phytotoxins. These data indicate that BcAtf1 controls a diversity of cellular processes and has broad regulatory functions.

Short-read sequencing for genomic analysis of the brown rot fungus Fibroporia radiculosa

The feasibility of short-read sequencing for genomic analysis was demonstrated for Fibroporia radiculosa, a copper-tolerant fungus that causes brown rot decay of wood. The effect of read quality on genomic assembly was assessed by filtering Illumina GAIIx reads from a single run of a paired-end library (75-nucleotide read length and 300-bp fragment size) at three different stringency levels and then assembling each data set with Velvet. A simple approach was devised to determine which filter stringency was "best." Venn diagrams identified the regions containing reads that were used in an assembly but were of a low-enough quality to be removed by a filter. By plotting base quality histograms of reads in this region, we judged whether a filter was too stringent or not stringent enough. Our best assembly had a genome size of 33.6 Mb, an N50 of 65.8 kb for a k-mer of 51, and a maximum contig length of 347 kb. Using GeneMark, 9,262 genes were predicted. TargetP and SignalP analyses showed that among the 1,213 genes with secreted products, 986 had motifs for signal peptides and 227 had motifs for signal anchors. Blast2GO analysis provided functional annotation for 5,407 genes. We identified 29 genes with putative roles in copper tolerance and 73 genes for lignocellulose degradation. A search for homologs of these 102 genes showed that F. radiculosa exhibited more similarity to Postia placenta than Serpula lacrymans. Notable differences were found, however, and their involvements in copper tolerance and wood decay are discussed.

The aspartic proteinase family of three Phytophthora species

BACKGROUND

Phytophthora species are oomycete plant pathogens with such major social and economic impact that genome sequences have been determined for Phytophthora infestans, P. sojae and P. ramorum. Pepsin-like aspartic proteinases (APs) are produced in a wide variety of species (from bacteria to humans) and contain conserved motifs and landmark residues. APs fulfil critical roles in infectious organisms and their host cells. Annotation of Phytophthora APs would provide invaluable information for studies into their roles in the physiology of Phytophthora species and interactions with their hosts.

RESULTS

Genomes of Phytophthora infestans, P. sojae and P. ramorum contain 11-12 genes encoding APs. Nine of the original gene models in the P. infestans database and several in P. sojae and P. ramorum (three and four, respectively) were erroneous. Gene models were corrected on the basis of EST data, consistent positioning of introns between orthologues and conservation of hallmark motifs. Phylogenetic analysis resolved the Phytophthora APs into 5 clades. Of the 12 sub-families, several contained an unconventional architecture, as they either lacked a signal peptide or a propart region. Remarkably, almost all APs are predicted to be membrane-bound.

CONCLUSIONS

One of the twelve Phytophthora APs is an unprecedented fusion protein with a putative G-protein coupled receptor as the C-terminal partner. The others appear to be related to well-documented enzymes from other species, including a vacuolar enzyme that is encoded in every fungal genome sequenced to date. Unexpectedly, however, the oomycetes were found to have both active and probably-inactive forms of an AP similar to vertebrate BACE, the enzyme responsible for initiating the processing cascade that generates the Aβ peptide central to Alzheimer's Disease. The oomycetes also encode enzymes similar to plasmepsin V, a membrane-bound AP that cleaves effector proteins of the malaria parasite Plasmodium falciparum during their translocation into the host red blood cell. Since the translocation of Phytophthora effector proteins is currently a topic of intense research activity, the identification in Phytophthora of potential functional homologues of plasmepsin V would appear worthy of investigation. Indeed, elucidation of the physiological roles of the APs identified here offers areas for future study. The significant revision of gene models and detailed annotation presented here should significantly facilitate experimental design.

Oxalate decarboxylase of the white-rot fungus Dichomitus squalens demonstrates a novel enzyme primary structure and non-induced expression on wood and in liquid cultures

Oxalate decarboxylase (ODC) catalyses the conversion of oxalic acid to formic acid and CO2 in bacteria and fungi. In wood-decaying fungi the enzyme has been linked to the regulation of intra- and extracellular quantities of oxalic acid, which is one of the key components in biological decomposition of wood. ODC enzymes are biotechnologically interesting for their potential in diagnostics, agriculture and environmental applications, e.g. removal of oxalic acid from industrial wastewaters. We identified a novel ODC in mycelial extracts of two wild-type isolates of Dichomitus squalens, and cloned the corresponding Ds-odc gene. The primary structure of the Ds-ODC protein contains two conserved Mn-binding cupin motifs, but at the N-terminus, a unique, approximately 60 aa alanine-serine-rich region is found. Real-time quantitative RT-PCR analysis confirmed gene expression when the fungus was cultivated on wood and in liquid medium. However, addition of oxalic acid in liquid cultures caused no increase in transcript amounts, thereby indicating a constitutive rather than inducible expression of Ds-odc. The detected stimulation of ODC activity by oxalic acid is more likely due to enzyme activation than to transcriptional upregulation of the Ds-odc gene. Our results support involvement of ODC in primary rather than secondary metabolism in fungi.

The Aspergillus nidulans stress response transcription factor StzA is ascomycete-specific and shows species-specific polymorphisms in the C-terminal region

Orthologues of the Aspergillus nidulans gene stzA were identified and characterised in an additional 19 fungi. These orthologues were restricted to, and found within all the Pezizomycotina subphyla of the Ascomycota, for which data are available, but not the Saccharomycotina or Taphrinomycotina subphyla. Intron analysis indicated that both intron loss and gain have occurred in this gene. The orthologous proteins demonstrate considerable size variation (between 663 and 897 amino acids); with almost all this variability accounted for by a hyper-variable region that is carboxy terminal to the zinc finger region. The Hypocrea jecorina orthologue (ACE1) has the binding site 5'AGGCA. There is evidence of competition, or interaction, between the ACE1/StzA and AreA binding sites in promoters of stzA and its orthologues, as well as genes involved in the metabolism of amino acids. The A. nidulans and A. fumigatus cpcA promoters have seven potential ACE1/StzA binding sites, six of which are highly conserved in position. Two very closely positioned sites are conserved across 14 of the 19 fungi analysed. Potential CpcA binding sites (5'TGAC/GTCA) have been identified between -50 and -170bp of the ATG start in the promoters of 16 of the stzA orthologues.

Purifying selection and birth-and-death evolution in the class II hydrophobin gene families of the ascomycete Trichoderma/Hypocrea

BACKGROUND

Hydrophobins are proteins containing eight conserved cysteine residues that occur uniquely in mycelial fungi. Their main function is to confer hydrophobicity to fungal surfaces in contact with air or during attachment of hyphae to hydrophobic surfaces of hosts, symbiotic partners or themselves resulting in morphogenetic signals. Based on their hydropathy patterns and solubility characteristics, hydrophobins are divided into two classes (I and II), the latter being found only in ascomycetes.

RESULTS

We have investigated the mechanisms driving the evolution of the class II hydrophobins in nine species of the mycoparasitic ascomycetous genus Trichoderma/Hypocrea, using three draft sequenced genomes (H. jecorina = T. reesei, H. atroviridis = T. atroviride; H. virens = T. virens) an additional 14,000 ESTs from six other Trichoderma spp. (T. asperellum, H. lixii = T. harzianum, T. aggressivum var. europeae, T. longibrachiatum, T. cf. viride). The former three contained six, ten and nine members, respectively. Ten is the highest number found in any ascomycete so far. All the hydrophobins we examined had the conserved four beta-strands/one helix structure, which is stabilized by four disulfide bonds. In addition, a small number of these hydrophobins (HFBs)contained an extended N-terminus rich in either proline and aspartate, or glycine-asparagine. Phylogenetic analysis reveals a mosaic of terminal clades containing duplicated genes and shows only three reasonably supported clades. Calculation of the ratio of differences in synonymous vs. non-synonymous nucleotide substitutions provides evidence for strong purifying selection (KS/Ka >> 1). A genome database search for class II HFBs from other ascomycetes retrieved a much smaller number of hydrophobins (2-4) from each species, and most were from Sordariomycetes. A combined phylogeny of these sequences with those of Trichoderma showed that the Trichoderma HFBs mostly formed their own clades, whereas those of other Sordariomycetes occurred in shared clades.

CONCLUSION

Our study shows that the genus Trichoderma/Hypocrea has a proliferated arsenal of class II hydrophobins which arose by birth-and-death evolution followed by purifying selection.

Fungal CSL transcription factors

Background

The CSL (CBF1/RBP-Jκ/Suppressor of Hairless/LAG-1) transcription factor family members are well-known components of the transmembrane receptor Notch signaling pathway, which plays a critical role in metazoan development. They function as context-dependent activators or repressors of transcription of their responsive genes, the promoters of which harbor the GTG(G/A)GAA consensus elements. Recently, several studies described Notch-independent activities of the CSL proteins.

Conclusion

Our findings support the evolutionary origin of the CSL transcription factor family in the last common ancestor of fungi and metazoans. We hypothesize that the ancestral CSL function involved DNA binding and Notch-independent regulation of transcription and that this function may still be shared, to a certain degree, by the present CSL family members from both fungi and metazoans.

A mutation in the gene encoding cytochrome c1 leads to a decreased ROS content and to a long-lived phenotype in the filamentous fungus Podospora anserina

We present here the properties of a complex III loss-of-function mutant of the filamentous fungus Podospora anserina. The mutation corresponds to a single substitution in the second intron of the gene cyc1 encoding cytochrome c(1), leading to a splicing defect. The cyc1-1 mutant is long-lived, exhibits a defect in ascospore pigmentation, has a reduced growth rate and a reduced ROS production associated with a stabilisation of its mitochondrial DNA. We also show that increased longevity is linked with morphologically modified mitochondria and an increased number of mitochondrial genomes. Overexpression of the alternative oxidase rescues all these phenotypes and restores aging. Interestingly, the absence of complex III in this mutant is not paralleled with a deficiency in complex I activity as reported in mammals although the respiratory chain of P. anserina has recently been demonstrated to be organized according to the "respirasome" model.