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

A genomic compendium of cultivated human gut fungi characterizes the gut mycobiome and its relevance to common diseases

The gut fungal community represents an essential element of human health, yet its functional and metabolic potential remains insufficiently elucidated, largely due to the limited availability of reference genomes. To address this gap, we presented the cultivated gut fungi (CGF) catalog, encompassing 760 fungal genomes derived from the feces of healthy individuals. This catalog comprises 206 species spanning 48 families, including 69 species previously unidentified. We explored the functional and metabolic attributes of the CGF species and utilized this catalog to construct a phylogenetic representation of the gut mycobiome by analyzing over 11,000 fecal metagenomes from Chinese and non-Chinese populations. Moreover, we identified significant common disease-related variations in gut mycobiome composition and corroborated the associations between fungal signatures and inflammatory bowel disease (IBD) through animal experimentation. These resources and findings substantially enrich our understanding of the biological diversity and disease relevance of the human gut mycobiome.

Evolutionary history of the cytochrome P450s from Colletotrichum species and prediction of their putative functional roles during host-pathogen interactions

The genomes of species belonging to the genus Colletotrichum harbor a substantial number of cytochrome P450 monooxygenases (CYPs) encoded by a broad diversity of gene families. However, the biological role of their CYP complement (CYPome) has not been elucidated. Here, we investigated the putative evolutionary scenarios that occurred during the evolution of the CYPome belonging to the Colletotrichum Graminicola species complex (s.c.) and their biological implications. The study revealed that most of the CYPome gene families belonging to the Graminicola s.c. experienced gene contractions. The reductive evolution resulted in species restricted CYPs are predominant in each CYPome of members from the Graminicola s.c., whereas only 18 families are absolutely conserved among these species. However, members of CYP families displayed a notably different phylogenetic relationship at the tertiary structure level, suggesting a putative convergent evolution scenario. Most of the CYP enzymes of the Graminicola s.c. share redundant functions in secondary metabolite biosynthesis and xenobiotic metabolism. Hence, this current work suggests that the presence of a broad CYPome in the genus Colletotrichum plays a critical role in the optimization of the colonization capability and virulence.

Mycoparasitism related targets of Tmk1 indicate stimulating regulatory functions of this MAP kinase in Trichoderma atroviride

Mycoparasitism is a key feature of Trichoderma (Hypocreales, Ascomycota) biocontrol agents. Recent studies of intracellular signal transduction pathways of the potent mycoparasite Trichoderma atroviride revealed the involvement of Tmk1, a mitogen-activated protein kinase (MAPK), in triggering the mycoparasitic response. We previously showed that mutants missing Tmk1 exhibit reduced mycoparasitic activity against several plant pathogenic fungi. In this study, we identified the most robustly regulated targets that were governed by Tmk1 during mycoparasitism using transcriptome and proteome profiling. Tmk1 mainly exerts a stimulating function for T. atroviride during its mycoparasitic interaction with the fungal plant pathogen Rhizoctonia solani, as reflected by 89% of strongly differently responding genes in the ∆tmk1 mutant compared to the wild type. Specifically, 54% of these genes showed strong downregulation in the response with a deletion of the tmk1 gene, whereas in the wild type the same genes were strongly upregulated during the interaction with the fungal host. These included the gene encoding the mycoparasitism-related proteinase Prb1; genes involved in signal transduction pathways such as a candidate coding for a conserved 14-3-3 protein, and a gene coding for Tmk2, the T. atroviride cell-wall integrity MAP kinase; genes encoding a specific siderophore synthetase, and multiple FAD-dependent oxidoreductases and aminotransferases. Due to the phosphorylating activity of Tmk1, different (phospho-)proteomics approaches were applied and identified proteins associated with cellular metabolism, energy production, protein synthesis and fate, and cell organization. Members of FAD- and NAD/NADP-binding-domain proteins, vesicular trafficking of molecules between cellular organelles, fungal translational, as well as protein folding apparatus were among others found to be phosphorylated by Tmk1 during mycoparasitism. Outstanding downregulation in the response of the ∆tmk1 mutant to the fungal host compared to the wild type at both the transcriptome and the proteome levels was observed for nitrilase, indicating that its defense and detoxification functions might be greatly dependent on Tmk1 during T. atroviride mycoparasitism. An intersection network analysis between the identified transcripts and proteins revealed a strong involvement of Tmk1 in molecular functions with GTPase and oxidoreductase activity. These data suggest that during T. atroviride mycoparasitism this MAPK mainly governs processes regulating cell responses to extracellular signals and those involved in reactive oxygen stress.

Identification and Characterization of Novel Candidate Effector Proteins from Magnaporthe oryzae

The fungal pathogen Magnaporthe oryzae secretes a large number of effector proteins to facilitate infection, most of which are not functionally characterized. We selected potential candidate effector genes from the genome of M. oryzae, field isolate P131, and cloned 69 putative effector genes for functional screening. Utilizing a rice protoplast transient expression system, we identified that four candidate effector genes, GAS1, BAS2, MoCEP1 and MoCEP2 induced cell death in rice. In particular, MoCEP2 also induced cell death in Nicotiana benthamiana leaves through Agrobacteria-mediated transient gene expression. We further identified that six candidate effector genes, MoCEP3 to MoCEP8, suppress flg22-induced ROS burst in N. benthamiana leaves upon transient expression. These effector genes were highly expressed at a different stage after M. oryzae infection. We successfully knocked out five genes in M. oryzae, MoCEP1, MoCEP2, MoCEP3, MoCEP5 and MoCEP7. The virulence tests suggested that the deletion mutants of MoCEP2, MoCEP3 and MoCEP5 showed reduced virulence on rice and barley plants. Therefore, those genes play an important role in pathogenicity.

Ancient and pervasive expansion of adaptin-related vesicle coat machinery across Parabasalia

The eukaryotic phylum Parabasalia is composed primarily of anaerobic, endobiotic organisms such as the veterinary parasite Tritrichomonas foetus and the human parasite Trichomonas vaginalis, the latter causing the most prevalent, non-viral, sexually transmitted disease world-wide. Although a parasitic lifestyle is generally associated with a reduction in cell biology, T. vaginalis provides a striking counter-example. The 2007 T. vaginalis genome paper reported a massive and selective expansion of encoded proteins involved in vesicle trafficking, particularly those implicated in the late secretory and endocytic systems. Chief amongst these were the hetero-tetrameric adaptor proteins or 'adaptins', with T. vaginalis encoding ∼3.5 times more such proteins than do humans. The provenance of such a complement, and how it relates to the transition from a free-living or endobiotic state to parasitism, remains unclear. In this study, we performed a comprehensive bioinformatic and molecular evolutionary investigation of the heterotetrameric cargo adaptor-derived coats, comparing the molecular complement and evolution of these proteins between T. vaginalis, T. foetus and the available diversity of endobiotic parabasalids. Notably, with the recent discovery of Anaeramoeba spp. as the free-living sister lineage to all parabasalids, we were able to delve back to time points earlier in the lineage's history than ever before. We found that, although T. vaginalis still encodes the most HTAC subunits amongst parabasalids, the duplications giving rise to the complement took place more deeply and at various stages across the lineage. While some duplications appear to have convergently shaped the parasitic lineages, the largest jump is in the transition from free-living to endobiotic lifestyle with both gains and losses shaping the encoded complement. This work details the evolution of a cellular system across an important lineage of parasites and provides insight into the evolutionary dynamics of an example of expansion of protein machinery, counter to the more common trends observed in many parasitic systems.

Role of Cytochrome P450 Enzyme in Plant Microorganisms' Communication: A Focus on Grapevine

Cytochromes P450 are ancient enzymes diffused in organisms belonging to all kingdoms of life, including viruses, with the largest number of P450 genes found in plants. The functional characterization of cytochromes P450 has been extensively investigated in mammals, where these enzymes are involved in the metabolism of drugs and in the detoxification of pollutants and toxic chemicals. The aim of this work is to present an overview of the often disregarded role of the cytochrome P450 enzymes in mediating the interaction between plants and microorganisms. Quite recently, several research groups have started to investigate the role of P450 enzymes in the interactions between plants and (micro)organisms, focusing on the holobiont Vitis vinifera. Grapevines live in close association with large numbers of microorganisms and interact with each other, regulating several vine physiological functions, from biotic and abiotic stress tolerance to fruit quality at harvest.

Virulence and pathogenicity determinants in whole genome sequence of Fusarium udum causing wilt of pigeon pea

The present study deals with whole genome analysis of Fusarium udum, a wilt causing pathogen of pigeon pea. The de novo assembly identified a total of 16,179 protein-coding genes, of which 11,892 genes (73.50%) were annotated using BlastP and 8,928 genes (55.18%) from KOG annotation. In addition, 5,134 unique InterPro domains were detected in the annotated genes. Apart from this, we also analyzed genome sequence for key pathogenic genes involved in virulence, and identified 1,060 genes (6.55%) as virulence genes as per the PHI-BASE database. The secretome profiling of these virulence genes indicated the presence of 1,439 secretory proteins. Of those, an annotation of 506 predicted secretory proteins through CAZyme database indicated maximum abundance of Glycosyl hydrolase (GH, 45%) family proteins followed by auxiliary activity (AA) family proteins. Interestingly, the presence of effectors for cell wall degradation, pectin degradation, and host cell death was found. The genome comprised approximately 895,132 bp of repetitive elements, which includes 128 long terminal repeats (LTRs), and 4,921 simple sequence repeats (SSRs) of 80,875 bp length. The comparative mining of effector genes among different Fusarium species revealed five common and two specific effectors in F. udum that are related to host cell death. Furthermore, wet lab experiment validated the presence of effector genes like SIX (for Secreted in Xylem). We conclude that deciphering the whole genome of F. udum would be instrumental in understanding evolution, virulence determinants, host-pathogen interaction, possible control strategies, ecological behavior, and many other complexities of the pathogen.

Genome-Wide Analysis of Cytochrome P450s of Alternaria Species: Evolutionary Origin, Family Expansion and Putative Functions

Cytochrome P450s are a group of monooxygenase enzymes involved in primary, secondary and xenobiotic metabolisms. They have a wide application in the agriculture sector where they could serve as a target for herbicides or fungicides, while they could function in the pharmaceutical industry as drugs or drugs structures or for bioconversions. Alternaria species are among the most commonly encountered fungal genera, with most of them living as saprophytes in different habitats, while others are parasites of plants and animals. This study was conducted to elucidate the diversity and abundance, evolutionary relationships and cellular localization of 372 cytochrome P450 in 13 Alternaria species. The 372 CYP proteins were phylogenetically clustered into ten clades. Forty (40) clans and seventy-one (71) cyp families were identified, of which eleven (11) families were found to appear in one species each. The majority of the CYP proteins were located in the endomembrane system. Polyketide synthase (PKS) gene cluster was the predominant secondary metabolic-related gene cluster in all the Alternaria species studied, except in A. porriof, where non-ribosomal peptide synthetase genes were dominant. This study reveals the expansion of cyps in these fungal genera, evident in the family and clan expansions, which is usually associated with the evolution of fungal characteristics, especially their lifestyle either as parasites or saprophytes, with the ability to metabolize a wide spectrum of substrates. This study can be used to understand the biology, physiology and toxigenic potentials of P450 in these fungal genera.

Phenolic degradation by catechol dioxygenases is associated with pathogenic fungi with a necrotrophic lifestyle in the Ceratocystidaceae

Fungal species of the Ceratocystidaceae grow on their host plants using a variety of different lifestyles, from saprophytic to highly pathogenic. Although many genomes of fungi in the Ceratocystidaceae are publicly available, it is not known how the genes that encode catechol dioxygenases (CDOs), enzymes involved in the degradation of phenolic plant defense compounds, differ among members of the Ceratocystidaceae. The aim of this study was therefore to identify and characterize the genes encoding CDOs in the genomes of Ceratocystidaceae representatives. We found that genes encoding CDOs are more abundant in pathogenic necrotrophic species of the Ceratocystidaceae and less abundant in saprophytic species. The loss of the CDO genes and the associated 3-oxoadipate catabolic pathway appears to have occurred in a lineage-specific manner. Taken together, this study revealed a positive association between CDO gene copy number and fungal lifestyle in Ceratocystidaceae representatives.

Expanding the Biological Role of Lipo-Chitooligosaccharides and Chitooligosaccharides in Laccaria bicolor Growth and Development

The role of lipo-chitooligosaccharides (LCOs) as signaling molecules that mediate the establishment of symbiotic relationships between fungi and plants is being redefined. New evidence suggests that the production of these molecular signals may be more of a common trait in fungi than what was previously thought. LCOs affect different aspects of growth and development in fungi. For the ectomycorrhizal forming fungi, Laccaria bicolor, the production and effects of LCOs have always been studied with a symbiotic plant partner; however, there is still no scientific evidence describing the effects that these molecules have on this organism. Here, we explored the physiological, molecular, and metabolomic changes in L. bicolor when grown in the presence of exogenous sulfated and non-sulfated LCOs, as well as the chitooligomers, chitotetraose (CO4), and chitooctaose (CO8). Physiological data from 21 days post-induction showed reduced fungal growth in response to CO and LCO treatments compared to solvent controls. The underlying molecular changes were interrogated by proteomics, which revealed substantial alterations to biological processes related to growth and development. Moreover, metabolite data showed that LCOs and COs caused a downregulation of organic acids, sugars, and fatty acids. At the same time, exposure to LCOs resulted in the overproduction of lactic acid in L. bicolor. Altogether, these results suggest that these signals might be fungistatic compounds and contribute to current research efforts investigating the emerging impacts of these molecules on fungal growth and development.

The role of filamentous fungi in advancing the development of a sustainable circular bioeconomy

Human activities generate enormous amounts of organic wastes and residues. Filamentous fungi (FF) are able to grow on a broad range of substrates and survive over a wide spectrum of growth conditions. These characteristics enable FF to be exploited in biorefineries for various waste streams. Valorization of food industry byproducts into biomass and various arrays of value-added products using FF creates promising pathways toward a sustainable circular economy. This approach might also contribute to reaching the sustainable development goals set by the United Nations, particularly for zero hunger as well as affordable and clean energy. This paper presents the application of filamentous fungi in food, feeds, fuels, biochemicals, and biopolymers. The nutritional values, health benefits, and safety of foods derived from byproducts of food industries are also addressed. The technoeconomical feasibilities, sustainability aspects and challenges and future perspectives for biorefineries using filamentous fungi are discussed.

Genetic response to nitrogen starvation in the aggressive Eucalyptus foliar pathogen Teratosphaeria destructans

Teratosphaeria destructans is one of the most aggressive foliar pathogens of Eucalyptus. The biological factors underpinning T. destructans infections, which include shoot and leaf blight on young trees, have never been interrogated. Thus, the means by which the pathogen modifies its host environment to overcome host defences remain unknown. By applying transcriptome sequencing, the aim of this study was to compare gene expression in a South African isolate of T. destructans grown on nitrogen-deficient and complete media. This made it possible to identify upregulated genes in a nitrogen-starved environment, often linked to the pathogenicity of the fungus. The results support the hypothesis that nitrogen starvation in T. destructans likely mirrors an in planta genetic response. This is because 45% of genes that were highly upregulated under nitrogen starvation have previously been reported to be associated with infection in other pathogen systems. These included several CAZymes, fungal effector proteins, peptidases, kinases, toxins, lipases and proteins associated with detoxification of toxic compounds. Twenty-five secondary metabolites were identified and expressed in both nitrogen-deficient and complete conditions. Additionally, the most highly expressed genes in both growth conditions had pathogenicity-related functions. This study highlights the large number of expressed genes associated with pathogenicity and overcoming plant defences. As such, the generated baseline knowledge regarding pathogenicity and aggressiveness in T. destructans is a valuable reference for future in planta work.

Molecular and comparative genomic analyses reveal evolutionarily conserved and unique features of the Schizosaccharomyces japonicus mycelial growth and the underlying genomic changes

Fungal pathogens, from phytopathogenic fungus to human pathogens, are able to alternate between the yeast-like form and filamentous forms. This morphological transition (dimorphism) is in close connection with their pathogenic lifestyles and with their responses to changing environmental conditions. The mechanisms governing these morphogenetic conversions are still not fully understood. Therefore, we studied the filamentous growth of the less-known, non-pathogenic dimorphic fission yeast, S. japonicus, which belongs to an ancient and early evolved branch of the Ascomycota. Its RNA sequencing revealed that several hundred genes were up- or down-regulated in the hyphae compared to the yeast-phase cells. These genes belonged to different GO categories, confirming that mycelial growth is a rather complex process. The genes of transport- and metabolic processes appeared especially in high numbers among them. High expression of genes involved in glycolysis and ethanol production was found in the hyphae, while other results pointed to the regulatory role of the protein kinase A (PKA) pathway. The homologues of 49 S. japonicus filament-associated genes were found by sequence alignments also in seven distantly related dimorphic and filamentous species. The comparative genomic analyses between S. japonicus and the closely related but non-dimorphic S. pombe shed some light on the differences in their genomes. All these data can contribute to a better understanding of hyphal growth and those genomic rearrangements that underlie it.

Genomic characterization of three marine fungi, including Emericellopsis atlantica sp. nov. with signatures of a generalist lifestyle and marine biomass degradation

Marine fungi remain poorly covered in global genome sequencing campaigns; the 1000 fungal genomes (1KFG) project attempts to shed light on the diversity, ecology and potential industrial use of overlooked and poorly resolved fungal taxa. This study characterizes the genomes of three marine fungi: Emericellopsis sp. TS7, wood-associated Amylocarpus encephaloides and algae-associated Calycina marina. These species were genome sequenced to study their genomic features, biosynthetic potential and phylogenetic placement using multilocus data. Amylocarpus encephaloides and C. marina were placed in the Helotiaceae and Pezizellaceae (Helotiales), respectively, based on a 15-gene phylogenetic analysis. These two genomes had fewer biosynthetic gene clusters (BGCs) and carbohydrate active enzymes (CAZymes) than Emericellopsis sp. TS7 isolate. Emericellopsis sp. TS7 (Hypocreales, Ascomycota) was isolated from the sponge Stelletta normani. A six-gene phylogenetic analysis placed the isolate in the marine Emericellopsis clade and morphological examination confirmed that the isolate represents a new species, which is described here as E. atlantica. Analysis of its CAZyme repertoire and a culturing experiment on three marine and one terrestrial substrates indicated that E. atlantica is a psychrotrophic generalist fungus that is able to degrade several types of marine biomass. FungiSMASH analysis revealed the presence of 35 BGCs including, eight non-ribosomal peptide synthases (NRPSs), six NRPS-like, six polyketide synthases, nine terpenes and six hybrid, mixed or other clusters. Of these BGCs, only five were homologous with characterized BGCs. The presence of unknown BGCs sets and large CAZyme repertoire set stage for further investigations of E. atlantica. The Pezizellaceae genome and the genome of the monotypic Amylocarpus genus represent the first published genomes of filamentous fungi that are restricted in their occurrence to the marine habitat and form thus a valuable resource for the community that can be used in studying ecological adaptions of fungi using comparative genomics.

Diversity and Function of Appressoria

Endophytic, saprobic, and pathogenic fungi have evolved elaborate strategies to obtain nutrients from plants. Among the diverse plant-fungi interactions, the most crucial event is the attachment and penetration of the plant surface. Appressoria, specialized infection structures, have been evolved to facilitate this purpose. In this review, we describe the diversity of these appressoria and classify them into two main groups: single-celled appressoria (proto-appressoria, hyaline appressoria, melanized (dark) appressoria) and compound appressoria. The ultrastructure of appressoria, their initiation, their formation, and their function in fungi are discussed. We reviewed the molecular mechanisms regulating the formation and function of appressoria, their strategies to evade host defenses, and the related genomics and transcriptomics. The current review provides a foundation for comprehensive studies regarding their evolution and diversity in different fungal groups.

Historical genomics reveals the evolutionary mechanisms behind multiple outbreaks of the host-specific coffee wilt pathogen Fusarium xylarioides

BACKGROUND

Nearly 50% of crop yields are lost to pests and disease, with plants and pathogens locked in an amplified co-evolutionary process of disease outbreaks. Coffee wilt disease, caused by Fusarium xylarioides, decimated coffee production in west and central Africa following its initial outbreak in the 1920s. After successful management, it later re-emerged and by the 2000s comprised two separate epidemics on arabica coffee in Ethiopia and robusta coffee in east and central Africa.

RESULTS

Here, we use genome sequencing of six historical culture collection strains spanning 52 years to identify the evolutionary processes behind these repeated outbreaks. Phylogenomic reconstruction using 13,782 single copy orthologs shows that the robusta population arose from the initial outbreak, whilst the arabica population is a divergent sister clade to the other strains. A screen for putative effector genes involved in pathogenesis shows that the populations have diverged in gene content and sequence mainly by vertical processes within lineages. However, 15 putative effector genes show evidence of horizontal acquisition, with close homology to genes from F. oxysporum. Most occupy small regions of homology within wider scaffolds, whereas a cluster of four genes occupy a 20Kb scaffold with strong homology to a region on a mobile pathogenicity chromosome in F. oxysporum that houses known effector genes. Lacking a match to the whole mobile chromosome, we nonetheless found close associations with DNA transposons, especially the miniature impala type previously proposed to facilitate horizontal transfer of pathogenicity genes in F. oxysporum. These findings support a working hypothesis that the arabica and robusta populations partly acquired distinct effector genes via transposition-mediated horizontal transfer from F. oxysporum, which shares coffee as a host and lives on other plants intercropped with coffee.

CONCLUSION

Our results show how historical genomics can help reveal mechanisms that allow fungal pathogens to keep pace with our efforts to resist them. Our list of putative effector genes identifies possible future targets for fungal control. In turn, knowledge of horizontal transfer mechanisms and putative donor taxa might help to design future intercropping strategies that minimize the risk of transfer of effector genes between closely-related Fusarium taxa.

Comparative Genomics: Insights on the Pathogenicity and Lifestyle of Rhizoctonia solani

Proper management of agricultural disease is important to ensure sustainable food security. Staple food crops like rice, wheat, cereals, and other cash crops hold great export value for countries. Ensuring proper supply is critical; hence any biotic or abiotic factors contributing to the shortfall in yield of these crops should be alleviated. Rhizoctonia solani is a major biotic factor that results in yield losses in many agriculturally important crops. This paper focuses on genome informatics of our Malaysian Draft R. solani AG1-IA, and the comparative genomics (inter- and intra- AG) with four AGs including China AG1-IA (AG1-IA_KB317705.1), AG1-IB, AG3, and AG8. The genomic content of repeat elements, transposable elements (TEs), syntenic genomic blocks, functions of protein-coding genes as well as core orthologous genic information that underlies R. solani’s pathogenicity strategy were investigated. Our analyses show that all studied AGs have low content and varying profiles of TEs. All AGs were dominant for Class I TE, much like other basidiomycete pathogens. All AGs demonstrate dominance in Glycoside Hydrolase protein-coding gene assignments suggesting its importance in infiltration and infection of host. Our profiling also provides a basis for further investigation on lack of correlation observed between number of pathogenicity and enzyme-related genes with host range. Despite being grouped within the same AG with China AG1-IA, our Draft AG1-IA exhibits differences in terms of protein-coding gene proportions and classifications. This implies that strains from similar AG do not necessarily have to retain similar proportions and classification of TE but must have the necessary arsenal to enable successful infiltration and colonization of host. In a larger perspective, all the studied AGs essentially share core genes that are generally involved in adhesion, penetration, and host colonization. However, the different infiltration strategies will depend on the level of host resilience where this is clearly exhibited by the gene sets encoded for the process of infiltration, infection, and protection from host.

Chorismate mutase and isochorismatase, two potential effectors of the migratory nematode Hirschmanniella oryzae, increase host susceptibility by manipulating secondary metabolite content of rice

Hirschmanniella oryzae is one of the most devastating nematodes on rice, leading to substantial yield losses. Effector proteins aid the nematode during the infection process by subduing plant defence responses. In this research we characterized two potential H. oryzae effector proteins, chorismate mutase (HoCM) and isochorismatase (HoICM), and investigated their enzymatic activity and their role in plant immunity. Both HoCM and HoICM proved to be enzymatically active in complementation tests in mutant Escherichia coli strains. Infection success by the migratory nematode H. oryzae was significantly higher in transgenic rice lines constitutively expressing HoCM or HoICM. Expression of HoCM, but not HoICM, increased rice susceptibility against the sedentary nematode Meloidogyne graminicola also. Transcriptome and metabolome analyses indicated reductions in secondary metabolites in the transgenic rice plants expressing the potential nematode effectors. The results presented here demonstrate that both HoCM and HoICM suppress the host immune system and that this may be accomplished by lowering secondary metabolite levels in the plant.

"CATAStrophy," a Genome-Informed Trophic Classification of Filamentous Plant Pathogens - How Many Different Types of Filamentous Plant Pathogens Are There?

The traditional classification of fungal and oomycete phytopathogens into three classes – biotrophs, hemibiotrophs, or necrotrophs – is unsustainable. This study highlights multiple phytopathogen species for which these labels have been inappropriately applied. We propose a novel and reproducible classification based solely on genome-derived analysis of carbohydrate-active enzyme (CAZyme) gene content called CAZyme-Assisted Training And Sorting of -trophy (CATAStrophy). CATAStrophy defines four major divisions for species associated with living plants. These are monomertrophs (Mo) (corresponding to biotrophs), polymertrophs (P) (corresponding to necrotrophs), mesotrophs (Me) (corresponding to hemibiotrophs), and vasculartrophs (including species commonly described as wilts, rots, or anthracnoses). The Mo class encompasses symbiont, haustorial, and non-haustorial species. Me are divided into the subclasses intracellular and extracellular Me, and the P into broad and narrow host sub-classes. This gives a total of seven discrete plant-pathogenic classes. The classification provides insight into the properties of these species and offers a facile route to develop control measures for newly recognized diseases. Software for CATAStrophy is available online at https://github.com/ccdmb/catastrophy. We present the CATAStrophy method for the prediction of trophic phenotypes based on CAZyme gene content, as a complementary method to the traditional tripartite “biotroph–hemibiotroph–necrotroph” classifications that may encourage renewed investigation and revision within the fungal biology community.

Expression and role of CYP505A1 in pathogenicity of Fusarium oxysporum f. sp. lactucae

BACKGROUND

Cytochrome P450 enzymes (CYPs) are monooxygenases present in every domain of life. In fungi CYPs are involved in virulence. Fusarium wilt of lettuce, caused by F. oxysporum f. sp. lactucae, is the most serious disease of lettuce. F. oxysporum f.sp. lactucae MSA35 is an antagonistic fungus. Pathogenic formae specialis of F. oxysporum possess a CYP belonging to the new family CYP505. This enzyme hydroxylates saturated fatty acids that play a role in plant defence.

METHODS

Molecular tools were adopted to search for cyp505 gene in MSA35 genome. cyp505 gene expression analysis in pathogenic and antagonistic Fusarium was performed. The enzyme was expressed in its recombinant form and used for catalytic reactions with fatty acids, the products of which were characterized by mass spectrometry analysis.

RESULTS

A novel MSA35 self-sufficient CYP505 is differentially expressed in antagonistic and pathogenic F. oxysporum. Its expression is induced by the host plant lettuce in both pathogenesis and antagonism during the early phase of the interaction, while it is silenced during the late phase only in antagonistic Fusarium. Mass-spectrometry investigations proved that CYP505A1 mono-hydroxylates lauric, palmitic and stearic acids.

CONCLUSIONS

The ability of CYP505A1 to oxidize fatty acids present in the cortical cell membranes together with its differential expression in its Fusarium antagonistic form point out to the possibility that this enzyme is associated with Fusarium pathogenicity in lettuce.

GENERAL SIGNIFICANCE

The CYP505 clan is present in pathogenic fungal phyla, making CYP505A1 enzyme a putative candidate as a new target for the development of novel antifungal molecules.

Detection of putative pathogenicity and virulence genes of Erysiphe pisi using genome-wide in-silico search and their suppression by er2 mediated resistance in garden pea

The biotrophic fungus, Erysiphe pisi is the chief causal agent of powdery mildew disease of garden pea. A genome-wide search using in-silico approach was carried to detect putative pathogenicity and virulence genes of E. pisi, since information about these genes and their interaction with pea is limited. Nineteen putative pathogenicity gene sequences were detected through genome-wide pathogenicity gene-search and confirmed them to be conserved in E. pisi through genomic PCRs. Fifteen of these genes expressed through reverse transcriptase-polymerase chain reaction (RT-PCR) amplifying expected band size along with fungal and plant specific internal controls. Gene sequencing and annotation revealed them to be Erysiphe-specific. A time course study was carried to monitor expression of nine of these genes through real-time quantitative (qRT)-PCR in Erysiphe-challenged plants of powdery mildew resistant pea genotype, JI-2480 carrying er2 gene and susceptible pea cultivar, Arkel. Expression of these genes was differentially and temporally regulated. They were found mostly related to signaling; cAMP-PKA (cPKA, CRP and AC) and MAPK (MST7) pathways along with MFP, TRE and PEX which are reported pathogenicity factors in other ascomycete members indicating that similar conserved pathways function in E. pisi also. These genes expressed at higher level at initial hours post inoculation (hpi) as early as 6 hpi in Arkel compared to JI-2480 implying them as pathogenicity factors. The elevated level of expression of MFP, TRE, CRP and cPKA gene sequences in E. pisi-challenged JI-2480 genotype at 12 hpi alone suggests these genes to possess a role in avirulence in JI-2480, conferring er2 mediated resistance.

Chitin-Binding Protein of Verticillium nonalfalfae Disguises Fungus from Plant Chitinases and Suppresses Chitin-Triggered Host Immunity

During fungal infections, plant cells secrete chitinases, which digest chitin in the fungal cell walls. The recognition of released chitin oligomers via lysin motif (LysM)-containing immune host receptors results in the activation of defense signaling pathways. We report here that Verticillium nonalfalfae, a hemibiotrophic xylem-invading fungus, prevents these digestion and recognition processes by secreting a carbohydrate-binding motif 18 (CBM18)-chitin-binding protein, VnaChtBP, which is transcriptionally activated specifically during the parasitic life stages. VnaChtBP is encoded by the Vna8.213 gene, which is highly conserved within the species, suggesting high evolutionary stability and importance for the fungal lifestyle. In a pathogenicity assay, however, Vna8.213 knockout mutants exhibited wilting symptoms similar to the wild-type fungus, suggesting that Vna8.213 activity is functionally redundant during fungal infection of hop. In a binding assay, recombinant VnaChtBP bound chitin and chitin oligomers in vitro with submicromolar affinity and protected fungal hyphae from degradation by plant chitinases. Moreover, the chitin-triggered production of reactive oxygen species from hop suspension cells was abolished in the presence of VnaChtBP, indicating that VnaChtBP also acts as a suppressor of chitin-triggered immunity. Using a yeast-two-hybrid assay, circular dichroism, homology modeling, and molecular docking, we demonstrated that VnaChtBP forms dimers in the absence of ligands and that this interaction is stabilized by the binding of chitin hexamers with a similar preference in the two binding sites. Our data suggest that, in addition to chitin-binding LysM (CBM50) and Avr4 (CBM14) fungal effectors, structurally unrelated CBM18 effectors have convergently evolved to prevent hydrolysis of the fungal cell wall against plant chitinases and to interfere with chitin-triggered host immunity.

Responses to Hydric Stress in the Seed-Borne Necrotrophic Fungus Alternaria brassicicola

Alternaria brassicicola is a necrotrophic fungus causing black spot disease and is an economically important seed-borne pathogen of cultivated brassicas. Seed transmission is a crucial component of its parasitic cycle as it promotes long-term survival and dispersal. Recent studies, conducted with the Arabidopsis thaliana/A. brassicicola pathosystem, showed that the level of susceptibility of the fungus to water stress strongly influenced its seed transmission ability. In this study, we gained further insights into the mechanisms involved in the seed infection process by analyzing the transcriptomic and metabolomic responses of germinated spores of A. brassicicola exposed to water stress. Then, the repertoire of putative hydrophilins, a group of proteins that are assumed to be involved in cellular dehydration tolerance, was established in A. brassicicola based on the expression data and additional structural and biochemical criteria. Phenotyping of single deletion mutants deficient for fungal hydrophilin-like proteins showed that they were affected in their transmission to A. thaliana seeds, although their aggressiveness on host vegetative tissues remained intact.

Blast Fungal Genomes Show Frequent Chromosomal Changes, Gene Gains and Losses, and Effector Gene Turnover

Pyricularia is a fungal genus comprising several pathogenic species causing the blast disease in monocots. Pyricularia oryzae, the best-known species, infects rice, wheat, finger millet, and other crops. As past comparative and population genomics studies mainly focused on isolates of P. oryzae, the genomes of the other Pyricularia species have not been well explored. In this study, we obtained a chromosomal-level genome assembly of the finger millet isolate P. oryzae MZ5-1-6 and also highly contiguous assemblies of Pyricularia sp. LS, P. grisea, and P. pennisetigena. The differences in the genomic content of repetitive DNA sequences could largely explain the variation in genome size among these new genomes. Moreover, we found extensive gene gains and losses and structural changes among Pyricularia genomes, including a large interchromosomal translocation. We searched for homologs of known blast effectors across fungal taxa and found that most avirulence effectors are specific to Pyricularia, whereas many other effectors share homologs with distant fungal taxa. In particular, we discovered a novel effector family with metalloprotease activity, distinct from the well-known AVR-Pita family. We predicted 751 gene families containing putative effectors in 7 Pyricularia genomes and found that 60 of them showed differential expression in the P. oryzae MZ5-1-6 transcriptomes obtained under experimental conditions mimicking the pathogen infection process. In summary, this study increased our understanding of the structural, functional, and evolutionary genomics of the blast pathogen and identified new potential effector genes, providing useful data for developing crops with durable resistance.

Brassica napus Infected with Leptosphaeria maculans

Alternative splicing (AS) is a post-transcriptional regulatory process that enhances transcriptome diversity, thereby affecting plant growth, development, and stress responses. To identify the new transcripts and changes in the isoform-level AS landscape of rapeseed (Brassica napus) infected with the fungal pathogen Leptosphaeria maculans, we compared eight RNA-seq libraries prepared from mock-inoculated and inoculated B. napus cotyledons and stems. The AS events that occurred in stems were almost the same as those in cotyledons, with intron retention representing the most common AS pattern. We identified 1892 differentially spliced genes between inoculated and uninoculated plants. We performed a weighted gene co-expression network analysis (WGCNA) to identify eight co-expression modules and their Hub genes, which are the genes most connected with other genes within each module. There are nine Hub genes, encoding nine transcription factors, which represent key regulators of each module, including members of the NAC, WRKY, TRAF, AP2/ERF-ERF, C2H2, C2C2-GATA, HMG, bHLH, and C2C2-CO-like families. Finally, 52 and 117 alternatively spliced genes in cotyledons and stems were also differentially expressed between mock-infected and infected materials, such as HMG and C2C2-Dof; which have dual regulatory mechanisms in response to L. maculans. The splicing of the candidate genes identified in this study could be exploited to improve resistance to L. maculans.

Complementary Proteomics, Genomics approaches identifies potential pathogenicity/virulence factors in Tilletia indica induced under the influence of host factor

Karnal bunt disease of wheat is incited by quarantine fungal pathogen T. indica. Till date, there is little information on the pathogenic mechanisms involved in Karnal bunt. In order to understand the molecular mechanisms of disease pathogenesis, highly aggressive T. indica TiK isolate was cultured in the presence of host factor extracted from developing spikes of wheat variety WH-542. Modulation in protein profile of mycelial proteins and secretome from TiK cultured in the absence and presence of host factor was analyzed by 2-DE. Fifteen and twenty nine protein spots were up-regulated/differentially regulated in the proteome of mycelial and secreted proteins, respectively and identified using MALDI-TOF/TOF. Identified proteins are involved in suppression of host defense responses, lignin degradation of plant cell wall, penetration, adhesion of pathogen to host tissues, pathogen mediated reactive oxygen species generation, hydrolytic enzymes, detoxification of host generated reactive oxygen species. Further, integration of proteomic and genomic analysis has led to candidate pathogenicity/virulence factors identification. They were functionally annotated by sequence as well as structure based analysis. In this study, complementation of proteomics and genomics approaches resulted in novel pathogenicity/virulence factor(s) identification in T. indica.

A positive-charged patch and stabilized hydrophobic core are essential for avirulence function of AvrPib in the rice blast fungus

Fungal avirulence effectors, a key weapon utilized by pathogens to promote their infection, are recognized by immune receptors to boost host R gene-mediated resistance. Many avirulence effectors share sparse sequence homology to proteins with known functions, and their molecular and biochemical functions together with the evolutionary relationship among different members remain largely unknown. Here, the crystal structure of AvrPib, an avirulence effector from Magnaporthe oryzae, was determined and showed a high degree of similarity to the M. oryzae Avrs and ToxB (MAX) effectors. Compared with other MAX effectors, AvrPib has a distinct positive-charge patch formed by five positive-charged residues (K29, K30, R50, K52 and K70) on the surface. These five key residues were essential to avirulence function of AvrPib and affected its nuclear localization into host cells. Moreover, residues V39 and V58, which locate in the hydrophobic core of the structure, cause loss of function of AvrPib by single-point mutation in natural isolates. In comparison with the wild-type AvrPib, the V39A or V58A mutations resulted in a partial or entire loss of secondary structure elements. Taken together, our results suggest that differences in the surface charge distribution of avirulence proteins could be one of the major bases for the variation in effector-receptor specificity, and that destabilization of the hydrophobic core is one of the major mechanisms employed by AvrPib for the fungus to evade recognition by resistance factors in the host cell.

Genome-wide analysis of cytochrome P450s of Trichoderma spp.: annotation and evolutionary relationships

Background

Cytochrome P450s form an important group of enzymes involved in xenobiotics degradation and metabolism, both primary and secondary. These enzymes are also useful in industry as biotechnological tools for bioconversion and a few are reported to be involved in pathogenicity. Trichoderma spp. are widely used in industry and agriculture and are known for their biosynthetic potential of a large number of secondary metabolites. For realising the full biosynthetic potential of an organism, it is important to do a genome-wide annotation and cataloguing of these enzymes.

Results

Here, we have studied the genomes of seven species (T. asperellum, T. atroviride, T. citrinoviride, T. longibrachiatum, T. reesei , T. harzianum and T. virens) and identified a total of 477 cytochrome P450s. We present here the classification, evolution and structure as well as predicted function of these proteins. This study would pave the way for functional characterization of these groups of enzymes and will also help in realization of their full economic potential.

Conclusion

Our CYPome annotation and evolutionary studies of the seven Trichoderma species now provides opportunities for exploration of research-driven strategies to select Trichoderma species for various applications especially in relation to secondary metabolism and degradation of environmental pollutants.

How a fungus shapes biotechnology: 100 years of Aspergillus niger research

In 1917, a food chemist named James Currie made a promising discovery: any strain of the filamentous mould Aspergillus niger would produce high concentrations of citric acid when grown in sugar medium. This tricarboxylic acid, which we now know is an intermediate of the Krebs cycle, had previously been extracted from citrus fruits for applications in food and beverage production. Two years after Currie’s discovery, industrial-level production using A. niger began, the biochemical fermentation industry started to flourish, and industrial biotechnology was born. A century later, citric acid production using this mould is a multi-billion dollar industry, with A. niger additionally producing a diverse range of proteins, enzymes and secondary metabolites. In this review, we assess main developments in the field of A. niger biology over the last 100 years and highlight scientific breakthroughs and discoveries which were influential for both basic and applied fungal research in and outside the A. niger community. We give special focus to two developments of the last decade: systems biology and genome editing. We also summarize the current international A. niger research community, and end by speculating on the future of fundamental research on this fascinating fungus and its exploitation in industrial biotechnology.

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

The first eukaryotic genome to be sequenced was fungal, and there continue to be more sequenced genomes in the kingdom Fungi than in any other eukaryotic kingdom. Comparison of these genomes reveals many sources of genetic variation, from single nucleotide polymorphisms to horizontal gene transfer and on to changes in the arrangement and number of chromosomes, not to mention endofungal bacteria and viruses. Population genomics shows that all sources generate variation all the time and implicate natural selection as the force maintaining genome stability. Variation in wild populations is a rich resource for associating genetic variation with phenotypic variation, whether through quantitative trait locus mapping, genome-wide association studies, or reverse ecology. Subjects of studies associating genetic and phenotypic variation include model fungi, e.g., Saccharomyces and Neurospora, but pioneering studies have also been made with fungi pathogenic to plants, e.g., Pyricularia (= Magnaporthe), Zymoseptoria, and Fusarium, and to humans, e.g., Coccidioides, Cryptococcus, and Candida.

Genome sequencing and comparative genomics reveal the potential pathogenic mechanism of Cercospora sojina Hara on soybean

Frogeye leaf spot, caused by Cercospora sojina Hara, is a common disease of soybean in most soybean-growing countries of the world. In this study, we report a high-quality genome sequence of C. sojina by Single Molecule Real-Time sequencing method. The 40.8-Mb genome encodes 11,655 predicated genes, and 8,474 genes are revealed by RNA sequencing. Cercospora sojina genome contains large numbers of gene clusters that are involved in synthesis of secondary metabolites, including mycotoxins and pigments. However, much less carbohydrate-binding module protein encoding genes are identified in C. sojina genome, when compared with other phytopathogenic fungi. Bioinformatics analysis reveals that C. sojina harbours about 752 secreted proteins, and 233 of them are effectors. During early infection, the genes for metabolite biosynthesis and effectors are significantly enriched, suggesting that they may play essential roles in pathogenicity. We further identify 13 effectors that can inhibit BAX-induced cell death. Taken together, our results provide insights into the infection mechanisms of C. sojina on soybean.

A spatial temporal analysis of the Fusarium graminearum transcriptome during symptomless and symptomatic wheat infection

Fusarium head blight of wheat is one of the most serious and hazardous crop diseases worldwide. Here, a transcriptomic investigation of Fusarium graminearum reveals a new model for symptomless and symptomatic wheat infection. The predicted metabolic state and secretome of F. graminearum were distinct within symptomless and symptomatic wheat tissues. Transcripts for genes involved in the biosynthesis of the mycotoxin, deoxynivalenol, plus other characterized and putative secondary metabolite clusters increased in abundance in symptomless tissue. Transcripts encoding for genes of distinct groups of putative secreted effectors increased within either symptomless or symptomatic tissue. Numerous pathogenicity-associated gene transcripts and transcripts representing PHI-base mutations that impacted on virulence increased in symptomless tissue. In contrast, hydrolytic carbohydrate-active enzyme (CAZyme) and lipase gene transcripts exhibited a different pattern of expression, resulting in elevated transcript abundance during the development of disease symptoms. Genome-wide comparisons with existing datasets confirmed that, within the wheat floral tissue, at a single time point, different phases of infection co-exist, which are spatially separated and reminiscent of both early and late infection. This study provides novel insights into the combined spatial temporal coordination of functionally characterized and hypothesized virulence strategies.

Evolutionary Analysis of Pectin Lyases of the Genus Colletotrichum

Pectin lyases (PNLs) are important enzymes that are involved in plant cell wall degradation during the infection process. Colletotrichum is a diverse genus of fungi, which allows the study of the evolution of PNLs and their possible role in pathogen-host interactions and lifestyle adaptations. The phylogenetic reconstruction of PNLs from Colletotrichum and analysis of selection pressures showed the formation of protein lineages by groups of species with different selection pressures and specific patterns. The analysis of positive selection at individual sites using different methods allowed for the identification of three codons with evidence of positive selection in the oligosaccharide-binding region and two codons on the antiparallel sheet, which may influence the interaction with the substrate. Seven codons on the surface of the protein, mainly in the peripheral helices of the PNLs, could have an important function in evasion of plant defenses, as has been proposed in other enzymes. According to our results, it is possible that events of genetic duplication occurred in ancestral lines, followed by episodes of genetic diversification and gene loss, probably influenced by differences in the composition of the host cell wall. Additionally, different patterns of evolution in Colletotrichum appear to be molded by a strong purifying selection and positive selection episodes that forged the observed evolutionary patterns, possibly influenced by host interaction or substrate specificity. This work represents a starting point for the study of sites that may be important for evasion of plant defenses and biotechnological purposes.

Plant Pathogenic Fungi

Fungi are among the dominant causal agents of plant diseases. To colonize plants and cause disease, pathogenic fungi use diverse strategies. Some fungi kill their hosts and feed on dead material (necrotrophs), while others colonize the living tissue (biotrophs). For successful invasion of plant organs, pathogenic development is tightly regulated and specialized infection structures are formed. To further colonize hosts and establish disease, fungal pathogens deploy a plethora of virulence factors. Depending on the infection strategy, virulence factors perform different functions. While basically all pathogens interfere with primary plant defense, necrotrophs secrete toxins to kill plant tissue. In contrast, biotrophs utilize effector molecules to suppress plant cell death and manipulate plant metabolism in favor of the pathogen. This article provides an overview of plant pathogenic fungal species and the strategies they use to cause disease.

The Nep1-like protein family of Magnaporthe oryzae is dispensable for the infection of rice plants

The necrosis- and ethylene-inducing protein 1 (Nep1)-like proteins (NLPs) are a class of microbe-associated molecular patterns widely distributed across diverse groups of plant-associated microorganisms. In spite of the cytotoxic activity in dicot plants, the role of most NLPs in the virulence of plant pathogens is still largely unknown. We showed that the MoNLP family of rice blast fungus varied very little in amino acid sequence, transient expression of three MoNLPs induced cell death and the production of reactive oxygen species in Nicotiana benthamiana, and the expression of MoNLPs was induced during infection of susceptible rice plants. To further investigate the biological role of the MoNLP family, a marker-free gene replacement vector was developed and used to knock out the whole family in Magnaporthe oryzae. Results showed no significant difference in disease levels caused by wild type and the quadruple ΔMoNLP mutant strains. Likewise, the sporulation and radial growth of the two strains were similar under various unfavorable cultural conditions including malnutrition and abiotic stresses. These observations demonstrated that the MoNLP family is dispensable for the fungal tolerance to the tested adverse cultural conditions, and more importantly, for the virulence of blast fungus on susceptible rice plants.

Verticillium dahliae's Isochorismatase Hydrolase Is a Virulence Factor That Contributes to Interference With Potato's Salicylate and Jasmonate Defense Signaling

This study aimed to dissect the function of the Isochorismatase Hydrolase (ICSH1) gene in Verticillium dahliae's pathogenesis on potato. VdICSH1 was up-regulated in V. dahliae after induction with extracts from potato tissues. Its expression increased more in response to root extracts than to leaf and stem extracts. However, such expression in response to root extracts was not significantly different in the highly and weakly aggressive isolates tested. During infection of detached potato leaves, VdICSH1 expression increased significantly in the highly aggressive isolate compared to the weakly aggressive one. We generated icsh1 mutants from a highly aggressive isolate of V. dahliae and compared their pathogenicity with that of the original wild type strain. The analysis showed that this gene is required for full virulence of V. dahliae on potatoes. When we previously found differential accumulation of ICSH1 protein in favor of the highly aggressive isolate, as opposed to the weakly aggressive one, we had hypothesized that ICSH would interfere with the host's defense SA-based signaling. Here, we measured the accumulation of both salicylic acid (SA) and jasmonic acid (JA) in potato plants inoculated with an icsh1 mutant in comparison with the wild type strain. The higher accumulation of bound SA in the leaves in response to the icsh1 mutant compared to the wild type confirms the hypothesis that ICSH1 interferes with SA. However, the different trends in SA and JA accumulation in potato in the roots and in the stems at the early infection stages compared to the leaves at later stages indicate that they are both associated to potato defenses against V. dahliae. The expression of members of the isochorismatase family in the icsh1 mutants compensate that of ICSH1 transcripts, but this compensation disappears in presence of the potato leaf extracts. This study indicates ICSH1's involvement in V. dahliae's pathogenicity and provides more insight into its alteration of the SA/JA defense signaling's networking.

A HOPS Protein, MoVps41, Is Crucially Important for Vacuolar Morphogenesis, Vegetative Growth, Reproduction and Virulence in Magnaporthe oryzae

The homotypic fusion and protein sorting protein complex (HOPS) is the first known tether complex identified in the endocytic system that plays a key role in promoting homotypic vacuolar fusion, vacuolar biogenesis and trafficking in a wide range of organisms, including plant and fungi. However, the exact influence of the HOPS complex on growth, reproduction and pathogenicity of the economically destructive rice blast fungus has not been investigated. In this study, we identified M. oryzae vacuolar protein sorting 41 (MoVps41) an accessory subunit of HOPS complex and used targeted gene deletion approach to evaluate its contribution to growth, reproduction and infectious life cycle of the rice blast fungus. Corresponding results obtained from this study showed that MoVps41 is required for optimum vegetative development of M. oryzae and observed that MoVps41 deletion mutant displayed defective vegetative growth. Our investigation further showed that MoVps41 deletion triggered vacuolar fragmentation, compromised membrane integrity and pathogenesis of the ΔMovps41 mutant. Our studies also showed for the first time that MoVps41 plays an essential role in the regulation of sexual and asexual reproduction of M. oryzae. In summary, our study provides insight into how MoVps41 mediated vacuolar fusion and biogenesis influences reproduction, pathogenesis, and vacuolar integrity in M. oryzae and also underscores the need to holistically investigate the HOPS complex in rice blast pathogen.

Using Population and Comparative Genomics to Understand the Genetic Basis of Effector-Driven Fungal Pathogen Evolution

Epidemics caused by fungal plant pathogens pose a major threat to agro-ecosystems and impact global food security. High-throughput sequencing enabled major advances in understanding how pathogens cause disease on crops. Hundreds of fungal genomes are now available and analyzing these genomes highlighted the key role of effector genes in disease. Effectors are small secreted proteins that enhance infection by manipulating host metabolism. Fungal genomes carry 100s of putative effector genes, but the lack of homology among effector genes, even for closely related species, challenges evolutionary and functional analyses. Furthermore, effector genes are often found in rapidly evolving chromosome compartments which are difficult to assemble. We review how population and comparative genomics toolsets can be combined to address these challenges. We highlight studies that associated genome-scale polymorphisms with pathogen lifestyles and adaptation to different environments. We show how genome-wide association studies can be used to identify effectors and other pathogenicity-related genes underlying rapid adaptation. We also discuss how the compartmentalization of fungal genomes into core and accessory regions shapes the evolution of effector genes. We argue that an understanding of genome evolution provides important insight into the trajectory of host-pathogen co-evolution.

Globally distributed root endophyte Phialocephala subalpina links pathogenic and saprophytic lifestyles

BACKGROUND

Whereas an increasing number of pathogenic and mutualistic ascomycetous species were sequenced in the past decade, species showing a seemingly neutral association such as root endophytes received less attention. In the present study, the genome of Phialocephala subalpina, the most frequent species of the Phialocephala fortinii s.l. - Acephala applanata species complex, was sequenced for insight in the genome structure and gene inventory of these wide-spread root endophytes.

RESULTS

The genome of P. subalpina was sequenced using Roche/454 GS FLX technology and a whole genome shotgun strategy. The assembly resulted in 205 scaffolds and a genome size of 69.7 Mb. The expanded genome size in P. subalpina was not due to the proliferation of transposable elements or other repeats, as is the case with other ascomycetous genomes. Instead, P. subalpina revealed an expanded gene inventory that includes 20,173 gene models. Comparative genome analysis of P. subalpina with 13 ascomycetes shows that P. subalpina uses a versatile gene inventory including genes specific for pathogens and saprophytes. Moreover, the gene inventory for carbohydrate active enzymes (CAZymes) was expanded including genes involved in degradation of biopolymers, such as pectin, hemicellulose, cellulose and lignin.

CONCLUSIONS

The analysis of a globally distributed root endophyte allowed detailed insights in the gene inventory and genome organization of a yet largely neglected group of organisms. We showed that the ubiquitous root endophyte P. subalpina has a broad gene inventory that links pathogenic and saprophytic lifestyles.

Endoplasmic reticulum membrane-bound MoSec62 is involved in the suppression of rice immunity and is essential for the pathogenicity of Magnaporthe oryzae

Pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) constitutes the first line of plant inducible immunity. As an important step of plant colonization, phytopathogens have to suppress PTI, and secreted effectors are therefore co-evolved and deployed. In this study, we characterized the function of MoSec62 of Magnaporthe oryzae, the causal agent of the destructive rice blast. MoSec62 encodes a homologue of Sec62p, a yeast endoplasmic reticulum (ER) membrane transporter for precursors of secretory proteins. We showed that a T-DNA insertion into the promoter region of MoSec62, causing a disturbance to the up-regulation of MoSec62 expression during blast invasion, resulted in a complete loss of blast virulence of the mutant, M1575. Both 3,3'-diaminobenzidine (DAB) staining of the infected rice leaves and expression analysis revealed that the infectious attempt by the mutant led to strong defence responses of rice. Consistently, in transcriptomic analysis of rice leaves subject to blast inoculation, a battery of defence responses was found to be induced exclusively on M1575 challenge. For further exploration, we tested the pathogenicity on a highly susceptible rice variety and detected the accumulation of Slp1, a known PTI suppressor. Both results suggested that the mutant most likely failed to overcome rice PTI. In addition, we showed that MoSec62 was able to rescue the thermosensitivity of a yeast Δsec62, and the MoSec62-GFP fusion was co-localized to the ER membrane, both suggesting the conservation of Sec62 homologues. In conclusion, our data indicate that MoSec62, probably as an ER membrane transporter, plays an essential role in antagonizing rice defence at the early stages of blast invasion.

In planta Identification of Putative Pathogenicity Factors from the Chickpea Pathogen Ascochyta rabiei by De novo Transcriptome Sequencing Using RNA-Seq and Massive Analysis of cDNA Ends

The most important foliar diseases in legumes worldwide are ascochyta blights. Up to now, in the Ascochyta-legume pathosystem most studies focused on the identification of resistance genes in the host, while very little is known about the pathogenicity factors of the fungal pathogen. Moreover, available data were often obtained from fungi growing under artificial conditions. Therefore, in this study we aimed at the identification of the pathogenicity factors of Ascochyta rabiei, causing ascochyta blight in chickpea. To identify potential fungal pathogenicity factors, we employed RNA-seq and Massive Analysis of cDNA Ends (MACE) to produce comprehensive expression profiles of A. rabiei genes isolated either from the fungus growing in absence of its host or from fungi infecting chickpea leaves. We further provide a comprehensive de novo assembly of the A. rabiei transcriptome comprising 22,725 contigs with an average length of 1178 bp. Since pathogenicity factors are usually secreted, we predicted the A. rabiei secretome, yielding 550 putatively secreted proteins. MACE identified 596 transcripts that were up-regulated during infection. An analysis of these genes identified a collection of candidate pathogenicity factors and unraveled the pathogen's strategy for infecting its host.

Structure Analysis Uncovers a Highly Diverse but Structurally Conserved Effector Family in Phytopathogenic Fungi

Phytopathogenic ascomycete fungi possess huge effector repertoires that are dominated by hundreds of sequence-unrelated small secreted proteins. The molecular function of these effectors and the evolutionary mechanisms that generate this tremendous number of singleton genes are largely unknown. To get a deeper understanding of fungal effectors, we determined by NMR spectroscopy the 3-dimensional structures of the Magnaporthe oryzae effectors AVR1-CO39 and AVR-Pia. Despite a lack of sequence similarity, both proteins have very similar 6 β-sandwich structures that are stabilized in both cases by a disulfide bridge between 2 conserved cysteins located in similar positions of the proteins. Structural similarity searches revealed that AvrPiz-t, another effector from M. oryzae, and ToxB, an effector of the wheat tan spot pathogen Pyrenophora tritici-repentis have the same structures suggesting the existence of a family of sequence-unrelated but structurally conserved fungal effectors that we named MAX-effectors (MagnaportheAvrs and ToxB like). Structure-informed pattern searches strengthened this hypothesis by identifying MAX-effector candidates in a broad range of ascomycete phytopathogens. Strong expansion of the MAX-effector family was detected in M. oryzae and M. grisea where they seem to be particularly important since they account for 5–10% of the effector repertoire and 50% of the cloned avirulence effectors. Expression analysis indicated that the majority of M. oryzae MAX-effectors are expressed specifically during early infection suggesting important functions during biotrophic host colonization. We hypothesize that the scenario observed for MAX-effectors can serve as a paradigm for ascomycete effector diversity and that the enormous number of sequence-unrelated ascomycete effectors may in fact belong to a restricted set of structurally conserved effector families.

Fungal plant pathogens are of outstanding economic and ecological importance and cause destructive diseases on many cultivated and wild plants. Effector proteins that are secreted during infection to manipulate the host and to promote disease are a key element in fungal virulence. Phytopathogenic fungi possess huge effector repertoires that are dominated by hundreds of sequence-unrelated small secreted proteins. The molecular functions of this most important class of fungal effectors and the evolutionary mechanisms that generate this tremendous numbers of apparently unrelated proteins are largely unknown. By investigating the 3-dimensional structures of effectors from the rice blast fungus M. oryzae, we discovered an effector family comprising structurally conserved but sequence-unrelated effectors from M. oryzae and the phylogenetically distant wheat pathogen Pyrenophora tritici-repentis that we named MAX-effectors (M. oryzaeAvrs and ToxB). Structure-informed searches of whole genome sequence databases suggest that MAX-effectors are present at low frequencies and with a patchy phylogenetic distribution in many ascomycete phytopathogens. They underwent strong lineage-specific expansion in fungi of the Pyriculariae family that contains M. oryzae where they seem particularly important during biotrophic plant colonization and account for 50% of the cloned Avr effectors and 5–10% of the effector repertoire. Based on our results on the MAX-effectors and the widely accepted concept that fungal effectors evolve according to a birth-and-death model we propose the hypothesis that the majority of the immense numbers of different ascomycete effectors could in fact belong to a limited set of structurally defined families whose members are phylogenetically related.

Deciphering Genome Content and Evolutionary Relationships of Isolates from the Fungus Magnaporthe oryzae Attacking Different Host Plants

Deciphering the genetic bases of pathogen adaptation to its host is a key question in ecology and evolution. To understand how the fungus Magnaporthe oryzae adapts to different plants, we sequenced eight M. oryzae isolates differing in host specificity (rice, foxtail millet, wheat, and goosegrass), and one Magnaporthe grisea isolate specific of crabgrass. Analysis of Magnaporthe genomes revealed small variation in genome sizes (39–43 Mb) and gene content (12,283–14,781 genes) between isolates. The whole set of Magnaporthe genes comprised 14,966 shared families, 63% of which included genes present in all the nine M. oryzae genomes. The evolutionary relationships among Magnaporthe isolates were inferred using 6,878 single-copy orthologs. The resulting genealogy was mostly bifurcating among the different host-specific lineages, but was reticulate inside the rice lineage. We detected traces of introgression from a nonrice genome in the rice reference 70-15 genome. Among M. oryzae isolates and host-specific lineages, the genome composition in terms of frequencies of genes putatively involved in pathogenicity (effectors, secondary metabolism, cazome) was conserved. However, 529 shared families were found only in nonrice lineages, whereas the rice lineage possessed 86 specific families absent from the nonrice genomes. Our results confirmed that the host specificity of M. oryzae isolates was associated with a divergence between lineages without major gene flow and that, despite the strong conservation of gene families between lineages, adaptation to different hosts, especially to rice, was associated with the presence of a small number of specific gene families. All information was gathered in a public database (http://genome.jouy.inra.fr/gemo).

The Plasmodiophora brassicae genome reveals insights in its life cycle and ancestry of chitin synthases

Plasmodiophora brassicae causes clubroot, a major disease of Brassica oil and vegetable crops worldwide. P. brassicae is a Plasmodiophorid, obligate biotrophic protist in the eukaryotic kingdom of Rhizaria. Here we present the 25.5 Mb genome draft of P. brassicae, developmental stage-specific transcriptomes and a transcriptome of Spongospora subterranea, the Plasmodiophorid causing powdery scab on potato. Like other biotrophic pathogens both Plasmodiophorids are reduced in metabolic pathways. Phytohormones contribute to the gall phenotypes of infected roots. We report a protein (PbGH3) that can modify auxin and jasmonic acid. Plasmodiophorids contain chitin in cell walls of the resilient resting spores. If recognized, chitin can trigger defense responses in plants. Interestingly, chitin-related enzymes of Plasmodiophorids built specific families and the carbohydrate/chitin binding (CBM18) domain is enriched in the Plasmodiophorid secretome. Plasmodiophorids chitin synthases belong to two families, which were present before the split of the eukaryotic Stramenopiles/Alveolates/Rhizaria/Plantae and Metazoa/Fungi/Amoebozoa megagroups, suggesting chitin synthesis to be an ancient feature of eukaryotes. This exemplifies the importance of genomic data from unexplored eukaryotic groups, such as the Plasmodiophorids, to decipher evolutionary relationships and gene diversification of early eukaryotes.

Horizontal transfer of carbohydrate metabolism genes into ectomycorrhizal Amanita

The genus Amanita encompasses both symbiotic, ectomycorrhizal fungi and asymbiotic litter decomposers; all species are derived from asymbiotic ancestors. Symbiotic species are no longer able to degrade plant cell walls. The carbohydrate esterases family 1 (CE1s) is a diverse group of enzymes involved in carbon metabolism, including decomposition and carbon storage. CE1 genes of the ectomycorrhizal A. muscaria appear diverged from all other fungal homologues, and more similar to CE1s of bacteria, suggesting a horizontal gene transfer (HGT) event. In order to test whether AmanitaCE1s were acquired horizontally, we built a phylogeny of CE1s collected from across the tree of life, and describe the evolution of CE1 genes among Amanita and relevant lineages of bacteria. CE1s of symbiotic Amanita were very different from CE1s of asymbiotic Amanita, and are more similar to bacterial CE1s. The protein structure of one CE1 gene of A. muscaria matched a depolymerase that degrades the carbon storage molecule poly((R)-3-hydroxybutyrate) (PHB). Asymbiotic Amanita do not carry sequence or structural homologues of these genes. The CE1s acquired through HGT may enable novel metabolisms, or play roles in signaling or defense. This is the first evidence for the horizontal transfer of carbohydrate metabolism genes into ectomycorrhizal fungi.