Genome-Wide Analysis in Three Fusarium Pathogens Identifies Rapidly Evolving Chromosomes and Genes Associated with Pathogenicity

New insights into mycotoxin risk management through fungal population genetics and genomics

Mycotoxin contamination of food and feed is a major global concern. Chronic or acute dietary exposure to contaminated food and feed can negatively affect both human and animal health. Contamination occurs through plant infection by toxigenic fungi, primarily Aspergillus and Fusarium spp., either before or after harvest. Despite the application of various management strategies, controlling these pathogens remains a major challenge primarily because of their ability to adapt to environmental changes and selection pressures. Understanding the genetic structure of plant pathogen populations is pivotal for gaining new insights into their biology and epidemiology, as well as for understanding the mechanisms behind their adaptability. Such deeper understanding is crucial for developing effective and preemptive management strategies tailored to the evolving nature of pathogenic populations. This review focuses on the population-level variations within the two most economically significant toxigenic fungal genera according to space, host, and pathogenicity. Outcomes in terms of migration patterns, gene flow within populations, mating abilities, and the potential for host jumps are examined. We also discuss effective yet often underutilized applications of population genetics and genomics to address practical challenges in the epidemiology and disease control of toxigenic fungi.

A roadmap to understanding diversity and function of coral reef-associated fungi

Tropical coral reefs are hotspots of marine productivity, owing to the association of reef-building corals with endosymbiotic algae and metabolically diverse bacterial communities. However, the functional importance of fungi, well-known for their contribution to shaping terrestrial ecosystems and global nutrient cycles, remains underexplored on coral reefs. We here conceptualize how fungal functional traits may have facilitated the spread, diversification, and ecological adaptation of marine fungi on coral reefs. We propose that functions of reef-associated fungi may be diverse and go beyond their hitherto described roles of pathogens and bioeroders, including but not limited to reef-scale biogeochemical cycles and the structuring of coral-associated and environmental microbiomes via chemical mediation. Recent technological and conceptual advances will allow the elucidation of the physiological, ecological, and chemical contributions of understudied marine fungi to coral holobiont and reef ecosystem functioning and health and may help provide an outlook for reef management actions.

Profile of the in silico secretome of the palm dieback pathogen, Fusarium oxysporum f. sp. albedinis, a fungus that puts natural oases at risk

Understanding biotic changes that occur alongside climate change constitute a research priority of global significance. Here, we address a plant pathogen that poses a serious threat to life on natural oases, where climate change is already taking a toll and severely impacting human subsistence. Fusarium oxysporum f. sp. albedinis is a pathogen that causes dieback disease on date palms, a tree that provides several critical ecosystem services in natural oases; and consequently, of major importance in this vulnerable habitat. Here, we assess the current state of global pathogen spread, we annotate the genome of a sequenced pathogen strain isolated from the native range and we analyse its in silico secretome. The palm dieback pathogen secretes a large arsenal of effector candidates including a variety of toxins, a distinguished profile of secreted in xylem proteins (SIX) as well as an expanded protein family with an N-terminal conserved motif [SG]PC[KR]P that could be involved in interactions with host membranes. Using agrobiodiversity as a strategy to decrease pathogen infectivity, while providing short term resilient solutions, seems to be widely overcome by the pathogen. Hence, the urgent need for future mechanistic research on the palm dieback disease and a better understanding of pathogen genetic diversity.

Transcriptome Analysis of Fusarium Root-Rot-Resistant and -Susceptible Alfalfa (Medicago sativa L.) Plants during Plant–Pathogen Interactions

Alfalfa (Medicago sativa L.) is a perennial leguminous forage cultivated globally. Fusarium spp.-induced root rot is a chronic and devastating disease affecting alfalfa that occurs in most production fields. Studying the disease resistance regulatory network and investigating the key genes involved in plant–pathogen resistance can provide vital information for breeding alfalfa that are resistant to Fusarium spp. In this study, a resistant and susceptible clonal line of alfalfa was inoculated with Fusarium proliferatum L1 and sampled at 24 h, 48 h, 72 h, and 7 d post-inoculation for RNA-seq analysis. Among the differentially expressed genes (DEGs) detected between the two clonal lines at the four time points after inoculation, approximately 81.8% were detected at 24 h and 7 d after inoculation. Many DEGs in the two inoculated clonal lines participated in PAMP-triggered immunity (PTI) and effector-triggered immunity (ETI) mechanisms. In addition, transcription factor families such as bHLH, SBP, AP2, WRKY, and MYB were detected in response to infection. These results are an important supplement to the few existing studies on the resistance regulatory network of alfalfa against Fusarium root rot and will help to understand the evolution of host–pathogen interactions.

Identification and Characterization of Fusarium nirenbergiae Associated with Saffron Corm Rot Disease

Saffron (Crocus sativus L.) is the most expensive spice plant and is distributed widely around the world. However, its production is limited by corm rot, a disastrous disease, attributed to Fusarium oxysporum in many regions of the world. In 2020, extensive surveys were carried out in Zhejiang, Shanghai, Anhui, and Guizhou provinces as saffron growing areas of China. Fourteen single-spore isolates were obtained from rotted corms and identified as F. nirenbergiae according to morphological appearance and multilocus phylogenetic analysis with translation elongation factor 1-α (tef1), DNA-directed RNA polymerase II largest subunit (rpb2), and β-tubulin (tub2). Results of the pathogenicity assay supported the conclusion that F. nirenbergiae is the pathogen responsible for corm rot. In this study, we obtained the whole genome sequence of two highly virulent F. nirenbergiae strains via the Illumina HiSeq platform. Genome sequence assemblies of approximately 52.7 and 52.2 Mb were generated for isolates WY5 and SH1, respectively. To the best of our knowledge, this is the first report of F. nirenbergiae causing C. sativus corm rot in China and indeed worldwide. Results from this research contribute to our understanding of genetic diversity, genomic information, and host determination, which will enable researchers to design appropriate management measures for this hazardous disease.

Defence response in plants and animals against a common fungal pathogen, Fusarium oxysporum

•

Fusarium oxysporum species complex (FOSC) is considered one of the most devastating plant pathogen.

•

FOSC is an emerging pathogen of immunocompromised individuals.

•

Mycotoxins produced by FOSC predisposes the host to other pathogens.

•

Comparative immune reactions in plant and invertebrate show that several antimicrobial peptides (AMPs) and secondary metabolites maybe used as control against Fusarium infection.

Plant pathogens emerging as threat to human and animal health has been a matter of concern within the scientific community. Fusarium oxysporum, predominantly a phytopathogen, can infect both plants and animals. As a plant pathogen, F. oxysporum is one of the most economically damaging pathogen. In humans, F. oxysporum can infect immunocompromised individuals and is increasingly being considered as a problematic pathogen. Mycotoxins produced by F. oxysporum supress the innate immune pathways in both plants and animals. Hence, F. oxysporum is the perfect example for studying similarities and differences between defence strategies adopted by plants and animals. In this review we will discuss the innate immune response of plant and animal hosts for protecting against F. oxysporum infection. Such studies will be helpful for identifying genes, protein and metabolites with antifungal properties suitable for protecting humans.

A single amino acid polymorphism in a conserved effector of the multihost blast fungus pathogen expands host-target binding spectrum

Accelerated gene evolution is a hallmark of pathogen adaptation and specialization following host-jumps. However, the molecular processes associated with adaptive evolution between host-specific lineages of a multihost plant pathogen remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), host specialization on different grass hosts is generally associated with dynamic patterns of gain and loss of virulence effector genes that tend to define the distinct genetic lineages of this pathogen. Here, we unravelled the biochemical and structural basis of adaptive evolution of APikL2, an exceptionally conserved paralog of the well-studied rice-lineage specific effector AVR-Pik. Whereas AVR-Pik and other members of the six-gene AVR-Pik family show specific patterns of presence/absence polymorphisms between grass-specific lineages of M. oryzae, APikL2 stands out by being ubiquitously present in all blast fungus lineages from 13 different host species. Using biochemical, biophysical and structural biology methods, we show that a single aspartate to asparagine polymorphism expands the binding spectrum of APikL2 to host proteins of the heavy-metal associated (HMA) domain family. This mutation maps to one of the APikL2-HMA binding interfaces and contributes to an altered hydrogen-bonding network. By combining phylogenetic ancestral reconstruction with an analysis of the structural consequences of allelic diversification, we revealed a common mechanism of effector specialization in the AVR-Pik/APikL2 family that involves two major HMA-binding interfaces. Together, our findings provide a detailed molecular evolution and structural biology framework for diversification and adaptation of a fungal pathogen effector family following host-jumps.

Genome-wide association study for deoxynivalenol production and aggressiveness in wheat and rye head blight by resequencing 92 isolates of Fusarium culmorum

Background

Fusarium culmorum is an important pathogen causing head blight of cereals in Europe. This disease is of worldwide importance leading to reduced yield, grain quality, and contamination by mycotoxins. These mycotoxins are harmful for livestock and humans; therefore, many countries have strict regulatory limits for raw materials and processed food. Extensive genetic diversity is described among field populations of F. culmorum isolates for aggressiveness and production of the trichothecene mycotoxin deoxynivalenol (DON). However, the causes for this quantitative variation are not clear, yet. We analyzed 92 isolates sampled from different field populations in Germany, Russia, and Syria together with an international collection for aggressiveness and DON production in replicated field experiments at two locations in two years with two hosts, wheat and rye. The 30x coverage whole-genome resequencing of all isolates resulted in the identification of 130,389 high quality single nucleotide polymorphisms (SNPs) that were used for the first genome-wide association study in this phytopathogenic fungus.

Results

In wheat, 20 and 27 SNPs were detected for aggressiveness and DON content, respectively, of which 10 overlapped. Additionally, two different SNPs were significantly associated with aggressiveness in rye that were among those SNPs being associated with DON production in wheat. Most of the SNPs explained only a small proportion of genotypic variance (p_G), however, four SNPs were associated with major quantitative trait loci (QTLs) with p_G ranging from 12 to 48%. The QTL with the highest p_G was involved in DON production and associated with a SNP most probably located within the Tri4 gene.

Conclusions

The diversity of 92 isolates of F. culmorum were captured using a heuristic approach. Key phenotypic traits, SNPs, and candidate genes underlying aggressiveness and DON production were identified. Clearly, many QTLs are responsible for aggressiveness and DON content in wheat, both traits following a quantitative inheritance. Several SNPs involved in DON metabolism, among them the Tri4 gene of the trichothecene pathway, were inferred as important source of variation in fungal aggressiveness. Using this information underlying the phenotypic variation will be of paramount importance in evaluating strategies for successful resistance breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07931-5.

Rapid sequence evolution driven by transposable elements at a virulence locus in a fungal wheat pathogen

BACKGROUND

Plant pathogens cause substantial crop losses in agriculture production and threaten food security. Plants evolved the ability to recognize virulence factors and pathogens have repeatedly escaped recognition due rapid evolutionary change at pathogen virulence loci (i.e. effector genes). The presence of transposable elements (TEs) in close physical proximity of effector genes can have important consequences for gene regulation and sequence evolution. Species-wide investigations of effector gene loci remain rare hindering our ability to predict pathogen evolvability.

RESULTS

Here, we performed genome-wide association studies (GWAS) on a highly polymorphic mapping population of 120 isolates of Zymoseptoria tritici, the most damaging pathogen of wheat in Europe. We identified a major locus underlying significant variation in reproductive success of the pathogen and damage caused on the wheat cultivar Claro. The most strongly associated locus is intergenic and flanked by genes encoding a predicted effector and a serine-type endopeptidase. The center of the locus contained a highly dynamic region consisting of multiple families of TEs. Based on a large global collection of assembled genomes, we show that the virulence locus has undergone substantial recent sequence evolution. Large insertion and deletion events generated length variation between the flanking genes by a factor of seven (5-35 kb). The locus showed also strong signatures of genomic defenses against TEs (i.e. RIP) contributing to the rapid diversification of the locus.

CONCLUSIONS

In conjunction, our work highlights the power of combining GWAS and population-scale genome analyses to investigate major effect loci in pathogens.

Whole-Genome Comparisons of Ergot Fungi Reveals the Divergence and Evolution of Species within the Genus Claviceps Are the Result of Varying Mechanisms Driving Genome Evolution and Host Range Expansion

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

Genome Sequence of Fusarium oxysporum f. sp. conglutinans, the Etiological Agent of Cabbage Fusarium Wilt

Fusarium oxysporum f. sp. conglutinans is the causal agent of Fusarium wilt of cabbage (Brassica oleracea var. capitata L.), which results in severe yield loss. Here, we report a high-quality genome sequence of a race 1 strain (IVC-1) of F. oxysporum f. sp. conglutinans, which was assembled using a combination of PacBio long-read and Illumina short-read sequences. The assembled IVC-1 genome has a total size of 71.18 Mb, with a contig N50 length of 4.59 Mb, and encodes 23,374 predicted protein-coding genes. The high-quality genome of IVC-1 provides a valuable resource for facilitating our understanding of F. oxysporum f. sp. conglutinans-cabbage interaction.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Divergence and Gene Flow Between Fusarium subglutinans and F. temperatum Isolated from Maize in Argentina

Fusarium subglutinans and F. temperatum are two important fungal pathogens of maize whose distinctness as separate species has been difficult to assess. We isolated strains of these species from commercial and native maize varieties in Argentina and sequenced >28,000 loci to estimate genetic variation in the sample. Our objectives were to measure genetic divergence between the species, infer demographic parameters related to their split, and describe the population structure of the sample. When analyzed together, over 30% of each species' polymorphic sites (>2,500 sites) segregate as polymorphisms in the other. Demographic modeling confirmed the species split predated maize domestication, but subsequent between-species gene flow has occurred, with gene flow from F. subglutinans into F. temperatum greater than gene flow in the reverse direction. In F. subglutinans, little evidence exists for substructure or recent selective sweeps, but there is evidence for limited sexual reproduction. In F. temperatum, there is clear evidence for population substructure and signals of abundant recent selective sweeps, with sexual reproduction probably less common than in F. subglutinans. Both genetic variation and the relative number of polymorphisms shared between species increase near the telomeres of all 12 chromosomes, where genes related to plant-pathogen interactions often are located. Our results suggest that species boundaries between closely related Fusarium species can be semipermeable and merit further study. Such semipermeability could facilitate unanticipated genetic exchange between species and enable quicker permanent responses to changes in the agro-ecosystem, e.g., pathogen-resistant host varieties, new chemical and biological control agents, and agronomic practices.

Whole genome analysis of the koa wilt pathogen (Fusarium oxysporum f. sp. koae) and the development of molecular tools for early detection and monitoring

BACKGROUND

Development and application of DNA-based methods to distinguish highly virulent isolates of Fusarium oxysporum f. sp. koae [Fo koae; cause of koa wilt disease on Acacia koa (koa)] will help disease management through early detection, enhanced monitoring, and improved disease resistance-breeding programs.

RESULTS

This study presents whole genome analyses of one highly virulent Fo koae isolate and one non-pathogenic F. oxysporum (Fo) isolate. These analyses allowed for the identification of putative lineage-specific DNA and predicted genes necessary for disease development on koa. Using putative chromosomes and predicted gene comparisons, Fo koae-exclusive, virulence genes were identified. The putative lineage-specific DNA included identified genes encoding products secreted in xylem (e. g., SIX1 and SIX6) that may be necessary for disease development on koa. Unique genes from Fo koae were used to develop pathogen-specific PCR primers. These diagnostic primers allowed target amplification in the characterized highly virulent Fo koae isolates but did not allow product amplification in low-virulence or non-pathogenic isolates of Fo. Thus, primers developed in this study will be useful for early detection and monitoring of highly virulent strains of Fo koae. Isolate verification is also important for disease resistance-breeding programs that require a diverse set of highly virulent Fo koae isolates for their disease-screening assays to develop disease-resistant koa.

CONCLUSIONS

These results provide the framework for understanding the pathogen genes necessary for koa wilt disease and the genetic variation of Fo koae populations across the Hawaiian Islands.

De novo sequencing, assembly and functional annotation of Armillaria borealis genome

BACKGROUND

Massive forest decline has been observed almost everywhere as a result of negative anthropogenic and climatic effects, which can interact with pests, fungi and other phytopathogens and aggravate their effects. Climatic changes can weaken trees and make fungi, such as Armillaria more destructive. Armillaria borealis (Marxm. & Korhonen) is a fungus from the Physalacriaceae family (Basidiomycota) widely distributed in Eurasia, including Siberia and the Far East. Species from this genus cause the root white rot disease that weakens and often kills woody plants. However, little is known about ecological behavior and genetics of A. borealis. According to field research data, A. borealis is less pathogenic than A. ostoyae, and its aggressive behavior is quite rare. Mainly A. borealis behaves as a secondary pathogen killing trees already weakened by other factors. However, changing environment might cause unpredictable effects in fungus behavior.

RESULTS

The de novo genome assembly and annotation were performed for the A. borealis species for the first time and presented in this study. The A. borealis genome assembly contained ~ 68 Mbp and was comparable with ~ 60 and ~ 79.5 Mbp for the A. ostoyae and A. mellea genomes, respectively. The N50 for contigs equaled 50,544 bp. Functional annotation analysis revealed 21,969 protein coding genes and provided data for further comparative analysis. Repetitive sequences were also identified. The main focus for further study and comparative analysis will be on the enzymes and regulatory factors associated with pathogenicity.

CONCLUSIONS

Pathogenic fungi such as Armillaria are currently one of the main problems in forest conservation. A comprehensive study of these species and their pathogenicity is of great importance and needs good genomic resources. The assembled genome of A. borealis presented in this study is of sufficiently good quality for further detailed comparative study on the composition of enzymes in other Armillaria species. There is also a fundamental problem with the identification and classification of species of the Armillaria genus, where the study of repetitive sequences in the genomes of basidiomycetes and their comparative analysis will help us identify more accurately taxonomy of these species and reveal their evolutionary relationships.

Pathogenomics and Management of Fusarium Diseases in Plants

There is an urgency to supplant the heavy reliance on chemical control of Fusarium diseases in different economically important, staple food crops due to development of resistance in the pathogen population, the high cost of production to the risk-averse grower, and the concomitant environmental impacts. Pathogenomics has enabled (i) the creation of genetic inventories which identify those putative genes, regulators, and effectors that are associated with virulence, pathogenicity, and primary and secondary metabolism; (ii) comparison of such genes among related pathogens; (iii) identification of potential genetic targets for chemical control; and (iv) better characterization of the complex dynamics of host–microbe interactions that lead to disease. This type of genomic data serves to inform host-induced gene silencing (HIGS) technology for targeted disruption of transcription of select genes for the control of Fusarium diseases. This review discusses the various repositories and browser access points for comparison of genomic data, the strategies for identification and selection of pathogenicity- and virulence-associated genes and effectors in different Fusarium species, HIGS and successful Fusarium disease control trials with a consideration of loss of RNAi, off-target effects, and future challenges in applying HIGS for management of Fusarium diseases.

Three-Locus Sequence Identification and Differential Tebuconazole Sensitivity Suggest Novel Fusarium equiseti Haplotype from Trinidad

The Fusarium incarnatum-equiseti species complex (FIESC) consists of 33 phylogenetic species according to multi-locus sequence typing (MLST) and Genealogical Concordance Phylogenetic Species Recognition (GCPSR). A multi-locus dataset consisting of nucleotide sequences of the translation elongation factor (EF-1α), calmodulin (CAM), partial RNA polymerase largest subunit (RPB1), and partial RNA polymerase second largest subunit (RPB2), was generated to distinguish among phylogenetic species within the FIESC isolates infecting bell pepper in Trinidad. Three phylogenetic species belonged to the Incarnatum clade (FIESC-15, FIESC-16, and FIESC-26), and one species belonged to the Equiseti clade (FIESC-14). Specific MLST types were sensitive to 10 µg/mL of tebuconazole fungicide as a discriminatory dose. The EC50 values were significantly different among the four MLST groups, which were separated into two homogeneous groups: FIESC-26a and FIESC-14a, demonstrating the “sensitive” azole phenotype and FIESC-15a and FIESC-16a as the “less sensitive” azole phenotype. CYP51C sequences of the Trinidad isolates, although under positive selection, were without any signatures of recombination, were highly conserved, and were not correlated with these azole phenotypes. CYP51C sequences were unable to resolve the FIESC isolates as phylogenetic inference indicated polytomic branching for these sequences. This data is important to different research communities, including those studying Fusarium phytopathology, mycotoxins, and public health impacts.

Pathogenicity and Virulence Factors of Fusarium graminearum Including Factors Discovered Using Next Generation Sequencing Technologies and Proteomics

Fusarium graminearum is a devasting mycotoxin-producing pathogen of grain crops. F. graminearum has been extensively studied to understand its pathogenicity and virulence factors. These studies gained momentum with the advent of next-generation sequencing (NGS) technologies and proteomics. NGS and proteomics have enabled the discovery of a multitude of pathogenicity and virulence factors of F. graminearum. This current review aimed to trace progress made in discovering F. graminearum pathogenicity and virulence factors in general, as well as pathogenicity and virulence factors discovered using NGS, and to some extent, using proteomics. We present more than 100 discovered pathogenicity or virulence factors and conclude that although a multitude of pathogenicity and virulence factors have already been discovered, more work needs to be done to take advantage of NGS and its companion applications of proteomics.

A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici

BACKGROUND

The gene content of a species largely governs its ecological interactions and adaptive potential. A species is therefore defined by both core genes shared between all individuals and accessory genes segregating presence-absence variation. There is growing evidence that eukaryotes, similar to bacteria, show intra-specific variability in gene content. However, it remains largely unknown how functionally relevant such a pangenome structure is for eukaryotes and what mechanisms underlie the emergence of highly polymorphic genome structures.

RESULTS

Here, we establish a reference-quality pangenome of a fungal pathogen of wheat based on 19 complete genomes from isolates sampled across six continents. Zymoseptoria tritici causes substantial worldwide losses to wheat production due to rapidly evolved tolerance to fungicides and evasion of host resistance. We performed transcriptome-assisted annotations of each genome to construct a global pangenome. Major chromosomal rearrangements are segregating within the species and underlie extensive gene presence-absence variation. Conserved orthogroups account for only ~ 60% of the species pangenome. Investigating gene functions, we find that the accessory genome is enriched for pathogenesis-related functions and encodes genes involved in metabolite production, host tissue degradation and manipulation of the immune system. De novo transposon annotation of the 19 complete genomes shows that the highly diverse chromosomal structure is tightly associated with transposable element content. Furthermore, transposable element expansions likely underlie recent genome expansions within the species.

CONCLUSIONS

Taken together, our work establishes a highly complex eukaryotic pangenome providing an unprecedented toolbox to study how pangenome structure impacts crop-pathogen interactions.

A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici

BACKGROUND

The gene content of a species largely governs its ecological interactions and adaptive potential. A species is therefore defined by both core genes shared between all individuals and accessory genes segregating presence-absence variation. There is growing evidence that eukaryotes, similar to bacteria, show intra-specific variability in gene content. However, it remains largely unknown how functionally relevant such a pangenome structure is for eukaryotes and what mechanisms underlie the emergence of highly polymorphic genome structures.

RESULTS

Here, we establish a reference-quality pangenome of a fungal pathogen of wheat based on 19 complete genomes from isolates sampled across six continents. Zymoseptoria tritici causes substantial worldwide losses to wheat production due to rapidly evolved tolerance to fungicides and evasion of host resistance. We performed transcriptome-assisted annotations of each genome to construct a global pangenome. Major chromosomal rearrangements are segregating within the species and underlie extensive gene presence-absence variation. Conserved orthogroups account for only ~ 60% of the species pangenome. Investigating gene functions, we find that the accessory genome is enriched for pathogenesis-related functions and encodes genes involved in metabolite production, host tissue degradation and manipulation of the immune system. De novo transposon annotation of the 19 complete genomes shows that the highly diverse chromosomal structure is tightly associated with transposable element content. Furthermore, transposable element expansions likely underlie recent genome expansions within the species.

CONCLUSIONS

Taken together, our work establishes a highly complex eukaryotic pangenome providing an unprecedented toolbox to study how pangenome structure impacts crop-pathogen interactions.

Dynamic virulence-related regions of the plant pathogenic fungus Verticillium dahliae display enhanced sequence conservation

Plant pathogens continuously evolve to evade host immune responses. During host colonization, many fungal pathogens secrete effectors to perturb such responses, but these in turn may become recognized by host immune receptors. To facilitate the evolution of effector repertoires, such as the elimination of recognized effectors, effector genes often reside in genomic regions that display increased plasticity, a phenomenon that is captured in the two-speed genome hypothesis. The genome of the vascular wilt fungus Verticillium dahliae displays regions with extensive presence/absence polymorphisms, so-called lineage-specific regions, that are enriched in in planta-induced putative effector genes. As expected, comparative genomics reveals differential degrees of sequence divergence between lineage-specific regions and the core genome. Unanticipated, lineage-specific regions display markedly higher sequence conservation in coding as well as noncoding regions than the core genome. We provide evidence that disqualifies horizontal transfer to explain the observed sequence conservation and conclude that sequence divergence occurs at a slower pace in lineage-specific regions of the V. dahliae genome. We hypothesize that differences in chromatin organisation may explain lower nucleotide substitution rates in the plastic, lineage-specific regions of V. dahliae.

Local and systemic regulation of PSII efficiency in triticale infected by the hemibiotrophic pathogen Microdochium nivale

Microdochium nivale is a fungal pathogen that causes yield losses of cereals during winter. Cold hardening under light conditions induces genotype-dependent resistance of a plant to infection. We aim to show how photosystem II (PSII) regulation contributes to plant resistance. Using mapping population of triticale doubled haploid lines, three M. nivale strains and different infection assays, we demonstrate that plants that maintain a higher maximum quantum efficiency of PSII show less leaf damage upon infection. The fungus can establish necrotrophic or biotrophic interactions with susceptible or resistant genotypes, respectively. It is suggested that local inhibition of photosynthesis during the infection of sensitive genotypes is not balanced by a supply of energy from the tissue surrounding the infected cells as efficiently as in resistant genotypes. Thus, defence is limited, which in turn results in extensive necrotic damage. Quantitative trait loci regions, involved in the control of both PSII functioning and resistance, were located on chromosomes 4 and 6, similar to a wide range of PSII- and resistance-related genes. A meta-analysis of microarray experiments showed that the expression of genes involved in the repair and de novo assembly of PSII was maintained at a stable level. However, to establish a favourable energy balance for defence, genes encoding PSII proteins resistant to oxidative degradation were downregulated to compensate for the upregulation of defence-related pathways. Finally, we demonstrate that the structural and functional integrity of the plant is a factor required to meet the energy demand of infected cells, photosynthesis-dependent systemic signalling and defence responses.

Transposable elements contribute to fungal genes and impact fungal lifestyle

The last decade brought a still growing experimental evidence of mobilome impact on host's gene expression. We systematically analysed genomic location of transposable elements (TEs) in 625 publicly available fungal genomes from the NCBI database in order to explore their potential roles in genome evolution and correlation with species' lifestyle. We found that non-autonomous TEs and remnant copies are evenly distributed across genomes. In consequence, they also massively overlap with regions annotated as genes, which suggests a great contribution of TE-derived sequences to host's coding genome. Younger and potentially active TEs cluster with one another away from genic regions. This non-randomness is a sign of either selection against insertion of TEs in gene proximity or target site preference among some types of TEs. Proteins encoded by genes with old transposable elements insertions have significantly less repeat and protein-protein interaction motifs but are richer in enzymatic domains. However, genes only proximal to TEs do not display any functional enrichment. Our findings show that adaptive cases of TE insertion remain a marginal phenomenon, and the overwhelming majority of TEs are evolving neutrally. Eventually, animal-related and pathogenic fungi have more TEs inserted into genes than fungi with other lifestyles. This is the first systematic, kingdom-wide study concerning mobile elements and their genomic neighbourhood. The obtained results should inspire further research concerning the roles TEs played in evolution and how they shape the life we know today.

Arms race: diverse effector proteins with conserved motifs

Effector proteins play important roles in the infection by pathogenic oomycetes and fungi or the colonization by endophytic and mycorrhizal fungi. They are either translocated into the host plant cells via specific translocation mechanisms and function in the host's cytoplasm or nucleus, or they reside in the apoplast of the plant cells and act at the extracellular host-microbe interface. Many effector proteins possess conserved motifs (such as the RXLR, CRN, LysM, RGD, DELD, EAR, RYWT, Y/F/WXC or CFEM motifs) localized in their N- or C-terminal regions. Analysis of the functions of effector proteins, especially so-called "core effectors", is crucial for the understanding of pathogenicity/symbiosis mechanisms and plant defense strategies, and helps to develop breeding strategies for pathogen-resistant cultivars, and to increase crop yield and quality as well as abiotic stress resistance. This review summarizes current knowledge about these effector proteins with the conversed motifs and their involvement in pathogenic or mutualistic plant/fungal interactions.

Multiple independent origins for a subtelomeric locus associated with growth rate in Fusarium circinatum

Fusarium is a diverse assemblage that includes a large number of species of considerable medical and agricultural importance. Not surprisingly, whole genome sequences for many Fusarium species have been published or are in the process of being determined, the availability of which is invaluable for deciphering the genetic basis of key phenotypic traits. Here we investigated the distribution, genic composition, and evolutionary history of a locus potentially determining growth rate in the pitch canker pathogen F. circinatum. We found that the genomic region underlying this locus is highly conserved amongst F. circinatum and its close relatives, except for the presence of a 12 000 base pair insertion in all of the examined isolates of F. circinatum. This insertion encodes for five genes and our phylogenetic analyses revealed that each was most likely acquired through horizontal gene transfer from polyphyletic origins. Our data further showed that this region is located in a region low in G+C content and enriched for repetitive sequences and transposable elements, which is situated near the telomere of Chromosome 3 of F. circinatum. As have been shown for other fungi, these findings thus suggest that the emergence of the unique 12 000 bp region in F. circinatum is linked to the dynamic evolutionary processes associated with subtelomeres that, in turn, have been implicated in the ecological adaptation of fungal pathogens.

Inter-genome comparison of the Quorn fungus Fusarium venenatum and the closely related plant infecting pathogen Fusarium graminearum

BACKGROUND

The soil dwelling saprotrophic non-pathogenic fungus Fusarium venenatum, routinely used in the commercial fermentation industry, is phylogenetically closely related to the globally important cereal and non-cereal infecting pathogen F. graminearum. This study aimed to sequence, assemble and annotate the F. venenatum (strain A3/5) genome, and compare this genome with F. graminearum.

RESULTS

Using shotgun sequencing, a 38,660,329 bp F. venenatum genome was assembled into four chromosomes, and a 78,618 bp mitochondrial genome. In comparison to F. graminearum, the predicted gene count of 13,946 was slightly lower. The F. venenatum centromeres were found to be 25% smaller compared to F. graminearum. Chromosome length was 2.8% greater in F. venenatum, primarily due to an increased abundance of repetitive elements and transposons, but not transposon diversity. On chromosome 3 a major sequence rearrangement was found, but its overall gene content was relatively unchanged. Unlike homothallic F. graminearum, heterothallic F. venenatum possessed the MAT1-1 type locus, but lacked the MAT1-2 locus. The F. venenatum genome has the type A trichothecene mycotoxin TRI5 cluster, whereas F. graminearum has type B. From the F. venenatum gene set, 786 predicted proteins were species-specific versus NCBI. The annotated F. venenatum genome was predicted to possess more genes coding for hydrolytic enzymes and species-specific genes involved in the breakdown of polysaccharides than F. graminearum. Comparison of the two genomes reduced the previously defined F. graminearum-specific gene set from 741 to 692 genes. A comparison of the F. graminearum versus F. venenatum proteomes identified 15 putative secondary metabolite gene clusters (SMC), 109 secreted proteins and 38 candidate effectors not found in F. venenatum. Five of the 15 F. graminearum-specific SMCs that were either absent or highly divergent in the F. venenatum genome showed increased in planta expression. In addition, two predicted F. graminearum transcription factors previously shown to be required for fungal virulence on wheat plants were absent or exhibited high sequence divergence.

CONCLUSIONS

This study identifies differences between the F. venenatum and F. graminearum genomes that may contribute to contrasting lifestyles, and highlights the repertoire of F. graminearum-specific candidate genes and SMCs potentially required for pathogenesis.

Fungal Genomes and Insights into the Evolution of the Kingdom

The kingdom Fungi comprises species that inhabit nearly all ecosystems. Fungi exist as both free-living and symbiotic unicellular and multicellular organisms with diverse morphologies. The genomes of fungi encode genes that enable them to thrive in diverse environments, invade plant and animal cells, and participate in nutrient cycling in terrestrial and aquatic ecosystems. The continuously expanding databases of fungal genome sequences have been generated by individual and large-scale efforts such as Génolevures, Broad Institute's Fungal Genome Initiative, and the 1000 Fungal Genomes Project (http://1000.fungalgenomes.org). These efforts have produced a catalog of fungal genes and genomic organization. The genomic datasets can be utilized to better understand how fungi have adapted to their lifestyles and ecological niches. Large datasets of fungal genomic and transcriptomic data have enabled the use of novel methodologies and improved the study of fungal evolution from a molecular sequence perspective. Combined with microscopes, petri dishes, and woodland forays, genome sequencing supports bioinformatics and comparative genomics approaches as important tools in the study of the biology and evolution of fungi.

FcRav2, a gene with a ROGDI domain involved in Fusarium head blight and crown rot on durum wheat caused by Fusarium culmorum

Fusarium culmorum is a soil-borne fungal pathogen which causes foot and root rot and Fusarium head blight on small-grain cereals, in particular wheat and barley. It causes significant yield and quality losses and results in the contamination of kernels with type B trichothecene mycotoxins. Our knowledge of the pathogenicity factors of this fungus is still limited. A transposon tagging approach based on the mimp1/impala double-component system has allowed us to select a mutant altered in multiple metabolic and morphological processes, trichothecene production and virulence. The flanking regions of mimp1 were used to seek homologies in the F. culmorum genome, and revealed that mimp1 had reinserted within the last exon of a gene encoding a hypothetical protein of 318 amino acids which contains a ROGDI-like leucine zipper domain, supposedly playing a protein-protein interaction or regulatory role. By functional complementation and bioinformatic analysis, we characterized the gene as the yeast Rav2 homologue, confirming the high level of divergence in multicellular fungi. Deletion of FcRav2 or its orthologous gene in F. graminearum highlighted its ability to influence a number of functions, including virulence, trichothecene type B biosynthesis, resistance to azoles and resistance to osmotic and oxidative stress. Our results indicate that the FcRav2 protein (and possibly the RAVE complex as a whole) may become a suitable target for new antifungal drug development or the plant-mediated resistance response in filamentous fungi of agricultural interest.

Fusarium verticillioides: Advancements in Understanding the Toxicity, Virulence, and Niche Adaptations of a Model Mycotoxigenic Pathogen of Maize

The importance of understanding the biology of the mycotoxigenic fungus Fusarium verticillioides and its various microbial and plant host interactions is critical given its threat to maize, one of the world's most valuable food crops. Disease outbreaks and mycotoxin contamination of grain threaten economic returns and have grave implications for human and animal health and food security. Furthermore, F. verticillioides is a member of a genus of significant phytopathogens and, thus, data regarding its host association, biosynthesis of secondary metabolites, and other metabolic (degradative) capabilities are consequential to both basic and applied research efforts across multiple pathosystems. Notorious among its secondary metabolites are the fumonisin mycotoxins, which cause severe animal diseases and are implicated in human disease. Additionally, studies of these mycotoxins have led to new understandings of F. verticillioides plant pathogenicity and provide tools for research into cellular processes and host-pathogen interaction strategies. This review presents current knowledge regarding several significant lines of F. verticillioides research, including facets of toxin production, virulence, and novel fitness strategies exhibited by this fungus across rhizosphere and plant environments.

Positively Selected Effector Genes and Their Contribution to Virulence in the Smut Fungus Sporisorium reilianum

Plants and fungi display a broad range of interactions in natural and agricultural ecosystems ranging from symbiosis to parasitism. These ecological interactions result in coevolution between genes belonging to different partners. A well-understood example is secreted fungal effector proteins and their host targets, which play an important role in pathogenic interactions. Biotrophic smut fungi (Basidiomycota) are well-suited to investigate the evolution of plant pathogens, because several reference genomes and genetic tools are available for these species. Here, we used the genomes of Sporisorium reilianum f. sp. zeae and S. reilianum f. sp. reilianum, two closely related formae speciales infecting maize and sorghum, respectively, together with the genomes of Ustilago hordei, Ustilago maydis, and Sporisorium scitamineum to identify and characterize genes displaying signatures of positive selection. We identified 154 gene families having undergone positive selection during species divergence in at least one lineage, among which 77% were identified in the two investigated formae speciales of S. reilianum. Remarkably, only 29% of positively selected genes encode predicted secreted proteins. We assessed the contribution to virulence of nine of these candidate effector genes in S. reilianum f. sp. zeae by deleting individual genes, including a homologue of the effector gene pit2 previously characterized in U. maydis. Only the pit2 deletion mutant was found to be strongly reduced in virulence. Additional experiments are required to understand the molecular mechanisms underlying the selection forces acting on the other candidate effector genes, as well as the large fraction of positively selected genes encoding predicted cytoplasmic proteins.

High-resolution mapping of the recombination landscape of the phytopathogen Fusarium graminearum suggests two-speed genome evolution

Recombination is a major evolutionary force, increasing genetic diversity and permitting efficient coevolution of fungal pathogen(s) with their host(s). The ascomycete Fusarium graminearum is a devastating pathogen of cereal crops, and can contaminate food and feed with harmful mycotoxins. Previous studies have suggested a high adaptive potential of this pathogen, illustrated by an increase in pathogenicity and resistance to fungicides. In this study, we provide the first detailed picture of the crossover events occurring during meiosis and discuss the role of recombination in pathogen evolution. An experimental recombinant population (n = 88) was created and genotyped using 1306 polymorphic markers obtained from restriction site-associated DNA sequencing (RAD-seq) and aligned to the reference genome. The construction of a high-density linkage map, anchoring 99% of the total length of the reference genome, allowed the identification of 1451 putative crossovers, positioned at a median resolution of 24 kb. The majority of crossovers (87.2%) occurred in a relatively small portion of the genome (30%). All chromosomes demonstrated recombination-active sections, which had a near 15-fold higher crossover rate than non-active recombinant sections. The recombination rate showed a strong positive correlation with nucleotide diversity, and recombination-active regions were enriched for genes with a putative role in host-pathogen interaction, as well as putative diversifying genes. Our results confirm the preliminary analysis observed in other F. graminearum strains and suggest a conserved 'two-speed' recombination landscape. The consequences with regard to the evolutionary potential of this major fungal pathogen are also discussed.

A high-resolution genetic map of the cereal crown rot pathogen Fusarium pseudograminearum provides a near-complete genome assembly

Fusarium pseudograminearum is an important pathogen of wheat and barley, particularly in semi-arid environments. Previous genome assemblies for this organism were based entirely on short read data and are highly fragmented. In this work, a genetic map of F. pseudograminearum has been constructed for the first time based on a mapping population of 178 individuals. The genetic map, together with long read scaffolding of a short read-based genome assembly, was used to give a near-complete assembly of the four F. pseudograminearum chromosomes. Large regions of synteny between F. pseudograminearum and F. graminearum, the related pathogen that is the primary causal agent of cereal head blight disease, were previously proposed in the core conserved genome, but the construction of a genetic map to order and orient contigs is critical to the validation of synteny and the placing of species-specific regions. Indeed, our comparative analyses of the genomes of these two related pathogens suggest that rearrangements in the F. pseudograminearum genome have occurred in the chromosome ends. One of these rearrangements includes the transposition of an entire gene cluster involved in the detoxification of the benzoxazolinone (BOA) class of plant phytoalexins. This work provides an important genomic and genetic resource for F. pseudograminearum, which is less well characterized than F. graminearum. In addition, this study provides new insights into a better understanding of the sexual reproduction process in F. pseudograminearum, which informs us of the potential of this pathogen to evolve.

Transcriptome Analysis of Tomato Leaf Spot Pathogen Fusarium proliferatum: De novo Assembly, Expression Profiling, and Identification of Candidate Effectors

Leaf spot disease caused by the fungus Fusarium proliferatum (Matsushima) Nirenberg is a destructive disease of tomato plants in China. Typical symptoms of infected tomato plants are softened and wilted stems and leaves, leading to the eventual death of the entire plant. In this study, we resorted to transcriptional profile analysis to gain insight into the repertoire of effectors involved in F. proliferatum–tomato interactions. A total of 61,544,598 clean reads were de novo assembled to provide a F. proliferatum reference transcriptome. From these, 75,044 unigenes were obtained, with 19.46% of the unigenes being assigned to 276 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, with 22.3% having a homology with genes from F. fujikuroi. A total of 18,075 differentially expressed genes (DEGs) were identified, 720 of which were found to code for secreted proteins. Of these, 184 were identified as candidate effectors, while 79.89% had an upregulated expression. Moreover, 17 genes that were differentially expressed in RNA-seq studies were randomly selected for validation by quantitative real-time polymerase chain reaction (qRT–PCR). The study demonstrates that transcriptome analysis could be an effective method for identifying the repertoire of candidate effectors and may provide an invaluable resource for future functional analyses of F. proliferatum pathogenicity in F. proliferatum and tomato plant–host interactions.

Fungi that Infect Humans

Fungi must meet four criteria to infect humans: growth at human body temperatures, circumvention or penetration of surface barriers, lysis and absorption of tissue, and resistance to immune defenses, including elevated body temperatures. Morphogenesis between small round, detachable cells and long, connected cells is the mechanism by which fungi solve problems of locomotion around or through host barriers. Secretion of lytic enzymes, and uptake systems for the released nutrients, are necessary if a fungus is to nutritionally utilize human tissue. Last, the potent human immune system evolved in the interaction with potential fungal pathogens, so few fungi meet all four conditions for a healthy human host. Paradoxically, the advances of modern medicine have made millions of people newly susceptible to fungal infections by disrupting immune defenses. This article explores how different members of four fungal phyla use different strategies to fulfill the four criteria to infect humans: the Entomophthorales, the Mucorales, the Ascomycota, and the Basidiomycota. Unique traits confer human pathogenic potential on various important members of these phyla: pathogenic Onygenales comprising thermal dimorphs such as Histoplasma and Coccidioides; the Cryptococcus spp. that infect immunocompromised as well as healthy humans; and important pathogens of immunocompromised patients-Candida, Pneumocystis, and Aspergillus spp. Also discussed are agents of neglected tropical diseases important in global health such as mycetoma and paracoccidiomycosis and common pathogens rarely implicated in serious illness such as dermatophytes. Commensalism is considered, as well as parasitism, in shaping genomes and physiological systems of hosts and fungi during evolution.

Secretome Characterization and Correlation Analysis Reveal Putative Pathogenicity Mechanisms and Identify Candidate Avirulence Genes in the Wheat Stripe Rust Fungus Puccinia striiformis f. sp. tritici

Stripe (yellow) rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most destructive diseases of wheat worldwide. Planting resistant cultivars is an effective way to control this disease, but race-specific resistance can be overcome quickly due to the rapid evolving Pst population. Studying the pathogenicity mechanisms is critical for understanding how Pst virulence changes and how to develop wheat cultivars with durable resistance to stripe rust. We re-sequenced 7 Pst isolates and included additional 7 previously sequenced isolates to represent balanced virulence/avirulence profiles for several avirulence loci in seretome analyses. We observed an uneven distribution of heterozygosity among the isolates. Secretome comparison of Pst with other rust fungi identified a large portion of species-specific secreted proteins, suggesting that they may have specific roles when interacting with the wheat host. Thirty-two effectors of Pst were identified from its secretome. We identified candidates for Avr genes corresponding to six Yr genes by correlating polymorphisms for effector genes to the virulence/avirulence profiles of the 14 Pst isolates. The putative AvYr76 was present in the avirulent isolates, but absent in the virulent isolates, suggesting that deleting the coding region of the candidate avirulence gene has produced races virulent to resistance gene Yr76. We conclude that incorporating avirulence/virulence phenotypes into correlation analysis with variations in genomic structure and secretome, particularly presence/absence polymorphisms of effectors, is an efficient way to identify candidate Avr genes in Pst. The candidate effector genes provide a rich resource for further studies to determine the evolutionary history of Pst populations and the co-evolutionary arms race between Pst and wheat. The Avr candidates identified in this study will lead to cloning avirulence genes in Pst, which will enable us to understand molecular mechanisms underlying Pst-wheat interactions, to determine the effectiveness of resistance genes and further to develop durable resistance to stripe rust.

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.

Emerging Fungal Threats to Plants and Animals Challenge Agriculture and Ecosystem Resilience

While fungi can make positive contributions to ecosystems and agro-ecosystems, for example, in mycorrhizal associations, they can also have devastating impacts as pathogens of plants and animals. In undisturbed ecosystems, most such negative interactions will be limited through the coevolution of fungi with their hosts. In this article, we explore what happens when pathogenic fungi spread beyond their natural ecological range and become invasive on naïve hosts in new ecosystems. We will see that such invasive pathogens have been problematic to humans and their domesticated plant and animal species throughout history, and we will discuss some of the most pressing fungal threats of today.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

NOXious gases and the unpredictability of emerging plant pathogens under climate change

Emerging pathogens of crops threaten food security and are increasingly problematic due to intensive agriculture and high volumes of trade and transport in plants and plant products. The ability to predict pathogen risk to agricultural regions would therefore be valuable. However, predictions are complicated by multi-faceted relationships between crops, their pathogens, and climate change. Climate change is related to industrialization, which has brought not only a rise in greenhouse gas emissions but also an increase in other atmospheric pollutants. Here, we consider the implications of rising levels of reactive nitrogen gases and their manifold interactions with crops and crop diseases.

Influence of phylogenetic conservatism and trait convergence on the interactions between fungal root endophytes and plants

Plants associate through their roots with fungal assemblages that impact their abundance and productivity. Non-mycorrhizal endophytes constitute an important component of such fungal diversity, but their implication in ecosystem processes is little known. Using a selection of 128 root-endophytic strains, we defined functional groups based on their traits and plant interactions with potential to predict community assembly and symbiotic association processes. In vitro tests of the strains' interactions with Arabidopsis thaliana, Microthlaspi erraticum and Hordeum vulgare showed a net negative effect of fungal colonization on plant growth. The effects partly depended on the phylogenetic affiliation of strains, but also varied considerably depending on the plant-strain combination. The variation was partly explained by fungal traits shared by different lineages, like growth rates or melanization. The origin of strains also affected their symbioses, with endophytes isolated from Microthlaspi spp. populations being more detrimental to M. erraticum than strains from other sources. Our findings suggest that plant-endophyte associations are subject to local processes of selection, in which particular combinations of symbionts are favored across landscapes. We also show that different common endophytic taxa have differential sets of traits found to affect interactions, hinting to a functional complementarity that can explain their frequent co-existence in natural communities.

Landscape of genomic diversity and host adaptation in Fusarium graminearum

BACKGROUND

Fusarium graminearum is one of the main causal agents of the Fusarium Head Blight, a worldwide disease affecting cereal cultures, whose presence can lead to contaminated grains with chemically stable and harmful mycotoxins. Resistant cultivars and fungicides are frequently used to control this pathogen, and several observations suggest an adaptation of F. graminearum that raises concerns regarding the future of current plant disease management strategies. To understand the genetic basis as well as the extent of its adaptive potential, we investigated the landscape of genomic diversity among six French isolates of F. graminearum, at single-nucleotide resolution using whole-genome re-sequencing.

RESULTS

A total of 242,756 high-confidence genetic variants were detected when compared to the reference genome, among which 96% are single nucleotides polymorphisms. One third of these variants were observed in all isolates. Seventy-seven percent of the total polymorphism is located in 32% of the total length of the genome, comprising telomeric/subtelomeric regions as well as discrete interstitial sections, delineating clear variant enriched genomic regions- 7.5 times in average. About 80% of all the F. graminearum protein-coding genes were found polymorphic. Biological functions are not equally affected: genes potentially involved in host adaptation are preferentially located within polymorphic islands and show greater diversification rate than genes fulfilling basal functions. We further identified 29 putative effector genes enriched with non-synonymous effect mutation.

CONCLUSIONS

Our results highlight a remarkable level of polymorphism in the genome of F. graminearum distributed in a specific pattern. Indeed, the landscape of genomic diversity follows a bi-partite organization of the genome according to polymorphism and biological functions. We measured, for the first time, the level of sequence diversity for the entire gene repertoire of F. graminearum and revealed that the majority are polymorphic. Those assumed to play a role in host-pathogen interaction are discussed, in the light of the subsequent consequences for host adaptation. The annotated genetic variants discovered for this major pathogen are valuable resources for further genetic and genomic studies.

Diverse data supports the transition of filamentous fungal model organisms into the post-genomics era

Filamentous fungi have been important as model organisms since the beginning of modern biological inquiry and have benefitted from open data since the earliest genetic maps were shared. From early origins in simple Mendelian genetics of mating types, parasexual genetics of colony colour, and the foundational demonstration of the segregation of a nutritional requirement, the contribution of research systems utilising filamentous fungi has spanned the biochemical genetics era, through the molecular genetics era, and now are at the very foundation of diverse omics approaches to research and development. Fungal model organisms have come from most major taxonomic groups although Ascomycete filamentous fungi have seen the most major sustained effort. In addition to the published material about filamentous fungi, shared molecular tools have found application in every area of fungal biology. Similarly, shared data has contributed to the success of model systems. The scale of data supporting research with filamentous fungi has grown by 10 to 12 orders of magnitude. From genetic to molecular maps, expression databases, and finally genome resources, the open and collaborative nature of the research communities has assured that the rising tide of data has lifted all of the research systems together.

Characterization of the Two-Speed Subgenomes of Fusarium graminearum Reveals the Fast-Speed Subgenome Specialized for Adaption and Infection

Fusarium head blight, caused by Fusarium graminearum, is one of the most severe diseases on wheat and barley worldwide. Although the genomic data of several strains were published, the intragenomic variation of F. graminearum was not well characterized. Here, we sequenced three Chinese strains and conducted genome-wide comparisons. Our data revealed that all the sequenced strains were distinct from each other and over 350 genes were functionally lost in each of them. Variants of each strain were unevenly distributed in a highly conserved pattern along the chromosomes, resulting in a conserved two-speed genome. The fast subgenome has a lower GC content, shorter gene length, and higher variation of exon numbers than the slow subgenome. Genes related to interaction and pathogenicity, under positive selection, and up-regulated in planta were all significantly enriched in the fast subgenome. Furthermore, we found that the fast subgenome coincided with facultative heterochromatin regions that were repressed in vegetative stage but activated during infection as measured by RNA-seq and ChIP-seq data, suggesting that the fast subgenome is epigenetically regulated. Taken together, our data demonstrated that F. graminearum has a highly conserved two-speed genome and the fast subgenome responsible for adaption and infection is under the control of heterochromatin.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Transcriptome analysis of the fungal pathogen Fusarium oxysporum f. sp. medicaginis during colonisation of resistant and susceptible Medicago truncatula hosts identifies differential pathogenicity profiles and novel candidate effectors

BACKGROUND

Pathogenic members of the Fusarium oxysporum species complex are responsible for vascular wilt disease on many important crops including legumes, where they can be one of the most destructive disease causing necrotrophic fungi. We previously developed a model legume-infecting pathosystem based on the reference legume Medicago truncatula and a pathogenic F. oxysporum forma specialis (f. sp.) medicaginis (Fom). To dissect the molecular pathogenicity arsenal used by this root-infecting pathogen, we sequenced its transcriptome during infection of a susceptible and resistant host accession.

RESULTS

High coverage RNA-Seq of Fom infected root samples harvested from susceptible (DZA315) or resistant (A17) M. truncatula seedlings at early or later stages of infection (2 or 7 days post infection (dpi)) and from vegetative (in vitro) samples facilitated the identification of unique and overlapping sets of in planta differentially expressed genes. This included enrichment, particularly in DZA315 in planta up-regulated datasets, for proteins associated with sugar, protein and plant cell wall metabolism, membrane transport, nutrient uptake and oxidative processes. Genes encoding effector-like proteins were identified, including homologues of the F. oxysporum f. sp. lycopersici Secreted In Xylem (SIX) proteins, and several novel candidate effectors based on predicted secretion, small protein size and high in-planta induced expression. The majority of the effector candidates contain no known protein domains but do share high similarity to predicted proteins predominantly from other F. oxysporum ff. spp. as well as other Fusaria (F. solani, F. fujikori, F. verticilloides, F. graminearum and F. pseudograminearum), and from another wilt pathogen of the same class, a Verticillium species. Overall, this suggests these novel effector candidates may play important roles in Fusaria and wilt pathogen virulence.

CONCLUSION

Combining high coverage in planta RNA-Seq with knowledge of fungal pathogenicity protein features facilitated the identification of differentially expressed pathogenicity associated genes and novel effector candidates expressed during infection of a resistant or susceptible M. truncatula host. The knowledge from this first in depth in planta transcriptome sequencing of any F. oxysporum ff. spp. pathogenic on legumes will facilitate the dissection of Fusarium wilt pathogenicity mechanisms on many important legume crops.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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