A target enrichment approach for enhanced recovery of Synchytrium endobioticum nuclear genome sequences

Potato wart disease is caused by the obligate fungal pathogen Synchytrium endobioticum. DNA extraction from compost, purified spores and crude wart tissue derived from tuber galls of infected potatoes often results in low S. endobioticum DNA concentration or highly contaminated with DNA coming from other microorganisms and the potato host. Therefore, Illumina sequencing of these samples generally results in suboptimal recovery of the nuclear genome sequences of S. endobioticum. A hybridization-based target enrichment protocol was developed to strongly enhance the recovery of S. endobioticum DNA while off-target organisms DNA remains uncaptured. The design strategy involved creating a set of 180,000 molecular baits targeting both gene and non-gene regions of S. endobioticum. The baits were applied to whole genome amplified DNA samples of various S. endobioticum pathotypes (races) in compost, from purified spores and crude wart tissue samples. This was followed by Illumina sequencing and bioinformatic analyses. Compared to non-enriched samples, target enriched samples: 1) showed a significant increase in the proportion of sequenced bases mapped to the S. endobioticum nuclear genome, especially for crude wart tissue samples; 2) yielded sequencing data with higher and better nuclear genome coverage; 3) biased genome assembly towards S. endobioticum sequences, yielding smaller assembly sizes but higher representation of putative S. endobioticum contigs; 4) showed an increase in the number of S. endobioticum genes detected in the genome assemblies. Our hybridization-based target enrichment protocol offers a valuable tool for enhancing genome sequencing and NGS-based molecular detection of S. endobioticum, especially in difficult samples.

An Efficient Triplex TaqMan Quantitative PCR to Detect a Blackleg-Causing Lineage of Pectobacterium brasiliense in Potato Based on a Pangenome Analysis

P. brasiliense is an important bacterial pathogen causing blackleg (BL) in potatoes. Nevertheless, P. brasiliense is often detected in seed lots that do not develop any of the typical blackleg symptoms in the potato crop when planted. Field bioassays identified that P. brasiliense strains can be categorized into two distinct classes, some able to cause blackleg symptoms and some unable to do it. A comparative pangenomic approach was performed on 116 P. brasiliense strains, of which 15 were characterized as BL-causing strains and 25 as non-causative. In a genetically homogeneous clade comprising all BL-causing P. brasiliense strains, two genes only present in the BL-causing strains were identified, one encoding a predicted lysozyme inhibitor Lprl (LZI) and one encoding a putative Toll/interleukin-1 receptor (TIR) domain-containing protein. TaqMan assays for the specific detection of BL-causing P. brasiliense were developed and integrated with the previously developed generic P. brasiliense assay into a triplex TaqMan assay. This simultaneous detection makes the scoring more efficient as only a single tube is needed, and it is more robust as BL-causing strains of P. brasiliense should be positive for all three assays. Individual P. brasiliense strains were found to be either positive for all three assays or only for the P. brasiliense assay. In potato samples, the mixed presence of BL-causing and not BL-causing P. brasiliense strains was observed as shown by the difference in Ct value of the TaqMan assays. However, upon extension of the number of strains, it became clear that in recent years additional BL-causing lineages of P. brasiliense were detected for which additional assays must be developed.

PanTools v3: functional annotation, classification and phylogenomics

SUMMARY

The ever-increasing number of sequenced genomes necessitates the development of pangenomic approaches for comparative genomics. Introduced in 2016, PanTools is a platform that allows pangenome construction, homology grouping and pangenomic read mapping. The use of graph database technology makes PanTools versatile, applicable from small viral genomes like SARS-CoV-2 up to large plant or animal genomes like tomato or human. Here, we present our third major update to PanTools that enables the integration of functional annotations and provides both gene-level analyses and phylogenetics.

AVAILABILITY AND IMPLEMENTATION

PanTools is implemented in Java 8 and released under the GNU GPLv3 license. Software and documentation are available at https://git.wur.nl/bioinformatics/pantools.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Synchytrium endobioticum, the potato wart disease pathogen

Potato wart disease is considered one of the most important quarantine pests for cultivated potato and is caused by the obligate biotrophic chytrid fungus Synchytrium endobioticum. This review integrates observations from early potato wart research and recent molecular, genetic, and genomic studies of the pathogen and its host potato. Taxonomy, epidemiology, pathology, and formation of new pathotypes are discussed, and a model for molecular S. endobioticum-potato interaction is proposed.

TAXONOMY

Currently classified as kingdom: Fungi, phylum: Chytridiomycota, class: Chytridiomycetes, order: Chytridiales, family: Synchytriaceae, genus: Synchytrium, species: Synchytrium endobioticum, there is strong molecular support for Synchytriaceae to be transferred to the order Synchytriales.

HOSTS AND DISEASE SYMPTOMS

Solanum tuberosum is the main host for S. endobioticum but other solanaceous species have been reported as alternative hosts. It is not known if these alternative hosts play a role in the survival of the pathogen in (borders of) infested fields. Disease symptoms on potato tubers are characterized by the warty cauliflower-like malformations that are the result of cell enlargement and cell multiplication induced by the pathogen. Meristematic tissue on tubers, stolons, eyes, sprouts, and inflorescences can be infected while the potato root system seems to be immune.

PATHOTYPES

For S. endobioticum over 40 pathotypes, which are defined as groups of isolates with a similar response to a set of differential potato varieties, are described. Pathotypes 1(D1), 2(G1), 6(O1), and 18(T1) are currently regarded to be most widespread. However, with the current differential set other pathogen diversity largely remains undetected.

PATHOGEN-HOST INTERACTION

A single effector has been described for S. endobioticum (AvrSen1), which is recognized by the potato Sen1 resistance gene product. This is also the first effector that has been described in Chytridiomycota, showing that in this fungal division resistance also fits the gene-for-gene concept. Although significant progress was made in the last decade in mapping wart disease resistance loci, not all resistances present in potato breeding germplasm could be identified. The use of resistant varieties plays an essential role in disease management.

Phylogenomic Analysis of a 55.1-kb 19-Gene Dataset Resolves a Monophyletic Fusarium that Includes the Fusarium solani Species Complex

Scientific communication is facilitated by a data-driven, scientifically sound taxonomy that considers the end-user's needs and established successful practice. In 2013, the Fusarium community voiced near unanimous support for a concept of Fusarium that represented a clade comprising all agriculturally and clinically important Fusarium species, including the F. solani species complex (FSSC). Subsequently, this concept was challenged in 2015 by one research group who proposed dividing the genus Fusarium into seven genera, including the FSSC described as members of the genus Neocosmospora, with subsequent justification in 2018 based on claims that the 2013 concept of Fusarium is polyphyletic. Here, we test this claim and provide a phylogeny based on exonic nucleotide sequences of 19 orthologous protein-coding genes that strongly support the monophyly of Fusarium including the FSSC. We reassert the practical and scientific argument in support of a genus Fusarium that includes the FSSC and several other basal lineages, consistent with the longstanding use of this name among plant pathologists, medical mycologists, quarantine officials, regulatory agencies, students, and researchers with a stake in its taxonomy. In recognition of this monophyly, 40 species described as genus Neocosmospora were recombined in genus Fusarium, and nine others were renamed Fusarium. Here the global Fusarium community voices strong support for the inclusion of the FSSC in Fusarium, as it remains the best scientific, nomenclatural, and practical taxonomic option available.

The Pectobacterium pangenome, with a focus on Pectobacterium brasiliense, shows a robust core and extensive exchange of genes from a shared gene pool

BACKGROUND

Bacterial plant pathogens of the Pectobacterium genus are responsible for a wide spectrum of diseases in plants, including important crops such as potato, tomato, lettuce, and banana. Investigation of the genetic diversity underlying virulence and host specificity can be performed at genome level by using a comprehensive comparative approach called pangenomics. A pangenomic approach, using newly developed functionalities in PanTools, was applied to analyze the complex phylogeny of the Pectobacterium genus. We specifically used the pangenome to investigate genetic differences between virulent and avirulent strains of P. brasiliense, a potato blackleg causing species dominantly present in Western Europe.

RESULTS

Here we generated a multilevel pangenome for Pectobacterium, comprising 197 strains across 19 species, including type strains, with a focus on P. brasiliense. The extensive phylogenetic analysis of the Pectobacterium genus showed robust distinct clades, with most detail provided by 452,388 parsimony-informative single-nucleotide polymorphisms identified in single-copy orthologs. The average Pectobacterium genome consists of 47% core genes, 1% unique genes, and 52% accessory genes. Using the pangenome, we zoomed in on differences between virulent and avirulent P. brasiliense strains and identified 86 genes associated to virulent strains. We found that the organization of genes is highly structured and linked with gene conservation, function, and transcriptional orientation.

CONCLUSION

The pangenome analysis demonstrates that evolution in Pectobacteria is a highly dynamic process, including gene acquisitions partly in clusters, genome rearrangements, and loss of genes. Pectobacterium species are typically not characterized by a set of species-specific genes, but instead present themselves using new gene combinations from the shared gene pool. A multilevel pangenomic approach, fusing DNA, protein, biological function, taxonomic group, and phenotypes, facilitates studies in a flexible taxonomic context.

Detecting Introgression Between Members of the Fusarium fujikuroi and F. oxysporum Species Complexes by Comparative Mitogenomics

The Fusarium fujikuroi species complex (FFSC) and F. oxysporum species complex (FOSC) are two related groups of plant pathogens causing a wide diversity of diseases in agricultural crops world wide. The aims of this study are (1) to clarify the phylogeny of the FFSC, (2) to identify potential deviation from tree-like evolution, (3) to explore the value of using mitogenomes for these kinds of analyses, and (4) to better understand mitogenome evolution. In total, we have sequenced 24 species from the FFSC and a representative set of recently analyzed FOSC strains was chosen, while F. redolens was used as outgroup for the two species complexes. A species tree was constructed based on the concatenated alignment of seven nuclear genes and the mitogenome, which was contrasted to individual gene trees to identify potential conflicts. These comparisons indicated conflicts especially within the previously described African clade of the FFSC. Furthermore, the analysis of the mitogenomes revealed the presence of a variant of the large variable (LV) region in FFSC which was previously only reported for FOSC. The distribution of this variant and the results of sequence comparisons indicate horizontal genetic transfer between members of the two species complexes, most probably through introgression. In addition, a duplication of atp9 was found inside an intron of cob, which suggests that even highly conserved mitochondrial genes can have paralogs. Paralogization in turn may lead to inaccurate single gene phylogenies. In conclusion, mitochondrial genomes provide a robust basis for phylogeny. Comparative phylogenetic analysis indicated that gene flow among and between members of FFSC and FOSC has played an important role in the evolutionary history of these two groups. Since mitogenomes show greater levels of conservation and synteny than nuclear regions, they are more likely to be compatible for recombination than nuclear regions. Therefore, mitogenomes can be used as indicators to detect interspecies gene flow.

Population Genomic Analysis Reveals a Highly Conserved Mitochondrial Genome in Fusarium asiaticum

Fusarium asiaticum is one of the pivotal members of the Fusarium graminearum species complex (FGSC) causing Fusarium head blight (FHB) on wheat, barley and rice in large parts of Asia. Besides resulting in yield losses, FHB also causes the accumulation of mycotoxins such as nivalenol (NIV) and deoxynivalenol (DON). The aim of this study was to conduct population studies on F. asiaticum from Southern China through mitochondrial genome analyses. All strains were isolated from wheat or rice from several geographic areas in seven provinces in Southern China. Based on geographic location and host, 210 isolates were selected for next generation sequencing, and their mitogenomes were assembled by GRAbB and annotated to explore the mitochondrial genome variability of F. asiaticum. The F. asiaticum mitogenome proves extremely conserved and variation is mainly caused by absence/presence of introns harboring homing endonuclease genes. These variations could be utilized to develop molecular markers for track and trace of migrations within and between populations. This study illustrates how mitochondrial introns can be used as markers for population genetic analysis. SNP analysis demonstrate the occurrence of mitochondrial recombination in F. asiaticum as was previously found for F. oxysporum and implied for F. graminearum. Furthermore, varying degrees of genetic diversity and recombination showed a high association with different geographic regions as well as with cropping systems. The mitogenome of F. graminearum showed a much higher SNP diversity while the interspecies intron variation showed no evidence of gene flow between the two closely related and sexual compatible species.

The Synchytrium endobioticum AvrSen1 Triggers a Hypersensitive Response in Sen1 Potatoes While Natural Variants Evade Detection

Synchytrium endobioticum is an obligate biotrophic fungus of division Chytridiomycota. It causes potato wart disease, has a worldwide quarantine status and is included on the Health and Human Services and United States Department of Agriculture Select Agent list. S. endobioticum isolates are grouped in pathotypes based on their ability to evade host resistance in a set of differential potato varieties. Thus far, 39 pathotypes are reported. A single dominant gene (Sen1) governs pathotype 1 (D1) resistance and we anticipated that the underlying molecular model would involve a pathogen effector (AvrSen1) that is recognized by the host. The S. endobioticum-specific secretome of 14 isolates representing six different pathotypes was screened for effectors specifically present in pathotype 1 (D1) isolates but absent in others. We identified a single AvrSen1 candidate. Expression of this candidate in potato Sen1 plants showed a specific hypersensitive response (HR), which cosegregated with the Sen1 resistance in potato populations. No HR was obtained with truncated genes found in pathotypes that evaded recognition by Sen1. These findings established that our candidate gene was indeed Avrsen1. The S. endobioticum AvrSen1 is a single-copy gene and encodes a 376-amino-acid protein without predicted function or functional domains, and is the first effector gene identified in Chytridiomycota, an extremely diverse yet underrepresented basal lineage of fungi.

The Synchytrium endobioticum AvrSen1 Triggers a Hypersensitive Response in Sen1 Potatoes While Natural Variants Evade Detection

Synchytrium endobioticum is an obligate biotrophic fungus of division Chytridiomycota. It causes potato wart disease, has a worldwide quarantine status and is included on the Health and Human Services and United States Department of Agriculture Select Agent list. S. endobioticum isolates are grouped in pathotypes based on their ability to evade host resistance in a set of differential potato varieties. Thus far, 39 pathotypes are reported. A single dominant gene (Sen1) governs pathotype 1 (D1) resistance and we anticipated that the underlying molecular model would involve a pathogen effector (AvrSen1) that is recognized by the host. The S. endobioticum-specific secretome of 14 isolates representing six different pathotypes was screened for effectors specifically present in pathotype 1 (D1) isolates but absent in others. We identified a single AvrSen1 candidate. Expression of this candidate in potato Sen1 plants showed a specific hypersensitive response (HR), which cosegregated with the Sen1 resistance in potato populations. No HR was obtained with truncated genes found in pathotypes that evaded recognition by Sen1. These findings established that our candidate gene was indeed Avrsen1. The S. endobioticum AvrSen1 is a single-copy gene and encodes a 376-amino-acid protein without predicted function or functional domains, and is the first effector gene identified in Chytridiomycota, an extremely diverse yet underrepresented basal lineage of fungi.

Comparative genomics of chytrid fungi reveal insights into the obligate biotrophic and pathogenic lifestyle of Synchytrium endobioticum

Synchytrium endobioticum is an obligate biotrophic soilborne Chytridiomycota (chytrid) species that causes potato wart disease, and represents the most basal lineage among the fungal plant pathogens. We have chosen a functional genomics approach exploiting knowledge acquired from other fungal taxa and compared this to several saprobic and pathogenic chytrid species. Observations linked to obligate biotrophy, genome plasticity and pathogenicity are reported. Essential purine pathway genes were found uniquely absent in S. endobioticum, suggesting that it relies on scavenging guanine from its host for survival. The small gene-dense and intron-rich chytrid genomes were not protected for genome duplications by repeat-induced point mutation. Both pathogenic chytrids Batrachochytrium dendrobatidis and S. endobioticum contained the largest amounts of repeats, and we identified S. endobioticum specific candidate effectors that are associated with repeat-rich regions. These candidate effectors share a highly conserved motif, and show isolate specific duplications. A reduced set of cell wall degrading enzymes, and LysM protein expansions were found in S. endobioticum, which may prevent triggering plant defense responses. Our study underlines the high diversity in chytrids compared to the well-studied Ascomycota and Basidiomycota, reflects characteristic biological differences between the phyla, and shows commonalities in genomic features among pathogenic fungi.

Development and Evaluation of a Triplex TaqMan Assay and Next-Generation Sequence Analysis for Improved Detection of Xylella in Plant Material

Xylella fastidiosa is a heterogenous gram-negative bacterial plant pathogen with a wide host range covering over 300 plant species. Since 2013, in Europe, the presence of the pathogen is increasing in a part of the Mediterranean area, but it causes in particular severe disease problems in olive orchards in the Southern part of Italy. Various subspecies of the pathogen were also diagnosed in natural outbreaks and intercepted ornamental plants in Europe, among them Olea europaea, Coffea arabica, and Nerium oleander. The host range of the pathogen can vary, depending on the subspecies and even the strain. The availability of fast and reliable diagnostic tools is indispensable in management strategies to control diseases caused by X. fastidiosa. To improve the reliability of the TaqMan assay, currently widely used in surveys, a triplex TaqMan assay was developed in which two specific and sensitive TaqMan assays, previously designed for X. fastidiosa, were combined with an internal control. The triplex assay exhibited the same diagnostic sensitivity as the simplex assays. In addition, the usefulness of a metagenomic approach using next-generation sequencing (NGS) was demonstrated, in which total DNA extracted from plant material was sequenced. DNA extracts from plant material free of X. fastidiosa, from artificially inoculated hosts plants or from naturally infected plants sampled in France, Spain, and Italy, or intercepted in Austria and the Netherlands, were analyzed for the presence of X. fastidiosa using the metagenomic approach. In all samples, even in samples with a low infection level, but not in the pathogen-free samples, DNA reads were detected specific for X. fastidiosa. In most cases, the pathogen could be identified to the subspecies level, and for one sample even the whole genome could be assembled and the sequence type could be determined. All results of NGS-analyzed samples were confirmed with the triplex TaqMan polymerase chain reaction and loop-mediated isothermal amplification.

The linear mitochondrial genome of the quarantine chytrid Synchytrium endobioticum; insights into the evolution and recent history of an obligate biotrophic plant pathogen

BACKGROUND

Chytridiomycota species (chytrids) belong to a basal lineage in the fungal kingdom. Inhabiting terrestrial and aquatic environments, most are free-living saprophytes but several species cause important diseases: e.g. Batrachochytrium dendrobatidis, responsible for worldwide amphibian decline; and Synchytrium endobioticum, causing potato wart disease. S. endobioticum has an obligate biotrophic lifestyle and isolates can be further characterized as pathotypes based on their virulence on a differential set of potato cultivars. Quarantine measures have been implemented globally to control the disease and prevent its spread. We used a comparative approach using chytrid mitogenomes to determine taxonomical relationships and to gain insights into the evolution and recent history of introductions of this plant pathogen.

RESULTS

We assembled and annotated the complete mitochondrial genome of 30 S. endobioticum isolates and generated mitochondrial genomes for five additional chytrid species. The mitochondrial genome of S. endobioticum is linear with terminal inverted repeats which was validated by tailing and PCR amplifying the telomeric ends. Surprisingly, no conservation in organisation and orientation of mitochondrial genes was observed among the Chytridiomycota except for S. endobioticum and its sister species Synchytrium microbalum. However, the mitochondrial genome of S. microbalum is circular and comprises only a third of the 72.9 Kbp found for S. endobioticum suggesting recent linearization and expansion. Four mitochondrial lineages were identified in the S. endobioticum mitochondrial genomes. Several pathotypes occur in different lineages, suggesting that these have emerged independently. In addition, variations for polymorphic sites in the mitochondrial genome of individual isolates were observed demonstrating that S. endobioticum isolates represent a community of different genotypes. Such communities were shown to be complex and stable over time, but we also demonstrate that the use of semi-resistant potato cultivars triggers a rapid shift in the mitochondrial haplotype associated with increased virulence.

CONCLUSIONS

Mitochondrial genomic variation shows that S. endobioticum has been introduced into Europe multiple times, that several pathotypes emerged multiple times, and that isolates represent communities of different genotypes. Our study represents the most comprehensive dataset of chytrid mitogenomes, which provides new insights into the extraordinary dynamics and evolution of mitochondrial genomes involving linearization, expansion and reshuffling.

Evolution and Diversity of Biosynthetic Gene Clusters in Fusarium

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

The Fusarium graminearum Histone Acetyltransferases Are Important for Morphogenesis, DON Biosynthesis, and Pathogenicity

Post-translational modifications of chromatin structure by histone acetyltransferase (HATs) play a central role in the regulation of gene expression and various biological processes in eukaryotes. Although HAT genes have been studied in many fungi, few of them have been functionally characterized. In this study, we identified and characterized four putative HATs (FgGCN5, FgRTT109, FgSAS2, FgSAS3) in the plant pathogenic ascomycete Fusarium graminearum, the causal agent of Fusarium head blight of wheat and barley. We replaced the genes and all mutant strains showed reduced growth of F. graminearum. The ΔFgSAS3 and ΔFgGCN5 mutant increased sensitivity to oxidative and osmotic stresses. Additionally, ΔFgSAS3 showed reduced conidia sporulation and perithecium formation. Mutant ΔFgGCN5 was unable to generate any conidia and lost its ability to form perithecia. Our data showed also that FgSAS3 and FgGCN5 are pathogenicity factors required for infecting wheat heads as well as tomato fruits. Importantly, almost no Deoxynivalenol (DON) was produced either in ΔFgSAS3 or ΔFgGCN5 mutants, which was consistent with a significant downregulation of TRI genes expression. Furthermore, we discovered for the first time that FgSAS3 is indispensable for the acetylation of histone site H3K4, while FgGCN5 is essential for the acetylation of H3K9, H3K18, and H3K27. H3K14 can be completely acetylated when FgSAS3 and FgGCN5 were both present. The RNA-seq analyses of the two mutant strains provide insight into their functions in development and metabolism. Results from this study clarify the functional divergence of HATs in F. graminearum, and may provide novel targeted strategies to control secondary metabolite expression and infections of F. graminearum.

Stress and sexual reproduction affect the dynamics of the wheat pathogen effector AvrStb6 and strobilurin resistance

Host resistance and fungicide treatments are cornerstones of plant-disease control. Here, we show that these treatments allow sex and modulate parenthood in the fungal wheat pathogen Zymoseptoria tritici. We demonstrate that the Z. tritici-wheat interaction complies with the gene-for-gene model by identifying the effector AvrStb6, which is recognized by the wheat resistance protein Stb6. Recognition triggers host resistance, thus implying removal of avirulent strains from pathogen populations. However, Z. tritici crosses on wheat show that sex occurs even with an avirulent parent, and avirulence alleles are thereby retained in subsequent populations. Crossing fungicide-sensitive and fungicide-resistant isolates under fungicide pressure results in a rapid increase in resistance-allele frequency. Isolates under selection always act as male donors, and thus disease control modulates parenthood. Modeling these observations for agricultural and natural environments reveals extended durability of host resistance and rapid emergence of fungicide resistance. Therefore, fungal sex has major implications for disease control.

Mitochondrial genomes reveal recombination in the presumed asexual Fusarium oxysporum species complex

BACKGROUND

The Fusarium oxysporum species complex (FOSC) contains several phylogenetic lineages. Phylogenetic studies identified two to three major clades within the FOSC. The mitochondrial sequences are highly informative phylogenetic markers, but have been mostly neglected due to technical difficulties.

RESULTS

A total of 61 complete mitogenomes of FOSC strains were de novo assembled and annotated. Length variations and intron patterns support the separation of three phylogenetic species. The variable region of the mitogenome that is typical for the genus Fusarium shows two new variants in the FOSC. The variant typical for Fusarium is found in members of all three clades, while variant 2 is found in clades 2 and 3 and variant 3 only in clade 2. The extended set of loci analyzed using a new implementation of the genealogical concordance species recognition method support the identification of three phylogenetic species within the FOSC. Comparative analysis of the mitogenomes in the FOSC revealed ongoing mitochondrial recombination within, but not between phylogenetic species.

CONCLUSIONS

The recombination indicates the presence of a parasexual cycle in F. oxysporum. The obstacles hindering the usage of the mitogenomes are resolved by using next generation sequencing and selective genome assemblers, such as GRAbB. Complete mitogenome sequences offer a stable basis and reference point for phylogenetic and population genetic studies.

GRAbB: Selective Assembly of Genomic Regions, a New Niche for Genomic Research

GRAbB (Genomic Region Assembly by Baiting) is a new program that is dedicated to assemble specific genomic regions from NGS data. This approach is especially useful when dealing with multi copy regions, such as mitochondrial genome and the rDNA repeat region, parts of the genome that are often neglected or poorly assembled, although they contain interesting information from phylogenetic or epidemiologic perspectives, but also single copy regions can be assembled. The program is capable of targeting multiple regions within a single run. Furthermore, GRAbB can be used to extract specific loci from NGS data, based on homology, like sequences that are used for barcoding. To make the assembly specific, a known part of the region, such as the sequence of a PCR amplicon or a homologous sequence from a related species must be specified. By assembling only the region of interest, the assembly process is computationally much less demanding and may lead to assemblies of better quality. In this study the different applications and functionalities of the program are demonstrated such as: exhaustive assembly (rDNA region and mitochondrial genome), extracting homologous regions or genes (IGS, RPB1, RPB2 and TEF1a), as well as extracting multiple regions within a single run. The program is also compared with MITObim, which is meant for the exhaustive assembly of a single target based on a similar query sequence. GRAbB is shown to be more efficient than MITObim in terms of speed, memory and disk usage. The other functionalities (handling multiple targets simultaneously and extracting homologous regions) of the new program are not matched by other programs. The program is available with explanatory documentation at https://github.com/b-brankovics/grabb. GRAbB has been tested on Ubuntu (12.04 and 14.04), Fedora (23), CentOS (7.1.1503) and Mac OS X (10.7). Furthermore, GRAbB is available as a docker repository: brankovics/grabb (https://hub.docker.com/r/brankovics/grabb/).

Development of Polymorphic Microsatellite Loci for Potato Wart from Next-Generation Sequence Data

Synchytrium endobioticum is the fungal agent causing potato wart disease. Because of its severity and persistence, quarantine measures are enforced worldwide to avoid the spread of this disease. Molecular markers exist for species-specific detection of this pathogen, yet markers to study the intraspecific genetic diversity of S. endobioticum were not available. Whole-genome sequence data from Dutch pathotype 1 isolate MB42 of S. endobioticum were mined for perfect microsatellite motifs. Of the 62 selected microsatellites, 21 could be amplified successfully and displayed moderate levels of polymorphism in 22 S. endobioticum isolates from different countries. Nineteen multilocus genotypes were observed, with only three isolates from Canada displaying identical profiles. The majority of isolates from Canada clustered genetically. In contrast, most isolates collected in Europe show no genetic clustering associated with their geographic origin. S. endobioticum isolates with the same pathotype displayed highly variable genotypes and none of the microsatellite markers correlated with a specific pathotype. The markers developed in this study can be used to assess intraspecific genetic diversity of S. endobioticum and allow track and trace of genotypes that will generate a better understanding of the migration and spread of this important fungal pathogen and support management of this disease.

Proteome catalog of Zymoseptoria tritici captured during pathogenesis in wheat

Zymoseptoria tritici is an economically important pathogen of wheat. However, the molecular basis of pathogenicity on wheat is still poorly understood. Here, we present a global survey of the proteins secreted by this fungus in the apoplast of resistant (cv. Shafir) and susceptible (cv. Obelisk) wheat cultivars after inoculation with reference Z. tritici strain IPO323. The fungal proteins present in apoplastic fluids were analyzed by gel electrophoresis and by data-independent acquisition liquid chromatography/mass spectrometry (LC/MS(E)) combined with data-dependent acquisition LC-MS/MS. Subsequent mapping mass spectrometry-derived peptide sequence data against the genome sequence of strain IPO323 identified 665 peptides in the MS(E) and 93 in the LC-MS/MS mode that matched to 85 proteins. The identified fungal proteins, including cell-wall degrading enzymes and proteases, might function in pathogenicity, but the functions of many remain unknown. Most fungal proteins accumulated in cv. Obelisk at the onset of necrotrophy. This inventory provides an excellent basis for future detailed studies on the role of these genes and their encoded proteins during pathogenesis in wheat.

Computational strategies for genome-based natural product discovery and engineering in fungi

Fungal natural products possess biological activities that are of great value to medicine, agriculture and manufacturing. Recent metagenomic studies accentuate the vastness of fungal taxonomic diversity, and the accompanying specialized metabolic diversity offers a great and still largely untapped resource for natural product discovery. Although fungal natural products show an impressive variation in chemical structures and biological activities, their biosynthetic pathways share a number of key characteristics. First, genes encoding successive steps of a biosynthetic pathway tend to be located adjacently on the chromosome in biosynthetic gene clusters (BGCs). Second, these BGCs are often are located on specific regions of the genome and show a discontinuous distribution among evolutionarily related species and isolates. Third, the same enzyme (super)families are often involved in the production of widely different compounds. Fourth, genes that function in the same pathway are often co-regulated, and therefore co-expressed across various growth conditions. In this mini-review, we describe how these partly interlinked characteristics can be exploited to computationally identify BGCs in fungal genomes and to connect them to their products. Particular attention will be given to novel algorithms to identify unusual classes of BGCs, as well as integrative pan-genomic approaches that use a combination of genomic and metabolomic data for parallelized natural product discovery across multiple strains. Such novel technologies will not only expedite the natural product discovery process, but will also allow the assembly of a high-quality toolbox for the re-design or even de novo design of biosynthetic pathways using synthetic biology approaches.

Effector discovery in the fungal wheat pathogen Zymoseptoria tritici

Fungal plant pathogens, such as Zymoseptoria tritici (formerly known as Mycosphaerella graminicola), secrete repertoires of effectors to facilitate infection or trigger host defence mechanisms. The discovery and functional characterization of effectors provides valuable knowledge that can contribute to the design of new and effective disease management strategies. Here, we combined bioinformatics approaches with expression profiling during pathogenesis to identify candidate effectors of Z. tritici. In addition, a genetic approach was conducted to map quantitative trait loci (QTLs) carrying putative effectors, enabling the validation of both complementary strategies for effector discovery. In planta expression profiling revealed that candidate effectors were up-regulated in successive waves corresponding to consecutive stages of pathogenesis, contrary to candidates identified by QTL mapping that were, overall, expressed at low levels. Functional analyses of two top candidate effectors (SSP15 and SSP18) showed their dispensability for Z. tritici pathogenesis. These analyses reveal that generally adopted criteria, such as protein size, cysteine residues and expression during pathogenesis, may preclude an unbiased effector discovery. Indeed, genetic mapping of genomic regions involved in specificity render alternative effector candidates that do not match the aforementioned criteria, but should nevertheless be considered as promising new leads for effectors that are crucial for the Z. tritici-wheat pathosystem.

Diagnosis of Fusarium Infections: Approaches to Identification by the Clinical Mycology Laboratory

Infections caused by the genus Fusarium have emerged over the past decades and range from onychomycosis and keratitis in healthy individuals to deep and disseminated infections with high mortality rates in immune-compromised patients. As antifungal susceptibility can differ between the different Fusarium species, identification at species level is recommended. Several clinical observations as hyaline hyphae in tissue, necrotic lesions in the skin and positive blood tests with fungal growth or presence of fungal cell wall components may be the first hints for fusariosis. Many laboratories rely on morphological identification, but especially multi-locus sequencing proves better to discriminate among members of the species complexes involved in human infection. DNA-based diagnostic tools have best discriminatory power when based on translation elongation factor 1-α or the RNA polymerase II second largest subunit. However, assays based on the detection of other fusarial cell compounds such as peptides and cell wall components may also be used for identification. The purpose of this review is to provide an overview and a comparison of the different tools currently available for the diagnosis of fusariosis.

Rapid identification of Fusarium graminearum species complex using Rolling Circle Amplification (RCA)

Rolling Circle Amplification (RCA) of DNA is a sensitive and cost effective method for the rapid identification of pathogenic fungi without the need for sequencing. Amplification products can be visualized on 1% agarose gel to verify the specificity of probe-template binding or directly by adding fluorescent dyes. Fusarium Head Blight (FHB) is currently the world's largest threat to the production of cereal crops with the production of a range of mycotoxins as an additional risk. We designed sets of RCA padlock probes based on polymorphisms in the elongation factor 1-α (EF-1α) gene to detect the dominant FHB species, comprising lineages of the Fusarium graminearum species complex (FGSC). The method also enabled the identification of species of the Fusarium oxysporum (FOSC), the Fusarium incarnatum-equiseti (FIESC), and the Fusarium tricinctum (FTSC) species complexes, and used strains from the CBS culture collection as reference. Subsequently probes were applied to characterize isolates from wheat and wild grasses, and inoculated wheat kernels. The RCA assays successfully amplified DNA of the target fungi, both in environmental samples and in the contaminated wheat samples, while no cross reactivity was observed with uncontaminated wheat or related Fusarium species. As RCA does not require expensive instrumentation, the technique has a good potential for local and point of care screening for toxigenic Fusarium species in cereals.

New broad-spectrum resistance to septoria tritici blotch derived from synthetic hexaploid wheat

Septoria tritici blotch (STB), caused by the ascomycete Mycosphaerella graminicola, is one of the most devastating foliar diseases of wheat. We screened five synthetic hexaploid wheats (SHs), 13 wheat varieties that represent the differential set of cultivars and two susceptible checks with a global set of 20 isolates and discovered exceptionally broad STB resistance in SHs. Subsequent development and analyses of recombinant inbred lines (RILs) from a cross between the SH M3 and the highly susceptible bread wheat cv. Kulm revealed two novel resistance loci on chromosomes 3D and 5A. The 3D resistance was expressed in the seedling and adult plant stages, and it controlled necrosis (N) and pycnidia (P) development as well as the latency periods of these parameters. This locus, which is closely linked to the microsatellite marker Xgwm494, was tentatively designated Stb16q and explained from 41 to 71% of the phenotypic variation at seedling stage and 28-31% in mature plants. The resistance locus on chromosome 5A was specifically expressed in the adult plant stage, associated with SSR marker Xhbg247, explained 12-32% of the variation in disease, was designated Stb17, and is the first unambiguously identified and named QTL for adult plant resistance to M. graminicola. Our results confirm that common wheat progenitors might be a rich source of new Stb resistance genes/QTLs that can be deployed in commercial breeding programs.

Genetic analysis of resistance to septoria tritici blotch in the French winter wheat cultivars Balance and Apache

The ascomycete Mycosphaerella graminicola is the causal agent of septoria tritici blotch (STB), one of the most destructive foliar diseases of bread and durum wheat globally, particularly in temperate humid areas. A screening of the French bread wheat cultivars Apache and Balance with 30 M. graminicola isolates revealed a pattern of resistant responses that suggested the presence of new genes for STB resistance. Quantitative trait loci (QTL) analysis of a doubled haploid (DH) population with five M. graminicola isolates in the seedling stage identified four QTLs on chromosomes 3AS, 1BS, 6DS and 7DS, and occasionally on 7DL. The QTL on chromosome 6DS flanked by SSR markers Xgpw5176 and Xgpw3087 is a novel QTL that now can be designated as Stb18. The QTLs on chromosomes 3AS and 1BS most likely represent Stb6 and Stb11, respectively, and the QTLs on chromosome 7DS are most probably identical with Stb4 and Stb5. However, the QTL identified on chromosome 7DL is expected to be a new Stb gene that still needs further characterization. Multiple isolates were used and show that not all isolates identify all QTLs, which clearly demonstrates the specificity in the M. graminicola-wheat pathosystem. QTL analyses were performed with various disease parameters. The development of asexual fructifications (pycnidia) in the characteristic necrotic blotches of STB, designated as parameter P, identified the maximum number of QTLs. All other parameters identified fewer but not different QTLs. The segregation of multiple QTLs in the Apache/Balance DH population enabled the identification of DH lines with single QTLs and multiple QTL combinations. Analyses of the marker data of these DH lines clearly demonstrated the positive effect of pyramiding QTLs to broaden resistance spectra as well as epistatic and additive interactions between these QTLs. Phenotyping of the Apache/Balance DH population in the field confirmed the presence of the QTLs that were identified in the seedling stage, but Stb18 was inconsistently expressed and might be particularly effective in young plants. In contrast, an additional QTL for STB resistance was identified on chromosome 2DS that is exclusively and consistently expressed in mature plants over locations and time, but it was also strongly related with earliness, tallness as well as resistance to Fusarium head blight. Although to date no Stb gene has been reported on chromosome 2D, the data provide evidence that this QTL is only indirectly related to STB resistance. This study shows that detailed genetic analysis of contemporary commercial bread wheat cultivars can unveil novel Stb genes that can be readily applied in marker-assisted breeding programs.

The het-c heterokaryon incompatibility gene in Aspergillus niger

Heterokaryon incompatibility among Aspergillus niger strains is a widespread phenomenon that is observed as the inability to form stable heterokaryons. The genetic basis of heterokaryon incompatibility reactions is well established in some sexual filamentous fungi but largely unknown in presumed asexual species, such as A. niger. To test whether the genes that determine heterokaryon incompatibility in Neurospora crassa, such as het-c, vib-1 and pin-c, have a similar function in A. niger, we performed a short in silico search for homologues of these genes in the A. niger and several related genomes. For het-c, pin-c and vib-1 we did indeed identify putative orthologues. We then screened a genetically diverse worldwide collection of incompatible black Aspergilli for polymorphisms in the het-c orthologue. No size variation was observed in the variable het-c indel region that determines the specificity in N. crassa. Sequence comparison showed only minor variation in the number of glutamine coding triplets. However, introduction of one of the three N. crassa alleles (het-c2) in A. niger by transformation resulted in an abortive phenotype, reminiscent of the heterokaryon incompatibility in N. crassa. We conclude that although the genes required are present and the het-c homologue could potentially function as a heterokaryon incompatibility gene, het-c has no direct function in heterokaryon incompatibility in A. niger because the necessary allelic variation is absent.

Large-scale gene discovery in the septoria tritici blotch fungus Mycosphaerella graminicola with a focus on in planta expression

The foliar disease septoria tritici blotch, caused by the fungus Mycosphaerella graminicola, is currently the most important wheat disease in Europe. Gene expression was examined under highly different conditions, using 10 expressed sequence tag libraries generated from M. graminicola isolate IPO323 using seven in vitro and three in planta growth conditions. To identify fungal clones in the interaction libraries, we developed a selection method based on hybridization with the entire genomic DNA of M. graminicola, to selectively enrich these libraries for fungal genes. Assembly of the 27,007 expressed sequence tags resulted in 9,190 unigenes, representing 5.2 Mb of the estimated 39-Mb genome size of M. graminicola. All libraries contributed significantly to the number of unigenes, especially the in planta libraries representing different stages of pathogenesis, which covered 15% of the library-specific unigenes. Even under presymptomatic conditions (5 days postinoculation), when fungal biomass is less than 5%, this method enabled us to efficiently capture fungal genes expressed during pathogenesis. Many of these genes were uniquely expressed in planta, indicating that in planta gene expression significantly differed from in vitro expression. Examples of gene discovery included a number of cell wall-degrading enzymes, a broad set of genes involved in signal transduction (n=11) and a range of ATP-binding cassette (n=20) and major facilitator superfamily transporter genes (n=12) potentially involved in protection against antifungal compounds or the secretion of pathogenicity factors. In addition, evidence is provided for a mycovirus in M. graminicola that is highly expressed under various stress conditions, in particular, under nitrogen starvation. Our analyses provide a unique window on in vitro and in planta gene expression of M. graminicola.

Identification and genetic mapping of highly polymorphic microsatellite loci from an EST database of the septoria tritici blotch pathogen Mycosphaerella graminicola

A database of 30,137 EST sequences from Mycosphaerella graminicola, the septoria tritici blotch fungus of wheat, was scanned with a custom software pipeline for di- and trinucleotide units repeated tandemly six or more times. The bioinformatics analysis identified 109 putative SSR loci, and for 99 of them, flanking primers were developed successfully and tested for amplification and polymorphism by PCR on five field isolates of diverse origin, including the parents of the standard M. graminicola mapping population. Seventy-seven of the 99 primer pairs generated an easily scored banding pattern and 51 were polymorphic, with up to four alleles per locus, among the isolates tested. Among these 51 loci, 23 were polymorphic between the parents of the mapping population. Twenty-one of these as well as two previously published microsatellite loci were positioned on the existing genetic linkage map of M. graminicola on 13 of the 24 linkage groups. Most (66%) of the primer pairs also amplified bands in the closely related barley pathogen Septoria passerinii, but only six were polymorphic among four isolates tested. A subset of the primer pairs also revealed polymorphisms when tested with DNA from the related banana black leaf streak (Black Sigatoka) pathogen, M. fijiensis. The EST database provided an excellent source of new, highly polymorphic microsatellite markers that can be multiplexed for high-throughput genetic analyses of M. graminicola and related species.