A novel mode of chromosomal evolution peculiar to filamentous Ascomycete fungi

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

Retracing the evolution of Pneumocystis species, with a focus on the human pathogen Pneumocystis jirovecii

SUMMARYEvery human being is presumed to be infected by the fungus Pneumocystis jirovecii at least once in his or her lifetime. This fungus belongs to a large group of species that appear to exclusively infect mammals, with P. jirovecii being the only one known to cause disease in humans. The mystery of P. jirovecii origin and speciation is just beginning to unravel. Here, we provide a review of the major steps of P. jirovecii evolution. The Pneumocystis genus likely originated from soil or plant-associated organisms during the period of Cretaceous ~165 million years ago and successfully shifted to mammals. The transition coincided with a substantial loss of genes, many of which are related to the synthesis of nutrients that can be scavenged from hosts or cell wall components that could be targeted by the mammalian immune system. Following the transition, the Pneumocystis genus cospeciated with mammals. Each species specialized at infecting its own host. Host specialization is presumably built at least partially upon surface glycoproteins, whose protogene was acquired prior to the genus formation. P. jirovecii appeared at ~65 million years ago, overlapping with the emergence of the first primates. P. jirovecii and its sister species P. macacae, which infects macaques nowadays, may have had overlapping host ranges in the distant past. Clues from molecular clocks suggest that P. jirovecii did not cospeciate with humans. Molecular evidence suggests that Pneumocystis speciation involved chromosomal rearrangements and the mounting of genetic barriers that inhibit gene flow among species.

Template-Based Modelling of the Structure of Fungal Effector Proteins

The discovery of new fungal effector proteins is necessary to enable the screening of cultivars for disease resistance. Sequence-based bioinformatics methods have been used for this purpose, but only a limited number of functional effector proteins have been successfully predicted and subsequently validated experimentally. A significant obstacle is that many fungal effector proteins discovered so far lack sequence similarity or conserved sequence motifs. The availability of experimentally determined three-dimensional (3D) structures of a number of effector proteins has recently highlighted structural similarities amongst groups of sequence-dissimilar fungal effectors, enabling the search for similar structural folds amongst effector sequence candidates. We have applied template-based modelling to predict the 3D structures of candidate effector sequences obtained from bioinformatics predictions and the PHI-BASE database. Structural matches were found not only with ToxA- and MAX-like effector candidates but also with non-fungal effector-like proteins-including plant defensins and animal venoms-suggesting the broad conservation of ancestral structural folds amongst cytotoxic peptides from a diverse range of distant species. Accurate modelling of fungal effectors were achieved using RaptorX. The utility of predicted structures of effector proteins lies in the prediction of their interactions with plant receptors through molecular docking, which will improve the understanding of effector-plant interactions.

het-B allorecognition in Podospora anserina is determined by pseudo-allelic interaction of genes encoding a HET and lectin fold domain protein and a PII-like protein

Filamentous fungi display allorecognition genes that trigger regulated cell death (RCD) when strains of unlike genotype fuse. Podospora anserina is one of several model species for the study of this allorecognition process termed heterokaryon or vegetative incompatibility. Incompatibility restricts transmission of mycoviruses between isolates. In P. anserina, genetic analyses have identified nine incompatibility loci, termed het loci. Here we set out to clone the genes controlling het-B incompatibility. het-B displays two incompatible alleles, het-B1 and het-B2. We find that the het-B locus encompasses two adjacent genes, Bh and Bp that exist as highly divergent allelic variants (Bh1/Bh2 and Bp1/Bp2) in the incompatible haplotypes. Bh encodes a protein with an N-terminal HET domain, a cell death inducing domain bearing homology to Toll/interleukin-1 receptor (TIR) domains and a C-terminal domain with a predicted lectin fold. The Bp product is homologous to PII-like proteins, a family of small trimeric proteins acting as sensors of adenine nucleotides in bacteria. We show that although the het-B system appears genetically allelic, incompatibility is in fact determined by the non-allelic Bh1/Bp2 interaction while the reciprocal Bh2/Bp1 interaction plays no role in incompatibility. The highly divergent C-terminal lectin fold domain of BH determines recognition specificity. Population studies and genome analyses indicate that het-B is under balancing selection with trans-species polymorphism, highlighting the evolutionary significance of the two incompatible haplotypes. In addition to emphasizing anew the central role of TIR-like HET domains in fungal RCD, this study identifies novel players in fungal allorecognition and completes the characterization of the entire het gene set in that species.

Ancient gene linkages support ctenophores as sister to other animals

A central question in evolutionary biology is whether sponges or ctenophores (comb jellies) are the sister group to all other animals. These alternative phylogenetic hypotheses imply different scenarios for the evolution of complex neural systems and other animal-specific traits1-6. Conventional phylogenetic approaches based on morphological characters and increasingly extensive gene sequence collections have not been able to definitively answer this question7-11. Here we develop chromosome-scale gene linkage, also known as synteny, as a phylogenetic character for resolving this question12. We report new chromosome-scale genomes for a ctenophore and two marine sponges, and for three unicellular relatives of animals (a choanoflagellate, a filasterean amoeba and an ichthyosporean) that serve as outgroups for phylogenetic analysis. We find ancient syntenies that are conserved between animals and their close unicellular relatives. Ctenophores and unicellular eukaryotes share ancestral metazoan patterns, whereas sponges, bilaterians, and cnidarians share derived chromosomal rearrangements. Conserved syntenic characters unite sponges with bilaterians, cnidarians, and placozoans in a monophyletic clade to the exclusion of ctenophores, placing ctenophores as the sister group to all other animals. The patterns of synteny shared by sponges, bilaterians, and cnidarians are the result of rare and irreversible chromosome fusion-and-mixing events that provide robust and unambiguous phylogenetic support for the ctenophore-sister hypothesis. These findings provide a new framework for resolving deep, recalcitrant phylogenetic problems and have implications for our understanding of animal evolution.

Stepwise recombination suppression around the mating-type locus in an ascomycete fungus with self-fertile spores

Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.

Comparative genomics of five Valsa species gives insights on their pathogenicity evolution

Valsa is a genus of ascomycetes within the Valsaceae family. This family includes many wood destructive pathogens such as the well known Valsa mali and Valsa pyri which cause canker diseases in fruit trees and threaten the global fruit production. Lack of genomic information of this family is impeding our understandings about their evolution and genetic basis of their pathogenicity divergence. Here, we report genome assemblies of Valsa malicola, Valsa persoonii, and Valsa sordida which represent close relatives of Valsa mali and Valsa pyri with different host preferences. Comparative genomics analysis revealed that segmental rearrangements, inversions, and translocations frequently occurred among Valsa spp. genomes. Gene families that exhibited gene copy expansions tended to be associated with secondary metabolism, transmembrane transport, and pyrophosphatase activities. Orthologous genes in regions lost synteny exhibited significantly higher rate of synonymous substitution (KS) than those in regions retained synteny. Moreover, among these genes, membrane transporter families associated with antidrug (MFS, DHA) activities and nutrient transportation (SP and APCs) activities were significantly over-represented. Lineage specific synonymous substitution (KS) and nonsynonymous substitution (KA) analysis based on the phylogeny constructed from 11 fungal species identified a set of genes with selection signatures in Valsa clade and these genes were significantly enriched in functions associated with fatty acid beta-oxidation, DNA helicase activity, and ATPase activity. Furthermore, unique genes that possessed or retained by each of the five Valsa species are more likely part of the secondary metabolic (SM) gene clusters. SM gene clusters conserved across five Valsa species showed various degrees of diversification in both identity and completeness. All 11 syntenically conserved SM clusters showed differential expression during the infection of apple branch with Valsa mali suggesting involvements of secondary metabolism in the pathogenicity of Valsa species.

Differential Colonization of the Plant Vasculature Between Endophytic Versus Pathogenic Fusarium oxysporum Strains

Plant xylem colonization is the hallmark of vascular wilt diseases caused by phytopathogens within the Fusarium oxysporum species complex. Recently, xylem colonization has also been reported among endophytic F. oxysporum strains, resulting in some uncertainty. This study compares xylem colonization processes by pathogenic versus endophytic strains in Arabidopsis thaliana and Solanum lycopersicum, using Arabidopsis pathogen Fo5176, tomato pathogen Fol4287, and the endophyte Fo47, which can colonize both plant hosts. We observed that all strains were able to advance from epidermis to endodermis within 3 days postinoculation (dpi) and reached the root xylem at 4 dpi. However, this shared progression was restricted to lateral roots and the elongation zone of the primary root. Only pathogens reached the xylem above the primary-root maturation zone (PMZ). Related to the distinct colonization patterns, we also observed stronger induction of callose at the PMZ and lignin deposition at primary-lateral root junctions by the endophyte in both plants. This observation was further supported by stronger induction of Arabidopsis genes involved in callose and lignin biosynthesis during the endophytic colonization (Fo47) compared with the pathogenic interaction (Fo5176). Moreover, both pathogens encode more plant cell wall-degrading enzymes than the endophyte Fo47. Therefore, observed differences in callose and lignin deposition could be the combination of host production and the subsequent fungal degradation. In summary, this study demonstrates spatial differences between endophytic and pathogenic colonization, strongly suggesting that further investigations of molecular arm-races are needed to understand how plants differentiate friend from foe. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates,1,2,3 bilaterian animals,4,5 and filamentous fungi6). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ∼400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization.7 We found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ∼100 million years, the microsynteny of many genes-especially genes in metabolic clusters (e.g., in the GAL gene cluster8)-was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

Mycorrhiza-induced mycocypins of Laccaria bicolor are potent protease inhibitors with nematotoxic and collembola antifeedant activity

Fungivory of mycorrhizal hyphae has a significant impact on fungal fitness and, by extension, on nutrient transfer between fungi and host plants in natural ecosystems. Mycorrhizal fungi have therefore evolved an arsenal of chemical compounds that are hypothesized to protect the hyphal tissues from being eaten, such as the protease inhibitors mycocypins. The genome of the ectomycorrhizal fungus Laccaria bicolor has an unusually high number of mycocypin-encoding genes. We have characterized the evolution of this class of proteins, identified those induced by symbiosis with a host plant and characterized the biochemical properties of two upregulated L. bicolor mycocypins. More than half of L. bicolor mycocypin-encoding genes are differentially expressed during symbiosis or fruiting body formation. We show that two L. bicolor mycocypins that are strongly induced during symbiosis are cysteine protease inhibitors and exhibit similar but distinct localization in fungal tissues at different developmental stages and during interaction with a host plant. Moreover, we show that these L. bicolor mycocypins have toxic and feeding deterrent effect on nematodes and collembolans, respectively. Therefore, L. bicolor mycocypins may be part of a mechanism by which this species deters grazing by different members of the soil food web.

Telomere-to-Telomere Genome Sequences across a Single Genus Reveal Highly Variable Chromosome Rearrangement Rates but Absolute Stasis of Chromosome Number

Genome rearrangements in filamentous fungi are prevalent but little is known about the modalities of their evolution, in part because few complete genomes are available within a single genus. To address this, we have generated and compared 15 complete telomere-to-telomere genomes across the phylogeny of a single genus of filamentous fungi, Epichloë. We find that the striking distinction between gene-rich and repeat-rich regions previously reported for isolated species is ubiquitous across the Epichloë genus. We built a species phylogeny from single-copy gene orthologs to provide a comparative framing to study chromosome composition and structural change through evolutionary time. All Epichloë genomes have exactly seven nuclear chromosomes, but despite this conserved ploidy, analyses reveal low synteny and substantial rearrangement of gene content across the genus. These rearrangements are highly lineage-dependent, with most occurring over short evolutionary distances, with long periods of structural stasis. Quantification of chromosomal rearrangements shows they are uncorrelated with numbers of substitutions and evolutionary distances, suggesting that different modes of evolution are acting to create nucleotide and chromosome-scale changes.

A chromosome-scale genome assembly of the tomato pathogen Cladosporium fulvum reveals a compartmentalized genome architecture and the presence of a dispensable chromosome

Cladosporium fulvum is a fungal pathogen that causes leaf mould of tomato. The reference genome of this pathogen was released in 2012 but its high repetitive DNA content prevented a contiguous assembly and further prohibited the analysis of its genome architecture. In this study, we combined third generation sequencing technology with the Hi-C chromatin conformation capture technique, to produce a high-quality and near complete genome assembly and gene annotation of a Race 5 isolate of C. fulvum. The resulting genome assembly contained 67.17 Mb organized into 14 chromosomes (Chr1-to-Chr14), all of which were assembled telomere-to-telomere. The smallest of the chromosomes, Chr14, is only 460 kb in size and contains 25 genes that all encode hypothetical proteins. Notably, PCR assays revealed that Chr14 was absent in 19 out of 24 isolates of a world-wide collection of C. fulvum, indicating that Chr14 is dispensable. Thus, C. fulvum is currently the second species of Capnodiales shown to harbour dispensable chromosomes. The genome of C. fulvum Race 5 is 49.7 % repetitive and contains 14 690 predicted genes with an estimated completeness of 98.9%, currently one of the highest among the Capnodiales. Genome structure analysis revealed a compartmentalized architecture composed of gene-dense and repeat-poor regions interspersed with gene-sparse and repeat-rich regions. Nearly 39.2 % of the C. fulvum Race 5 genome is affected by Repeat-Induced Point (RIP) mutations and evidence of RIP leakage toward non-repetitive regions was observed in all chromosomes, indicating the RIP plays an important role in the evolution of this pathogen. Finally, 345 genes encoding candidate effectors were identified in C. fulvum Race 5, with a significant enrichment of their location in gene-sparse regions, in accordance with the 'two-speed genome' model of evolution. Overall, the new reference genome of C. fulvum presents several notable features and is a valuable resource for studies in plant pathogens.

An automated and combinative method for the predictive ranking of candidate effector proteins of fungal plant pathogens

Fungal plant-pathogens promote infection of their hosts through the release of 'effectors'-a broad class of cytotoxic or virulence-promoting molecules. Effectors may be recognised by resistance or sensitivity receptors in the host, which can determine disease outcomes. Accurate prediction of effectors remains a major challenge in plant pathology, but if achieved will facilitate rapid improvements to host disease resistance. This study presents a novel tool and pipeline for the ranking of predicted effector candidates-Predector-which interfaces with multiple software tools and methods, aggregates disparate features that are relevant to fungal effector proteins, and applies a pairwise learning to rank approach. Predector outperformed a typical combination of secretion and effector prediction methods in terms of ranking performance when applied to a curated set of confirmed effectors derived from multiple species. We present Predector ( https://github.com/ccdmb/predector ) as a useful tool for the ranking of predicted effector candidates, which also aggregates and reports additional supporting information relevant to effector and secretome prediction in a simple, efficient, and reproducible manner.

Remote homology clustering identifies lowly conserved families of effector proteins in plant-pathogenic fungi

Plant diseases caused by fungal pathogens are typically initiated by molecular interactions between 'effector' molecules released by a pathogen and receptor molecules on or within the plant host cell. In many cases these effector-receptor interactions directly determine host resistance or susceptibility. The search for fungal effector proteins is a developing area in fungal-plant pathology, with more than 165 distinct confirmed fungal effector proteins in the public domain. For a small number of these, novel effectors can be rapidly discovered across multiple fungal species through the identification of known effector homologues. However, many have no detectable homology by standard sequence-based search methods. This study employs a novel comparison method (RemEff) that is capable of identifying protein families with greater sensitivity than traditional homology-inference methods, leveraging a growing pool of confirmed fungal effector data to enable the prediction of novel fungal effector candidates by protein family association. Resources relating to the RemEff method and data used in this study are available from https://figshare.com/projects/Effector_protein_remote_homology/87965.

Metatranscriptomic Comparison of Endophytic and Pathogenic Fusarium-Arabidopsis Interactions Reveals Plant Transcriptional Plasticity

Plants are continuously exposed to beneficial and pathogenic microbes, but how plants recognize and respond to friends versus foes remains poorly understood. Here, we compared the molecular response of Arabidopsis thaliana independently challenged with a Fusarium oxysporum endophyte Fo47 versus a pathogen Fo5176. These two F. oxysporum strains share a core genome of about 46 Mb, in addition to 1,229 and 5,415 unique accessory genes. Metatranscriptomic data reveal a shared pattern of expression for most plant genes (about 80%) in responding to both fungal inoculums at all timepoints from 12 to 96 h postinoculation (HPI). However, the distinct responding genes depict transcriptional plasticity, as the pathogenic interaction activates plant stress responses and suppresses functions related to plant growth and development, while the endophytic interaction attenuates host immunity but activates plant nitrogen assimilation. The differences in reprogramming of the plant transcriptome are most obvious in 12 HPI, the earliest timepoint sampled, and are linked to accessory genes in both fungal genomes. Collectively, our results indicate that the A. thaliana and F. oxysporum interaction displays both transcriptome conservation and plasticity in the early stages of infection, providing insights into the fine-tuning of gene regulation underlying plant differential responses to fungal endophytes and pathogens.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Chromosome-level genome assembly and manually-curated proteome of model necrotroph Parastagonospora nodorum Sn15 reveals a genome-wide trove of candidate effector homologs, and redundancy of virulence-related functions within an accessory chromosome

Background

The fungus Parastagonospora nodorum causes septoria nodorum blotch (SNB) of wheat (Triticum aestivum) and is a model species for necrotrophic plant pathogens. The genome assembly of reference isolate Sn15 was first reported in 2007. P. nodorum infection is promoted by its production of proteinaceous necrotrophic effectors, three of which are characterised – ToxA, Tox1 and Tox3.

Results

A chromosome-scale genome assembly of P. nodorum Australian reference isolate Sn15, which combined long read sequencing, optical mapping and manual curation, produced 23 chromosomes with 21 chromosomes possessing both telomeres. New transcriptome data were combined with fungal-specific gene prediction techniques and manual curation to produce a high-quality predicted gene annotation dataset, which comprises 13,869 high confidence genes, and an additional 2534 lower confidence genes retained to assist pathogenicity effector discovery. Comparison to a panel of 31 internationally-sourced isolates identified multiple hotspots within the Sn15 genome for mutation or presence-absence variation, which was used to enhance subsequent effector prediction. Effector prediction resulted in 257 candidates, of which 98 higher-ranked candidates were selected for in-depth analysis and revealed a wealth of functions related to pathogenicity. Additionally, 11 out of the 98 candidates also exhibited orthology conservation patterns that suggested lateral gene transfer with other cereal-pathogenic fungal species. Analysis of the pan-genome indicated the smallest chromosome of 0.4 Mbp length to be an accessory chromosome (AC23). AC23 was notably absent from an avirulent isolate and is predominated by mutation hotspots with an increase in non-synonymous mutations relative to other chromosomes. Surprisingly, AC23 was deficient in effector candidates, but contained several predicted genes with redundant pathogenicity-related functions.

Conclusions

We present an updated series of genomic resources for P. nodorum Sn15 – an important reference isolate and model necrotroph – with a comprehensive survey of its predicted pathogenicity content.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07699-8.

The Architecture of Metabolism Maximizes Biosynthetic Diversity in the Largest Class of Fungi

Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or “specialized” metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

Fungal natural products galaxy: Biochemistry and molecular genetics toward blockbuster drugs discovery

Secondary metabolites synthesized by fungi have become a precious source of inspiration for the design of novel drugs. Indeed, fungi are prolific producers of fascinating, diverse, structurally complex, and low-molecular-mass natural products with high therapeutic leads, such as novel antimicrobial compounds, anticancer compounds, immunosuppressive agents, among others. Given that these microorganisms possess the extraordinary capacity to secrete diverse chemical scaffolds, they have been highly exploited by the giant pharma companies to generate small molecules. This has been made possible because the isolation of metabolites from fungal natural sources is feasible and surpasses the organic synthesis of compounds, which otherwise remains a significant bottleneck in the drug discovery process. Here in this comprehensive review, we have discussed recent studies on different fungi (pathogenic, non-pathogenic, commensal, and endophytic/symbiotic) from different habitats (terrestrial and marines), the specialized metabolites they biosynthesize, and the drugs derived from these specialized metabolites. Moreover, we have unveiled the logic behind the biosynthesis of vital chemical scaffolds, such as NRPS, PKS, PKS-NRPS hybrid, RiPPS, terpenoids, indole alkaloids, and their genetic mechanisms. Besides, we have provided a glimpse of the concept behind mycotoxins, virulence factor, and host immune response based on fungal infections.

Telomere length de novo assembly of all 7 chromosomes and mitogenome sequencing of the model entomopathogenic fungus, Metarhizium brunneum, by means of a novel assembly pipeline

Background

More accurate and complete reference genomes have improved understanding of gene function, biology, and evolutionary mechanisms. Hybrid genome assembly approaches leverage benefits of both long, relatively error-prone reads from third-generation sequencing technologies and short, accurate reads from second-generation sequencing technologies, to produce more accurate and contiguous de novo genome assemblies in comparison to using either technology independently. In this study, we present a novel hybrid assembly pipeline that allowed for both mitogenome de novo assembly and telomere length de novo assembly of all 7 chromosomes of the model entomopathogenic fungus, Metarhizium brunneum.

Results

The improved assembly allowed for better ab initio gene prediction and a more BUSCO complete proteome set has been generated in comparison to the eight current NCBI reference Metarhizium spp. genomes. Remarkably, we note that including the mitogenome in ab initio gene prediction training improved overall gene prediction. The assembly was further validated by comparing contig assembly agreement across various assemblers, assessing the assembly performance of each tool. Genomic synteny and orthologous protein clusters were compared between Metarhizium brunneum and three other Hypocreales species with complete genomes, identifying core proteins, and listing orthologous protein clusters shared uniquely between the two entomopathogenic fungal species, so as to further facilitate the understanding of molecular mechanisms underpinning fungal-insect pathogenesis.

Conclusions

The novel assembly pipeline may be used for other haploid fungal species, facilitating the need to produce high-quality reference fungal genomes, leading to better understanding of fungal genomic evolution, chromosome structuring and gene regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07390-y.

Comparative genomic analysis of iprodione-degrading Paenarthrobacter strains reveals the iprodione catabolic molecular mechanism in Paenarthrobacter sp. strain YJN-5

Degradation of the fungicide iprodione by the Paenarthrobacter sp. strain YJN-5 is initiated via hydrolysis of its N1 amide bond to form N-(3,5-dichlorophenyl)-2,4-dioxoimidazolidine. In this study, another iprodione-degrading strain, Paenarthrobacter sp. YJN-D, which harbours the same metabolic pathway as strain YJN-5 was isolated and characterized. The genes that encode the conserved iprodione catabolic pathway were identified based on comparative analysis of the genomes of the two iprodione-degrading Paenarthrobacter sp. and subsequent experimental validation. These genes include an amidase gene, ipaH (previously reported in AEM e01150-18); a deacetylase gene, ddaH, which is responsible for hydantoin ring cleavage of N-(3,5-dichlorophenyl)-2,4-dioxoimidazolidine, and a hydrolase gene, duaH, which is responsible for cleavage of the urea side chain of (3,5-dichlorophenylurea)acetic acid, thus yielding 3,5-dichloroaniline as the end product. These iprodione-catabolic genes are distributed on three plasmids in strain YJN-5 and are highly conserved between the two iprodione-degrading Paenarthrobacter strains. However, only the ipaH gene is flanked by a mobile genetic element. Two iprodione degradation cassettes bearing ipaH-ddaH-duaH were constructed and expressed in strains Pseudomonas putida KT2440 and Bacillus subtilis SCK6 respectively. Our findings enhance the current understanding of the microbial degradation mechanism of iprodione.

Fibrillarin evolution through the Tree of Life: Comparative genomics and microsynteny network analyses provide new insights into the evolutionary history of Fibrillarin

Fibrillarin (FIB), a methyltransferase essential for life in the vast majority of eukaryotes, is involved in methylation of rRNA required for proper ribosome assembly, as well as methylation of histone H2A of promoter regions of rRNA genes. RNA viral progression that affects both plants and animals requires FIB proteins. Despite the importance and high conservation of fibrillarins, there little is known about the evolutionary dynamics of this small gene family. We applied a phylogenomic microsynteny-network approach to elucidate the evolutionary history of FIB proteins across the Tree of Life. We identified 1063 non-redundant FIB sequences across 1049 completely sequenced genomes from Viruses, Bacteria, Archaea, and Eukarya. FIB is a highly conserved single-copy gene through Archaea and Eukarya lineages, except for plants, which have a gene family expansion due to paleopolyploidy and tandem duplications. We found a high conservation of the FIB genomic context during plant evolution. Surprisingly, FIB in mammals duplicated after the Eutheria split (e.g., ruminants, felines, primates) from therian mammals (e.g., marsupials) to form two main groups of sequences, the FIB and FIB-like groups. The FIB-like group transposed to another genomic context and remained syntenic in all the eutherian mammals. This transposition correlates with differences in the expression patterns of FIB-like proteins and with elevated Ks values potentially due to reduced evolutionary constraints of the duplicated copy. Our results point to a unique evolutionary event in mammals, between FIB and FIB-like genes, that led to non-redundant roles of the vital processes in which this protein is involved.

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.

Accessory Chromosomes in Fusarium oxysporum

Most genomes within the species complex of Fusarium oxysporum are organized into two compartments: the core chromosomes (CCs) and accessory chromosomes (ACs). As opposed to CCs, which are conserved and vertically transmitted to carry out essential housekeeping functions, lineage- or strain-specific ACs are believed to be initially horizontally acquired through unclear mechanisms. These two genomic compartments are different in terms of gene density, the distribution of transposable elements, and epigenetic markers. Although common in eukaryotes, the functional importance of ACs is uniquely emphasized among fungal species, specifically in relationship to fungal pathogenicity and their adaptation to diverse hosts. With a focus on the cross-kingdom fungal pathogen F. oxysporum, this review provides a summary of the differences between CCs and ACs based on current knowledge of gene functions, genome structures, and epigenetic signatures, and explores the transcriptional crosstalk between the core and accessory genomes.

CSBFinder: discovery of colinear syntenic blocks across thousands of prokaryotic genomes

MOTIVATION

Identification of conserved syntenic blocks across microbial genomes is important for several problems in comparative genomics such as gene annotation, study of genome organization and evolution and prediction of gene interactions. Current tools for syntenic block discovery do not scale up to the large quantity of prokaryotic genomes available today.

RESULTS

We present a novel methodology for the discovery, ranking and taxonomic distribution analysis of colinear syntenic blocks (CSBs)-groups of genes that are consistently located close to each other, in the same order, across a wide range of taxa. We present an efficient algorithm that identifies CSBs in large genomic datasets. The algorithm is implemented and incorporated in a novel tool with a graphical user interface, denoted CSBFinder, that ranks the discovered CSBs according to a probabilistic score and clusters them to families according to their gene content similarity. We apply CSBFinder to data mine 1487 prokaryotic genomes including chromosomes and plasmids. For post-processing analysis, we generate heatmaps for visualizing the distribution of CSB family members across various taxa. We exemplify the utility of CSBFinder in operon prediction, in deciphering unknown gene function and in taxonomic analysis of colinear syntenic blocks.

AVAILABILITY AND IMPLEMENTATION

CSBFinder software and code are publicly available at https://github.com/dinasv/CSBFinder.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A chromosome-scale assembly of the smallest Dothideomycete genome reveals a unique genome compaction mechanism in filamentous fungi

Background

The wide variation in the size of fungal genomes is well known, but the reasons for this size variation are less certain. Here, we present a chromosome-scale assembly of ectophytic Peltaster fructicola, a surface-dwelling extremophile, based on long-read DNA sequencing technology, to assess possible mechanisms associated with genome compaction.

Results

At 18.99 million bases (Mb), P. fructicola possesses one of the smallest known genomes sequence among filamentous fungi. The genome is highly compact relative to other fungi, with substantial reductions in repeat content, ribosomal DNA copies, tRNA gene quantity, and intron sizes, as well as intergenic lengths and the size of gene families. Transposons take up just 0.05% of the entire genome, and no full-length transposon was found. We concluded that reduced genome sizes in filamentous fungi such as P. fructicola, Taphrina deformans and Pneumocystis jirovecii occurred through reduction in ribosomal DNA copy number and reduced intron sizes. These dual mechanisms contrast with genome reduction in the yeast fungus Saccharomyces cerevisiae, whose small and compact genome is associated solely with intron loss.

Conclusions

Our results reveal a unique genomic compaction architecture of filamentous fungi inhabiting plant surfaces, and broaden the understanding of the mechanisms associated with compaction of fungal genomes.

Fungal evolution: diversity, taxonomy and phylogeny of the Fungi

The fungal kingdom comprises a hyperdiverse clade of heterotrophic eukaryotes characterized by the presence of a chitinous cell wall, the loss of phagotrophic capabilities and cell organizations that range from completely unicellular monopolar organisms to highly complex syncitial filaments that may form macroscopic structures. Fungi emerged as a 'Third Kingdom', embracing organisms that were outside the classical dichotomy of animals versus vegetals. The taxonomy of this group has a turbulent history that is only now starting to be settled with the advent of genomics and phylogenomics. We here review the current status of the phylogeny and taxonomy of fungi, providing an overview of the main defined groups. Based on current knowledge, nine phylum-level clades can be defined: Opisthosporidia, Chytridiomycota, Neocallimastigomycota, Blastocladiomycota, Zoopagomycota, Mucoromycota, Glomeromycota, Basidiomycota and Ascomycota. For each group, we discuss their main traits and their diversity, focusing on the evolutionary relationships among the main fungal clades. We also explore the diversity and phylogeny of several groups of uncertain affinities and the main phylogenetic and taxonomical controversies and hypotheses in the field.

Evolutionary and functional patterns of shared gene neighbourhood in fungi

Gene clusters comprise genomically co-localized and potentially co-regulated genes that tend to be conserved across species. In eukaryotes, multiple examples of metabolic gene clusters are known, particularly among fungi and plants. However, little is known about how gene clustering patterns vary among taxa or with respect to functional roles. Furthermore, mechanisms of the formation, maintenance and evolution of gene clusters remain unknown. We surveyed 341 fungal genomes to discover gene clusters shared by different species, independently of their functions. We inferred 12,120 cluster families, which comprised roughly one third of the gene space and were enriched in genes associated with diverse cellular functions. Additionally, most clusters did not encode transcription factors, suggesting that they are regulated distally. We used phylogenomics to characterize the evolutionary history of these clusters. We found that most clusters originated once and were transmitted vertically, coupled to differential loss. However, convergent evolution-that is, independent appearance of the same cluster-was more prevalent than anticipated. Finally, horizontal gene transfer of entire clusters was somewhat restricted, with the exception of those associated with secondary metabolism. Altogether, our results provide insights on the evolution of gene clustering as well as a broad catalogue of evolutionarily conserved gene clusters whose function remains to be elucidated.

Variable genome evolution in fungi after transposon-mediated amplification of a housekeeping gene

Background

Transposable elements (TEs) can be key drivers of evolution, but the mechanisms and scope of how they impact gene and genome function are largely unknown. Previous analyses revealed that TE-mediated gene amplifications can have variable effects on fungal genomes, from inactivation of function to production of multiple active copies. For example, a DNA methyltransferase gene in the wheat pathogen Zymoseptoria tritici (synonym Mycosphaerella graminicola) was amplified to tens of copies, all of which were inactivated by Repeat-Induced Point mutation (RIP) including the original, resulting in loss of cytosine methylation. In another wheat pathogen, Pyrenophora tritici-repentis, a histone H3 gene was amplified to tens of copies with little evidence of RIP, leading to many potentially active copies. To further test the effects of transposon-aided gene amplifications on genome evolution and architecture, the repetitive fraction of the significantly expanded genome of the banana pathogen, Pseudocercospora fijiensis, was analyzed in greater detail.

Results

These analyses identified a housekeeping gene, histone H3, which was captured and amplified to hundreds of copies by a hAT DNA transposon, all of which were inactivated by RIP, except for the original. In P. fijiensis the original H3 gene probably was not protected from RIP, but most likely was maintained intact due to strong purifying selection. Comparative analyses revealed that a similar event occurred in five additional genomes representing the fungal genera Cercospora, Pseudocercospora and Sphaerulina.

Conclusions

These results indicate that the interplay of TEs and RIP can result in different and unpredictable fates of amplified genes, with variable effects on gene and genome evolution.

Electronic supplementary material

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

Global population genomics of the forest pathogen Dothistroma septosporum reveal chromosome duplications in high dothistromin-producing strains

Dothistroma needle blight is one of the most devastating pine tree diseases worldwide. New and emerging epidemics have been frequent over the last 25 years, particularly in the Northern Hemisphere, where they are in part associated with changing weather patterns. One of the main Dothistroma needle blight pathogens, Dothistroma septosporum, has a global distribution but most molecular plant pathology research has been confined to Southern Hemisphere populations that have limited genetic diversity. Extensive genomic and transcriptomic data are available for a D. septosporum reference strain from New Zealand, where an introduced clonal population of the pathogen predominates. Due to the global importance of this pathogen, we determined whether the genome of this reference strain is representative of the species worldwide by sequencing the genomes of 18 strains sampled globally from different pine hosts. Genomic polymorphism shows substantial variation within the species, clustered into two distinct groups of strains with centres of diversity in Central and South America. A reciprocal chromosome translocation uniquely identifies the New Zealand strains. Globally, strains differ in their production of the virulence factor dothistromin, with extremely high production levels in strain ALP3 from Germany. Comparisons with the New Zealand reference revealed that several strains are aneuploids; for example, ALP3 has duplications of three chromosomes. Increased gene copy numbers therefore appear to contribute to increased production of dothistromin, emphasizing that studies of population structure are a necessary adjunct to functional analyses of genetic polymorphisms to identify the molecular basis of virulence in this important forest pathogen.

Physiological and comparative genomic analysis of new isolated yeasts Spathaspora sp. JA1 and Meyerozyma caribbica JA9 reveal insights into xylitol production

Xylitol is a five-carbon polyol of economic interest that can be produced by microbial xylose reduction from renewable resources. The current study sought to investigate the potential of two yeast strains, isolated from Brazilian Cerrado biome, in the production of xylitol as well as the genomic characteristics that may impact this process. Xylose conversion capacity by the new isolates Spathaspora sp. JA1 and Meyerozyma caribbica JA9 was evaluated and compared with control strains on xylose and sugarcane biomass hydrolysate. Among the evaluated strains, Spathaspora sp. JA1 was the strongest xylitol producer, reaching product yield and productivity as high as 0.74 g/g and 0.20 g/(L.h) on xylose, and 0.58 g/g and 0.44 g/(L.h) on non-detoxified hydrolysate. Genome sequences of Spathaspora sp. JA1 and M. caribbica JA9 were obtained and annotated. Comparative genomic analysis revealed that the predicted xylose metabolic pathway is conserved among the xylitol-producing yeasts Spathaspora sp. JA1, M. caribbica JA9 and Meyerozyma guilliermondii, but not in Spathaspora passalidarum, an efficient ethanol-producing yeast. Xylitol-producing yeasts showed strictly NADPH-dependent xylose reductase and NAD+-dependent xylitol-dehydrogenase activities. This imbalance of cofactors favors the high xylitol yield shown by Spathaspora sp. JA1, which is similar to the most efficient xylitol producers described so far.

Genomic overview of closely related fungi with different Protea host ranges

Genome comparisons of species with distinctive ecological traits can elucidate genetic divergence that influenced their differentiation. The interaction of a microorganism with its biotic environment is largely regulated by secreted compounds, and these can be predicted from genome sequences. In this study, we considered Knoxdaviesia capensis and Knoxdaviesia proteae, two closely related saprotrophic fungi found exclusively in Protea plants. We investigated their genome structure to compare their potential inter-specific interactions based on gene content. Their genomes displayed macrosynteny and were approximately 10 % repetitive. Both species had fewer secreted proteins than pathogens and other saprotrophs, reflecting their specialized habitat. The bulk of the predicted species-specific and secreted proteins coded for carbohydrate metabolism, with a slightly higher number of unique carbohydrate-degrading proteins in the broad host-range K. capensis. These fungi have few secondary metabolite gene clusters, suggesting minimal competition with other microbes and symbiosis with antibiotic-producing bacteria common in this niche. Secreted proteins associated with detoxification and iron sequestration likely enable these Knoxdaviesia species to tolerate antifungal compounds and compete for resources, facilitating their unusual dominance. This study confirms the genetic cohesion between Protea-associated Knoxdaviesia species and reveals aspects of their ecology that have likely evolved in response to their specialist niche.

Assembly of Schizosaccharomyces cryophilus chromosomes and their comparative genomic analyses revealed principles of genome evolution of the haploid fission yeasts

The fission yeast clade, which has a distinct life history from other yeasts, can provide important clues about evolutionary changes. To reveal these changes the large S. cryophilus supercontigs were assembled into chromosomes using synteny relationships and the conserved pericentromeric, subtelomeric genes. Togetherness of the supercontigs was confirmed by PCR. Investigation of the gene order revealed localisation of the rDNA arrays, more than 300 new conserved orthologues and proved that S. cryophilus supercontigs were mosaics of collinear blocks. PFGE analysis showed that size of the S. cryophilus chromosomes differ from the S. pombe chromosomes. Comparative genomic analyses of the newly assembled chromosomes confirmed that the closest relative of S. cryophilus was S. octosporus not just in sequence similarity but also in a structural way, and revealed that preservation of the conserved regions did not arise from the lower number of chromosomal rearrangements. Translocations were more typical in the closely related species, while the number of inversions increased with the phylogenetic distances. Our data suggested that sites of the chromosomal rearrangements were not random and often associated with repetitive sequences, structural- and nucleotide evolution might correlate. Chromosomal rearrangements of the fission yeasts compared to other lineages were also discussed.

Ceratocystidaceae exhibit high levels of recombination at the mating-type (MAT) locus

Mating is central to many fungal life cycles and is controlled by genes at the mating-type (MAT) locus. These genes determine whether the fungus will be self-sterile (heterothallic) or self-fertile (homothallic). Species in the ascomycete family Ceratocystidaceae have different mating strategies, making them interesting to consider with regards to their MAT loci. The aim of this study was to compare the composition of the MAT locus flanking regions in 11 species of Ceratocystidaceae representing four genera. Genome assemblies for each species were examined to identify the MAT locus and determine the structure of the flanking regions. Large contigs containing the MAT locus were then functionally annotated and analysed for the presence of transposable elements. Genes typically flanking the MAT locus in sordariomycetes were found to be highly conserved in the Ceratocystidaceae. The different genera in the Ceratocystidaceae displayed little synteny outside of the immediate MAT locus flanking genes. Even though species ofCeratocystis did not show much synteny outside of the immediate MAT locus flanking genes, species of Huntiella and Endoconidiophora were comparatively syntenic. Due to the high number of transposable elements present in Ceratocystis MAT flanking regions, we hypothesise that Ceratocystis species may have undergone recombination in this region.

Pan-Parastagonospora Comparative Genome Analysis-Effector Prediction and Genome Evolution

We report a fungal pan-genome study involving Parastagonospora spp., including 21 isolates of the wheat (Triticum aestivum) pathogen Parastagonospora nodorum, 10 of the grass-infecting Parastagonospora avenae, and 2 of a closely related undefined sister species. We observed substantial variation in the distribution of polymorphisms across the pan-genome, including repeat-induced point mutations, diversifying selection and gene gains and losses. We also discovered chromosome-scale inter and intraspecific presence/absence variation of some sequences, suggesting the occurrence of one or more accessory chromosomes or regions that may play a role in host-pathogen interactions. The presence of known pathogenicity effector loci SnToxA, SnTox1, and SnTox3 varied substantially among isolates. Three P. nodorum isolates lacked functional versions for all three loci, whereas three P. avenae isolates carried one or both of the SnTox1 and SnTox3 genes, indicating previously unrecognized potential for discovering additional effectors in the P. nodorum-wheat pathosystem. We utilized the pan-genomic comparative analysis to improve the prediction of pathogenicity effector candidates, recovering the three confirmed effectors among our top-ranked candidates. We propose applying this pan-genomic approach to identify the effector repertoire involved in other host-microbe interactions involving necrotrophic pathogens in the Pezizomycotina.

Genomic Analysis of the Insect-Killing Fungus Beauveria bassiana JEF-007 as a Biopesticide

Insect-killing fungi have high potential in pest management. A deeper insight into the fungal genes at the whole genome level is necessary to understand the inter-species or intra-species genetic diversity of fungal genes, and to select excellent isolates. In this work, we conducted a whole genome sequencing of Beauveria bassiana (Bb) JEF-007 and characterized pathogenesis-related features and compared with other isolates including Bb ARSEF2860. A large number of Bb JEF-007 genes showed high identity with Bb ARSEF2860, but some genes showed moderate or low identity. The two Bb isolates showed a significant difference in vegetative growth, antibiotic-susceptibility, and virulence against Tenebrio molitor larvae. When highly identical genes between the two Bb isolates were subjected to real-time PCR, their transcription levels were different, particularly in heat shock protein 30 (hsp30) gene which is related to conidial thermotolerance. In several B. bassiana isolates, chitinases and trypsin-like protease genes involved in pathogenesis were highly conserved, but other genes showed noticeable sequence variation within the same species. Given the transcriptional and genetic diversity in B. bassiana, a selection of virulent isolates with industrial advantages is a pre-requisite, and this genetic approach could support the development of excellent biopesticides with intellectual property protection.

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

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

Genomic evidence for intraspecific hybridization in a clonal and extremely halotolerant yeast

Background

The black yeast Hortaea werneckii (Dothideomycetes, Ascomycota) is one of the most extremely halotolerant fungi, capable of growth at NaCl concentrations close to saturation. Although dothideomycetous fungi are typically haploid, the reference H. werneckii strain has a diploid genome consisting of two subgenomes with a high level of heterozygosity.

Results

In order to explain the origin of the H. werneckii diploid genome we here report the genome sequencing of eleven strains isolated from different habitats and geographic locations. Comparison of nine diploid and two haploid strains showed that the reference genome was likely formed by hybridization between two haploids and not by endoreduplication as suggested previously. Results also support additional hybridization events in the evolutionary history of investigated strains, however exchange of genetic material in the species otherwise appears to be rare. Possible links between such unusual reproduction and the extremotolerance of H. werneckii remain to be investigated.

Conclusions

H. werneckii appears to be able to form persistent haploid as well as diploid strains, is capable of occasional hybridization between relatively heterozygous haploids, but is otherwise limited to clonal reproduction. The reported data and the first identification of haploid H. werneckii strains establish this species as a good model for studying the effects of ploidy and hybridization in an extremotolerant system unperturbed by frequent genetic recombination.

Electronic supplementary material

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

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

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

The Fungal Tree of Life: from Molecular Systematics to Genome-Scale Phylogenies

The kingdom Fungi is one of the more diverse clades of eukaryotes in terrestrial ecosystems, where they provide numerous ecological services ranging from decomposition of organic matter and nutrient cycling to beneficial and antagonistic associations with plants and animals. The evolutionary relationships of the kingdom have represented some of the more recalcitrant problems in systematics and phylogenetics. The advent of molecular phylogenetics, and more recently phylogenomics, has greatly advanced our understanding of the patterns and processes associated with fungal evolution, however. In this article, we review the major phyla, subphyla, and classes of the kingdom Fungi and provide brief summaries of ecologies, morphologies, and exemplar taxa. We also provide examples of how molecular phylogenetics and evolutionary genomics have advanced our understanding of fungal evolution within each of the phyla and some of the major classes. In the current classification we recognize 8 phyla, 12 subphyla, and 46 classes within the kingdom. The ancestor of fungi is inferred to be zoosporic, and zoosporic fungi comprise three lineages that are paraphyletic to the remainder of fungi. Fungi historically classified as zygomycetes do not form a monophyletic group and are paraphyletic to Ascomycota and Basidiomycota. Ascomycota and Basidiomycota are each monophyletic and collectively form the subkingdom Dikarya.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

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

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

Specialized plant biochemistry drives gene clustering in fungi

The fitness and evolution of prokaryotes and eukaryotes are affected by the organization of their genomes. In particular, the physical clustering of genes can coordinate gene expression and can prevent the breakup of co-adapted alleles. Although clustering may thus result from selection for phenotype optimization and persistence, the impact of environmental selection pressures on eukaryotic genome organization has rarely been systematically explored. Here, we investigated the organization of fungal genes involved in the degradation of phenylpropanoids, a class of plant-produced secondary metabolites that mediate many ecological interactions between plants and fungi. Using a novel gene cluster detection method, we identified 1110 gene clusters and many conserved combinations of clusters in a diverse set of fungi. We demonstrate that congruence in genome organization over small spatial scales is often associated with similarities in ecological lifestyle. Additionally, we find that while clusters are often structured as independent modules with little overlap in content, certain gene families merge multiple modules into a common network, suggesting they are important components of phenylpropanoid degradation strategies. Together, our results suggest that phenylpropanoids have repeatedly selected for gene clustering in fungi, and highlight the interplay between genome organization and ecological evolution in this ancient eukaryotic lineage.

Whole-genome assembly of the coral reef Pearlscale Pygmy Angelfish (Centropyge vrolikii)

The diversity of DNA sequencing methods and algorithms for genome assemblies presents scientists with a bewildering array of choices. Here, we construct and compare eight candidate assemblies combining overlapping shotgun read data, mate-pair and Chicago libraries and four different genome assemblers to produce a high-quality draft genome of the iconic coral reef Pearlscale Pygmy Angelfish, Centropyge vrolikii (family Pomacanthidae). The best candidate assembly combined all four data types and had a scaffold N50 127.5 times higher than the candidate assembly obtained from shotgun data only. Our best candidate assembly had a scaffold N50 of 8.97 Mb, contig N50 of 189,827, and 97.4% complete for BUSCO v2 (Actinopterygii set) and 95.6% complete for CEGMA matches. These contiguity and accuracy scores are higher than those of any other fish assembly released to date that did not apply linkage map information, including those based on more expensive long-read sequencing data. Our analysis of how different data types improve assembly quality will help others choose the most appropriate de novo genome sequencing strategy based on resources and target applications. Furthermore, the draft genome of the Pearlscale Pygmy angelfish will play an important role in future studies of coral reef fish evolution, diversity and conservation.

Prospects for Fungal Bioremediation of Acidic Radioactive Waste Sites: Characterization and Genome Sequence of Rhodotorula taiwanensis MD1149

Highly concentrated radionuclide waste produced during the Cold War era is stored at US Department of Energy (DOE) production sites. This radioactive waste was often highly acidic and mixed with heavy metals, and has been leaking into the environment since the 1950s. Because of the danger and expense of cleanup of such radioactive sites by physicochemical processes, in situ bioremediation methods are being developed for cleanup of contaminated ground and groundwater. To date, the most developed microbial treatment proposed for high-level radioactive sites employs the radiation-resistant bacterium Deinococcus radiodurans. However, the use of Deinococcus spp. and other bacteria is limited by their sensitivity to low pH. We report the characterization of 27 diverse environmental yeasts for their resistance to ionizing radiation (chronic and acute), heavy metals, pH minima, temperature maxima and optima, and their ability to form biofilms. Remarkably, many yeasts are extremely resistant to ionizing radiation and heavy metals. They also excrete carboxylic acids and are exceptionally tolerant to low pH. A special focus is placed on Rhodotorula taiwanensis MD1149, which was the most resistant to acid and gamma radiation. MD1149 is capable of growing under 66 Gy/h at pH 2.3 and in the presence of high concentrations of mercury and chromium compounds, and forming biofilms under high-level chronic radiation and low pH. We present the whole genome sequence and annotation of R. taiwanensis strain MD1149, with a comparison to other Rhodotorula species. This survey elevates yeasts to the frontier of biology's most radiation-resistant representatives, presenting a strong rationale for a role of fungi in bioremediation of acidic radioactive waste sites.

Evolution within the fungal genus Verticillium is characterized by chromosomal rearrangement and gene loss

The fungal genus Verticillium contains ten species, some of which are notorious plant pathogens causing vascular wilt diseases in host plants, while others are known as saprophytes and opportunistic plant pathogens. Whereas the genome of V. dahliae, the most notorious plant pathogen of the genus, has been well characterized, evolution and speciation of other members of the genus received little attention thus far. Here, we sequenced the genomes of the nine haploid Verticillium spp. to study evolutionary trajectories of their divergence from a last common ancestor. Frequent occurrence of chromosomal rearrangement and gene family loss was identified. In addition to ∼11 000 genes that are shared at least between two species, only 200-600 species-specific genes occur. Intriguingly, these species-specific genes show different features than the shared genes.