Huiquan Liu

Functional analysis of the kinome of the wheat scab fungus Fusarium graminearum.

As in other eukaryotes, protein kinases play major regulatory roles in filamentous fungi. Although the genomes of many plant pathogenic fungi have been sequenced, systematic characterization of their kinomes has not been reported. The wheat scab fungus Fusarium graminearum has 116 protein kinases (PK) genes. Although twenty of them appeared to be essential, we generated deletion mutants for the other 96 PK genes, including 12 orthologs of essential genes in yeast. All of the PK mutants were assayed for changes in 17 phenotypes, including growth, conidiation, pathogenesis, stress responses, and sexual reproduction. Overall, deletion of 64 PK genes resulted in at least one of the phenotypes examined, including three mutants blocked in conidiation and five mutants with increased tolerance to hyperosmotic stress. In total, 42 PK mutants were significantly reduced in virulence or non-pathogenic, including mutants deleted of key components of the cAMP signaling and three MAPK pathways. A number of these PK genes, including Fg03146 and Fg04770 that are unique to filamentous fungi, are dispensable for hyphal growth and likely encode novel fungal virulence factors. Ascospores play a critical role in the initiation of wheat scab. Twenty-six PK mutants were blocked in perithecia formation or aborted in ascosporogenesis. Additional 19 mutants were defective in ascospore release or morphology. Interestingly, F. graminearum contains two aurora kinase genes with distinct functions, which has not been reported in fungi. In addition, we used the interlog approach to predict the PK-PK and PK-protein interaction networks of F. graminearum. Several predicted interactions were verified with yeast two-hybrid or co-immunoprecipitation assays. To our knowledge, this is the first functional characterization of the kinome in plant pathogenic fungi. Protein kinase genes important for various aspects of growth, developmental, and infection processes in F. graminearum were identified in this study.

High genome heterozygosity and endemic genetic recombination in the wheat stripe rust fungus

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most destructive diseases of wheat. Here we report a 110-Mb draft sequence of Pst isolate CY32, obtained using a ‘fosmid-to-fosmid’ strategy, to better understand its race evolution and pathogenesis. The Pst genome is highly heterozygous and contains 25,288 protein-coding genes. Compared with non-obligate fungal pathogens, Pst has a more diverse gene composition and more genes encoding secreted proteins. Re-sequencing analysis indicates significant genetic variation among six isolates collected from different continents. Approximately 35% of SNPs are in the coding sequence regions, and half of them are non-synonymous. High genetic diversity in Pst suggests that sexual reproduction has an important role in the origin of different regional races. Our results show the effectiveness of the ‘fosmid-to-fosmid’ strategy for sequencing dikaryotic genomes and the feasibility of genome analysis to understand race evolution in Pst and other obligate pathogens.

Genome-wide A-to-I RNA editing in fungi independent of ADAR enzymes.

Yeasts and filamentous fungi do not have adenosine deaminase acting on RNA (ADAR) orthologs and are believed to lack A-to-I RNA editing, which is the most prevalent editing of mRNA in animals. However, during this study with the PUK1(FGRRES_01058) pseudokinase gene important for sexual reproduction in Fusarium graminearum, we found that two tandem stop codons, UA(1831)GUA(1834)G, in its kinase domain were changed to UG(1831)GUG(1834)G by RNA editing in perithecia. To confirm A-to-I editing of PUK1 transcripts, strand-specific RNA-seq data were generated with RNA isolated from conidia, hyphae, and perithecia. PUK1 was almost specifically expressed in perithecia, and 90% of transcripts were edited to UG(1831)GUG(1834)G. Genome-wide analysis identified 26,056 perithecium-specific A-to-I editing sites. Unlike those in animals, 70.5% of A-to-I editing sites inF. graminearum occur in coding regions, and more than two-thirds of them result in amino acid changes, including editing of 69PUK1-like pseudogenes with stop codons in ORFs.PUK1orthologs and other pseudogenes also displayed stage-specific expression and editing in Neurospora crassa and F. verticillioides Furthermore,F. graminearum differs from animals in the sequence preference and structure selectivity of A-to-I editing sites. Whereas As embedded in RNA stems are targeted by ADARs, RNA editing inF. graminearum preferentially targets As in hairpin loops, which is similar to the anticodon loop of tRNA targeted by adenosine deaminases acting on tRNA (ADATs). Overall, our results showed that A-to-I RNA editing occurs specifically during sexual reproduction and mainly in the coding regions in filamentous ascomycetes, involving adenosine deamination mechanisms distinct from metazoan ADARs.

Genome sequence of Valsa canker pathogens uncovers a potential adaptation of colonization of woody bark

Canker caused by ascomycetous Valsa species are among the most destructive diseases of woody plants worldwide. These pathogens are distinct from other pathogens because they only effectively attack tree bark in the field. To unravel the potential adaptation mechanism of bark colonization, we examined the genomes of Valsa mali and Valsa pyri that preferentially infect apple and pear, respectively. We reported the 44.7 and 35.7 Mb genomes of V. mali and V. pyri, respectively. We also identified the potential genomic determinants of wood colonization by comparing them with related cereal pathogens. Both genomes encode a plethora of pathogenicity-related genes involved in plant cell wall degradation and secondary metabolite biosynthesis. In order to adapt to the nutrient limitation and low pH environment in bark, they seem to employ membrane transporters associated with nitrogen uptake and secrete proteases predominantly with acidic pH optima. Remarkably, both Valsa genomes are especially suited for pectin decomposition, but are limited in lignocellulose and cutin degradation. Besides many similarities, the two genomes show distinct variations in many secondary metabolism gene clusters. Our results show a potential adaptation of Valsa canker pathogens to colonize woody bark. Secondary metabolism gene clusters are probably responsible for this host specificity.

Two Cdc2 Kinase Genes with Distinct Functions in Vegetative and Infectious Hyphae in Fusarium graminearum

Eukaryotic cell cycle involves a number of protein kinases important for the onset and progression through mitosis, most of which are well characterized in the budding and fission yeasts and conserved in other fungi. However, unlike the model yeast and filamentous fungi that have a single Cdc2 essential for cell cycle progression, the wheat scab fungus Fusarium graminearum contains two CDC2 orthologs. The cdc2A and cdc2B mutants had no obvious defects in growth rate and conidiation but deletion of both of them is lethal, indicating that these two CDC2 orthologs have redundant functions during vegetative growth and asexual reproduction. However, whereas the cdc2B mutant was normal, the cdc2A mutant was significantly reduced in virulence and rarely produced ascospores. Although deletion of CDC2A had no obvious effect on the formation of penetration branches or hyphopodia, the cdc2A mutant was limited in the differentiation and growth of infectious growth in wheat tissues. Therefore, CDC2A plays stage-specific roles in cell cycle regulation during infectious growth and sexual reproduction. Both CDC2A and CDC2B are constitutively expressed but only CDC2A was up-regulated during plant infection and ascosporogenesis. Localization of Cdc2A- GFP to the nucleus but not Cdc2B-GFP was observed in vegetative hyphae, ascospores, and infectious hyphae. Complementation assays with chimeric fusion constructs showed that both the N- and C-terminal regions of Cdc2A are important for its functions in pathogenesis and ascosporogenesis but only the N-terminal region is important for its subcellular localization. Among the Sordariomycetes, only three Fusarium species closely related to F. graminearum have two CDC2 genes. Furthermore, F. graminearum uniquely has two Aurora kinase genes and one additional putative cyclin gene, and its orthologs of CAK1 and other four essential mitotic kinases in the budding yeast are dispensable for viability. Overall, our data indicate that cell cycle regulation is different between vegetative and infectious hyphae in F. graminearum and Cdc2A, possibly by interacting with a stage-specific cyclin, plays a more important role than Cdc2B during ascosporogenesis and plant infection.

Molecular evolution and functional divergence of tubulin superfamily in the fungal tree of life

Microtubules are essential for various cellular activities and β-tubulins are the target of benzimidazole fungicides. However, the evolution and molecular mechanisms driving functional diversification in fungal tubulins are not clear. In this study, we systematically identified tubulin genes from 59 representative fungi across the fungal kingdom. Phylogenetic analysis showed that α-/β-tubulin genes underwent multiple independent duplications and losses in different fungal lineages and formed distinct paralogous/orthologous clades. The last common ancestor of basidiomycetes and ascomycetes likely possessed two paralogs of α-tubulin (α1/α2) and β-tubulin (β1/β2) genes but α2-tubulin genes were lost in basidiomycetes and β2-tubulin genes were lost in most ascomycetes. Molecular evolutionary analysis indicated that α1, α2, and β2-tubulins have been under strong divergent selection and adaptive positive selection. Many positively selected sites are at or adjacent to important functional sites and likely contribute to functional diversification. We further experimentally confirmed functional divergence of two β-tubulins in Fusarium and identified type II variations in FgTub2 responsible for function shifts. In this study, we also identified δ-/ε-/η-tubulins in Chytridiomycetes. Overall, our results illustrated that different evolutionary mechanisms drive functional diversification of α-/β-tubulin genes in different fungal lineages, and residues under positive selection could provide targets for further experimental study.

PKA activity is essential for relieving the suppression of hyphal growth and appressorium formation by MoSfl1 in Magnaporthe oryzae

In the rice blast fungus Magnaporthe oryzae, the cAMP-PKA pathway regulates surface recognition, appressorium turgor generation, and invasive growth. However, deletion of CPKA failed to block appressorium formation and responses to exogenous cAMP. In this study, we generated and characterized the cpk2 and cpkA cpk2 mutants and spontaneous suppressors of cpkA cpk2 in M. oryzae. Our results demonstrate that CPKA and CPK2 have specific and overlapping functions, and PKA activity is essential for appressorium formation and plant infection. Unlike the single mutants, the cpkA cpk2 mutant was significantly reduced in growth and rarely produced conidia. It failed to form appressoria although the intracellular cAMP level and phosphorylation of Pmk1 MAP kinase were increased. The double mutant also was defective in plant penetration and Mps1 activation. Interestingly, it often produced fast-growing spontaneous suppressors that formed appressoria but were still non-pathogenic. Two suppressor strains of cpkA cpk2 had deletion and insertion mutations in the MoSFL1 transcription factor gene. Deletion of MoSFL1 or its C-terminal 93-aa (MoSFL1ΔCT) was confirmed to suppress the defects of cpkA cpk2 in hyphal growth but not appressorium formation or pathogenesis. We also isolated 30 spontaneous suppressors of the cpkA cpk2 mutant in Fusarium graminearum and identified mutations in 29 of them in FgSFL1. Affinity purification and co-IP assays showed that this C-terminal region of MoSfl1 was essential for its interaction with the conserved Cyc8-Tup1 transcriptional co-repressor, which was reduced by cAMP treatment. Furthermore, the S211D mutation at the conserved PKA-phosphorylation site in MoSFL1 partially suppressed the defects of cpkA cpk2. Overall, our results indicate that PKA activity is essential for appressorium formation and proper activation of Pmk1 or Mps1 in M. oryzae, and phosphorylation of MoSfl1 by PKA relieves its interaction with the Cyc8-Tup1 co-repressor and suppression of genes important for hyphal growth.

A-to-I RNA editing independent of ADARs in filamentous fungi

ABSTRACT ADAR mediated A-to-I RNA editing is thought to be unique to animals and occurs mainly in the non-coding regions. Recently filamentous fungi such as Fusarium graminearum were found to lack orthologs of animal ADARs but have stage-specific A-to-I editing during sexual reproduction. Unlike animals, majority of editing sites are in the coding regions and often result in missense and stop loss changes in fungi. Furthermore, whereas As in RNA stems are targeted by animal ADARs, RNA editing in fungi preferentially targets As in hairpin loops, implying that fungal RNA editing involves mechanisms related to editing of the anticodon loop by ADATs. Identification and characterization of fungal adenosine deaminases and their stage-specific co-factors may be helpful to understand the evolution of human ADARs. Fungi also can be used to study biological functions of missense and stop loss RNA editing events in eukaryotic organisms.

Identification of a fungi-specific lineage of protein kinases closely related to tyrosine kinases.

Tyrosine kinases (TKs) specifically catalyze the phosphorylation of tyrosine residues in proteins and play essential roles in many cellular processes. Although TKs mainly exist in animals, recent studies revealed that some organisms outside the Opisthokont clade also contain TKs. The fungi, as the sister group to animals, are thought to lack TKs. To better understand the origin and evolution of TKs, it is important to investigate if fungi have TK or TK-related genes. We therefore systematically identified possible TKs across the fungal kingdom by using the profile hidden Markov Models searches and phylogenetic analyses. Our results confirmed that fungi lack the orthologs of animal TKs. We identified a fungi-specific lineage of protein kinases (FslK) that appears to be a sister group closely related to TKs. Sequence analysis revealed that members of the FslK clade contain all the conserved protein kinase sub-domains and thus are likely enzymatically active. However, they lack key amino acid residues that determine TK-specific activities, indicating that they are not true TKs. Phylogenetic analysis indicated that the last common ancestor of fungi may have possessed numerous members of FslK. The ancestral FslK genes were lost in Ascomycota and Ustilaginomycotina and Pucciniomycotina of Basidiomycota during evolution. Most of these ancestral genes, however, were retained and expanded in Agaricomycetes. The discovery of the fungi-specific lineage of protein kinases closely related to TKs helps shed light on the origin and evolution of TKs and also has potential implications for the importance of these kinases in mushroom fungi.

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

Significance This study systematically identified adenosine to inosine (A-to-I) editing sites in Neurospora crassa and showed the existence of stage-specific editing events at different sexual stages. Unlike in humans, fungal A-to-I editing mainly occurred in coding regions and caused nonsynonymous changes that significantly increased proteome complexity. In general, nonsynonymous editing sites in Neurospora are adaptive and favored by positive selection. RNA editing enables stage-specific functions or expression of proteins important for different sexual developmental processes. Some editing events are well conserved and may affect genes important for other genetic and epigenetic phenomena occurring during sexual reproduction. Overall, our results provide insights into the complex regulation of sexual development and reveal the role of A-to-I editing for adaptive evolution in Neurospora. Although fungi lack adenosine deaminase acting on RNA (ADAR) enzymes, adenosine to inosine (A-to-I) RNA editing was reported recently in Fusarium graminearum during sexual reproduction. In this study, we profiled the A-to-I editing landscape and characterized its functional and adaptive properties in the model filamentous fungus Neurospora crassa. A total of 40,677 A-to-I editing sites were identified, and approximately half of them displayed stage-specific editing or editing levels at different sexual stages. RNA-sequencing analysis with the Δstc-1 and Δsad-1 mutants confirmed A-to-I editing occurred before ascus development but became more prevalent during ascosporogenesis. Besides fungal-specific sequence and secondary structure preference, 63.5% of A-to-I editing sites were in the coding regions and 81.3% of them resulted in nonsynonymous recoding, resulting in a significant increase in the proteome complexity. Many genes involved in RNA silencing, DNA methylation, and histone modifications had extensive recoding, including sad-1, sms-3, qde-1, and dim-2. Fifty pseudogenes harbor premature stop codons that require A-to-I editing to encode full-length proteins. Unlike in humans, nonsynonymous editing events in N. crassa are generally beneficial and favored by positive selection. Almost half of the nonsynonymous editing sites in N. crassa are conserved and edited in Neurospora tetrasperma. Furthermore, hundreds of them are conserved in F. graminearum and had higher editing levels. Two unknown genes with editing sites conserved between Neurospora and Fusarium were experimentally shown to be important for ascosporogenesis. This study comprehensively analyzed A-to-I editing in N. crassa and showed that RNA editing is stage-specific and generally adaptive, and may be functionally related to repeat induced point mutation and meiotic silencing by unpaired DNA.