Genomic surveillance uncovers a pandemic clonal lineage of the wheat blast fungus

Wheat, one of the most important food crops, is threatened by a blast disease pandemic. Here, we show that a clonal lineage of the wheat blast fungus recently spread to Asia and Africa following two independent introductions from South America. Through a combination of genome analyses and laboratory experiments, we show that the decade-old blast pandemic lineage can be controlled by the Rmg8 disease resistance gene and is sensitive to strobilurin fungicides. However, we also highlight the potential of the pandemic clone to evolve fungicide-insensitive variants and sexually recombine with African lineages. This underscores the urgent need for genomic surveillance to track and mitigate the spread of wheat blast outside of South America and to guide preemptive wheat breeding for blast resistance.

Genomic surveillance uncovers a pandemic clonal lineage of the wheat blast fungus

Wheat, one of the most important food crops, is threatened by a blast disease pandemic. Here, we show that a clonal lineage of the wheat blast fungus recently spread to Asia and Africa following two independent introductions from South America. Through a combination of genome analyses and laboratory experiments, we show that the decade-old blast pandemic lineage can be controlled by the Rmg8 disease resistance gene and is sensitive to strobilurin fungicides. However, we also highlight the potential of the pandemic clone to evolve fungicide-insensitive variants and sexually recombine with African lineages. This underscores the urgent need for genomic surveillance to track and mitigate the spread of wheat blast outside of South America and to guide preemptive wheat breeding for blast resistance.

Genome-wide characterization of mitochondrial DNA methylation in human brain

Background

There is growing interest in the role of DNA methylation in regulating the transcription of mitochondrial genes, particularly in brain disorders characterized by mitochondrial dysfunction. Here, we present a novel approach to interrogate the mitochondrial DNA methylome at single base resolution using targeted bisulfite sequencing. We applied this method to investigate mitochondrial DNA methylation patterns in post-mortem superior temporal gyrus and cerebellum brain tissue from seven human donors.

Results

We show that mitochondrial DNA methylation patterns are relatively low but conserved, with peaks in DNA methylation at several sites, such as within the D-LOOP and the genes MT-ND2, MT-ATP6, MT-ND4, MT-ND5 and MT-ND6, predominantly in a non-CpG context. The elevated DNA methylation we observe in the D-LOOP we validate using pyrosequencing. We identify loci that show differential DNA methylation patterns associated with age, sex and brain region. Finally, we replicate previously reported differentially methylated regions between brain regions from a methylated DNA immunoprecipitation sequencing study.

Conclusions

We have annotated patterns of DNA methylation at single base resolution across the mitochondrial genome in human brain samples. Looking to the future this approach could be utilized to investigate the role of mitochondrial epigenetic mechanisms in disorders that display mitochondrial dysfunction.

Appressorium-mediated plant infection by Magnaporthe oryzae is regulated by a Pmk1-dependent hierarchical transcriptional network

Rice blast is a devastating disease caused by the fungal pathogen Magnaporthe oryzae that threatens rice production around the world. The fungus produces a specialized infection cell, called the appressorium, that enables penetration through the plant cell wall in response to surface signals from the rice leaf. The underlying biology of plant infection, including the regulation of appressorium formation, is not completely understood. Here we report the identification of a network of temporally coregulated transcription factors that act downstream of the Pmk1 mitogen-activated protein kinase pathway to regulate gene expression during appressorium-mediated plant infection. We show that this tiered regulatory mechanism involves Pmk1-dependent phosphorylation of the Hox7 homeobox transcription factor, which regulates genes associated with induction of major physiological changes required for appressorium development-including cell-cycle control, autophagic cell death, turgor generation and melanin biosynthesis-as well as controlling a additional set of virulence-associated transcription factor-encoding genes. Pmk1-dependent phosphorylation of Mst12 then regulates gene functions involved in septin-dependent cytoskeletal re-organization, polarized exocytosis and effector gene expression, which are necessary for plant tissue invasion. Identification of this regulatory cascade provides new potential targets for disease intervention.

Cautionary Notes on Use of the MoT3 Diagnostic Assay for Magnaporthe oryzae Wheat and Rice Blast Isolates

The blast fungus Magnaporthe oryzae is comprised of lineages that exhibit varying degrees of specificity on about 50 grass hosts, including rice, wheat, and barley. Reliable diagnostic tools are essential given that the pathogen has a propensity to jump to new hosts and spread to new geographic regions. Of particular concern is wheat blast, which has suddenly appeared in Bangladesh in 2016 before spreading to neighboring India. In these Asian countries, wheat blast strains are now co-occurring with the destructive rice blast pathogen raising the possibility of genetic exchange between these destructive pathogens. We assessed the recently described MoT3 diagnostic assay and found that it did not distinguish between wheat and rice blast isolates from Bangladesh. The assay is based on primers matching the WB12 sequence corresponding to a fragment of the M. oryzae MGG_02337 gene annotated as a short chain dehydrogenase. These primers could not reliably distinguish between wheat and rice blast isolates from Bangladesh based on DNA amplification experiments performed in separate laboratories in Bangladesh and in the United Kingdom. Specifically, all eight rice blast isolates tested in this study produced the WB12 amplicon. In addition, comparative genomics of the WB12 nucleotide sequence revealed a complex underlying genetic structure with related sequences across M. oryzae strains and in both rice and wheat blast isolates. We, therefore, caution against the indiscriminate use of this assay to identify wheat blast and encourage further development of the assay to ensure its value in diagnosis.

Extreme genome diversity in the hyper-prevalent parasitic eukaryote Blastocystis

Blastocystis is the most prevalent eukaryotic microbe colonizing the human gut, infecting approximately 1 billion individuals worldwide. Although Blastocystis has been linked to intestinal disorders, its pathogenicity remains controversial because most carriers are asymptomatic. Here, the genome sequence of Blastocystis subtype (ST) 1 is presented and compared to previously published sequences for ST4 and ST7. Despite a conserved core of genes, there is unexpected diversity between these STs in terms of their genome sizes, guanine-cytosine (GC) content, intron numbers, and gene content. ST1 has 6,544 protein-coding genes, which is several hundred more than reported for ST4 and ST7. The percentage of proteins unique to each ST ranges from 6.2% to 20.5%, greatly exceeding the differences observed within parasite genera. Orthologous proteins also display extreme divergence in amino acid sequence identity between STs (i.e., 59%-61% median identity), on par with observations of the most distantly related species pairs of parasite genera. The STs also display substantial variation in gene family distributions and sizes, especially for protein kinase and protease gene families, which could reflect differences in virulence. It remains to be seen to what extent these inter-ST differences persist at the intra-ST level. A full 26% of genes in ST1 have stop codons that are created on the mRNA level by a novel polyadenylation mechanism found only in Blastocystis. Reconstructions of pathways and organellar systems revealed that ST1 has a relatively complete membrane-trafficking system and a near-complete meiotic toolkit, possibly indicating a sexual cycle. Unlike some intestinal protistan parasites, Blastocystis ST1 has near-complete de novo pyrimidine, purine, and thiamine biosynthesis pathways and is unique amongst studied stramenopiles in being able to metabolize α-glucans rather than β-glucans. It lacks all genes encoding heme-containing cytochrome P450 proteins. Predictions of the mitochondrion-related organelle (MRO) proteome reveal an expanded repertoire of functions, including lipid, cofactor, and vitamin biosynthesis, as well as proteins that may be involved in regulating mitochondrial morphology and MRO/endoplasmic reticulum (ER) interactions. In sharp contrast, genes for peroxisome-associated functions are absent, suggesting Blastocystis STs lack this organelle. Overall, this study provides an important window into the biology of Blastocystis, showcasing significant differences between STs that can guide future experimental investigations into differences in their virulence and clarifying the roles of these organisms in gut health and disease.

Tpc1 is an important Zn(II)2Cys6 transcriptional regulator required for polarized growth and virulence in the rice blast fungus

The establishment of polarity is a critical process in pathogenic fungi, mediating infection-related morphogenesis and host tissue invasion. Here, we report the identification of TPC1 (Transcription factor for Polarity Control 1), which regulates invasive polarized growth in the rice blast fungus Magnaporthe oryzae. TPC1 encodes a putative transcription factor of the fungal Zn(II)₂Cys₆ family, exclusive to filamentous fungi. Tpc1-deficient mutants show severe defects in conidiogenesis, infection-associated autophagy, glycogen and lipid metabolism, and plant tissue colonisation. By tracking actin-binding proteins, septin-5 and autophagosome components, we show that Tpc1 regulates cytoskeletal dynamics and infection-associated autophagy during appressorium-mediated plant penetration. We found that Tpc1 interacts with Mst12 and modulates its DNA-binding activity, while Tpc1 nuclear localisation also depends on the MAP kinase Pmk1, consistent with the involvement of Tpc1 in this signalling pathway, which is critical for appressorium development. Importantly, Tpc1 directly regulates NOXD expression, the p22^phox subunit of the fungal NADPH oxidase complex via an interaction with Mst12. Tpc1 therefore controls spatial and temporal regulation of cortical F-actin through regulation of the NADPH oxidase complex during appressorium re-polarisation. Consequently, Tpc1 is a core developmental regulator in filamentous fungi, linking the regulated synthesis of reactive oxygen species and the Pmk1 pathway, with polarity control during host invasion.

Cellular polarity is an intrinsic feature of filamentous fungal growth and pathogenesis. In this study, we identified a gene required for fungal polar growth and virulence in the rice blast fungus Magnaporthe oryzae. This gene has been named TPC1 (Transcription factor for Polarity Control 1). The Tpc1 protein belongs to the fungal Zn(II)₂Cys₆ binuclear cluster family. This DNA-binding motif is present exclusively in the fungal kingdom. We have characterised defects associated with lack of Tpc1 in M. oryzae. We show that Tpc1 is involved in polarised growth and virulence. The M. oryzae Δtpc1 mutant shows a delay in glycogen and lipid metabolism, and infection-associated autophagy–processes that regulate appressorium-mediated M. oryzae plant infection. The saprophytic fungus Neurospora crassa contains a Tpc1 homolog (NcTpc1) involved in vegetative growth and sustained tip elongation, suggesting that Tpc1-like proteins act as core regulators of polarised growth and development in filamentous fungi. A comparative transcriptome analysis has allowed us to identify genes regulated by Tpc1 in M. oryzae including NoxD, an important component of the fungal NADPH complex. Significantly, Tpc1 interacts with Mst12, a component of the Pmk1 signalling pathway essential for appressorium development, and modulates Mst12 binding affinity to NOXD promoter region. We conclude that Tpc1 is a key regulator of polarity in M. oryzae that regulates growth, autophagy and septin-mediated reorientation of the F-actin cytoskeleton to facilitate plant cell invasion.

The β-1,3-glucanosyltransferases (Gels) affect the structure of the rice blast fungal cell wall during appressorium-mediated plant infection

The fungal wall is pivotal for cell shape and function, and in interfacial protection during host infection and environmental challenge. Here, we provide the first description of the carbohydrate composition and structure of the cell wall of the rice blast fungus Magnaporthe oryzae. We focus on the family of glucan elongation proteins (Gels) and characterize five putative β‐1,3‐glucan glucanosyltransferases that each carry the Glycoside Hydrolase 72 signature. We generated targeted deletion mutants of all Gel isoforms, that is, the GH72⁺, which carry a putative carbohydrate‐binding module, and the GH72⁻ Gels, without this motif. We reveal that M. oryzaeGH72⁺GELs are expressed in spores and during both infective and vegetative growth, but each individual Gel enzymes are dispensable for pathogenicity. Further, we demonstrated that a Δgel1Δgel3Δgel4 null mutant has a modified cell wall in which 1,3‐glucans have a higher degree of polymerization and are less branched than the wild‐type strain. The mutant showed significant differences in global patterns of gene expression, a hyper‐branching phenotype and no sporulation, and thus was unable to cause rice blast lesions (except via wounded tissues). We conclude that Gel proteins play significant roles in structural modification of the fungal cell wall during appressorium‐mediated plant infection.

Emergence of wheat blast in Bangladesh was caused by a South American lineage of Magnaporthe oryzae

BACKGROUND

In February 2016, a new fungal disease was spotted in wheat fields across eight districts in Bangladesh. The epidemic spread to an estimated 15,000 hectares, about 16 % of the cultivated wheat area in Bangladesh, with yield losses reaching up to 100 %. Within weeks of the onset of the epidemic, we performed transcriptome sequencing of symptomatic leaf samples collected directly from Bangladeshi fields.

RESULTS

Reinoculation of seedlings with strains isolated from infected wheat grains showed wheat blast symptoms on leaves of wheat but not rice. Our phylogenomic and population genomic analyses revealed that the wheat blast outbreak in Bangladesh was most likely caused by a wheat-infecting South American lineage of the blast fungus Magnaporthe oryzae.

CONCLUSION

Our findings suggest that genomic surveillance can be rapidly applied to monitor plant disease outbreaks and provide valuable information regarding the identity and origin of the infectious agent.

Protein kinase C is essential for viability of the rice blast fungus Magnaporthe oryzae

Protein kinase C constitutes a family of serine–threonine kinases found in all eukaryotes and implicated in a wide range of cellular functions, including regulation of cell growth, cellular differentiation and immunity. Here, we present three independent lines of evidence which indicate that protein kinase C is essential for viability of Magnaporthe oryzae. First, all attempts to generate a target deletion of PKC1, the single copy protein kinase C‐encoding gene, proved unsuccessful. Secondly, conditional gene silencing of PKC1 by RNA interference led to severely reduced growth of the fungus, which was reversed by targeted deletion of the Dicer2‐encoding gene, MDL2. Finally, selective kinase inhibition of protein kinase C by targeted allelic replacement with an analogue‐sensitive PKC1^AS allele led to specific loss of fungal viability in the presence of the PP1 inhibitor. Global transcriptional profiling following selective PKC inhibition identified significant changes in gene expression associated with cell wall re‐modelling, autophagy, signal transduction and secondary metabolism. When considered together, these results suggest protein kinase C is essential for growth and development of M. oryzae with extensive downstream targets in addition to the cell integrity pathway. Targeting protein kinase C signalling may therefore prove an effective means of controlling rice blast disease.

FAR1 and FAR2 regulate the expression of genes associated with lipid metabolism in the rice blast fungus Magnaporthe oryzae

The rice blast fungus Magnaporthe oryzae causes plant disease via specialised infection structures called appressoria. These dome-shaped cells are able to generate enormous internal pressure, which enables penetration of rice tissue by invasive hyphae. Previous studies have shown that mobilisation of lipid bodies and subsequent lipid metabolism are essential pre-requisites for successful appressorium-mediated plant infection, which requires autophagic recycling of the contents of germinated spores and germ tubes to the developing appressorium. Here, we set out to identify putative regulators of lipid metabolism in the rice blast fungus. We report the identification of FAR1 and FAR2, which encode highly conserved members of the Zn2-Cys6 family of transcriptional regulators. We generated Δfar1, Δfar2 and Δfar1Δfar2 double mutants in M. oryzae and show that these deletion mutants are deficient in growth on long chain fatty acids. In addition, Δfar2 mutants are also unable to grow on acetate and short chain fatty acids. FAR1 and FAR2 are necessary for differential expression of genes involved in fatty acid β-oxidation, acetyl-CoA translocation, peroxisomal biogenesis, and the glyoxylate cycle in response to the presence of lipids. Furthermore, FAR2 is necessary for expression of genes associated with acetyl-CoA synthesis. Interestingly, Δfar1, Δfar2 and Δfar1Δfar2 mutants show no observable delay or reduction in lipid body mobilisation during plant infection, suggesting that these transcriptional regulators control lipid substrate utilization by the fungus but not the mobilisation of intracellular lipid reserves during infection-related morphogenesis.

Hook is an adapter that coordinates kinesin-3 and dynein cargo attachment on early endosomes

The Ustilago maydis Hook protein Hok1 is part of an evolutionarily conserved protein complex that regulates bidirectional early endosome trafficking by controlling attachment of both kinesin-3 and dynein.

Bidirectional membrane trafficking along microtubules is mediated by kinesin-1, kinesin-3, and dynein. Several organelle-bound adapters for kinesin-1 and dynein have been reported that orchestrate their opposing activity. However, the coordination of kinesin-3/dynein-mediated transport is not understood. In this paper, we report that a Hook protein, Hok1, is essential for kinesin-3– and dynein-dependent early endosome (EE) motility in the fungus Ustilago maydis. Hok1 binds to EEs via its C-terminal region, where it forms a complex with homologues of human fused toes (FTS) and its interactor FTS- and Hook-interacting protein. A highly conserved N-terminal region is required to bind dynein and kinesin-3 to EEs. To change the direction of EE transport, kinesin-3 is released from organelles, and dynein binds subsequently. A chimaera of human Hook3 and Hok1 rescues the hok1 mutant phenotype, suggesting functional conservation between humans and fungi. We conclude that Hok1 is part of an evolutionarily conserved protein complex that regulates bidirectional EE trafficking by controlling attachment of both kinesin-3 and dynein.

A SUF Fe-S cluster biogenesis system in the mitochondrion-related organelles of the anaerobic protist Pygsuia

BACKGROUND

Many microbial eukaryotes have evolved anaerobic alternatives to mitochondria known as mitochondrion-related organelles (MROs). Yet, only a few of these have been experimentally investigated. Here we report an RNA-seq-based reconstruction of the MRO proteome of Pygsuia biforma, an anaerobic representative of an unexplored deep-branching eukaryotic lineage.

RESULTS

Pygsuia's MRO has a completely novel suite of functions, defying existing "function-based" organelle classifications. Most notable is the replacement of the mitochondrial iron-sulfur cluster machinery by an archaeal sulfur mobilization (SUF) system acquired via lateral gene transfer (LGT). Using immunolocalization in Pygsuia and heterologous expression in yeast, we show that the SUF system does indeed localize to the MRO. The Pygsuia MRO also possesses a unique assemblage of features, including: cardiolipin, phosphonolipid, amino acid, and fatty acid metabolism; a partial Kreb's cycle; a reduced respiratory chain; and a laterally acquired rhodoquinone (RQ) biosynthesis enzyme. The latter observation suggests that RQ is an electron carrier of a fumarate reductase-type complex II in this MRO.

CONCLUSIONS

The unique functional profile of this MRO underscores the tremendous plasticity of mitochondrial function within eukaryotes and showcases the role of LGT in forging metabolic mosaics of ancestral and newly acquired organellar pathways.

Identifying the emerging human pathogen Scedosporium prolificans by using a species-specific monoclonal antibody that binds to the melanin biosynthetic enzyme tetrahydroxynaphthalene reductase

The dematiaceous (melanized) fungus Scedosporium prolificans is an emerging and frequently fatal pathogen of immunocompromised humans and which, along with the closely related fungi Pseudallescheria boydii, Scedosporium apiospermum and S. aurantiacum in the Pseudallescheria-Scedosporium complex, is a contributing aetiology to tsunami lung and central nervous system infections in near-drowning victims who have aspirated water laden with spores. At present, the natural habitat of the fungus is largely unknown, and accurate detection methods are needed to identify environmental reservoirs of infectious propagules. In this study, we report the development of a monoclonal antibody (mAb) (CA4) specific to S. prolificans, which does not cross-react with closely related fungi in the Pseudallescheria-Scedosporium complex or with a wide range of mould and yeast species pathogenic to humans. Using genome sequencing of a soil isolate and targeted gene disruption of the CA4 antigen-encoding gene, we show that mAb CA4 binds to the melanin-biosynthetic enzyme tetrahydroxynaphthalene reductase. Enzyme-deficient mutants produce orange-brown or green-brown spore suspensions compared with the black spore suspension of the wild-type strain. Using mAb CA4 and a mAb (HG12) specific to the related fungi P. boydii, P. apiosperma, S. apiospermum and S. aurantiacum, we demonstrate how the mAbs can be used in combination with a semiselective isolation procedure to track these opportunistic pathogens in environmental samples containing mixed populations of human pathogenic fungi. Specificity of mAb CA4 was confirmed by sequencing of the internally transcribed spacer 1 (ITS1)-5.8S-ITS2 rRNA-encoding regions of fungi isolated from estuarine muds.

Genome-wide transcriptional profiling of appressorium development by the rice blast fungus Magnaporthe oryzae

The rice blast fungus Magnaporthe oryzae is one of the most significant pathogens affecting global food security. To cause rice blast disease the fungus elaborates a specialised infection structure called an appressorium. Here, we report genome wide transcriptional profile analysis of appressorium development using next generation sequencing (NGS). We performed both RNA-Seq and High-Throughput SuperSAGE analysis to compare the utility of these procedures for identifying differential gene expression in M. oryzae. We then analysed global patterns of gene expression during appressorium development. We show evidence for large-scale gene expression changes, highlighting the role of autophagy, lipid metabolism and melanin biosynthesis in appressorium differentiation. We reveal the role of the Pmk1 MAP kinase as a key global regulator of appressorium-associated gene expression. We also provide evidence for differential expression of transporter-encoding gene families and specific high level expression of genes involved in quinate uptake and utilization, consistent with pathogen-mediated perturbation of host metabolism during plant infection. When considered together, these data provide a comprehensive high-resolution analysis of gene expression changes associated with cellular differentiation that will provide a key resource for understanding the biology of rice blast disease.

Saprotrophic competitiveness and biocontrol fitness of a genetically modified strain of the plant-growth-promoting fungus Trichoderma hamatum GD12

Trichoderma species are ubiquitous soil fungi that hold enormous potential for the development of credible alternatives to agrochemicals and synthetic fertilizers in sustainable crop production. In this paper, we show that substantial improvements in plant productivity can be met by genetic modification of a plant-growth-promoting and biocontrol strain of Trichoderma hamatum, but that these improvements are obtained in the absence of disease pressure only. Using a quantitative monoclonal antibody-based ELISA, we show that an N-acetyl-β-d-glucosaminidase-deficient mutant of T. hamatum, generated by insertional mutagenesis of the corresponding gene, has impaired saprotrophic competitiveness during antagonistic interactions with Rhizoctonia solani in soil. Furthermore, its fitness as a biocontrol agent of the pre-emergence damping-off pathogen Sclerotinia sclerotiorum is significantly reduced, and its ability to promote plant growth is constrained by the presence of both pathogens. This work shows that while gains in T. hamatum-mediated plant-growth-promotion can be met through genetic manipulation of a single beneficial trait, such a modification has negative impacts on other aspects of its biology and ecology that contribute to its success as a saprotrophic competitor and antagonist of soil-borne pathogens. The work has important implications for fungal morphogenesis, demonstrating a clear link between hyphal architecture and secretory potential. Furthermore, it highlights the need for a holistic approach to the development of genetically modified Trichoderma strains for use as crop stimulants and biocontrol agents in plant agriculture.

Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea

Sclerotinia sclerotiorum and Botrytis cinerea are closely related necrotrophic plant pathogenic fungi notable for their wide host ranges and environmental persistence. These attributes have made these species models for understanding the complexity of necrotrophic, broad host-range pathogenicity. Despite their similarities, the two species differ in mating behaviour and the ability to produce asexual spores. We have sequenced the genomes of one strain of S. sclerotiorum and two strains of B. cinerea. The comparative analysis of these genomes relative to one another and to other sequenced fungal genomes is provided here. Their 38-39 Mb genomes include 11,860-14,270 predicted genes, which share 83% amino acid identity on average between the two species. We have mapped the S. sclerotiorum assembly to 16 chromosomes and found large-scale co-linearity with the B. cinerea genomes. Seven percent of the S. sclerotiorum genome comprises transposable elements compared to <1% of B. cinerea. The arsenal of genes associated with necrotrophic processes is similar between the species, including genes involved in plant cell wall degradation and oxalic acid production. Analysis of secondary metabolism gene clusters revealed an expansion in number and diversity of B. cinerea-specific secondary metabolites relative to S. sclerotiorum. The potential diversity in secondary metabolism might be involved in adaptation to specific ecological niches. Comparative genome analysis revealed the basis of differing sexual mating compatibility systems between S. sclerotiorum and B. cinerea. The organization of the mating-type loci differs, and their structures provide evidence for the evolution of heterothallism from homothallism. These data shed light on the evolutionary and mechanistic bases of the genetically complex traits of necrotrophic pathogenicity and sexual mating. This resource should facilitate the functional studies designed to better understand what makes these fungi such successful and persistent pathogens of agronomic crops.

Resolving the question of trypanosome monophyly: a comparative genomics approach using whole genome data sets with low taxon sampling

Since the first attempts to classify the evolutionary history of trypanosomes, there have been conflicting reports regarding their true phylogenetic relationships and, in particular, their relationships with other vertebrate trypanosomatids, e.g. Leishmania sp., as well as with the many insect parasitising trypanosomatids. Perhaps the issue that has provided most debate is that concerning the monophyly (or otherwise) of genus Trypanosoma and, even with the advent of molecular methods, the findings of numerous studies have varied significantly depending on the gene sequences analysed, number of taxa included, choice of outgroup and phylogenetic methodology. While of arguably limited applied importance, resolution of the question as to whether or not trypanosomes are monophyletic is critical to accurate evaluation of competing, mutually exclusive evolutionary scenarios for these parasites, namely the 'vertebrate-first' or 'insect-first' hypotheses. Therefore, a new approach, which could overcome previous limitations was needed. At its most simple, the problem can be defined within the framework of a trifurcated tree with three hypothetical positions at which the root can be placed. Using BLASTp and whole-genome gene-by-gene phylogenetic analyses of Trypanosoma brucei, Trypanosoma cruzi, Leishmania major and Naegleria gruberi, we have identified 599 gene markers--putative homologues--that were shared between the genomes of these four taxa. Of these, 75 homologous gene families that demonstrate monophyly of the kinetoplastids were identified. We then used these data sets in combination with an additional outgroup, Euglena gracilis, coupled with large-scale gene concatenation and diverse phylogenetic techniques to investigate the relative branching order of T. brucei, T. cruzi and L. major. Our findings confirm the monophyly of genus Trypanosoma and demonstrate that <1% of the analysed gene markers shared between the genomes of T. brucei, T. cruzi and L. major reject the hypothesis that the trypanosomes form a monophyletic group.

Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism

Powdery mildews are phytopathogens whose growth and reproduction are entirely dependent on living plant cells. The molecular basis of this life-style, obligate biotrophy, remains unknown. We present the genome analysis of barley powdery mildew, Blumeria graminis f.sp. hordei (Blumeria), as well as a comparison with the analysis of two powdery mildews pathogenic on dicotyledonous plants. These genomes display massive retrotransposon proliferation, genome-size expansion, and gene losses. The missing genes encode enzymes of primary and secondary metabolism, carbohydrate-active enzymes, and transporters, probably reflecting their redundancy in an exclusively biotrophic life-style. Among the 248 candidate effectors of pathogenesis identified in the Blumeria genome, very few (less than 10) define a core set conserved in all three mildews, suggesting that most effectors represent species-specific adaptations.

Characterization of MoLDB1 required for vegetative growth, infection-related morphogenesis, and pathogenicity in the rice blast fungus Magnaporthe oryzae

An insertional mutagenesis screen in the rice blast fungus, Magnaporthe oryzae, identified a novel mutant, A2-12-3, which is defective in infection-related morphogenesis and pathogenicity. Analysis of the mutation confirmed an insertion into MoLDB1, which putatively encodes an 806-amino-acid protein with a predicted LIM binding domain. Targeted gene deletion mutants of MoLDB1 were unable to produce asexual or sexual spores and were significantly impaired in vegetative growth and fungal virulence. The Δmoldb1 mutants also showed reduced expression of genes coding hydrophobic proteins (e.g. MPG1 and MHP1), resulting in an easily wettable phenotype in vegetative culture. Moreover, the expression of four genes encoding LIM proteins predicted from the M. oryzae genome was significantly downregulated by deletion of MoLDB1. Analysis of an M. oryzae strain expressing a MoLbd1-green fluorescent protein gene fusion was consistent with the protein being nuclear localized. When considered together, MoLdb1 appears to be involved in regulation of cell wall proteins, including hydrophobins and LIM proteins, and is essential for conidiation, sexual development, appressorium formation, and pathogenicity in M. oryzae.

Comparative genome analysis reveals an absence of leucine-rich repeat pattern-recognition receptor proteins in the kingdom Fungi

Background

In plants and animals innate immunity is the first line of defence against attack by microbial pathogens. Specific molecular features of bacteria and fungi are recognised by pattern recognition receptors that have extracellular domains containing leucine rich repeats. Recognition of microbes by these receptors induces defence responses that protect hosts against potential microbial attack.

Methodology/Principal Findings

A survey of genome sequences from 101 species, representing a broad cross-section of the eukaryotic phylogenetic tree, reveals an absence of leucine rich repeat-domain containing receptors in the fungal kingdom. Uniquely, however, fungi possess adenylate cyclases that contain distinct leucine rich repeat-domains, which have been demonstrated to act as an alternative means of perceiving the presence of bacteria by at least one fungal species. Interestingly, the morphologically similar osmotrophic oomycetes, which are taxonomically distant members of the stramenopiles, possess pattern recognition receptors with similar domain structures to those found in plants.

Conclusions

The absence of pattern recognition receptors suggests that fungi may possess novel classes of pattern-recognition receptor, such as the modified adenylate cyclase, or instead rely on secretion of anti-microbial secondary metabolites for protection from microbial attack. The absence of pattern recognition receptors in fungi, coupled with their abundance in oomycetes, suggests this may be a unique characteristic of the fungal kingdom rather than a consequence of the osmotrophic growth form.

A relational database for the discovery of genes encoding amino acid biosynthetic enzymes in pathogenic fungi

Fungal phytopathogens continue to cause major economic impact, either directly, through crop losses, or due to the costs of fungicide application. Attempts to understand these organisms are hampered by a lack of fungal genome sequence data. A need exists, however, to develop specific bioinformatics tools to collate and analyse the sequence data that currently is available. A web-accessible gene discovery database (http://cogeme.ex.ac.uk/biosynthesis.html) was developed as a demonstration tool for the analysis of metabolic and signal transduction pathways in pathogenic fungi using incomplete gene inventories. Using Bayesian probability to analyse the currently available gene information from pathogenic fungi, we provide evidence that the obligate pathogen Blumeria graminis possesses all amino acid biosynthetic pathways found in free-living fungi, such as Saccharomyces cerevisiae. Phylogenetic analysis was also used to deduce a gene history of succinate-semialdehyde dehydrogenase, an enzyme in the glutamate and lysine biosynthesis pathways. The database provides a tool and methodology to researchers to direct experimentation towards predicting pathway conservation in pathogenic microorganisms.

Phylogenomic analysis demonstrates a pattern of rare and ancient horizontal gene transfer between plants and fungi

Horizontal gene transfer (HGT) describes the transmission of genetic material across species boundaries and is an important evolutionary phenomenon in the ancestry of many microbes. The role of HGT in plant evolutionary history is, however, largely unexplored. Here, we compare the genomes of six plant species with those of 159 prokaryotic and eukaryotic species and identify 1689 genes that show the highest similarity to corresponding genes from fungi. We constructed a phylogeny for all 1689 genes identified and all homolog groups available from the rice (Oryza sativa) genome (3177 gene families) and used these to define 14 candidate plant-fungi HGT events. Comprehensive phylogenetic analyses of these 14 data sets, using methods that account for site rate heterogeneity, demonstrated support for nine HGT events, demonstrating an infrequent pattern of HGT between plants and fungi. Five HGTs were fungi-to-plant transfers and four were plant-to-fungi HGTs. None of the fungal-to-plant HGTs involved angiosperm recipients. These results alter the current view of organismal barriers to HGT, suggesting that phagotrophy, the consumption of a whole cell by another, is not necessarily a prerequisite for HGT between eukaryotes. Putative functional annotation of the HGT candidate genes suggests that two fungi-to-plant transfers have added phenotypes important for life in a soil environment. Our study suggests that genetic exchange between plants and fungi is exceedingly rare, particularly among the angiosperms, but has occurred during their evolutionary history and added important metabolic traits to plant lineages.

Comparative genomics of MAP kinase and calcium-calcineurin signalling components in plant and human pathogenic fungi

Mitogen-activated protein kinase (MAPK) cascades and the calcium-calcineurin pathway control fundamental aspects of fungal growth, development and reproduction. Core elements of these signalling pathways are required for virulence in a wide array of fungal pathogens of plants and mammals. In this review, we have used the available genome databases to explore the structural conservation of three MAPK cascades and the calcium-calcineurin pathway in ten different fungal species, including model organisms, plant pathogens and human pathogens. While most known pathway components from the model yeast Saccharomyces cerevisiae appear to be widely conserved among taxonomically and biologically diverse fungi, some of them were found to be restricted to the Saccharomycotina. The presence of multiple paralogues in certain species such as the zygomycete Rhizopus oryzae and the incorporation of new functional domains that are lacking in S. cerevisiae signalling proteins, most likely reflect functional diversification or adaptation as filamentous fungi have evolved to occupy distinct ecological niches.

Moving targets: rapid evolution of oomycete effectors

Plant pathogenic microbes secrete proteins known as effectors, which enter the cytoplasm of plant cells and suppress host defences. Known effectors in oomycete pathogens possess an RXLR-EER motif in their amino acid sequence that is necessary for transport of the effector into a host plant cell. A large number of putative effectors have now been identified in oomycete genomes, the sequences of which show evidence of diversifying selection at their C terminus. Here, we describe recent progress in characterizing RXLR-EER effectors and discuss why so many of these rapidly evolving proteins are encoded by the genomes of plant pathogenic oomycetes.

Comparative genome analysis of filamentous fungi reveals gene family expansions associated with fungal pathogenesis

Fungi and oomycetes are the causal agents of many of the most serious diseases of plants. Here we report a detailed comparative analysis of the genome sequences of thirty-six species of fungi and oomycetes, including seven plant pathogenic species, that aims to explore the common genetic features associated with plant disease-causing species. The predicted translational products of each genome have been clustered into groups of potential orthologues using Markov Chain Clustering and the data integrated into the e-Fungi object-oriented data warehouse (http://www.e-fungi.org.uk/). Analysis of the species distribution of members of these clusters has identified proteins that are specific to filamentous fungal species and a group of proteins found only in plant pathogens. By comparing the gene inventories of filamentous, ascomycetous phytopathogenic and free-living species of fungi, we have identified a set of gene families that appear to have expanded during the evolution of phytopathogens and may therefore serve important roles in plant disease. We have also characterised the predicted set of secreted proteins encoded by each genome and identified a set of protein families which are significantly over-represented in the secretomes of plant pathogenic fungi, including putative effector proteins that might perturb host cell biology during plant infection. The results demonstrate the potential of comparative genome analysis for exploring the evolution of eukaryotic microbial pathogenesis.

A Methodology for Comparative Functional Genomics

The continuing and rapid increase in the number of fully sequenced genomes is creating new opportunities for comparative studies. However, although many genomic databases store data from multiple organisms, for the most part they provide limited support for comparative genomics. We argue that refocusing genomic data management to provide more direct support for comparative studies enables systematic identification of important relationships between species, thereby increasing the value that can be obtained from sequenced genomes. The principal result of the paper is a methodology, in which comparative analyses are constructed over a foundation based on sequence clusters and evolutionary relationships. This methodology has been applied in a systematic study of the fungi, and we describe how comparative analyses have been implemented as an analysis library over the e-Fungi data warehouse.

Comparative genome analysis across a kingdom of eukaryotic organisms: specialization and diversification in the fungi

The recent proliferation of genome sequencing in diverse fungal species has provided the first opportunity for comparative genome analysis across a eukaryotic kingdom. Here, we report a comparative study of 34 complete fungal genome sequences, representing a broad diversity of Ascomycete, Basidiomycete, and Zygomycete species. We have clustered all predicted protein-encoding gene sequences from these species to provide a means of investigating gene innovations, gene family expansions, protein family diversification, and the conservation of essential gene functions-empirically determined in Saccharomyces cerevisiae-among the fungi. The results are presented with reference to a phylogeny of the 34 fungal species, based on 29 universally conserved protein-encoding gene sequences. We contrast this phylogeny with one based on gene presence and absence and show that, while the two phylogenies are largely in agreement, there are differences in the positioning of some species. We have investigated levels of gene duplication and demonstrate that this varies greatly between fungal species, although there are instances of coduplication in distantly related fungi. We have also investigated the extent of orthology for protein families and demonstrate unexpectedly high levels of diversity among genes involved in lipid metabolism. These analyses have been collated in the e-Fungi data warehouse, providing an online resource for comparative genomic analysis of the fungi.

e-Fungi: a data resource for comparative analysis of fungal genomes

BACKGROUND

The number of sequenced fungal genomes is ever increasing, with about 200 genomes already fully sequenced or in progress. Only a small percentage of those genomes have been comprehensively studied, for example using techniques from functional genomics. Comparative analysis has proven to be a useful strategy for enhancing our understanding of evolutionary biology and of the less well understood genomes. However, the data required for these analyses tends to be distributed in various heterogeneous data sources, making systematic comparative studies a cumbersome task. Furthermore, comparative analyses benefit from close integration of derived data sets that cluster genes or organisms in a way that eases the expression of requests that clarify points of similarity or difference between species.

DESCRIPTION

To support systematic comparative analyses of fungal genomes we have developed the e-Fungi database, which integrates a variety of data for more than 30 fungal genomes. Publicly available genome data, functional annotations, and pathway information has been integrated into a single data repository and complemented with results of comparative analyses, such as MCL and OrthoMCL cluster analysis, and predictions of signaling proteins and the sub-cellular localisation of proteins. To access the data, a library of analysis tasks is available through a web interface. The analysis tasks are motivated by recent comparative genomics studies, and aim to support the study of evolutionary biology as well as community efforts for improving the annotation of genomes. Web services for each query are also available, enabling the tasks to be incorporated into workflows.

CONCLUSION

The e-Fungi database provides fungal biologists with a resource for comparative studies of a large range of fungal genomes. Its analysis library supports the comparative study of genome data, functional annotation, and results of large scale analyses over all the genomes stored in the database. The database is accessible at http://www.e-fungi.org.uk, as is the WSDL for the web services.

Functional analysis of lipid metabolism in Magnaporthe grisea reveals a requirement for peroxisomal fatty acid beta-oxidation during appressorium-mediated plant infection

The rice blast fungus Magnaporthe grisea infects plants by means of specialized infection structures known as appressoria. Turgor generated in the appressorium provides the invasive force that allows the fungus to breach the leaf cuticle with a narrow-penetration hypha gaining entry to the underlying epidermal cell. Appressorium maturation in M. grisea involves mass transfer of lipid bodies to the developing appressorium, coupled to autophagic cell death in the conidium and rapid lipolysis at the onset of appressorial turgor generation. Here, we report identification of the principal components of lipid metabolism in M. grisea based on genome sequence analysis. We show that deletion of any of the eight putative intracellular triacylglycerol lipase-encoding genes from the fungus is insufficient to prevent plant infection, highlighting the complexity and redundancy associated with appressorial lipolysis. In contrast, we demonstrate that a peroxisomally located multifunctional, fatty acid beta-oxidation enzyme is critical to appressorium physiology, and blocking peroxisomal biogenesis prevents plant infection. Taken together, our results indicate that, although triacylglycerol breakdown in the appressorium involves the concerted action of several lipases, fatty acid metabolism and consequent generation of acetyl CoA are necessary for M. grisea to complete its prepenetration phase of development and enter the host plant.

Peroxisomal carnitine acetyl transferase is required for elaboration of penetration hyphae during plant infection by Magnaporthe grisea

Magnaporthe grisea, the causal agent of rice blast disease, invades plant tissue due to the action of specialized infection structures called appressoria, which are used to breach the leaf cuticle and allow development of intracellular, infectious hyphae. In this report we demonstrate that peroxisomal carnitine acetyl transferase (CAT) activity is necessary for appressorium function, and in particular, for the elaboration of primary penetration hyphae. The major CAT activity in M. grisea is encoded by the PTH2 gene, which shows elevated expression in response to acetate and lipid, and is regulated by the cyclic AMP response pathway. Furthermore, a Pth2-GFP fusion protein colocalizes with a peroxisomal marker protein. Targeted deletion of PTH2, generated mutants that were completely non-pathogenic, lacked CAT activity and were unable to utilize a range of lipid substrates. The impairment of appressorium function in Deltapth2 was associated with a delay in lipid reserve mobilization from germ tubes into developing infection cells, and abnormal chitin distribution in infection structures. Addition of glucose to Deltapth2 mutants partially restored the ability to cause rice blast disease and lipid reserve mobilization. Taken together, our findings provide evidence that Pth2 plays a role in the generation of acetyl CoA pools necessary for appressorium function and rapid elaboration of penetration hyphae during host infection.

Comparative genomic analysis of phytopathogenic fungi using expressed sequence tag (EST) collections

SUMMARY We describe the analysis of 57 727 unique expressed sequence tags (ESTs) from 15 species of phytopathogenic and three species of saprophytic fungi. This resource is held within the COGEME phytopathogen EST database (http://cogeme.ex.ac.uk/). Comparative analysis was performed to investigate the differences between pathogenic and free-living fungi based on a substantial collection of expressed gene sequences and available, completed fungal genome sequences. We report that the expressed gene inventories of pathogenic fungi were not significantly more similar to each other than to those of free-living filamentous fungi. As expected, however, filamentous fungi as a group share more sequences in common than with the free-living yeast species Saccharomyces cerevisiae. Interestingly, ESTs of the obligate biotrophic fungus Blumeria graminis f. sp. hordei were more dissimilar to those of all other fungal species assessed, having a lower number of sequences in common with filamentous ascomycetes studied to date and also possessing a larger proportion of unisequences of unknown function. Our analysis of ESTs in the COGEME database enabled identification of a set of functional groups of genes that are more highly represented in the genomes of pathogenic fungi than non-pathogenic species.

A bioinformatic tool for analysis of EST transcript abundance during infection-related development by Magnaporthe grisea

SUMMARY Information regarding the levels of mRNA transcript abundance under different conditions, or in specific tissue types, can be obtained by analysis of the frequency of EST sequences in randomly sequenced cDNA libraries. Here we report a bioinformatics tool, which provides a means of identifying genes that are differentially expressed during pathogenesis-related development by the rice blast fungus Magnaporthe grisea. A total of 31 534 M. grisea ESTs were obtained from dbEST at NCBI, clustered into 8821 unique sequences (unisequences) and manually annotated. Transcript profiles were then calculated for 958 unigenes identified from eight different cDNA libraries. The data were integrated into the Consortium for Functional Genomics of Microbial Eukaryotes (COGEME) database (http://cogeme.ex.ac.uk/) and a web-based front end was designed to allow users to access and interrogate the generated datasets.

The molecular biology of appressorium turgor generation by the rice blast fungus Magnaporthe grisea

The rice blast fungus Magnaporthe grisea develops specialized infection structures known as appressoria, which develop enormous turgor pressure to bring about plant infection. Turgor is generated by accumulation of compatible solutes, including glycerol, which is synthesized in large quantities in the appressorium. Glycogen, trehalose and lipids represent the most abundant storage products in M. grisea conidia. Trehalose and glycogen are rapidly degraded during conidial germination and it is known that trehalose synthesis is required for virulence of the fungus. Lipid bodies are transported to the developing appressoria and degraded at the onset of turgor generation, in a process that is cAMP-dependent. A combined biochemical and genetic approach is being used to dissect the process of turgor generation in the rice blast fungus.

The genome sequence of the rice blast fungus Magnaporthe grisea

Magnaporthe grisea is the most destructive pathogen of rice worldwide and the principal model organism for elucidating the molecular basis of fungal disease of plants. Here, we report the draft sequence of the M. grisea genome. Analysis of the gene set provides an insight into the adaptations required by a fungus to cause disease. The genome encodes a large and diverse set of secreted proteins, including those defined by unusual carbohydrate-binding domains. This fungus also possesses an expanded family of G-protein-coupled receptors, several new virulence-associated genes and large suites of enzymes involved in secondary metabolism. Consistent with a role in fungal pathogenesis, the expression of several of these genes is upregulated during the early stages of infection-related development. The M. grisea genome has been subject to invasion and proliferation of active transposable elements, reflecting the clonal nature of this fungus imposed by widespread rice cultivation.

Regulation of the MPG1 hydrophobin gene in the rice blast fungus Magnaporthe grisea

The hydrophobin-encoding gene MPG1 of the rice blast fungus Magnaporthe grisea is highly expressed during the initial stages of host plant infection and targeted deletion of the gene results in a mutant strain that is reduced in virulence, conidiation, and appressorium formation. The green fluorescent protein-encoding allele sGFP was used as a reporter to investigate regulatory genes that control MPG1 expression. The MAP kinase-encoding gene PMK1 and the wide domain regulators of nitrogen source utilization, NPR1 and NUT1, were required for full expression of MPG1 in response to starvation stress. The CPKA gene, encoding the catalytic subunit of protein kinase A, was required for repression of MPG1 during growth in rich nutrient conditions. During appressorium morphogenesis, high-level MPG1 expression was found to require the CPKA and NPR1 genes. Expression of a destabilized GFP allele indicated that de novo MPG1 expression occurs during appressorium formation. Three regions of the MPG1 promoter were identified which are required for high-level expression of MPG1 during appressorium formation and are necessary for the biological activity of the MPG1 hydrophobin during spore formation and plant infection.

Genomics of phytopathogenic fungi and the development of bioinformatic resources

Genomic resources available to researchers studying phytopathogenic fungi are limited. Here, we briefly review the genomic and bioinformatic resources available and the current status of fungal genomics. We also describe a relational database containing sequences of expressed sequence tags (ESTs) from three phytopathogenic fungi, Blumeria graminis, Magnaporthe grisea, and Mycosphaerella graminicola, and the methods and underlying principles required for its construction. The database contains significant annotation for each EST sequence and is accessible at http://cogeme.ex.ac.uk. An easy-to-use interface allows the user to identify gene sequences by using simple text queries or homology searches. New querying functions and large sequence sets from a variety of phytopathogenic species will be incorporated in due course.