Karyotype Variability in Plant-Pathogenic Fungi

The dynamics of fungal genome organization and its impact on host adaptation and antifungal resistance

Fungi are a diverse kingdom characterized by remarkable genomic plasticity that facilitates pathogenicity and adaptation to adverse environmental conditions. In this review, we delve into the dynamic organization of fungal genomes and its implications for host adaptation and antifungal resistance. We examine key features and the heterogeneity of genomes across different fungal species, including but not limited to their chromosome content, DNA composition, distribution and arrangement of their content across chromosomes, and other major traits. We further highlight how this variability in genomic traits influences their virulence and adaptation to adverse conditions. Fungal genomes exhibit large variations in size, gene content, and structural features, such as the abundance of transposable elements (TEs), compartmentalization into gene-rich and TE-rich regions, and the presence or absence of dispensable chromosomes. Genomic structural variations are equally diverse in fungi, ranging from whole-chromosome duplications that may enhance tolerance to antifungal compounds, to targeted deletion of effector encoding genes that may promote virulence. Finally, the often-overlooked fungal mitochondrial genomes can also affect virulence and resistance to fungicides. Such and other features of fungal genome organization are reviewed and discussed in the context of host-microbe interactions and antifungal resistance.

High-quality genome assembly and comparative analysis reveal extensive genomic variation in Talaromyces marneffei

Talaromyces marneffei is a dimorphic fungus that transitions from a filamentous form at 25 °C to a pathogenic yeast form at 37 °C, demonstrating pathogenicity mostly in immunocompromised individuals, such as those with human immunodeficiency virus/AIDS. Though it is one of the most severe infectious fungi in Southeast Asia, the lack of comprehensive genomic analysis has hindered advancement in strain differentiation, diagnosis and treatment. In this study, we assembled a high-quality genome of T. marneffei ATCC 18224, resulting in a 28.9 Mb genome distributed across 11 contigs, using third-generation Oxford Nanopore Technologies sequencing reads. Notably, we identified a strain-specific 740-kb segmental duplication in strain ATCC 18224, potentially mediated by inserting a Ty1/Copia long terminal repeat (LTR) retrotransposon. This segmental duplication includes various functional genes, with 75 differentially expressed during its dimorphic transition. Comparative genomic analysis revealed large-scale rearrangements in strains PM1 and 11CN-20-091, which were inconsistent with the phylogenomic trees of six T. marneffei strains and required further investigation. Additionally, we observed substantial genetic structural variations in LTR retrotransposons, particularly within the Ty1/Copia family, including two significant recent expansions in strain ATCC 18224. In summary, the identification and characterization of these extensive genomic structural variations in T. marneffei contribute to a deep understanding of its genetic diversity and will facilitate improvements in genotyping, classification and genomic surveillance.

From Lesions to Lessons: Two Decades of Filamentous Plant Pathogen Genomics

Many filamentous microorganisms, such as fungi and oomycetes, have evolved the ability to colonize plants and cause devastating crop diseases. Coevolutionary conflicts with their hosts have shaped the genomes of these plant pathogens. Over the past 20 years, genomics and genomics-enabled technologies have revealed remarkable diversity in genome size, architecture, and gene regulatory mechanisms. Technical and conceptual advances continue to provide novel insights into evolutionary dynamics, diversification of distinct genomic compartments, and facilitated molecular disease diagnostics. In this review, we discuss how genomics has advanced our understanding of genome organization and plant-pathogen coevolution and provide a perspective on future developments in the field. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Recombination landscape and karyotypic variations revealed by linkage mapping in the grapevine downy mildew pathogen Plasmopara viticola

Plasmopara viticola, the causal agent of grapevine downy mildew, is a biotrophic oomycete engaged in a tight coevolutionary relationship with its host. Rapid adaptation of the pathogen is favored by annual sexual reproduction that generates genotypic diversity. With the aim of studying the recombination landscape across the P. viticola genome, we generated 2 half-sibling F1 progenies (N = 189 and 162). Using targeted SNP sequencing, between 1,405 and 1,894 markers were included in parental linkage maps, and a consensus map was obtained by integrating 4,509 markers. The reference genome could be assembled into 17 pseudochromosomes, anchoring 88% of its physical length. We observed a strong collinearity between parental genomes and extensive synteny with the downy mildew Peronospora effusa. In the consensus map, the median recombination rate was 13.8 cM/Mb. The local recombination rate was highly variable along chromosomes, and recombination was suppressed in putative centromeric regions. Recombination rate was found negatively correlated with repeats' coverage and positively correlated with gene coverage. However, genes encoding secreted proteins and putative effectors were underrepresented in highly recombining regions. In both progenies, about 5% of the individuals presented karyotypic anomalies. Aneuploidies and triploidies almost exclusively originated from the male-transmitted chromosomes. Triploids resulted from fertilization by diploid gametes, but also from dispermy. Obligatory sexual reproduction each year may explain the lower level of karyotypic variation in P. viticola compared to other oomycetes. The linkage maps will be useful to guide future de novo chromosome-scale assemblies of P. viticola genomes and to perform forward genetics.

Complete Annotated Genome Assembly of Flax Pathogen Colletotrichum lini

Colletotrichum lini is a fungal pathogen of flax that can cause significant yield and quality losses. In this work, we obtained the first complete annotated genome assembly of the highly virulent C. lini strain #394-2. The nuclear genome consisted of ten core and two accessory chromosomes and had a length of 53.7 Mb. The mitochondrial genome was 39.1 kb. The assembly was obtained by the Canu-Racon ×2-Medaka-Polca algorithm using Oxford Nanopore Technologies and Illumina data. As a result of the annotation with the Illumina RNA-Seq data, 12,449 genes were identified. Potential signaling proteins were tested for effector functions and 550 effector proteins were predicted using EffectorP. The visualization of the effector protein localization revealed that the presence of effector proteins was associated with repeat-rich regions. A comparison of the genomic structure of C. lini with chromosome-level and complete assemblies of the genus Colletotrichum representatives revealed that the genomes of Colletotrichum species differed by the presence of chromosomal rearrangements. The obtained assembly expands the knowledge of the genomic structure of Colletotrichum species and provides the basis for further studies of C. lini, which will help to understand the virulence mechanisms and protect flax from anthracnose.

The landscape and predicted roles of structural variants in Fusarium graminearum genomes

Structural rearrangements, such as inversions, translocations, duplications, and large insertions and deletions, are large-scale genomic variants that can play an important role in shaping phenotypic variation and in genome adaptation and evolution. We used chromosomal-level assemblies from eight Fusarium graminearum isolates to study structural variants and their role in fungal evolution. We generated the assemblies of four of these genomes after Oxford Nanopore sequencing. A total of 87 inversions, 159 translocations, 245 duplications, 58,489 insertions, and 34,102 deletions were detected. Regions of high recombination rate are associated with structural rearrangements, and a significant proportion of inversions, translocations, and duplications overlap with the repeat content of the genome, suggesting recombination and repeat elements are major factors in the origin of structural rearrangements in F. graminearum. Large insertions and deletions introduce presence-absence polymorphisms for many genes, including secondary metabolite biosynthesis cluster genes and predicted effectors genes. Translocation events were found to be shuffling predicted effector-rich regions of the genomes and are likely contributing to the gain and loss of effectors facilitated by recombination. Breakpoints of some structural rearrangements fall within coding sequences and are likely altering the protein products. Structural rearrangements in F. graminearum thus have an important role to play in shaping pathogen-host interactions and broader evolution through genome reorganization, the introduction of presence-absence polymorphisms, and changing protein products and gene regulation.

Comparative genomics of the closely related fungal genera Cryptococcus and Kwoniella reveals karyotype dynamics and suggests evolutionary mechanisms of pathogenesis

In exploring the evolutionary trajectories of both pathogenesis and karyotype dynamics in fungi, we conducted a large-scale comparative genomic analysis spanning the Cryptococcus genus, encompassing both global human fungal pathogens and nonpathogenic species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species, covering virtually all known diversity within these genera. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at preadaptive pathogenic potential, our analysis found evidence of gene gain (via horizontal gene transfer) and gene loss in pathogenic Cryptococcus species, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the 2 genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5, or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes showed reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Overall, our findings advance our understanding of genetic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Frequent horizontal chromosome transfer between asexual fungal insect pathogens

Entire chromosomes are typically only transmitted vertically from one generation to the next. The horizontal transfer of such chromosomes has long been considered improbable, yet gained recent support in several pathogenic fungi where it may affect the fitness or host specificity. To date, it is unknown how these transfers occur, how common they are, and whether they can occur between different species. In this study, we show multiple independent instances of horizontal transfers of the same accessory chromosome between two distinct strains of the asexual entomopathogenic fungus Metarhizium robertsii during experimental co-infection of its insect host, the Argentine ant. Notably, only the one chromosome-but no other-was transferred from the donor to the recipient strain. The recipient strain, now harboring the accessory chromosome, exhibited a competitive advantage under certain host conditions. By phylogenetic analysis, we further demonstrate that the same accessory chromosome was horizontally transferred in a natural environment between M. robertsii and another congeneric insect pathogen, Metarhizium guizhouense. Hence, horizontal chromosome transfer is not limited to the observed frequent events within species during experimental infections but also occurs naturally across species. The accessory chromosome that was transferred contains genes that may be involved in its preferential horizontal transfer or support its establishment. These genes encode putative histones and histone-modifying enzymes, as well as putative virulence factors. Our study reveals that both intra- and interspecies horizontal transfer of entire chromosomes is more frequent than previously assumed, likely representing a not uncommon mechanism for gene exchange.

Comparative genomics of Cryptococcus and Kwoniella reveals pathogenesis evolution and contrasting karyotype dynamics via intercentromeric recombination or chromosome fusion

A large-scale comparative genomic analysis was conducted for the global human fungal pathogens within the Cryptococcus genus, compared to non-pathogenic Cryptococcus species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species of both genera, resulting in a dataset encompassing virtually all of their known diversity. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at pre-adaptive pathogenic potential, our analysis found evidence in pathogenic Cryptococcus species of specific examples of gene gain (via horizontal gene transfer) and gene loss, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the two genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5 or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes, underwent chromosome reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Taken together, our findings advance our understanding of genomic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Gapless Genome Assembly of Puccinia triticina Provides Insights into Chromosome Evolution in Pucciniales

Chromosome evolution drives species evolution, speciation, and adaptive radiation. Accurate genome assembly is crucial to understanding chromosome evolution of species, such as dikaryotic fungi. Rust fungi (Pucciniales) in dikaryons represent the largest group of plant pathogens, but the evolutionary process of adaptive radiation in Pucciniales remains poorly understood. Here, we report a gapless genome for the wheat leaf rust fungus Puccinia triticina determined using PacBio high-fidelity (HiFi) sequencing. This gapless assembly contains two sets of chromosomes, showing that one contig represents one chromosome. Comparisons of homologous chromosomes between the phased haplotypes revealed that highly frequent small-scale sequence divergence shapes haplotypic variation. Genome analyses of Puccinia triticina along with other rusts revealed that recent transposable element bursts and extensive segmental gene duplications synergistically highlight the evolution of chromosome structures. Comparative analysis of chromosomes indicated that frequent chromosomal rearrangements may act as a major contributor to rapid radiation of Pucciniales. This study presents the first gapless, phased assembly for a dikaryotic rust fungus and provides insights into adaptive evolution and species radiation in Pucciniales. IMPORTANCE Rust fungi (Pucciniales) are the largest group of plant pathogens. Adaptive radiation is a predominant feature in Pucciniales evolution. Chromosome evolution plays an important role in adaptive evolution. Accurate chromosome-scale assembly is required to understand the role of chromosome evolution in Pucciniales. We took advantage of HiFi sequencing to construct a gapless, phased genome for Puccinia triticina. Further analyses revealed that the evolution of chromosome structures in rust lineage is shaped by the combination of transposable element bursts and segmental gene duplications. Chromosome comparisons of Puccinia triticina and other rusts suggested that frequent chromosomal arrangements may make remarkable contributions to high species diversity of rust fungi. Our results present the first gapless genome for Pucciniales and shed light on the feature of chromosome evolution in Pucciniales.

Comprehending the dynamism of B chromosomes in their journey towards becoming unselfish

Investigated for more than a century now, B chromosomes (Bs) research has come a long way from Bs being considered parasitic or neutral to becoming unselfish and bringing benefits to their hosts. B chromosomes exist as accessory chromosomes along with the standard A chromosomes (As) across eukaryotic taxa. Represented singly or in multiple copies, B chromosomes are largely heterochromatic but also contain euchromatic and organellar segments. Although B chromosomes are derived entities, they follow their species-specific evolutionary pattern. B chromosomes fail to pair with the standard chromosomes during meiosis and vary in their number, size, composition and structure across taxa and ensure their successful transmission through non-mendelian mechanisms like mitotic, pre-meiotic, meiotic or post-meiotic drives, unique non-disjunction, self-pairing or even imparting benefits to the host when they lack drive. B chromosomes have been associated with cellular processes like sex determination, pathogenicity, resistance to pathogens, phenotypic effects, and differential gene expression. With the advancements in B-omics research, novel insights have been gleaned on their functions, some of which have been associated with the regulation of gene expression of A chromosomes through increased expression of miRNAs or differential expression of transposable elements located on them. The next-generation sequencing and emerging technologies will further likely unravel the cellular, molecular and functional behaviour of these enigmatic entities. Amidst the extensive fluidity shown by B chromosomes in their structural and functional attributes, we perceive that the existence and survival of B chromosomes in the populations most likely seem to be a trade-off between the drive efficiency and adaptive significance versus their adverse effects on reproduction.

Karyon: a computational framework for the diagnosis of hybrids, aneuploids, and other nonstandard architectures in genome assemblies

BACKGROUND

Recent technological developments have made genome sequencing and assembly highly accessible and widely used. However, the presence in sequenced organisms of certain genomic features such as high heterozygosity, polyploidy, aneuploidy, heterokaryosis, or extreme compositional biases can challenge current standard assembly procedures and result in highly fragmented assemblies. Hence, we hypothesized that genome databases must contain a nonnegligible fraction of low-quality assemblies that result from such type of intrinsic genomic factors.

FINDINGS

Here we present Karyon, a Python-based toolkit that uses raw sequencing data and de novo genome assembly to assess several parameters and generate informative plots to assist in the identification of nonchanonical genomic traits. Karyon includes automated de novo genome assembly and variant calling pipelines. We tested Karyon by diagnosing 35 highly fragmented publicly available assemblies from 19 different Mucorales (Fungi) species.

CONCLUSIONS

Our results show that 10 (28.57%) of the assemblies presented signs of unusual genomic configurations, suggesting that these are common, at least for some lineages within the Fungi.

Chromosomal polymorphism of the Ceratocystis fimbriata species complex in Brazil

Ceratocystis fimbriata is an important pathogen that causes wilt in several plant species. Despite the importance of this pathogen, knowledge about its karyotypic polymorphism and genomic architecture is limited. The main objective of this study was to investigate the karyotype of isolates of the C. fimbriata species complex from different host plants and geographical origins in Brazil. First, the identity of the isolates was confirmed conducting multilocus sequence analysis (MLSA) phylogeny using β-tubulin (TUBB), translation elongation factor 1α (TEF-1α) and mating-type (MAT1 and MAT2) gene sequences. To investigate the chromosomal polymorphism, two conditions of pulsed-field gel electrophoresis (PFGE) were used and the karyotypes of the isolates obtained. The retrotransposon-microsatellite amplified polymorphism (REMAP) molecular marker was utilized to assess the genetic variability among isolates. In the MLSA utilizing the concatenated gene sequences, Ceratocystis cacaofunesta and C. fimbriata formed separate clades, but considerable variation among C. fimbriata isolates was observed. Polymorphism in chromosome number and size was found, indicating the existence of genomic differences among isolates and occurrence of chromosomal rearrangements in the species complex. The number of chromosomes varied from seven to nine and the estimated minimum chromosome sizes were estimated to be between 2.7 to 6.0 Mbp. Small polymorphic chromosomes ranging from 1.2 to 1.8 Mbp were observed in all isolates, raising the hypothesis that they could be supernumerary chromosomes. REMAP analysis revealed a high genetic variability and that isolates from the same host tend to group together in a same cluster. Our results bring new insights into the chromosomal diversity and genome organization of the C. fimbriata complex.

Non-Mendelian transmission of accessory chromosomes in fungi

Non-Mendelian transmission has been reported for various genetic elements, ranging from small transposons to entire chromosomes. One prime example of such a transmission pattern are B chromosomes in plants and animals. Accessory chromosomes in fungi are similar to B chromosomes in showing presence/absence polymorphism and being non-essential. How these chromosomes are transmitted during meiosis is however poorly understood—despite their often high impact on the fitness of the host. For several fungal organisms, a non-Mendelian transmission or a mechanistically unique meiotic drive of accessory chromosomes have been reported. In this review, we provide an overview of the possible mechanisms that can cause the non-Mendelian transmission or meiotic drives of fungal accessory chromosomes. We compare processes responsible for the non-Mendelian transmission of accessory chromosomes for different fungal eukaryotes and discuss the structural traits of fungal accessory chromosomes affecting their meiotic transmission. We conclude that research on fungal accessory chromosomes, due to their small size, ease of sequencing, and epigenetic profiling, can complement the study of B chromosomes in deciphering factors that influence and regulate the non-Mendelian transmission of entire chromosomes.

Genetic Characteristics and Metabolic Interactions between Pseudocercospora fijiensis and Banana: Progress toward Controlling Black Sigatoka

The international importance of banana and severity of black Sigatoka disease have led to extensive investigations into the genetic characteristics and metabolic interactions between the Dothideomycete Pseudocercospora fijiensis and its banana host. P. fijiensis was shown to have a greatly expanded genome compared to other Dothideomycetes, due to the proliferation of retrotransposons. Genome analysis suggests the presence of dispensable chromosomes that may aid in fungal adaptation as well as pathogenicity. Genomic research has led to the characterization of genes and metabolic pathways involved in pathogenicity, including: secondary metabolism genes such as PKS10-2, genes for mitogen-activated protein kinases such as Fus3 and Slt2, and genes for cell wall proteins such as glucosyl phosphatidylinositol (GPI) and glycophospholipid surface (Gas) proteins. Studies conducted on resistance mechanisms in banana have documented the role of jasmonic acid and ethylene pathways. With the development of banana transformation protocols, strategies for engineering resistance include transgenes expressing antimicrobial peptides or hydrolytic enzymes as well as host-induced gene silencing (HIGS) targeting pathogenicity genes. Pseudocercospora fijiensis has been identified as having high evolutionary potential, given its large genome size, ability to reproduce both sexually and asexually, and long-distance spore dispersal. Thus, multiple control measures are needed for the sustainable control of black Sigatoka disease.

Parasexual recombination enables Aspergillus fumigatus to persist in cystic fibrosis

Aspergillus fumigatus is a saprobic fungus that causes a range of pulmonary diseases, some of which are characterised by fungal persistence such as is observed in cystic fibrosis (CF) patients. Creation of genetic variation is critical for A. fumigatus to adapt to the lung environment, but biofilm formation, especially in CF patients, may preclude mutational supply in A. fumigatus due to its confinement to the hyphal morphotype. We tested our hypothesis that genetic variation is created through parasexual recombination in chronic biofilms by phenotypic and genetic analysis of A. fumigatus isolates cultured from different origins.

As diploids are the hallmark of parasex, we screened 799 A. fumigatus isolates obtained from patients with CF, chronic pulmonary lung disease and acute invasive aspergillosis, and from the environment for spore size. Benomyl sensitivity, nuclear content measurements through fluorescence-activated cell sorting and scanning electron microscopy were used to confirm the diploid state of large size spores. Whole genome sequencing was used to characterise diploid-associated genetic variation.

We identified 11 diploids in isolates recovered from six of 11 (55%) CF patients and from one of 24 (4%) chronic aspergillosis patients, but not in 368 isolates from patients with acute Aspergillus infection and the environment. Diploid formation was associated with accumulation of mutations and variable haploid offspring including a voriconazole-resistant isolate.

Parasexual recombination allows A. fumigatus to adapt and persist in CF patients, and plays a role in azole resistance development. Our findings are highly significant for understanding the genetics and biology of A. fumigatus in the human lung.

Aspergillus fumigatus can undergo parasexual recombination in the lungs of chronically colonised individuals. This is highly relevant, as the parasexual cycle can facilitate genetic variation and, consequently, in-host adaptation.https://bit.ly/3kAtrXd

Whole genome analysis of the koa wilt pathogen (Fusarium oxysporum f. sp. koae) and the development of molecular tools for early detection and monitoring

BACKGROUND

Development and application of DNA-based methods to distinguish highly virulent isolates of Fusarium oxysporum f. sp. koae [Fo koae; cause of koa wilt disease on Acacia koa (koa)] will help disease management through early detection, enhanced monitoring, and improved disease resistance-breeding programs.

RESULTS

This study presents whole genome analyses of one highly virulent Fo koae isolate and one non-pathogenic F. oxysporum (Fo) isolate. These analyses allowed for the identification of putative lineage-specific DNA and predicted genes necessary for disease development on koa. Using putative chromosomes and predicted gene comparisons, Fo koae-exclusive, virulence genes were identified. The putative lineage-specific DNA included identified genes encoding products secreted in xylem (e. g., SIX1 and SIX6) that may be necessary for disease development on koa. Unique genes from Fo koae were used to develop pathogen-specific PCR primers. These diagnostic primers allowed target amplification in the characterized highly virulent Fo koae isolates but did not allow product amplification in low-virulence or non-pathogenic isolates of Fo. Thus, primers developed in this study will be useful for early detection and monitoring of highly virulent strains of Fo koae. Isolate verification is also important for disease resistance-breeding programs that require a diverse set of highly virulent Fo koae isolates for their disease-screening assays to develop disease-resistant koa.

CONCLUSIONS

These results provide the framework for understanding the pathogen genes necessary for koa wilt disease and the genetic variation of Fo koae populations across the Hawaiian Islands.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

Molecular analysis of mating type loci from the mycophenolic acid producer Penicillium brevicompactum: Phylogeny and MAT protein characterization suggest a cryptic sexual life cycle

The mycophenolic acid producing ascomycete Penicillium brevicompactum is considered to be an anamorphic (asexual) species, for which a sexual cycle was never observed. However, since recent reports of otherwise asexually propagating filamentous fungi have demonstrated a sexual cycle controlled by mating type loci, we carried out a molecular analysis of mating type loci from P. brevicompactum. Using data from extensive DNA sequencing analysis, we determined the mating type loci from 22 strains derived from various type culture collections. We found 8 strains carrying a MAT1-1 locus encoding a 362 amino acid alpha domain transcription factor. The other 14 possessed a MAT1-2 locus encoding a 298 amino acid HMG domain transcription factor. cDNA analysis confirmed that both mating type loci are transcriptionally expressed. The karyotype of six selected strains, determined using contour-clamped homogeneous electric field (CHEF) electrophoresis, demonstrated distinct differences in size and numbers of chromosomes between the strains investigated. Interestingly, our phylogenetic survey of 72 strains from 11 different Penicillium species revealed that MAT genes serve as excellent molecular markers to determine phylogenetic relationships among species closely related to P. brevicompactum. Based on our sequencing results, we constructed transformation vectors for site-specific deletion of mating type loci from two selected strains of opposite mating type. Complementation strains were constructed containing both the mating type locus deletion cassette and a MAT-egfp fusion gene. These strains were used for comparative phenotypic analyses between strains containing or lacking the mating type gene. Whereas all MAT1-2 strains were indistinguishable, the MAT1-1 and MAT1-1-1 deletion strains differed distinctly. The MAT1-1-1 deletion strain produced more conidiospores on solid media, but smaller pellets in liquid media. This is probably the consequence of fewer conidial germ tubes than with the wild type mating type strain. Finally, we showed that the MAT-EGPF fusion protein is localized to the nuclei and detectable in protein samples by Western analysis. Together, our results suggest that the asexually propagating fungus P. brevicompactum might be a heterothallic species with a cryptic sexual life cycle.

Genome compartmentalization predates species divergence in the plant pathogen genus Zymoseptoria

Background

Antagonistic co-evolution can drive rapid adaptation in pathogens and shape genome architecture. Comparative genome analyses of several fungal pathogens revealed highly variable genomes, for many species characterized by specific repeat-rich genome compartments with exceptionally high sequence variability. Dynamic genome structure may enable fast adaptation to host genetics. The wheat pathogen Zymoseptoria tritici with its highly variable genome, has emerged as a model organism to study genome evolution of plant pathogens. Here, we compared genomes of Z. tritici isolates and of sister species infecting wild grasses to address the evolution of genome composition and structure.

Results

Using long-read technology, we sequenced and assembled genomes of Z. ardabiliae, Z. brevis, Z. pseudotritici and Z. passerinii, together with two isolates of Z. tritici. We report a high extent of genome collinearity among Zymoseptoria species and high conservation of genomic, transcriptomic and epigenomic signatures of compartmentalization. We identify high gene content variability both within and between species. In addition, such variability is mainly limited to the accessory chromosomes and accessory compartments. Despite strong host specificity and non-overlapping host-range between species, predicted effectors are mainly shared among Zymoseptoria species, yet exhibiting a high level of presence-absence polymorphism within Z. tritici. Using in planta transcriptomic data from Z. tritici, we suggest different roles for the shared orthologs and for the accessory genes during infection of their hosts.

Conclusion

Despite previous reports of high genomic plasticity in Z. tritici, we describe here a high level of conservation in genomic, epigenomic and transcriptomic composition and structure across the genus Zymoseptoria. The compartmentalized genome allows the maintenance of a functional core genome co-occurring with a highly variable accessory genome.

Advances in understanding the evolution of fungal genome architecture

Diversity within the fungal kingdom is evident from the wide range of morphologies fungi display as well as the various ecological roles and industrial purposes they serve. Technological advances, particularly in long-read sequencing, coupled with the increasing efficiency and decreasing costs across sequencing platforms have enabled robust characterization of fungal genomes. These sequencing efforts continue to reveal the rampant diversity in fungi at the genome level. Here, we discuss studies that have furthered our understanding of fungal genetic diversity and genomic evolution. These studies revealed the presence of both small-scale and large-scale genomic changes. In fungi, research has recently focused on many small-scale changes, such as how hypermutation and allelic transmission impact genome evolution as well as how and why a few specific genomic regions are more susceptible to rapid evolution than others. High-throughput sequencing of a diverse set of fungal genomes has also illuminated the frequency, mechanisms, and impacts of large-scale changes, which include chromosome structural variation and changes in chromosome number, such as aneuploidy, polyploidy, and the presence of supernumerary chromosomes. The studies discussed herein have provided great insight into how the architecture of the fungal genome varies within species and across the kingdom and how modern fungi may have evolved from the last common fungal ancestor and might also pave the way for understanding how genomic diversity has evolved in all domains of life.

Genomic Instability in Fungal Plant Pathogens

Fungi and fungal-like organisms (oomycetes) that cause diseases in plants have impacted human communities for centuries and probably from the dawn of agriculture. In modern agriculture, there is a constant race between new strategies to manage fungal plant pathogens and their ability to adapt. An important component in this race is fungal genetic diversity. Mechanisms such as sexual and parasexual recombination that contribute to the creation of novel allele combinations in fungal plant pathogens are briefly discussed in the first part of this review. Advances in genomics have enabled the investigation of chromosomal aberrations of agriculturally important fungal isolates at the nucleotide level. Some of these cases are summarized in the second part of this review; it is claimed that the effect of chromosomal aberrations on pathogenicity should be studied mechanistically. More data on the effect of gene copy number variations on phenotypes that are relevant to agriculture are especially needed. Genome rearrangements through translocations have shaped the genome of fungal plant pathogens by creating lineage-specific chromosome territories encoding for genes participating in plant diseases. Pathogenicity chromosomes are unique cases of such lineage-specific genetic elements, interestingly these chromosomes can be transferred horizontally and thus transforming a non-pathogenic strain to a pathogenic one. The third part of this review describes our attempts to reveal mutators in fungal plant pathogens by identifying fungi that lack important DNA repair genes or respond to DNA damage in an unconventional way. We found that a group of fungal plant pathogens lack conserved genes that are needed for an important Holliday junction resolution pathway. In addition, in Fusarium oxysporum, the rate-limiting step in dNTP production is not induced under DNA replication stress. This is very different from organisms from bacteria to humans. It remains to be seen if these mechanisms promote genetic instability in fungal plant pathogens.

Karyotypic instability of endoprophase and mitotic cells of Amoeba sp. strain Cont from the "proteus-type" group (Amoebozoa, Euamoebida, Amoebidae)

We performed karyotyping of Amoeba sp. strain Cont. Based on the results of a cytological analysis, we concluded that the chromosome number of Amoeba sp. strain Cont in mitosis was unstable. In all cases they appeared to be hypergaploid (the basic chromosome number is 30), with monosomy of all chromosomes except four shortest ones. The presence of "extrachromosomes" in the nucleus could prolong until the beginning of the anaphase. It was only then that they were ejected from the nucleus and the euploidy (haploidy) was restored. The stage of endoprophase nucleus was revealed in the cell cycle of Amoeba sp. strain Cont. This stage has not yet been found in other amoebae from the "proteus-type" group that had been previously studied (A. proteus strain B and A. borokensis). The maximum number of endoreplication rounds in the strain Cont amoebae nuclear cycle was 4 or 5. The regular extrusion of chromosomes from the nucleus into the cytoplasm occurred in each of the endoreplication rounds. Comparative cytological analysis of A. proteus strain B, A. borokensis and Amoeba sp. strain Cont karyotypes indicated that strain Cont, though rather close to the former two amoebae, is actually a distinct species.

The architecture of the Plasmodiophora brassicae nuclear and mitochondrial genomes

Plasmodiophora brassicae is a soil-borne pathogen that attacks roots of cruciferous plants causing clubroot disease. The pathogen belongs to the Plasmodiophorida order in Phytomyxea. Here we used long-read SMRT technology to clarify the P. brassicae e3 genomic constituents along with comparative and phylogenetic analyses. Twenty contigs representing the nuclear genome and one mitochondrial (mt) contig were generated, together comprising 25.1 Mbp. Thirteen of the 20 nuclear contigs represented chromosomes from telomere to telomere characterized by [TTTTAGGG] sequences. Seven active gene candidates encoding synaptonemal complex-associated and meiotic-related protein homologs were identified, a finding that argues for possible genetic recombination events. The circular mt genome is large (114,663 bp), gene dense and intron rich. It shares high synteny with the mt genome of Spongospora subterranea, except in a unique 12 kb region delimited by shifts in GC content and containing tandem minisatellite- and microsatellite repeats with partially palindromic sequences. De novo annotation identified 32 protein-coding genes, 28 structural RNA genes and 19 ORFs. ORFs predicted in the repeat-rich region showed similarities to diverse organisms suggesting possible evolutionary connections. The data generated here form a refined platform for the next step involving functional analysis, all to clarify the complex biology of P. brassicae.

Fungal clones win the battle, but recombination wins the war

Clonal reproduction is common in fungi and fungal-like organisms during epidemics and invasion events. The success of clonal fungi shaped systems for their classification and some pathogens are tacitly treated as asexual. We argue that genetic recombination driven by sexual reproduction must be a starting hypothesis when dealing with fungi for two reasons: (1) Clones eventually crash because they lack adaptability; and (2) fungi find a way to exchange genetic material through recombination, whether sexual, parasexual, or hybridisation. Successful clones may prevail over space and time, but they are the product of recombination and the next successful clone will inevitably appear. Fungal pathogen populations are dynamic rather than static, and they need genetic recombination to adapt to a changing environment.

Effector gene reshuffling involves dispensable mini-chromosomes in the wheat blast fungus

Newly emerged wheat blast disease is a serious threat to global wheat production. Wheat blast is caused by a distinct, exceptionally diverse lineage of the fungus causing rice blast disease. Through sequencing a recent field isolate, we report a reference genome that includes seven core chromosomes and mini-chromosome sequences that harbor effector genes normally found on ends of core chromosomes in other strains. No mini-chromosomes were observed in an early field strain, and at least two from another isolate each contain different effector genes and core chromosome end sequences. The mini-chromosome is enriched in transposons occurring most frequently at core chromosome ends. Additionally, transposons in mini-chromosomes lack the characteristic signature for inactivation by repeat-induced point (RIP) mutation genome defenses. Our results, collectively, indicate that dispensable mini-chromosomes and core chromosomes undergo divergent evolutionary trajectories, and mini-chromosomes and core chromosome ends are coupled as a mobile, fast-evolving effector compartment in the wheat pathogen genome.

Effector gene reshuffling involves dispensable mini-chromosomes in the wheat blast fungus

Newly emerged wheat blast disease is a serious threat to global wheat production. Wheat blast is caused by a distinct, exceptionally diverse lineage of the fungus causing rice blast disease. Through sequencing a recent field isolate, we report a reference genome that includes seven core chromosomes and mini-chromosome sequences that harbor effector genes normally found on ends of core chromosomes in other strains. No mini-chromosomes were observed in an early field strain, and at least two from another isolate each contain different effector genes and core chromosome end sequences. The mini-chromosome is enriched in transposons occurring most frequently at core chromosome ends. Additionally, transposons in mini-chromosomes lack the characteristic signature for inactivation by repeat-induced point (RIP) mutation genome defenses. Our results, collectively, indicate that dispensable mini-chromosomes and core chromosomes undergo divergent evolutionary trajectories, and mini-chromosomes and core chromosome ends are coupled as a mobile, fast-evolving effector compartment in the wheat pathogen genome.

Fungal evolution: major ecological adaptations and evolutionary transitions

Fungi are a highly diverse group of heterotrophic eukaryotes characterized by the absence of phagotrophy and the presence of a chitinous cell wall. While unicellular fungi are far from rare, part of the evolutionary success of the group resides in their ability to grow indefinitely as a cylindrical multinucleated cell (hypha). Armed with these morphological traits and with an extremely high metabolical diversity, fungi have conquered numerous ecological niches and have shaped a whole world of interactions with other living organisms. Herein we survey the main evolutionary and ecological processes that have guided fungal diversity. We will first review the ecology and evolution of the zoosporic lineages and the process of terrestrialization, as one of the major evolutionary transitions in this kingdom. Several plausible scenarios have been proposed for fungal terrestralization and we here propose a new scenario, which considers icy environments as a transitory niche between water and emerged land. We then focus on exploring the main ecological relationships of Fungi with other organisms (other fungi, protozoans, animals and plants), as well as the origin of adaptations to certain specialized ecological niches within the group (lichens, black fungi and yeasts). Throughout this review we use an evolutionary and comparative-genomics perspective to understand fungal ecological diversity. Finally, we highlight the importance of genome-enabled inferences to envision plausible narratives and scenarios for important transitions.

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.

Karyotype evolution in Fusarium

The germ tube burst method (GTBM) was employed to examine karyotypes of 33 Fusarium species representative of 11 species complexes that span the phylogenetic breadth of the genus. The karyotypes revealed that the nucleolar organizing region (NOR), which includes the ribosomal rDNA region, was telomeric in the species where it was discernible. Variable karyotypes were detected in eight species due to variation in numbers of putative core and/or supernumerary chromosomes. The putative core chromosome number (CN) was most variable in the F. solani (CN = 9‒12) and F. buharicum (CN = 9+1 and 18-20) species complexes. Quantitative real-time PCR and genome sequence analysis rejected the hypothesis that the latter variation in CN was due to diploidization. The core CN in six other species complexes where two or more karyotypes were obtained was less variable or fixed. Karyotypes of 10 species in the sambucinum species complex, which is the most derived lineage of Fusarium, revealed that members of this complex possess the lowest CN in the genus. When viewed in context of the species phylogeny, karyotype evolution in Fusarium appears to have been dominated by a reduction in core CN in five closely related complexes that share a most recent common ancestor (tricinctum and incarnatum-equiseti CN = 8-9, chlamydosporum CN = 8, heterosporum CN = 7, sambucinum CN = 4-5) but not in the sister to these complexes (nisikadoi CN = 11, oxysporum CN = 11 and fujikuroi CN = 10-12). CN stability is best illustrated by the F. sambucinum subclade, where the only changes observed since it diverged from other fusaria appear to have involved two independent putative telomere to telomere fusions that reduced the core CN from five to four, once each in the sambucinum and graminearum subclades. Results of the present study indicate a core CN of 4 may be fixed in the latter subclade, which is further distinguished by the absence of putative supernumerary chromosomes. Karyotyping of fusaria in the not too distant future will be done by whole-genome sequencing such that each scaffold represents a complete chromosome from telomere to telomere. The CN data presented here should be of value to assist such full genome assembling.