BACKGROUND
genome plasticity is one of the most important drivers of evolutions in eukaryotes.
this plasticity includes large scale genome duplications, rearrangements, deletions and compartmentalization .
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  <cit> .
indeed, the division of the genome into a core and accessory part evolving at different speeds, has been described in multiple pathogens, particularly fungi .
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  <cit> .
recent reviews have noted on the exceptional genome plasticity and particularly the widespread presence of supernumerary chromosomes as a hallmark of pathogenic fungi  <cit> .
for this reason, they represent excellent model organisms to investigate eukaryotic genome evolution  <cit> .
in some cases, the supernumerary chromosomes contain genes involved in pathogenicity and/or efficient host colonization like in alternaria alternata  <cit> , fusarium solani   <cit> , and f. oxysporum f. sp.
lycopersici  <cit> .
in other species, the role of supernumerary chromosomes is less clear as they do not show obvious pathogenicity related functions, such as zymoseptoria tritici   <cit> .
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  <cit> .
different evolutionary pressure on supernumerary chromosomes has also been proposed as an explanation for the detected differences  <cit> .
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  <cit> .
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  <cit> .
transposable elements  play an important role in fungal genome diversity and the evolutionary success of some pathogens  <cit> .
examples are the vast differences in genome sizes of fusarium and phytophthora species  <cit> , and the shaping of pathogenicity in l. maculans and pyrenophora tritici-repentis  <cit> .
the link between te-mediated genome expansion and the evolution of virulence factors has been reviewed extensively  <cit> .
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 , that efficiently inactivates tes by introducing cytosine to thymidine mutations  <cit> .
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  <cit>  and neurospora crassa  <cit> .
this finding has been termed the evolutionary cost of genome defense  <cit> .
rip functions on repetitive sequences with greater than 80 % identity and exceeding +/- 800 bp in length  <cit> .
the rip process occurs only during the di-karyotic pre-meiotic phase and is therefore intricately associated with sexual cycle  <cit> .
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- <dig> or the mat1- <dig> 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  <cit> .
the presumption of active meiosis becomes substantiated when markers for a sexual lifestyle are considered, including recombination  <cit> , rip  <cit> , distributions of the mating idiomorphs in the population  <cit> , and functional constraint on the genes implicated for meiosis  <cit> .
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  <cit> .
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  <cit> .
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  <cit> .
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  <cit> .
individuals of this species were shown to contain a highly variable set of supernumerary chromosomes  <cit> .
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  <cit> .
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.
RESULTS
the genome is composed of a core and supernumerary part
the genome of f. poae isolate  <dig> was assembled from the single-molecule real-time  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  were queried against two smrt assemblies based on different parameters .
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  <dig> chromosomes was obtained, originating from  <dig> contigs , accounting for a total  <dig>  mb of sequence.
the merged assembly of  <dig> chromosomes was supplemented with the remaining contigs and the degenerate unitigs from assembly a, and only the  <dig> non-redundant contigs were kept .
the general statistics of this final assembly can be found in table  <dig>  and are compared with a de novo assembly of the hiseq reads for this isolate.
the smrt assembly is larger by  <dig>  mb, significantly reduces the total number of contigs and has a much larger representation of bases in large contigs .
the  <dig> merged chromosomes contain one, two, two and one telomeres, respectively.
the long arm of chromosome  <dig> misses the telomere in this assembly, while the short arm of chromosome  <dig> ends in the ribosomal dna tandem repeat.
the long arm of chromosome  <dig> contains a 5000 n placeholder at 150 kb from the telomere.
at this junction a 150 kb contig was joined to the rest of the assembly on the basis of its collinearity with f. graminearum and other f. poae isolates.table  <dig> comparison of the smrt and hiseq assemblies of isolate  <dig>  the statistics for the smrt assembly were extracted from the final version of the assembly: four core chromosomes and  <dig> supernumerary contigs
the base quality of the assembly of the four core chromosomes was checked by mapping the hiseq reads of isolate  <dig> to the reference assembly.
over these  <dig>  mb, only one single nucleotide polymorphism  was detected between the hiseq and smrt reads.
two hundred and twenty-two variants were detected in homopolymeric stretches of nucleotides and low complexity  regions .
for these variants, read mapping was inconclusive for both hiseq and smrt reads, with both read batches giving support for different nucleotide calls.
figure  <dig> shows the result of a whole genome alignment of the four core chromosomes of f. poae isolate  <dig> and f. graminearum reference isolate ph- <dig>  the latter was recently assembled to the chromosomal level  <cit> .
aside from two major chromosomal inversions  and several smaller ones, the core chromosomes show extensive macrosynteny between the two species.
moreover, the entire f. graminearum sequence complement is present in our assembly, with the exception of  <dig> copies of the rdna repeat at the end of chromosome  <dig>  two blocks of 204 kb and 464 kb on chromosome  <dig> of f. poae  <dig> do not show synteny with f. graminearum  and are described in detail further below.fig.
<dig> chromosome alignments between f. graminearum  and f. poae .
the best 1: <dig> alignment is shown between the four chromosomes of f. graminearum and the four core chromosomes of f. poae.
red indicates best hits in the same orientation while blue indicates inversions.
the short arm of f. graminearum chromosome four ends in ~ <dig>  mb of rdna repeats that are not assembled in f. poae.
all f. graminearum telomeres except the telomere of the short arm of chromosome  <dig> are assembled.
for f. poae, the same telomere is lacking as well as the one on the long arm of chromosome  <dig>  telomeres that are assembled are shown with green bars on the arms of the chromosomes.
two insertions into f. poae chromosome  <dig> are denoted with black arrows.
approximate locations of the centromeres are shown with black circles
the remaining  <dig> contigs  do not show any synteny with sister species f. graminearum.
these  <dig> contigs contain eight copies of the zit <dig> transposable element  described earlier as a specific marker for supernumerary chromosomes in f. poae  <cit> ; this element was not found on the four core chromosomes.
not only zit <dig>  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  that is the most abundant repetitive element in the genome.
all  <dig> copies  were found dispersed over the core chromosomes.
in sharp contrast no mite was found on the  <dig> 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  was identical to at least one other isolate.
taken together, the  <dig> 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  <dig>  mb that the  <dig> 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  double assembly .
besides the four core chromosomes and  <dig> 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
isolate  <dig> was grown in six diverse conditions, to stimulate transcription of as many genes as possible.
rna was extracted and sequenced, and  <dig>  <dig>  <dig> rnaseq sequence reads were obtained.
these were quality trimmed and the resulting  <dig>  <dig>  <dig> reads were used in the braker <dig> pipeline  <cit> .
this is a novel annotation method that uses rnaseq reads as extrinsic evidence, to annotate the genome in a rapid and automated way without any manual curation steps.
in total  <dig> genes were predicted for isolate  <dig>  table  <dig> lists some core features of the machine annotation of f. poae  <dig> compared to the annotation of f. graminearum ph- <dig>  <cit> .table  <dig> general features of the machine annotation of f. poae  <dig>  these were compared to the features of the published annotation of f. graminearum ph- <dig>  <cit>
f. poae
f. graminearum
the busco  data set for fungi was used to assess whether the annotation can be considered complete  <cit> .
this set comprises proteins that are very likely to be present in a queried genome, based on an analysis of other genomes within a particular kingdom.
the predicted proteins from f. poae  <dig> as well as the proteins from the most recent annotation of f. graminearum ph- <dig>  <cit>  were analyzed by comparing them to the busco data set.
table  <dig> shows the output for both species.
the f. poae protein set is assessed at equally high quality as the f. graminearum set, indicating that the f. poae genome annotation is as accurate and complete as the f. graminearum annotation.table  <dig> busco analyses of f. poae and f. graminearum
f. poae
f. graminearum
fragmented buscos are proteins that are only partially recovered  <cit> .
these were analyzed manually .
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  <dig>  for f. graminearum an rnaseq data set described before was used  <cit> .
the busco analysis suggests that the annotation of f. poae  <dig> did not miss any conserved genes, and within the conserved genes, < <dig>  % 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  <dig>  the mat1- <dig> locus was extracted from the assembly, and its architecture is presented in additional file  <dig>  as all four isolates in this study have the mat1- <dig> mating type, the architecture of the mat1- <dig> locus could not be investigated.
the number, order and direction of the genes occupying the mat1- <dig> locus is identical to that in other fusarium species  <cit> .
the mat1-1- <dig>  mat1-1- <dig> and mat1-1- <dig> genes have a predicted orf with high similarity to those found in related species .
the mat1-1- <dig> gene was previously identified for f. poae  <cit>  and has 99 % protein similarity with the gene model in this study.
transcription of mat1- <dig> idiomorph was noted in the rnaseq data, and the predicted splice forms lead to functional proteins .
a collection of  <dig> isolates was screened for the presence of mat1- <dig> and mat1- <dig> and both idiomorphs were detected, albeit in heavily skewed distribution .
the kegg pathway for meiosis in f. graminearum  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 fgr <dig> give best reciprocal protein hits with f. poae at expect values below 10- <dig>  indicating that all ingredients of the meiotic toolbox are present in f. poae .
only a limited number of genes have been identified that belong to the machinery for rip.
a homolog of the rid  gene, shown to be vital for rip, of f. graminearum is present, intact and transcribed in f. poae .
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
table  <dig> lists the distribution of tes throughout the genome.
figure  <dig> visualizes the chromosomal distribution of tes.
many te families only contain copies on the supernumerary part of the genome.
repeatmasker analysis with the identified tes classified  <dig>  % of the core genome and  <dig>  % of the supernumerary genome as tes.table  <dig> classification and key characteristics of te families in the genome of f. poae  <dig>  elements below the length threshold for rip are not included .
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  <cit> .
r retrotransposon, d dna transposon, l long terminal repeat , t terminal inverted repeat , 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
rlg_maggy
rlg_skippy
rlc_ghost
rxx_marsu
dtf_fot4
dtf_fot8
dtf_fot2
dtf_fot3-a
dtf_fot3-b
dtf_fot3-c
dtf_fot5-a
dtf_fot5-c
dtf_esp4-a
dtf_esp4-b
dtf_drogon
dtf_viserion
dtf_rhaegal
dtf_balerion
dta_rlt1
dta_rlt2
dta_rlt3
dta_hornet1
dta_hornet2
dta_hornet3
dta_tfo1
dta_tfo2
dta_drifter
dta_obara
dta_nymeria
dta_obella
dta_sarella
dtm_hop7
dtm_hop4
fig.
<dig> circos plot showing differences between the core and supernumerary parts of the genome.
outer circle: blue lines denote the distribution of a mite, red triangles denote zit <dig> copies.
second circle: core chromosomes and supernumerary contigs are colored, blue blocks on the chromosomes indicate the centromeres, black blocks show the two insertions of supernumerary sequence into the core chromosomes.
third circle: black lines represent intact  copies of tes.
fourth circle: red lines represent ripped copies of tes.
at the center of the plot, black lines connect gene duplications between the core genome and the supernumerary genome.
only protein hits larger than  <dig> amino acids are shown as their corresponding genes are supposed to be above the length threshold for rip.
duplications within the supernumerary genome are not mapped
transcription and splicing of the predicted introns were noted for many intact elements.
single nucleotide polymorphisms  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  <dig> and annotated in additional file  <dig>  the functional and structural features of every te were used for the classification into superfamilies .
additional file  <dig> 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_hop <dig>  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  <dig> transposable element families, separately for the core and supernumerary genome.
the results are presented in additional file  <dig>  on the core genome, there are  <dig> families that show ripcal patterns that are typical for rip .
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 .
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 .
transposable elements copy number is dynamic between isolates of the same species
figure  <dig> shows the te copy number variation between the isolates in this study.
based on genome coverage the most abundant element in isolate  <dig>  occurs only once in isolate bfb <dig>  a strain originating from china.
as its orf and tirs are intact, it remains unknown why this element has not proliferated in isolate bfb <dig>  isolates  <dig> and  <dig> were isolated from the same belgian field at the same time, but show sharp differences in te copy number.
two families that contain multiple copies in the belgian isolates, are not present as intact copies in the chinese isolate .
however, ripped copies present in the genome of bfb <dig> indicate that during the evolution of the lineage that isolate bfb <dig> belongs to, intact copies of these families have been present but were effectively eradicated from the genome.fig.
<dig> estimation of te numbers in the different f. poae isolates used in the study, as determined by a coverage-based method.
repeat families are classified in decreasing order of incidence in the genome of f. poae 2516; only class i and ii transposable elements that are intact in f. poae  <dig> are included, therefore elements such as the rdna tandem and two families of telomere linked recq helicases are not in the table.
x denotes families for which rip was detected.
it should be noted that average read coverage does not account for possible truncations and therefore the numbers in this table should be considered an estimate
interestingly, rip of certain elements seems isolate-specific, such as dta_rlt <dig> in isolate bfb <dig> and rlc_ghost in isolate  <dig>  a process similar to the loss of dtm_hop <dig> and rlg_maggy in isolate bfb <dig> may have occurred species-wide, as ripped elements in isolate  <dig> were detected of up to  <dig> 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  <dig> 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  <dig>  figure  <dig> 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  <dig> elements on the core chromosomes of isolate  <dig> are present in the same location in isolates  <dig>   <dig> and bfb <dig> .
the integration of these elements therefore seems to have happened recently.
in contrast, on the supernumerary genome of isolate  <dig> elements can be found that show identical integration in isolate  <dig> and one or more of the other isolates.
figure  <dig> illustrates this for supernumerary contig  <dig>  the four tracks show the te presence  and genome coverage  for every isolate.
several elements have identical flanks in all isolates, indicating that they are ancestral integrations .
additional file  <dig> shows the profile for supernumerary contigs  <dig>   <dig> and  <dig>  together these four contigs are the largest supernumerary contigs, totaling  <dig>  mb.
whole blocks of sequence are absent from some isolates: most of contig  <dig> in isolate bfb <dig>  parts of contig  <dig> in isolates  <dig> and bfb <dig>  this absence/presence of sequence on the supernumerary genome is not cumulative for any one isolate or contig investigated.
moreover, the integration of tes on the supernumerary genome is also not concordant with vertical inheritance.
this is illustrated in additional file  <dig>  the recombination-like picture of sequence absence/presence and te integration, may reflect the dynamics the supernumerary chromosomes undergo during crossing.fig.
<dig> integration of intact tes on supernumerary contig  <dig>  the graphs shows in a sliding 1 kb window the fraction of bases from the reference contig that is covered by hiseq reads of every isolate .
the upper track shows all tes on contig  <dig> of isolate  <dig> that are >1 kb and >90 % identity to the element prototype  with yellow dots.
this te landscape was used for comparison with isolates  <dig>   <dig> and bfb <dig>  dots for these three isolates indicate elements for which there is read mapping that an element has integrated in the exact same location as the element in isolate  <dig> .
dots that align vertically are conserved in multiple isolates
a comparison was made between the estimated divergence time of tes in the core genome and those in the supernumerary genome of isolate  <dig>  this is based on the principle that when tes are present at a certain location for a longer period of time, they gradually accumulate more snps, which can be used to calculate the time elapsed since their insertion.
as can be seen in fig.
<dig>  tes in the supernumerary genome are more divergent and are therefore presumed to result from more ancient transposition events.fig.
<dig> divergence estimation of intact  te copies on the core  and supernumerary  genomes.
copies were aligned and branch lengths extracted from a maximum-likelihood phylogenetic tree.
branch lengths were used to calculate divergence times with a fixed substitution rate .
y axis scale was cut off at  <dig> mya, but for the supernumerary genome many outliers are above this value.
additional file  <dig> shows the boxplot with outliers for the supernumerary genome.
the boxes for every te show the lower and upper quartile of the divergence estimates and the median .
the whiskers represent the minimum and maximum values.
circles and asterisks are outliers and extreme values which fall respectively outside of one-and-a-half additional box lengths and three additional box lengths counted from the upper quartile limit
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  <dig>  this te has recently integrated into the genome of isolate  <dig> and contains  <dig> snps compared to the other copies on the core chromosomes of this isolate.
these  <dig> 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  <dig> has sometimes occurred within the coding region of genes.
the environment of all  <dig> recent integrations in the core genome was investigated .
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  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  <dig>  two sequence blocks do not show any synteny with f. graminearum, with coordinates 115073- <dig>  and 883738- <dig> .
analysis of the flanking sequences of these two regions shows that they are continuous in isolates  <dig>   <dig> and bfb <dig>  therefore, these regions represent translocations of supernumerary sequence to the core genome of isolate  <dig>  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  <dig> 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  <dig> .
isolate  <dig> is one of these three isolates.
isolates  <dig>   <dig> and bfb <dig> 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  <cit> .
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  <dig> 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  and represent genuine gene duplications that have not been inactivated by rip.
a total of  <dig> genes on the core chromosomes have one or more duplicates in the supernumerary genome, totaling  <dig> hits on the latter.
figure  <dig> visualizes these gene duplications as lines connecting both the gene on one of the four core chromosomes and its paralog on the supernumerary genome.
additional file  <dig> lists the functional annotation of these duplicated genes.
notable instances include the transcription factor ebr <dig>  <cit> , key component of the rna silencing pathway dicer <dig>  <cit>  and the secondary metabolite backbone gene pks <dig>  <cit> .
in the case of ebr <dig>  up to  <dig> 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  <dig> duplications were identified as likely double assemblies.
the remaining  <dig> instances are likely genuine gene duplications .
additional file  <dig> 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  <dig> shows that some duplications are unique for isolate  <dig>  while others are not cumulative over the different isolates.
indeed, isolates  <dig>   <dig> and bfb <dig> each share unique duplications with isolate  <dig> that are absent in the other isolates.
discussion
chromosomes that vary in both size and number, and have an uneven distribution among individuals of the same species, have been described in animals  <cit> , plants  <cit>  and fungi  <cit> .
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  <cit> , have different gene density and function  <cit>  and/or gc-content  <cit> , are epigenetically very dissimilar  <cit> , can be transmitted by horizontal transfer  <cit>  and are unstable in meiosis  <cit> .
this sharp contrast between the two sets of chromosomes can be explained by different evolutionary pressure, different origins, or a combination of both  <cit> .
the functions of the supernumerary chromosomes are unclear.
speculations range from selfish dna fragments without benefit to the host  <cit>  to components that are critical for pathogenicity and survival  <cit> .
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  <cit> , and mitochondrial dna insertions into the nuclear genome of many eukaryotes  <cit> .
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  <cit> .
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  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  <dig> of f. oxysporum f. sp.
lycopersici has all the characteristics of supernumerary sequence  <cit> .
it has been shown that this region is highly syntenic with one of its supernumerary chromosomes  <cit> .
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  <cit> .
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 ebr <dig> is present as a single copy gene in f. graminearum  <cit> .
in f. oxysporum f. sp.
lycopersici it has undergone gene duplication on the supernumerary chromosomes  <cit> , 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  <cit> , and the advantages of a temporary rip relaxation or arrest may include accelerated evolution and divergence between related species  <cit> .
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  <cit> .
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  <cit> .
specific exclusion of duplicated genes from rip has been observed before in f. solani , a species distantly related to f. poae  <cit> .
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  <cit> , 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.
the regional variability of rip extends beyond supernumerary chromosomes however.
nucleolus organizer regions  contain the rdna tandem repeats in many fungal species, and within this region they are protected from rip  <cit> .
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  <cit> .
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 g <dig> or near the replication fork during the s phase  <cit> .
it functions only in the nucleus, or nuclei, that contain dna duplications and does so on a single dna strand.
mis-pairing of duplicated dna has been hypothesized to deliver the substrate for rip  <cit> .
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  <cit>  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  <cit> .
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  <cit> .
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  <cit> , transcription  <cit>  and gene content  <cit> .
CONCLUSIONS
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.

 genome plasticity is one of the most important drivers of evolutions in eukaryotes. this plasticity includes large scale genome duplications, rearrangements, deletions and compartmentalization. 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. 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. 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, 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. 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. different evolutionary pressure on supernumerary chromosomes has also been proposed as an explanation for the detected differences. 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. 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 play an important role in fungal genome diversity and the evolutionary success of some pathogens. examples are the vast differences in genome sizes of fusarium and phytophthora species, and the shaping of pathogenicity in l. maculans and pyrenophora tritici-repentis. 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, 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- or the mat1- 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. 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 was assembled from the single-molecule real-time 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 were queried against two smrt assemblies based on different parameters. 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 chromosomes was obtained, originating from contigs, accounting for a total  mb of sequence. the merged assembly of chromosomes was supplemented with the remaining contigs and the degenerate unitigs from assembly a, and only the non-redundant contigs were kept. the general statistics of this final assembly can be found in table  and are compared with a de novo assembly of the hiseq reads for this isolate. the smrt assembly is larger by  mb, significantly reduces the total number of contigs and has a much larger representation of bases in large contigs. the merged chromosomes contain one, two, two and one telomeres, respectively. the long arm of chromosome misses the telomere in this assembly, while the short arm of chromosome ends in the ribosomal dna tandem repeat. the long arm of chromosome contains a 5000 n placeholder at 150 kb from the telomere. at this junction a 150 kb contig was joined to the rest of the assembly on the basis of its collinearity with f. graminearum and other f. poae isolates.table comparison of the smrt and hiseq assemblies of isolate the statistics for the smrt assembly were extracted from the final version of the assembly: four core chromosomes and supernumerary contigs the base quality of the assembly of the four core chromosomes was checked by mapping the hiseq reads of isolate to the reference assembly. over these  mb, only one single nucleotide polymorphism was detected between the hiseq and smrt reads. two hundred and twenty-two variants were detected in homopolymeric stretches of nucleotides and low complexity regions. for these variants, read mapping was inconclusive for both hiseq and smrt reads, with both read batches giving support for different nucleotide calls. figure  shows the result of a whole genome alignment of the four core chromosomes of f. poae isolate and f. graminearum reference isolate ph- the latter was recently assembled to the chromosomal level. aside from two major chromosomal inversions and several smaller ones, the core chromosomes show extensive macrosynteny between the two species. moreover, the entire f. graminearum sequence complement is present in our assembly, with the exception of copies of the rdna repeat at the end of chromosome two blocks of 204 kb and 464 kb on chromosome of f. poae do not show synteny with f. graminearum and are described in detail further below.fig. chromosome alignments between f. graminearum and f. poae. the best 1: alignment is shown between the four chromosomes of f. graminearum and the four core chromosomes of f. poae. red indicates best hits in the same orientation while blue indicates inversions. the short arm of f. graminearum chromosome four ends in ~  mb of rdna repeats that are not assembled in f. poae. all f. graminearum telomeres except the telomere of the short arm of chromosome are assembled. for f. poae, the same telomere is lacking as well as the one on the long arm of chromosome telomeres that are assembled are shown with green bars on the arms of the chromosomes. two insertions into f. poae chromosome are denoted with black arrows. approximate locations of the centromeres are shown with black circles the remaining contigs do not show any synteny with sister species f. graminearum. these contigs contain eight copies of the zit transposable element described earlier as a specific marker for supernumerary chromosomes in f. poae; this element was not found on the four core chromosomes. not only zit 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 that is the most abundant repetitive element in the genome. all copies were found dispersed over the core chromosomes. in sharp contrast no mite was found on the 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 was identical to at least one other isolate. taken together, the 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  mb that the 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 double assembly. besides the four core chromosomes and 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 isolate was grown in six diverse conditions, to stimulate transcription of as many genes as possible. rna was extracted and sequenced, and rnaseq sequence reads were obtained. these were quality trimmed and the resulting reads were used in the braker pipeline. this is a novel annotation method that uses rnaseq reads as extrinsic evidence, to annotate the genome in a rapid and automated way without any manual curation steps. in total genes were predicted for isolate table  lists some core features of the machine annotation of f. poae compared to the annotation of f. graminearum ph-.table general features of the machine annotation of f. poae these were compared to the features of the published annotation of f. graminearum ph- f. poae f. graminearum the busco data set for fungi was used to assess whether the annotation can be considered complete. this set comprises proteins that are very likely to be present in a queried genome, based on an analysis of other genomes within a particular kingdom. the predicted proteins from f. poae as well as the proteins from the most recent annotation of f. graminearum ph- were analyzed by comparing them to the busco data set. table  shows the output for both species. the f. poae protein set is assessed at equally high quality as the f. graminearum set, indicating that the f. poae genome annotation is as accurate and complete as the f. graminearum annotation.table busco analyses of f. poae and f. graminearum f. poae f. graminearum fragmented buscos are proteins that are only partially recovered. these were analyzed manually. 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 for f. graminearum an rnaseq data set described before was used. the busco analysis suggests that the annotation of f. poae did not miss any conserved genes, and within the conserved genes, <  % 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 the mat1- locus was extracted from the assembly, and its architecture is presented in additional file as all four isolates in this study have the mat1- mating type, the architecture of the mat1- locus could not be investigated. the number, order and direction of the genes occupying the mat1- locus is identical to that in other fusarium species. the mat1-1- mat1-1- and mat1-1- genes have a predicted orf with high similarity to those found in related species. the mat1-1- gene was previously identified for f. poae and has 99 % protein similarity with the gene model in this study. transcription of mat1- idiomorph was noted in the rnaseq data, and the predicted splice forms lead to functional proteins. a collection of isolates was screened for the presence of mat1- and mat1- and both idiomorphs were detected, albeit in heavily skewed distribution. the kegg pathway for meiosis in f. graminearum 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 fgr give best reciprocal protein hits with f. poae at expect values below 10- indicating that all ingredients of the meiotic toolbox are present in f. poae. only a limited number of genes have been identified that belong to the machinery for rip. a homolog of the rid gene, shown to be vital for rip, of f. graminearum is present, intact and transcribed in f. poae. 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 table  lists the distribution of tes throughout the genome. figure  visualizes the chromosomal distribution of tes. many te families only contain copies on the supernumerary part of the genome. repeatmasker analysis with the identified tes classified  % of the core genome and  % of the supernumerary genome as tes.table classification and key characteristics of te families in the genome of f. poae elements below the length threshold for rip are not included. 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, t terminal inverted repeat, 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  rlg_maggy  rlg_skippy  rlc_ghost  rxx_marsu  dtf_fot4  dtf_fot8  dtf_fot2  dtf_fot3-a  dtf_fot3-b  dtf_fot3-c  dtf_fot5-a  dtf_fot5-c  dtf_esp4-a  dtf_esp4-b  dtf_drogon  dtf_viserion  dtf_rhaegal  dtf_balerion  dta_rlt1  dta_rlt2  dta_rlt3  dta_hornet1  dta_hornet2  dta_hornet3  dta_tfo1  dta_tfo2  dta_drifter  dta_obara  dta_nymeria  dta_obella  dta_sarella  dtm_hop7  dtm_hop4 fig. circos plot showing differences between the core and supernumerary parts of the genome. outer circle: blue lines denote the distribution of a mite, red triangles denote zit copies. second circle: core chromosomes and supernumerary contigs are colored, blue blocks on the chromosomes indicate the centromeres, black blocks show the two insertions of supernumerary sequence into the core chromosomes. third circle: black lines represent intact copies of tes. fourth circle: red lines represent ripped copies of tes. at the center of the plot, black lines connect gene duplications between the core genome and the supernumerary genome. only protein hits larger than amino acids are shown as their corresponding genes are supposed to be above the length threshold for rip. duplications within the supernumerary genome are not mapped transcription and splicing of the predicted introns were noted for many intact elements. single nucleotide polymorphisms 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 and annotated in additional file the functional and structural features of every te were used for the classification into superfamilies. additional file 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_hop 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 transposable element families, separately for the core and supernumerary genome. the results are presented in additional file on the core genome, there are families that show ripcal patterns that are typical for rip. 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. 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. transposable elements copy number is dynamic between isolates of the same species figure  shows the te copy number variation between the isolates in this study. based on genome coverage the most abundant element in isolate occurs only once in isolate bfb a strain originating from china. as its orf and tirs are intact, it remains unknown why this element has not proliferated in isolate bfb isolates and were isolated from the same belgian field at the same time, but show sharp differences in te copy number. two families that contain multiple copies in the belgian isolates, are not present as intact copies in the chinese isolate. however, ripped copies present in the genome of bfb indicate that during the evolution of the lineage that isolate bfb belongs to, intact copies of these families have been present but were effectively eradicated from the genome.fig. estimation of te numbers in the different f. poae isolates used in the study, as determined by a coverage-based method. repeat families are classified in decreasing order of incidence in the genome of f. poae 2516; only class i and ii transposable elements that are intact in f. poae are included, therefore elements such as the rdna tandem and two families of telomere linked recq helicases are not in the table. x denotes families for which rip was detected. it should be noted that average read coverage does not account for possible truncations and therefore the numbers in this table should be considered an estimate interestingly, rip of certain elements seems isolate-specific, such as dta_rlt in isolate bfb and rlc_ghost in isolate a process similar to the loss of dtm_hop and rlg_maggy in isolate bfb may have occurred species-wide, as ripped elements in isolate were detected of up to 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 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 figure  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 elements on the core chromosomes of isolate are present in the same location in isolates and bfb. the integration of these elements therefore seems to have happened recently. in contrast, on the supernumerary genome of isolate elements can be found that show identical integration in isolate and one or more of the other isolates. figure  illustrates this for supernumerary contig the four tracks show the te presence and genome coverage for every isolate. several elements have identical flanks in all isolates, indicating that they are ancestral integrations. additional file shows the profile for supernumerary contigs and together these four contigs are the largest supernumerary contigs, totaling  mb. whole blocks of sequence are absent from some isolates: most of contig in isolate bfb parts of contig in isolates and bfb this absence/presence of sequence on the supernumerary genome is not cumulative for any one isolate or contig investigated. moreover, the integration of tes on the supernumerary genome is also not concordant with vertical inheritance. this is illustrated in additional file the recombination-like picture of sequence absence/presence and te integration, may reflect the dynamics the supernumerary chromosomes undergo during crossing.fig. integration of intact tes on supernumerary contig the graphs shows in a sliding 1 kb window the fraction of bases from the reference contig that is covered by hiseq reads of every isolate. the upper track shows all tes on contig of isolate that are >1 kb and >90 % identity to the element prototype with yellow dots. this te landscape was used for comparison with isolates and bfb dots for these three isolates indicate elements for which there is read mapping that an element has integrated in the exact same location as the element in isolate. dots that align vertically are conserved in multiple isolates a comparison was made between the estimated divergence time of tes in the core genome and those in the supernumerary genome of isolate this is based on the principle that when tes are present at a certain location for a longer period of time, they gradually accumulate more snps, which can be used to calculate the time elapsed since their insertion. as can be seen in fig.  tes in the supernumerary genome are more divergent and are therefore presumed to result from more ancient transposition events.fig. divergence estimation of intact te copies on the core and supernumerary genomes. copies were aligned and branch lengths extracted from a maximum-likelihood phylogenetic tree. branch lengths were used to calculate divergence times with a fixed substitution rate. y axis scale was cut off at mya, but for the supernumerary genome many outliers are above this value. additional file shows the boxplot with outliers for the supernumerary genome. the boxes for every te show the lower and upper quartile of the divergence estimates and the median. the whiskers represent the minimum and maximum values. circles and asterisks are outliers and extreme values which fall respectively outside of one-and-a-half additional box lengths and three additional box lengths counted from the upper quartile limit 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 this te has recently integrated into the genome of isolate and contains snps compared to the other copies on the core chromosomes of this isolate. these 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 has sometimes occurred within the coding region of genes. the environment of all recent integrations in the core genome was investigated. 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 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 two sequence blocks do not show any synteny with f. graminearum, with coordinates 115073- and 883738-. analysis of the flanking sequences of these two regions shows that they are continuous in isolates and bfb therefore, these regions represent translocations of supernumerary sequence to the core genome of isolate 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 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. isolate is one of these three isolates. isolates and bfb 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 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 and represent genuine gene duplications that have not been inactivated by rip. a total of genes on the core chromosomes have one or more duplicates in the supernumerary genome, totaling hits on the latter. figure  visualizes these gene duplications as lines connecting both the gene on one of the four core chromosomes and its paralog on the supernumerary genome. additional file lists the functional annotation of these duplicated genes. notable instances include the transcription factor ebr, key component of the rna silencing pathway dicer and the secondary metabolite backbone gene pks. in the case of ebr up to 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 duplications were identified as likely double assemblies. the remaining instances are likely genuine gene duplications. additional file 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 shows that some duplications are unique for isolate while others are not cumulative over the different isolates. indeed, isolates and bfb each share unique duplications with isolate that are absent in the other isolates. discussion 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. 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 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 of f. oxysporum f. sp. lycopersici has all the characteristics of supernumerary sequence. 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. 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 ebr 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. specific exclusion of duplicated genes from rip has been observed before in f. solani, 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. the regional variability of rip extends beyond supernumerary chromosomes however. nucleolus organizer regions 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 g or near the replication fork during the s phase. it functions only in the nucleus, or nuclei, that contain 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. 
the supernumerary genome is a reservoir for tes that migrate to the core genome, and even large blocks of supernumerary sequence have recently translocated to the core. a “living apart together” crosstalk exists between the core and supernumerary genome, accelerating chromosomal and organismal evolution. vice versa, the supernumerary genome acts as a refuge for genes that are duplicated from the core genome. the core genome contains limited transposable elements and no gene duplications, while the supernumerary genome holds up to 25 % tes and multiple gene duplications. single-molecule real-time sequencingsupernumerary chromosomesfusariumrepeat-induced point mutationtranslocationgene duplicationstransposable elementsghent university special research fund01di2314vanheule adriaan fonds wetenschappelijk onderzoek k147115nvanheule adriaan dutch ministry of economic affairs, agriculture and innovationdutch ministry of economic affairs, agriculture and innovationissue-copyright-statement© the author 2016 the absence of rip on the supernumerary genome accounts for the differences between the two genomes, and results in a functional crosstalk between them. here we use single-molecule real-time sequencing to assemble the genome of a prominent fungal wheat pathogen, fusarium poae, including at least one supernumerary chromosome. for the first time, a mechanism was determined that explains the differences that exist between the core and supernumerary genome in fungi. 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. the core genome shows all hallmarks of repeat-induced point mutation, a defense mechanism against tes, specific for fungi.
 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 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 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, a defense mechanism against tes, specific for fungi. the absence of rip on the supernumerary genome accounts for the differences between the two 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 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. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords single-molecule real-time sequencingsupernumerary chromosomesfusariumrepeat-induced point mutationtranslocationgene duplicationstransposable elementsghent university special research fund01di2314vanheule adriaan fonds wetenschappelijk onderzoek k147115nvanheule adriaan dutch ministry of economic affairs, agriculture and innovationdutch ministry of economic affairs, agriculture and innovationissue-copyright-statement© the author 2016 