Introns and splicing elements of five diverse fungi

Genes divided according to the relative position of the longest intron show increased representation in different KEGG pathways

Despite the fact that introns mean an energy and time burden for eukaryotic cells, they play an irreplaceable role in the diversification and regulation of protein production. As a common feature of eukaryotic genomes, it has been reported that in protein-coding genes, the longest intron is usually one of the first introns. The goal of our work was to find a possible difference in the biological function of genes that fulfill this common feature compared to genes that do not. Data on the lengths of all introns in genes were extracted from the genomes of six vertebrates (human, mouse, koala, chicken, zebrafish and fugu) and two other model organisms (nematode worm and arabidopsis). We showed that more than 40% of protein-coding genes have the relative position of the longest intron located in the second or third tertile of all introns. Genes divided according to the relative position of the longest intron were found to be significantly increased in different KEGG pathways. Genes with the longest intron in the first tertile predominate in a range of pathways for amino acid and lipid metabolism, various signaling, cell junctions or ABC transporters. Genes with the longest intron in the second or third tertile show increased representation in pathways associated with the formation and function of the spliceosome and ribosomes. In the two groups of genes defined in this way, we further demonstrated the difference in the length of the longest introns and the distribution of their absolute positions. We also pointed out other characteristics, namely the positive correlation between the length of the longest intron and the sum of the lengths of all other introns in the gene and the preservation of the exact same absolute and relative position of the longest intron between orthologous genes.

Unique features of conventional and nonconventional introns in Euglena gracilis

BACKGROUND

Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns.

RESULTS

A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans.

CONCLUSIONS

Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

A systematic screen identifies Saf5 as a link between splicing and transcription in fission yeast

Splicing is an important step of gene expression regulation in eukaryotes, as there are many mRNA precursors that can be alternatively spliced in different tissues, at different cell cycle phases or under different external stimuli. We have developed several integrated fluorescence-based in vivo splicing reporter constructs that allow the quantification of fission yeast splicing in vivo on intact cells, and we have compared their splicing efficiency in a wild type strain and in a prp2-1 (U2AF65) genetic background, showing a clear dependency between Prp2 and a consensus signal at 5' splicing site (5'SS). To isolate novel genes involved in regulated splicing, we have crossed the reporter showing more intron retention with the Schizosaccharomyces pombe knock out collection. Among the candidate genes involved in the regulation of splicing, we have detected strong splicing defects in two of the mutants -Δcwf12, a member of the NineTeen Complex (NTC) and Δsaf5, a methylosome subunit that acts together with the survival motor neuron (SMN) complex in small nuclear ribonucleoproteins (snRNP) biogenesis. We have identified that strains with mutations in cwf12 have inefficient splicing, mainly when the 5'SS differs from the consensus. However, although Δsaf5 cells also have some dependency on 5'SS sequence, we noticed that when one intron of a given pre-mRNA was affected, the rest of the introns of the same pre-mRNA had high probabilities of being also affected. This observation points Saf5 as a link between transcription rate and splicing.

An interpretable model of pre-mRNA splicing for animal and plant genes

Pre-mRNA splicing is a fundamental step in gene expression, conserved across eukaryotes, in which the spliceosome recognizes motifs at the 3' and 5' splice sites (SSs), excises introns, and ligates exons. SS recognition and pairing is often influenced by protein splicing factors (SFs) that bind to splicing regulatory elements (SREs). Here, we describe SMsplice, a fully interpretable model of pre-mRNA splicing that combines models of core SS motifs, SREs, and exonic and intronic length preferences. We learn models that predict SS locations with 83 to 86% accuracy in fish, insects, and plants and about 70% in mammals. Learned SRE motifs include both known SF binding motifs and unfamiliar motifs, and both motif classes are supported by genetic analyses. Our comparisons across species highlight similarities between non-mammals, increased reliance on intronic SREs in plant splicing, and a greater reliance on SREs in mammalian splicing.

Unveiling a Microexon Switch: Novel Regulation of the Activities of Sugar Assimilation and Plant-Cell-Wall-Degrading Xylanases and Cellulases by Xlr2 in Trichoderma virens

Functional microexons have not previously been described in filamentous fungi. Here, we describe a novel mechanism of transcriptional regulation in Trichoderma requiring the inclusion of a microexon from the Xlr2 gene. In low-glucose environments, a long mRNA including the microexon encodes a protein with a GAL4-like DNA-binding domain (Xlr2-α), whereas in high-glucose environments, a short mRNA that is produced encodes a protein lacking this DNA-binding domain (Xlr2-β). Interestingly, the protein isoforms differ in their impact on cellulase and xylanase activity. Deleting the Xlr2 gene reduced both xylanase and cellulase activity and growth on different carbon sources, such as carboxymethylcellulose, xylan, glucose, and arabinose. The overexpression of either Xlr2-α or Xlr2-β in T. virens showed that the short isoform (Xlr2-β) caused higher xylanase activity than the wild types or the long isoform (Xlr2-α). Conversely, cellulase activity did not increase when overexpressing Xlr2-β but was increased with the overexpression of Xlr2-α. This is the first report of a novel transcriptional regulation mechanism of plant-cell-wall-degrading enzyme activity in T. virens. This involves the differential expression of a microexon from a gene encoding a transcriptional regulator.

Insights in inulin binding and inulin oligosaccharide formation by novel multi domain endo-inulinases from Botrytis cinerea

World-wide, pathogenic fungi such as Botrytis cinerea cause tremendous yield losses in terms of food production and post-harvest food decay. Many fungi produce inulin-type oligosaccharides (IOSs) from inulin through endo-inulinases which typically show a two domain structure. B.cinerea lacks a two domain endo-inulinase but contains a three domain structure instead. Genome mining revealed three and four domain (d4) enzymes in the fungal kingdom. Here, three and two domain enzymes were compared in their capacity to produce IOSs from inulin. Hill kinetics were observed in three domain enzymes as compared to Michaelis-Menten kinetics in two domain enzymes, suggesting that the N-terminal extension functions as a carbohydrate binding module. Analysis of the IOS product profiles generated from purified GF6, GF12, GF16 and GF18 inulins and extensive sugar docking approaches led to enhanced insights in the active site functioning, revealing subtle differences between the endo-inulinases from Aspergillus niger and B. cinerea. Improved insights in structure-function relationships in fungal endo-inulinases offer opportunities to develop superior enzymes for the production of specific IOS formulations to improve plant and animal health (priming agents, prebiotics).

Near chromosome-level genome assembly of the microsporidium Hamiltosporidium tvaerminnensis

Microsporidia are intracellular parasitic fungi whose genomes rank among the smallest of all known eukaryotes. A number of outstanding questions remain concerning the evolution of their large-scale variation in genome architecture, responsible for genome size variation of more than an order of magnitude. This genome report presents the first near-chromosomal assembly of a large-genome microsporidium, Hamiltosporidium tvaerminnensis. Combined Oxford Nanopore, Pacific Biosciences (PacBio), and Illumina sequencing led to a genome assembly of 17 contigs, 11 of which represent complete chromosomes. Our assembly is 21.64 Mb in length, has an N50 of 1.44 Mb, and consists of 39.56% interspersed repeats. We introduce a novel approach in microsporidia, PacBio Iso-Seq, as part of a larger annotation pipeline for obtaining high-quality annotations of 3,573 protein-coding genes. Based on direct evidence from the full-length Iso-Seq transcripts, we present evidence for alternative polyadenylation and variation in splicing efficiency, which are potential regulation mechanisms for gene expression in microsporidia. The generated high-quality genome assembly is a necessary resource for comparative genomics that will help elucidate the evolution of genome architecture in response to intracellular parasitism.

The Biological Action and Structural Characterization of Eryngitin 3 and 4, Ribotoxin-like Proteins from Pleurotus eryngii Fruiting Bodies

Ribotoxin-like proteins (RL-Ps) are specific ribonucleases found in mushrooms that are able to cleave a single phosphodiester bond located in the sarcin-ricin loop (SRL) of the large rRNA. The cleaved SRL interacts differently with some ribosomal proteins (P-stalk). This action blocks protein synthesis because the damaged ribosomes are unable to interact with elongation factors. Here, the amino acid sequences of eryngitin 3 and 4, RL-Ps isolated from Pleurotus eryngii fruiting bodies, were determined to (i) obtain structural information on this specific ribonuclease family from edible mushrooms and (ii) explore the structural determinants which justify their different biological and antipathogenic activities. Indeed, eryngitin 3 exhibited higher toxicity with respect to eryngitin 4 against tumoral cell lines and model fungi. Structurally, eryngitin 3 and 4 consist of 132 amino acids, most of them identical and exhibiting a single free cysteinyl residue. The amino acidic differences between the two toxins are (i) an additional phenylalanyl residue at the N-terminus of eryngitin 3, not retrieved in eryngitin 4, and (ii) an additional arginyl residue at the C-terminus of eryngitin 4, not retrieved in eryngitin 3. The 3D models of eryngitins show slight differences at the N- and C-terminal regions. In particular, the positive electrostatic surface at the C-terminal of eryngitin 4 is due to the additional arginyl residue not retrieved in eryngitin 3. This additional positive charge could interfere with the binding to the SRL (substrate) or with some ribosomal proteins (P-stalk structure) during substrate recognition.

Genome-wide identification of the opsin protein in Leptosphaeria maculans and comparison with other fungi (pathogens of Brassica napus)

The largest family of transmembrane receptors are G-protein-coupled receptors (GPCRs). These receptors respond to perceived environmental signals and infect their host plants. Family A of the GPCR includes opsin. However, there is little known about the roles of GPCRs in phytopathogenic fungi. We studied opsin in Leptosphaeria maculans, an important pathogen of oilseed rape (Brassica napus) that causes blackleg disease, and compared it with six other fungal pathogens of oilseed rape. A phylogenetic tree analysis of 31 isoforms of the opsin protein showed six major groups and six subgroups. All three opsin isoforms of L. maculans are grouped in the same clade in the phylogenetic tree. Physicochemical analysis revealed that all studied opsin proteins are stable and hydrophobic. Subcellular localization revealed that most isoforms were localized in the endoplasmic reticulum membrane except for several isoforms in Verticillium species, which were localized in the mitochondrial membrane. Most isoforms comprise two conserved domains. One conserved motif was observed across all isoforms, consisting of the BACTERIAL_OPSIN_1 domain, which has been hypothesized to have an identical sensory function. Most studied isoforms showed seven transmembrane helices, except for one isoform of V. longisporum and four isoforms of Fusarium oxysporum. Tertiary structure prediction displayed a conformational change in four isoforms of F. oxysporum that presumed differences in binding to other proteins and sensing signals, thereby resulting in various pathogenicity strategies. Protein-protein interactions and binding site analyses demonstrated a variety of numbers of ligands and pockets across all isoforms, ranging between 0 and 13 ligands and 4 and 10 pockets. According to the phylogenetic analysis in this study and considerable physiochemically and structurally differences of opsin proteins among all studied fungi hypothesized that this protein acts in the pathogenicity, growth, sporulation, and mating of these fungi differently.

Quantitative and qualitative plant-pathogen interactions call upon similar pathogenicity genes with a spectrum of effects

Septoria leaf blotch is a foliar wheat disease controlled by a combination of plant genetic resistances and fungicides use. R-gene-based qualitative resistance durability is limited due to gene-for-gene interactions with fungal avirulence (Avr) genes. Quantitative resistance is considered more durable but the mechanisms involved are not well documented. We hypothesize that genes involved in quantitative and qualitative plant-pathogen interactions are similar. A bi-parental population of Zymoseptoria tritici was inoculated on wheat cultivar 'Renan' and a linkage analysis performed to map QTL. Three pathogenicity QTL, Qzt-I05-1, Qzt-I05-6 and Qzt-I07-13, were mapped on chromosomes 1, 6 and 13 in Z. tritici, and a candidate pathogenicity gene on chromosome 6 was selected based on its effector-like characteristics. The candidate gene was cloned by Agrobacterium tumefaciens-mediated transformation, and a pathology test assessed the effect of the mutant strains on 'Renan'. This gene was demonstrated to be involved in quantitative pathogenicity. By cloning a newly annotated quantitative-effect gene in Z. tritici that is effector-like, we demonstrated that genes underlying pathogenicity QTL can be similar to Avr genes. This opens up the previously probed possibility that 'gene-for-gene' underlies not only qualitative but also quantitative plant-pathogen interactions in this pathosystem.

A genetic tool to express long fungal biosynthetic genes

BACKGROUND

Secondary metabolites (SMs) from mushroom-forming fungi (Basidiomycota) and early diverging fungi (EDF) such as Mucoromycota are scarcely investigated. In many cases, production of SMs is induced by unknown stress factors or is accompanied by seasonable developmental changes on fungal morphology. Moreover, many of these fungi are considered as non-culturable under laboratory conditions which impedes investigation into SM. In the post-genomic era, numerous novel SM genes have been identified especially from EDF. As most of them encode multi-module enzymes, these genes are usually long which limits cloning and heterologous expression in traditional hosts.

RESULTS

An expression system in Aspergillus niger is presented that is suitable for the production of SMs from both Basidiomycota and EDF. The akuB gene was deleted in the expression host A. niger ATNT∆pyrG, resulting in a deficient nonhomologous end-joining repair mechanism which in turn facilitates the targeted gene deletion via homologous recombination. The ∆akuB mutant tLK01 served as a platform to integrate overlapping DNA fragments of long SM genes into the fwnA locus required for the black pigmentation of conidia. This enables an easy discrimination of correct transformants by screening the transformation plates for fawn-colored colonies. Expression of the gene of interest (GOI) is induced dose-dependently by addition of doxycycline and is enhanced by the dual TetON/terrein synthase promoter system (ATNT) from Aspergillus terreus. We show that the 8 kb polyketide synthase gene lpaA from the basidiomycete Laetiporus sulphureus is correctly assembled from five overlapping DNA fragments and laetiporic acids are produced. In a second approach, we expressed the yet uncharacterized > 20 kb nonribosomal peptide synthetase gene calA from the EDF Mortierella alpina. Gene expression and subsequent LC-MS/MS analysis of mycelial extracts revealed the production of the antimycobacterial compound calpinactam. This is the first report on the heterologous production of a full-length SM multidomain enzyme from EDF.

CONCLUSIONS

The system allows the assembly, targeted integration and expression of genes of > 20 kb size in A. niger in one single step. The system is suitable for evolutionary distantly related SM genes from both Basidiomycota and EDF. This uncovers new SM resources including genetically intractable or non-culturable fungi.

A broadly conserved fungal alcohol oxidase (AOX) facilitates fungal invasion of plants

Alcohol oxidases (AOXs) are ecologically important enzymes that facilitate a number of plant-fungal interactions. Within Ascomycota they are primarily associated with methylotrophy, as a peroxisomal AOX catalysing the conversion of methanol to formaldehyde in methylotrophic yeast. In this study we demonstrate that AOX orthologues are phylogenetically conserved proteins that are common in the genomes of nonmethylotrophic, plant-associating fungi. Additionally, AOX orthologues are highly expressed during infection in a range of diverse pathosystems. To study the role of AOX in plant colonization, AOX knockout mutants were generated in the broad host range pathogen Sclerotinia sclerotiorum. Disease assays in soybean showed that these mutants had a significant virulence defect as evidenced by markedly reduced stem lesions and mortality rates. Chemical genomics suggested that SsAOX may function as an aromatic AOX, and growth assays demonstrated that ΔSsAOX is incapable of properly utilizing plant extract as a nutrient source. Profiling of known aromatic alcohols pointed towards the monolignol coniferyl alcohol (CA) as a possible substrate for SsAOX. As CA and other monolignols are ubiquitous among land plants, the presence of highly conserved AOX orthologues throughout Ascomycota implies that this is a broadly conserved protein used by ascomycete fungi during plant colonization.

Biogenesis of telomerase RNA from a protein-coding mRNA precursor

Telomerase is a eukaryotic ribonucleoprotein (RNP) enzyme that adds DNA repeats onto chromosome ends to maintain genomic stability and confer cellular immortality in cancer and stem cells. The telomerase RNA (TER) component is essential for telomerase catalytic activity and provides the template for telomeric DNA synthesis. The biogenesis of TERs is extremely divergent across eukaryotic kingdoms, employing distinct types of transcription machinery and processing pathways. In ciliates and plants, TERs are transcribed by RNA polymerase III (Pol III), while animal and ascomycete fungal TERs are transcribed by RNA Pol II and share biogenesis pathways with small nucleolar RNA (snoRNA) and small nuclear RNA (snRNA), respectively. Here, we report an unprecedented messenger RNA (mRNA)-derived biogenesis pathway for the 1,291 nucleotide TER from the basidiomycete fungus Ustilago maydis. The U. maydis TER (UmTER) contains a 5'-monophosphate, distinct from the 5' 2,2,7-trimethylguanosine (TMG) cap common to animal and ascomycete fungal TERs. The mature UmTER is processed from the 3'-untranslated region (3'-UTR) of a larger RNA precursor that possesses characteristics of mRNA including a 5' 7-methyl-guanosine (m7G) cap, alternative splicing of introns, and a poly(A) tail. Moreover, this mRNA transcript encodes a protein called Early meiotic induction protein 1 (Emi1) that is conserved across dikaryotic fungi. A recombinant UmTER precursor expressed from an mRNA promoter is processed correctly to yield mature UmTER, confirming an mRNA-processing pathway for producing TER. Our findings expand the plethora of TER biogenesis mechanisms and demonstrate a pathway for producing a functional long noncoding RNA from a protein-coding mRNA precursor.

Splicing of branchpoint-distant exons is promoted by Cactin, Tls1 and the ubiquitin-fold-activated Sde2

Intron diversity facilitates regulated gene expression and alternative splicing. Spliceosomes excise introns after recognizing their splicing signals: the 5'-splice site (5'ss), branchpoint (BP) and 3'-splice site (3'ss). The latter two signals are recognized by U2 small nuclear ribonucleoprotein (snRNP) and its accessory factors (U2AFs), but longer spacings between them result in weaker splicing. Here, we show that excision of introns with a BP-distant 3'ss (e.g. rap1 intron 2) requires the ubiquitin-fold-activated splicing regulator Sde2 in Schizosaccharomyces pombe. By monitoring splicing-specific ura4 reporters in a collection of S. pombe mutants, Cay1 and Tls1 were identified as additional regulators of this process. The role of Sde2, Cay1 and Tls1 was further confirmed by increasing BP-3'ss spacings in a canonical tho5 intron. We also examined BP-distant exons spliced independently of these factors and observed that RNA secondary structures possibly bridged the gap between the two signals. These proteins may guide the 3'ss towards the spliceosome's catalytic centre by folding the RNA between the BP and 3'ss. Orthologues of Sde2, Cay1 and Tls1, although missing in the intron-poor Saccharomyces cerevisiae, are present in intron-rich eukaryotes, including humans. This type of intron-specific pre-mRNA splicing appears to have evolved for regulated gene expression and alternative splicing of key heterochromatin factors.

Transcript-specific determinants of pre-mRNA splicing revealed through in vivo kinetic analyses of the 1st and 2nd chemical steps

We generate high-precision measurements of the in vivo rates of both chemical steps of pre-mRNA splicing across the genome-wide complement of substrates in yeast by coupling metabolic labeling, multiplexed primer-extension sequencing, and kinetic modeling. We demonstrate that the rates of intron removal vary widely, splice-site sequences are primary determinants of 1st step but have little apparent impact on 2nd step rates, and the 2nd step is generally faster than the 1st step. Ribosomal protein genes (RPGs) are spliced faster than non-RPGs at each step, and RPGs share evolutionarily conserved properties that may contribute to their faster splicing. A genetic variant defective in the 1st step of the pathway reveals a genome-wide defect in the 1st step but an unexpected, transcript-specific change in the 2nd step. Our work demonstrates that extended co-transcriptional association is an important determinant of splicing rate, a conclusion at odds with recent claims of ultra-fast splicing.

Unique and Repeated Stwintrons (Spliceosomal Twin Introns) in the Hypoxylaceae

Introns are usually non-coding sequences interrupting open reading frames in pre-mRNAs [D1,2]. Stwintrons are nested spliceosomal introns, where an internal intron splits a second donor sequence into two consecutive splicing reactions leading to mature mRNA. In Hypoxylon sp. CO27-5, 36 highly sequence-similar [D1,2] stwintrons are extant (sister stwintrons). An additional 81 [D1,2] sequence-unrelated stwintrons are described here. Most of them are located at conserved gene positions rooted deep in the Hypoxylaceae. Absence of exonic sequence bias at the exon–stwintron junctions and a very similar phase distribution were noted for both groups. The presence of an underlying sequence symmetry in all 117 stwintrons was striking. This symmetry, more pronounced near the termini of most of the full-length sister stwintrons, may lead to a secondary structure that brings into close proximity the most distal splice sites, the donor of the internal and the acceptor of the external intron. The Hypoxylon stwintrons were overwhelmingly excised by consecutive splicing reactions precisely removing the whole intervening sequence, whereas one excision involving the distal splice sites led to a frameshift. Alternative (mis)splicing took place for both sister and uniquely occurring stwintrons. The extraordinary symmetry of the sister stwintrons thus seems dispensable for the infrequent, direct utilisation of the distal splice sites.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

Characterization of cross-species transcription and splicing from Penicillium to Saccharomyces cerevisiae

Heterologous expression of eukaryotic gene clusters in yeast has been widely used for producing high-value chemicals and bioactive secondary metabolites. However, eukaryotic transcription cis-elements are still undercharacterized, and the cross-species expression mechanism remains poorly understood. Here we used the whole expression unit (including original promoter, terminator, and open reading frame with introns) of orotidine 5'-monophosphate decarboxylases from 14 Penicillium species as a showcase, and analyzed their cross-species expression in Saccharomyces cerevisiae. We found that pyrG promoters from the Penicillium species could drive URA3 expression in yeast, and that inefficient cross-species splicing of Penicillium introns might result in weak cross-species expression. Thus, this study demonstrates cross-species expression from Penicillium to yeast, and sheds light on the opportunities and challenges of cross-species expression of fungi expression units and gene clusters in yeast without refactoring for novel natural product discovery.

The long non-coding RNA landscape of Candida yeast pathogens

Long non-coding RNAs (lncRNAs) constitute a poorly studied class of transcripts with emerging roles in key cellular processes. Despite efforts to characterize lncRNAs across a wide range of species, these molecules remain largely unexplored in most eukaryotic microbes, including yeast pathogens of the Candida clade. Here, we analyze thousands of publicly available sequencing datasets to infer and characterize the lncRNA repertoires of five major Candida pathogens: Candida albicans, Candida tropicalis, Candida parapsilosis, Candida auris and Candida glabrata. Our results indicate that genomes of these species encode hundreds of lncRNAs that show levels of evolutionary constraint intermediate between those of intergenic genomic regions and protein-coding genes. Despite their low sequence conservation across the studied species, some lncRNAs are syntenic and are enriched in shared sequence motifs. We find co-expression of lncRNAs with certain protein-coding transcripts, hinting at potential functional associations. Finally, we identify lncRNAs that are differentially expressed during infection of human epithelial cells for four of the studied species. Our comprehensive bioinformatic analyses of Candida lncRNAs pave the way for future functional characterization of these transcripts.

Long non-coding RNAs (lncRNAs) play roles in key cellular processes, but remain largely unexplored in fungal pathogens such as Candida. Here, Hovhannisyan and Gabaldón analyze thousands of sequencing datasets to infer and characterize the lncRNA repertoires of five Candida species, paving the way for their future functional characterization.

Evolution of the Early Spliceosomal Complex-From Constitutive to Regulated Splicing

Pre-mRNA splicing is a major process in the regulated expression of genes in eukaryotes, and alternative splicing is used to generate different proteins from the same coding gene. Splicing is a catalytic process that removes introns and ligates exons to create the RNA sequence that codifies the final protein. While this is achieved in an autocatalytic process in ancestral group II introns in prokaryotes, the spliceosome has evolved during eukaryogenesis to assist in this process and to finally provide the opportunity for intron-specific splicing. In the early stage of splicing, the RNA 5' and 3' splice sites must be brought within proximity to correctly assemble the active spliceosome and perform the excision and ligation reactions. The assembly of this first complex, termed E-complex, is currently the least understood process. We focused in this review on the formation of the E-complex and compared its composition and function in three different organisms. We highlight the common ancestral mechanisms in S. cerevisiae, S. pombe, and mammals and conclude with a unifying model for intron definition in constitutive and regulated co-transcriptional splicing.

Cross talk between the upstream exon-intron junction and Prp2 facilitates splicing of non-consensus introns

Splicing of mRNA precursors is essential in the regulation of gene expression. U2AF65 recognizes the poly-pyrimidine tract and helps in the recognition of the branch point. Inactivation of fission yeast U2AF65 (Prp2) blocks splicing of most, but not all, pre-mRNAs, for reasons that are not understood. Here, we have determined genome-wide the splicing efficiency of fission yeast cells as they progress into synchronous meiosis in the presence or absence of functional Prp2. Our data indicate that in addition to the splicing elements at the 3' end of any intron, the nucleotides immediately upstream the intron will determine whether Prp2 is required or dispensable for splicing. By changing those nucleotides in any given intron, we regulate its Prp2 dependency. Our results suggest a model in which Prp2 is required for the coordinated recognition of both intronic ends, placing Prp2 as a key regulatory element in the determination of the exon-intron boundaries.

Intron distribution and emerging role of alternative splicing in fungi

Spliceosomal introns are noncoding sequences that are spliced from pre-mRNA. They are ubiquitous in eukaryotic genomes, although the average number of introns per gene varies considerably between different eukaryotic species. Fungi are diverse in terms of intron numbers ranging from 4% to 99% genes with introns. Alternative splicing is one of the most common modes of posttranscriptional regulation in eukaryotes, giving rise to multiple transcripts from a single pre-mRNA and is widespread in metazoans and drives extensive proteome diversity. Earlier, alternative splicing was considered to be rare in fungi, but recently, increasing numbers of studies have revealed that alternative splicing is also widespread in fungi and has been implicated in the regulation of fungal growth and development, protein localization and the improvement of survivability, likely underlying their unique capacity to adapt to changing environmental conditions. However, the role of alternative splicing in pathogenicity and development of drug resistance is only recently gaining attention. In this review, we describe the intronic landscape in fungi. We also present in detail the newly discovered functions of alternative splicing in various cellular processes and outline areas particularly in pathogenesis and clinical drug resistance for future studies that could lead to the development of much needed new therapeutics.

Modeling in yeast how rDNA introns slow growth and increase desiccation tolerance in lichens

We connect ribosome biogenesis to desiccation tolerance in lichens, widespread symbioses between specialized fungi (mycobionts) and unicellular phototrophs. We test whether the introns present in the nuclear ribosomal DNA of lichen mycobionts contribute to their anhydrobiosis. Self-splicing introns are found in the rDNA of several eukaryotic microorganisms, but most introns populating lichen rDNA are unable to self-splice, being either catalytically impaired group I introns, or spliceosomal introns ectopically present in rDNA. Although the mycobiont’s splicing machinery removes all introns from rRNA, Northern analysis indicates delayed post-transcriptional removal during rRNA processing, suggesting interference with ribosome assembly. To study the effects of lichen introns in a model system, we used CRISPR to introduce a spliceosomal rDNA intron from the lichen fungus Cladonia grayi into all nuclear rDNA copies of Saccharomyces cerevisiae, which lacks rDNA introns. Three intron-bearing yeast mutants were constructed with the intron inserted either in the 18S rRNA genes, the 25S rRNA genes, or in both. The mutants removed the introns correctly but had half the rDNA genes of the wildtype, grew 4.4–6 times slower, and were 40–1700 times more desiccation tolerant depending on intron position and number. Intracellular trehalose, a disaccharide implicated in desiccation tolerance, was detected at low concentration. Our data suggest that the interference of the splicing machinery with ribosome assembly leads to fewer ribosomes and proteins and to slow growth and increased desiccation tolerance in the yeast mutants. The relevance of these findings for slow growth and desiccation tolerance in lichens is discussed.

Internally Symmetrical Stwintrons and Related Canonical Introns in Hypoxylaceae Species

Spliceosomal introns are pervasive in eukaryotes. Intron gains and losses have occurred throughout evolution, but the origin of new introns is unclear. Stwintrons are complex intervening sequences where one of the sequence elements (5′-donor, lariat branch point element or 3′-acceptor) necessary for excision of a U2 intron (external intron) is itself interrupted by a second (internal) U2 intron. In Hypoxylaceae, a family of endophytic fungi, we uncovered scores of donor-disrupted stwintrons with striking sequence similarity among themselves and also with canonical introns. Intron–exon structure comparisons suggest that these stwintrons have proliferated within diverging taxa but also give rise to proliferating canonical introns in some genomes. The proliferated (stw)introns have integrated seamlessly at novel gene positions. The recently proliferated (stw)introns appear to originate from a conserved ancestral stwintron characterised by terminal inverted repeats (45–55 nucleotides), a highly symmetrical structure that may allow the formation of a double-stranded intron RNA molecule. No short tandem duplications flank the putatively inserted intervening sequences, which excludes a DNA transposition-based mechanism of proliferation. It is tempting to suggest that this highly symmetrical structure may have a role in intron proliferation by (an)other mechanism(s).

Evidence for the Role of CYP51A and Xenobiotic Detoxification in Differential Sensitivity to Azole Fungicides in Boxwood Blight Pathogens

Boxwood blight, a fungal disease of ornamental plants (Buxus spp.), is caused by two sister species, Calonectria pseudonaviculata (Cps) and C. henricotiae (Che). Compared to Cps, Che is documented to display reduced sensitivity to fungicides, including the azole class of antifungals, which block synthesis of a key fungal membrane component, ergosterol. A previous study reported an ergosterol biosynthesis gene in Cps, CYP51A, to be a pseudogene, and RNA-Seq data confirm that a functional CYP51A is expressed only in Che. The lack of additional ergosterol biosynthesis genes showing significant differential expression suggests that the functional CYP51A in Che could contribute to reduced azole sensitivity when compared to Cps. RNA-Seq and bioinformatic analyses found that following azole treatment, 55 genes in Cps, belonging to diverse pathways, displayed a significant decrease in expression. Putative xenobiotic detoxification genes overexpressed in tetraconazole-treated Che encoded predicted monooxygenase and oxidoreductase enzymes. In summary, expression of a functional CYP51A gene and overexpression of predicted xenobiotic detoxification genes appear likely to contribute to differential fungicide sensitivity in these two sister taxa.

Identification and Expression Analysis of Sugar Transporter Gene Family in Aspergillus oryzae

Sugar transporter (SUT) genes are associated with multiple physiological and biochemical processes in filamentous fungi, such as the response to various stresses. However, limited systematic analysis and functional information of SUT gene family have been available on Aspergillus oryzae (A. oryzae). To investigate the potential roles of SUTs in A. oryzae, we performed an integrative analysis of the SUT gene family in this study. Based on the conserved protein domain search, 127 putative SUT genes were identified in A. oryzae and further categorized into eight distinct subfamilies. The result of gene structure and conserved motif analysis illustrated functional similarities among the AoSUT proteins within the same subfamily. Additionally, expression profiles of the AoSUT genes at different growth stages elucidated that most of AoSUT genes have high expression levels at the stationary phase while low in the adaptive phase. Furthermore, expression profiles of AoSUT genes under salt stress showed that AoSUT genes may be closely linked to salt tolerance and involved in sophisticated transcriptional process. The protein-protein interaction network of AoSUT propounded some potentially interacting proteins. A comprehensive overview of the AoSUT gene family will offer new insights into the structural and functional features as well as facilitate further research on the roles of AoSUT genes in response to abiotic stresses.

Gene Organization, Expression, and Localization of Ribotoxin-Like Protein Ageritin in Fruiting Body and Mycelium of Agrocybe aegerita

The edible mushroom Agrocybe aegerita produces a ribotoxin-like protein known as Ageritin. In this work, the gene encoding Ageritin was characterized by sequence analysis. It contains several typical features of fungal genes such as three short introns (60, 55 and 69 bp) located at the 5' region of the coding sequence and typical splice junctions. This sequence codes for a precursor of 156 amino acids (~17-kDa) containing an additional N-terminal peptide of 21 amino acid residues, absent in the purified toxin (135 amino acid residues; ~15-kDa). The presence of 17-kDa and 15-kDa forms was investigated by Western blot in specific parts of fruiting body and in mycelia of A. aegerita. Data show that the 15-kDa Ageritin is the only form retrieved in the fruiting body and the principal form in mycelium. The immunolocalization by confocal laser scanning microscopy and transmission electron microscopy proves that Ageritin has vacuolar localization in hyphae. Coupling these data with a bioinformatics approach, we suggest that the N-terminal peptide of Ageritin (not found in the purified toxin) is a new signal peptide in fungi involved in intracellular routing from endoplasmic reticulum to vacuole, necessary for self-defense of A. aegerita ribosomes from Ageritin toxicity.

Modification of transcriptional factor ACE3 enhances protein production in Trichoderma reesei in the absence of cellulase gene inducer

BACKGROUND

Trichoderma reesei is one of the best-known cellulolytic organisms, producing large quantities of a complete set of extracellular cellulases and hemicellulases for the degradation of lignocellulosic substances. Hence, T. reesei is a biotechnically important host and it is used commercially in enzyme production, of both native and foreign origin. Many strategies for producing enzymes in T. reesei rely on the cbh1 and other cellulase gene promoters for high-level expression and these promoters require induction by sophorose, lactose or other inducers for high productivity during manufacturing.

RESULTS

We described an approach for producing high levels of secreted proteins by overexpression of a transcription factor ACE3 in T. reesei. We refined the ace3 gene structure and identified specific ACE3 variants that enable production of secreted cellulases and hemicellulases on glucose as a sole carbon source (i.e., in the absence of an inducer). These specific ACE3 variants contain a full-length Zn₂Cys₆ binuclear cluster domain at the N-terminus and a defined length of truncations at the C-terminus. When expressed at a moderate level in the fungal cells, the ACE3 variants can induce high-level expression of cellulases and hemicellulases on glucose (i.e., in the absence of an inducer), and further improve expression on lactose or glucose/sophorose (i.e., in the presence of an inducer). Finally, we demonstrated that this method is applicable to industrial strains and fermentation conditions, improving protein production both in the absence and in the presence of an inducer.

CONCLUSIONS

This study demonstrates that overexpression of ACE3 variants enables a high level of protein production in the absence of an inducer, and boosts protein production in the presence of an inducer. It is an efficient approach to increase protein productivity and to reduce manufacturing costs.

Genome-Wide Characterization and Comparative Analysis of MYB Transcription Factors in Ganoderma Species

Numerous studies in plants have shown the vital roles of MYB transcription factors in signal transduction, developmental regulation, biotic/abiotic stress responses and secondary metabolism regulation. However, less is known about the functions of MYBs in Ganoderma. In this study, five medicinal macrofungi of genus Ganoderma were subjected to a genome-wide comparative analysis of MYB genes. A total of 75 MYB genes were identified and classified into four types: 1R-MYBs (52), 2R-MYBs (19), 3R-MYBs (2) and 4R-MYBs (2). Gene structure analysis revealed varying exon numbers (3-14) and intron lengths (7-1058 bp), and noncanonical GC-AG introns were detected in G. lucidum and G. sinense. In a phylogenetic analysis, 69 out of 75 MYB genes were clustered into 15 subgroups, and both single-copy orthologous genes and duplicated genes were identified. The promoters of the MYB genes harbored multiple cis-elements, and specific genes were co-expressed with the G. lucidum MYB genes, indicating the potential roles of these MYB genes in stress response, development and metabolism. This comprehensive and systematic study of MYB family members provides a reference and solid foundation for further functional analysis of MYB genes in Ganoderma species.

The functional study of human proteins using humanized yeast

The functional and optimal expression of genes is crucial for survival of all living organisms. Numerous experiments and efforts have been performed to reveal the mechanisms required for the functional and optimal expression of human genes. The yeast Saccharomyces cerevisiae has evolved independently of humans for billions of years. Nevertheless, S. cerevisiae has many conserved genes and expression mechanisms that are similar to those in humans. Yeast is the most commonly used model organism for studying the function and expression mechanisms of human genes because it has a relatively simple genome structure, which is easy to manipulate. Many previous studies have focused on understanding the functions and mechanisms of human proteins using orthologous genes and biological systems of yeast. In this review, we mainly introduce two recent studies that replaced human genes and nucleosomes with those of yeast. Here, we suggest that, although yeast is a relatively small eukaryotic cell, its humanization is useful for the direct study of human proteins. In addition, yeast can be used as a model organism in a broader range of studies, including drug screening.

Complex intron generation in the yeast genus Lipomyces

In primary transcripts of eukaryotic nuclear genes, coding sequences are often interrupted by U2-type introns. Such intervening sequences can constitute complex introns excised by consecutive splicing reactions. The origin of spliceosomal introns is a vexing problem. Sequence variation existent across fungal taxa provides means to study their structure and evolution. In one class of complex introns called [D] stwintrons, an (internal) U2 intron is nested within the 5'-donor element of another (external) U2 intron. In the gene for a reticulon-like protein in species of the ascomycete yeast genus Lipomyces, the most 5' terminal intron position is occupied by one of three complex intervening sequences consistent of differently nested U2 intron units, as demonstrated in L. lipofer, L. suomiensis, and L. starkeyi. In L. starkeyi, the donor elements of the constituent introns are abutting and the complex intervening sequence can be excised alternatively either with one standard splicing reaction or, as a [D] stwintron, by two consecutive reactions. Our work suggests how [D] stwintrons could emerge by the appearance of new functional splice sites within an extant intron. The stepwise stwintronisation mechanism may involve duplication of the functional intron donor element of the ancestor intron.

A Silent Exonic Mutation in a Rice Integrin-α FG-GAP Repeat-Containing Gene Causes Male-Sterility by Affecting mRNA Splicing

Pollen development plays crucial roles in the life cycle of higher plants. Here we characterized a rice mutant with complete male-sterile phenotype, pollen-less 1 (pl1). pl1 exhibited smaller anthers with arrested pollen development, absent Ubisch bodies, necrosis-like tapetal hypertrophy, and smooth anther cuticular surface. Molecular mapping revealed a synonymous mutation in the fourth exon of PL1 co-segregated with the mutant phenotype. This mutation disrupts the exon-intron splice junction in PL1, generating aberrant mRNA species and truncated proteins. PL1 is highly expressed in the tapetal cells of developing anther, and its protein is co-localized with plasma membrane (PM) and endoplasmic reticulum (ER) signal. PL1 encodes an integrin-α FG-GAP repeat-containing protein, which has seven β-sheets and putative Ca²⁺-binding motifs and is broadly conserved in terrestrial plants. Our findings therefore provide insights into both the role of integrin-α FG-GAP repeat-containing protein in rice male fertility and the influence of exonic mutation on intronic splice donor site selection.

Reconstitution of biosynthetic machinery of fungal natural products in heterologous hosts

Ascomycota and basidiomycota fungi are prolific sources of biologically active natural products. Recent genomic data and bioinformatic analysis indicate that fungi possess a large number of biosynthetic gene clusters for bioactive natural products but more than 90% are silent. Heterologous expression in the filamentous fungi as hosts is one of the powerful tools to expression of the silent gene clusters. This review introduces recent studies on the total biosynthesis of representative family members via common platform intermediates, genome mining of novel di- and sesterterpenoids including detailed cyclization pathway, and development of expression host for basidiomycota genes with efficient genome editing method. In addition, this review will discuss the several strategies, for the generation of structural diversity, which are found through these studies.

Animal, Fungi, and Plant Genome Sequences Harbor Different Non-Canonical Splice Sites

Most protein-encoding genes in eukaryotes contain introns, which are interwoven with exons. Introns need to be removed from initial transcripts in order to generate the final messenger RNA (mRNA), which can be translated into an amino acid sequence. Precise excision of introns by the spliceosome requires conserved dinucleotides, which mark the splice sites. However, there are variations of the highly conserved combination of GT at the 5' end and AG at the 3' end of an intron in the genome. GC-AG and AT-AC are two major non-canonical splice site combinations, which have been known for years. Recently, various minor non-canonical splice site combinations were detected with numerous dinucleotide permutations. Here, we expand systematic investigations of non-canonical splice site combinations in plants across eukaryotes by analyzing fungal and animal genome sequences. Comparisons of splice site combinations between these three kingdoms revealed several differences, such as an apparently increased CT-AC frequency in fungal genome sequences. Canonical GT-AG splice site combinations in antisense transcripts are a likely explanation for this observation, thus indicating annotation errors. In addition, high numbers of GA-AG splice site combinations were observed in Eurytemoraaffinis and Oikopleuradioica. A variant in one U1 small nuclear RNA (snRNA) isoform might allow the recognition of GA as a 5' splice site. In depth investigation of splice site usage based on RNA-Seq read mappings indicates a generally higher flexibility of the 3' splice site compared to the 5' splice site across animals, fungi, and plants.

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.

A dispensable paralog of succinate dehydrogenase subunit C mediates standing resistance towards a subclass of SDHI fungicides in Zymoseptoria tritici

Succinate dehydrogenase inhibitor (SDHI) fungicides are widely used for the control of a broad range of fungal diseases. This has been the most rapidly expanding fungicide group in terms of new molecules discovered and introduced for agricultural use over the past fifteen years. A particular pattern of differential sensitivity (resistance) to the stretched heterocycle amide SDHIs (SHA-SDHIs), a subclass of chemically-related SDHIs, was observed in naïve Zymoseptoria tritici populations not previously exposed to these chemicals. Subclass-specific resistance was confirmed at the enzyme level but did not correlate with the genotypes of the succinate dehydrogenase (SDH) encoding genes. Mapping and characterization of the molecular mechanisms responsible for standing SHA-SDHI resistance in natural field isolates identified a gene paralog of SDHC, termed ZtSDHC3, which encodes for an alternative C subunit of succinate dehydrogenase, named alt-SDHC. Using reverse genetics, we showed that alt-SDHC associates with the three other SDH subunits, leading to a fully functional enzyme and that a unique Qp-site residue within the alt-SDHC protein confers SHA-SDHI resistance. Enzymatic assays, computational modelling and docking simulations for the two SQR enzymes (altC-SQR, WT_SQR) enabled us to describe enzyme-inhibitor interactions at an atomistic level and to propose rational explanations for differential potency and resistance across SHA-SDHIs. European Z. tritici populations displayed a presence (20–30%) / absence polymorphism of ZtSDHC3, as well as differences in ZtSDHC3 expression levels and splicing efficiency. These polymorphisms have a strong impact on SHA-SDHI resistance phenotypes. Characterization of the ZtSDHC3 promoter in European Z. tritici populations suggests that transposon insertions are associated with the strongest resistance phenotypes. These results establish that a dispensable paralogous gene determines SHA-SDHIs fungicide resistance in natural populations of Z. tritici. This study paves the way to an increased awareness of the role of fungicidal target paralogs in resistance to fungicides and demonstrates the paramount importance of population genomics in fungicide discovery.

Zymoseptoria tritici is the causal agent of Septoria tritici leaf blotch (STB) of wheat, the most devastating disease for cereal production in Europe. Multiple succinate dehydrogenase inhibitor (SDHI) fungicides have been developed and introduced for the control of STB. We report the discovery and detailed characterization of a paralog of the C subunit of the SDH enzyme conferring standing resistance towards the SHA-SDHIs, a particular chemical subclass of the SDHIs. The SDHC paralog is characterized by its presence/absence, expression and alternative splicing polymorphisms, which in turn influence resistance levels. The identified mechanisms exemplify the importance of population genomics for the discovery and rational design of the most adapted solutions.

Genome-wide evolution of wobble base-pairing nucleotides of branchpoint motifs with increasing organismal complexity

How have the branchpoint motifs evolved in organisms of different complexity? Here we identified and examined the consensus motifs (R₁C₂T₃R₄A₅Y₆, R: A or G, Y: C or T) of 898 fungal genomes. In Ascomycota unicellular yeasts, the G₄/A₄ ratio is mostly (98%) below 0.125 but increases sharply in multicellular species by about 40 times on average, and in the more complex Basidiomycota, it increases further by about 7 times. The global G₄ increase is consistent with A₄ to G₄ transitions in evolution. Of the G₄/A₄-interacting amino acids of the branchpoint binding protein MSL5 (SF1) and the HSH155 (SF3B1), as well as the 5' splice sites (SS) and U2 snRNA genes, the 5' SS G³/A³ co-vary with the G₄ to some extent. However, corresponding increase of the G₄-complementary GCAGTA-U2 gene is rare, suggesting wobble-base pairing between the G₄-containing branchpoint motif and GTAGTA-U2 in most of these species. Interestingly, the G₄/A₄ ratio correlates well with the abundance of alternative splicing in the two phyla, and G₄ enriched significantly at the alternative 3' SS of genes in RNA metabolism, kinases and membrane proteins. Similar wobble nucleotides also enriched at the 3' SS of multicellular fungi with only thousands of protein-coding genes. Thus, branchpoint motifs have evolved U2-complementarity in unicellular Ascomycota yeasts, but have gradually gained more wobble base-pairing nucleotides in fungi of higher complexity, likely to destabilize branchpoint motif-U2 interaction and/or branchpoint A protrusion for alternative splicing. This implies an important role of relaxing the branchpoint signals in the multicellularity and further complexity of fungi.

Ascomycete Aspergillus oryzae Is an Efficient Expression Host for Production of Basidiomycete Terpenes by Using Genomic DNA Sequences

Basidiomycete fungi are an attractive resource for biologically active natural products for use in pharmaceutically relevant compounds. Recently, genome projects on mushroom fungi have provided a great deal of biosynthetic gene cluster information. However, functional analyses of the gene clusters for natural products were largely unexplored because of the difficulty of cDNA preparation and lack of gene manipulation tools for basidiomycete fungi. To develop a versatile host for basidiomycete genes, we examined gene expression using genomic DNA sequences in the robust ascomycete host Aspergillus oryzae, which is frequently used for the production of metabolites from filamentous fungi. Exhaustive expression of 30 terpene synthase genes from the basidiomycetes Clitopilus pseudo-pinsitus and Stereum hirsutum showed two splicing patterns, i.e., completely spliced cDNAs giving terpenes (15 cases) and mostly spliced cDNAs, indicating that A. oryzae correctly spliced most introns at the predicted positions and lengths. The mostly spliced cDNAs were expressed after PCR-based removal of introns, resulting in the successful production of terpenes (14 cases). During this study, we observed relatively frequent mispredictions in the automated program. Hence, the complementary use of A. oryzae expression and automated prediction will be a powerful tool for genome mining.IMPORTANCE The recent large influx of genome sequences from basidiomycetes, which are prolific producers of bioactive natural products, may provide opportunities to develop novel drug candidates. The development of a reliable expression system is essential for the genome mining of natural products because of the lack of a tractable host for heterologous expression of basidiomycete genes. For this purpose, we applied the ascomycete Aspergillus oryzae system for the direct expression of fungal natural product biosynthetic genes from genomic DNA. Using this system, 29 sesquiterpene synthase genes and diterpene biosynthetic genes for bioactive pleuromutilin were successfully expressed. Together with the use of computational tools for intron prediction, this Aspergillus oryzae system represents a practical method for the production of basidiomycete natural products.

PUF60-activated exons uncover altered 3' splice-site selection by germline missense mutations in a single RRM

PUF60 is a splicing factor that binds uridine (U)-rich tracts and facilitates association of the U2 small nuclear ribonucleoprotein with primary transcripts. PUF60 deficiency (PD) causes a developmental delay coupled with intellectual disability and spinal, cardiac, ocular and renal defects, but PD pathogenesis is not understood. Using RNA-Seq, we identify human PUF60-regulated exons and show that PUF60 preferentially acts as their activator. PUF60-activated internal exons are enriched for Us upstream of their 3′ splice sites (3′ss), are preceded by longer AG dinucleotide exclusion zones and more distant branch sites, with a higher probability of unpaired interactions across a typical branch site location as compared to control exons. In contrast, PUF60-repressed exons show U-depletion with lower estimates of RNA single-strandedness. We also describe PUF60-regulated, alternatively spliced isoforms encoding other U-bound splicing factors, including PUF60 partners, suggesting that they are co-regulated in the cell, and identify PUF60-regulated exons derived from transposed elements. PD-associated amino-acid substitutions, even within a single RNA recognition motif (RRM), altered selection of competing 3′ss and branch points of a PUF60-dependent exon and the 3′ss choice was also influenced by alternative splicing of PUF60. Finally, we propose that differential distribution of RNA processing steps detected in cells lacking PUF60 and the PUF60-paralog RBM39 is due to the RBM39 RS domain interactions. Together, these results provide new insights into regulation of exon usage by the 3′ss organization and reveal that germline mutation heterogeneity in RRMs can enhance phenotypic variability at the level of splice-site and branch-site selection.

A spliceosomal twin intron (stwintron) participates in both exon skipping and evolutionary exon loss

Spliceosomal twin introns (stwintrons) are introns where any of the three consensus sequences involved in splicing is interrupted by another intron (internal intron). In Aspergillus nidulans, a donor-disrupted stwintron (intron-1) is extant in the transcript encoding a reticulon-like protein. The orthologous transcript of Aspergillus niger can be alternatively spliced; the exon downstream the stwintron could be skipped by excising a sequence that comprises this stwintron, the neighbouring intron-2, and the exon bounded by these. This process involves the use of alternative 3' splice sites for the internal intron, the resulting alternative intervening sequence being a longer 3'-extended stwintron. In 29 species of Onygenales, a multi-step splicing process occurs in the orthologous transcript, in which a complex intervening sequence including the stwintron and neigbouring intron-2, generates by three splicing reactions a "second order intron" which must then be excised with a fourth splicing event. The gene model in two species can be envisaged as one canonical intron (intron-1) evolved from this complex intervening sequence of nested canonical introns found elsewhere in Onygenales. Postulated splicing intermediates were experimentally verified in one or more species. This work illustrates a role of stwintrons in both alternative splicing and the evolution of intron structure.

Identification and Functional Analysis of NLP-Encoding Genes from the Postharvest Pathogen Penicillium expansum

Penicillium expansum is a major postharvest pathogen that infects different fruits, mainly through injuries inflicted during harvest or subsequent handling after harvest. Several effectors were suggested to mediate pathogenicity of P. expansum in fruit tissue. Among these effectors Nep1-like proteins (NLPs), produced by various microorganisms with different lifestyles, are known for their ability to induce necrosis in dicot plants and were shown to be involved in virulence of several plant-related pathogens. This study was aimed at the identification and functional characterization of two NLP genes found in the genome of P. expansum. The genes were designated Penlp1 and Penlp2 and were found to code type1 and type3 NLP respectively. Necrosis-inducing activity of the two proteins was demonstrated by transient expression in Nicotiana benthamiana leaves. While Penlp1 expression was induced during apple infection and in liquid culture, the highest level of Penlp2 expression was found in ungerminated spores. Deletion of Penlp1, but not Penlp2, resulted in reduced virulence on apples manifested by reduced rate of lesion development (disease severity).