RNA-Seq analysis reveals new gene models and alternative splicing in the fungal pathogen Fusarium graminearum

Alternative splicing of the conserved drug-resistant orthologue FpNcb2 is associated with its nuclear accumulation of products and full virulence of Fusarium pseudograminearum

BACKGROUND

Negative cofactor 2 NC2β (Ncb2 or Dr1) is the beta subunit of a conserved heterodimeric regulator of transcription negative cofactor 2 (NC2) complex that has been identified as key regulator of drug resistance in model fungi. However, its role in plant pathogens is still unclear.

RESULTS

We identified an NC2β orthologue, FpNcb2, in Fusarium pseudograminearum, which is not only a significant regulatory function in drug resistance, but also essential for growth, conidiation and penetration. Moreover, FpNcb2 undergoes alternative splicing which creates two mRNA isoforms. As a putative CCAAT binding protein, FpNcb2 concentrates in the nuclei, contributing to the expression of two spliced mRNA of FpNcb2 in hypha, conidiophores and conidia, with exception of FpNcb2ISOA in germlings. Expression of each spliced mRNA of FpNcb2 in Δfpncb2 mutant could full complement the defects on growth, conidiation and fungicides sensitivity to that of wild type. However, FpNcb2ISOA and FpNcb2ISOB have different effects on virulence. FpNcb2 acts as a regulator for the transcription of some genes encoding drug efflux and hydrolases.

CONCLUSION

Our analysis showed the existence of alternative mRNA splicing in the NC2β orthologue, which is associated with protein subcellular localization and fungal virulence. The further elucidation of the target genes of NC2β will provide insights into the potential regulation mechanisms in the antifungal resistance and pathogenesis of F. pseudograminearum. © 2024 Society of Chemical Industry.

Secreted in Xylem 6 (SIX6) Mediates Fusarium oxysporum f. sp. fragariae Race 1 Avirulence on FW1-Resistant Strawberry Cultivars

Fusarium oxysporum f. sp. fragariae (Fof) race 1 is avirulent on cultivars with the dominant resistance gene FW1, while Fof race 2 is virulent on FW1-resistant cultivars. We hypothesized there was a gene-for-gene interaction between a gene at the FW1 locus and an avirulence gene (AvrFW1) in Fof race 1. To identify a candidate AvrFW1, we compared genomes of 24 Fof race 1 and three Fof race 2 isolates. We found one candidate gene that was present in race 1, was absent in race 2, was highly expressed in planta, and was homologous to a known effector, secreted in xylem 6 (SIX6). We knocked out SIX6 in two Fof race 1 isolates by homologous recombination. All SIX6 knockout transformants (ΔSIX6) gained virulence on FW1/fw1 cultivars, whereas ectopic transformants and the wildtype isolates remained avirulent. ΔSIX6 isolates were quantitatively less virulent on FW1/fw1 cultivars Fronteras and San Andreas than fw1/fw1 cultivars. Seedlings from an FW1/fw1 × fw1/fw1 population were genotyped for FW1 and tested for susceptibility to a SIX6 knockout isolate. Results suggested that additional minor-effect quantitative resistance genes could be present at the FW1 locus. This work demonstrates that SIX6 acts as an avirulence factor interacting with a resistance gene at the FW1 locus. The identification of AvrFW1 enables surveillance for Fof race 2 and provides insight into the mechanisms of FW1-mediated resistance. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Unanticipated Large-Scale Deletion in Fusarium graminearum Genome Using CRISPR/Cas9 and Its Impact on Growth and Virulence

Fusarium graminearum, a filamentous fungus, and causal agent of Fusarium head blight (FHB) in wheat and other cereals, leads to significant economic losses globally. This study aimed to investigate the roles of specific genes in F. graminearum virulence using CRISPR/Cas9-mediated gene deletions. Illumina sequencing was used to characterize the genomic changes due to editing. Unexpectedly, a large-scale deletion of 525,223 base pairs on chromosome 2, comprising over 222 genes, occurred in two isolates. Many of the deleted genes were predicted to be involved in essential molecular functions, such as oxidoreductase activity, transmembrane transporter activity, hydrolase activity, as well as biological processes, such as carbohydrate metabolism and transmembrane transport. Despite the substantial loss of genetic material, the mutant isolate exhibited normal growth rates and virulence on wheat under most conditions. However, growth rates were significantly reduced under high temperatures and on some media. Additionally, wheat inoculation assays using clip dipping, seed inoculation, and head point inoculation methods were performed. No significant differences in virulence were observed, suggesting that these genes were not involved in infection or alternative compensatory pathways, and allow the fungi to maintain pathogenicity despite the extensive genomic deletion.

New insights into genome annotation in Podospora anserina through re-exploiting multiple RNA-seq data

BACKGROUND

Publicly available RNA-seq datasets are often underused although being helpful to improve functional annotation of eukaryotic genomes. This is especially true for filamentous fungi genomes which structure differs from most well annotated yeast genomes. Podospora anserina is a filamentous fungal model, which genome has been sequenced and annotated in 2008. Still, the current annotation lacks information about cis-regulatory elements, including promoters, transcription starting sites and terminators, which are instrumental to integrate epigenomic features into global gene regulation strategies.

RESULTS

Here we took advantage of 37 RNA-seq experiments that were obtained in contrasted developmental and physiological conditions, to complete the functional annotation of P. anserina genome. Out of the 10,800 previously annotated genes, 5'UTR and 3'UTR were defined for 7554, among which, 3328 showed differential transcriptional signal starts and/or transcriptional end sites. In addition, alternative splicing events were detected for 2350 genes, mostly due alternative 3'splice sites and 1732 novel transcriptionally active regions (nTARs) in unannotated regions were identified.

CONCLUSIONS

Our study provides a comprehensive genome-wide functional annotation of P. anserina genome, including chromatin features, cis-acting elements such as UTRs, alternative splicing events and transcription of non-coding regions. These new findings will likely improve our understanding of gene regulation strategies in compact genomes, such as those of filamentous fungi. Characterization of alternative transcripts and nTARs paves the way to the discovery of putative new genes, alternative peptides or regulatory non-coding RNAs.

Relation between CarS expression and activation of carotenogenesis by stress in Fusarium fujikuroi

Fusarium fujikuroi, a model organism for secondary metabolism in fungi, produces carotenoids, terpenoid pigments with antioxidant activity. Previous results indicate that carotenoid synthesis in F. fujikuroi is stimulated by light or by different stress conditions and downregulated by a RING finger protein encoded by carS gene. Here, we have analyzed the effects of three stressors, nitrogen scarcity, heat shock, and oxidative stress. We compared them with the effect of light in the wild type, a carS mutant that overproduces carotenoids, and its complemented strain. The assayed stressors increase the synthesis of carotenoids in the three strains, but mRNA levels of structural genes of carotenogenesis, carRA and carB, are only enhanced in the presence of a functional carS gene. In the wild-type strain, the four conditions affect in different manners the mRNA levels of carS: greater in the presence of light, without significant changes in nitrogen starvation, and with patent decreases after heat shock or oxidative stress, suggesting different activation mechanisms. The spores of the carS mutant are more resistant to H₂O₂ than those of the wild type; however, the mutant shows a greater H₂O₂ sensitivity at the growth level, which may be due to the participation of CarS in the regulation of genes with catalase domains, formerly described. A possible mechanism of regulation by heat stress has been found in the alternative splicing of the intron of the carS gene, located close to its 3′ end, giving rise to the formation of a shorter protein. This action could explain the inducing effect of the heat shock, but not of the other inducing conditions, which may involve other mechanisms of action on the CarS regulator, either transcriptionally or post-transcriptionally.

Landscape and regulation of alternative splicing and alternative polyadenylation in a plant pathogenic fungus

Alternative splicing (AS) and alternative polyadenylation (APA) contribute significantly to the regulation of gene expression in higher eukaryotes. Their biological impact in filamentous fungi, however, is largely unknown. Here we combine PacBio Isoform-Sequencing and strand-specific RNA-sequencing of multiple tissues and mutant characterization to reveal the landscape and regulation of AS and APA in Fusarium graminearum. We generated a transcript annotation comprising 51 617 isoforms from 17 189 genes. In total, 4997 and 11 133 genes are alternatively spliced and polyadenylated, respectively. Majority of the AS events alter coding sequences. Unexpectedly, the AS transcripts containing premature-termination codons are not sensitive to nonsense-mediated messenger RNA decay. Unlike in yeasts and animals, distal APA sites have strong signals, but proximal APA isoforms are highly expressed in F. graminearum. The 3'-end processing factors FgRNA15, FgHRP1, and FgFIP1 play roles in promoting proximal APA site usage and intron splicing. A genome-wide increase in intron inclusion and distal APA site usage and downregulation of the spliceosomal and 3'-end processing factors were observed in older and quiescent tissues, indicating intron inclusion and 3'-untranslated region lengthening as novel mechanisms in regulating aging and dormancy in fungi. This study provides new insights into the complexity and regulation of AS and APA in filamentous fungi.

TuRLK1, a leucine-rich repeat receptor-like kinase, is indispensable for stripe rust resistance of YrU1 and confers broad resistance to multiple pathogens

BACKGROUND

YrU1 is a nucleotide-binding site (NBS) and leucine-rich repeat (LRR) protein (NLR), with additional ankyrin-repeat and WRKY domains and confers effective resistance to stripe rust fungus Puccinia striiformis f. sp. Tritici (Pst). YrU1 was positionally cloned in the progenitor species of the A genome of bread wheat, Tricicum urartu, recently. However, the molecular mechanism and components involved in YrU1-mediated resistance are not clear.

RESULTS

In this study, we found that the transcript level of TuRLK1, which encodes a novel leucine-rich repeat receptor-like kinase, was up-regulated after inoculation with Pst in the presence of YrU1, through RNA-seq analysis in T. urartu accession PI428309. TuRLK1 contained only a small number of LRR motifs, and was localized in the plasma-membrane. Transient expression of TuRLK1 induced hypersensitive cell death response in N. benthamiana leaves. Silencing of TuRLK1, using barley stripe mosaic virus (BSMV)-induced gene silencing (VIGS) system in PI428309 that contains YrU1, compromised the resistance against stripe rust caused by Pst CY33, indicating that TuRLK1 was required for YrU1-activated plant immunity. Furthermore, overexpression of TuRLK1 could enhance powdery mildew resistance in bread wheat and Arabidopsis thaliana after inoculating with the corresponding pathogens.

CONCLUSIONS

Our study indicates that TuRLK1 is required for immune response mediated by the unique NLR protein YrU1, and likely plays an important role in disease resistance to other pathogens.

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.

Transcriptional Landscapes of Long Non-coding RNAs and Alternative Splicing in Pyricularia oryzae Revealed by RNA-Seq

Pyricularia oryzae causes the rice blast, which is one of the most devastating crop diseases worldwide, and is a model fungal pathogen widely used for dissecting the molecular mechanisms underlying fungal virulence/pathogenicity. Although the whole genome sequence of P. oryzae is publicly available, its current transcriptomes remain incomplete, lacking the information on non-protein coding genes and alternative splicing. Here, we performed and analyzed RNA-Seq of conidia and hyphae, resulting in the identification of 3,374 novel genes. Interestingly, the vast majority of these novel genes likely transcribed long non-coding RNAs (lncRNAs), and most of them were localized in the intergenic regions. Notably, their expressions were concomitant with the transcription of neighboring genes thereof in conidia and hyphae. In addition, 2,358 genes were found to undergo alternative splicing events. Furthermore, we exemplified that a lncRNA was important for hyphal growth likely by regulating the neighboring protein-coding gene and that alternative splicing of the transcription factor gene CON7 was required for appressorium formation. In summary, results from this study indicate that lncRNA transcripts and alternative splicing events are two important mechanisms for regulating the expression of genes important for conidiation, hyphal growth, and pathogenesis, and provide new insights into transcriptomes and gene regulation in the rice blast fungus.

Identification of Putative Virulence Genes by DNA Methylation Studies in the Cereal Pathogen Fusarium graminearum

DNA methylation mediates organisms’ adaptations to environmental changes in a wide range of species. We investigated if a such a strategy is also adopted by Fusarium graminearum in regulating virulence toward its natural hosts. A virulent strain of this fungus was consecutively sub-cultured for 50 times (once a week) on potato dextrose agar. To assess the effect of subculturing on virulence, wheat seedlings and heads (cv. A416) were inoculated with subcultures (SC) 1, 23, and 50. SC50 was also used to re-infect (three times) wheat heads (SC50×3) to restore virulence. In vitro conidia production, colonies growth and secondary metabolites production were also determined for SC1, SC23, SC50, and SC50×3. Seedling stem base and head assays revealed a virulence decline of all subcultures, whereas virulence was restored in SC50×3. The same trend was observed in conidia production. The DNA isolated from SC50 and SC50×3 was subject to a methylation content-sensitive enzyme and double-digest, restriction-site-associated DNA technique (ddRAD-MCSeEd). DNA methylation analysis indicated 1024 genes, whose methylation levels changed in response to the inoculation on a healthy host after subculturing. Several of these genes are already known to be involved in virulence by functional analysis. These results demonstrate that the physiological shifts following sub-culturing have an impact on genomic DNA methylation levels and suggest that the ddRAD-MCSeEd approach can be an important tool for detecting genes potentially related to fungal virulence.

The RNA binding protein FgRbp1 regulates specific pre-mRNA splicing via interacting with U2AF23 in Fusarium

Precursor messenger RNA (pre-mRNA) splicing is an essential and tightly regulated process in eukaryotic cells; however, the regulatory mechanisms for the splicing are not well understood. Here, we characterize a RNA binding protein named FgRbp1 in Fusarium graminearum, a fungal pathogen of cereal crops worldwide. Deletion of FgRbp1 leads to reduced splicing efficiency in 47% of the F. graminearum intron-containing gene transcripts that are involved in various cellular processes including vegetative growth, development, and virulence. The human ortholog RBM42 is able to fully rescue the growth defects of ΔFgRbp1. FgRbp1 binds to the motif CAAGR in its target mRNAs, and interacts with the splicing factor FgU2AF23, a highly conserved protein involved in 3' splice site recognition, leading to enhanced recruitment of FgU2AF23 to the target mRNAs. This study demonstrates that FgRbp1 is a splicing regulator and regulates the pre-mRNA splicing in a sequence-dependent manner in F. graminearum.

A conserved Zn2Cys6 transcription factor, identified in a spontaneous mutant from in vitro passaging, is involved in pathogenicity of the blackleg fungus Leptosphaeria maculans

Continuous passaging in vitro can lead to the accumulation of changes in DNA sequence that potentially affect the properties of microbes, making them different from the original isolates. The identification of such genetic alterations is rare in fungi. A set of insertional mutants in the plant pathogenic fungus Leptosphaeria maculans, all derived from the same transformation experiment, had independent Agrobacterium T-DNA insertions and reduced pathogenicity on canola (Brassica napus). None of the insertions co-segregated in progeny from crosses with the reduction in pathogenicity. Genome sequences of three strains were analysed, and a mutation identified in a gene (ptf1, for pathogenicity-associated transcription factor 1) encoding a putative Zn₂(II)Cys₆ transcription factor. Homologs are found in other ascomycetes, and are required for pathogenicity by Fusarium graminearum, Fusarium oxysporum and Magnaporthe oryzae. The mutation in the L. maculans ptf1 gene co-segregates in progeny from crosses with the reduction in pathogenicity, a strain with an independent mutant allele isolated using CRISPR-Cas9 editing has reduced pathogenicity, and addition of wild type copies of the gene restores pathogenicity. Thus, this work defines a base pair substitution that occurred during in vitro passaging of a fungus that contributed to an attenuation of pathogenicity.

Genome-wide alternative splicing profiling in the fungal plant pathogen Sclerotinia sclerotiorum during the colonization of diverse host families

Sclerotinia sclerotiorum is a notorious generalist plant pathogen that threatens more than 600 host plants, including wild and cultivated species. The molecular bases underlying the broad compatibility of S. sclerotiorum with its hosts is not fully elucidated. In contrast to higher plants and animals, alternative splicing (AS) is not well studied in plant-pathogenic fungi. AS is a common regulated cellular process that increases cell protein and RNA diversity. In this study, we annotated spliceosome genes in the genome of S. sclerotiorum and characterized their expression in vitro and during the colonization of six host species. Several spliceosome genes were differentially expressed in planta, suggesting that AS was altered during infection. Using stringent parameters, we identified 1,487 S. sclerotiorum genes differentially expressed in planta and exhibiting alternative transcripts. The most common AS events during the colonization of all plants were retained introns and the alternative 3' receiver site. We identified S. sclerotiorum genes expressed in planta for which (a) the relative accumulation of alternative transcripts varies according to the host being colonized and (b) alternative transcripts harbour distinct protein domains. This notably included 42 genes encoding predicted secreted proteins showing high-confidence AS events. This study indicates that AS events are taking place in the plant pathogenic fungus S. sclerotiorum during the colonization of host plants and could generate functional diversity in the repertoire of proteins secreted by S. sclerotiorum during infection.

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.

Transcriptomes of the interaction between Fusarium verticillioides and a Streptomyces strain reveal the fungal defense strategy under the pressure of a potential biocontrol agent

The actinobacteria Streptomyces sp. AV05 appears to be a potential biocontrol agent (BCA) against mycotoxigenic fungi. It was found to significantly inhibit F. verticillioides growth and mycotoxin production during their co-cultivation. F. verticillioides growth was durably affected while the decrease of the toxin production levels was reversible, suggesting different BCA actions. The study of both transcriptomes brought useful information on the microbial interaction. RNA-seq data indicated that the dual interaction modified genetic expression of both microorganisms as 18.5 % of the genes were differentially expressed for the fungus against 3.8 % for the actinobacteria. Fungal differentially expressed genes (DEGs) were equally up and down regulated while bacterial ones were mainly upregulated. We especially focused the analysis of DEGs on fungal defense reaction to bacterial attack. For example, if this potential BCA implements a strategy of antibiosis with the over expression of 'siderophore-interacting protein' linked to the production of bacteriocins, the fungus in a state of stress is able to adapt its metabolism by up-regulation of amidase. It could correspond to the induction of resistance gene clusters and suggest a detoxification process. Moreover fumonisins-related pathway appears underexpressed in the presence of Streptomyces that explain the reduction of fumonisin accumulation observed.

Global identification of alternative splicing in Shiraia bambusicola and analysis of its regulation in hypocrellin biosynthesis

Hypocrellins, as natural pigments from Shiraia bambusicola, have extensive applications in the agricultural, cosmetic, food, and feed industries, and play a vital role in photodynamic therapy for anticancer and antiviral treatments. However, environmental stresses are always the bottlenecks for increasing hypocrellin yield during the process of fermentation. Pre-mRNA alternative splicing (AS) is an essential mechanism in the defense of abiotic stresses in the animal and plant kingdom, but is seldom involved in fungi. In this study, AS from genome-wide sequencing and RNA-seq data for S. bambusicola was analyzed for the first time. Interestingly, the proportion of AS in S. bambusicola was 38.44% (most of them participated in metabolic processes, covering pigmentation and response to stimulus), a much higher ratio than seen in that of other fungal species (1.3-18%). Here, we identified the relationship of AS and secondary metabolic (SM) biosynthesis under a series of abiotic stresses. Suitable fungicides suppressed hypocrellin production significantly, and AS occurred in key functional genes (sbFLO, sbMFS, sbPKS) of hypocrellin biosynthesis. In contrast, H₂O₂ improved the yield of hypocrellins, but AS were not found in the corresponding gene cluster. A further study showed that overexpressing an isoform of sbPKS (sbPKSa) in Shiraia bambusicola could dramatically down-regulate the expression of the original gene sbPKS and nearly inhibit the production of hypocrellins. Altogether, our study strongly supported the hypothesis that AS had a vital role in the regulation of hypocrellin biosynthesis under stresses, and initially explored whether SM functional genes were relevant for fungi.

Alternative Splicing in Heat Shock Protein Transcripts as a Mechanism of Cell Adaptation in Trichophyton rubrum

Heat shock proteins (HSPs) are involved in critical processes like host tissue invasion, resistance, and pathogenicity in dermatophytes. RNA-Seq analysis of Trichophyton rubrum exposed to undecanoic acid (UDA) revealed intron retention events in HSP transcripts. Because HSPs are modulated in response to various stimuli and as alternative splicing (AS) can result in a broad diversity in the proteome of eukaryotic cells, our objective was to confirm the aforementioned retention events, investigating their consequences and extent. Furthermore, we aimed to determine: (1) the expression profile of HSP genes in an infection-like scenario and (2) the importance of Hsp90 for the keratinolytic potential of T. rubrum. RT and qPCR analyses comparing the exposure to UDA and terbinafine (TRB) confirmed the presence of two mRNA isoforms of the hsp7-like gene, with distinct expression patterns in response to UDA and TRB. The HSP expression profile revealed two upregulated, three downregulated, and four unmodulated transcripts; Hsp90 inhibition by 17-AAG resulted in a significant decrease in keratinolytic potential at 37 °C. Altogether, these results broaden the current knowledge on the importance of HSP-mediated pathways for cell adaptation and other aspects of dermatophyte biology, indicating that HSP network proteins can be potential targets for antifungal therapy.

Biochemical Characterization of the Fusarium graminearum Candidate ACC-Deaminases and Virulence Testing of Knockout Mutant Strains

Fusarium graminearum is a plant pathogenic fungus which is able to infect wheat and other economically important cereal crop species. The role of ethylene in the interaction with host plants is unclear and controversial. We have analyzed the inventory of genes with a putative function in ethylene production or degradation of the ethylene precursor 1-aminocyclopropane carboxylic acid (ACC). F. graminearum, in contrast to other species, does not contain a candidate gene encoding ethylene-forming enzyme. Three genes with similarity to ACC synthases exist; heterologous expression of these did not reveal enzymatic activity. The F. graminearum genome contains in addition two ACC deaminase candidate genes. We have expressed both genes in E. coli and characterized the enzymatic properties of the affinity-purified products. One of the proteins had indeed ACC deaminase activity, with kinetic properties similar to ethylene-stress reducing enzymes of plant growth promoting bacteria. The other candidate was inactive with ACC but turned out to be a d-cysteine desulfhydrase. Since it had been reported that ethylene insensitivity in transgenic wheat increased Fusarium resistance and reduced the content of the mycotoxin deoxynivalenol (DON) in infected wheat, we generated single and double knockout mutants of both genes in the F. graminearum strain PH-1. No statistically significant effect of the gene disruptions on fungal spread or mycotoxin content was detected, indicating that the ability of the fungus to manipulate the production of the gaseous plant hormones ethylene and H₂S is dispensable for full virulence.

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).

Functional characterization of host toxic EcdB transcription factor protein of echinocandin B biosynthetic gene cluster

The ecdB is a transcription factor, located in the echinocandin B biosynthetic gene cluster of Emericella rugulosa NRRL11440. Here, we validated the ecdB mRNA sequence for functional expression and to explore the role of EcdB protein in the echinocandin B regulation. The sequence alignment study revealed that the ecdB coding sequence was found 75 bp shorter than the reference mRNA sequence. This coding sequence encodes for EcdB protein and comprises three conserved domains; DNA binding domain (DBD), coiled-coil domain, and signature middle homology region. The full-length and DBD (truncated) DNA sequences were expressed in Escherichia coli BL21(DE3) under different tested conditions. The expression of EcdB protein was found to be toxic, which curbs the cell growth. In contrast to truncated protein (GST:EcdB1-54), the full-length (GST:EcdB) protein was expressed at very low titer and not detectable in SDS-PAGE under the varying isopropyl β-d-1-thiogalactopyranoside (IPTG), temperature, and media conditions. However, GST:EcdB1-54 was successfully purified under standard conditions (0.5 mM IPTG at 0.5OD) with 33 kDa expected size. The functionality of GST:EcdB1-54 was attained by electrophoretic mobility shift assay study as a clear band shifting showed with ecdA promoter. Taken together, we conclude that EcdB interacts with the ecdA promoter that reflected to require for echinocandin B regulation.

Study of the whole genome, methylome and transcriptome of Cordyceps militaris

The complete genome of Cordyceps militaris was sequenced using single-molecule real-time (SMRT) sequencing technology at a coverage over 300×. The genome size was 32.57 Mb, and 14 contigs ranging from 0.35 to 4.58 Mb with an N50 of 2.86 Mb were assembled, including 4 contigs with telomeric sequences on both ends and an additional 8 contigs with telomeric sequences on either the 5′ or 3′ end. A methylome database of the genome was constructed using SMRT and m4C and m6A methylated nucleotides, and many unknown modification types were identified. The major m6A methylation motif is GA and GGAG, and the major m4C methylation motif is GC or CG/GC. In the C. militaris genome DNA, there were four types of methylated nucleotides that we confirmed using high-resolution LCMS-IT-TOF. Using PacBio Iso-Seq, a total of 31,133 complete cDNA sequences were obtained in the fruiting body. The conserved domains of the nontranscribed regions of the genome include TATA boxes, which are the initial regions of genome replication. There were 406 structural variants between the HN and CM01 strains, and there were 1,114 structural variants between the HN and ATCC strains.

Opposite Polarity Monospore Genome De Novo Sequencing and Comparative Analysis Reveal the Possible Heterothallic Life Cycle of Morchella importuna

Morchella is a popular edible fungus worldwide due to its rich nutrition and unique flavor. Many research efforts were made on the domestication and cultivation of Morchella all over the world. In recent years, the cultivation of Morchella was successfully commercialized in China. However, the biology is not well understood, which restricts the further development of the morel fungus cultivation industry. In this paper, we performed de novo sequencing and assembly of the genomes of two monospores with a different mating type (M04M24 and M04M26) isolated from the commercially cultivated strain M04. Gene annotation and comparative genome analysis were performed to study differences in CAZyme (Carbohydrate-active enzyme) enzyme content, transcription factors, duplicated sequences, structure of mating type sites, and differences at the gene and functional levels between the two monospore strains of M. importuna. Results showed that the de novo assembled haploid M04M24 and M04M26 genomes were 48.98 and 51.07 Mb, respectively. A complete fine physical map of M. importuna was obtained from genome coverage and gene completeness evaluation. A total of 10,852 and 10,902 common genes and 667 and 868 endemic genes were identified from the two monospore strains, respectively. The Gene Ontology (GO) and KAAS (KEGG Automatic Annotation Serve) enrichment analyses showed that the endemic genes performed different functions. The two monospore strains had 99.22% collinearity with each other, accompanied with certain position and rearrangement events. Analysis of complete mating-type loci revealed that the two monospore M. importuna strains contained an independent mating-type structure and remained conserved in sequence and location. The phylogenetic and divergence time of M. importuna was analyzed at the whole-genome level for the first time. The bifurcation time of morel and tuber was estimated to be 201.14 million years ago (Mya); the two monospore strains with a different mating type represented the evolution of different nuclei, and the single copy homologous genes between them were also different due to a genetic differentiation distance about 0.65 Mya. Compared with truffles, M. importuna had an extension of 28 clusters of orthologous genes (COGs) and a contraction of two COGs. The two different polar nuclei with different degrees of contraction and expansion suggested that they might have undergone different evolutionary processes. The different mating-type structures, together with the functional clustering and enrichment analysis results of the endemic genes of the two different polar nuclei, imply that M. importuna might be a heterothallic fungus and the interaction between the endemic genes may be necessary for its complete life history. Studies on the genome of M. importuna facilitate a better understanding of morel biology and evolution.

The novel cyst nematode effector protein 30D08 targets host nuclear functions to alter gene expression in feeding sites

Cyst nematodes deliver effector proteins into host cells to manipulate cellular processes and establish a metabolically hyperactive feeding site. The novel 30D08 effector protein is produced in the dorsal gland of parasitic juveniles, but its function has remained unknown. We demonstrate that expression of 30D08 contributes to nematode parasitism, the protein is packaged into secretory granules and it is targeted to the plant nucleus where it interacts with SMU2 (homolog of suppressor of mec-8 and unc-52 2), an auxiliary spliceosomal protein. We show that SMU2 is expressed in feeding sites and an smu2 mutant is less susceptible to nematode infection. In Arabidopsis expressing 30D08 under the SMU2 promoter, several genes were found to be alternatively spliced and the most abundant functional classes represented among differentially expressed genes were involved in RNA processing, transcription and binding, as well as in development, and hormone and secondary metabolism, representing key cellular processes known to be important for feeding site formation. In conclusion, we demonstrated that the 30D08 effector is secreted from the nematode and targeted to the plant nucleus where its interaction with a host auxiliary spliceosomal protein may alter the pre-mRNA splicing and expression of a subset of genes important for feeding site formation.

Evolution and Diversity of Biosynthetic Gene Clusters in Fusarium

Plant pathogenic fungi in the Fusarium genus cause severe damage to crops, resulting in great financial losses and health hazards. Specialized metabolites synthesized by these fungi are known to play key roles in the infection process, and to provide survival advantages inside and outside the host. However, systematic studies of the evolution of specialized metabolite-coding potential across Fusarium have been scarce. Here, we apply a combination of bioinformatic approaches to identify biosynthetic gene clusters (BGCs) across publicly available genomes from Fusarium, to group them into annotated families and to study gain/loss events of BGC families throughout the history of the genus. Comparison with MIBiG reference BGCs allowed assignment of 29 gene cluster families (GCFs) to pathways responsible for the production of known compounds, while for 57 GCFs, the molecular products remain unknown. Comparative analysis of BGC repertoires using ancestral state reconstruction raised several new hypotheses on how BGCs contribute to Fusarium pathogenicity or host specificity, sometimes surprisingly so: for example, a gene cluster for the biosynthesis of hexadehydro-astechrome was identified in the genome of the biocontrol strain Fusarium oxysporum Fo47, while being absent in that of the tomato pathogen F. oxysporum f.sp. lycopersici. Several BGCs were also identified on supernumerary chromosomes; heterologous expression of genes for three terpene synthases encoded on the Fusarium poae supernumerary chromosome and subsequent GC/MS analysis showed that these genes are functional and encode enzymes that each are able to synthesize koraiol; this observed functional redundancy supports the hypothesis that localization of copies of BGCs on supernumerary chromosomes provides freedom for evolutionary innovations to occur, while the original function remains conserved. Altogether, this systematic overview of biosynthetic diversity in Fusarium paves the way for targeted natural product discovery based on automated identification of species-specific pathways as well as for connecting species ecology to the taxonomic distributions of BGCs.

The NB-LRR gene Pm60 confers powdery mildew resistance in wheat

Powdery mildew is one of the most devastating diseases of wheat. To date, few powdery mildew resistance genes have been cloned from wheat due to the size and complexity of the wheat genome. Triticum urartu is the progenitor of the A genome of wheat and is an important source for powdery mildew resistance genes. Using molecular markers designed from scaffolds of the sequenced T. urartu accession and standard map-based cloning, a powdery mildew resistance locus was mapped to a 356-kb region, which contains two nucleotide-binding and leucine-rich repeat domain (NB-LRR) protein-encoding genes. Virus-induced gene silencing, single-cell transient expression, and stable transformation assays demonstrated that one of these two genes, designated Pm60, confers resistance to powdery mildew. Overexpression of full-length Pm60 and two allelic variants in Nicotiana benthamiana leaves induced hypersensitive cell death response, but expression of the coiled-coil domain alone was insufficient to induce hypersensitive response. Yeast two-hybrid, bimolecular fluorescence complementation and luciferase complementation imaging assays showed that Pm60 protein interacts with its neighboring NB-containing protein, suggesting that they might be functionally related. The identification and cloning of this novel wheat powdery mildew resistance gene will facilitate breeding for disease resistance in wheat.

The matrices and constraints of GT/AG splice sites of more than 1000 species/lineages

To provide a resource for the splice sites (SS) of different species, we calculated the matrices of nucleotide compositions of about 38 million splice sites from >1000 species/lineages. The matrices are enriched of aGGTAAGT (5'SS) or (Y)₆N(C/t)AG(g/a)t (3'SS) overall; however, they are quite diverse among hundreds of species. The diverse matrices remain prominent even under sequence selection pressures, suggesting the existence of diverse constraints as well as U snRNAs and other spliceosomal factors and/or their interactions with the splice sites. Using an algorithm to measure and compare the splice site constraints across all species, we demonstrate their distinct differences quantitatively. As an example of the resource's application to answering specific questions, we confirm that high constraints of particular positions are significantly associated with transcriptome-wide, increased occurrences of alternative splicing when uncommon nucleotides are present. More interestingly, the abundance of alternative splicing in 16 species correlates with the average constraint index of splice sites in a bell curve. This resource will allow users to assess specific sequences/splice sites against the consensus of every Ensembl-annotated species, and to explore the evolutionary changes or relationship to alternative splicing and transcriptome diversity. Web-search or update features are also included.

Transcriptome analysis of avian reovirus-mediated changes in gene expression of normal chicken fibroblast DF-1 cells

Background

Avian reovirus (ARV) is an important poultry pathogen that can cause immunosuppression. In this study, RNA-Seq technology was applied to investigate the transcriptome-wide changes of DF-1 cells upon ARV infection at the middle stage.

Results

Total RNA of ARV-infected or mock-infected samples at 10 and 18 h post infection (hpi) was extracted to build RNA-Seq datasets. Analysis of the sequencing data revealed that the expressions of numerous genes were altered, and a panel of differentially expressed genes were confirmed with RT-qPCR. At 10 hpi, 104 genes were down-regulated and 64 were up-regulated, while the expressions of 47 genes were increased and only one was down-regulated, which may play a role in retinoic acid biosynthesis, at 18 hpi in the ARV-infected cells. The similar profiles of up-regulated genes between the two groups of infected cells suggest that ARV infection activated a prolonged antiviral response of host cells. Alternative splicing analysis found no significantly changed events altered by ARV infection.

Conclusions

Overall, the differential expression profile presented in this study can be used to expand our understanding of the comprehensive interactions between ARV and the host cells, and may be helpful for us to reveal the pathogenic mechanism on the molecular level.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4310-5) contains supplementary material, which is available to authorized users.

Architecture and Distribution of Introns in Core Genes of Four Fusarium Species

Removal of introns from transcribed RNA represents a crucial step during the production of mRNA in eukaryotes. Available whole-genome sequences and expressed sequence tags (ESTs) have increased our knowledge of this process and revealed various commonalities among eukaryotes. However, certain aspects of intron structure and diversity are taxon-specific, which can complicate the accuracy of in silico gene prediction methods. Using core genes, we evaluated the distribution and architecture of Fusarium circinatum spliceosomal introns, and linked these characteristics to the accuracy of the predicted gene models of the genome of this fungus. We also evaluated intron distribution and architecture in F. verticillioides, F. oxysporum, and F. graminearum, and made comparisons with F. circinatum. Results indicated that F. circinatum and the three other Fusarium species have canonical 5′ and 3′ splice sites, but with subtle differences that are apparently not shared with those of other fungal genera. The polypyrimidine tract of Fusarium introns was also found to be highly divergent among species and genes. Furthermore, the conserved adenosine nucleoside required during the first step of splicing is contained within unique branch site motifs in certain Fusarium introns. Data generated here show that introns of F. circinatum, as well as F. verticillioides, F. oxysporum, and F. graminearum, are characterized by a number of unique features such as the CTHAH and ACCAT motifs of the branch site. Incorporation of such information into genome annotation software will undoubtedly improve the accuracy of gene prediction methods used for Fusarium species and related fungi.

Transcriptome analysis reveals the host selection fitness mechanisms of the Rhizoctonia solani AG1IA pathogen

Rhizoctonia solani AG1IA is a major generalist pathogen that causes sheath blight. Its genome, which was the first to be sequenced from the Rhizoctonia genus, may serve as a model for studying pathogenic mechanisms. To explore the pathogen-host fitness mechanism of sheath-blight fungus, a comprehensive comparative transcriptome ecotype analysis of R. solani AG1IA isolated from rice, soybean and corn during infection was performed. Special characteristics in gene expression, gene ontology terms and expression of pathogenesis-associated genes, including genes encoding secreted proteins, candidate effectors, hydrolases, and proteins involved in secondary metabolite production and the MAPK pathway, were revealed. Furthermore, as an important means of pathogenic modulation, diverse alternative splicing of key pathogenic genes in Rhizoctonia solani AG1IA during infections of the abovementioned hosts was uncovered for the first time. These important findings of key factors in the pathogenicity of R. solani AG1IA ecotypes during infection of various hosts explain host preference and provide novel insights into the pathogenic mechanisms and host-pathogen selection. Furthermore, they provide information on the fitness of Rhizoctonia, a severe pathogen with a wide host range.

Transcriptional Analysis Allows Genome Reannotation and Reveals that Cryptococcus gattii VGII Undergoes Nutrient Restriction during Infection

Cryptococcus gattii is a human and animal pathogen that infects healthy hosts and caused the Pacific Northwest outbreak of cryptococcosis. The inhalation of infectious propagules can lead to internalization of cryptococcal cells by alveolar macrophages, a niche in which C. gattii cells can survive and proliferate. Although the nutrient composition of macrophages is relatively unknown, the high induction of amino acid transporter genes inside the phagosome indicates a preference for amino acid uptake instead of synthesis. However, the presence of countable errors in the R265 genome annotation indicates significant inhibition of transcriptomic analysis in this hypervirulent strain. Thus, we analyzed RNA-Seq data from in vivo and in vitro cultures of C. gattii R265 to perform the reannotation of the genome. In addition, based on in vivo transcriptomic data, we identified highly expressed genes and pathways of amino acid metabolism that would enable C. gattii to survive and proliferate in vivo. Importantly, we identified high expression in three APC amino acid transporters as well as the GABA permease. The use of amino acids as carbon and nitrogen sources, releasing ammonium and generating carbohydrate metabolism intermediaries, also explains the high expression of components of several degradative pathways, since glucose starvation is an important host defense mechanism.

Fungal Plant Pathogenesis Mediated by Effectors

The interactions between fungi and plants encompass a spectrum of ecologies ranging from saprotrophy (growth on dead plant material) through pathogenesis (growth of the fungus accompanied by disease on the plant) to symbiosis (growth of the fungus with growth enhancement of the plant). We consider pathogenesis in this article and the key roles played by a range of pathogen-encoded molecules that have collectively become known as effectors.

Candidate gene based association mapping in Fusarium culmorum for field quantitative pathogenicity and mycotoxin production in wheat

Background

Quantitative traits are common in nature, but quantitative pathogenicity has received only little attention in phytopathology. In this study, we used 100 Fusarium culmorum isolates collected from natural field environments to assess their variation for two quantitative traits, aggressiveness and deoxynivalenol (DON) production on wheat plants grown in four different field environments (location-year combinations). Seventeen Fusarium graminearum pathogenicity candidate genes were assessed for their effect on the aggressiveness and DON production of F. culmorum under field conditions.

Results

For both traits, genotypic variance among isolates was high and significant while the isolate-by-environment interaction was also significant, amounting to approximately half of the genotypic variance. Among the studied candidate genes, the mitogen-activated protein kinase (MAPK) HOG1 was found to be significantly associated with aggressiveness and deoxynivalenol (DON) production, explaining 10.29 and 6.05% of the genotypic variance, respectively.

Conclusions

To the best of our knowledge, this is the first report of a protein kinase regulator explaining differences in field aggressiveness and mycotoxin production among individuals from natural populations of a plant pathogen.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0511-9) contains supplementary material, which is available to authorized users.

Transcriptome analysis of smut fungi reveals widespread intergenic transcription and conserved antisense transcript expression

Background

Biotrophic fungal plant pathogens cause billions of dollars in losses to North American crops annually. The model for functional investigation of these fungi is Ustilago maydis. Its 20.5 Mb annotated genome sequence has been an excellent resource for investigating biotrophic plant pathogenesis. Expressed-sequence tag libraries and microarray hybridizations have provided insight regarding the type of transcripts produced by U. maydis but these analyses were not comprehensive and there were insufficient data for transcriptome comparison to other smut fungi. To improve transcriptome annotation and enable comparative analyses, comprehensive strand-specific RNA-seq was performed on cell-types of three related smut species: U. maydis (common smut of corn), Ustilago hordei (covered smut of barley), and Sporisorium reilianum (head smut of corn).

Results

In total, >1 billion paired-end sequence reads were obtained from haploid cell, dikaryon and teliospore RNA of U. maydis, haploid cell RNA of U. hordei, and haploid and dikaryon cell RNA of S. reilianum. The sequences were assembled into transfrags using Trinity, and updated gene models were created using PASA and categorized with Cufflinks Cuffcompare. Representative genes that were predicted for the first time with these RNA-seq analyses and genes with novel annotation features were independently assessed by reverse transcriptase PCR. The analyses indicate hundreds more predicted proteins, relative to the previous genome annotation, could be produced by U. maydis from altered transcript forms, and that the number of non-coding RNAs produced, including transcribed intergenic sequences and natural antisense transcripts, approximately equals the number of mRNAs. This high representation of non-coding RNAs appears to be a conserved feature of the smut fungi regardless of whether they have RNA interference machinery. Approximately 50% of the identified NATs were conserved among the smut fungi.

Conclusions

Overall, these analyses revealed: 1) smut genomes encode a number of transcriptional units that is twice the number of annotated protein-coding genes, 2) a small number of intergenic transcripts may encode proteins with characteristics of fungal effectors, 3) the vast majority of intergenic and antisense transcripts do not contain ORFs, 4) a large proportion of the identified antisense transcripts were detected at orthologous loci among the smut fungi, and 5) there is an enrichment of functional categories among orthologous loci that suggests antisense RNAs could have a genome-wide, non-RNAi-mediated, influence on gene expression in smut fungi.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3720-8) contains supplementary material, which is available to authorized users.

Transcriptome analysis reveals the complexity of alternative splicing regulation in the fungus Verticillium dahliae

Background

Alternative splicing (AS) regulation is extensive and shapes the functional complexity of higher organisms. However, the contribution of alternative splicing to fungal biology is not well studied.

Results

This study provides sequences of the transcriptomes of the plant wilt pathogen Verticillium dahliae, using two different strains and multiple methods for cDNA library preparations. We identified alternatively spliced mRNA isoforms in over a half of the multi-exonic fungal genes. Over one-thousand isoforms involve TopHat novel splice junction; multiple types of combinatory alternative splicing patterns were identified. We showed that one Verticillium gene could use four different 5′ splice sites and two different 3′ donor sites to produce up to five mature mRNAs, representing one of the most sophisticated alternative splicing model in eukaryotes other than animals. Hundreds of novel intron types involving a pair of new splice sites were identified in the V. dahliae genome. All the types of AS events were validated by using RT-PCR. Functional enrichment analysis showed that AS genes are involved in most known biological functions and enriched in ATP biosynthesis, sexual/asexual reproduction, morphogenesis, signal transduction etc., predicting that the AS regulation modulates mRNA isoform output and shapes the V. dahliae proteome plasticity of the pathogen in response to the environmental and developmental changes.

Conclusions

These findings demonstrate the comprehensive alternative splicing mechanisms in a fungal plant pathogen, which argues the importance of this fungus in developing complicate genome regulation strategies in eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3507-y) contains supplementary material, which is available to authorized users.

Whole genome sequencing and comparative genomics of closely related Fusarium Head Blight fungi: Fusarium graminearum, F. meridionale and F. asiaticum

BACKGROUND

The Fusarium graminearum species complex is composed of many distinct fungal species that cause several diseases in economically important crops, including Fusarium Head Blight of wheat. Despite being closely related, these species and individuals within species have distinct phenotypic differences in toxin production and pathogenicity, with some isolates reported as non-pathogenic on certain hosts. In this report, we compare genomes and gene content of six new isolates from the species complex, including the first available genomes of F. asiaticum and F. meridionale, with four other genomes reported in previous studies.

RESULTS

A comparison of genome structure and gene content revealed a 93-99% overlap across all ten genomes. We identified more than 700 k base pairs (kb) of single nucleotide polymorphisms (SNPs), insertions, and deletions (indels) within common regions of the genome, which validated the species and genetic populations reported within species. We constructed a non-redundant pan gene list containing 15,297 genes from the ten genomes and among them 1827 genes or 12% were absent in at least one genome. These genes were co-localized in telomeric regions and select regions within chromosomes with a corresponding increase in SNPs and indels. Many are also predicted to encode for proteins involved in secondary metabolism and other functions associated with disease. Genes that were common between isolates contained high levels of nucleotide variation and may be pseudogenes, allelic, or under diversifying selection.

CONCLUSIONS

The genomic resources we have contributed will be useful for the identification of genes that contribute to the phenotypic variation and niche specialization that have been reported among members of the F. graminearum species complex.

Schizophyllum commune has an extensive and functional alternative splicing repertoire

Recent genome-wide studies have demonstrated that fungi possess the machinery to alternatively splice pre-mRNA. However, there has not been a systematic categorization of the functional impact of alternative splicing in a fungus. We investigate alternative splicing and its functional consequences in the model mushroom forming fungus Schizophyllum commune. Alternative splicing was demonstrated for 2,285 out of 12,988 expressed genes, resulting in 20% additional transcripts. Intron retentions were the most common alternative splicing events, accounting for 33% of all splicing events, and 43% of the events in coding regions. On the other hand, exon skipping events were rare in coding regions (1%) but enriched in UTRs where they accounted for 57% of the events. Specific functional groups, including transcription factors, contained alternatively spliced genes. Alternatively spliced transcripts were regulated differently throughout development in 19% of the 2,285 alternatively spliced genes. Notably, 69% of alternatively spliced genes have predicted alternative functionality by loss or gain of functional domains, or by acquiring alternative subcellular locations. S. commune exhibits more alternative splicing than any other studied fungus. Taken together, alternative splicing increases the complexity of the S. commune proteome considerably and provides it with a rich repertoire of alternative functionality that is exploited dynamically.

Living apart together: crosstalk between the core and supernumerary genomes in a fungal plant pathogen

BACKGROUND

Eukaryotes display remarkable genome plasticity, which can include supernumerary chromosomes that differ markedly from the core chromosomes. Despite the widespread occurrence of supernumerary chromosomes in fungi, their origin, relation to the core genome and the reason for their divergent characteristics are still largely unknown. The complexity of genome assembly due to the presence of repetitive DNA partially accounts for this.

RESULTS

Here we use single-molecule real-time (SMRT) sequencing to assemble the genome of a prominent fungal wheat pathogen, Fusarium poae, including at least one supernumerary chromosome. The core genome contains limited transposable elements (TEs) and no gene duplications, while the supernumerary genome holds up to 25 % TEs and multiple gene duplications. The core genome shows all hallmarks of repeat-induced point mutation (RIP), a defense mechanism against TEs, specific for fungi. The absence of RIP on the supernumerary genome accounts for the differences between the two (sub)genomes, and results in a functional crosstalk between them. The supernumerary genome is a reservoir for TEs that migrate to the core genome, and even large blocks of supernumerary sequence (>200 kb) have recently translocated to the core. Vice versa, the supernumerary genome acts as a refuge for genes that are duplicated from the core genome.

CONCLUSIONS

For the first time, a mechanism was determined that explains the differences that exist between the core and supernumerary genome in fungi. Different biology rather than origin was shown to be responsible. A "living apart together" crosstalk exists between the core and supernumerary genome, accelerating chromosomal and organismal evolution.

Genome-Wide Association Study Identifies Novel Candidate Genes for Aggressiveness, Deoxynivalenol Production, and Azole Sensitivity in Natural Field Populations of Fusarium graminearum

Genome-wide association studies can identify novel genomic regions and genes that affect quantitative traits. Fusarium head blight is a destructive disease caused by Fusarium graminearum that exhibits several quantitative traits, including aggressiveness, mycotoxin production, and fungicide resistance. Restriction site-associated DNA sequencing was performed for 220 isolates of F. graminearum. A total of 119 isolates were phenotyped for aggressiveness and deoxynivalenol (DON) production under natural field conditions across four environments. The effective concentration of propiconazole that inhibits isolate growth in vitro by 50% was calculated for 220 strains. Approximately 29,000 single nucleotide polymorphism markers were associated to each trait, resulting in 50, 29, and 74 quantitative trait nucleotides (QTNs) that were significantly associated to aggressiveness, DON production, and propiconazole sensitivity, respectively. Approximately 41% of these QTNs caused nonsynonymous substitutions in predicted exons, while the remainder were synonymous substitutions or located in intergenic regions. Three QTNs associated with propiconazole sensitivity were significant after Bonferroni correction. These QTNs were located in genes not previously associated with azole sensitivity. The majority of the detected QTNs were located in genes with predicted regulatory functions, suggesting that nucleotide variation in regulatory genes plays a major role in the corresponding quantitative trait variation.

Alternatively spliced, spliceosomal twin introns in Helminthosporium solani

Spliceosomal twin introns, "stwintrons", have been defined as complex intervening sequences that carry a second intron ("internal intron") interrupting one of the conserved sequence domains necessary for their correct splicing via consecutive excision events. Previously, we have described and experimentally verified stwintrons in species of Sordariomycetes, where an "internal intron" interrupted the donor sequence of an "external intron". Here we describe and experimentally verify two novel stwintrons of the potato pathogen Helminthosporium solani. One instance involves alternative splicing of an internal intron interrupting the donor domain of an external intron and a second one interrupting the acceptor domain of an overlapping external intron, both events leading to identical mature mRNAs. In the second case, an internal intron interrupts the donor domain of the external intron, while an alternatively spliced intron leads to an mRNA carrying a premature chain termination codon. We thus extend the stwintron concept to the acceptor domain and establish a link of the occurrence of stwintrons with that of alternative splicing.

Transcriptome analysis of the white pine blister rust pathogen Cronartium ribicola: de novo assembly, expression profiling, and identification of candidate effectors

Background

The fungus Cronartium ribicola (Cri) is an economically and ecologically important forest pathogen that causes white pine blister rust (WPBR) disease on five-needle pines. To cause stem cankers and kill white pine trees the fungus elaborates a life cycle with five stages of spore development on five-needle pines and the alternate host Ribes plants. To increase our understanding of molecular WP-BR interactions, here we report genome-wide transcriptional profile analysis of C. ribicola using RNA-seq.

Results

cDNA libraries were constructed from aeciospore, urediniospore, and western white pine (Pinus monticola) tissues post Cri infection. Over 200 million RNA-seq 100-bp paired-end (PE) reads from rust fungal spores were de novo assembled and a reference transcriptome was generated with 17,880 transcripts that were expressed from 13,629 unigenes. A total of 734 unique proteins were predicted as a part of the Cri secretome from complete open reading frames (ORFs), and 41 % of them were Cronartium-specific. This study further identified a repertoire of candidate effectors and other pathogenicity determinants. Differentially expressed genes (DEGs) were identified to gain an understanding of molecular events important during the WPBR fungus life cycle by comparing Cri transcriptomes at different infection stages. Large-scale changes of in planta gene expression profiles were observed, revealing that multiple fungal biosynthetic pathways were enhanced during mycelium growth inside infected pine stem tissues. Conversely, many fungal genes that were up-regulated at the urediniospore stage appeared to be signalling components and transporters. The secreted fungal protein genes that were up-regulated in pine needle tissues during early infection were primarily associated with cell wall modifications, possibly to mask the rust pathogen from plant defenses.

Conclusion

This comprehensive transcriptome profiling substantially improves our current understanding of molecular WP-BR interactions. The repertoire of candidate effectors and other putative pathogenicity determinants identified here are valuable for future functional analysis of Cri virulence and pathogenicity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1861-1) contains supplementary material, which is available to authorized users.

RNA-seq Profiles of Immune Related Genes in the Spleen of Necrotic Enteritis-afflicted Chicken Lines

The study aimed to compare the necrotic enteritis (NE)-induced transcriptome differences between the spleens of Marek’s disease resistant chicken line 6.3 and susceptible line 7.2 co-infected with Eimeria maxima/Clostridium perfringens using RNA-Seq. Total RNA from the spleens of two chicken lines were used to make libraries, generating 42,736,296 and 42,617,720 usable reads, which were assembled into groups of 29,897 and 29,833 mRNA genes, respectively. The transcriptome changes were investigated using the differentially expressed genes (DEGs) package, which indicated 3,255, 2,468 and 2,234 DEGs of line 6.3, line 7.2, and comparison between two lines, respectively (fold change ≥2, p<0.01). The transcription levels of 14 genes identified were further examined using qRT-PCR. The results of qRT-PCR were consistent with the RNA-seq data. All of the DEGs were analysed using gene ontology terms, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database and the DEGs in each term were found to be more highly expressed in line 6.3 than in line 7.2. RNA-seq analysis indicated 139 immune related genes, 44 CD molecular genes and 150 cytokines genes which were differentially expressed among chicken lines 6.3 and 7.2 (fold change ≥2, p<0.01). Novel mRNA analysis indicated 15,518 novel genes, for which the expression was shown to be higher in line 6.3 than in line 7.2 including some immune-related targets. These findings will help to understand host-pathogen interaction in the spleen and elucidate the mechanism of host genetic control of NE, and provide basis for future studies that can lead to the development of marker-based selection of highly disease-resistant chickens.

The exosome controls alternative splicing by mediating the gene expression and assembly of the spliceosome complex

The exosome is a complex with exoribonuclease activity that regulates RNA surveillance and turnover. The exosome also plays a role in regulating the degradation of precursor mRNAs to maintain the expression of splicing variants. In Neurospora, the silencing of rrp44, which encodes the catalytic subunit of the exosome, changed the expression of a set of spliceosomal snRNA, snRNP genes and SR protein related genes. The knockdown of rrp44 also affected the assembly of the spliceosome. RNA-seq analysis revealed a global change in bulk splicing events. Exosome-mediated splicing may regulate alternative splicing of NCU05290, NCU07421 and the circadian clock gene frequency (frq). The knockdown of rrp44 led to an increased ratio of splicing variants without intron 6 (I-6) and shorter protein isoform small FRQ (s-FRQ) as a consequence. These findings suggest that the exosome controls splicing events by regulating the degradation of precursor mRNAs and the gene expression, assembly and function of the spliceosome.

CodingQuarry: highly accurate hidden Markov model gene prediction in fungal genomes using RNA-seq transcripts

BACKGROUND

The impact of gene annotation quality on functional and comparative genomics makes gene prediction an important process, particularly in non-model species, including many fungi. Sets of homologous protein sequences are rarely complete with respect to the fungal species of interest and are often small or unreliable, especially when closely related species have not been sequenced or annotated in detail. In these cases, protein homology-based evidence fails to correctly annotate many genes, or significantly improve ab initio predictions. Generalised hidden Markov models (GHMM) have proven to be invaluable tools in gene annotation and, recently, RNA-seq has emerged as a cost-effective means to significantly improve the quality of automated gene annotation. As these methods do not require sets of homologous proteins, improving gene prediction from these resources is of benefit to fungal researchers. While many pipelines now incorporate RNA-seq data in training GHMMs, there has been relatively little investigation into additionally combining RNA-seq data at the point of prediction, and room for improvement in this area motivates this study.

RESULTS

CodingQuarry is a highly accurate, self-training GHMM fungal gene predictor designed to work with assembled, aligned RNA-seq transcripts. RNA-seq data informs annotations both during gene-model training and in prediction. Our approach capitalises on the high quality of fungal transcript assemblies by incorporating predictions made directly from transcript sequences. Correct predictions are made despite transcript assembly problems, including those caused by overlap between the transcripts of adjacent gene loci. Stringent benchmarking against high-confidence annotation subsets showed CodingQuarry predicted 91.3% of Schizosaccharomyces pombe genes and 90.4% of Saccharomyces cerevisiae genes perfectly. These results are 4-5% better than those of AUGUSTUS, the next best performing RNA-seq driven gene predictor tested. Comparisons against whole genome Sc. pombe and S. cerevisiae annotations further substantiate a 4-5% improvement in the number of correctly predicted genes.

CONCLUSIONS

We demonstrate the success of a novel method of incorporating RNA-seq data into GHMM fungal gene prediction. This shows that a high quality annotation can be achieved without relying on protein homology or a training set of genes. CodingQuarry is freely available ( https://sourceforge.net/projects/codingquarry/ ), and suitable for incorporation into genome annotation pipelines.

RNA-Seq-based transcriptome analysis of aflatoxigenic Aspergillus flavus in response to water activity

Aspergillus flavus is one of the most important producers of carcinogenic aflatoxins in crops, and the effect of water activity (a_w) on growth and aflatoxin production of A. flavus has been previously studied. Here we found the strains under 0.93 a_w exhibited decreased conidiation and aflatoxin biosynthesis compared to that under 0.99 a_w. When RNA-Seq was used to delineate gene expression profile under different water activities, 23,320 non-redundant unigenes, with an average length of 1297 bp, were yielded. By database comparisons, 19,838 unigenes were matched well (e-value < 10⁻⁵) with known gene sequences, and another 6767 novel unigenes were obtained by comparison to the current genome annotation of A. flavus. Based on the RPKM equation, 5362 differentially expressed unigenes (with |log₂Ratio| ≥ 1) were identified between 0.99 a_w and 0.93 a_w treatments, including 3156 up-regulated and 2206 down-regulated unigenes, suggesting that A. flavus underwent an extensive transcriptome response during water activity variation. Furthermore, we found that the expression of 16 aflatoxin producing-related genes decreased obviously when water activity decreased, and the expression of 11 development-related genes increased after 0.99 a_w treatment. Our data corroborate a model where water activity affects aflatoxin biosynthesis through increasing the expression of aflatoxin producing-related genes and regulating development-related genes.

New gene models and alternative splicing in the maize pathogen Colletotrichum graminicola revealed by RNA-Seq analysis

BACKGROUND

An annotated genomic sequence of the corn anthracnose fungus Colletotrichum graminicola has been published previously, but correct identification of gene models by means of automated gene annotation remains a challenge. RNA-Seq offers the potential for substantially improved gene annotations and for the identification of posttranscriptional RNA modifications, such as alternative splicing and RNA editing.

RESULTS

Based on the nucleotide sequence information of transcripts, we identified 819 novel transcriptionally active regions (nTARs) and revised 906 incorrectly predicted gene models, including revisions of exon-intron structure, gene orientation and sequencing errors. Among the nTARs, 146 share significant similarity with proteins that have been identified in other species suggesting that they are hitherto unidentified genes in C. graminicola. Moreover, 5'- and 3'-UTR sequences of 4378 genes have been retrieved and alternatively spliced variants of 69 genes have been identified. Comparative analysis of RNA-Seq data and the genome sequence did not provide evidence for RNA editing in C. graminicola.

CONCLUSIONS

We successfully employed deep sequencing RNA-Seq data in combination with an elaborate bioinformatics strategy in order to identify novel genes, incorrect gene models and mechanisms of transcript processing in the corn anthracnose fungus C. graminicola. Sequence data of the revised genome annotation including several hundreds of novel transcripts, improved gene models and candidate genes for alternative splicing have been made accessible in a comprehensive database. Our results significantly contribute to both routine laboratory experiments and large-scale genomics or transcriptomic studies in C. graminicola.