The landscape of transposable elements in the finished genome of the fungal wheat pathogen Mycosphaerella graminicola

Defense Pathways of Wheat Plants Inoculated with Zymoseptoria tritici under NaCl Stress Conditions: An Overview

Due to being sessile, plants develop a broad range of defense pathways when they face abiotic or biotic stress factors. Although plants are subjected to more than one type of stress at a time in nature, the combined effects of either multiple stresses of one kind (abiotic or biotic) or more kinds (abiotic and biotic) have now been realized in agricultural lands due to increases in global warming and environmental pollution, along with population increases. Soil-borne pathogens, or pathogens infecting aerial parts, can have devastating effects on plants when combined with other stressors. Obtaining yields or crops from sensitive or moderately resistant plants could be impossible, and it could be very difficult from resistant plants. The mechanisms of combined stress in many plants have previously been studied and elucidated. Recent studies proposed new defense pathways and mechanisms through signaling cascades. In light of these mechanisms, it is now time to develop appropriate strategies for crop protection under multiple stress conditions. This may involve using disease-resistant or stress-tolerant plant varieties, implementing proper irrigation and drainage practices, and improving soil quality. However, generation of both stress-tolerant and disease-resistant crop plants is of crucial importance. The establishment of a database and understanding of the defense mechanisms under combined stress conditions would be meaningful for the development of resistant and tolerant plants. It is clear that leaf pathogens show great tolerance to salinity stress and result in pathogenicity in crop plants. We noticed that regulation of the stomata through biochemical applications and some effort with the upregulation of the minor gene expressions indirectly involved with the defense mechanisms could be a great way to increase the defense metabolites without interfering with quality parameters. In this review, we selected wheat as a model plant and Zymoseptoria tritici as a model leaf pathogen to evaluate the defense mechanisms under saline conditions through physiological, biochemical, and molecular pathways and suggested various ways to generate tolerant and resistant cereal plants.

IRAPs in Combination with Highly Informative ISSRs Confer Effective Potentials for Genetic Diversity and Fidelity Assessment in Rhododendron

The species belonging to the Rhododendron genus are well-known for their colorful corolla. Molecular marker systems have the potential to elucidate genetic diversity as well as to assess genetic fidelity in rhododendrons. In the present study, the reverse transcription domains of long terminal repeat retrotransposons were cloned from rhododendrons and used to develop an inter-retrotransposon amplified polymorphism (IRAP) marker system. Subsequently, 198 polymorphic loci were generated from the IRAP and inter-simple sequence repeat (ISSR) markers, of which 119 were derived from the IRAP markers. It was shown that in rhododendrons, IRAP markers were superior to the ISSRs in some polymorphic parameters, such as the average number of polymorphic loci (14.88 versus 13.17). The combination of the IRAP and ISSR systems was more discriminative for detecting 46 rhododendron accessions than each of the systems on their own. Furthermore, IRAP markers demonstrated more efficiency in genetic fidelity detection of in-vitro-grown R. bailiense Y.P.Ma, C.Q.Zhang and D.F.Chamb, an endangered species recently recorded in Guizhzhou Province, China. The available evidence revealed the distinct properties of IRAP and ISSR markers in the rhododendron-associated applications, and highlighted the availability of highly informative ISSR and IRAP markers in the evaluation of genetic diversity and genetic fidelity of rhododendrons, which may facilitate preservation and genetic breeding of rhododendron plants.

Major proliferation of transposable elements shaped the genome of the soybean rust pathogen Phakopsora pachyrhizi

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

The activity of transposable elements (TEs) contributes significantly to pathogen genome evolution. TEs often destabilize genome integrity but may also confer adaptive variation in pathogenicity or resistance traits. De-repression of epigenetically silenced TEs often initiates bursts of transposition activity that may be counteracted by purifying selection and genome defenses. However, how these forces interact to determine the expansion routes of TEs within a pathogen species remains largely unknown. Here, we analyzed a set of 19 telomere-to-telomere genomes of the fungal wheat pathogen Zymoseptoria tritici. Phylogenetic reconstruction and ancestral state estimates of individual TE families revealed that TEs have undergone distinct activation and repression periods resulting in highly uneven copy numbers between genomes of the same species. Most TEs are clustered in gene poor niches, indicating strong purifying selection against insertions near coding sequences, or as a consequence of insertion site preferences. TE families with high copy numbers have low sequence divergence and strong signatures of defense mechanisms (i.e., RIP). In contrast, small non-autonomous TEs (i.e., MITEs) are less impacted by defense mechanisms and are often located in close proximity to genes. Individual TE families have experienced multiple distinct burst events that generated many nearly identical copies. We found that a Copia element burst was initiated from recent copies inserted substantially closer to genes compared to older copies. Overall, TE bursts tended to initiate from copies in GC-rich niches that escaped inactivation by genomic defenses. Our work shows how specific genomic environments features provide triggers for TE proliferation in pathogen genomes.

The pangenome of the wheat pathogen Pyrenophora tritici-repentis reveals novel transposons associated with necrotrophic effectors ToxA and ToxB

BACKGROUND

In fungal plant pathogens, genome rearrangements followed by selection pressure for adaptive traits have facilitated the co-evolutionary arms race between hosts and their pathogens. Pyrenophora tritici-repentis (Ptr) has emerged recently as a foliar pathogen of wheat worldwide and its populations consist of isolates that vary in their ability to produce combinations of different necrotrophic effectors. These effectors play vital roles in disease development. Here, we sequenced the genomes of a global collection (40 isolates) of Ptr to gain insights into its gene content and genome rearrangements.

RESULTS

A comparative genome analysis revealed an open pangenome, with an abundance of accessory genes (~ 57%) reflecting Ptr's adaptability. A clear distinction between pathogenic and non-pathogenic genomes was observed in size, gene content, and phylogenetic relatedness. Chromosomal rearrangements and structural organization, specifically around effector coding genes, were detailed using long-read assemblies (PacBio RS II) generated in this work in addition to previously assembled genomes. We also discovered the involvement of large mobile elements associated with Ptr's effectors: ToxA, the gene encoding for the necrosis effector, was found as a single copy within a 143-kb 'Starship' transposon (dubbed 'Horizon') with a clearly defined target site and target site duplications. 'Horizon' was located on different chromosomes in different isolates, indicating mobility, and the previously described ToxhAT transposon (responsible for horizontal transfer of ToxA) was nested within this newly identified Starship. Additionally, ToxB, the gene encoding the chlorosis effector, was clustered as three copies on a 294-kb element, which is likely a different putative 'Starship' (dubbed 'Icarus') in a ToxB-producing isolate. ToxB and its putative transposon were missing from the ToxB non-coding reference isolate, but the homolog toxb and 'Icarus' were both present in a different non-coding isolate. This suggests that ToxB may have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption of ToxB vertical inheritance. Finally, the genome architecture of Ptr was defined as 'one-compartment' based on calculated gene distances and evolutionary rates.

CONCLUSIONS

These findings together reflect on the highly plastic nature of the Ptr genome which has likely helped to drive its worldwide adaptation and has illuminated the involvement of giant transposons in facilitating the evolution of virulence in Ptr.

Comparative genomics reveals a dynamic genome evolution in the ectomycorrhizal milk-cap (Lactarius) mushrooms

Ectomycorrhizal fungi play a key role in forests by establishing mutualistic symbioses with woody plants. Genome analyses have identified conserved symbiosis-related traits among ectomycorrhizal fungal species, but the molecular mechanisms underlying host specificity remain poorly known. We sequenced and compared the genomes of seven species of milk-cap fungi (Lactarius, Russulales) with contrasting host specificity. We also compared these genomes with those of symbiotic and saprotrophic Russulales species, aiming to identify genes involved in their ecology and host specificity. The size of Lactarius genomes is significantly larger than other Russulales species, owing to a massive accumulation of transposable elements and duplication of dispensable genes. As expected, their repertoire of genes coding for plant cell wall-degrading enzymes is restricted, but they retained a substantial set of genes involved in microbial cell wall degradation. Notably, Lactarius species showed a striking expansion of genes encoding proteases, such as secreted ectomycorrhiza-induced sedolisins. A high copy number of genes coding for small secreted LysM proteins and Lactarius-specific lectins were detected, which may be linked to host specificity. This study revealed a large diversity in the genome landscapes and gene repertoires within Russulaceae. The known host specificity of Lactarius symbionts may be related to mycorrhiza-induced species-specific genes, including secreted sedolisins.

Genomic Variations and Mutational Events Associated with Plant-Pathogen Interactions

Simple Summary

Plants, unlike animals, do not have defender cells or an adaptive immune system. Instead, plants rely on each cell’s innate immunity and systemic signals emitted from infection sites. On the other hand, not all plants, even within the same species, are genetically identical, and their genetic backgrounds determine how well they respond to stress factors. Through evolution, plants have acquired various defense mechanisms that play important roles in the never-ending fight between plants and pathogens. Genetic variation in relation to plant disease resistance can thus be contextualized to provide new insights into these defense mechanisms and evolutionary processes that lead to resistance to pathogens. By focusing on genetic variations and mutational events linked with plant–pathogen interactions, the paper explores how genome compartments facilitate plant and pathogen evolutionary processes.

Abstract

Phytopathologists are actively researching the molecular basis of plant–pathogen interactions. The mechanisms of responses to pathogens have been studied extensively in model crop plant species and natural populations. Today, with the rapid expansion of genomic technologies such as DNA sequencing, transcriptomics, proteomics, and metabolomics, as well as the development of new methods and protocols, data analysis, and bioinformatics, it is now possible to assess the role of genetic variation in plant–microbe interactions and to understand the underlying molecular mechanisms of plant defense and microbe pathogenicity with ever-greater resolution and accuracy. Genetic variation is an important force in evolution that enables organisms to survive in stressful environments. Moreover, understanding the role of genetic variation and mutational events is essential for crop breeders to produce improved cultivars. This review focuses on genetic variations and mutational events associated with plant–pathogen interactions and discusses how these genome compartments enhance plants’ and pathogens’ evolutionary processes.

Population genomics of transposable element activation in the highly repressive genome of an agricultural pathogen

The activity of transposable elements (TEs) can be an important driver of genetic diversity with TE-mediated mutations having a wide range of fitness consequences. To avoid deleterious effects of TE activity, some fungi have evolved highly sophisticated genomic defences to reduce TE proliferation across the genome. Repeat-induced point mutation (RIP) is a fungal-specific TE defence mechanism efficiently targeting duplicated sequences. The rapid accumulation of RIPs is expected to deactivate TEs over the course of a few generations. The evolutionary dynamics of TEs at the population level in a species with highly repressive genome defences is poorly understood. Here, we analyse 366 whole-genome sequences of Parastagonospora nodorum, a fungal pathogen of wheat with efficient RIP. A global population genomics analysis revealed high levels of genetic diversity and signs of frequent sexual recombination. Contrary to expectations for a species with RIP, we identified recent TE activity in multiple populations. The TE composition and copy numbers showed little divergence among global populations regardless of the demographic history. Miniature inverted-repeat transposable elements (MITEs) and terminal repeat retrotransposons in miniature (TRIMs) were largely underlying recent intra-species TE expansions. We inferred RIP footprints in individual TE families and found that recently active, high-copy TEs have possibly evaded genomic defences. We find no evidence that recent positive selection acted on TE-mediated mutations rather that purifying selection maintained new TE insertions at low insertion frequencies in populations. Our findings highlight the complex evolutionary equilibria established by the joint action of TE activity, selection and genomic repression.

Dynamics of transposable elements in recently diverged fungal pathogens: lineage-specific transposable element content and efficiency of genome defenses

Transposable elements (TEs) impact genome plasticity, architecture, and evolution in fungal plant pathogens. The wide range of TE content observed in fungal genomes reflects diverse efficacy of host-genome defense mechanisms that can counter-balance TE expansion and spread. Closely related species can harbor drastically different TE repertoires. The evolution of fungal effectors, which are crucial determinants of pathogenicity, has been linked to the activity of TEs in pathogen genomes. Here, we describe how TEs have shaped genome evolution of the fungal wheat pathogen Zymoseptoria tritici and four closely related species. We compared de novo TE annotations and repeat-induced point mutation signatures in 26 genomes from the Zymoseptoria species-complex. Then, we assessed the relative insertion ages of TEs using a comparative genomics approach. Finally, we explored the impact of TE insertions on genome architecture and plasticity. The 26 genomes of Zymoseptoria species reflect different TE dynamics with a majority of recent insertions. TEs associate with accessory genome compartments, with chromosomal rearrangements, with gene presence/absence variation, and with effectors in all Zymoseptoria species. We find that the extent of RIP-like signatures varies among Z. tritici genomes compared to genomes of the sister species. The detection of a reduction of RIP-like signatures and TE recent insertions in Z. tritici reflects ongoing but still moderate TE mobility.

Comparative Genomics: Insights on the Pathogenicity and Lifestyle of Rhizoctonia solani

Proper management of agricultural disease is important to ensure sustainable food security. Staple food crops like rice, wheat, cereals, and other cash crops hold great export value for countries. Ensuring proper supply is critical; hence any biotic or abiotic factors contributing to the shortfall in yield of these crops should be alleviated. Rhizoctonia solani is a major biotic factor that results in yield losses in many agriculturally important crops. This paper focuses on genome informatics of our Malaysian Draft R. solani AG1-IA, and the comparative genomics (inter- and intra- AG) with four AGs including China AG1-IA (AG1-IA_KB317705.1), AG1-IB, AG3, and AG8. The genomic content of repeat elements, transposable elements (TEs), syntenic genomic blocks, functions of protein-coding genes as well as core orthologous genic information that underlies R. solani’s pathogenicity strategy were investigated. Our analyses show that all studied AGs have low content and varying profiles of TEs. All AGs were dominant for Class I TE, much like other basidiomycete pathogens. All AGs demonstrate dominance in Glycoside Hydrolase protein-coding gene assignments suggesting its importance in infiltration and infection of host. Our profiling also provides a basis for further investigation on lack of correlation observed between number of pathogenicity and enzyme-related genes with host range. Despite being grouped within the same AG with China AG1-IA, our Draft AG1-IA exhibits differences in terms of protein-coding gene proportions and classifications. This implies that strains from similar AG do not necessarily have to retain similar proportions and classification of TE but must have the necessary arsenal to enable successful infiltration and colonization of host. In a larger perspective, all the studied AGs essentially share core genes that are generally involved in adhesion, penetration, and host colonization. However, the different infiltration strategies will depend on the level of host resilience where this is clearly exhibited by the gene sets encoded for the process of infiltration, infection, and protection from host.

Conversion of DNA Sequences: From a Transposable Element to a Tandem Repeat or to a Gene

Eukaryotic genomes are rich in repetitive DNA sequences grouped in two classes regarding their genomic organization: tandem repeats and dispersed repeats. In tandem repeats, copies of a short DNA sequence are positioned one after another within the genome, while in dispersed repeats, these copies are randomly distributed. In this review we provide evidence that both tandem and dispersed repeats can have a similar organization, which leads us to suggest an update to their classification based on the sequence features, concretely regarding the presence or absence of retrotransposons/transposon specific domains. In addition, we analyze several studies that show that a repetitive element can be remodeled into repetitive non-coding or coding sequences, suggesting (1) an evolutionary relationship among DNA sequences, and (2) that the evolution of the genomes involved frequent repetitive sequence reshuffling, a process that we have designated as a "DNA remodeling mechanism". The alternative classification of the repetitive DNA sequences here proposed will provide a novel theoretical framework that recognizes the importance of DNA remodeling for the evolution and plasticity of eukaryotic genomes.

Transposable Elements Adaptive Role in Genome Plasticity, Pathogenicity and Evolution in Fungal Phytopathogens

Transposable elements (TEs) are agents of genetic variability in phytopathogens as they are a source of adaptive evolution through genome diversification. Although many studies have uncovered information on TEs, the exact mechanism behind TE-induced changes within the genome remains poorly understood. Furthermore, convergent trends towards bigger genomes, emergence of novel genes and gain or loss of genes implicate a TE-regulated genome plasticity of fungal phytopathogens. TEs are able to alter gene expression by revamping the cis-regulatory elements or recruiting epigenetic control. Recent findings show that TEs recruit epigenetic control on the expression of effector genes as part of the coordinated infection strategy. In addition to genome plasticity and diversity, fungal pathogenicity is an area of economic concern. A survey of TE distribution suggests that their proximity to pathogenicity genes TEs may act as sites for emergence of novel pathogenicity factors via nucleotide changes and expansion or reduction of the gene family. Through a systematic survey of literature, we were able to conclude that the role of TEs in fungi is wide: ranging from genome plasticity, pathogenicity to adaptive behavior in evolution. This review also identifies the gaps in knowledge that requires further elucidation for a better understanding of TEs' contribution to genome architecture and versatility.

Genome and evolution of the arbuscular mycorrhizal fungus Diversispora epigaea (formerly Glomus versiforme) and its bacterial endosymbionts

Arbuscular mycorrhizal (AM) fungi form endosymbioses with most plants, and they themselves are hosts for Mollicutes/Mycoplasma-related endobacteria (MRE). Despite their significance, genomic information for AM fungi and their MRE are relatively sparse, which hinders our understanding of their biology and evolution. We assembled the genomes of the AM fungus Diversispora epigaea (formerly Glomus versiforme) and its MRE and performed comparative genomics and evolutionary analyses. The D. epigaea genome showed a pattern of substantial gene duplication and differential evolution of gene families, including glycosyltransferase family 25, whose activities are exclusively lipopolysaccharide biosynthesis. Genes acquired by horizontal transfer from bacteria possibly function in defense against foreign DNA or viruses. The MRE population was diverse, with multiple genomes displaying characteristics of differential evolution and encoding many MRE-specific genes as well as genes of AM fungal origin. Gene family expansion in D. epigaea may enhance adaptation to both external and internal environments, such as expansion of kinases for signal transduction upon external stimuli and expansion of nucleoside salvage pathway genes potentially for competition with MRE, whose genomes lack purine and pyrimidine biosynthetic pathways. Collectively, this metagenome provides high-quality references and begins to reveal the diversity within AM fungi and their MRE.

Meiotic drive of female-inherited supernumerary chromosomes in a pathogenic fungus

Meiosis is a key cellular process of sexual reproduction that includes pairing of homologous sequences. In many species however, meiosis can also involve the segregation of supernumerary chromosomes, which can lack a homolog. How these unpaired chromosomes undergo meiosis is largely unknown. In this study we investigated chromosome segregation during meiosis in the haploid fungus Zymoseptoria tritici that possesses a large complement of supernumerary chromosomes. We used isogenic whole chromosome deletion strains to compare meiotic transmission of chromosomes when paired and unpaired. Unpaired chromosomes inherited from the male parent as well as paired supernumerary chromosomes in general showed Mendelian inheritance. In contrast, unpaired chromosomes inherited from the female parent showed non-Mendelian inheritance but were amplified and transmitted to all meiotic products. We concluded that the supernumerary chromosomes of Z. tritici show a meiotic drive and propose an additional feedback mechanism during meiosis, which initiates amplification of unpaired female-inherited chromosomes.

The Landscape of Repetitive Elements in the Refined Genome of Chilli Anthracnose Fungus Colletotrichum truncatum

The ascomycete fungus Colletotrichum truncatum is a major phytopathogen with a broad host range which causes anthracnose disease of chilli. The genome sequencing of this fungus led to the discovery of functional categories of genes that may play important roles in fungal pathogenicity. However, the presence of gaps in C. truncatum draft assembly prevented the accurate prediction of repetitive elements, which are the key players to determine the genome architecture and drive evolution and host adaptation. We re-sequenced its genome using single-molecule real-time (SMRT) sequencing technology to obtain a refined assembly with lesser and smaller gaps and ambiguities. This enabled us to study its genome architecture by characterising the repetitive sequences like transposable elements (TEs) and simple sequence repeats (SSRs), which constituted 4.9 and 0.38% of the assembled genome, respectively. The comparative analysis among different Colletotrichum species revealed the extensive repeat rich regions, dominated by Gypsy superfamily of long terminal repeats (LTRs), and the differential composition of SSRs in their genomes. Our study revealed a recent burst of LTR amplification in C. truncatum, C. higginsianum, and C. scovillei. TEs in C. truncatum were significantly associated with secretome, effectors and genes in secondary metabolism clusters. Some of the TE families in C. truncatum showed cytosine to thymine transitions indicative of repeat-induced point mutation (RIP). C. orbiculare and C. graminicola showed strong signatures of RIP across their genomes and "two-speed" genomes with extensive AT-rich and gene-sparse regions. Comparative genomic analyses of Colletotrichum species provided an insight into the species-specific SSR profiles. The SSRs in the coding and non-coding regions of the genome revealed the composition of trinucleotide repeat motifs in exons with potential to alter the translated protein structure through amino acid repeats. This is the first genome-wide study of TEs and SSRs in C. truncatum and their comparative analysis with six other Colletotrichum species, which would serve as a useful resource for future research to get insights into the potential role of TEs in genome expansion and evolution of Colletotrichum fungi and for development of SSR-based molecular markers for population genomic studies.

Extraordinary Genome Instability and Widespread Chromosome Rearrangements During Vegetative Growth

The haploid genome of the pathogenic fungus Zymoseptoria tritici is contained on "core" and "accessory" chromosomes. While 13 core chromosomes are found in all strains, as many as eight accessory chromosomes show presence/absence variation and rearrangements among field isolates. The factors influencing these presence/absence polymorphisms are so far unknown. We investigated chromosome stability using experimental evolution, karyotyping, and genome sequencing. We report extremely high and variable rates of accessory chromosome loss during mitotic propagation in vitro and in planta Spontaneous chromosome loss was observed in 2 to >50% of cells during 4 weeks of incubation. Similar rates of chromosome loss in the closely related Zymoseptoria ardabiliae suggest that this extreme chromosome dynamic is a conserved phenomenon in the genus. Elevating the incubation temperature greatly increases instability of accessory and even core chromosomes, causing severe rearrangements involving telomere fusion and chromosome breakage. Chromosome losses do not affect the fitness of Zymoseptoria tritici in vitro, but some lead to increased virulence, suggesting an adaptive role of this extraordinary chromosome instability.

The evolution of genomic and epigenomic features in two Pleurotus fungi

Pleurotus tuoliensis (Bailinggu, designated Pt) and P. eryngii var. eryngii (Xingbaogu, designated Pe) are highly valued edible mushrooms. We report de novo assemblies of high-quality genomes for both mushrooms based on PacBio RS II sequencing and annotation of all identified genes. A comparative genomics analysis between Pt and Pe with P. ostreatus as an outgroup taxon revealed extensive genomic divergence between the two mushroom genomes primarily due to the rapid gain of taxon-specific genes and disruption of synteny in either taxon. The re-appraised phylogenetic relationship between Pt and Pe at the genome-wide level validates earlier proposals to designate Pt as an independent species. Variation of the identified wood-decay-related gene content can largely explain the variable adaptation and host specificity of the two mushrooms. On the basis of the two assembled genome sequences, methylomes and the regulatory roles of DNA methylation in gene expression were characterized and compared. The genome, methylome and transcriptome data of these two important mushrooms will provide valuable information for advancing our understanding of the evolution of Pleurotus and related genera and for facilitating genome- and epigenome-based strategies for mushroom breeding.

Analysis of Transposable Elements in Coccidioides Species

Coccidioides immitis and C. posadasii are primary pathogenic fungi that cause disease in immunologically-normal animals and people. The organism is found exclusively in arid regions of the Southwestern United States, Mexico, and South America, but not in other parts of the world. This study is a detailed analysis of the transposable elements (TE) in Coccidioides spp. As is common in most fungi, Class I and Class II transposons were identified and the LTR Gypsy superfamily is the most common. The minority of Coccidioides Gypsy transposons contained regions highly homologous to polyprotein domains. Phylogenetic analysis of the integrase and reverse transcriptase sequences revealed that many, but not all, of the Gypsy reverse transcriptase and integrase domains clustered by species suggesting extensive transposition after speciation of the two Coccidiodies spp. The TEs were clustered and the distribution is enriched for the ends on contigs. Analysis of gene expression data from C. immitis found that protein-coding genes within 1 kB of hAT or Gypsy TEs were poorly expressed. The expression of C. posadasii genes within 1 kB of Gypsy TEs was also significantly lower compared to all genes but the difference in expression was smaller than C. immitis. C. posadasii orthologs of C. immitis Gyspsy-associated genes were also likely to be TE-associated. In both C. immitis and C. posadasii the TEs were preferentially associated with genes annotated with protein kinase gene ontology terms. These observations suggest that TE may play a role in influencing gene expression in Coccidioides spp. Our hope is that these bioinformatic studies of the potential TE influence on expression and evolution of Coccidioides will prompt the development of testable hypotheses to better understand the role of TEs in the biology and gene regulation of Coccidioides spp.

Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici

We conducted a comprehensive analysis of virulence in the fungal wheat pathogen Zymoseptoria tritici using quantitative trait locus (QTL) mapping. High-throughput phenotyping based on automated image analysis allowed the measurement of pathogen virulence on a scale and with a precision that was not previously possible. Across two mapping populations encompassing more than 520 progeny, 540 710 pycnidia were counted and their sizes and grey values were measured. A significant correlation was found between pycnidia size and both spore size and number. Precise measurements of percentage leaf area covered by lesions provided a quantitative measure of host damage. Combining these large and accurate phenotypic datasets with a dense panel of restriction site-associated DNA sequencing (RADseq) genetic markers enabled us to genetically dissect pathogen virulence into components related to host damage and those related to pathogen reproduction. We showed that different components of virulence can be under separate genetic control. Large- and small-effect QTLs were identified for all traits, with some QTLs specific to mapping populations, cultivars and traits and other QTLs shared among traits within the same mapping population. We associated the presence of four accessory chromosomes with small, but significant, increases in several virulence traits, providing the first evidence for a meaningful function associated with accessory chromosomes in this organism. A large-effect QTL involved in host specialization was identified on chromosome 7, leading to the identification of candidate genes having a large effect on virulence.

Comparative genomics of Coniophora olivacea reveals different patterns of genome expansion in Boletales

Background

Coniophora olivacea is a basidiomycete fungus belonging to the order Boletales that produces brown-rot decay on dead wood of conifers. The Boletales order comprises a diverse group of species including saprotrophs and ectomycorrhizal fungi that show important differences in genome size.

Results

In this study we report the 39.07-megabase (Mb) draft genome assembly and annotation of C. olivacea. A total of 14,928 genes were annotated, including 470 putatively secreted proteins enriched in functions involved in lignocellulose degradation. Using similarity clustering and protein structure prediction we identified a new family of 10 putative lytic polysaccharide monooxygenase genes. This family is conserved in basidiomycota and lacks of previous functional annotation. Further analyses showed that C. olivacea has a low repetitive genome, with 2.91% of repeats and a restrained content of transposable elements (TEs). The annotation of TEs in four related Boletales yielded important differences in repeat content, ranging from 3.94 to 41.17% of the genome size. The distribution of insertion ages of LTR-retrotransposons showed that differential expansions of these repetitive elements have shaped the genome architecture of Boletales over the last 60 million years.

Conclusions

Coniophora olivacea has a small, compact genome that shows macrosynteny with Coniophora puteana. The functional annotation revealed the enzymatic signature of a canonical brown-rot. The annotation and comparative genomics of transposable elements uncovered their particular contraction in the Coniophora genera, highlighting their role in the differential genome expansions found in Boletales species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4243-z) contains supplementary material, which is available to authorized users.

Plasticity of the MFS1 Promoter Leads to Multidrug Resistance in the Wheat Pathogen Zymoseptoria tritici

The ascomycete Zymoseptoria tritici is the causal agent of Septoria leaf blotch on wheat. Disease control relies mainly on resistant wheat cultivars and on fungicide applications. The fungus displays a high potential to circumvent both methods. Resistance against all unisite fungicides has been observed over decades. A different type of resistance has emerged among wild populations with multidrug-resistant (MDR) strains. Active fungicide efflux through overexpression of the major facilitator gene MFS1 explains this emerging resistance mechanism. Applying a bulk-progeny sequencing approach, we identified in this study a 519-bp long terminal repeat (LTR) insert in the MFS1 promoter, a relic of a retrotransposon cosegregating with the MDR phenotype. Through gene replacement, we show the insert as a mutation responsible for MFS1 overexpression and the MDR phenotype. Besides this type I insert, we found two different types of promoter inserts in more recent MDR strains. Type I and type II inserts harbor potential transcription factor binding sites, but not the type III insert. Interestingly, all three inserts correspond to repeated elements present at different genomic locations in either IPO323 or other Z. tritici strains. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici and which contribute to its adaptive potential. IMPORTANCE Disease control through fungicides remains an important means to protect crops from fungal diseases and to secure the harvest. Plant-pathogenic fungi, especially Zymoseptoria tritici, have developed resistance against most currently used active ingredients, reducing or abolishing their efficacy. While target site modification is the most common resistance mechanism against single modes of action, active efflux of multiple drugs is an emerging phenomenon in fungal populations reducing additionally fungicides' efficacy in multidrug-resistant strains. We have investigated the mutations responsible for increased drug efflux in Z. tritici field strains. Our study reveals that three different insertions of repeated elements in the same promoter lead to multidrug resistance in Z. tritici. The target gene encodes the membrane transporter MFS1 responsible for drug efflux, with the promoter inserts inducing its overexpression. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici.

Gapless genome assembly of Colletotrichum higginsianum reveals chromosome structure and association of transposable elements with secondary metabolite gene clusters

Background

The ascomycete fungus Colletotrichum higginsianum causes anthracnose disease of brassica crops and the model plant Arabidopsis thaliana. Previous versions of the genome sequence were highly fragmented, causing errors in the prediction of protein-coding genes and preventing the analysis of repetitive sequences and genome architecture.

Results

Here, we re-sequenced the genome using single-molecule real-time (SMRT) sequencing technology and, in combination with optical map data, this provided a gapless assembly of all twelve chromosomes except for the ribosomal DNA repeat cluster on chromosome 7. The more accurate gene annotation made possible by this new assembly revealed a large repertoire of secondary metabolism (SM) key genes (89) and putative biosynthetic pathways (77 SM gene clusters). The two mini-chromosomes differed from the ten core chromosomes in being repeat- and AT-rich and gene-poor but were significantly enriched with genes encoding putative secreted effector proteins. Transposable elements (TEs) were found to occupy 7% of the genome by length. Certain TE families showed a statistically significant association with effector genes and SM cluster genes and were transcriptionally active at particular stages of fungal development. All 24 subtelomeres were found to contain one of three highly-conserved repeat elements which, by providing sites for homologous recombination, were probably instrumental in four segmental duplications.

Conclusion

The gapless genome of C. higginsianum provides access to repeat-rich regions that were previously poorly assembled, notably the mini-chromosomes and subtelomeres, and allowed prediction of the complete SM gene repertoire. It also provides insights into the potential role of TEs in gene and genome evolution and host adaptation in this asexual pathogen.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4083-x) contains supplementary material, which is available to authorized users.

Karyotype Variability in Plant-Pathogenic Fungi

Recent advances in genetic and molecular technologies gradually paved the way for the transition from traditional fungal karyotyping to more comprehensive chromosome biology studies. Extensive chromosomal polymorphisms largely resulting from chromosomal rearrangements (CRs) are widely documented in fungal genomes. These extraordinary CRs in fungi generate substantial genome plasticity compared to other eukaryotic organisms. Here, we review the most recent findings on fungal CRs and their underlying mechanisms and discuss the functional consequences of CRs for adaptation, fungal evolution, host range, and pathogenicity of fungal plant pathogens in the context of chromosome biology. In addition to a complement of permanent chromosomes called core chromosomes, the genomes of many fungal pathogens comprise distinct unstable chromosomes called dispensable chromosomes (DCs) that also contribute to chromosome polymorphisms. Compared to the core chromosomes, the structural features of DCs usually differ for gene density, GC content, housekeeping genes, and recombination frequency. Despite their dispensability for normal growth and development, DCs have important biological roles with respect to pathogenicity in some fungi but not in others. Therefore, their evolutionary origin is also reviewed in relation to overall fungal physiology and pathogenicity.

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

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

Trichoderma reesei complete genome sequence, repeat-induced point mutation, and partitioning of CAZyme gene clusters

BACKGROUND

Trichoderma reesei (Ascomycota, Pezizomycotina) QM6a is a model fungus for a broad spectrum of physiological phenomena, including plant cell wall degradation, industrial production of enzymes, light responses, conidiation, sexual development, polyketide biosynthesis, and plant-fungal interactions. The genomes of QM6a and its high enzyme-producing mutants have been sequenced by second-generation-sequencing methods and are publicly available from the Joint Genome Institute. While these genome sequences have offered useful information for genomic and transcriptomic studies, their limitations and especially their short read lengths make them poorly suited for some particular biological problems, including assembly, genome-wide determination of chromosome architecture, and genetic modification or engineering.

RESULTS

We integrated Pacific Biosciences and Illumina sequencing platforms for the highest-quality genome assembly yet achieved, revealing seven telomere-to-telomere chromosomes (34,922,528 bp; 10877 genes) with 1630 newly predicted genes and >1.5 Mb of new sequences. Most new sequences are located on AT-rich blocks, including 7 centromeres, 14 subtelomeres, and 2329 interspersed AT-rich blocks. The seven QM6a centromeres separately consist of 24 conserved repeats and 37 putative centromere-encoded genes. These findings open up a new perspective for future centromere and chromosome architecture studies. Next, we demonstrate that sexual crossing readily induced cytosine-to-thymine point mutations on both tandem and unlinked duplicated sequences. We also show by bioinformatic analysis that T. reesei has evolved a robust repeat-induced point mutation (RIP) system to accumulate AT-rich sequences, with longer AT-rich blocks having more RIP mutations. The widespread distribution of AT-rich blocks correlates genome-wide partitions with gene clusters, explaining why clustering of genes has been reported to not influence gene expression in T. reesei.

CONCLUSION

Compartmentation of ancestral gene clusters by AT-rich blocks might promote flexibilities that are evolutionarily advantageous in this fungus' soil habitats and other natural environments. Our analyses, together with the complete genome sequence, provide a better blueprint for biotechnological and industrial applications.

de Jong T, de Vries I, Dietrich R, Farmer AD, Fortes Fereira C, Garcia S, Guzman M, Hamelin RC, Lindquist EA, Mehrabi R, Quiros O, Schmutz J, Shapiro H, Reynolds E, Scalliet G, Souza M, Stergiopoulos I, Van der Lee TAJ, De Wit PJGM, Zapater MF, Zwiers LH, Grigoriev IV, Goodwin SB, Kema GHJ. Combating a Global Threat to a Clonal Crop: Banana Black Sigatoka Pathogen Pseudocercospora fijiensis (Synonym Mycosphaerella fijiensis) Genomes Reveal Clues for Disease Control

Black Sigatoka or black leaf streak disease, caused by the Dothideomycete fungus Pseudocercospora fijiensis (previously: Mycosphaerella fijiensis), is the most significant foliar disease of banana worldwide. Due to the lack of effective host resistance, management of this disease requires frequent fungicide applications, which greatly increase the economic and environmental costs to produce banana. Weekly applications in most banana plantations lead to rapid evolution of fungicide-resistant strains within populations causing disease-control failures throughout the world. Given its extremely high economic importance, two strains of P. fijiensis were sequenced and assembled with the aid of a new genetic linkage map. The 74-Mb genome of P. fijiensis is massively expanded by LTR retrotransposons, making it the largest genome within the Dothideomycetes. Melting-curve assays suggest that the genomes of two closely related members of the Sigatoka disease complex, P. eumusae and P. musae, also are expanded. Electrophoretic karyotyping and analyses of molecular markers in P. fijiensis field populations showed chromosome-length polymorphisms and high genetic diversity. Genetic differentiation was also detected using neutral markers, suggesting strong selection with limited gene flow at the studied geographic scale. Frequencies of fungicide resistance in fungicide-treated plantations were much higher than those in untreated wild-type P. fijiensis populations. A homologue of the Cladosporium fulvum Avr4 effector, PfAvr4, was identified in the P. fijiensis genome. Infiltration of the purified PfAVR4 protein into leaves of the resistant banana variety Calcutta 4 resulted in a hypersensitive-like response. This result suggests that Calcutta 4 could carry an unknown resistance gene recognizing PfAVR4. Besides adding to our understanding of the overall Dothideomycete genome structures, the P. fijiensis genome will aid in developing fungicide treatment schedules to combat this pathogen and in improving the efficiency of banana breeding programs.

Transposable Elements versus the Fungal Genome: Impact on Whole-Genome Architecture and Transcriptional Profiles

Transposable elements (TEs) are exceptional contributors to eukaryotic genome diversity. Their ubiquitous presence impacts the genomes of nearly all species and mediates genome evolution by causing mutations and chromosomal rearrangements and by modulating gene expression. We performed an exhaustive analysis of the TE content in 18 fungal genomes, including strains of the same species and species of the same genera. Our results depicted a scenario of exceptional variability, with species having 0.02 to 29.8% of their genome consisting of transposable elements. A detailed analysis performed on two strains of Pleurotus ostreatus uncovered a genome that is populated mainly by Class I elements, especially LTR-retrotransposons amplified in recent bursts from 0 to 2 million years (My) ago. The preferential accumulation of TEs in clusters led to the presence of genomic regions that lacked intra- and inter-specific conservation. In addition, we investigated the effect of TE insertions on the expression of their nearby upstream and downstream genes. Our results showed that an important number of genes under TE influence are significantly repressed, with stronger repression when genes are localized within transposon clusters. Our transcriptional analysis performed in four additional fungal models revealed that this TE-mediated silencing was present only in species with active cytosine methylation machinery. We hypothesize that this phenomenon is related to epigenetic defense mechanisms that are aimed to suppress TE expression and control their proliferation.

Transposable elements (TEs) are enigmatic genetic units that have played important roles in the evolution of eukaryotic genomes. Since their discovery in the 1950s, they have gained increasing attention and are known today as active genome modelers in multiple species. Although these elements have been widely studied in plants, much less is known about their occurrence and impact on the fungal kingdom. Using a diverse set of basidiomycete and ascomycete fungi, we quantified and characterized a huge diversity of DNA and RNA transposable elements, and we identified species that had 0.02 to 29.8% of their genomes occupied by transposable elements. In addition, using our basidiomycete model Pleurotus ostreatus, we demonstrated how TE insertions produced detrimental effects on the expression of upstream and downstream genes, which were downregulated compared with the control groups. This silencing mechanism was present in the basidiomycetes tested but exhibited a patchy distribution in ascomycetes, and might be related to specific genome defense mechanisms that control transposon proliferation. This finding reveals the broader impact of transposable elements in fungi. In addition to their importance as long-term evolutionary forces, they play major roles in the more dynamic transcriptome regulation of certain species.

Next-generation re-sequencing as a tool for rapid bioinformatic screening of presence and absence of genes and accessory chromosomes across isolates of Zymoseptoria tritici

The wheat pathogen Zymoseptoria tritici possesses a large number of accessory chromosomes that may be present or absent in its genome. The genome of the reference isolate IPO323 has been assembled to a very high standard and contains 21 full length chromosome sequences, 8 of which represent accessory chromosomes. The IPO323 reference, when combined with low-cost next-generation sequencing and bioinformatics, can be used as a powerful tool to assess the presence or absence of accessory chromosomes. We present an outline of a range of bioinformatics techniques that can be applied to the analysis of presence-absence variation among accessory chromosomes across 13 novel isolates of Z. tritici.

The impact of Septoria tritici Blotch disease on wheat: An EU perspective

•

Zymospetoria tritici is a threat to wheat production in the EU Z.t.’s plastic genome increases the potential severity of this threat in the future.

•

Climate change may also affect the risk from Z.t.

•

We estimate the spore numbers produced by Z.t. during each infection cycle.

•

We calculate 1) the economic value of wheat in the three main EU producers 2) the cost of and economic return for fungicide treatment of wheat vs Z.t.

Zymoseptoria tritici is the causal agent of one of the European Union’s most devastating foliar diseases of wheat: Septoria tritici Blotch (STB). It is also a notable pathogen of wheat grown in temperate climates throughout the world. In this commentary, we highlight the importance of STB on wheat in the EU. To better understand STB, it is necessary to consider the host crop, the fungal pathogen and their shared environment. Here, we consider the fungus per se and its interaction with its host and then focus on a more agricultural overview of the impact STB on wheat. We consider the climatic and weather factors which influence its spread and severity, allude to the agricultural practices which may mitigate or enhance its impact on crop yields, and evaluate the economic importance of wheat as a food and animal feed crop in the UK and EU. Finally, we estimate the cost of STB disease to EU agriculture.

Cell biology of Zymoseptoria tritici: Pathogen cell organization and wheat infection

•

Cell biology of the infection begins to shed light on the host–pathogen interaction.

•

The cell biology of the fungal pathogen is highly understudied.

•

Intensified cell biology research promises new fungicide targets and will help mode-of-action studies.

Cell biological research in the wheat pathogen Zymoseptoria tritici (formerly Mycosphaerella graminicola) has led to a good understanding of the histology of the infection process. Expression profiling and bioinformatic approaches, combined with molecular studies on signaling pathways, effectors and potential necrosis factors provides first insight into the complex interplay between the host and the pathogen. Cell biological studies will help to further our understanding of the infection strategy of the fungus. The cellular organization and intracellular dynamics of the fungus itself is largely unexplored. Insight into essential cellular processes within the pathogen will expand our knowledge of the basic biology of Z. tritici, thereby providing putative new anti-fungal targets.

Comparative genome analysis and genome evolution of members of the magnaporthaceae family of fungi

BACKGROUND

Magnaporthaceae, a family of ascomycetes, includes three fungi of great economic importance that cause disease in cereal and turf grasses: Magnaporthe oryzae (rice blast), Gaeumannomyces graminis var. tritici (take-all disease), and Magnaporthe poae (summer patch disease). Recently, the sequenced and assembled genomes for these three fungi were reported. Here, the genomes were compared for orthologous genes in order to identified genes that are unique to the Magnaporthaceae family of fungi. In addition, ortholog clustering was used to identify a core proteome for the Magnaporthaceae, which was examined for diversifying and purifying selection and evidence of two-speed genome evolution.

RESULTS

A genome-scale comparative study was conducted across 74 fungal genomes to identify clusters of orthologous genes unique to the three Magnaporthaceae species as well as species specific genes. We found 1149 clusters that were unique to the Magnaporthaceae family of fungi with 295 of those containing genes from all three species. Gene clusters involved in metabolic and enzymatic activities were highly represented in the Magnaporthaceae specific clusters. Also highly represented in the Magnaporthaceae specific clusters as well as in the species specific genes were transcriptional regulators. In addition, we examined the relationship between gene evolution and distance to repetitive elements found in the genome. No correlations between diversifying or purifying selection and distance to repetitive elements or an increased rate of evolution in secreted and small secreted proteins were observed.

CONCLUSIONS

Taken together, these data show that at the genome level, there is no evidence to suggest multi-speed genome evolution or that proximity to repetitive elements play a role in diversification of genes.

Microsatellites in Pursuit of Microbial Genome Evolution

Microsatellites or short sequence repeats are widespread genetic markers which are hypermutable 1-6 bp long short nucleotide motifs. Significantly, their applications in genetics are extensive due to their ceaseless mutational degree, widespread length variations and hypermutability skills. These features make them useful in determining the driving forces of evolution by using powerful molecular techniques. Consequently, revealing important questions, for example, what is the significance of these abundant sequences in DNA, what are their roles in genomic evolution? The answers of these important questions are hidden in the ways these short motifs contributed in altering the microbial genomes since the origin of life. Even though their size ranges from 1 -to- 6 bases, these repeats are becoming one of the most popular genetic probes in determining their associations and phylogenetic relationships in closely related genomes. Currently, they have been widely used in molecular genetics, biotechnology and evolutionary biology. However, due to limited knowledge; there is a significant gap in research and lack of information concerning hypermutational mechanisms. These mechanisms play a key role in microsatellite loci point mutations and phase variations. This review will extend the understandings of impacts and contributions of microsatellite in genomic evolution and their universal applications in microbiology.

Vaginal Microbiome Characterization of Nellore Cattle Using Metagenomic Analysis

Understanding of microbial communities inhabiting cattle vaginal tract may lead to a better comprehension of bovine physiology and reproductive health being of great economic interest. Up to date, studies involving cattle microbiota are focused on the gastrointestinal tract, and little is known about the vaginal microbiota. This study aimed to investigate the vaginal microbiome in Nellore cattle, heifers and cows, pregnant and non-pregnant, using a culture independent approach. The main bacterial phyla found were Firmicutes (~40–50%), Bacteroidetes (~15–25%) and Proteobacteria (~5–25%), in addition to ~10–20% of non-classified bacteria. 45–55% of the samples were represented by only ten OTUs: Aeribacillus, Bacteroides, Clostridium, Ruminococcus, Rikenella, Alistipes, Bacillus, Eubacterium, Prevotella and non-classified bacteria. Interestingly, microbiota from all 20 animals could be grouped according to the respiratory metabolism of the main OTUs found, creating three groups of vaginal microbiota in cattle. Archaeal samples were dominated by the Methanobrevibacter genus (Euryarchaeota, ~55–70%). Ascomycota was the main fungal phylum (~80–95%) and Mycosphaerella the most abundant genus (~70–85%). Hormonal influence was not clear, but a tendency for the reduction of bacterial and increase of archaeal populations in pregnant animals was observed. Eukaryotes did not vary significantly between pregnant and non-pregnant animals, but tended to be more abundant on cows than on heifers. The present work describes a great microbial variability in the vaginal community among the evaluated animals and groups (heifers and cows, pregnant and non-pregnant), which is significantly different from the findings previously reported using culture dependent methods, pointing out the need for further studies on this issue. The microbiome found also indicates that the vaginal colonization appears to be influenced by the gastrointestinal community.

Histone modifications rather than the novel regional centromeres of Zymoseptoria tritici distinguish core and accessory chromosomes

BACKGROUND

Supernumerary chromosomes have been found in many organisms. In fungi, these "accessory" or "dispensable" chromosomes are present at different frequencies in populations and are usually characterized by higher repetitive DNA content and lower gene density when compared to the core chromosomes. In the reference strain of the wheat pathogen, Zymoseptoria tritici, eight discrete accessory chromosomes have been found. So far, no functional role has been assigned to these chromosomes; however, they have existed as separate entities in the karyotypes of Zymoseptoria species over evolutionary time. In this study, we addressed what-if anything-distinguishes the chromatin of accessory chromosomes from core chromosomes. We used chromatin immunoprecipitation combined with high-throughput sequencing ("ChIP-seq") of DNA associated with the centromere-specific histone H3, CENP-A (CenH3), to identify centromeric DNA, and ChIP-seq with antibodies against dimethylated H3K4, trimethylated H3K9 and trimethylated H3K27 to determine the relative distribution and proportion of euchromatin, obligate and facultative heterochromatin, respectively.

RESULTS

Centromeres of the eight accessory chromosomes have the same sequence composition and structure as centromeres of the 13 core chromosomes and they are of similar length. Unlike those of most other fungi, Z. tritici centromeres are not composed entirely of repetitive DNA; some centromeres contain only unique DNA sequences, and bona fide expressed genes are located in regions enriched with CenH3. By fluorescence microscopy, we showed that centromeres of Z. tritici do not cluster into a single chromocenter during interphase. We found dramatically higher enrichment of H3K9me3 and H3K27me3 on the accessory chromosomes, consistent with the twofold higher proportion of repetitive DNA and poorly transcribed genes. In contrast, no single histone modification tested here correlated with the distribution of centromeric nucleosomes.

CONCLUSIONS

All centromeres are similar in length and composed of a mixture of unique and repeat DNA, and most contain actively transcribed genes. Centromeres, subtelomeric regions or telomere repeat length cannot account for the differences in transfer fidelity between core and accessory chromosomes, but accessory chromosomes are greatly enriched in nucleosomes with H3K27 trimethylation. Genes on accessory chromosomes appear to be silenced by trimethylation of H3K9 and H3K27.

Genome Sequences of Three Phytopathogenic Species of the Magnaporthaceae Family of Fungi

Magnaporthaceae is a family of ascomycetes that includes three fungi of great economic importance: Magnaporthe oryzae, Gaeumannomyces graminis var. tritici, and Magnaporthe poae. These three fungi cause widespread disease and loss in cereal and grass crops, including rice blast disease (M. oryzae), take-all disease in wheat and other grasses (G. graminis), and summer patch disease in turf grasses (M. poae). Here, we present the finished genome sequence for M. oryzae and draft sequences for M. poae and G. graminis var. tritici. We used multiple technologies to sequence and annotate the genomes of M. oryzae, M. poae, and G. graminis var. tritici. The M. oryzae genome is now finished to seven chromosomes whereas M. poae and G. graminis var. tritici are sequenced to 40.0× and 25.0× coverage respectively. Gene models were developed by the use of multiple computational techniques and further supported by RNAseq data. In addition, we performed preliminary analysis of genome architecture and repetitive element DNA.

Single-Molecule Real-Time Sequencing Combined with Optical Mapping Yields Completely Finished Fungal Genome

Next-generation sequencing (NGS) technologies have increased the scalability, speed, and resolution of genomic sequencing and, thus, have revolutionized genomic studies. However, eukaryotic genome sequencing initiatives typically yield considerably fragmented genome assemblies. Here, we assessed various state-of-the-art sequencing and assembly strategies in order to produce a contiguous and complete eukaryotic genome assembly, focusing on the filamentous fungus Verticillium dahliae. Compared with Illumina-based assemblies of the V. dahliae genome, hybrid assemblies that also include PacBio-generated long reads establish superior contiguity. Intriguingly, provided that sufficient sequence depth is reached, assemblies solely based on PacBio reads outperform hybrid assemblies and even result in fully assembled chromosomes. Furthermore, the addition of optical map data allowed us to produce a gapless and complete V. dahliae genome assembly of the expected eight chromosomes from telomere to telomere. Consequently, we can now study genomic regions that were previously not assembled or poorly assembled, including regions that are populated by repetitive sequences, such as transposons, allowing us to fully appreciate an organism’s biological complexity. Our data show that a combination of PacBio-generated long reads and optical mapping can be used to generate complete and gapless assemblies of fungal genomes.

Studying whole-genome sequences has become an important aspect of biological research. The advent of next-generation sequencing (NGS) technologies has nowadays brought genomic science within reach of most research laboratories, including those that study nonmodel organisms. However, most genome sequencing initiatives typically yield (highly) fragmented genome assemblies. Nevertheless, considerable relevant information related to genome structure and evolution is likely hidden in those nonassembled regions. Here, we investigated a diverse set of strategies to obtain gapless genome assemblies, using the genome of a typical ascomycete fungus as the template. Eventually, we were able to show that a combination of PacBio-generated long reads and optical mapping yields a gapless telomere-to-telomere genome assembly, allowing in-depth genome analyses to facilitate functional studies into an organism’s biology.

Single-Molecule Real-Time Sequencing Combined with Optical Mapping Yields Completely Finished Fungal Genome

Next-generation sequencing (NGS) technologies have increased the scalability, speed, and resolution of genomic sequencing and, thus, have revolutionized genomic studies. However, eukaryotic genome sequencing initiatives typically yield considerably fragmented genome assemblies. Here, we assessed various state-of-the-art sequencing and assembly strategies in order to produce a contiguous and complete eukaryotic genome assembly, focusing on the filamentous fungus Verticillium dahliae. Compared with Illumina-based assemblies of the V. dahliae genome, hybrid assemblies that also include PacBio-generated long reads establish superior contiguity. Intriguingly, provided that sufficient sequence depth is reached, assemblies solely based on PacBio reads outperform hybrid assemblies and even result in fully assembled chromosomes. Furthermore, the addition of optical map data allowed us to produce a gapless and complete V. dahliae genome assembly of the expected eight chromosomes from telomere to telomere. Consequently, we can now study genomic regions that were previously not assembled or poorly assembled, including regions that are populated by repetitive sequences, such as transposons, allowing us to fully appreciate an organism's biological complexity. Our data show that a combination of PacBio-generated long reads and optical mapping can be used to generate complete and gapless assemblies of fungal genomes.

IMPORTANCE

Studying whole-genome sequences has become an important aspect of biological research. The advent of next-generation sequencing (NGS) technologies has nowadays brought genomic science within reach of most research laboratories, including those that study nonmodel organisms. However, most genome sequencing initiatives typically yield (highly) fragmented genome assemblies. Nevertheless, considerable relevant information related to genome structure and evolution is likely hidden in those nonassembled regions. Here, we investigated a diverse set of strategies to obtain gapless genome assemblies, using the genome of a typical ascomycete fungus as the template. Eventually, we were able to show that a combination of PacBio-generated long reads and optical mapping yields a gapless telomere-to-telomere genome assembly, allowing in-depth genome analyses to facilitate functional studies into an organism's biology.

Chromatin-based control of effector gene expression in plant-associated fungi

Plant-associated fungi often present in their genome areas enriched in repeat sequences and effector genes, the latter being specifically induced in planta. The location of effector genes in regions enriched in repeats has been shown to have an impact on adaptability of fungi but could also provide for tight control of effector gene expression through chromatin-based regulation. The distribution of two repressive histone marks was shown to be an important regulatory layer in two fungal species with different lifestyles. Chromatin-based control of effector gene expression is likely to provide an evolutionary advantage by preventing the expression of genes not needed during vegetative growth and allow for a massive concerted expression at particular time-points of plant infection.

RNA-seq-Based Gene Annotation and Comparative Genomics of Four Fungal Grass Pathogens in the Genus Zymoseptoria Identify Novel Orphan Genes and Species-Specific Invasions of Transposable Elements

The fungal pathogen Zymoseptoria tritici (synonym Mycosphaerella graminicola) is a prominent pathogen of wheat. The reference genome of the isolate IPO323 is one of the best-assembled eukaryotic genomes and encodes more than 10,000 predicted genes. However, a large proportion of the previously annotated gene models are incomplete, with either no start or no stop codons. The availability of RNA-seq data allows better predictions of gene structure. We here used two different RNA-seq datasets, de novo transcriptome assemblies, homology-based comparisons, and trained ab initio gene callers to generate a new gene annotation of Z. tritici IPO323. The annotation pipeline was also applied to re-sequenced genomes of three closely related species of Z. tritici: Z. pseudotritici, Z. ardabiliae, and Z. brevis. Comparative analyses of the predicted gene models using the four Zymoseptoria species revealed sets of species-specific orphan genes enriched with putative pathogenicity-related genes encoding small secreted proteins that may play essential roles in virulence and host specificity. De novo repeat identification allowed us to show that few families of transposable elements are shared between Zymoseptoria species while we observe many species-specific invasions and expansions. The annotation data presented here provide a high-quality resource for future studies of Z. tritici and its sister species and provide detailed insight into gene and genome evolution of fungal plant pathogens.

Chromatin analyses of Zymoseptoria tritici: Methods for chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq)

The presence or absence of specific transcription factors, chromatin remodeling machineries, chromatin modification enzymes, post-translational histone modifications and histone variants all play crucial roles in the regulation of pathogenicity genes. Chromatin immunoprecipitation (ChIP) followed by high-throughput sequencing (ChIP-seq) provides an important tool to study genome-wide protein-DNA interactions to help understand gene regulation in the context of native chromatin. ChIP-seq is a convenient in vivo technique to identify, map and characterize occupancy of specific DNA fragments with proteins against which specific antibodies exist or which can be epitope-tagged in vivo. We optimized existing ChIP protocols for use in the wheat pathogen Zymoseptoria tritici and closely related sister species. Here, we provide a detailed method, underscoring which aspects of the technique are organism-specific. Library preparation for Illumina sequencing is described, as this is currently the most widely used ChIP-seq method. One approach for the analysis and visualization of representative sequence is described; improved tools for these analyses are constantly being developed. Using ChIP-seq with antibodies against H3K4me2, which is considered a mark for euchromatin or H3K9me3 and H3K27me3, which are considered marks for heterochromatin, the overall distribution of euchromatin and heterochromatin in the genome of Z. tritici can be determined. Our ChIP-seq protocol was also successfully applied to Z. tritici strains with high levels of melanization or aberrant colony morphology, and to different species of the genus (Z. ardabiliae and Z. pseudotritici), suggesting that our technique is robust. The methods described here provide a powerful framework to study new aspects of chromatin biology and gene regulation in this prominent wheat pathogen.