Genome evolution in filamentous plant pathogens: why bigger can be better

Advances in understanding obligate biotrophy in rust fungi

Contents Summary 1190 I. Introduction 1190 II. Rust fungi: a diverse and serious threat to agriculture 1191 III. The different facets of rust life cycles and unresolved questions about their evolution 1191 IV. The biology of rust infection 1192 V. Rusts in the genomics era: the ever-expanding list of candidate effector genes 1195 VI. Functional characterization of rust effectors 1197 VII. Putting rusts to sleep: Pucciniales research outlooks 1201 Acknowledgements 1202 References 1202 SUMMARY: Rust fungi (Pucciniales) are the largest group of plant pathogens and represent one of the most devastating threats to agricultural crops worldwide. Despite the economic importance of these highly specialized pathogens, many aspects of their biology remain obscure, largely because rust fungi are obligate biotrophs. The rise of genomics and advances in high-throughput sequencing technology have presented new options for identifying candidate effector genes involved in pathogenicity mechanisms of rust fungi. Transcriptome analysis and integrated bioinformatics tools have led to the identification of key genetic determinants of host susceptibility to infection by rusts. Thousands of genes encoding secreted proteins highly expressed during host infection have been reported for different rust species, which represents significant potential towards understanding rust effector function. Recent high-throughput in planta expression screen approaches (effectoromics) have pushed the field ahead even further towards predicting high-priority effectors and identifying avirulence genes. These new insights into rust effector biology promise to inform future research and spur the development of effective and sustainable strategies for managing rust diseases.

Smut infection of perennial hosts: the genome and the transcriptome of the Brassicaceae smut fungus Thecaphora thlaspeos reveal functionally conserved and novel effectors

Biotrophic fungal plant pathogens can balance their virulence and form intricate relationships with their hosts. Sometimes, this leads to systemic host colonization over long time scales without macroscopic symptoms. However, how plant-pathogenic endophytes manage to establish their sustained systemic infection remains largely unknown. Here, we present a genomic and transcriptomic analysis of Thecaphora thlaspeos. This relative of the well studied grass smut Ustilago maydis is the only smut fungus adapted to Brassicaceae hosts. Its ability to overwinter with perennial hosts and its systemic plant infection including roots are unique characteristics among smut fungi. The T. thlaspeos genome was assembled to the chromosome level. It is a typical smut genome in terms of size and genome characteristics. In silico prediction of candidate effector genes revealed common smut effector proteins and unique members. For three candidates, we have functionally demonstrated effector activity. One of these, TtTue1, suggests a potential link to cold acclimation. On the plant side, we found evidence for a typical immune response as it is present in other infection systems, despite the absence of any macroscopic symptoms during infection. Our findings suggest that T. thlaspeos distinctly balances its virulence during biotrophic growth ultimately allowing for long-lived infection of its perennial hosts.

A Novel Effector Gene SCRE2 Contributes to Full Virulence of Ustilaginoidea virens to Rice

Ustilaginoidea virens, the causal agent of rice false smut (RFS), has become one of the most devastating rice pathogens worldwide. As a group of essential virulence factors, the effectors in the filamentous fungus might play central roles in the interaction between plants and pathogens. However, little is known about the roles of individual effectors in U. virens virulence. In this study, we identified and characterized a small secreted cysteine-rich effector, SCRE2, in U. virens. SCRE2 was first confirmed as an effector through yeast secretion, protein localization and translocation assays, as well as its expression pattern during U. virens infection. Transient expression of SCRE2 in Nicotiana benthamiana suppressed necrosis-like defense symptoms triggered by the mammalian BAX and oomycete elicitin INF1 proteins. The ability of SCRE2 to inhibit immunity-associated responses in N. benthamiana, including elicitor-triggered cell death and oxidative burst, is further defined to a small peptide region SCRE2^68-85 through expressing a series of truncated proteins. Convincingly, ectopic expression of SCRE2 in the transgenic rice cells significantly inhibited pathogen-associated molecular pattern-triggered immunity including flg22- and chitin-induced defense gene expression and oxidative burst. Furthermore, the scre2 knockout mutant generated by the CRISPR/Cas9 system greatly attenuated in U. virens virulence to rice. Collectively, this study indicates that the effector SCRE2 is able to inhibit plant immunity and is required for full virulence of U. virens.

Evolution of Mutator transposable elements across eukaryotic diversity

Background

Mutator-like elements (MULEs) are a significant superfamily of DNA transposons on account of their: (i) great transpositional activity and propensity for insertion in or near gene sequences, (ii) their consequent high mutagenic capacity, and, (iii) their tendency to acquire host gene fragments. Consequently, MULEs are important genetic tools and represent a key study system for research into host-transposon interactions. Yet, while several studies have focused on the impacts of MULEs on crop and fungus genomes, their evolution remains poorly explored.

Results

We perform comprehensive bioinformatic and phylogenetic analyses to address currently available MULE diversity and reconstruct evolution for the group. For this, we mine MULEs from online databases, and combine search results with available transposase sequences retrieved from previously published studies. Our analyses uncover two entirely new MULE clades that contain elements almost entirely restricted to arthropod hosts, considerably expanding the set of MULEs known from this group, suggesting that many additional MULEs may await discovery from further arthropod genomes. In several cases, close relationships occur between MULEs recovered from distantly related host organisms, suggesting that horizontal transfer events may have played an important role in the evolution of the group. However, it is apparent that MULEs from plants remain separate from MULEs identified from other host groups. MULE structure varies considerably across phylogeny, and TIR length is shown to vary greatly both within and between MULE groups. Our phylogeny suggests that MULE diversity is clustered in well-supported groups, typically according to host taxonomy. With reference to this, we make suggestions on how MULE diversity can be partitioned to provide a robust taxonomic framework.

Conclusions

Our study represents a considerable advance in the understanding of MULE diversity, host range and evolution, and provides a taxonomic framework for the classification of further MULE elements that await discovery. Our findings also raise a number of questions relating to MULE biology, suggesting that this group will provide a rich avenue for future study.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0153-8) contains supplementary material, which is available to authorized users.

Comparative genomics of plant pathogenic Botrytis species with distinct host specificity

BACKGROUND

Fungi of the genus Botrytis (presently containing ~ 35 species) are able to infect more than 1400 different plant species and cause losses in a wide range of crops of economic importance. The best studied species is B. cinerea, which has a broad host range and is one of the best studied necrotrophic plant pathogenic fungi. Most other Botrytis spp. have a narrow host range and have been studied in less detail. To characterize genomic variation among different representatives of Botrytis spp., we sequenced and annotated the draft genomes of nine Botrytis species: B. calthae, B. convoluta, B. elliptica, B. galanthina, B. hyacinthi, B. narcissicola, B. paeoniae, B. porri and B. tulipae.

RESULTS

Bioinformatics and comparative genomics tools were applied to determine a core of 7668 shared protein families in all Botrytis species, which grouped them in two distinct phylogenetic clades. The secretome of all nine Botrytis spp. was similar in number (ranging from 716 to 784 predicted proteins). A detailed analysis of the molecular functions of the secretome revealed that shared activities were highly similar. Orthologs to effectors functionally studied in B. cinerea were also present in the other Botrytis species. A complex pattern of presence/absence of secondary metabolite biosynthetic key enzymes was observed.

CONCLUSIONS

Comparative genomics of Botrytis show that overall, species share the main signatures and protein families in the secreted proteins, and of known effectors. Our study provides leads to study host range determinants in the genus Botrytis and provides a stepping stone to elucidate the roles of effector candidates in the infection process of these species.

Mitotic Recombination and Rapid Genome Evolution in the Invasive Forest Pathogen Phytophthora ramorum

Invasive alien species often have reduced genetic diversity and must adapt to new environments. Given the success of many invasions, this is sometimes called the genetic paradox of invasion. Phytophthora ramorum is invasive, limited to asexual reproduction within four lineages, and presumed clonal. It is responsible for sudden oak death in the United States, sudden larch death in Europe, and ramorum blight in North America and Europe. We sequenced the genomes of 107 isolates to determine how this pathogen can overcome the invasion paradox. Mitotic recombination (MR) associated with transposons and low gene density has generated runs of homozygosity (ROH) affecting 2,698 genes, resulting in novel genotypic diversity within the lineages. One ROH enriched in effectors was fixed in the NA1 lineage. An independent ROH affected the same scaffold in the EU1 lineage, suggesting an MR hot spot and a selection target. Differences in host infection between EU1 isolates with and without the ROH suggest that they may differ in aggressiveness. Non-core regions (not shared by all lineages) had signatures of accelerated evolution and were enriched in putative pathogenicity genes and transposons. There was a striking pattern of gene loss, including all effectors, in the non-core EU2 genome. Positive selection was observed in 8.0% of RxLR and 18.8% of Crinkler effector genes compared with 0.9% of the core eukaryotic gene set. We conclude that the P. ramorum lineages are diverging via a rapidly evolving non-core genome and that the invasive asexual lineages are not clonal, but display genotypic diversity caused by MR.IMPORTANCE Alien species are often successful invaders in new environments, despite the introduction of a few isolates with a reduced genetic pool. This is called the genetic paradox of invasion. We found two mechanisms by which the invasive forest pathogen causing sudden oak and sudden larch death can evolve. Extensive mitotic recombination producing runs of homozygosity generates genotypic diversity even in the absence of sexual reproduction, and rapid turnover of genes in the non-core, or nonessential portion of genome not shared by all isolates, allows pathogenicity genes to evolve rapidly or be eliminated while retaining essential genes. Mitotic recombination events occur in genomic hot spots, resulting in similar ROH patterns in different isolates or groups; one ROH, independently generated in two different groups, was enriched in pathogenicity genes and may be a target for selection. This provides important insights into the evolution of invasive alien pathogens and their potential for adaptation and future persistence.

Intercellular cooperation in a fungal plant pathogen facilitates host colonization

Cooperation is associated with major transitions in evolution such as the emergence of multicellularity. It is central to the evolution of many complex traits in nature, including growth and virulence in pathogenic bacteria. Whether cells of multicellular parasites function cooperatively during infection remains, however, largely unknown. Here, we show that hyphal cells of the fungal pathogen Sclerotinia sclerotiorum reprogram toward division of labor to facilitate the colonization of host plants. Using global transcriptome sequencing, we reveal that gene expression patterns diverge markedly in cells at the center and apex of hyphae during Arabidopsis thaliana colonization compared with in vitro growth. We reconstructed a genome-scale metabolic model for S. sclerotiorum and used flux balance analysis to demonstrate metabolic heterogeneity supporting division of labor between hyphal cells. Accordingly, continuity between the central and apical compartments of invasive hyphae was required for optimal growth in planta Using a multicell model of fungal hyphae, we show that this cooperative functioning enhances fungal growth predominantly during host colonization. Our work identifies cooperation in fungal hyphae as a mechanism emerging at the multicellular level to support host colonization and virulence.

Analysis of Epichloë festucae small secreted proteins in the interaction with Lolium perenne

Epichloë festucae is an endophyte of the agriculturally important perennial ryegrass. This species systemically colonises the aerial tissues of this host where its growth is tightly regulated thereby maintaining a mutualistic symbiotic interaction. Recent studies have suggested that small secreted proteins, termed effectors, play a vital role in the suppression of host defence responses. To date only a few effectors with important roles in mutualistic interactions have been described. Here we make use of the fully assembled E. festucae genome and EffectorP to generate a suite of 141 effector candidates. These were analysed with respect to their genome location and expression profiles in planta and in several symbiosis-defective mutants. We found an association between effector candidates and a class of transposable elements known as MITEs, but no correlation with other dynamic features of the E. festucae genome, such as transposable element-rich regions. Three effector candidates and a small GPI-anchored protein were chosen for functional analysis based on their high expression in planta compared to in culture and their differential regulation in symbiosis defective E. festucae mutants. All three candidate effector proteins were shown to possess a functional signal peptide and two could be detected in the extracellular medium by western blotting. Localization of the effector candidates in planta suggests that they are not translocated into the plant cell, but rather, are localized in the apoplastic space or are attached to the cell wall. Deletion and overexpression of the effector candidates, as well as the putative GPI-anchored protein, did not affect the plant growth phenotype or restrict growth of E. festucae mutants in planta. These results indicate that these proteins are either not required for the interaction at the observed life stages or that there is redundancy between effectors expressed by E. festucae.

Genome-wide DNA methylation and transcriptomic profiles in the lifestyle strategies and asexual development of the forest fungal pathogen Heterobasidion parviporum

Heterobasidion parviporum is the most devastating fungal pathogen of conifer forests in Northern Europe. The fungus has dual life strategies, necrotrophy on living trees and saprotrophy on dead woods. DNA cytosine methylation is an important epigenetic modification in eukaryotic organisms. Our presumption is that the lifestyle transition and asexual development in H. parviporum could be driven by epigenetic effects. Involvements of DNA methylation in the regulation of aforementioned processes have never been studied thus far. RNA-seq identified lists of highly induced genes enriched in carbohydrate-active enzymes during necrotrophic interaction with host trees and saprotrophic sawdust growth. It also highlighted signaling- and transcription factor-related genes potentially associated with the transition of saprotrophic to necrotrophic lifestyle and groups of primary cellular activities throughout asexual development. Whole-genome bisulfite sequencing revealed that DNA methylation displayed pronounced preference in CpG dinucleotide context across the genome and mostly targeted transposable element (TE)-rich regions. TE methylation level demonstrated a strong negative correlation with TE expression, reinforcing the protective function of DNA methylation in fungal genome stability. Small groups of genes putatively subject to methylation transcriptional regulation in response to saprotrophic and necrotrophic growth in comparison with free-living mycelia were also explored. Our study reported on the first methylome map of a forest pathogen. Analysis of transcriptome and methylome variations associated with asexual development and different lifestyle strategies provided further understanding of basic biological processes in H. parviporum. More importantly, our work raised additional potential roles of DNA methylation in fungi apart from controlling the proliferation of TEs.

Comparative genomic analysis revealed rapid differentiation in the pathogenicity-related gene repertoires between Pyricularia oryzae and Pyricularia penniseti isolated from a Pennisetum grass

BACKGROUND

A number of Pyricularia species are known to infect different grass species. In the case of Pyricularia oryzae (syn. Magnaporthe oryzae), distinct populations are known to be adapted to a wide variety of grass hosts, including rice, wheat and many other grasses. The genome sizes of Pyricularia species are typical for filamentous ascomycete fungi [~ 40 Mbp for P. oryzae, and ~ 45 Mbp for P. grisea]. Genome plasticity, mediated in part by deletions promoted by recombination between repetitive elements [Genome Res 26:1091-1100, 2016, Nat Rev Microbiol 10:417-430,2012] and transposable elements [Annu Rev Phytopathol 55:483-503,2017] contributes to host adaptation. Therefore, comparisons of genome structure of individual species will provide insight into the evolution of host specificity. However, except for the P. oryzae subgroup, little is known about the gene content or genome organization of other Pyricularia species, such as those infecting Pennisetum grasses.

RESULTS

Here, we report the genome sequence of P. penniseti strain P1609 isolated from a Pennisetum grass (JUJUNCAO) using PacBio SMRT sequencing technology. Phylogenomic analysis of 28 Magnaporthales species and 5 non-Magnaporthales species indicated that P1609 belongs to a Pyricularia subclade, which is genetically distant from P. oryzae. Comparative genomic analysis revealed that the pathogenicity-related gene repertoires had diverged between P1609 and the P. oryzae strain 70-15, including the known avirulence genes, other putative secreted proteins, as well as some other predicted Pathogen-Host Interaction (PHI) genes. Genomic sequence comparison also identified many genomic rearrangements relative to P. oryzae.

CONCLUSION

Our results suggested that the genomic sequence of the P. penniseti P1609 could be a useful resource for the genetic study of the Pennisetum-infecting Pyricularia species and provide new insight into evolution of pathogen genomes during host adaptation.

A downy mildew effector evades recognition by polymorphism of expression and subcellular localization

Pathogen co-evolution with plants involves selection for evasion of host surveillance systems. The oomycete Hyaloperonospora arabidopsidis (Hpa) causes downy mildew on Arabidopsis, and race-specific interactions between Arabidopsis accessions and Hpa isolates fit the gene-for-gene model in which host resistance or susceptibility are determined by matching pairs of plant Resistance (R) genes and pathogen Avirulence (AVR) genes. Arabidopsis Col-0 carries R gene RPP4 that confers resistance to Hpa isolates Emoy2 and Emwa1, but its cognate recognized effector(s) were unknown. We report here the identification of the Emoy2 AVR effector gene recognized by RPP4 and show resistance-breaking isolates of Hpa on RPP4-containing Arabidopsis carry the alleles that either are not expressed, or show cytoplasmic instead of nuclear subcellular localization.

Plant pathogens have evolved to evade detection by their hosts. Here, Asai et al. show that virulent isolates of the oomycete Hyaloperonospora arabidopsidis can break resistance conferred by the Arabidopsis RPP4 resistance gene via variation in effector expression or subcellular localization.

A New Reference Genome Shows the One-Speed Genome Structure of the Barley Pathogen Ramularia collo-cygni

Ramularia leaf spot has recently emerged as a major threat to barley production world-wide, causing 25% yield loss in many barley growing regions. Here, we provide a new reference genome of the causal agent, the Dothideomycete Ramularia collo-cygni. The assembly of 32 Mb consists of 78 scaffolds. We used RNA-seq to identify 11,622 genes of which 1,303 and 282 are coding for predicted secreted proteins and putative effectors respectively.

The pathogen separated from its nearest sequenced relative, Zymoseptoria tritici ∼27 Ma. We calculated the divergence of the two species on protein level and see remarkably high synonymous and nonsynonymous divergence. Unlike in many other plant pathogens, the comparisons of transposable elements and gene distributions, show a very homogeneous genome for R. collo-cygni. We see no evidence for higher selective pressure on putative effectors or other secreted proteins and repetitive sequences are spread evenly across the scaffolds. These findings could be associated to the predominantly endophytic life-style of the pathogen. We hypothesize that R. collo-cygni only recently became pathogenic and that therefore its genome does not yet show the typical pathogen characteristics. Because of its high scaffold length and improved CDS annotations, our new reference sequence provides a valuable resource for the community for future comparative genomics and population genetics studies.

Arbuscular mycorrhizal fungi: intraspecific diversity and pangenomes

Contents Summary 1129 I. Introduction 1129 II. Intraspecific phenotypic variation and the plant host 1130 III. High inter-isolate genetic diversity in model AMF 1130 IV. Genome diversity within the model AM fungus Rhizophagus irregularis 1131 V. Pangenomes and the future of AMF ecological genomics 1131 Acknowledgements 1133 References 1133 SUMMARY: Arbuscular mycorrhizal fungi (AMF) are ubiquitous plant symbionts with an intriguing population biology. Conspecific AMF strains can vary substantially at the genetic and phenotypic levels, leading to direct and quantifiable variation in plant growth. Recent studies have shown that high intraspecific diversity is very common in AMF, and not only found in model species. Studies have also revealed how the phenotype of conspecific isolates varies depending on the plant host, highlighting the functional relevance of intraspecific phenotypic plasticity for the AMF ecology and mycorrhizal symbiosis. Recent work has also demonstrated that conspecific isolates of the model AMF Rhizophagus irregularis harbor large and highly variable pangenomes, highlighting the potential role of intraspecific genome diversity for the ecological adaptation of these symbionts.

Phytophthora methylomes are modulated by 6mA methyltransferases and associated with adaptive genome regions

BACKGROUND

Filamentous plant pathogen genomes often display a bipartite architecture with gene-sparse, repeat-rich compartments serving as a cradle for adaptive evolution. The extent to which this two-speed genome architecture is associated with genome-wide DNA modifications is unknown.

RESULTS

We show that the oomycetes Phytophthora infestans and Phytophthora sojae possess functional adenine N6-methylation (6mA) methyltransferases that modulate patterns of 6mA marks across the genome. In contrast, 5-methylcytosine could not be detected in these species. Methylated DNA IP sequencing (MeDIP-seq) of each species reveals 6mA is depleted around the transcription start sites (TSSs) and is associated with lowly expressed genes, particularly transposable elements. Genes occupying the gene-sparse regions have higher levels of 6mA in both genomes, possibly implicating the methylome in adaptive evolution. All six putative adenine methyltransferases from P. infestans and P. sojae, except PsDAMT2, display robust enzymatic activities. Surprisingly, single knockouts in P. sojae significantly reduce in vivo 6mA levels, indicating that the three enzymes are not fully redundant. MeDIP-seq of the psdamt3 mutant reveals uneven 6mA methylation reduction across genes, suggesting that PsDAMT3 may have a preference for gene body methylation after the TSS. Furthermore, transposable elements such as DNA elements are more active in the psdamt3 mutant. A large number of genes, particularly those from the adaptive genomic compartment, are differentially expressed.

CONCLUSIONS

Our findings provide evidence that 6mA modification is potentially an epigenetic mark in Phytophthora genomes, and complex patterns of 6mA methylation may be associated with adaptive evolution in these important plant pathogens.

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.

Evolutionary Genomics of Metchnikovella incurvata (Metchnikovellidae): An Early Branching Microsporidium

Metchnikovellids are highly specialized hyperparasites, which infect and reproduce inside gregarines (Apicomplexa) inhabiting marine invertebrates. Their phylogenetic affiliation was under constant discussion until recently, when analysis of the first near-complete metchnikovellid genome, that of Amphiamblys sp., placed it in a basal position with respect to most Microsporidia. Microsporidia are a highly diversified lineage of extremely reduced parasites related to Rozellida (Rozellosporidia = Rozellomycota = Cryptomycota) within the Holomycota clade of Opisthokonta. By sequencing DNA from a single-isolated infected gregarine cell we obtained an almost complete genome of a second metchnikovellid species, and the first one of a taxonomically described and well-documented species, Metchnikovella incurvata. Our phylogenomic analyses show that, despite being considerably divergent from each other, M. incurvata forms a monophyletic group with Amphiamplys sp., and confirm that metchnikovellids are one of the deep branches of Microsporidia. Comparative genomic analysis demonstrates that, like most Microsporidia, metchnikovellids lack mitochondrial genes involved in energy transduction and are thus incapable of synthesizing their own ATP via mitochondrial oxidative phosphorylation. They also lack the horizontally acquired ATP transporters widespread in most Microsporidia. We hypothesize that a family of mitochondrial carrier proteins evolved to transport ATP from the host into the metchnikovellid cell. We observe the progressive reduction of genes involved in DNA repair pathways along the evolutionary path of Microsporidia, which might explain, at least partly, the extremely high evolutionary rate of the most derived species. Our data also suggest that genome reduction and acquisition of novel genes co-occurred during the adaptation of Microsporidia to their hosts.

Human Pathogenic Entomophthorales

The pathogenic entomophthoralean fungi cause infection in insects and mammalian hosts. Basidiobolus and Conidiobolus species can be found in soil and insect, reptile, and amphibian droppings in tropical and subtropical areas. The life cycles of these fungi occur in these environments where infecting sticky conidia are developed. The infection is acquired by insect bite or contact with contaminated environments through open skin. Conidiobolus coronatus typically causes chronic rhinofacial disease in immunocompetent hosts, whereas some Conidiobolus species can be found in immunocompromised patients. Basidiobolus ranarum infection is restricted to subcutaneous tissues but may be involved in intestinal and disseminated infections. Its early diagnosis remains challenging due to clinical similarities to other intestinal diseases. Infected tissues characteristically display eosinophilic granulomas with the Splendore-Höeppli phenomenon. However, in immunocompromised patients, the above-mentioned inflammatory reaction is absent. Laboratory diagnosis includes wet mount, culture serological assays, and molecular methodologies. The management of entomophthoralean fungi relies on traditional antifungal therapies, such as potassium iodide (KI), amphotericin B, itraconazole, and ketoconazole, and surgery. These species are intrinsically resistant to some antifungals, prompting physicians to experiment with combinations of therapies. Research is needed to investigate the immunology of entomophthoralean fungi in infected hosts. The absence of an animal model and lack of funding severely limit research on these fungi.

Comparative genome analyses reveal sequence features reflecting distinct modes of host-adaptation between dicot and monocot powdery mildew

BACKGROUND

Powdery mildew (PM) is one of the most important and widespread plant diseases caused by biotrophic fungi. Notably, while monocot (grass) PM fungi exhibit high-level of host-specialization, many dicot PM fungi display a broad host range. To understand such distinct modes of host-adaptation, we sequenced the genomes of four dicot PM biotypes belonging to Golovinomyces cichoracearum or Oidium neolycopersici.

RESULTS

We compared genomes of the four dicot PM together with those of Blumeria graminis f.sp. hordei (both DH14 and RACE1 isolates), B. graminis f.sp. tritici, and Erysiphe necator infectious on barley, wheat and grapevine, respectively. We found that despite having a similar gene number (6620-6961), the PM genomes vary from 120 to 222 Mb in size. This high-level of genome size variation is indicative of highly differential transposon activities in the PM genomes. While the total number of genes in any given PM genome is only about half of that in the genomes of closely related ascomycete fungi, most (~ 93%) of the ascomycete core genes (ACGs) can be found in the PM genomes. Yet, 186 ACGs were found absent in at least two of the eight PM genomes, of which 35 are missing in some dicot PM biotypes, but present in the three monocot PM genomes, indicating remarkable, independent and perhaps ongoing gene loss in different PM lineages. Consistent with this, we found that only 4192 (3819 singleton) genes are shared by all the eight PM genomes, the remaining genes are lineage- or biotype-specific. Strikingly, whereas the three monocot PM genomes possess up to 661 genes encoding candidate secreted effector proteins (CSEPs) with families containing up to 38 members, all the five dicot PM fungi have only 116-175 genes encoding CSEPs with limited gene amplification.

CONCLUSIONS

Compared to monocot (grass) PM fungi, dicot PM fungi have a much smaller effectorome. This is consistent with their contrasting modes of host-adaption: while the monocot PM fungi show a high-level of host specialization, which may reflect an advanced host-pathogen arms race, the dicot PM fungi tend to practice polyphagy, which might have lessened selective pressure for escalating an with a particular host.

Genomic insights into host adaptation between the wheat stripe rust pathogen (Puccinia striiformis f. sp. tritici) and the barley stripe rust pathogen (Puccinia striiformis f. sp. hordei)

BACKGROUND

Plant fungal pathogens can rapidly evolve and adapt to new environmental conditions in response to sudden changes of host populations in agro-ecosystems. However, the genomic basis of their host adaptation, especially at the forma specialis level, remains unclear.

RESULTS

We sequenced two isolates each representing Puccinia striiformis f. sp. tritici (Pst) and P. striiformis f. sp. hordei (Psh), different formae speciales of the stripe rust fungus P. striiformis highly adapted to wheat and barley, respectively. The divergence of Pst and Psh, estimated to start 8.12 million years ago, has been driven by high nucleotide mutation rates. The high genomic variation within dikaryotic urediniospores of P. striiformis has provided raw genetic materials for genome evolution. No specific gene families have enriched in either isolate, but extensive gene loss events have occurred in both Pst and Psh after the divergence from their most recent common ancestor. A large number of isolate-specific genes were identified, with unique genomic features compared to the conserved genes, including 1) significantly shorter in length; 2) significantly less expressed; 3) significantly closer to transposable elements; and 4) redundant in pathways. The presence of specific genes in one isolate (or forma specialis) was resulted from the loss of the homologues in the other isolate (or forma specialis) by the replacements of transposable elements or losses of genomic fragments. In addition, different patterns and numbers of telomeric repeats were observed between the isolates.

CONCLUSIONS

Host adaptation of P. striiformis at the forma specialis level is a complex pathogenic trait, involving not only virulence-related genes but also other genes. Gene loss, which might be adaptive and driven by transposable element activities, provides genomic basis for host adaptation of different formae speciales of P. striiformis.

Interspecific Gene Exchange as a Driver of Adaptive Evolution in Fungi

Throughout evolutionary history in the kingdom Fungi, taxa have exchanged genetic information among species, as revealed in particular by analyses of genome sequences. In fungi, hybridization can occur by sexual mating or by fusion of vegetative structures giving rise to new species or leaving traces of introgression in the genome. Furthermore, gene exchange can occur by horizontal gene transfer between species and can even include organisms outside the kingdom Fungi. In several cases, interspecific gene exchange has been instrumental in rapid adaptive evolution of fungal species and has notably played a role in the emergence of new pathogens. Here we summarize mechanisms and examples of gene exchange in fungi with a particular focus on the genomic context. We emphasize the need for and potential of applying population genetic approaches to better understand the processes and the impact of interspecific gene exchange in rapid adaptive evolution and species diversification. The broad occurrence of gene exchange among fungal species challenges our species concepts in the kingdom Fungi.

Evidence for intragenic recombination and selective sweep in an effector gene of Phytophthora infestans

Effectors, a group of small proteins secreted by pathogens, play a critical role in the antagonistic interaction between plant hosts and pathogens through their dual functions in regulating host immune systems and pathogen infection capability. In this study, evolution in effector genes was investigated through population genetic analysis of Avr3a sequences generated from 96 Phytophthora infestans isolates collected from six locations representing a range of thermal variation and cropping systems in China. We found high genetic variation in the Avr3a gene resulting from diverse mechanisms extending beyond point mutations, frameshift, and defeated start and stop codons to intragenic recombination. A total of 51 nucleotide haplotypes encoding 38 amino acid isoforms were detected in the 96 full sequences with nucleotide diversity in the pathogen populations ranging from 0.007 to 0.023 (mean = 0.017). Although haplotype and nucleotide diversity were high, the effector gene was dominated by only three haplotypes. Evidence for a selective sweep was provided by (i) the population genetic differentiation (GST) of haplotypes being lower than the population differentiation (FST) of SSR marker loci; and (ii) negative values of Tajima's D and Fu's FS. Annual mean temperature in the collection sites was negatively correlated with the frequency of the virulent form (Avr3aEM), indicating Avr3a may be regulated by temperature. These results suggest that elevated air temperature due to global warming may hamper the development of pathogenicity traits in P. infestans and further study under confined thermal regimes may be required to confirm the hypothesis.

The Genome Biology of Effector Gene Evolution in Filamentous Plant Pathogens

Filamentous pathogens, including fungi and oomycetes, pose major threats to global food security. Crop pathogens cause damage by secreting effectors that manipulate the host to the pathogen's advantage. Genes encoding such effectors are among the most rapidly evolving genes in pathogen genomes. Here, we review how the major characteristics of the emergence, function, and regulation of effector genes are tightly linked to the genomic compartments where these genes are located in pathogen genomes. The presence of repetitive elements in these compartments is associated with elevated rates of point mutations and sequence rearrangements with a major impact on effector diversification. The expression of many effectors converges on an epigenetic control mediated by the presence of repetitive elements. Population genomics analyses showed that rapidly evolving pathogens show high rates of turnover at effector loci and display a mosaic in effector presence-absence polymorphism among strains. We conclude that effective pathogen containment strategies require a thorough understanding of the effector genome biology and the pathogen's potential for rapid adaptation.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

The impact of transposable elements in adaptive evolution

The growing knowledge about the influence of transposable elements (TEs) on (a) long-term genome and transcriptome evolution; (b) genomic, transcriptomic and epigenetic variation within populations; and (c) patterns of somatic genetic differences in individuals continues to spur the interest of evolutionary biologists in the role of TEs in adaptive evolution. As TEs can trigger a broad range of molecular variation in a population with potentially severe fitness and phenotypic consequences for individuals, different mechanisms evolved to keep TE activity in check, allowing for a dynamic interplay between the host, its TEs and the environment in evolution. Here, we review evidence for adaptive phenotypic changes associated with TEs and the basic molecular mechanisms by which the underlying genetic changes arise: (a) domestication, (b) exaptation, (c) host gene regulation, (d) TE-mediated formation of intronless gene copies-so-called retrogenes and (e) overall increased genome plasticity. Furthermore, we review and discuss how the stress-dependent incapacitation of defence mechanisms against the activity of TEs might facilitate adaptive responses to environmental challenges and how such mechanisms might be particularly relevant in species frequently facing novel environments, such as invasive, pathogenic or parasitic species.

Accessories Make the Outfit: Accessory Chromosomes and Other Dispensable DNA Regions in Plant-Pathogenic Fungi

Fungal pathogen genomes can often be divided into core and accessory regions. Accessory regions ARs) may be comprised of either ARs (within core chromosomes (CCs) or wholly dispensable (accessory) chromosomes (ACs). Fungal ACs and ARs typically accumulate mutations and structural rearrangements more rapidly over time than CCs and many harbor genes relevant to host-pathogen interactions. These regions are of particular interest in plant pathology and include host-specific virulence factors and secondary metabolite synthesis gene clusters. This review outlines known ACs and ARs in fungal genomes, methods used for their detection, their common properties that differentiate them from the core genome, and what is currently known of their various roles in pathogenicity. Reports on the evolutionary processes generating and shaping AC and AR compartments are discussed, including repeat induced point mutation and breakage fusion bridge cycles. Previously ACs have been studied extensively within key genera, including Fusarium, Zymoseptoria, and Alternaria, but are growing in frequency of observation and perceived importance across a wider range of fungal species. Recent advances in sequencing technologies permit affordable genome assembly and resequencing of populations that will facilitate further discovery and routine screening of ACs.

A fungal avirulence factor encoded in a highly plastic genomic region triggers partial resistance to septoria tritici blotch

Cultivar-strain specificity in the wheat-Zymoseptoria tritici pathosystem determines the infection outcome and is controlled by resistance genes on the host side, many of which have been identified. On the pathogen side, however, the molecular determinants of specificity remain largely unknown. We used genetic mapping, targeted gene disruption and allele swapping to characterise the recognition of the new avirulence factor Avr3D1. We then combined population genetic and comparative genomic analyses to characterise the evolutionary trajectory of Avr3D1. Avr3D1 is specifically recognised by wheat cultivars harbouring the Stb7 resistance gene, triggering a strong defence response without preventing pathogen infection and reproduction. Avr3D1 resides in a cluster of putative effector genes located in a genome region populated by independent transposable element insertions. The gene was present in all 132 investigated strains and is highly polymorphic, with 30 different protein variants identified. We demonstrated that specific amino acid substitutions in Avr3D1 led to evasion of recognition. These results demonstrate that quantitative resistance and gene-for-gene interactions are not mutually exclusive. Localising avirulence genes in highly plastic genomic regions probably facilitates accelerated evolution that enables escape from recognition by resistance proteins.

The coming of age of EvoMPMI: evolutionary molecular plant-microbe interactions across multiple timescales

Plant-microbe interactions are great model systems to study co-evolutionary dynamics across multiple timescales. However, mechanistic research on plant-microbe interactions has often been conducted with little consideration of evolutionary concepts and methods. Conversely, evolutionary research has rarely integrated the range of mechanisms and models from the molecular plant-microbe interactions field. In recent years, the incipient field of evolutionary molecular plant-microbe interactions (EvoMPMI) has emerged to bridge this gap. Here, we report on some of the recent advances in EvoMPMI. In particular, we highlight new systems to study microbe interactions with early diverging land plants, and new findings from studies of adaptive evolution in pathogens and plants. By linking mechanistic and evolutionary research, EvoMPMI promises to expand our understanding of plant-microbe interactions.

Up-regulated RxLR effector genes of Plasmopara viticola in synchronized host-free stages and infected leaves of hosts with different susceptibility

Fast recognition of host signals and early activation of infection mechanisms in Plasmopara viticola are decisive for successful infestation of Vitis vinifera. To better understand interactive processes at the first front line of combat between the pathogen and its host, a specific pre-infective stage was generated in a host-free system. Zoospore encystment was triggered within minutes after treatment with CaCl₂. Subsequently, high rates of germ tube formation occurred in a synchronized manner. This method was employed to compare development-related gene expression in strains of different virulence. Soon after germination, spores showed strong up-regulation of two effector genes, PvRxLR18 and PvRxLR28, particularly in the high virulence strain. On infected grapevine leaf-discs of cultivars with different susceptibility, a similar up-regulation was found at 6 hours post inoculation (hpi). This effect was much more evident in the high virulence than in the low virulence strain and was significantly higher on leaves of the tolerant cultivar Regent than on Müller-Thurgau. In addition, PvRxLR67 was up-regulated 24 hpi in the high virulence strain indicating that different effectors are active in later infection stages. Differences in the expression pattern of RxLR effector genes between the two strains corroborated with infection symptoms visible by sporulation.

Transposable element insertions shape gene regulation and melanin production in a fungal pathogen of wheat

BACKGROUND

Fungal plant pathogens pose major threats to crop yield and sustainable food production if they are highly adapted to their host and the local environment. Variation in gene expression contributes to phenotypic diversity within fungal species and affects adaptation. However, very few cases of adaptive regulatory changes have been reported in fungi and the underlying mechanisms remain largely unexplored. Fungal pathogen genomes are highly plastic and harbor numerous insertions of transposable elements, which can potentially contribute to gene expression regulation. In this work, we elucidated how transposable elements contribute to variation in melanin accumulation, a quantitative trait in fungi that affects survival under stressful conditions.

RESULTS

We demonstrated that differential transcriptional regulation of the gene encoding the transcription factor Zmr1, which controls expression of the genes in the melanin biosynthetic gene cluster, is responsible for variation in melanin accumulation in the fungal plant pathogen Zymoseptoria tritici. We show that differences in melanin levels between two strains of Z. tritici are due to two levels of transcriptional regulation: (1) variation in the promoter sequence of Zmr1 and (2) an insertion of transposable elements upstream of the Zmr1 promoter. Remarkably, independent insertions of transposable elements upstream of Zmr1 occurred in 9% of Z. tritici strains from around the world and negatively regulated Zmr1 expression, contributing to variation in melanin accumulation.

CONCLUSIONS

Our studies identified two levels of transcriptional control that regulate the synthesis of melanin. We propose that these regulatory mechanisms evolved to balance the fitness costs associated with melanin production against its positive contribution to survival in stressful environments.

The genome assembly of the fungal pathogen Pyrenochaeta lycopersici from Single-Molecule Real-Time sequencing sheds new light on its biological complexity

The first draft genome sequencing of the non-model fungal pathogen Pyrenochaeta lycopersici showed an expansion of gene families associated with heterokaryon incompatibility and lacking of mating-type genes, providing insights into the genetic basis of this “imperfect” fungus which lost the ability to produce the sexual stage. However, due to the Illumina short-read technology, the draft genome was too fragmented to allow a comprehensive characterization of the genome, especially of the repetitive sequence fraction. In this work, the sequencing of another P. lycopersici isolate using long-read Single Molecule Real-Time sequencing technology was performed with the aim of obtaining a gapless genome. Indeed, a gapless genome assembly of 62.7 Mb was obtained, with a fraction of repetitive sequences representing 30% of the total bases. The gene content of the two P. lycopersici isolates was very similar, and the large difference in genome size (about 8 Mb) might be attributable to the high fraction of repetitive sequences detected for the new sequenced isolate. The role of repetitive elements, including transposable elements, in modulating virulence effectors is well established in fungal plant pathogens. Moreover, transposable elements are of fundamental importance in creating and re-modelling genes, especially in imperfect fungi. Their abundance in P. lycopersici, together with the large expansion of heterokaryon incompatibility genes in both sequenced isolates, suggest the presence of possible mechanisms alternative to gene re-assorting mediated by sexual recombination. A quite large fraction (~9%) of repetitive elements in P. lycopersici, has no homology with known classes, strengthening this hypothesis. The availability of a gapless genome of P. lycopersici allowed the in-depth analysis of its genome content, by annotating functional genes and TEs. This goal will be an important resource for shedding light on the evolution of the reproductive and pathogenic behaviour of this soilborne pathogen and the onset of a possible speciation within this species.

Nine draft genome sequences of Claviceps purpurea s.lat., including C. arundinis, C. humidiphila, and C. cf. spartinae, pseudomolecules for the pitch canker pathogen Fusarium circinatum, draft genome of Davidsoniella eucalypti, Grosmannia galeiformis, Quambalaria eucalypti, and Teratosphaeria destructans

This genome announcement includes draft genomes from Claviceps purpurea s.lat., including C. arundinis, C. humidiphila and C. cf. spartinae. The draft genomes of Davidsoniella eucalypti, Quambalaria eucalypti and Teratosphaeria destructans, all three important eucalyptus pathogens, are presented. The insect associate Grosmannia galeiformis is also described. The pine pathogen genome of Fusarium circinatum has been assembled into pseudomolecules, based on additional sequence data and by harnessing the known synteny within the Fusarium fujikuroi species complex. This new assembly of the F. circinatum genome provides 12 pseudomolecules that correspond to the haploid chromosome number of F. circinatum. These are comparable to other chromosomal assemblies within the FFSC and will enable more robust genomic comparisons within this species complex.

Novel jack-in-the-box effector of the barley powdery mildew pathogen?

This article comments on: Nottensteiner M, Zechmann B, McCollum C, Hückelhoven R. 2018. A barley powdery mildew fungus non-autonomous retrotransposon encodes a peptide that supports penetration success on barley. Journal of Experimental Botany 69, 3745–3758.

Multiple independent origins for a subtelomeric locus associated with growth rate in Fusarium circinatum

Fusarium is a diverse assemblage that includes a large number of species of considerable medical and agricultural importance. Not surprisingly, whole genome sequences for many Fusarium species have been published or are in the process of being determined, the availability of which is invaluable for deciphering the genetic basis of key phenotypic traits. Here we investigated the distribution, genic composition, and evolutionary history of a locus potentially determining growth rate in the pitch canker pathogen F. circinatum. We found that the genomic region underlying this locus is highly conserved amongst F. circinatum and its close relatives, except for the presence of a 12 000 base pair insertion in all of the examined isolates of F. circinatum. This insertion encodes for five genes and our phylogenetic analyses revealed that each was most likely acquired through horizontal gene transfer from polyphyletic origins. Our data further showed that this region is located in a region low in G+C content and enriched for repetitive sequences and transposable elements, which is situated near the telomere of Chromosome 3 of F. circinatum. As have been shown for other fungi, these findings thus suggest that the emergence of the unique 12 000 bp region in F. circinatum is linked to the dynamic evolutionary processes associated with subtelomeres that, in turn, have been implicated in the ecological adaptation of fungal pathogens.

Signatures of host specialization and a recent transposable element burst in the dynamic one-speed genome of the fungal barley powdery mildew pathogen

Background

Powdery mildews are biotrophic pathogenic fungi infecting a number of economically important plants. The grass powdery mildew, Blumeria graminis, has become a model organism to study host specialization of obligate biotrophic fungal pathogens. We resolved the large-scale genomic architecture of B. graminis forma specialis hordei (Bgh) to explore the potential influence of its genome organization on the co-evolutionary process with its host plant, barley (Hordeum vulgare).

Results

The near-chromosome level assemblies of the Bgh reference isolate DH14 and one of the most diversified isolates, RACE1, enabled a comparative analysis of these haploid genomes, which are highly enriched with transposable elements (TEs). We found largely retained genome synteny and gene repertoires, yet detected copy number variation (CNV) of secretion signal peptide-containing protein-coding genes (SPs) and locally disrupted synteny blocks. Genes coding for sequence-related SPs are often locally clustered, but neither the SPs nor the TEs reside preferentially in genomic regions with unique features. Extended comparative analysis with different host-specific B. graminis formae speciales revealed the existence of a core suite of SPs, but also isolate-specific SP sets as well as congruence of SP CNV and phylogenetic relationship. We further detected evidence for a recent, lineage-specific expansion of TEs in the Bgh genome.

Conclusions

The characteristics of the Bgh genome (largely retained synteny, CNV of SP genes, recently proliferated TEs and a lack of significant compartmentalization) are consistent with a “one-speed” genome that differs in its architecture and (co-)evolutionary pattern from the “two-speed” genomes reported for several other filamentous phytopathogens.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4750-6) contains supplementary material, which is available to authorized users.

Development of a Highly Sensitive FcMito qPCR Assay for the Quantification of the Toxigenic Fungal Plant Pathogen Fusarium culmorum

Fusarium culmorum is a ubiquitous, soil-borne fungus (ascomycete) causing foot and root rot and Fusarium head blight on cereals. It is responsible for yield and quality losses as well as grain contamination with mycotoxins, which are a potential health hazard. An extremely sensitive mitochondrial-based qPCR assay (FcMito qPCR) for quantification of F. culmorum was developed in this study. To provide specificity, the FcMito assay was successfully validated against 85 F. culmorum strains and 53 isolates of 30 other fungal species. The assay efficiency and sensitivity were evaluated against different F. culmorum strains with various amounts of pure fungal DNA and in the presence of background wheat DNA. The results demonstrated the high efficiency of the assay (97.2⁻106.0%, R²-values > 0.99). It was also shown that, in the presence of background DNA, 0.01 pg of fungal template could be reliably quantified. The FcMito assay was used to quantify F. culmorum DNA using 108 grain samples with different trichothecene levels. A significant positive correlation was found between fungal DNA quantity and the total trichothecene content. The obtained results showed that the sensitivity of the FcMito assay was much higher than the nuclear-based qPCR assay for F. culmorum.

A Dispensable Chromosome Is Required for Virulence in the Hemibiotrophic Plant Pathogen Colletotrichum higginsianum

The hemibiotrophic plant pathogen Colletotrichum higginsianum infects Brassicaceae and in combination with Arabidopsis thaliana, represents an important model system to investigate various ecologically important fungal pathogens and their infection strategies. After penetration of plant cells by appressoria, C. higginsianum establishes large biotrophic primary hyphae in the first infected cell. Shortly thereafter, a switch to necrotrophic growth occurs leading to the invasion of neighboring cells by secondary hyphae. In a forward genetic screen for virulence mutants by insertional mutagenesis, we identified mutants that penetrate the plant but show a defect in the passage from biotrophy to necrotrophy. Genome sequencing and pulsed-field gel electrophoresis revealed that two mutants were lacking chromosome 11, encoding potential pathogenicity genes. We established a chromosome loss assay to verify that strains lacking this small chromosome abort infection during biotrophy, while their ability to grow on artificial media was not affected. C. higginsianum harbors a second small chromosome, which can be lost without effects on virulence or growth on agar plates. Furthermore, we found that chromosome 11 is required to suppress Arabidopsis thaliana plant defense mechanisms dependent on tryptophan derived secondary metabolites.

Whole genome and transcriptome analysis reveal adaptive strategies and pathogenesis of Calonectria pseudoreteaudii to Eucalyptus

BACKGROUND

Leaf blight caused by Calonectria spp. is one of the most destructive diseases to affect Eucalyptus nurseries and plantations. These pathogens mainly attack Eucalyptus, a tree with a diversity of secondary metabolites employed as defense-related phytoalexins. To unravel the fungal adaptive mechanisms to various phytoalexins, we examined the genome of C. pseudoreteaudii, which is one of the most aggressive pathogens in southeast Asia.

RESULTS

A 63.7 Mb genome with 14,355 coding genes of C. pseudoreteaudii were assembled. Genomic comparisons identified 1785 species-specific gene families in C. pseudoreteaudii. Most of them were not annotated and those annotated genes were enriched in peptidase activity, pathogenesis, oxidoreductase activity, etc. RNA-seq showed that 4425 genes were differentially expressed on the eucalyptus(the resistant cultivar E. grandis×E.camaldulensis M1) tissue induced medium. The annotation of GO term and KEGG pathway indicated that some of the differential expression genes were involved in detoxification and transportation, such as genes encoding ABC transporters, degrading enzymes of aromatic compounds and so on.

CONCLUSIONS

Potential genomic determinants of phytoalexin detoxification were identified in C. pseudoreteaudii by comparison with 13 other fungi. This pathogen seems to employ membrane transporters and degradation enzymes to detoxify Eucalyptus phytoalexins. Remarkably, the Calonectria genome possesses a surprising number of secondary metabolism backbone enzyme genes involving toxin biosynthesis. It is also especially suited for cutin and lignin degradation. This indicates that toxin and cell wall degrading enzymes may act important roles in the establishment of Calonectria leaf blight. This study provides further understanding on the mechanism of pathogenesis in Calonectria.

Massive Expansion of Gypsy-Like Retrotransposons in Microbotryum Fungi

Transposable elements (TEs) are selfish, autonomously replicating DNA sequences that constitute a major component of eukaryotic genomes and contribute to genome evolution through their movement and amplification. Many fungal genomes, including the anther-smut fungi in the basidiomycete genus Microbotryum, have genome defense mechanisms, such as repeat-induced point mutation (RIP), which hypermutate repetitive DNA and limit TE activity. Little is known about how hypermutation affects the tempo of TE activity and their sequence evolution. Here we report the identification of a massive burst-like expansion of Gypsy-like retrotransposons in a strain of Microbotryum. This TE expansion evidently occurred in the face of RIP-like hypermutation activity. By examining the fitness of individual TE insertion variants, we found that RIP-like mutations impair TE fitness and limit proliferation. Our results provide evidence for a punctuated pattern of TE expansion in a fungal genome, similar to that observed in animals and plants. While targeted hypermutation is often thought of as an effective protection against mobile element activity, our findings suggest that active TEs can persist and undergo selection while they proliferate in genomes that have RIP-like defenses.

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

The first eukaryotic genome to be sequenced was fungal, and there continue to be more sequenced genomes in the kingdom Fungi than in any other eukaryotic kingdom. Comparison of these genomes reveals many sources of genetic variation, from single nucleotide polymorphisms to horizontal gene transfer and on to changes in the arrangement and number of chromosomes, not to mention endofungal bacteria and viruses. Population genomics shows that all sources generate variation all the time and implicate natural selection as the force maintaining genome stability. Variation in wild populations is a rich resource for associating genetic variation with phenotypic variation, whether through quantitative trait locus mapping, genome-wide association studies, or reverse ecology. Subjects of studies associating genetic and phenotypic variation include model fungi, e.g., Saccharomyces and Neurospora, but pioneering studies have also been made with fungi pathogenic to plants, e.g., Pyricularia (= Magnaporthe), Zymoseptoria, and Fusarium, and to humans, e.g., Coccidioides, Cryptococcus, and Candida.

The G-protein γ subunit of Phytophthora infestans is involved in sporangial development

The oomycete Phytophthora infestans is a notorious plant pathogen with potato and tomato as its primary hosts. Previous research showed that the heterotrimeric G-protein subunits Gα and Gβ have a role in zoospore motility and virulence, and sporangial development, respectively. Here, we present analyses of the gene encoding a Gγ subunit in P. infestans, Pigpg1. The overall similarity of PiGPG1 with non-oomycete Gγ subunits is low, with only the most conserved amino acids maintained, but similarity with its homologs in other oomycetes is high. Pigpg1 is expressed in all life stages and shows a similar expression profile as the gene encoding the Gβ subunit, Pigpb1. To elucidate its function, transformants were generated in which Pigpg1 is silenced or overexpressed and their phenotypes were analyzed. Pigpg1-silenced lines produce less sporangia, which are malformed. Altogether, the results show that PiGPG1 is crucial for proper sporangia development and zoosporogenesis. PiGPG1 is a functional Gγ, and likely forms a dimer with PiGPB1 that mediates signaling.

Heterogeneous rates of genome rearrangement contributed to the disparity of species richness in Ascomycota

Background

Chromosomal rearrangements have been shown to facilitate speciation through creating a barrier of gene flow. However, it is not known whether heterogeneous rates of chromosomal rearrangement at the genome scale contributed to the huge disparity of species richness among different groups of organisms, which is one of the most remarkable and pervasive patterns on Earth. The largest fungal phylum Ascomycota is an ideal study system to address this question because it comprises three subphyla (Saccharomycotina, Taphrinomycotina, and Pezizomycotina) whose species numbers differ by two orders of magnitude (59,000, 1000, and 150 respectively).

Results

We quantified rates of genome rearrangement for 71 Ascomycota species that have well-assembled genomes. The rates of inter-species genome rearrangement, which were inferred based on the divergence rates of gene order, are positively correlated with species richness at both ranks of subphylum and class in Ascomycota. This finding is further supported by our quantification of intra-species rearrangement rates based on paired-end genome sequencing data of 216 strains from three representative species, suggesting a difference of intrinsic genome instability among Ascomycota lineages. Our data also show that different rates of imbalanced rearrangements, such as deletions, are a major contributor to the heterogenous rearrangement rates.

Conclusions

Various lines of evidence in this study support that a higher rate of rearrangement at the genome scale might have accelerated the speciation process and increased species richness during the evolution of Ascomycota species. Our findings provide a plausible explanation for the species disparity among Ascomycota lineages, which will be valuable to unravel the underlying causes for the huge disparity of species richness in various taxonomic groups.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4683-0) contains supplementary material, which is available to authorized users.

The ash dieback invasion of Europe was founded by two genetically divergent individuals

Accelerating international trade and climate change make pathogen spread an increasing concern. Hymenoscyphus fraxineus, the causal agent of ash dieback, is a fungal pathogen that has been moving across continents and hosts from Asian to European ash. Most European common ash trees (Fraxinus excelsior) are highly susceptible to H. fraxineus, although a minority (~5%) have partial resistance to dieback. Here, we assemble and annotate a H. fraxineus draft genome which approaches chromosome scale. Pathogen genetic diversity across Europe and in Japan, reveals a strong bottleneck in Europe, though a signal of adaptive diversity remains in key host interaction genes. We find that the European population was founded by two divergent haploid individuals. Divergence between these haplotypes represents the ancestral polymorphism within a large source population. Subsequent introduction from this source would greatly increase adaptive potential of the pathogen. Thus, further introgression of H. fraxineus into Europe represents a potential threat and Europe-wide biological security measures are needed to manage this disease.

Comparative genomics of the wheat fungal pathogen Pyrenophora tritici-repentis reveals chromosomal variations and genome plasticity

Background

Pyrenophora tritici-repentis (Ptr) is a necrotrophic fungal pathogen that causes the major wheat disease, tan spot. We set out to provide essential genomics-based resources in order to better understand the pathogenicity mechanisms of this important pathogen.

Results

Here, we present eight new Ptr isolate genomes, assembled and annotated; representing races 1, 2 and 5, and a new race. We report a high quality Ptr reference genome, sequenced by PacBio technology with Illumina paired-end data support and optical mapping. An estimated 98% of the genome coverage was mapped to 10 chromosomal groups, using a two-enzyme hybrid approach. The final reference genome was 40.9 Mb and contained a total of 13,797 annotated genes, supported by transcriptomic and proteogenomics data sets.

Conclusions

Whole genome comparative analysis revealed major chromosomal segmental rearrangements and fusions, highlighting intraspecific genome plasticity in this species. Furthermore, the Ptr race classification was not supported at the whole genome level, as phylogenetic analysis did not cluster the ToxA producing isolates. This expansion of available Ptr genomics resources will directly facilitate research aimed at controlling tan spot disease.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4680-3) contains supplementary material, which is available to authorized users.

Contribution of transposable elements in the plant's genome

Plants maintain extensive growth flexibility under different environmental conditions, allowing them to continuously and rapidly adapt to alterations in their environment. A large portion of many plant genomes consists of transposable elements (TEs) that create new genetic variations within plant species. Different types of mutations may be created by TEs in plants. Many TEs can avoid the host's defense mechanisms and survive alterations in transposition activity, internal sequence and target site. Thus, plant genomes are expected to utilize a variety of mechanisms to tolerate TEs that are near or within genes. TEs affect the expression of not only nearby genes but also unlinked inserted genes. TEs can create new promoters, leading to novel expression patterns or alternative coding regions to generate alternate transcripts in plant species. TEs can also provide novel cis-acting regulatory elements that act as enhancers or inserts within original enhancers that are required for transcription. Thus, the regulation of plant gene expression is strongly managed by the insertion of TEs into nearby genes. TEs can also lead to chromatin modifications and thereby affect gene expression in plants. TEs are able to generate new genes and modify existing gene structures by duplicating, mobilizing and recombining gene fragments. They can also facilitate cellular functions by sharing their transposase-coding regions. Hence, TE insertions can not only act as simple mutagens but can also alter the elementary functions of the plant genome. Here, we review recent discoveries concerning the contribution of TEs to gene expression in plant genomes and discuss the different mechanisms by which TEs can affect plant gene expression and reduce host defense mechanisms.

Comparative genomics provides insights into the lifestyle and reveals functional heterogeneity of dark septate endophytic fungi

Dark septate endophytes (DSE) are a form-group of root endophytic fungi with elusive functions. Here, the genomes of two common DSE of semiarid areas, Cadophora sp. and Periconia macrospinosa were sequenced and analyzed with another 32 ascomycetes of different lifestyles. Cadophora sp. (Helotiales) and P. macrospinosa (Pleosporales) have genomes of 70.46 Mb and 54.99 Mb with 22,766 and 18,750 gene models, respectively. The majority of DSE-specific protein clusters lack functional annotation with no similarity to characterized proteins, implying that they have evolved unique genetic innovations. Both DSE possess an expanded number of carbohydrate active enzymes (CAZymes), including plant cell wall degrading enzymes (PCWDEs). Those were similar in three other DSE, and contributed a signal for the separation of root endophytes in principal component analyses of CAZymes, indicating shared genomic traits of DSE fungi. Number of secreted proteases and lipases, aquaporins, and genes linked to melanin synthesis were also relatively high in our fungi. In spite of certain similarities between our two DSE, we observed low levels of convergence in their gene family evolution. This suggests that, despite originating from the same habitat, these two fungi evolved along different evolutionary trajectories and display considerable functional differences within the endophytic lifestyle.

Transposable Element Genomic Fissuring in Pyrenophora teres Is Associated With Genome Expansion and Dynamics of Host-Pathogen Genetic Interactions

Pyrenophora teres, P. teres f. teres (PTT) and P. teres f. maculata (PTM) cause significant diseases in barley, but little is known about the large-scale genomic differences that may distinguish the two forms. Comprehensive genome assemblies were constructed from long DNA reads, optical and genetic maps. As repeat masking in fungal genomes influences the final gene annotations, an accurate and reproducible pipeline was developed to ensure comparability between isolates. The genomes of the two forms are highly collinear, each composed of 12 chromosomes. Genome evolution in P. teres is characterized by genome fissuring through the insertion and expansion of transposable elements (TEs), a process that isolates blocks of genic sequence. The phenomenon is particularly pronounced in PTT, which has a larger, more repetitive genome than PTM and more recent transposon activity measured by the frequency and size of genome fissures. PTT has a longer cultivated host association and, notably, a greater range of host-pathogen genetic interactions compared to other Pyrenophora spp., a property which associates better with genome size than pathogen lifestyle. The two forms possess similar complements of TE families with Tc1/Mariner and LINE-like Tad-1 elements more abundant in PTT. Tad-1 was only detectable as vestigial fragments in PTM and, within the forms, differences in genome sizes and the presence and absence of several TE families indicated recent lineage invasions. Gene differences between P. teres forms are mainly associated with gene-sparse regions near or within TE-rich regions, with many genes possessing characteristics of fungal effectors. Instances of gene interruption by transposons resulting in pseudogenization were detected in PTT. In addition, both forms have a large complement of secondary metabolite gene clusters indicating significant capacity to produce an array of different molecules. This study provides genomic resources for functional genetics to help dissect factors underlying the host-pathogen interactions.

Genomics analysis of Aphanomyces spp. identifies a new class of oomycete effector associated with host adaptation

Background

Oomycetes are a group of filamentous eukaryotic microorganisms that have colonized all terrestrial and oceanic ecosystems, and they include prominent plant pathogens. The Aphanomyces genus is unique in its ability to infect both plant and animal species, and as such exemplifies oomycete versatility in adapting to different hosts and environments. Dissecting the underpinnings of oomycete diversity provides insights into their specificity and pathogenic mechanisms.

Results

By carrying out genomic analyses of the plant pathogen A. euteiches and the crustacean pathogen A. astaci, we show that host specialization is correlated with specialized secretomes that are adapted to the deconstruction of the plant cell wall in A. euteiches and protein degradation in A. astaci. The A. euteiches genome is characterized by a large repertoire of small secreted protein (SSP)-encoding genes that are highly induced during plant infection, and are not detected in other oomycetes. Functional analysis revealed an SSP from A. euteiches containing a predicted nuclear-localization signal which shuttles to the plant nucleus and increases plant susceptibility to infection.

Conclusion

Collectively, our results show that Aphanomyces host adaptation is associated with evolution of specialized secretomes and identify SSPs as a new class of putative oomycete effectors.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0508-5) contains supplementary material, which is available to authorized users.