Sex or no sex: evolutionary adaptation occurs regardless

Functional analysis of the mating type genes in Verticillium dahliae

BACKGROUND

Populations of the plant pathogenic fungus Verticillium dahliae display a complex and rich genetic diversity, yet the existence of sexual reproduction in the fungus remains contested. As pivotal genes, MAT genes play a crucial role in regulating cell differentiation, morphological development, and mating of compatible cells. However, the functions of the two mating type genes in V. dahliae, VdMAT1-1-1, and VdMAT1-2-1, remain poorly understood.

RESULTS

In this study, we confirmed that the MAT loci in V. dahliae are highly conserved, including both VdMAT1-1-1 and VdMAT1-2-1 which share high collinearity. The conserved core transcription factor encoded by the two MAT loci may facilitate the regulation of pheromone precursor and pheromone receptor genes by directly binding to their promoter regions. Additionally, peptide activity assays demonstrated that the signal peptide of the pheromone VdPpg1 possessed secretory activity, while VdPpg2, lacked a predicted signal peptide. Chemotactic growth assays revealed that V. dahliae senses and grows towards the pheromones FO-a and FO-α of Fusarium oxysporum, as well as towards VdPpg2 of V. dahliae, but not in response to VdPpg1. The findings herein also revealed that VdMAT1-1-1 and VdMAT1-2-1 regulate vegetative growth, carbon source utilization, and resistance to stressors in V. dahliae, while negatively regulating virulence.

CONCLUSIONS

These findings underscore the potential roles of VdMAT1-1-1 and VdMAT1-2-1 in sexual reproduction and confirm their involvement in various asexual processes of V. dahliae, offering novel insights into the functions of mating type genes in this species.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BACKGROUND

Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

RESULTS

Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.

CONCLUSIONS

Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

Polyethyleneimine-coated MXene quantum dots improve cotton tolerance to Verticillium dahliae by maintaining ROS homeostasis

Verticillium dahliae is a soil-borne hemibiotrophic fungal pathogen that threatens cotton production worldwide. In this study, we assemble the genomes of two V. dahliae isolates: the more virulence and defoliating isolate V991 and nondefoliating isolate 1cd3-2. Transcriptome and comparative genomics analyses show that genes associated with pathogen virulence are mostly induced at the late stage of infection (Stage II), accompanied by a burst of reactive oxygen species (ROS), with upregulation of more genes involved in defense response in cotton. We identify the V991-specific virulence gene SP3 that is highly expressed during the infection Stage II. V. dahliae SP3 knock-out strain shows attenuated virulence and triggers less ROS production in cotton plants. To control the disease, we employ polyethyleneimine-coated MXene quantum dots (PEI-MQDs) that possess the ability to remove ROS. Cotton seedlings treated with PEI-MQDs are capable of maintaining ROS homeostasis with enhanced peroxidase, catalase, and glutathione peroxidase activities and exhibit improved tolerance to V. dahliae. These results suggest that V. dahliae trigger ROS production to promote infection and scavenging ROS is an effective way to manage this disease. This study reveals a virulence mechanism of V. dahliae and provides a means for V. dahliae resistance that benefits cotton production.

High Nucleotide Substitution Rates Associated with Retrotransposon Proliferation Drive Dynamic Secretome Evolution in Smut Pathogens

Transposable elements (TEs) play a pivotal role in shaping diversity in eukaryotic genomes. The covered smut pathogen on barley, Ustilago hordei, encountered a recent genome expansion. Using long reads, we assembled genomes of 6 U. hordei strains and 3 sister species, to study this genome expansion. We found that larger genome sizes can mainly be attributed to a higher genome fraction of long terminal repeat retrotransposons (LTR-RTs). In the studied smut genomes, LTR-RTs fractions are the largest in U. hordei and are positively correlated with the mating-type locus sizes, which is up to ~560 kb in U. hordei. Furthermore, LTR-RTs were found to be associated with higher nucleotide substitution levels, as these occur in specific genome regions of smut species with a recent LTR-RT proliferation. Moreover, genes in genome regions with higher nucleotide substitution levels generally reside closer to LTR-RTs than other genome regions. Genome regions with many nucleotide substitutions encountered an especially high fraction of CG substitutions, which is not observed for LTR-RT sequences. The high nucleotide substitution levels particularly accelerate the evolution of secretome genes, as their more accessory nature results in substitutions that often lead to amino acid alterations. IMPORTANCE Genomic alteration can be generated through various means, in which transposable elements (TEs) can play a pivotal role. Their mobility causes mutagenesis in itself and can disrupt the function of the sequences they insert into. They also impact genome evolution as their repetitive nature facilitates nonhomologous recombination. Furthermore, TEs have been linked to specific epigenetic genome organizations. We report a recent TE proliferation in the genome of the barley covered smut fungus, Ustilago hordei. This proliferation is associated with a distinct nucleotide substitution regime that has a higher rate and a higher fraction of CG substitutions. This different regime shapes the evolution of genes in subjected genome regions. We hypothesize that TEs may influence the error-rate of DNA polymerase in a hitherto unknown fashion.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga’s genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host–pathogen interactions.

Sexual reproduction contributes to the evolution of resistance breaking isolates of the spinach pathogen Peronospora effusa

Peronospora effusa causes downy mildew, the economically most important disease of cultivated spinach worldwide. To date, 19 P. effusa races have been denominated based on their capacity to break spinach resistances, but their genetic diversity and the evolutionary processes that contribute to race emergence are unknown. Here, we performed the first systematic analysis of P. effusa races showing that those emerge by both asexual and sexual reproduction. Specifically, we studied the diversity of 26 P. effusa isolates from 16 denominated races based on mitochondrial and nuclear comparative genomics. Mitochondrial genomes based on long‐read sequencing coupled with diversity assessment based on short‐read sequencing uncovered two mitochondrial haplogroups, each with distinct genome organization. Nuclear genome‐wide comparisons of the 26 isolates revealed that 10 isolates from six races could clearly be divided into three asexually evolving groups, in concordance with their mitochondrial phylogeny. The remaining isolates showed signals of reticulated evolution and discordance between nuclear and mitochondrial phylogenies, suggesting that these evolved through sexual reproduction. Increased understanding of this pathogen's reproductive modes will provide the framework for future studies into the molecular mechanisms underlying race emergence and into the P. effusa‐spinach interaction, thus assisting in sustainable production of spinach through knowledge‐driven resistance breeding.

Unraveling the Genome Sequence of Plant Growth Promoting Aspergillus niger (CSR3) Provides Insight into the Synthesis of Secondary Metabolites and Its Comparative Genomics

Aspergillus niger strain CSR3 is an endophytic fungus that regulates plant endogenous hormones, secondary metabolites, and promotes plant growth during abiotic stress conditions. In this study, we sequenced the genome of A. niger (CSR3) and compared it with previously available A. niger strains. The final genome assembly was 35.8 Mb in size, consisting of 23 scaffolds with N50 scaffold length of 2.4 Mb. A total of 12,442 protein coding genes, 270 tRNA, and 57 rRNA were predicted in the CSR3 genome. We used comparative genomic analysis to provide insights into the genome’s evolution and to elucidate the adaptive genomic signatures for bioactive secondary metabolite biosynthesis, hormones biosynthesis, and plant growth promoting activities. We also analyzed the transposable elements (TEs), simple sequence repeats (SSRs), CAZymes families, genes involved in gibberellin biosynthesis, and secondary metabolite clusters in the CSR3 genome. A total of 21 secondary metabolite biosynthesis gene clusters were detected, with 18 essential enzymes involved in the mevalonate pathway (MVA). The repeat analysis revealed about 3431 SSR, 274 TEs, and 205 inverted repeats (IR). Further gene family analysis revealed that 124 gene families were gained, whereas 125 gene families were lost in CSR3 genome, compared to A. niger ASM151534V and A. niger ASM285V2 genomes. The results improve our understanding of the CSR3 genome and will assist in future investigations on the genetic basis of A. niger CSR3, including the identification of CSR3 phytostimulant properties.

Chromatin profiling reveals heterogeneity in clinical isolates of the human pathogen Aspergillus fumigatus

Invasive Pulmonary Aspergillosis, which is caused by the filamentous fungus Aspergillus fumigatus, is a life-threatening infection for immunosuppressed patients. Chromatin structure regulation is important for genome stability maintenance and has the potential to drive genome rearrangements and affect virulence and pathogenesis of pathogens. Here, we performed the first A. fumigatus global chromatin profiling of two histone modifications, H3K4me3 and H3K9me3, focusing on the two most investigated A. fumigatus clinical isolates, Af293 and CEA17. In eukaryotes, H3K4me3 is associated with active transcription, while H3K9me3 often marks silent genes, DNA repeats, and transposons. We found that H3K4me3 deposition is similar between the two isolates, while H3K9me3 is more variable and does not always represent transcriptional silencing. Our work uncovered striking differences in the number, locations, and expression of transposable elements between Af293 and CEA17, and the differences are correlated with H3K9me3 modifications and higher genomic variations among strains of Af293 background. Moreover, we further showed that the Af293 strains from different laboratories actually differ in their genome contents and found a frequently lost region in chromosome VIII. For one such Af293 variant, we identified the chromosomal changes and demonstrated their impacts on its secondary metabolites production, growth and virulence. Overall, our findings not only emphasize the influence of genome heterogeneity on A. fumigatus fitness, but also caution about unnoticed chromosomal variations among common laboratory strains.

Phenotypic Characterization and Comparative Genomics of the Melanin-Producing Yeast Exophiala lecanii-corni Reveals a Distinct Stress Tolerance Profile and Reduced Ribosomal Genetic Content

The black yeast Exophiala lecanii-corni of the order Chaetothyriales is notable for its ability to produce abundant quantities of DHN-melanin. While many other Exophiala species are frequent causal agents of human infection, E. lecanii-corni CBS 102400 lacks the thermotolerance requirements that enable pathogenicity, making it appealing for use in targeted functional studies and biotechnological applications. Here, we report the stress tolerance characteristics of E. lecanii-corni, with an emphasis on the influence of melanin on its resistance to various forms of stress. We find that E. lecanii-corni has a distinct stress tolerance profile that includes variation in resistance to temperature, osmotic, and oxidative stress relative to the extremophilic and pathogenic black yeast Exophiala dermatitidis. Notably, the presence of melanin substantially impacts stress resistance in E. lecanii-corni, while this was not found to be the case in E. dermatitidis. The cellular context, therefore, influences the role of melanin in stress protection. In addition, we present a detailed analysis of the E. lecanii-corni genome, revealing key differences in functional genetic content relative to other ascomycetous species, including a significant decrease in abundance of genes encoding ribosomal proteins. In all, this study provides insight into how genetics and physiology may underlie stress tolerance and enhances understanding of the genetic diversity of black yeasts.

Genome Sequence Data of MAT1-1 and MAT1-2 Idiomorphs from Verticillium dahliae

Though Verticillium dahliae is an asexually reproducing fungus, it is considered heterothallic owing to the presence of only one of the two mating-type idiomorphs (MAT1-1 or MAT1-2) in individual isolates. But sexual reproduction has never been observed either in nature or in the laboratory. All of the genomic information in the literature thus far has therefore come from studies on isolates carrying only the MAT1-2 idiomorph. Herein, we sequenced and compared high-quality reference genomes of MAT1-1 strain S011 and MAT1-2 strain S023 obtained from the same sunflower field. The two genomic sequences displayed high synteny, and encoded similar number genes, a similarity especially notable among pathogenicity-related genes. Homolog analysis between these two genomes revealed that 80% of encoded genes are highly conserved (95% identity and coverage), but only 20% of the single copy genes were identical. These novel genome resources will support the analysis of the structure and function of the two idiomorphs and provide valuable tools to elucidate the evolution and potential mechanisms of sexual reproduction in V. dahliae.

Key Insights and Research Prospects at the Dawn of the Population Genomics Era for Verticillium dahliae

The genomics era has ushered in exciting possibilities to examine the genetic bases that undergird the characteristic features of Verticillium dahliae and other plant pathogens. In this review, we provide historical perspectives on some of the salient biological characteristics of V. dahliae, including its morphology, microsclerotia formation, host range, disease symptoms, vascular niche, reproduction, and population structure. The kaleidoscopic population structure of this pathogen is summarized, including different races of the pathogen, defoliating and nondefoliating phenotypes, vegetative compatibility groupings, and clonal populations. Where possible, we place the characteristic differences in the context of comparative and functional genomics analyses that have offered insights into population divergence within V. dahliae and the related species.Current challenges are highlighted along with some suggested future population genomics studies that will contribute to advancing our understanding of the population divergence in V. dahliae.

A 20-kb lineage-specific genomic region tames virulence in pathogenic amphidiploid Verticillium longisporum

Amphidiploid fungal Verticillium longisporum strains Vl43 and Vl32 colonize the plant host Brassica napus but differ in their ability to cause disease symptoms. These strains represent two V. longisporum lineages derived from different hybridization events of haploid parental Verticillium strains. Vl32 and Vl43 carry same-sex mating-type genes derived from both parental lineages. Vl32 and Vl43 similarly colonize and penetrate plant roots, but asymptomatic Vl32 proliferation in planta is lower than virulent Vl43. The highly conserved Vl43 and Vl32 genomes include less than 1% unique genes, and the karyotypes of 15 or 16 chromosomes display changed genetic synteny due to substantial genomic reshuffling. A 20 kb Vl43 lineage-specific (LS) region apparently originating from the Verticillium dahliae-related ancestor is specific for symptomatic Vl43 and encodes seven genes, including two putative transcription factors. Either partial or complete deletion of this LS region in Vl43 did not reduce virulence but led to induction of even more severe disease symptoms in rapeseed. This suggests that the LS insertion in the genome of symptomatic V. longisporum Vl43 mediates virulence-reducing functions, limits damage on the host plant, and therefore tames Vl43 from being even more virulent.

The Interspecific Fungal Hybrid Verticillium longisporum Displays Subgenome-Specific Gene Expression

Hybridization is an important evolutionary mechanism that can enable organisms to adapt to environmental challenges. It has previously been shown that the fungal allodiploid species Verticillium longisporum, the causal agent of verticillium stem striping in rapeseed, originated from at least three independent hybridization events between two haploid Verticillium species. To reveal the impact of genome duplication as a consequence of hybridization, we studied the genome and transcriptome dynamics upon two independent V. longisporum hybridization events, represented by the hybrid lineages “A1/D1” and “A1/D3.” We show that V. longisporum genomes are characterized by extensive chromosomal rearrangements, including between parental chromosomal sets. V. longisporum hybrids display signs of evolutionary dynamics that are typically associated with the aftermath of allodiploidization, such as haploidization and more relaxed gene evolution. The expression patterns of the two subgenomes within the two hybrid lineages are more similar than those of the shared A1 parent between the two lineages, showing that the expression patterns of the parental genomes homogenized within a lineage. However, as genes that display differential parental expression in planta do not typically display the same pattern in vitro, we conclude that subgenome-specific responses occur in both lineages. Overall, our study uncovers genomic and transcriptomic plasticity during the evolution of the filamentous fungal hybrid V. longisporum and illustrates its adaptive potential.

Transposable Elements Contribute to Genome Dynamics and Gene Expression Variation in the Fungal Plant Pathogen Verticillium dahliae

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization and adaptation. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as polymorphic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the total TE complement. As expected, the polyrmorphic TEs are enriched in the adaptive genomic regions. Besides, we observed an association of polymorphic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of adaptive genomic regions.

Integrated transcriptome and metabolome analyses of biochar-induced pathways in response to Fusarium wilt infestation in pepper

The present study used soils contaminated with Fusarium oxysporum f. sp. capsici (CCS) and CCS amended with bamboo biochar (CCS + BC) to grow the pepper variety Qujiao No.1. The physiological performance, and transcriptome and metabolome profiling in leaf (L) and fruit (F) of Qujiao No.1 were conducted. Application of biochar improved soil properties, pepper plant nutrition and increased activities of enzymes related to pest/disease resistance, leading to superior physiological performance and lesser F. wilt disease incidence than plants from CCS. Most of the differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were involved in protein processing in endoplasmic reticulum (fruit), plant pathogen interaction (fruit), photosynthesis (leaf), phenylpropanoid biosynthesis (both tissues) and metabolic pathways (both tissues). Biochar improved plant photosynthesis, enhanced the immune system, energy production and increased stress signaling pathways. Overall, our results provide evidence of a number of pathways induced by biochar in pepper regulating its response to F. wilt disease.

Evolutionary and genomic comparisons of hybrid uninucleate and nonhybrid Rhizoctonia fungi

The basidiomycetous fungal genus, Rhizoctonia, can cause severe damage to many plants and is composed of multinucleate, binucleate, and uninucleate species differing in pathogenicity. Here we generated chromosome-scale genome assemblies of the three nuclear types of Rhizoctonia isolates. The genomic comparisons revealed that the uninucleate JN strain likely arose by somatic hybridization of two binucleate isolates, and maintained a diploid nucleus. Homeolog gene pairs in the JN genome have experienced both decelerated or accelerated evolution. Homeolog expression dominance occurred between JN subgenomes, in which differentially expressed genes show potentially less evolutionary constraint than the genes without. Analysis of mating-type genes suggested that Rhizoctonia maintains the ancestral tetrapolarity of the Basidiomycota. Long terminal repeat-retrotransposons displayed a reciprocal correlation with the chromosomal GC content in the three chromosome-scale genomes. The more aggressive multinucleate XN strain had more genes encoding enzymes for host cell wall decomposition. These findings demonstrate some evolutionary changes of a recently derived hybrid and in multiple nuclear types of Rhizoctonia.

Lineage-Specific Genes and Cryptic Sex: Parallels and Differences between Arbuscular Mycorrhizal Fungi and Fungal Pathogens

Arbuscular mycorrhizal fungi (AMF) live as obligate root symbionts on almost all land plants. They have long been regarded as ancient asexuals that have propagated clonally for millions of years. However, genomic studies in Rhizophagus irregularis and other AMF revealed many features indicative of sex. Surprisingly, comparative genomics of conspecific isolates of R. irregularis revealed an unexpected interstrain diversity, suggesting that AMF carry a high number of lineage-specific (LS) genes. Intriguingly, cryptic sex and LS genomic regions have previously been reported in a number of fungal pathogens of plants and humans. Here, we discuss these genomic similarities and highlight their potential relevance for AMF adaptation to the environment and for symbiotic functioning.

A unique chromatin profile defines adaptive genomic regions in a fungal plant pathogen

Genomes store information at scales beyond the linear nucleotide sequence, which impacts genome function at the level of an individual, while influences on populations and long-term genome function remains unclear. Here, we addressed how physical and chemical DNA characteristics influence genome evolution in the plant pathogenic fungus Verticillium dahliae. We identified incomplete DNA methylation of repetitive elements, associated with specific genomic compartments originally defined as Lineage-Specific (LS) regions that contain genes involved in host adaptation. Further chromatin characterization revealed associations with features such as H3 Lys-27 methylated histones (H3K27me3) and accessible DNA. Machine learning trained on chromatin data identified twice as much LS DNA as previously recognized, which was validated through orthogonal analysis, and we propose to refer to this DNA as adaptive genomic regions. Our results provide evidence that specific chromatin profiles define adaptive genomic regions, and highlight how different epigenetic factors contribute to the organization of these regions.

Metarhizium: jack of all trades, master of many

The genus Metarhizium and Pochonia chlamydosporia comprise a monophyletic clade of highly abundant globally distributed fungi that can transition between long-term beneficial associations with plants to transitory pathogenic associations with frequently encountered protozoans, nematodes or insects. Some very common 'specialist generalist' species are adapted to particular soil and plant ecologies, but can overpower a wide spectrum of insects with numerous enzymes and toxins that result from extensive gene duplications made possible by loss of meiosis and associated genome defence mechanisms. These species use parasexuality instead of sex to combine beneficial mutations from separate clonal individuals into one genome (Vicar of Bray dynamics). More weakly endophytic species which kill a narrow range of insects retain sexuality to facilitate host-pathogen coevolution (Red Queen dynamics). Metarhizium species can fit into numerous environments because they are very flexible at the genetic, physiological and ecological levels, providing tractable models to address how new mechanisms for econutritional heterogeneity, host switching and virulence are acquired and relate to diverse sexual life histories and speciation. Many new molecules and functions have been discovered that underpin Metarhizium associations, and have furthered our understanding of the crucial ecology of these fungi in multiple habitats.

Repetitive Elements Contribute to the Diversity and Evolution of Centromeres in the Fungal Genus Verticillium

Centromeres are chromosomal regions that are crucial for chromosome segregation during mitosis and meiosis, and failed centromere formation can contribute to chromosomal anomalies. Despite this conserved function, centromeres differ significantly between and even within species. Thus far, systematic studies into the organization and evolution of fungal centromeres remain scarce. In this study, we identified the centromeres in each of the 10 species of the fungal genus Verticillium and characterized their organization and evolution. Chromatin immunoprecipitation of the centromere-specific histone CenH3 (ChIP-seq) and chromatin conformation capture (Hi-C) followed by high-throughput sequencing identified eight conserved, large (∼150-kb), AT-, and repeat-rich regional centromeres that are embedded in heterochromatin in the plant pathogen Verticillium dahliae Using Hi-C, we similarly identified repeat-rich centromeres in the other Verticillium species. Strikingly, a single degenerated long terminal repeat (LTR) retrotransposon is strongly associated with centromeric regions in some but not all Verticillium species. Extensive chromosomal rearrangements occurred during Verticillium evolution, of which some could be linked to centromeres, suggesting that centromeres contributed to chromosomal evolution. The size and organization of centromeres differ considerably between species, and centromere size was found to correlate with the genome-wide repeat content. Overall, our study highlights the contribution of repetitive elements to the diversity and rapid evolution of centromeres within the fungal genus VerticilliumIMPORTANCE The genus Verticillium contains 10 species of plant-associated fungi, some of which are notorious pathogens. Verticillium species evolved by frequent chromosomal rearrangements that contribute to genome plasticity. Centromeres are instrumental for separation of chromosomes during mitosis and meiosis, and failed centromere functionality can lead to chromosomal anomalies. Here, we used a combination of experimental techniques to identify and characterize centromeres in each of the Verticillium species. Intriguingly, we could strongly associate a single repetitive element to the centromeres of some of the Verticillium species. The presence of this element in the centromeres coincides with increased centromere sizes and genome-wide repeat expansions. Collectively, our findings signify a role of repetitive elements in the function, organization, and rapid evolution of centromeres in a set of closely related fungal species.

Comparative Pathogenicity and Host Ranges of Magnaporthe oryzae and Related Species

Host shifting and host expansion of fungal plant pathogens increases the rate of emergence of new pathogens and the incidence of disease in various crops, which threaten global food security. Magnaporthe species cause serious disease in rice, namely rice blast disease, as well as in many alternative hosts, including wheat, barley, and millet. A severe outbreak of wheat blast due to Magnaporthe oryzae occurred recently in Bangladesh, after the fungus was introduced from South America, causing great loss of yield. This outbreak of wheat blast is of growing concern, because it might spread to adjacent wheat-producing areas. Therefore, it is important to understand the host range and population structure of M. oryzae and related species for determining the evolutionary relationships among Magnaporthe species and for managing blast disease in the field. Here, we collected isolates of M. oryzae and related species from various Poaceae species, including crops and weeds surrounding rice fields, in Korea and determined their phylogenetic relationships and host species specificity. Internal transcribed spacer-mediated phylogenetic analysis revealed that M. oryzae and related species are classified into four groups primarily including isolates from rice, crabgrass, millet and tall fescue. Based on pathogenicity assays, M. oryzae and related species can infect different Poaceae hosts and move among hosts, suggesting the potential for host shifting and host expansion in nature. These results provide important information on the diversification of M. oryzae and related species with a broad range of Poaceae as hosts in crop fields.

Epigenetic Regulation of Verticillium dahliae Virulence: Does DNA Methylation Level Play A Role?

Verticillium dahliae is the etiological agent of Verticillium wilt of olive. The virulence of Defoliating V. dahliae isolates usually displays differences and high plasticity. This work studied whether an epigenetic mechanism was involved in this plasticity. An inverse correlation between virulence and DNA methylation of protein-coding genes was found. A set of 831 genes was selected for their highly consistent inverse methylation profile and virulence in the five studied isolates. Of these genes, ATP-synthesis was highly represented, which indicates that the more virulent D isolates are, the more energy requirements they may have. Furthermore, there were numerous genes in the protein biosynthesis process: genes coding for the chromatin structure, which suggests that epigenetic changes may also affect chromatin condensation; many transmembrane transporter genes, which is consistent with denser compounds, traffic through membranes in more virulent isolates; a fucose-specific lectin that may play a role in the attachment to plant cell walls during the host infection process; and pathogenic cutinases that facilitate plant invasion and sporulation genes for rapid spreading alongside plants. Our findings support the notion that differences in the virulence of the Defoliating V. dahliae isolates may be controlled, at least to some extent, by an epigenetic mechanism.

Asexual Evolution and Forest Conditions Drive Genetic Parallelism in Phytophthora ramorum

It is commonly assumed that asexual lineages are short-lived evolutionarily, yet many asexual organisms can generate genetic and phenotypic variation, providing an avenue for further evolution. Previous work on the asexual plant pathogen Phytophthora ramorum NA1 revealed considerable genetic variation in the form of Structural Variants (SVs). To better understand how SVs arise and their significance to the California NA1 population, we studied the evolutionary histories of SVs and the forest conditions associated with their emergence. Ancestral state reconstruction suggests that SVs arose by somatic mutations among multiple independent lineages, rather than by recombination. We asked if this unusual phenomenon of parallel evolution between isolated populations is transmitted to extant lineages and found that SVs persist longer in a population if their genetic background had a lower mutation load. Genetic parallelism was also found in geographically distant demes where forest conditions such as host density, solar radiation, and temperature, were similar. Parallel SVs overlap with genes involved in pathogenicity such as RXLRs and have the potential to change the course of an epidemic. By combining genomics and environmental data, we identified an unexpected pattern of repeated evolution in an asexual population and identified environmental factors potentially driving this phenomenon.

Emergence of a hybrid PKS-NRPS secondary metabolite cluster in a clonal population of the rice blast fungus Magnaporthe oryzae

Secondary metabolites (SMs) are crucial for fungi and vary in function from beneficial antibiotics to pathogenicity factors. To generate diversified SMs that enable different functions, SM-coding regions rapidly evolve in fungal genomes. However, the driving force and genetic mechanism of fungal SM diversification in the context of host-pathogen interactions remain largely unknown. Previously, we grouped field populations of the rice blast fungus Magnaporthe oryzae (syn: Pyricularia oryzae) into three major globally distributed clades based on population genomic analyses. Here, we characterize a recent duplication of an avirulent gene-containing SM cluster, ACE1, in a clonal M. oryzae population (Clade 2). We demonstrate that the ACE1 cluster is specifically duplicated in Clade 2, a dominant clade in indica rice-growing areas. With long-read sequencing, we obtained chromosome-level genome sequences of four Clade 2 isolates, which displayed differences in genomic organization of the ACE1 duplication process. Comparative genomic analyses suggested that the original ACE1 cluster experienced frequent rearrangement in Clade 2 isolates and revealed that the new ACE1 cluster is located in a newly formed and transposable element-rich region. Taken together, these results highlight the frequent mutation and expansion of an avirulent gene-containing SM cluster through transposable element-mediated whole-cluster duplication in the context of host-pathogen interactions.

Mitotic Recombination and Adaptive Genomic Changes in Human Pathogenic Fungi

Genome rearrangements and ploidy alterations are important for adaptive change in the pathogenic fungal species Candida and Cryptococcus, which propagate primarily through clonal, asexual reproduction. These changes can occur during mitotic growth and lead to enhanced virulence, drug resistance, and persistence in chronic infections. Examples of microevolution during the course of infection were described in both human infections and mouse models. Recent discoveries defining the role of sexual, parasexual, and unisexual cycles in the evolution of these pathogenic fungi further expanded our understanding of the diversity found in and between species. During mitotic growth, damage to DNA in the form of double-strand breaks (DSBs) is repaired, and genome integrity is restored by the homologous recombination and non-homologous end-joining pathways. In addition to faithful repair, these pathways can introduce minor sequence alterations at the break site or lead to more extensive genetic alterations that include loss of heterozygosity, inversions, duplications, deletions, and translocations. In particular, the prevalence of repetitive sequences in fungal genomes provides opportunities for structural rearrangements to be generated by non-allelic (ectopic) recombination. In this review, we describe DSB repair mechanisms and the types of resulting genome alterations that were documented in the model yeast Saccharomyces cerevisiae. The relevance of similar recombination events to stress- and drug-related adaptations and in generating species diversity are discussed for the human fungal pathogens Candida albicans and Cryptococcus neoformans.

Activity and Phylogenetics of the Broadly Occurring Family of Microbial Nep1-Like Proteins

Necrosis- and ethylene-inducing peptide 1 (Nep1)-like proteins (NLP) have an extremely broad taxonomic distribution; they occur in bacteria, fungi, and oomycetes. NLPs come in two forms, those that are cytotoxic to eudicot plants and those that are noncytotoxic. Cytotoxic NLPs bind to glycosyl inositol phosphoryl ceramide (GIPC) sphingolipids that are abundant in the outer leaflet of plant plasma membranes. Binding allows the NLP to become cytolytic in eudicots but not monocots. The function of noncytotoxic NLPs remains enigmatic, but the expansion of NLP genes in oomycete genomes suggests they are important. Several plant species have evolved the capacity to recognize NLPs as molecular patterns and trigger plant immunity, e.g., Arabidopsis thaliana detects nlp peptides via the receptor-like protein RLP23. In this review, we provide a historical perspective from discovery to understanding of molecular mechanisms and describe the latest developments in the NLP field to shed light on these fascinating microbial proteins.

Dynamic virulence-related regions of the plant pathogenic fungus Verticillium dahliae display enhanced sequence conservation

Plant pathogens continuously evolve to evade host immune responses. During host colonization, many fungal pathogens secrete effectors to perturb such responses, but these in turn may become recognized by host immune receptors. To facilitate the evolution of effector repertoires, such as the elimination of recognized effectors, effector genes often reside in genomic regions that display increased plasticity, a phenomenon that is captured in the two-speed genome hypothesis. The genome of the vascular wilt fungus Verticillium dahliae displays regions with extensive presence/absence polymorphisms, so-called lineage-specific regions, that are enriched in in planta-induced putative effector genes. As expected, comparative genomics reveals differential degrees of sequence divergence between lineage-specific regions and the core genome. Unanticipated, lineage-specific regions display markedly higher sequence conservation in coding as well as noncoding regions than the core genome. We provide evidence that disqualifies horizontal transfer to explain the observed sequence conservation and conclude that sequence divergence occurs at a slower pace in lineage-specific regions of the V. dahliae genome. We hypothesize that differences in chromatin organisation may explain lower nucleotide substitution rates in the plastic, lineage-specific regions of V. dahliae.

Link between epigenetic diversity and invasive status of south-eastern European populations of phytopathogenic fungus Cryphonectria parasitica

Epigenetic modifications may play an important role in invasion and adaptation of clonal and invasive populations to different environments. The aim of this study was to analyse epigenetic diversity and structure within and among populations of invasive pathogenic fungus Cryphonectria parasitica from south-eastern Europe, where one haplotype S12 dominates. The highest level of epigenetic diversity was found in haplotype S1, followed by S2, while the lowest level of epigenetic diversity was found in haplotype S12. Similar pattern of epigenetic diversity was detected in the control, genetically diverse Croatian population where S1 haplotype dominates. In four south-eastern European populations, the highest level of epigenetic diversity was observed in the Italian population, the oldest population in the studied area, while the lowest diversity was found in most recently established Bulgarian population. This relationship between epigenetic diversity and population age implies the important role of epigenetic modifications on the process of invasion. Our data suggest that epigenetic differences might affect the success of expansion of certain haplotype into new regions. Understanding the role of epigenetic processes in expansion and (pre)adaptation of fungal plant pathogens, besides fundamental knowledge, can contribute to development of strategies for control of fungal spread and pathogenesis.

Chromatin-dependent regulation of secondary metabolite biosynthesis in fungi: is the picture complete?

Fungal secondary metabolites are small molecules that exhibit diverse biological activities exploited in medicine, industry and agriculture. Their biosynthesis is governed by co-expressed genes that often co-localize in gene clusters. Most of these secondary metabolite gene clusters are inactive under laboratory conditions, which is due to a tight transcriptional regulation. Modifications of chromatin, the complex of DNA and histone proteins influencing DNA accessibility, play an important role in this regulation. However, tinkering with well-characterised chemical and genetic modifications that affect chromatin alters the expression of only few biosynthetic gene clusters, and thus the regulation of the vast majority of biosynthetic pathways remains enigmatic. In the past, attempts to activate silent gene clusters in fungi mainly focused on histone acetylation and methylation, while in other eukaryotes many other post-translational modifications are involved in transcription regulation. Thus, how chromatin regulates the expression of gene clusters remains a largely unexplored research field. In this review, we argue that focusing on only few well-characterised chromatin modifications is significantly hampering our understanding of the chromatin-based regulation of biosynthetic gene clusters. Research on underexplored chromatin modifications and on the interplay between different modifications is timely to fully explore the largely untapped reservoir of fungal secondary metabolites.

The genomic determinants of adaptive evolution in a fungal pathogen

Unravelling the strength, frequency, and distribution of selective variants along the genome as well as the underlying factors shaping this distribution are fundamental goals of evolutionary biology. Antagonistic host-pathogen coevolution is thought to be a major driver of genome evolution between interacting species. While rapid evolution of pathogens has been documented in several model organisms, the genetic mechanisms of their adaptation are still poorly understood and debated, particularly the role of sexual reproduction. Here, we apply a population genomic approach to infer genome-wide patterns of selection among 13 isolates of Zymoseptoria tritici, a fungal pathogen characterized by extremely high genetic diversity, gene density, and recombination rates. We report that the genome of Z. tritici undergoes a high rate of adaptive substitutions, with 44% of nonsynonymous substitutions being adaptive on average. This fraction reaches 68% in so-called effector genes encoding determinants of pathogenicity, and the distribution of fitness effects differs in this class of genes as they undergo adaptive mutations with stronger positive fitness effects, but also more slightly deleterious mutations. Besides the globally high rate of adaptive substitutions, we report a negative relationship between pN/pS and the fine-scale recombination rate and a strong positive correlation between the rate of adaptive nonsynonymous substitutions (ωa) and recombination rate. This result suggests a pervasive role of both background selection and Hill-Robertson interference even in a species with an exceptionally high recombination rate (60 cM/Mb on average). While transposable elements (TEs) have been suggested to contribute to adaptation by creating compartments of fast-evolving genomic regions, we do not find a significant effect of TEs on the rate of adaptive mutations. Overall our study suggests that sexual recombination is a significant driver of genome evolution, even in rapidly evolving organisms subject to recurrent mutations with large positive effects.

Genomic Plasticity Mediated by Transposable Elements in the Plant Pathogenic Fungus Colletotrichum higginsianum

Phytopathogen genomes are under constant pressure to change, as pathogens are locked in an evolutionary arms race with their hosts, where pathogens evolve effector genes to manipulate their hosts, whereas the hosts evolve immune components to recognize the products of these genes. Colletotrichum higginsianum (Ch), a fungal pathogen with no known sexual morph, infects Brassicaceae plants including Arabidopsis thaliana. Previous studies revealed that Ch differs in its virulence toward various Arabidopsis thaliana ecotypes, indicating the existence of coevolutionary selective pressures. However, between-strain genomic variations in Ch have not been studied. Here, we sequenced and assembled the genome of a Ch strain, resulting in a highly contiguous genome assembly, which was compared with the chromosome-level genome assembly of another strain to identify genomic variations between strains. We found that the two closely related strains vary in terms of large-scale rearrangements, the existence of strain-specific regions, and effector candidate gene sets and that these variations are frequently associated with transposable elements (TEs). Ch has a compartmentalized genome consisting of gene-sparse, TE-dense regions with more effector candidate genes and gene-dense, TE-sparse regions harboring conserved genes. Additionally, analysis of the conservation patterns and syntenic regions of effector candidate genes indicated that the two strains vary in their effector candidate gene sets because of de novo evolution, horizontal gene transfer, or gene loss after divergence. Our results reveal mechanisms for generating genomic diversity in this asexual pathogen, which are important for understanding its adaption to hosts.

The Genome of the Fungal Pathogen Verticillium dahliae Reveals Extensive Bacterial to Fungal Gene Transfer

Horizontal gene transfer (HGT) involves the transmission of genetic material between distinct evolutionary lineages and can be an important source of biological innovation. Reports of interkingdom HGT to eukaryotic microbial pathogens have accumulated over recent years. Verticillium dahliae is a notorious plant pathogen that causes vascular wilt disease on hundreds of plant species, resulting in high economic losses every year. Previously, the effector gene Ave1 and a glucosyltransferase-encoding gene were identified as virulence factor-encoding genes that were proposed to be horizontally acquired from a plant and a bacterial donor, respectively. However, to what extent HGT contributed to the overall genome composition of V. dahliae remained elusive. Here, we systematically searched for evidence of interkingdom HGT events in the genome of V. dahliae and provide evidence for extensive horizontal gene acquisition from bacterial origin.

Determining Whether Geographic Origin and Potato Genotypes Shape the Population Structure of Phytophthora infestans in the Central Region of Colombia

Pathogen variation plays an important role in the dynamics of infectious diseases. In this study, the genetic variation of 279 Phytophthora infestans isolates was assessed using a combination of 12 microsatellite simple-sequence repeat markers. Isolates were collected from 11 different potato cultivars in 11 different geographic localities of the central region of Colombia. The objective of this study was to determine whether populations were differentiated by host genotype or geographic origin. Within a single clonal lineage, EC-1, 76 genotypes were detected. An analysis of molecular variance attributed most of the variation to differences within host genotypes rather than among the host genotypes, suggesting that host cultivars do not structure the populations of the pathogen. Furthermore, the lack of a genetic population structure according to the host cultivar was confirmed by all of the analyses, including the Bayesian clustering analysis and the minimum spanning network that used the Bruvo genetic distance, which suggested that there are no significant barriers to gene flow for P. infestans among potato cultivars. According to the geographic origin, the populations of P. infestans were also not structured, and most of the variation among the isolates was attributed to differences within localities. Only some but not all localities in the north and west of the central region of Colombia showed some genetic differentiation from the other regions. The absence of sexual reproduction of this pathogen in Colombia was also demonstrated. Important insights are discussed regarding the genetic population dynamics of the P. infestans populations of the central region of Colombia that were provided by the results. In Colombia, there is a high genetic variation within the EC-1 clonal lineage with closely related genotypes, none dominant, that coexist in a wide geographic area and on several potato cultivars.

Comparative Genomics of Aspergillus flavus S and L Morphotypes Yield Insights into Niche Adaptation

Aspergillus flavus, the primary causal agent for aflatoxin contamination on crops, consists of isolates with two distinct morphologies: isolates of the S morphotype produce numerous small sclerotia and lower numbers of conidia while isolates of the L morphotype produce fewer large sclerotia and abundant conidia. The morphotypes also differ in aflatoxin production with S isolates consistently producing high concentrations of aflatoxin, whereas L isolates range from atoxigenic to highly toxigenic. The production of abundant sclerotia by the S morphotype suggests adaptation for long-term survival in the soil, whereas the production of abundant conidia by the L morphotype suggests adaptation for aerial dispersal to the phyllosphere. To identify genomic changes that support differential niche adaption, the sequences of three S and three L morphotype isolates were compared. Differences in genome structure and gene content were identified between the morphotypes. A >530 kb inversion between the morphotypes affect a secondary metabolite gene cluster and a cutinase gene. The morphotypes also differed in proteins predicted to be involved in carbon/nitrogen metabolism, iron acquisition, antimicrobial defense, and evasion of host immunity. The S morphotype genomes contained more intact secondary metabolite clusters indicating there is higher selection pressure to maintain secondary metabolism in the soil and that it is not limited to aflatoxin production. The L morphotype genomes were enriched in amino acid transporters, suggesting efficient nitrogen transport may be critical in the nutrient limited phyllosphere. These findings indicate the genomes of the two morphotypes differ beyond developmental genes and have diverged as they adapted to their respective niches.

Single nucleus sequencing reveals evidence of inter-nucleus recombination in arbuscular mycorrhizal fungi

Eukaryotes thought to have evolved clonally for millions of years are referred to as ancient asexuals. The oldest group among these are the arbuscular mycorrhizal fungi (AMF), which are plant symbionts harboring hundreds of nuclei within one continuous cytoplasm. Some AMF strains (dikaryons) harbor two co-existing nucleotypes but there is no direct evidence that such nuclei recombine in this life-stage, as is expected for sexual fungi. Here, we show that AMF nuclei with distinct genotypes can undergo recombination. Inter-nuclear genetic exchange varies in frequency among strains, and despite recombination all nuclear genomes have an average similarity of at least 99.8%. The present study demonstrates that AMF can generate genetic diversity via meiotic-like processes in the absence of observable mating. The AMF dikaryotic life-stage is a primary source of nuclear variability in these organisms, highlighting its potential for strain enhancement of these symbionts.

High intraspecific genome diversity in the model arbuscular mycorrhizal symbiont Rhizophagus irregularis

Arbuscular mycorrhizal fungi (AMF) are known to improve plant fitness through the establishment of mycorrhizal symbioses. Genetic and phenotypic variations among closely related AMF isolates can significantly affect plant growth, but the genomic changes underlying this variability are unclear. To address this issue, we improved the genome assembly and gene annotation of the model strain Rhizophagus irregularis DAOM197198, and compared its gene content with five isolates of R. irregularis sampled in the same field. All isolates harbor striking genome variations, with large numbers of isolate-specific genes, gene family expansions, and evidence of interisolate genetic exchange. The observed variability affects all gene ontology terms and PFAM protein domains, as well as putative mycorrhiza-induced small secreted effector-like proteins and other symbiosis differentially expressed genes. High variability is also found in active transposable elements. Overall, these findings indicate a substantial divergence in the functioning capacity of isolates harvested from the same field, and thus their genetic potential for adaptation to biotic and abiotic changes. Our data also provide a first glimpse into the genome diversity that resides within natural populations of these symbionts, and open avenues for future analyses of plant-AMF interactions that link AMF genome variation with plant phenotype and fitness.

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.

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.

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.

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.

Gene expression polymorphism underpins evasion of host immunity in an asexual lineage of the Irish potato famine pathogen

BACKGROUND

Outbreaks caused by asexual lineages of fungal and oomycete pathogens are a continuing threat to crops, wild animals and natural ecosystems (Fisher MC, Henk DA, Briggs CJ, Brownstein JS, Madoff LC, McCraw SL, Gurr SJ, Nature 484:186-194, 2012; Kupferschmidt K, Science 337:636-638, 2012). However, the mechanisms underlying genome evolution and phenotypic plasticity in asexual eukaryotic microbes remain poorly understood (Seidl MF, Thomma BP, BioEssays 36:335-345, 2014). Ever since the 19th century Irish famine, the oomycete Phytophthora infestans has caused recurrent outbreaks on potato and tomato crops that have been primarily caused by the successive rise and migration of pandemic asexual lineages (Goodwin SB, Cohen BA, Fry WE, Proc Natl Acad Sci USA 91:11591-11595, 1994; Yoshida K, Burbano HA, Krause J, Thines M, Weigel D, Kamoun S, PLoS Pathog 10:e1004028, 2014; Yoshida K, Schuenemann VJ, Cano LM, Pais M, Mishra B, Sharma R, Lanz C, Martin FN, Kamoun S, Krause J, et al. eLife 2:e00731, 2013; Cooke DEL, Cano LM, Raffaele S, Bain RA, Cooke LR, Etherington GJ, Deahl KL, Farrer RA, Gilroy EM, Goss EM, et al. PLoS Pathog 8:e1002940, 2012). However, the dynamics of genome evolution within these clonal lineages have not been determined. The objective of this study was to use a comparative genomics and transcriptomics approach to determine the molecular mechanisms that underpin phenotypic variation within a clonal lineage of P. infestans.

RESULTS

Here, we reveal patterns of genomic and gene expression variation within a P. infestans asexual lineage by comparing strains belonging to the South American EC-1 clone that has dominated Andean populations since the 1990s (Yoshida K, Burbano HA, Krause J, Thines M, Weigel D, Kamoun S, PLoS Pathog 10e1004028, 2014; Yoshida K, Schuenemann VJ, Cano LM, Pais M, Mishra B, Sharma R, Lanz C, Martin FN, Kamoun S, Krause J, et al. eLife 2:e00731, 2013; Delgado RA, Monteros-Altamirano AR, Li Y, Visser RGF, van der Lee TAJ, Vosman B, Plant Pathol 62:1081-1088, 2013; Forbes GA, Escobar XC, Ayala CC, Revelo J, Ordonez ME, Fry BA, Doucett K, Fry WE, Phytopathology 87:375-380, 1997; Oyarzun PJ, Pozo A, Ordonez ME, Doucett K, Forbes GA, Phytopathology 88:265-271, 1998). We detected numerous examples of structural variation, nucleotide polymorphisms and loss of heterozygosity within the EC-1 clone. Remarkably, 17 genes are not expressed in one of the two EC-1 isolates despite apparent absence of sequence polymorphisms. Among these, silencing of an effector gene was associated with evasion of disease resistance conferred by a potato immune receptor.

CONCLUSIONS

Our findings highlight the molecular changes underpinning the exceptional genetic and phenotypic plasticity associated with host adaptation in a pandemic clonal lineage of a eukaryotic plant pathogen. We observed that the asexual P. infestans lineage EC-1 can exhibit phenotypic plasticity in the absence of apparent genetic mutations resulting in virulence on a potato carrying the Rpi-vnt1.1 gene. Such variant alleles may be epialleles that arose through epigenetic changes in the underlying genes.

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.

Karyotype Variability and Inter-Population Genomic Differences in Freshwater Ostracods (Crustacea) Showing Geographical Parthenogenesis

Transitions from sexual to asexual reproduction are often associated with polyploidy and increased chromosomal plasticity in asexuals. We investigated chromosomes in the freshwater ostracod species Eucypris virens (Jurine, 1820), where sexual, asexual and mixed populations can be found. Our initial karyotyping of multiple populations from Europe and North Africa, both sexual and asexual, revealed a striking variability in chromosome numbers. This would suggest that chromosomal changes are likely to be accelerated in asexuals because the constraints of meiosis are removed. Hence, we employed comparative genomic hybridization (CGH) within and among sexual and asexual populations to get insights into E. virens genome arrangements. CGH disclosed substantial genomic imbalances among the populations analyzed, and three patterns of genome arrangement between these populations: 1. Only putative ribosomal DNA (rDNA)-bearing regions were conserved in the two populations compared indicating a high sequence divergence between these populations. This pattern is comparable with our findings at the interspecies level of comparison; 2. Chromosomal regions were shared by both populations to a varying extent with a distinct copy number variation in pericentromeric and presumable rDNA-bearing regions. This indicates a different rate of evolution in repetitive sequences; 3. A mosaic pattern of distribution of genomic material that can be explained as non-reciprocal genetic introgression and evidence of a hybrid origin of these individuals. We show an overall increased chromosomal dynamics in E. virens that is complementary with available phylogenetic and population genetic data reporting highly differentiated diploid sexual and asexual lineages with a wide variety of genetic backgrounds.

The birth and death of effectors in rapidly evolving filamentous pathogen genomes

Plant pathogenic fungi and oomycetes are major risks to food security due to their evolutionary success in overcoming plant defences. Pathogens produce effectors to interfere with host defences and metabolism. These effectors are often encoded in rapidly evolving compartments of the genome. We review how effector genes emerged and were lost in pathogen genomes drawing on the links between effector evolution and chromosomal rearrangements. Some new effectors entered pathogen genomes via horizontal transfer or introgression. However, new effector functions also arose through gene duplication or from previously non-coding sequences. The evolutionary success of an effector is tightly linked to its transcriptional regulation during host colonization. Some effectors converged on an epigenetic control of expression imposed by genomic defences against transposable elements. Transposable elements were also drivers of effector diversification and loss that led to mosaics in effector presence-absence variation. Such effector mosaics within species was the foundation for rapid pathogen adaptation.

Vegetative compatibility groups partition variation in the virulence of Verticillium dahliae on strawberry

Verticillium dahliae infection of strawberry (Fragaria x ananassa) is a major cause of disease-induced wilting in soil-grown strawberries across the world. To understand what components of the pathogen are affecting disease expression, the presence of the known effector VdAve1 was screened in a sample of Verticillium dahliae isolates. Isolates from strawberry were found to contain VdAve1 and were divided into two major clades, based upon their vegetative compatibility groups (VCG); no UK strawberry isolates contained VdAve1. VC clade was strongly related to their virulence levels. VdAve1-containing isolates pathogenic on strawberry were found in both clades, in contrast to some recently published findings. On strawberry, VdAve1-containing isolates had significantly higher virulence during early infection, which diminished in significance as the infection progressed. Transformation of a virulent non-VdAve1 containing isolate, with VdAve1 was found neither to increase nor decrease virulence when inoculated on a susceptible strawberry cultivar. There are therefore virulence factors that are epistatic to VdAve1 and potentially multiple independent routes to high virulence on strawberry in V. dahliae lineages. Genome sequencing a subset of isolates across the two VCGs revealed that isolates were differentiated at the whole genome level and contained multiple changes in putative effector content, indicating that different clonal VCGs may have evolved different strategies for infecting strawberry, leading to different virulence levels in pathogenicity tests. It is therefore important to consider both clonal lineage and effector complement as the adaptive potential of each lineage will differ, even if they contain the same race determining effector.

Stress and sexual reproduction affect the dynamics of the wheat pathogen effector AvrStb6 and strobilurin resistance

Host resistance and fungicide treatments are cornerstones of plant-disease control. Here, we show that these treatments allow sex and modulate parenthood in the fungal wheat pathogen Zymoseptoria tritici. We demonstrate that the Z. tritici-wheat interaction complies with the gene-for-gene model by identifying the effector AvrStb6, which is recognized by the wheat resistance protein Stb6. Recognition triggers host resistance, thus implying removal of avirulent strains from pathogen populations. However, Z. tritici crosses on wheat show that sex occurs even with an avirulent parent, and avirulence alleles are thereby retained in subsequent populations. Crossing fungicide-sensitive and fungicide-resistant isolates under fungicide pressure results in a rapid increase in resistance-allele frequency. Isolates under selection always act as male donors, and thus disease control modulates parenthood. Modeling these observations for agricultural and natural environments reveals extended durability of host resistance and rapid emergence of fungicide resistance. Therefore, fungal sex has major implications for disease control.