Living apart together: crosstalk between the core and supernumerary genomes in a fungal plant pathogen
© The Author(s). 2016
Received: 30 January 2016
Accepted: 14 July 2016
Published: 23 August 2016
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.
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.
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.
KeywordsSingle-molecule real-time sequencing Supernumerary chromosomes Fusarium Repeat-induced point mutation Translocation Gene duplications Transposable elements
Genome plasticity is one of the most important drivers of evolutions in eukaryotes. This plasticity includes large scale genome duplications, rearrangements, deletions and compartmentalization [1–3]. Fungi exhibit these dynamics better than any other kingdom, possible due to their occurrence in highly diverse niches. The organization of fungal genomes varies remarkably and in many cases facilitates rapid evolution and speciation . Indeed, the division of the genome into a core and accessory part evolving at different speeds, has been described in multiple pathogens, particularly fungi [5–7]. Supernumerary chromosomes are one of most radical extensions of genome plasticity in fungi. They represent chromosomal structures that vary in size and distribution among individuals of the same species and show presence/absence polymorphism . Recent reviews have noted on the exceptional genome plasticity and particularly the widespread presence of supernumerary chromosomes as a hallmark of pathogenic fungi [6, 8]. For this reason, they represent excellent model organisms to investigate eukaryotic genome evolution .
In some cases, the supernumerary chromosomes contain genes involved in pathogenicity and/or efficient host colonization like in Alternaria alternata , Fusarium solani (formerly known as Nectria haematococca) , and F. oxysporum f. sp. lycopersici . In other species, the role of supernumerary chromosomes is less clear as they do not show obvious pathogenicity related functions, such as Zymoseptoria tritici (formerly known as Mycosphaerella graminicola) . Regardless, in all cases described up to now supernumerary chromosomes differ markedly from the core chromosomes in characteristics such as gene content, codon usage and distribution of transposable elements.
The specific reasons for the these differences are unclear. A popular hypothesis proposes a different evolutionary origin for supernumerary chromosomes, that are subsequently acquired by horizontal chromosome transfer [11, 13]. Different evolutionary pressure on supernumerary chromosomes has also been proposed as an explanation for the detected differences [5, 6]. It has been argued that supernumerary chromosomes represent extreme cases of genome compartmentalization, as was demonstrated within the core chromosomes of Fusarium graminearum and Leptosphaeria maculans [14, 15]. These compartments may serve as evolutionary cradles, enriched in genes such as secondary metabolite clusters, often transcriptionally silent and only expressed under specific conditions. In accordance, degeneration from the core genome has been proposed as a potential origin of supernumerary chromosomes .
Transposable elements (TE) play an important role in fungal genome diversity and the evolutionary success of some pathogens [2, 16]. Examples are the vast differences in genome sizes of Fusarium and Phytophthora species [11, 17], and the shaping of pathogenicity in L. maculans and Pyrenophora tritici-repentis [15, 18]. The link between TE-mediated genome expansion and the evolution of virulence factors has been reviewed extensively . The possible deleterious effects of mobilization of TEs include gene disruption and intra- or inter-element recombination, potentially leading to gene loss. Fungi have evolved a specific genome defense mechanism against repetitive DNA, repeat-induced point mutation (RIP), that efficiently inactivates TEs by introducing cytosine to thymidine mutations . However, this process does not discriminate between TE proliferation and gene duplications, and therefore the near-absence of paralogs has been found to be a hallmark of a RIP-active species, e.g. in Fusarium graminearum  and Neurospora crassa . This finding has been termed the evolutionary cost of genome defense . RIP functions on repetitive sequences with greater than 80 % identity and exceeding +/- 800 bp in length .
The RIP process occurs only during the di-karyotic pre-meiotic phase and is therefore intricately associated with sexual cycle . Meiosis in fungi is partially regulated by the genes occupying the mating type locus. In heterothallic fungal species the locus is occupied by either the MAT1-1 or the MAT1-2 idiomorph, and isolates of these species require partners of the opposing mating type to enter into meiosis. For many species, a “cryptic” sexual cycle is presumed to occur in the field, that has never been witnessed nor simulated in the lab . The presumption of active meiosis becomes substantiated when markers for a sexual lifestyle are considered, including recombination , RIP , distributions of the mating idiomorphs in the population , and functional constraint on the genes implicated for meiosis .
Meiosis is one of the drivers of diversity in length and number of supernumerary chromosomes in fungi. It has been shown that during meiosis a process called nondisjunction is responsible for the loss of these chromosomes in Z. tritici offspring, even if both parents contained the supernumerary chromosome . Importantly, the offspring of these crosses are viable, underlining the conditionally dispensable nature of this part of the genome. The birth of a new supernumerary chromosome has been experimentally shown to occur through fusion of sister chromatids during meiosis, followed by breakage-fusion-bridge cycles . How the presence of supernumerary chromosomes influences the fate of the core chromosomes and whether a crosstalk between the two genome complements exists, has not been investigated.
The genus Fusarium comprises many agriculturally and medically important pathogens . As described above, species of this genus contain the hallmarks of fungal genome plasticity such as supernumerary chromosomes and compartmentalization of the core chromosomes. The Fusarium Head Blight disease on wheat and other small-grain cereals is caused by a number of species often co-occurring on the ear. Within this complex, Fusarium poae has been increasingly detected in a number of countries [32, 33]. Individuals of this species were shown to contain a highly variable set of supernumerary chromosomes . The genome of the related species F. graminearum is one of the best assembled fungal genomes and therefore represents an excellent model for comparative genome biology studies, moreover it lacks supernumerary chromosomes . In this study, we used SMRT sequencing to provide a high quality genome assembly of F. poae, and by a comparison with F. graminearum the cause of the sharp differences between the core and supernumerary genome was determined. This has led to a genetic crosstalk between the core and supernumerary genome, and the role of supernumerary chromosomes as evolutionary cradles that accelerate chromosomal and organismal evolution in fungi, was reaffirmed.
The genome is composed of a core and supernumerary part
The genome of F. poae isolate 2516 was assembled from the single-molecule real-time (SMRT) reads using a combined approach, based on macrosynteny with related Fusarium species as well as support from different assemblies using different parameters. In this combined approach, the largest contigs from one SMRT assembly (assembly A) were queried against two SMRT assemblies based on different parameters (see Methods). By detecting long collinear stretches at the edges of these contigs, they could be manually merged to longer contiguous sequences and finally to chromosomes. Correct merging of contigs was verified by inspecting the mapping of SMRT reads. Finally a merged assembly of 4 chromosomes was obtained, originating from 9 contigs (two, three, three and one respectively), accounting for a total 38.13 Mb of sequence. The merged assembly of 4 chromosomes was supplemented with the remaining contigs and the degenerate unitigs from assembly A, and only the 172 non-redundant contigs were kept (see Methods).
Comparison of the SMRT and HiSeq assemblies of isolate 2516. The statistics for the SMRT assembly were extracted from the final version of the assembly: four core chromosomes and 172 supernumerary contigs
Number of contigs
Total sequence length (bp)
Average sequence length (bp)
Minimum sequence length (bp)
Maximum sequence length (bp)
N50 sequence index (# of contigs)
N50 sequence length (bp)
The base quality of the assembly of the four core chromosomes was checked by mapping the HiSeq reads of isolate 2516 to the reference assembly. Over these 38.13 Mb, only one single nucleotide polymorphism (SNP) was detected between the HiSeq and SMRT reads. Two hundred and twenty-two variants were detected in homopolymeric stretches of nucleotides and low complexity (low GC%) regions (219 and 3 respectively). For these variants, read mapping was inconclusive for both HiSeq and SMRT reads, with both read batches giving support for different nucleotide calls.
The remaining 172 contigs (8.18 Mb of the total 46.3 Mb) do not show any synteny with sister species F. graminearum. These 172 contigs contain eight copies of the ZIT1 transposable element (TE) described earlier as a specific marker for supernumerary chromosomes in F. poae ; this element was not found on the four core chromosomes. Not only ZIT1, but all TEs show an unequal distribution between the core chromosomes and supernumerary contigs, which is described in detail further below. The most striking difference in TE distribution comes from a Miniature Inverted–Repeat Transposable Element (MITE) that is the most abundant repetitive element in the genome. All 712 copies (with expect value < e-10) were found dispersed over the core chromosomes. In sharp contrast no MITE was found on the 172 supernumerary contigs. Its positional conservation among F. poae isolates was investigated, and for all four isolates the localization of the vast majority of MITE instances (97.7-99.8 %) was identical to at least one other isolate.
Taken together, the 172 contigs are likely to make up one or more supernumerary chromosomes, and they are designated as the “supernumerary genome” for the purpose of this study. The 8.18 Mb that the 172 contigs contain is likely a slight overestimation, as in some instances the end of one contig is collinear and identical to the start of the contig that follows in the assembly, potentially indicating a (partial) double assembly (see further below for example). Besides the four core chromosomes and 172 supernumerary contigs, the assembly contains the mitochondrial genome of 138 kb, and three mitochondrial plasmids. These plasmids were confirmed to be mitochondrial as their ORFs encode a reverse transcriptase only with the mold mitochondrial genetic code.
A high quality machine annotation
General features of the machine annotation of F. poae 2516. These were compared to the features of the published annotation of F. graminearum PH-1 
Genome size (bp)
46 309 701
38 129 297
8 180 404
37 958 956
# of genes
Mean gene density (per Mb)
Median gene length (bp)
Median intron length (bp)
BUSCO analyses of F. poae and F. graminearum
Fragmented BUSCOs are proteins that are only partially recovered . These were analyzed manually (Additional file 1). Three of the fragmented BUSCOs were shared between F. poae and F. graminearum, and examination of the RNAseq data did not provide conclusive evidence that the genes are miss-annotated. The remaining four and three proteins of F. poae and F. graminearum represent gene models that are likely to be miss-annotated in these species. In all four cases of F. poae and in one case of F. graminearum, a hybrid gene model was built from two separate genes. The remaining two F. graminearum gene models respectively lack two exons and contain two exons in excess. The validation of all fragmented BUSCOs with RNAseq data can be found in Additional file 1. For F. graminearum an RNAseq data set described before was used . The BUSCO analysis suggests that the annotation of F. poae 2516 did not miss any conserved genes, and within the conserved genes, <0.5 % is miss-annotated.
The ingredients for meiosis and RIP are present in the genome
RIP only functions during the sexual cycle, which has not been definitively shown in F. poae. Therefore the conservation of the necessary ingredients for meiosis was investigated for isolate 2516. The MAT1-1 locus was extracted from the assembly, and its architecture is presented in Additional file 2. As all four isolates in this study have the MAT1-1 mating type, the architecture of the MAT1-2 locus could not be investigated. The number, order and direction of the genes occupying the MAT1-1 locus is identical to that in other Fusarium species . The MAT1-1-1, MAT1-1-2 and MAT1-1-3 genes have a predicted ORF with high similarity to those found in related species (85, 86 and 92 % similarity to proteins from F. graminearum). The MAT1-1-1 gene was previously identified for F. poae  and has 99 % protein similarity with the gene model in this study. Transcription of MAT1-1 idiomorph was noted in the RNAseq data, and the predicted splice forms lead to functional proteins (Additional file 2). A collection of 60 isolates was screened for the presence of MAT1-1 and MAT1-2 and both idiomorphs were detected, albeit in heavily skewed distribution (Additional file 3).
The KEGG pathway for meiosis in F. graminearum (fgr04113) was consulted to identify proteins involved in a putative sexual cycle. The conservation of this ‘meiotic toolbox’ was investigated in F. poae. All fifty-one entries in fgr01443 give best reciprocal protein hits with F. poae at expect values below 10-150, indicating that all ingredients of the meiotic toolbox are present in F. poae (Additional file 4).
Only a limited number of genes have been identified that belong to the machinery for RIP. A homolog of the rid (RIP defective) gene, shown to be vital for RIP, of F. graminearum is present, intact and transcribed in F. poae (Additional file 2). However, the expected intron was not spliced in the RNAseq data and the splice variant that was observed encodes a protein with a premature termination of translation.
Distribution of transposable elements differs markedly between core and supernumerary genome
Classification and key characteristics of TE families in the genome of F. poae 2516. Elements below the length threshold for RIP are not included (MITE, ZIT1). Repetitive elements such as the rDNA tandem and two families of telomere linked RecQ helicases are not included. Nomenclature of TEs is as recommended in literature . R retrotransposon, D DNA transposon, L long terminal repeat (LTR), T terminal inverted repeat (TIR), G Gypsy, C Copia, F Fot1/Pogo, T Tc1/mariner, M Mutator, A hAT, x unknown. n/a designates instances where a TIR/LTR could not be detected for a specific element
Transcription and splicing of the predicted introns were noted for many intact elements. Single nucleotide polymorphisms (SNPs) in the RNAseq reads allowed for the specific copy/copies of TEs that were transcribed, to be identified. These are included as TE prototypes in Additional file 5 and annotated in Additional file 6. The functional and structural features of every TE were used for the classification into superfamilies (Table 4). Additional file 6 also contains phylogenetic trees of every element for which a protein coding sequence could be determined. In most cases, TE phylogeny lines up well with species phylogeny. Exceptions are the TEs DTA_Nymeria and DTM_Hop7, that show higher similarity to elements from unrelated fungi than to elements from related Fusarium species.
Unbalanced RIP between core and supernumerary genome
RIPcal was used to analyze the occurrence of RIPped copies from the 33 transposable element families, separately for the core and supernumerary genome. The results are presented in Additional file 7. On the core genome, there are 13 families that show RIPcal patterns that are typical for RIP (dominance of CpA → TpA mutations; red trace). On the supernumerary genome, there are no such instances. The alignments that were used to perform RIPcal analysis were then manually examined in an attempt to quantify the number of RIPped copies per family, on the core and the supernumerary genome. Indeed, RIPped copies of TEs were only detected on the core genome (Table 4).
In a few cases, low complexity regions on the supernumerary genome resemble RIP of intact elements. These also contain most transversions when compared to genuinely RIPped copies (Additional file 8).
Transposable elements copy number is dynamic between isolates of the same species
Interestingly, RIP of certain elements seems isolate-specific, such as DTA_RLT1 in isolate bfb0173 and RLC_Ghost in isolate 2516. A process similar to the loss of DTM_Hop7 and RLG_Maggy in isolate bfb0173 may have occurred species-wide, as RIPped elements in isolate 2516 were detected of up to 14 families that no longer contain any intact copies in this isolate, or any other isolate in this study. For a retrotransposon of the Gypsy family, RIPped copies are present in all isolates, but only isolate 2548 contains intact copies.
Localization and divergence of transposable elements differs between the core and supernumerary genome
The localization and divergence of the intact TEs was investigated. One hundred and thirty-five intact TEs are present on the core chromosomes of isolate 2516. Figure 3 shows that elements of these families are often, but not always, also present in multiple copies in the genomes of the other isolates in this study. However, read mapping shows that none of the 135 elements on the core chromosomes of isolate 2516 are present in the same location in isolates 2548, 7555 and bfb0173 (as exemplified in Additional file 9). The integration of these elements therefore seems to have happened recently.
The core genome is invaded by transposable elements from the supernumerary genome
On the supernumerary genome, TEs are more divergent than TEs on the core chromosomes and some are present at identical sites in at least two isolates. We therefore wanted to test whether TEs in the core genome may originate from the supernumerary genome. This directionality was best illustrated by one element of RLG_Skippy located on the short arm of chromosome 3. This TE has recently integrated into the genome of isolate 2516 and contains 23 SNPs compared to the other copies on the core chromosomes of this isolate. These 23 SNPs, together with one additional SNP, are also present in a copy on the supernumerary genome, that is at the exact same location for all isolates and therefore is an ancestral insertion.
The recent integration of TEs into the core genome of isolate 2516 has sometimes occurred within the coding region of genes. The environment of all 135 recent integrations in the core genome was investigated (Additional file 12). Ten instances were found where integration disrupted a gene. Notably, DTF_Drogon integrations account for nine of these. Remarkable differences can be detected between the environments of the different TE families. DTF_Fot5-A elements consistently integrate within RIPped or low complexity (low GC%) environments, while DTF_Drogon elements have integrated within regions of average GC%.
An extreme case of core genome invasion is found near the telomere of the long arm of chromosome 3. Two sequence blocks do not show any synteny with F. graminearum, with coordinates 115073-319336 (204 kb) and 883738-1348064 (464 kb). Analysis of the flanking sequences of these two regions shows that they are continuous in isolates 2548, 7555 and bfb0173. Therefore, these regions represent translocations of supernumerary sequence to the core genome of isolate 2516. All parameters that were used to compare the core and supernumerary genome in this study, support the classification of these sequence blocks as supernumerary sequence. For the purpose of this study they have been regarded as part of the “supernumerary genome”. The 204 kb insertion is an underestimate as it holds a 5000 N placeholder, where presumably one or more of the 172 supernumerary contigs belong.
These insertions into the core genome, both of single TEs and whole blocks of supernumerary sequence, may have large implications for the biology of the organism such as respectively gene disruptions and hampered meiotic alignment. We investigated how common the two supernumerary sequence translocations are in a population of 60 F. poae isolates and found that seven isolates contain the first insertion, closest to the telomere. These seven isolates were isolated from three different locations in Belgium. Three of these seven isolates also contain the second insertion, at 883738 bp into chromosome 3 (Additional file 3). Isolate 2516 is one of these three isolates. Isolates 2548, 7555 and bfb0173 were confirmed not to have any of the insertions.
The supernumerary genome is a refuge for gene duplications
The absence of paralogs is a hallmark of a RIP-active species . In a blastp of all proteins encoded by genes on the core chromosomes against themselves, no hits with identity >80 % and length above the RIP length threshold can be found. This confirms that the core chromosomes are subjected to RIP. Additional file 13 shows the blastp output of all proteins encoded on the supernumerary genome queried against those encoded in the core genome. Many hits well above the RIP length threshold show high protein identity (>80 %) and represent genuine gene duplications that have not been inactivated by RIP. A total of 44 genes on the core chromosomes have one or more duplicates in the supernumerary genome, totaling 104 hits on the latter. Figure 2 visualizes these gene duplications as lines connecting both the gene on one of the four core chromosomes and its paralog(s) on the supernumerary genome.
Additional file 13 lists the functional annotation of these duplicated genes. Notable instances include the transcription factor EBR1 [41, 42], key component of the RNA silencing pathway Dicer2  and the secondary metabolite backbone gene PKS8 . In the case of EBR1, up to 11 duplications were found. To ascertain that these are not artefacts from the assembly, the duplicated sequences and their 500 bp flanking regions were aligned. Assembly artefacts are identified by sequence alignments with nearly 100 % sequence identity across the entire region. Five of the 104 duplications were identified as likely double assemblies. The remaining 99 instances are likely genuine gene duplications (Additional file 13). Additional file 14 shows an example of a potential assembly artefact as well as an example of an expected duplicate gene.
The presence of the duplicated genes was assessed between the different isolates of F. poae in this study. Additional file 13 shows that some duplications are unique for isolate 2516, while others are not cumulative over the different isolates. Indeed, isolates 2548, 7555 and bfb0173 each share unique duplications with isolate 2516 that are absent in the other isolates.
Chromosomes that vary in both size and number, and have an uneven distribution among individuals of the same species, have been described in animals , plants  and fungi . Throughout these kingdoms they are identified as supernumerary, accessory, dispensable or B chromosomes, in contrast to the core or A chromosomes. These supernumerary chromosomes show distinct features compared to the core chromosomes: they can be high in repeats and transposable elements , have different gene density and function  and/or GC-content , are epigenetically very dissimilar , can be transmitted by horizontal transfer  and are unstable in meiosis . This sharp contrast between the two sets of chromosomes can be explained by different evolutionary pressure, different origins, or a combination of both [6, 13]. The functions of the supernumerary chromosomes are unclear. Speculations range from selfish DNA fragments without benefit to the host  to components that are critical for pathogenicity and survival . As they occur in the same nucleus, this raises questions how the different sub-genomes are managed within the organism and if there are potential conflicts and interactions. The dynamics between the core genome of an organism and its extra-chromosomal DNA have been studied for plasmids in bacteria , and mitochondrial DNA insertions into the nuclear genome of many eukaryotes . In this study we aimed to provide a contiguous genome assembly of a fungal pathogen that contains supernumerary chromosomes, and to compare it with a sister species that contains none. This allowed us to determine differences in genome biology as the causal agent for the differences between the core and supernumerary genomes, and to observe a significant crosstalk between them as a result.
The assembly using SMRT long reads allowed the division of the genome into a core and supernumerary part, a feature that was not possible with the assembly using short reads due to the occurrence of highly repetitive DNA. The advantages of SMRT sequencing for fungal genome assembly was recently demonstrated for Verticillium dahliae . The core chromosomes of F. poae showed a high level of macrosynteny with F. graminearum and cover the entire F. graminearum sequence complement. They showed characteristics of chromosomes under the control of RIP, such as many inactivated TEs and no gene paralogs with high sequence identity. The opposite holds for the supernumerary chromosomes, and the absence of RIP on the supernumerary genome is responsible for the differences between the two core and supernumerary genome in F. poae. Indeed, on the supernumerary chromosomes, no RIPped TE copies are found, and many duplicated genes are present. This is the first time a definitive causal agent is identified for the sharp contrast between the core and supernumerary chromosomes in a fungal pathogen.
The different rules that govern the core and supernumerary chromosomes, lead to a genetic crosstalk between them. We found many cases of exchange of genetic material between the core and supernumerary genomes. Genes from the core chromosomes are duplicated to the supernumerary genome, where some of these genes have undergone further copy number expansion. Vice versa, transposable elements originating from the supernumerary genome, have integrated into the core chromosomes and, in some instances, have led to gene disruptions. Most drastically, large sequence blocks (>200 kb) have been translocated from the supernumerary genome to the core chromosomes. Translocation of whole regions from supernumerary chromosomes to core chromosomes is not restricted to F. poae. A region of core chromosome 1 of F. oxysporum f. sp. lycopersici has all the characteristics of supernumerary sequence [8, 11]. It has been shown that this region is highly syntenic with one of its supernumerary chromosomes . This situation may have arisen from ancient translocation followed by chromosome gain, or by duplication and integration of the supernumerary sequence. Our results show that unique events such as large insertions into the core chromosomes give rise to novel genotypes in F. poae, which may be able to rapidly spread as they were recovered from three different locations. Interestingly, both in F. poae and F. oxysporum f. sp. lycopersici, insertion of blocks of supernumerary sequence has occurred close to a telomere of core chromosomes, which supports the finding that core chromosomes in Fusarium species may be divided into distinct chromosomal regions on both a structural and functional level [14, 53].
The evolutionary advantage of a genome region not burdened by RIP defense is evident from the many gene duplications occurring specifically on that region. The primary versus secondary metabolism master regulator EBR1 is present as a single copy gene in F. graminearum . In F. oxysporum f. sp. lycopersici it has undergone gene duplication on the supernumerary chromosomes , which seemingly has independently occurred on the supernumerary chromosomes of F. poae. In Metarhizium, it has been hypothesized that an arrest of RIP was instrumental in the evolution to generalist infection agents , and the advantages of a temporary RIP relaxation or arrest may include accelerated evolution and divergence between related species . The mechanism that gave rise to the existence of supernumerary chromosomes in F. poae is unknown. In Z. tritici, it has been suggested that supernumerary chromosomes may have originated from core chromosomes and subsequently degenerated and evolved separately . A recent large scale duplication seems unlikely for F. poae, as an ancient MITE that characterizes the core chromosomes is completely absent from the supernumerary chromosomes. The contribution of horizontal chromosome transfer remains to be investigated, but has likely contributed to the supernumerary chromosome diversity in other species [11, 55].
Specific exclusion of duplicated genes from RIP has been observed before in F. solani (formerly known as N. haematococca), a species distantly related to F. poae . Similar to those in F. poae, its supernumerary chromosomes are rich in TEs and gene duplications. It was experimentally determined that progeny, of which one of the parents contained two copies of the hph marker gene, contained both an intact and a RIPped copy of that gene , contrary to the standard modus operandi of RIP wherein all copies are inactivated. Apparently, a region in the genome of N. haematococca is excluded from RIP, but it is not known whether this coincides with the supernumerary chromosome(s). The regional variability of RIP extends beyond supernumerary chromosomes however. Nucleolus organizer regions (NORs) contain the rDNA tandem repeats in many fungal species, and within this region they are protected from RIP . Either this is an active form of protection, which may be similar to what is happening on the supernumerary chromosomes of F. poae, or rDNA is mutated by RIP and mutated copies subsequently undergo reduced intrachromosomal recombination to give rise to a full-length rDNA tandem during meiosis .
It is thought that RIP functions after plasmogamy but before the final pre-meiotic DNA replication and karyogamy. It works multiple times during the rounds of nuclear division that occur at this point, presumably during G1 or near the replication fork during the S phase . It functions only in the nucleus, or nuclei, that contain(s) DNA duplications and does so on a single DNA strand. Mis-pairing of duplicated DNA has been hypothesized to deliver the substrate for RIP . Why supernumerary chromosomes in F. poae escape RIP is unknown. While physical alignment of duplicated copies is presumed to be important, the exact search mechanism for homology is unknown. Clutterbuck  proposed two hypotheses for the function of RIP that implicate the temporal or spatial proximity of haploid nuclei in dikaryotic cells, where RIP acts, to diploid cells undergoing meiotic pairing. This pairing was previously shown to be hampered for supernumerary chromosomes through their high variability . The splice form of the rid gene, shown to be vital for RIP, that was detected in this study does not lead to a functional protein. One explanation for this may be that the gene model constitutes a case of crucial alternative splicing, and the rid gene is only correctly spliced in the pre-meiotic phase.
Supernumerary or extra-chromosomal structures are considered to be evolutionary cradles for pathogenicity in viruses, bacteria and fungi . We showed that they shelter TEs and gene duplications in a fungal plant pathogen, and that large sequence blocks may translocate to the core genome, with a profound effect on genome biology. The “living apart together” dynamic between the core and supernumerary genome is explained by different rules that govern both genome compartments. Future investigation will be able to further characterize the core versus supernumerary chromosomes in parameters that have been shown to be relevant in other pathogens such as histone modification , transcription  and gene content .
The genome biology of the related plant pathogenic fungi F. poae and F. graminearum primarily differs by the presence of supernumerary chromosomes in F. poae, not present in F. graminearum. These supernumerary chromosomes differ markedly from the core chromosomes. We provided a high quality genome assembly for F. poae and determined the cause of these differences to be RIP, a mutational defense mechanism against TEs which functions on the core chromosomes but not the supernumerary chromosomes. This has led to a dynamic crosstalk between the core and supernumerary genome, which significantly affects chromosomal and organismal evolution.
Isolates used for whole genome sequencing
Nucleic acid manipulation, library preparation and sequencing
Detailed information on growth conditions and nucleic acid manipulation can be found in Additional file 15. In short, for HiSeq sequencing, DNA was extracted from isolates bfb0173, 2516, 2548 and 7555. DNA from isolate bfb0173 was used for random sheared shotgun library preparation using the NEXTflex ChIP-seq Library prep kit with adaptations for low DNA input according to the manufacturer’s instructions (Bioscientific). The library was loaded as (part of) one lane of an Illumina Paired End flowcell for cluster generation using a cBot. Sequencing was performed on an Illumina HiSeq2000 instrument using 101, 7 and 101 flow cycles for forward, index and reverse reads respectively. De-multiplexing of resulting data was carried out using Casava 1.8. Shotgun libraries were made for isolates 2516, 2548 and 7555 using the Illumina TruSeq LT DNA sample prep kit according to the manufacturer’s instructions (Illumina). Libraries were then pooled equimolarly and loaded on one flowcell lane for 2x100 nt paired end sequencing on an Illumina HiSeq2000 platform as described above.
For SMRT sequencing, DNA was isolated from isolate 2516. Genomic DNA was extracted with the Wizard Magnetic DNA Purification System for Food (Promega) according to the manufacturer’s instructions. Twenty μg DNA was used for a large insert (10 kb) library prep according to the manufacturer’s instructions (Pacific Biosciences) with small adaptations (Additional file 16). After library prep and SMRT bell adapter ligation, the SMRT bells were size selected with a 7000 bp minimum cutoff. Sequencing was done on a Pacbio RS II system using one cell per well, C4 chemistry and 240 min movie time. A total of 16 SMRT cells was run.
RNA samples were prepared from isolate 2516 grown under different conditions, designed to have maximal numbers of genes expressed. Details can be found in Additional file 16. In short, samples were designed to favor different metabolic stages. Primary metabolism was simulated in a rich medium. Secondary metabolism was induced in five distinct conditions, namely trichothecene biosynthesis induction, fungicide application, N-starvation, C-starvation, and conidiation under UV. RNA was extracted from the rich medium and the 5 “stress-inducing” treatments with TRIzol (Life Technologies). Subsequently, the crude RNA was purified with the RNA cleanup protocol included in the RNeasy Plant Mini kit (Qiagen) according to the manufacturer’s instructions. RNAseq library preparation was performed using the Illumina TruSeq total RNA sample preparation kit and guidelines. Libraries were pooled equimolarly and loaded on one Illumina HiSeq2000 flowcell lane for 2x100 nt paired end sequencing as described above.
From SMRT Portal version 18.104.22.168893, the HGAP2 (Hierarchical Genome Assembly Process) was initiated using data from 16 SMRT cells. Raw reads were filtered on Q = 0.83, polymerase read length >1000 bp and subread length >1200 bp. The seed read length for the error correction procedure was manually set to 6 kb. After the error correction step, the data was filtered for reads >9 kb. With this dataset the assembly was performed with Celera, using the default settings provided by HGAP2. The contigs from this assembly (assembly A) were taken as the basis for a merged assembly, for which the contigs were supplemented with those derived from two additional automatic HGAP2 assemblies (assemblies B and C). One assembly was performed with 10 SMRT cells using default settings, and one with 16 SMRT cells using only HQ input data (read quality > = 0.85, polymerase and subread length >4000 bp).
The largest contigs, corresponding to a major part of a hypothetical chromosome were used as queries for blastn searches against the other contigs in the assemblies. Contigs were joined on the basis of their colinearity (usually excluding the very end of one contig and the very beginning of another, where the assembler presumably stalled or followed a wrong seed for a particular assembly) and their macrosynteny with other Fusarium species. Rightful joining of contigs was evaluated by mapping SMRT long reads. The resulting merged assembly of four chromosomes was polished using Quiver (SMRT Portal resequencing protocol) for 2 times, using HQ reads (read quality > = 0.80, polymerase and subread length >3000 bp). The merged assembly was supplemented with the remaining contigs and the degenerate unitigs from assembly A. Nine contigs that only contained rDNA tandem repeats, and 66 contigs that contained mitochondrial sequence, were removed. The remaining contigs were added to the four-chromosome assembly, and the entire assembly was error corrected in 1 pass using quiver (read quality > = 0.84, polymerase and subread length >1000 bp). After this quiver run, 13 contigs were removed from the assembly that had overall base quality scores close to zero, compared to an average base quality of 50 for the rest of the assembly.
The mitochondrial genome was assembled with GRABb using standard settings and with the PH-1 mitochondrion (NCBI accession HG970331.1) as bait (Brankovics et al., submitted). One mitochondrial plasmid was assembled in SMRT assembly A. Two additional mitochondrial plasmids were taken from the HiSeq assembly (see below), that were not present in the SMRT assembly. This may result from the fragment size selection that was performed, as the plasmids are <3 kb in size. The final assembly therefore contains four chromosomes, 172 unplaced contigs, one mitochondrial genome and three mitochondrial plasmids. For the four chromosomes, the error rate of the SMRT assembly was checked by mapping the HiSeq reads of isolate 2516 to the SMRT assembly with CLC Genomics Workbench 7.5 (length and similarity fraction = 0.8). For the four core chromosomes, basic variant detection was run with minimum coverage = 50, minimum count = 10 and minimum frequency = 70 %, and other settings at standard value. A de novo assembly of Illumina HiSeq reads for isolate 2516 was perfomed with CLC Genomics Workbench 7.5 using standard settings.
Annotation of the reference genome
The Illumina HiSeq RNAseq paired-end reads were cleaned and trimmed using Trimmomatic . Tophat2  was used to map the trimmed reads to the SMRT assembly of isolate 2516. The mapping results were used in the genome annotation pipeline BRAKER1  for training GeneMark  and Augustus . BRAKER1 uses the introns parsed from the TopHat2 mapping as extrinsic evidence for the final gene models predicted by Augustus. The annotation is outputted as a GFF file with genes, introns, exons and the protein sequences predicted to be encoded by these genes. A short Python script was used to extract the protein sequences from the GFF.
Repeat identification, localization, structural and functional characterization
RepeatModeler  was run on the genome of isolate 2516 with standard settings. RepeatModeler output was manually curated to obtain complete elements. These elements were then subjected to functional and structural characterization. When possible, terminal inverted repeats (TIR) and long terminal repeats (LTR) were identified. Bowtie2/TopHat2 read mapping as well as related NCBI accessions were consulted to find intron/exon boundaries. The translations of the predicted ORF for every TE were used as blastp queries against the NCBI non-redundant protein sequences database. The 15 best hits were aligned with the TE query using ClustalO  implemented in CLC Genomics Workbench 7.5. The resulting neighbor-joining phylogenetic trees are included in the TE data sheets. The elements were divided into superfamilies based on their domain similarities to described TEs. For RIP analysis, a copy coding for a functional protein (determined as described above) was used as a query for blastn in CLC Genomics Workbench 7.5 against the core and supernumerary genomes separately, at expect value < 1e-10. All hits were aligned to the query with ClustalO , and the alignments were fed directly to RIPcal for alignment-based RIP analysis with the query as the defined reference .
Analysis of transposable element integration sites
Blastn was used to obtain genomic coordinates of all intact and RIPped copies of TEs. Using the getfasta utility of BEDtools  these hits, including their flanking regions, were extracted. HiSeq reads from all isolates in this study were subsequently mapped to the extracted reference sequences at high stringency (minimum length fraction 0.95, minimum similarity fraction 0.95) in CLC Genomics Workbench 7.5. Results were manually inspected to find identical genomic environments. Elements with read support for only one flank were also considered to be identically inserted.
For synteny of the MITE, a prototype of the element (see Additional file 5) was used as a blastn query against the entire genome (expect value < 1e-10), hits were extracted including 500 bp upstream and downstream flanking sequence with the getfasta utility of BEDtools . This was done for the assembly with long reads of isolate 2516 and the short read de novo assemblies of isolates 2548, 7555 and bfb0173. For every isolate, the resulting sequence list was queried against the genome assemblies of the other isolates with blastn, and the matches longer than 640 bp were counted (indicating instances where the localization of the MITE coincides between isolates: at least one flank of 500 bp and the 140 bp element are shared).
Divergence estimates of TE copies
Intact (not RIPped) copies were extracted as described above from the core and unplaced sequence separately. Only families containing five or more copies were retained. ClustalO alignments  were fed to PhyML  and maximum-likelihood phylogenies were built for every family with settings retrieved from literature . Specifically, a neighborhood joining tree was used as starting tree, the transition/transversion ratio was 4, the HKY85 evolution model was used and distribution parameters were allowed to optimize. In the resulting phylogenies, terminal branch lengths represent the relative age of every separate element. These branch lengths were extracted from the Newick files with Newick Utilities . Using the substitution rate determined for protein-coding genes in fungi (1.05 * 10-9 ), divergence time estimates were calculated from the branch lengths. These were then visualized as boxplots using SPSS.
Paralogs and gene duplications
Blastp was used to find paralogs for all proteins (initial blastp parameters at expect value < 1e-5). Results were filtered to > 80 % identity and a length above the RIP threshold (+/- 800 nt or 266 amino acids). When the proteins encoded by genes on the core genome were queried against themselves, no hits other than self-hits were found. When the proteins encoded by genes in the supernumerary genome were queried against those from the core genome, 104 hits were obtained. These were formatted for Circos visualization. Using the gff2sequence tool  the 104 genes in the supernumerary genome and their 44 paralogs on the core chromosomes were extracted, including 500 bp up- and downstream. These sequences were aligned all-vs-all with Smith-Waterman using a python application based on PaSWAS , which produced local alignments in SAM  format.
Circos  was used for circular genome visualization. Locations of TEs were extracted from blastn output (expect value < 1e-10). Gene duplications above the RIP threshold were parsed from blastp output.
To determine the MAT1-1 and MAT1-2 distribution in the F. poae population, primer pairs POA-1-F/POA-1-R and POA-2-F/POA-2-R were used . For the two insertions of supernumerary sequence into the core chromosomes, primers were designed flanking insertion site as well as covering the extremes of the inserted block (visualized in Additional file 16). All primers used in this study can be found in Additional file 17. PCR reactions were performed as described earlier . Gel electrophoresis of the diagnostic PCRs is shown in Additional file 17.
Comparative genomics intra- and inter-species
To estimate the number of intact copies for every family in the isolates that were sequenced with only Illumina technology, reads were mapped to the curated library of repeats (see above), and the resulting coverage was normalized against the mean coverage of the single-copy genome for each isolate. To estimate the coverage across the four largest supernumerary contigs, HiSeq reads from every isolate were mapped with CLC Genomics Workbench 7.5 (length and similarity fraction = 0.8) to the reference assembly which was masked for TE with RepeatMasker. BAM files were processed with the coverage utility of BEDtools  to find the fraction of bases covered by reads in a 1 kb sliding window. For whole-genome alignment, the genome of the reference isolate was masked using the curated repeat library with RepeatMasker. The masked genome was aligned with the completed genome of F. graminearum PH-1 using MUMmer .
BUSCO, benchmarking universal single-copy orthologs; LTR, long terminal repeat; MITE, miniature inverted-repeat transposable element; Rid, RIP defective; RIP, repeat-induced point mutation; SMRT, single-molecule real-time; SNP, single nucleotide polymorphism; TE, transposable element; TIR, terminal inverted repeat
We kindly acknowledge Els Verstappen, Marga van Gent-Pelzer, Thamara Hesselink, Bas te Lintel-Hekkert and Boris Bekaert for excellent technical support.
AV is supported by a grant from the Ghent University Special Research Fund (grant number 01DI2314) and a travel grant from the Research Foundation - Flanders (grant number K147115N). TvdL and CW are supported by a grant from the Dutch Ministry of Economic Affairs, Agriculture and Innovation. The funding bodies had no role in the design of the study, in collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The datasets generated and analyzed during the current study have been deposited with links to BioProject accession number PRJNA319914 in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA319914).
AV, KA, CW and TvdL conceived and designed the study; AV performed genome analysis; ES coordinated sequencing, HvdG performed assembly; SW performed annotation; AV, SW, CW and TvdL analyzed the data; AV, KA, SW, HvdG, ES, MH, SDS, GH, CW and TvdL wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Coghlan A, Eichler EE, Oliver SG, Paterson AH, Stein L. Chromosome evolution in eukaryotes: a multi-kingdom perspective. Trends Genet. 2005;21:673–82.View ArticlePubMedGoogle Scholar
- Stukenbrock EH, Croll D. The evolving fungal genome. Fungal Biol Rev. 2014;28:1–12.View ArticleGoogle Scholar
- Vanneste K, Baele G, Maere S, Van de Peer Y. Analysis of 41 plant genomes supports a wave of successful genome duplications in association with the Cretaceous–Paleogene boundary. Genome Res. 2014;24:1334–47.View ArticlePubMedPubMed CentralGoogle Scholar
- Stukenbrock EH. Evolution, selection and isolation: a genomic view of speciation in fungal plant pathogens. New Phytol. 2013;199:895–907.View ArticlePubMedGoogle Scholar
- Sperschneider J, Gardiner DM, Thatcher LF, Lyons R, Singh KB, Manners JM et al. Genome-wide analysis in three Fusarium pathogens identifies rapidly evolving chromosomes and genes associated with pathogenicity. Genome Biol Evol. 2015:evv092.Google Scholar
- Croll D, McDonald BA. The accessory genome as a cradle for adaptive evolution in pathogens. PLoS Pathog. 2012;8(4):e1002608.View ArticlePubMedPubMed CentralGoogle Scholar
- Dong S, Raffaele S, Kamoun S. The two-speed genomes of filamentous pathogens: waltz with plants. Curr Opin Genet Dev. 2015;35:57–65.View ArticlePubMedGoogle Scholar
- Galazka JM, Freitag M. Variability of chromosome structure in pathogenic fungi - of ‘ends and odds’. Curr Opin Microbiol. 2014;20:19–26.View ArticlePubMedGoogle Scholar
- Hatta R, Ito K, Hosaki Y, Tanaka T, Tanaka A, Yamamoto M, et al. A conditionally dispensable chromosome controls host-specific pathogenicity in the fungal plant pathogen Alternaria alternata. Genetics. 2002;161:59–70.PubMedPubMed CentralGoogle Scholar
- Coleman JJ, Rounsley SD, Rodriguez-Carres M, Kuo A, Wasmann CC, Grimwood J, et al. The genome of Nectria haematococca: contribution of supernumerary chromosomes to gene expansion. PLoS Genet. 2009;5(8):e1000618.View ArticlePubMedPubMed CentralGoogle Scholar
- Ma L-J, van der Does HC, Borkovich KA, Coleman JJ, Daboussi M-J, Di Pietro A, et al. Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium. Nature. 2010;464:367–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Goodwin SB, Ben M'Barek S, Dhillon B, Wittenberg AHJ, Crane CF, Hane JK, et al. Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis. PLoS Genet. 2011;7(6):e1002070.View ArticlePubMedPubMed CentralGoogle Scholar
- Mehrabi R, Bahkali AH, Abd-Elsalam KA, Moslem M, Ben M'Barek S, Gohari AM, et al. Horizontal gene and chromosome transfer in plant pathogenic fungi affecting host range. FEMS Microbiol Rev. 2011;35:542–54.View ArticlePubMedGoogle Scholar
- Zhao C, Waalwijk C, de Wit PJGM, Tang D, van der Lee T. Relocation of genes generates non-conserved chromosomal segments in Fusarium graminearum that show distinct and co-regulated gene expression patterns. BMC Genomics. 2014;15:191.View ArticlePubMedPubMed CentralGoogle Scholar
- Rouxel T, Grandaubert J, Hane JK, Hoede C, van de Wouw AP, Couloux A, et al. Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations. Nat Commun. 2011;2:202.View ArticlePubMedPubMed CentralGoogle Scholar
- Raffaele S, Kamoun S. Genome evolution in filamentous plant pathogens: why bigger can be better. Nat Rev Microbiol. 2012;10:417–30.PubMedGoogle Scholar
- Haas BJ, Kamoun S, Zody MC, Jiang RHY, Handsaker RE, Cano LM, et al. Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans. Nature. 2009;461:393–8.View ArticlePubMedGoogle Scholar
- Manning VA, Pandelova I, Dhillon B, Wilhelm LJ, Goodwin SB, Berlin AM, et al. Comparative genomics of a plant-pathogenic fungus, pyrenophora tritici-repentis, reveals transduplication and the impact of repeat elements on pathogenicity and population divergence. G3 (Bethesda). 2013;3:41–63.View ArticleGoogle Scholar
- Cambareri EB, Jensen BC, Schabtach E, Selker EU. Repeat-induced G-C to A-T mutations in Neurospora. Science. 1989;244:1571–5.View ArticlePubMedGoogle Scholar
- Cuomo CA, Gueldener U, Xu J-R, Trail F, Turgeon BG, Di Pietro A, et al. The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization. Science. 2007;317:1400–2.View ArticlePubMedGoogle Scholar
- Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003;422:859–68.View ArticlePubMedGoogle Scholar
- Galagan JE, Selker EU. RIP: the evolutionary cost of genome defense. Trends Genet. 2004;20:417–23.View ArticlePubMedGoogle Scholar
- Pomraning K. Characterization of Neurospora crassa and Fusarium graminearum Mutants Defective in Repeat-induced Point Mutation. Oregon State University; 2012. http://ir.library.oregonstate.edu/xmlui/bitstream/handle/1957/35582/PomraningKyleR2012.pdf?sequence=1.
- Turgeon BG. Application of mating type gene technology to problems in fungal biology. Annu Rev Phytopathol. 1998;36:115–37.View ArticlePubMedGoogle Scholar
- Tavanti A, Gow NAR, Maiden MCJ, Odds FC, Shaw DJ. Genetic evidence for recombination in Candida albicans based on haplotype analysis. Fungal Genet Biol. 2004;41:553–62.View ArticlePubMedGoogle Scholar
- Clutterbuck AJ. Genomic evidence of repeat-induced point mutation (RIP) in filamentous ascomycetes. Fungal Genet Biol. 2011;48:306–26.View ArticlePubMedGoogle Scholar
- Paoletti M, Rydholm C, Schwier EU, Anderson MJ, Szakacs G, Lutzoni F, et al. Evidence for sexuality in the opportunistic fungal pathogen Aspergillus fumigatus. Curr Biol. 2005;15:1242–8.View ArticlePubMedGoogle Scholar
- Ropars J, Dupont J, Fontanillas E, de la Vega RCR, Malagnac F, Coton M, et al. Sex in cheese: evidence for sexuality in the fungus Penicillium roqueforti. PLoS One. 2012;7(11):e49665.View ArticlePubMedPubMed CentralGoogle Scholar
- Wittenberg AHJ, van der Lee TAJ, Ben M'Barek S, Ware SB, Goodwin SB, Kilian A, et al. Meiosis drives extraordinary genome plasticity in the haploid fungal plant pathogen Mycosphaerella graminicola. PLoS One. 2009;4(6):e5863.View ArticlePubMedPubMed CentralGoogle Scholar
- Croll D, Zala M, McDonald BA. Breakage-fusion-bridge cycles and large insertions contribute to the rapid evolution of accessory chromosomes in a fungal pathogen. PLoS Genet. 2013;9(6):e1003567.View ArticlePubMedPubMed CentralGoogle Scholar
- O'Donnell K, Rooney AP, Proctor RH, Brown DW, McCormick SP, Ward TJ, et al. Phylogenetic analyses of RPB1 and RPB2 support a middle Cretaceous origin for a clade comprising all agriculturally and medically important fusaria. Fungal Genet Biol. 2013;52:20–31.View ArticlePubMedGoogle Scholar
- Xu XM, Parry DW, Nicholson P, Thomsett MA, Simpson D, Edwards SG, et al. Predominance and association of pathogenic fungi causing Fusarium ear blight in wheat in four European countries. Eur J Plant Pathol. 2005;112:143–54.View ArticleGoogle Scholar
- Vanheule A, Audenaert K, De Boevre M, Landschoot S, Bekaert B, Munaut F, et al. The compositional mosaic of Fusarium species and their mycotoxins in unprocessed cereals, food and feed products in Belgium. Int J Food Microbiol. 2014;181:28–36.View ArticlePubMedGoogle Scholar
- Fekete C, Hornok L. A repetitive DNA sequence associated with karyotype variability in Fusarium poae. Acta Phytopathol Entomol Hung. 1997;32:29–38.Google Scholar
- King R, Urban M, Hammond-Kosack MC, Hassani-Pak K, Hammond-Kosack KE. The completed genome sequence of the pathogenic ascomycete fungus Fusarium graminearum. BMC Genomics. 2015;16:544.View ArticlePubMedPubMed CentralGoogle Scholar
- Hoff KJ, Lange S, Lomsadze A, Borodovsky M, Stanke M. BRAKER1: Unsupervised RNA-Seq-based genome annotation with GeneMark-ET and AUGUSTUS. Bioinformatics. 2015:btv661.Google Scholar
- Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31:3210–2.View ArticlePubMedGoogle Scholar
- Zhao C, Waalwijk C, de Wit PJGM, Tang D, van der Lee T. RNA-Seq analysis reveals new gene models and alternative splicing in the fungal pathogen Fusarium graminearum. BMC Genomics. 2013;14:21.View ArticlePubMedPubMed CentralGoogle Scholar
- Ma L-J, Geiser DM, Proctor RH, Rooney AP, O'Donnell K, Trail F, et al. Fusarium Pathogenomics. Ann Rev Microbiol. 2013;67:399–416.View ArticleGoogle Scholar
- Kerenyi Z, Moretti A, Waalwijk C, Olah B, Hornok L. Mating type sequences in asexually reproducing Fusarium species. Appl Environ Microbiol. 2004;70:4419–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Jonkers W, Xayamongkhon H, Haas M, Olivain C, van der Does HC, Broz K, et al. EBR1 genomic expansion and its role in virulence of Fusarium species. Environ Microbiol. 2014;16:1982–2003.View ArticlePubMedGoogle Scholar
- Zhao C, Waalwijk C, de Wit PJGM, van der Lee T, Tang D. EBR1, a Novel Zn(2)Cys(6) transcription factor, affects virulence and apical dominance of the hyphal tip in Fusarium graminearum. Mol Plant Microbe Interact. 2011;24:1407–18.View ArticlePubMedGoogle Scholar
- Segers GC, Zhang X, Deng F, Sun Q, Nuss DL. Evidence that RNA silencing functions as an antiviral defense mechanism in fungi. Proc Natl Acad Sci U S A. 2007;104:12902–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Hansen FT, Gardiner DM, Lysoe E, Fuertes PR, Tudzynski B, Wiemann P, et al. An update to polyketide synthase and non-ribosomal synthetase genes and nomenclature in Fusarium. Fungal Genet Biol. 2015;75:20–9.View ArticlePubMedGoogle Scholar
- Vujosevic M, Blagojevic J. B chromosomes in populations of mammals. Cytogenet Genome Res. 2004;106:247–56.View ArticlePubMedGoogle Scholar
- Jones N, Houben A. B chromosomes in plants: escapees from the A chromosome genome? Trends Plant Sci. 2003;8:417–23.View ArticlePubMedGoogle Scholar
- Jones RN. Tansley review no 85 - B chromosomes in plants. New Phytol. 1995;131:411–34.View ArticleGoogle Scholar
- Schotanus K, Soyer JL, Connolly LR, Grandaubert J, Happel P, Smith KM, et al. Histone modifications rather than the novel regional centromeres of Zymoseptoria tritici distinguish core and accessory chromosomes. Epigenetics Chromatin. 2015;8:1–18.View ArticleGoogle Scholar
- Nur U, Werren JH, Eickbush DG, Burke WD, Eickbush TH. A selfish B-chromosome that enhances its transmission by eliminating the paternal genome. Science. 1988;240:512–4.View ArticlePubMedGoogle Scholar
- Godfrey SAC, Lovell HC, Mansfield JW, Corry DS, Jackson RW, Arnold DL. The stealth episome: suppression of gene expression on the excised genomic Island PPHGI-1 from Pseudomonas syringae pv. phaseolicola. PLoS Pathog. 2011;7(3):e1002010.View ArticlePubMedPubMed CentralGoogle Scholar
- Hazkani-Covo E, Zeller RM, Martin W. Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes. PLoS Genet. 2010;6(2):e1000834.View ArticlePubMedPubMed CentralGoogle Scholar
- Faino L, Seidl MF, Datema E, van den Berg GC, Janssen A, Wittenberg AH, et al. Single-molecule real-time sequencing combined with optical mapping yields completely finished fungal genome. MBio. 2015;6:e00936–00915.View ArticlePubMedPubMed CentralGoogle Scholar
- Connolly LR, Smith KM, Freitag M. The Fusarium graminearum Histone H3 K27 methyltransferase KMT6 regulates development and expression of secondary metabolite gene clusters. PLoS Genet. 2013;9(10):e1003916.View ArticlePubMedPubMed CentralGoogle Scholar
- Hu X, Xiao G, Zheng P, Shang Y, Su Y, Zhang X, et al. Trajectory and genomic determinants of fungal-pathogen speciation and host adaptation. Proc Natl Acad Sci U S A. 2014;111:16796–801.View ArticlePubMedPubMed CentralGoogle Scholar
- Hu J, Chen C, Peever T, Dang H, Lawrence C, Mitchell T. Genomic characterization of the conditionally dispensable chromosome in Alternaria arborescens provides evidence for horizontal gene transfer. BMC Genomics. 2012;13:171.View ArticlePubMedPubMed CentralGoogle Scholar
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.View ArticlePubMedPubMed CentralGoogle Scholar
- Lukashin AV, Borodovsky M. GeneMark.hmm: new solutions for gene finding. Nucleic Acids Res. 1998;26:1107–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Stanke M, Schoffmann O, Morgenstern B, Waack S. Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources. BMC Bioinformatics. 2006;7:62.View ArticlePubMedPubMed CentralGoogle Scholar
- RepeatModeler Open-1.0 [http://www.repeatmasker.org].
- Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539.View ArticlePubMedPubMed CentralGoogle Scholar
- Hane JK, Oliver RP. RIPCAL: a tool for alignment-based analysis of repeat-induced point mutations in fungal genomic sequences. BMC Bioinformatics. 2008;9:478.View ArticlePubMedPubMed CentralGoogle Scholar
- Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.View ArticlePubMedGoogle Scholar
- Junier T, Zdobnov EM. The Newick utilities: high-throughput phylogenetic tree processing in the Unix shell. Bioinformatics. 2010;26:1669–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Kasuga T, White TJ, Taylor JW. Estimation of nucleotide substitution rates in eurotiomycete fungi. Mol Biol Evol. 2002;19:2318–24.View ArticlePubMedGoogle Scholar
- Camiolo S, Porceddu A. gff2sequence, a new user friendly tool for the generation of genomic sequences. BioData Min. 2013;6(1):15.View ArticlePubMedPubMed CentralGoogle Scholar
- Warris S, Yalcin F, Jackson KJL, Nap JP. Flexible, fast and accurate sequence alignment profiling on GPGPU with PaSWAS. PLoS One. 2015;10(4):e0122524.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, et al. Circos: An information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Audenaert K, Van Broeck R, Bekaert B, De Witte F, Heremans B, Messens K, et al. Fusarium head blight (FHB) in Flanders: population diversity, inter-species associations and DON contamination in commercial winter wheat varieties. Eur J Plant Pathol. 2009;125:445–58.View ArticleGoogle Scholar
- Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, et al. Versatile and open software for comparing large genomes. Genome Biol. 2004;5(2):R12.View ArticlePubMedPubMed CentralGoogle Scholar
- Wicker T, Sabot F, Hua-Van A, Bennetzen JL, Capy P, Chalhoub B, et al. A unified classification system for eukaryotic transposable elements. Nat Rev Genet. 2007;8:973–82.View ArticlePubMedGoogle Scholar
- Yang L, van der Lee T, Yang X, Yu D, Waalwijk C. Fusarium populations on chinese barley show a dramatic gradient in mycotoxin profiles. Phytopathology. 2008;98:719–27.View ArticlePubMedGoogle Scholar