Genome-wide analysis of the Fusarium oxysporum mimp family of MITEs and mobilization of both native and de novo created mimps

Fusarium oxysporum effector clustering version 2: An updated pipeline to infer host range

The fungus Fusarium oxysporum is infamous for its devastating effects on economically important crops worldwide. F. oxysporum isolates are grouped into formae speciales based on their ability to cause disease on different hosts. Assigning F. oxysporum strains to formae speciales using non-experimental procedures has proven to be challenging due to their genetic heterogeneity and polyphyletic nature. However, genetically diverse isolates of the same forma specialis encode similar repertoires of effectors, proteins that are secreted by the fungus and contribute to the establishment of compatibility with the host. Based on this observation, we previously designed the F. oxysporum Effector Clustering (FoEC) pipeline which is able to classify F. oxysporum strains by forma specialis based on hierarchical clustering of the presence of predicted putative effector sequences, solely using genome assemblies as input. Here we present the updated FoEC2 pipeline which is more user friendly, customizable and, due to multithreading, has improved scalability. It is designed as a Snakemake pipeline and incorporates a new interactive visualization app. We showcase FoEC2 by clustering 537 publicly available F. oxysporum genomes and further analysis of putative effector families as multiple sequence alignments. We confirm classification of isolates into formae speciales and are able to further identify their subtypes. The pipeline is available on github: https://github.com/pvdam3/FoEC2.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Basal Rot of Narcissus: Understanding Pathogenicity in Fusarium oxysporum f. sp. narcissi

Fusarium oxysporum is a globally distributed soilborne fungal pathogen causing root rots, bulb rots, crown rots and vascular wilts on a range of horticultural plants. Pathogenic F. oxysporum isolates are highly host specific and are classified as formae speciales. Narcissus is an important ornamental crop and both the quality and yield of flowers and bulbs can be severely affected by a basal rot caused by F. oxysporum f. sp. narcissi (FON); 154 Fusarium isolates were obtained from different locations and Narcissus cultivars in the United Kingdom, representing a valuable resource. A subset of 30 F. oxysporum isolates were all found to be pathogenic and were therefore identified as FON. Molecular characterisation of isolates through sequencing of three housekeeping genes, suggested a monophyletic origin with little divergence. PCR detection of 14 Secreted in Xylem (SIX) genes, previously shown to be associated with pathogenicity in other F. oxysporum f. spp., revealed different complements of SIX7, SIX9, SIX10, SIX12 and SIX13 within FON isolates which may suggest a race structure. SIX gene sequences were unique to FON and SIX10 was present in all isolates, allowing for molecular identification of FON for the first time. The genome of a highly pathogenic isolate was sequenced and lineage specific (LS) regions identified which harboured putative effectors including the SIX genes. Real-time RT-PCR, showed that SIX genes and selected putative effectors were expressed in planta with many significantly upregulated during infection. This is the first study to characterise molecular variation in FON and provide an analysis of the FON genome. Identification of expressed genes potentially associated with virulence provides the basis for future functional studies and new targets for molecular diagnostics.

The Distribution of Miniature Impala Elements and SIX Genes in the Fusarium Genus is Suggestive of Horizontal Gene Transfer

The mimp family of miniature inverted-repeat transposable elements was previously found only in genomes of Fusarium oxysporum and is contextually associated with virulence genes in this species. Through extensive comparative analysis of 83 F. oxysporum and 52 other Fusarium genomes, we uncovered the distribution of different mimp families throughout the genus. We show that (i) mimps are not exclusive to F. oxysporum; (ii) pathogenic isolates generally possess more mimps than non-pathogenic strains and (iii) two isolates of F. hostae and one F. proliferatum isolate display evidence for horizontal transfer of genetic material to or from F. oxysporum. Multiple instances of mimp elements identical to F. oxysporum mimps were encountered in the genomes of these isolates. Moreover, homologs of effector genes (SIX1, 2, 6, 7, 11 and FomAVR2) were discovered here, several with very high (97-100%) pairwise nucleotide sequence identity scores. These three strains were isolated from infected flower bulbs (Hyacinthus and Lilium spp.). Their ancestors may thus have lived in close proximity to pathogenic strains of F. oxysporum f. sp. hyacinthi and f. sp. lilii. The Fo f. sp. lycopersici SIX2 effector gene was found to be widely distributed (15/18 isolates) throughout the F. fujikuroi species complex, exhibiting a predominantly vertical inheritance pattern. These findings shed light on the potential evolutionary mechanism underlying plant-pathogenicity in Fusarium and show that interspecies horizontal gene transfer may have occurred.

The characteristics and functions of a miniature inverted-repeat transposable element TaMITE81 in the 5' UTR of TaCHS7BL from Triticum aestivum

Miniature inverted-repeat transposable elements (MITEs) are truncated derivatives of autonomous DNA transposons, and are dispersed abundantly in eukaryotic and prokaryotic genomes. In this article, a MITE, TaMITE81, was isolated from the 5' untranslated region (UTR) of TaCHS7BL, chalcone synthase (CHS) catalyzing the first committed step of anthocyanin biosynthesis, in the wheat cultivar 'Opata' with white grain. TaMITE81 was only 81 nucleotides, including a terminal inverted repeat with 39 nucleotides and was flanked by two nucleotides, "TA", target site duplications that were typical features of stowaway-like MITEs. Compared with the wheat cultivar 'Gy115' with purple grain, which is without the insertion, the expression of TaCHS7BL was lower in several organs of 'Opata'. The insertion of TaMITE81 into the 5' UTR of the GUS gene also reduced the transient expression of GUS on the coleoptiles of 'Opata', which means the insertion of TaMITE81 was the reason for the low expression of TaCHS7BL in 'Opata'. But the genotype of TaCHS7BL was not linked to phenotype of grain color in the RILs derived from a cross 'Gy115' and 'Opata'. The TaMITE81 density of the hexaploid variety of T. aestivum was more than 10 times that of diploid relatives, which implies that polyploidization caused the amplification of TaMITE81 homologous sequences. Further research should be conducted on decoding the relationship between TaCHS7BL and other traits relative to anthocyanin biosynthesis in wheat, and discovering the mechanism of TaMITE81 transposon action.

Effector profiles distinguish formae speciales of Fusarium oxysporum

Formae speciales (ff.spp.) of the fungus Fusarium oxysporum are often polyphyletic within the species complex, making it impossible to identify them on the basis of conserved genes. However, sequences that determine host-specific pathogenicity may be expected to be similar between strains within the same forma specialis. Whole genome sequencing was performed on strains from five different ff.spp. (cucumerinum, niveum, melonis, radicis-cucumerinum and lycopersici). In each genome, genes for putative effectors were identified based on small size, secretion signal, and vicinity to a "miniature impala" transposable element. The candidate effector genes of all genomes were collected and the presence/absence patterns in each individual genome were clustered. Members of the same forma specialis turned out to group together, with cucurbit-infecting strains forming a supercluster separate from other ff.spp. Moreover, strains from different clonal lineages within the same forma specialis harbour identical effector gene sequences, supporting horizontal transfer of genetic material. These data offer new insight into the genetic basis of host specificity in the F. oxysporum species complex and show that (putative) effectors can be used to predict host specificity in F. oxysporum.

Evolutionary dynamics of hAT DNA transposon families in Saccharomycetaceae

Transposable elements (TEs) are widespread in eukaryotes but uncommon in yeasts of the Saccharomycotina subphylum, in terms of both host species and genome fraction. The class II elements are especially scarce, but the hAT element Rover is a noteworthy exception that deserves further investigation. Here, we conducted a genome-wide analysis of hAT elements in 40 ascomycota. A novel family, Roamer, was found in three species, whereas Rover was detected in 15 preduplicated species from Kluyveromyces, Eremothecium, and Lachancea genera, with up to 41 copies per genome. Rover acquisition seems to have occurred by horizontal transfer in a common ancestor of these genera. The detection of remote Rover copies in Naumovozyma dairenensis and in the sole Saccharomyces cerevisiae strain AWRI1631, without synteny, suggests that two additional independent horizontal transfers took place toward these genomes. Such patchy distribution of elements prevents any anticipation of TE presence in incoming sequenced genomes, even closely related ones. The presence of both putative autonomous and defective Rover copies, as well as their diversification into five families, indicate particular dynamics of Rover elements in the Lachancea genus. Especially, we discovered the first miniature inverted-repeat transposable elements (MITEs) to be described in yeasts, together with their parental autonomous copies. Evidence of MITE insertion polymorphism among Lachancea waltii strains suggests their recent activity. Moreover, 40% of Rover copies appeared to be involved in chromosome rearrangements, showing the large structural impact of TEs on yeast genome and opening the door to further investigations to understand their functional and evolutionary consequences.

Reverse transcriptase and intron number evolution

BACKGROUND

Introns are universal in eukaryotic genomes and play important roles in transcriptional regulation, mRNA export to the cytoplasm, nonsense-mediated decay as both a regulatory and a splicing quality control mechanism, R-loop avoidance, alternative splicing, chromatin structure, and evolution by exon-shuffling.

METHODS

Sixteen complete fungal genomes were used 13 of which were sequenced and annotated by JGI. Ustilago maydis, Cryptococcus neoformans, and Coprinus cinereus (also named Coprinopsis cinerea) were from the Broad Institute. Gene models from JGI-annotated genomes were taken from the GeneCatalog track that contained the best representative gene models. Varying fractions of the GeneCatalog were manually curated by external users. For clarity, we used the JGI unique database identifier.

RESULTS

The last common ancestor of eukaryotes (LECA) has an estimated 6.4 coding exons per gene (EPG) and evolved into the diverse eukaryotic life forms, which is recapitulated by the development of a stem cell. We found a parallel between the simulated reverse transcriptase (RT)-mediated intron loss and the comparative analysis of 16 fungal genomes that spanned a wide range of intron density. Although footprints of RT (RTF) were dynamic, relative intron location (RIL) to the 5'-end of mRNA faithfully traced RT-mediated intron loss and revealed 7.7 EPG for LECA. The mode of exon length distribution was conserved in simulated intron loss, which was exemplified by the shared mode of 75 nt between fungal and Chlamydomonas genomes. The dominant ancient exon length was corroborated by the average exon length of the most intron-rich genes in fungal genomes and consistent with ancient protein modules being ~25 aa. Combined with the conservation of a protein length of 400 aa, the earliest ancestor of eukaryotes could have 16 EPG. During earlier evolution, Ascomycota's ancestor had significantly more 3'-biased RT-mediated intron loss that was followed by dramatic RTF loss. There was a down trend of EPG from more conserved to less conserved genes. Moreover, species-specific genes have higher exon-densities, shorter exons, and longer introns when compared to genes conserved at the phylum level. However, intron length in species-specific genes became shorter than that of genes conserved in all species after genomes experiencing drastic intron loss. The estimated EPG from the most frequent exon length is more than double that from the RIL method.

CONCLUSIONS

This implies significant intron loss during the very early period of eukaryotic evolution. De novo gene-birth contributes to shorter exons, longer introns, and higher exon-density in species-specific genes relative to conserved genes.

Insights on genome size evolution from a miniature inverted repeat transposon driving a satellite DNA

The genome size in eukaryotes does not correlate well with the number of genes they contain. We can observe this so-called C-value paradox in amphibian species. By analyzing an amphibian genome we asked how repetitive DNA can impact genome size and architecture. We describe here our discovery of a Tc1/mariner miniature inverted-repeat transposon family present in Xenopus frogs. These transposons named miDNA4 are unique since they contain a satellite DNA motif. We found that miDNA4 measured 331 bp, contained 25 bp long inverted terminal repeat sequences and a sequence motif of 119 bp present as a unique copy or as an array of 2-47 copies. We characterized the structure, dynamics, impact and evolution of the miDNA4 family and its satellite DNA in Xenopus frog genomes. This led us to propose a model for the evolution of these two repeated sequences and how they can synergize to increase genome size.

Miniature inverted-repeat transposable elements: discovery, distribution, and activity

Eukaryotic organisms have dynamic genomes, with transposable elements (TEs) as a major contributing factor. Although the large autonomous TEs can significantly shape genomic structures during evolution, genomes often harbor more miniature nonautonomous TEs that can infest genomic niches where large TEs are rare. In spite of their cut-and-paste transposition mechanisms that do not inherently favor copy number increase, miniature inverted-repeat transposable elements (MITEs) are abundant in eukaryotic genomes and exist in high copy numbers. Based on the large number of MITE families revealed in previous studies, accurate annotation of MITEs, particularly in newly sequenced genomes, will identify more genomes highly rich in these elements. Novel families identified from these analyses, together with the currently known families, will further deepen our understanding of the origins, transposase sources, and dramatic amplification of these elements.

MITEs in the promoters of effector genes allow prediction of novel virulence genes in Fusarium oxysporum

BACKGROUND

The plant-pathogenic fungus Fusarium oxysporum f.sp.lycopersici (Fol) has accessory, lineage-specific (LS) chromosomes that can be transferred horizontally between strains. A single LS chromosome in the Fol4287 reference strain harbors all known Fol effector genes. Transfer of this pathogenicity chromosome confers virulence to a previously non-pathogenic recipient strain. We hypothesize that expression and evolution of effector genes is influenced by their genomic context.

RESULTS

To gain a better understanding of the genomic context of the effector genes, we manually curated the annotated genes on the pathogenicity chromosome and identified and classified transposable elements. Both retro- and DNA transposons are present with no particular overrepresented class. Retrotransposons appear evenly distributed over the chromosome, while DNA transposons tend to concentrate in large chromosomal subregions. In general, genes on the pathogenicity chromosome are dispersed within the repeat landscape. Effector genes are present within subregions enriched for DNA transposons. A miniature Impala (mimp) is always present in their promoters. Although promoter deletion studies of two effector gene loci did not reveal a direct function of the mimp for gene expression, we were able to use proximity to a mimp as a criterion to identify new effector gene candidates. Through xylem sap proteomics we confirmed that several of these candidates encode proteins secreted during plant infection.

CONCLUSIONS

Effector genes in Fol reside in characteristic subregions on a pathogenicity chromosome. Their genomic context allowed us to develop a method for the successful identification of novel effector genes. Since our approach is not based on effector gene similarity, but on unique genomic features, it can easily be extended to identify effector genes in Fo strains with different host specificities.

Genome-wide comparative analysis of pogo-like transposable elements in different Fusarium species

The recent availability of genome sequences of four different Fusarium species offers the opportunity to perform extensive comparative analyses, in particular of repeated sequences. In a recent work, the overall content of such sequences in the genomes of three phylogenetically related Fusarium species, F. graminearum, F. verticillioides, and F. oxysporum f. sp. lycopersici has been estimated. In this study, we present an exhaustive characterization of pogo-like elements, named Fots, in four Fusarium genomes. Overall 10 Fot and two Fot-related miniature inverted-repeat transposable element families were identified, revealing a diversification of multiple lineages of pogo-like elements, some of which accompanied by a gain of introns. This analysis also showed that such elements are present in an unusual high proportion in the genomes of F. oxysporum f. sp. lycopersici and Nectria haematococca (anamorph F. solani f. sp. pisi) in contrast with most other fungal genomes in which retroelements are the most represented. Interestingly, our analysis showed that the most numerous Fot families all contain potentially active or mobilisable copies, thus conferring a mutagenic potential of these transposable elements and consequently a role in strain adaptation and genome evolution. This role is strongly reinforced when examining their genomic distribution which is clearly biased with a high proportion (more than 80%) located on strain- or species-specific regions enriched in genes involved in pathogenicity and/or adaptation. Finally, the different reproductive characteristics of the four Fusarium species allowed us to investigate the impact of the process of repeat-induced point mutations on the expansion and diversification of Fot elements.

Abundant degenerate miniature inverted-repeat transposable elements in genomes of epichloid fungal endophytes of grasses

Miniature inverted-repeat transposable elements (MITEs) are abundant repeat elements in plant and animal genomes; however, there are few analyses of these elements in fungal genomes. Analysis of the draft genome sequence of the fungal endophyte Epichloë festucae revealed 13 MITE families that make up almost 1% of the E. festucae genome, and relics of putative autonomous parent elements were identified for three families. Sequence and DNA hybridization analyses suggest that at least some of the MITEs identified in the study were active early in the evolution of Epichloë but are not found in closely related genera. Analysis of MITE integration sites showed that these elements have a moderate integration site preference for 5' genic regions of the E. festucae genome and are particularly enriched near genes for secondary metabolism. Copies of the EFT-3m/Toru element appear to have mediated recombination events that may have abolished synthesis of two fungal alkaloids in different epichloae. This work provides insight into the potential impact of MITEs on epichloae evolution and provides a foundation for analysis in other fungal genomes.

New virulence groups in Fusarium oxysporum f. sp. phaseoli: the expression of the gene coding for the transcription factor ftf1 correlates with virulence

Fusarium oxysporum f. sp. phaseoli strains isolated from runner bean plants showing Fusarium wilt symptoms were characterized. The analysis of the genetic diversity of these strains and the comparison with strains formerly isolated from diseased common bean plants grown in the same region of Spain indicated a close genetic similarity among them. Pathogenicity assays carried out on runner bean plants showed virulence differences that allowed the classification of these strains into three groups: super virulent, highly virulent, and weakly virulent. However, all the analyzed strains behaved as highly virulent when inoculated on common bean plants, indicating that virulence is specific of the host-pathogen interaction. We also analyzed the number of copies and expression of the gene encoding the transcription factor ftf1, which has been shown to be specific of virulent F. oxysporum strains and highly up-regulated during plant infection. In planta real-time quantitative polymerase chain reaction expression analysis showed that expression of ftf1 was correlated with the degree of virulence. The comparative analysis of the polymorphic copies of ftf1 detected in the strains here characterized and those detected in the genome sequence of F. oxysporum f. sp. lycopersici strain 4287 indicates that some of the copies are likely nonfunctional.

Characterization of a group of MITEs with unusual features from two coral genomes

Background

Miniature inverted-repeat transposable elements (MITEs), which are common in eukaryotic genomes, are small non-coding elements that transpose by utilizing transposases encoded by autonomous transposons. Recent genome-wide analyses and cross-mobilization assays have greatly improved our knowledge on MITE proliferation, however, specific mechanisms for the origin and evolution of MITEs are still unclear.

Principal Findings

A group of coral MITEs called CMITE were identified from two corals, Acropora millepora and Acropora palmata. CMITEs conform to many common characteristics of MITEs, but also present several unusual features. The most unusual feature of CMITEs is conservation of the internal region, which is more conserved between MITE families than the TIRs. The origin of this internal region remains unknown, although we found one CMITE family that seems to be derived from a piggyBac-like transposon in A. millepora. CMITEs can form tandem arrays, suggesting an unconventional way for MITEs to increase copy numbers. We also describe a case in which a novel transposable element was created by a CMITE insertion event.

Conclusions

To our knowledge, this is the first report of identification of MITEs from coral genomes. Proliferation of CMITEs seems to be related to the transposition machinery of piggyBac-like autonomous transposons. The highly conserved internal region of CMITEs suggests a potential role for this region in their successful transposition. However, the origin of these unusual features in CMITEs remains unclear, and thus represents an intriguing topic for future investigations.