IMA Genome-F 6: Draft genome sequences of Armillaria fuscipes, Ceratocystiopsis minuta, Ceratocystis adiposa, Endoconidiophora laricicola, E. polonica and Penicillium freii DAOMC 242723

Unidirectional mating-type switching is underpinned by a conserved MAT1 locus architecture

Unidirectional mating-type switching is a form of homothallic reproduction known only in a small number of filamentous ascomycetes. Their ascospores can give rise to either self-sterile isolates that require compatible partners for subsequent sexual reproduction, or self-fertile individuals capable of completing this process in isolation. The limited studies previously conducted in these fungi suggest that the differences in mating specificity are determined by the architecture of the MAT1 locus. In self-fertile isolates that have not undergone unidirectional mating-type switching, the locus contains both MAT1-1 and MAT1-2 mating-type genes, typical of primary homothallism. In the self-sterile isolates produced after a switching event, the MAT1-2 genes are lacking from the locus, likely due to a recombination-mediated deletion of the MAT1-2 gene information. To determine whether these arrangements of the MAT1 locus support unidirectional mating-type switching in the Ceratocystidaceae, the largest known fungal assemblage capable of this reproduction strategy, a combination of genetic and genomic approaches were used. The MAT1 locus was annotated in representative species of Ceratocystis, Endoconidiophora, and Davidsoniella. In all cases, MAT1-2 genes interrupted the MAT1-1-1 gene in self-fertile isolates. The MAT1-2 genes were flanked by two copies of a direct repeat that accurately predicted the boundaries of the deletion event that would yield the MAT1 locus of self-sterile isolates. Although the relative position of the MAT1-2 gene region differed among species, it always disrupted the MAT1-1-1 gene and/or its expression in the self-fertile MAT1 locus. Following switching, this gene and/or its expression was restored in the self-sterile arrangement of the locus. This mirrors what has been reported in other species capable of unidirectional mating-type switching, providing the strongest support for a conserved MAT1 locus structure that is associated with this process. This study contributes to our understanding of the evolution of unidirectional mating-type switching.

A practical approach to genome assembly and annotation of Basidiomycota using the example of Armillaria

Technological advancements in genome sequencing, assembly and annotation platforms and algorithms that resulted in several genomic studies have created an opportunity to further our understanding of the biology of phytopathogens, including Armillaria species. Most Armillaria species are facultative necrotrophs that cause root- and stem-rot, usually on woody plants, significantly impacting agriculture and forestry worldwide. Genome sequencing, assembly and annotation in terms of samples used and methods applied in Armillaria genome projects are evaluated in this review. Infographic guidelines and a database of resources to facilitate future Armillaria genome projects were developed. Knowledge gained from genomic studies of Armillaria species is summarized and prospects for further research are provided. This guide can be applied to other diploid and dikaryotic fungal genomics.

Hybrid Genome Assembly of Berkeleyomyces rouxiae, an Emerging Cannabis Fungal Pathogen Causing Black Root Rot in an Aeroponic Facility

The resurged interest in cultivation of Cannabis sativa has presented an array of new challenges. Among them are the difficult-to-control pests and pathogens that infect cannabis plants. The limited methods for disease control available to cannabis growers necessitates early detection of plant pathogens, something that molecular techniques such as DNA sequencing has greatly improved. This study reports for the first time the fungal plant pathogen Berkeleyomyces rouxiae causing black root rot in high THC-containing cannabis. Aeroponically grown cannabis plants at a licenced production facility in Cranbrook BC, Canada, rapidly displayed root discoloration and rot symptoms despite testing negative for all commercially available pathogen tests. Developing sequencing-based disease diagnostics requires genomic information, so this study presents the first whole genome sequence of the multihost, widespread black root rot pathogen B. rouxiae. Hybrid genome assembly using Oxford Nanopore long-reads and Illumina short-reads yielded a genome size of 28.2 Mb represented over 404 contigs with an N50 of 267 kb. Genome annotation predicted 6,960 protein-coding genes with 59,477 functional annotations. The availability of this genome will assist in sequence-based diagnostic development, comparative genomics, and taxonomic resolution of this globally important plant pathogen.

Whole genome sequencing and analysis of Armillaria gallica Jzi34 symbiotic with Gastrodia elata

BACKGROUND

Armillaria species are plant pathogens, but a few Armillaria species can establish a symbiotic relationship with Gastrodia elata, a rootless and leafless orchid, that is used as a Chinese herbal medicine. Armillaria is a nutrient source for the growth of G. elata. However, there are few reports on the molecular mechanism of symbiosis between Armillaria species and G. elata. The genome sequencing and analysis of Armillaria symbiotic with G. elata would provide genomic information for further studying the molecular mechanism of symbiosis.

RESULTS

The de novo genome assembly was performed with the PacBio Sequel platform and Illumina NovaSeq PE150 for the A. gallica Jzi34 strain, which was symbiotic with G. elata. Its genome assembly contained ~ 79.9 Mbp and consisted of 60 contigs with an N50 of 2,535,910 bp. There were only 4.1% repetitive sequences in the genome assembly. Functional annotation analysis revealed a total of 16,280 protein coding genes. Compared with the other five genomes of Armillaria, the carbohydrate enzyme gene family of the genome was significantly contracted, while it had the largest set of glycosyl transferase (GT) genes. It also had an expansion of auxiliary activity enzymes AA3-2 gene subfamily and cytochrome P450 genes. The synteny analysis result of P450 genes reveals that the evolutionary relationship of P450 proteins between A. gallica Jzi34 and other four Armillaria was complex.

CONCLUSIONS

These characteristics may be beneficial for establishing a symbiotic relationship with G. elata. These results explore the characteristics of A. gallica Jzi34 from a genomic perspective and provide an important genomic resource for further detailed study of Armillaria. This will help to further study the symbiotic mechanism between A. gallica and G. elata.

An in Silico Approach to Identifying TF Binding Sites: Analysis of the Regulatory Regions of BUSCO Genes from Fungal Species in the Ceratocystidaceae Family

Transcriptional regulation controls gene expression through regulatory promoter regions that contain conserved sequence motifs. These motifs, also known as regulatory elements, are critically important to expression, which is driving research efforts to identify and characterize them. Yeasts have been the focus of such studies in fungi, including in several in silico approaches. This study aimed to determine whether in silico approaches could be used to identify motifs in the Ceratocystidaceae family, and if present, to evaluate whether these correspond to known transcription factors. This study targeted the 1000 base-pair region upstream of the start codon of 20 single-copy genes from the BUSCO dataset for motif discovery. Using the MEME and Tomtom analysis tools, conserved motifs at the family level were identified. The results show that such in silico approaches could identify known regulatory motifs in the Ceratocystidaceae and other unrelated species. This study provides support to ongoing efforts to use in silico analyses for motif discovery.

The compact mitogenome of Ceratocystiopsis pallidobrunnea

Ceratocystiopsis is a fungal genus that has been assigned to the Ophiostomatales, fungi known for their association with various bark beetles and other arthropods. The mitochondrial genome of Ceratocystiopsis pallidobrunnea has been characterized and compared with other members of the genus Ceratocystiopsis and Ophiostomatales. At 29 022 bp, the mitogenome of C. pallidobrunnea is the smallest reported so far for this genus. Gene arrangement was observed to be conserved for this group of fungi and mitogenome variation appears to be mostly due to the absence and presence of introns. The long-term goal is to apply mitogenomes to resolve taxonomic issues within the Ophiostomatales and within the various genera that comprise the Ophiostomataceae.

Phenolic degradation by catechol dioxygenases is associated with pathogenic fungi with a necrotrophic lifestyle in the Ceratocystidaceae

Fungal species of the Ceratocystidaceae grow on their host plants using a variety of different lifestyles, from saprophytic to highly pathogenic. Although many genomes of fungi in the Ceratocystidaceae are publicly available, it is not known how the genes that encode catechol dioxygenases (CDOs), enzymes involved in the degradation of phenolic plant defense compounds, differ among members of the Ceratocystidaceae. The aim of this study was therefore to identify and characterize the genes encoding CDOs in the genomes of Ceratocystidaceae representatives. We found that genes encoding CDOs are more abundant in pathogenic necrotrophic species of the Ceratocystidaceae and less abundant in saprophytic species. The loss of the CDO genes and the associated 3-oxoadipate catabolic pathway appears to have occurred in a lineage-specific manner. Taken together, this study revealed a positive association between CDO gene copy number and fungal lifestyle in Ceratocystidaceae representatives.

Phylogenetic and phylogenomic analyses reveal two new genera and three new species of ophiostomatalean fungi from termite fungus combs

The Ophiostomatales (Ascomycota) accommodates more than 300 species characterized by similar morphological adaptations to arthropod dispersal. Most species in this order are wood-inhabiting fungi associated with bark or ambrosia beetles. However, a smaller group of species occur in other niches such as in soil and Protea infructescences. Recent surveys of Termitomyces fungus gardens (fungus combs) of fungus-growing termites led to the discovery of characteristic ophiostomatalean-like fruiting structures. In this study, these ophiostomatalean-like fungi were identified using morphological characteristics, conventional molecular markers, and whole genome sequencing. In addition, the influence of the extracts derived from various parts of Termitomyces combs on the growth of these fungi in culture was considered. Based on phylogenomic analyses, two new genera (Intubia and Chrysosphaeria) were introduced to accommodate these ophiostomatalean species. Phylogenetic analyses revealed that the isolates resided in three well-supported lineages, and these were described as three new species (Intubia macrotermitinarum, I. oerlemansii, and Chrysosphaeria jan-nelii). Culture-based studies showed that these species do not depend on the Termitomyces comb material for growth.

Ophiostomatoid fungi associated with pine bark beetles and infested pines in south-eastern Australia, including Graphilbum ipis-grandicollis sp. nov

The ophiostomatoid fungi are an assemblage of ascomycetes which are arguably best-known for their associations with bark and ambrosia beetles (Curculonidae) and blue stain (sap stain) of many economically important tree species. These fungi are considered a significant threat to coniferous forests, which has resulted in numerous studies characterising the diversity of bark beetles and their ophiostomatoid associates globally. The diversity of ophiostomatoid fungi present in Australian pine plantations, however, remains largely undetermined. The aims of this study were therefore to reconsider the diversity of ophiostomatoid fungi associated with Pinus in Australia, and to establish the baseline of expected taxa found within these plantation ecosystems. To achieve this, we reviewed Australian plant pathogen reference collections, and analysed samples collected during forest health surveillance programs from the major pine growing regions in south-eastern Australia. In total, 135 ophiostomatoid isolates (15 from reference collections and 120 collected during the current study) were assessed using morphological identification and ITS screening which putatively distinguished 15 taxonomic groups. Whole genome sequencing (WGS) of representative isolates from each taxon was performed to obtain high-quality sequence data for multi-locus phylogenetic analysis. Our results revealed a greater than expected diversity, expanding the status of ophiostomatoid fungi associated with Pinus in Australia to include 14 species from six genera in the Ophiostomatales and a single species residing in the Microascales. While most of these were already known to science, our study includes seven first records for Australia and the description of one new species, Graphilbum ipis-grandicollis sp. nov.. This study also provides an early example of whole genome sequencing (WGS) approaches replacing traditional PCR-based methods for taxonomic surveys. This not only allowed for robust multi-locus sequence extraction during taxonomic assessment, but also permitted the rapid establishment of a curated genomic database for ophiostomatoid fungi which will continue to aid in the development of improved diagnostic resources and capabilities for Australian biosecurity.

Characterization of the Ergosterol Biosynthesis Pathway in Ceratocystidaceae

Terpenes represent the biggest group of natural compounds on earth. This large class of organic hydrocarbons is distributed among all cellular organisms, including fungi. The different classes of terpenes produced by fungi are mono, sesqui, di- and triterpenes, although triterpene ergosterol is the main sterol identified in cell membranes of these organisms. The availability of genomic data from members in the Ceratocystidaceae enabled the detection and characterization of the genes encoding the enzymes in the mevalonate and ergosterol biosynthetic pathways. Using a bioinformatics approach, fungal orthologs of sterol biosynthesis genes in nine different species of the Ceratocystidaceae were identified. Ergosterol and some of the intermediates in the pathway were also detected in seven species (Ceratocystis manginecans, C. adiposa, Huntiella moniliformis, Thielaviopsis punctulata, Bretziella fagacearum, Endoconidiophora polonica and Davidsoniella virescens), using gas chromatography-mass spectrometry analysis. The average ergosterol content differed among different genera of Ceratocystidaceae. We also identified all possible terpene related genes and possible biosynthetic clusters in the genomes used in this study. We found a highly conserved terpene biosynthesis gene cluster containing some genes encoding ergosterol biosynthesis enzymes in the analysed genomes. An additional possible terpene gene cluster was also identified in all of the Ceratocystidaceae. We also evaluated the sensitivity of the Ceratocystidaceae to a triazole fungicide that inhibits ergosterol synthesis. The results showed that different members of this family behave differently when exposed to different concentrations of triazole tebuconazole.

The Mitogenomes of Ophiostoma minus and Ophiostoma piliferum and Comparisons With Other Members of the Ophiostomatales

Fungi assigned to the Ophiostomatales are of economic concern as many are blue-stain fungi and some are plant pathogens. The mitogenomes of two blue-stain fungi, Ophiostoma minus and Ophiostoma piliferum, were sequenced and compared with currently available mitogenomes for other members of the Ophiostomatales. Species representing various genera within the Ophiostomatales have been examined for gene content, gene order, phylogenetic relationships, and the distribution of mobile elements. Gene synteny is conserved among the Ophiostomatales but some members were missing the atp9 gene. A genome wide intron landscape has been prepared to demonstrate the distribution of the mobile genetic elements (group I and II introns and homing endonucleases) and to provide insight into the evolutionary dynamics of introns among members of this group of fungi. Examples of complex introns or nested introns composed of two or three intron modules have been observed in some species. The size variation among the mitogenomes (from 23.7 kb to about 150 kb) is mostly due to the presence and absence of introns. Members of the genus Sporothrix sensu stricto appear to have the smallest mitogenomes due to loss of introns. The taxonomy of the Ophiostomatales has recently undergone considerable revisions; however, some lineages remain unresolved. The data showed that genera such as Raffaelea appear to be polyphyletic and the separation of Sporothrix sensu stricto from Ophiostoma is justified.

Epidemiology, Biotic Interactions and Biological Control of Armillarioids in the Northern Hemisphere

Armillarioids, including the genera Armillaria, Desarmillaria and Guyanagaster, represent white-rot specific fungal saprotrophs with soilborne pathogenic potentials on woody hosts. They propagate in the soil by root-like rhizomorphs, connecting between susceptible root sections of their hosts, and often forming extended colonies in native forests. Pathogenic abilities of Armillaria and Desarmillaria genets can readily manifest in compromised hosts, or hosts with full vigour can be invaded by virulent mycelia when exposed to a larger number of newly formed genets. Armillaria root rot-related symptoms are indicators of ecological imbalances in native forests and plantations at the rhizosphere levels, often related to abiotic environmental threats, and most likely unfavourable changes in the microbiome compositions in the interactive zone of the roots. The less-studied biotic impacts that contribute to armillarioid host infection include fungi and insects, as well as forest conditions. On the other hand, negative biotic impactors, like bacterial communities, antagonistic fungi, nematodes and plant-derived substances may find applications in the environment-friendly, biological control of armillarioid root diseases, which can be used instead of, or in combination with the classical, but frequently problematic silvicultural and chemical control measures.

Draft genomic sequence of Armillaria gallica 012m: insights into its symbiotic relationship with Gastrodia elata

Armillaria species (Basidiomycota, Physalacriaceae) are well known as plant pathogens related to serious root rot disease on various trees in forests and plantations. Interestingly, some Armillaria species are essential symbionts of the rare Chinese medicinal herb Gastrodia elata, a rootless and leafless orchid used for over 2000 years. In this work, an 87.3-M draft genome of Armillaria gallica 012m strain, which was symbiotic with G. elata, was assembled. The genome includes approximately 23.6% repetitive sequences and encodes 26,261 predicted genes. In comparison with other four genomes of Armillaria, the following gene families related to pathogenicity/saprophytic phase, including cytochrome P450 monooxygenases, carbohydrate-active enzyme AA3, and hydrophobins, were significantly contracted in A. gallica 012m. These characteristics may be beneficial for G. elata to get less injuries. The genome-guided analysis of differential expression between rhizomorph (RH) and vegetative mycelium (VM) showed that a total of 2549 genes were differentially expressed, including 632 downregulated genes and 1917 upregulated genes. In the RH, most differentially expressed genes (DEGs) related to pathogenicity were significantly upregulated. To further elucidate gene function, Gene Ontology enrichment analysis showed that the upregulated DEGs significantly grouped into monooxygenase activity, hydrolase activity, glucosidase activity, extracellular region, fungal cell wall, response to xenobiotic stimulus, response to toxic substance, etc. These phenomena indicate that RH had better infection ability than VM. The infection ability of RH may be beneficial for G. elata to obtain nutrition, because the rhizomorph constantly infected the nutritional stems of G. elata and formed the hyphae that can be digested by G. elata. These results clarified the characteristics of A. gallica 012m and the reason why the strain 012m can establish a symbiotic relationship with G. elata in some extent from the perspective of genomics.

Phylogenomic incongruence in Ceratocystis: a clue to speciation?

Background

The taxonomic history of Ceratocystis, a genus in the Ceratocystidaceae, has been beset with questions and debate. This is due to many of the commonly used species recognition concepts (e.g., morphological and biological species concepts) providing different bases for interpretation of taxonomic boundaries. Species delineation in Ceratocystis primarily relied on genealogical concordance phylogenetic species recognition (GCPSR) using multiple standard molecular markers.

Results

Questions have arisen regarding the utility of these markers e.g., ITS, BT and TEF1-α due to evidence of intragenomic variation in the ITS, as well as genealogical incongruence, especially for isolates residing in a group referred to as the Latin-American clade (LAC) of the species. This study applied a phylogenomics approach to investigate the extent of phylogenetic incongruence in Ceratocystis. Phylogenomic analyses of a total of 1121 shared BUSCO genes revealed widespread incongruence within Ceratocystis, particularly within the LAC, which was typified by three equally represented topologies. Comparative analyses of the individual gene trees revealed evolutionary patterns indicative of hybridization. The maximum likelihood phylogenetic tree generated from the concatenated dataset comprised of 1069 shared BUSCO genes provided improved phylogenetic resolution suggesting the need for multiple gene markers in the phylogeny of Ceratocystis.

Conclusion

The incongruence observed among single gene phylogenies in this study call into question the utility of single or a few molecular markers for species delineation. Although this study provides evidence of interspecific hybridization, the role of hybridization as the source of discordance will require further research because the results could also be explained by high levels of shared ancestral polymorphism in this recently diverged lineage. This study also highlights the utility of BUSCO genes as a set of multiple orthologous genes for phylogenomic studies.

Chromosomal assembly and analyses of genome-wide recombination rates in the forest pathogenic fungus Armillaria ostoyae

Recombination shapes the evolutionary trajectory of populations and plays an important role in the faithful transmission of chromosomes during meiosis. Levels of sexual reproduction and recombination are important properties of host-pathogen interactions because the speed of antagonistic co-evolution depends on the ability of hosts and pathogens to generate genetic variation. However, our understanding of the importance of recombination is limited because large taxonomic groups remain poorly investigated. Here, we analyze recombination rate variation in the basidiomycete fungus Armillaria ostoyae, which is an aggressive pathogen on a broad range of conifers and other trees. We analyzed a previously constructed, dense genetic map based on 198 single basidiospore progeny from a cross. Progeny were genotyped at a genome-wide set of single-nucleotide polymorphism (SNP) markers using double digest restriction site associated DNA sequencing. Based on a linkage map of on 11,700 SNPs spanning 1007.5 cM, we assembled genomic scaffolds into 11 putative chromosomes of a total genome size of 56.6 Mb. We identified 1984 crossover events among all progeny and found that recombination rates were highly variable along chromosomes. Recombination hotspots tended to be in regions close to the telomeres and were more gene-poor than the genomic background. Genes in proximity to recombination hotspots were encoding on average shorter proteins and were enriched for pectin degrading enzymes. Our analyses enable more powerful population and genome-scale studies of a major tree pathogen.

IMA Genome-F 11: Draft genome sequences of Fusarium xylarioides, Teratosphaeria gauchensis and T. zuluensis and genome annotation for Ceratocystis fimbriata

Draft genomes of the fungal species Fusarium xylarioides, Teratosphaeria gauchensis and T. zuluensis are presented. In addition an annotation of the genome of Ceratocystis fimbriata is presented. Overall these genomes provide a valuable resource for understanding the molecular processes underlying pathogenicity and potential management strategies of these economically important fungi.

Patterns of coevolution between ambrosia beetle mycangia and the Ceratocystidaceae, with five new fungal genera and seven new species

Ambrosia beetles farm specialised fungi in sapwood tunnels and use pocket-like organs called mycangia to carry propagules of the fungal cultivars. Ambrosia fungi selectively grow in mycangia, which is central to the symbiosis, but the history of coevolution between fungal cultivars and mycangia is poorly understood. The fungal family Ceratocystidaceae previously included three ambrosial genera (Ambrosiella, Meredithiella, and Phialophoropsis), each farmed by one of three distantly related tribes of ambrosia beetles with unique and relatively large mycangium types. Studies on the phylogenetic relationships and evolutionary histories of these three genera were expanded with the previously unstudied ambrosia fungi associated with a fourth mycangium type, that of the tribe Scolytoplatypodini. Using ITS rDNA barcoding and a concatenated dataset of six loci (28S rDNA, 18S rDNA, tef1-α, tub, mcm7, and rpl1), a comprehensive phylogeny of the family Ceratocystidaceae was developed, including Inodoromyces interjectus gen. & sp. nov., a non-ambrosial species that is closely related to the family. Three minor morphological variants of the pronotal disk mycangium of the Scolytoplatypodini were associated with ambrosia fungi in three respective clades of Ceratocystidaceae: Wolfgangiella gen. nov., Toshionella gen. nov., and Ambrosiella remansi sp. nov. Closely-related species that are not symbionts of ambrosia beetles are accommodated by Catunica adiposa gen. & comb. nov. and Solaloca norvegica gen. & comb. nov. The divergent morphology of the ambrosial genera and their phylogenetic placement among non-ambrosial genera suggest three domestication events in the Ceratocystidaceae. Estimated divergence dates for the ambrosia fungi and mycangia suggest that Scolytoplatypodini mycangia may have been the first to acquire Ceratocystidaceae symbionts and other ambrosial fungal genera emerged shortly after the evolution of new mycangium types. There is no evidence of reversion to a non-ambrosial lifestyle in the mycangial symbionts.

Genomic analysis of the aggressive tree pathogen Ceratocystis albifundus

The overall goal of this study was to determine whether the genome of an important plant pathogen in Africa, Ceratocystis albifundus, is structured into subgenomic compartments, and if so, to establish how these compartments are distributed across the genome. For this purpose, the publicly available genome of C. albifundus was complemented with the genome sequences for four additional isolates using the Illumina HiSeq platform. In addition, a reference genome for one of the individuals was assembled using both PacBio and Illumina HiSeq technologies. Our results showed a high degree of synteny between the five genomes, although several regions lacked detectable long-range synteny. These regions were associated with the presence of accessory genes, lower genetic similarity, variation in read-map depth, as well as transposable elements and genes associated with host-pathogen interactions (e.g. effectors and CAZymes). Such patterns are regarded as hallmarks of accelerated evolution, particularly of accessory subgenomic compartments in fungal pathogens. Our findings thus showed that the genome of C. albifundus is made-up of core and accessory subgenomic compartments, which is an important step towards characterizing its pangenome. This study also highlights the value of comparative genomics for understanding mechanisms that may underly and influence the biology and evolution of pathogens.

The mitochondrial genome of Endoconidiophora resinifera is intron rich

Endoconidiophora resinifera (=Ceratocystis resinifera) is a blue-stain fungus that occurs on conifers. The data showed that the Endoconidiophora resinifera mitochondrial genome is one of the largest mitochondrial genomes (>220 kb) so far reported among members of the Ascomycota. An exceptional large number of introns (81) were noted and differences among the four strains were restricted to minor variations in intron numbers and a few indels and single nucleotide polymorphisms. The major differences among the four strains examined are due to size polymorphisms generated by the absence or presence of mitochondrial introns. Also, these mitochondrial genomes encode the largest cytochrome oxidase subunit 1 gene (47.5 kb) reported so far among the fungi. The large size for this gene again can be attributed to the large number of intron insertions. This study reports the first mitochondrial genome for the genus Endoconidiophora, previously members of this genus were assigned to Ceratocystis. The latter genus has recently undergone extensive taxonomic revisions and the mitochondrial genome might provide loci that could be applied as molecular markers assisting in the identification of taxa within this group of economically important fungi. The large mitochondrial genome also may provide some insight on mechanisms that can lead to mitochondrial genome expansion.

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.

Heterothallism revealed in the root rot fungi Berkeleyomyces basicola and B. rouxiae

Berkeleyomyces basicola and Berkeleyomyces rouxiae, two sister species previously treated collectively as Thielaviopsis basicola, reside in the Ceratocystidaceae (Microascales, Ascomycota). Both species are important root pathogens of many important agricultural crops and ornamental plants. Although T. basicola has been known for more than 150y, a sexual state has never been found and it has been assumed to be an asexual pathogen. The aim of this study was to determine the mating strategy of the two Berkeleyomyces species. Investigation of the genome sequences of two B. basicola isolates allowed for the complete characterization of the MATlocus, revealing that it has a typical heterothallic mating system with the MAT1-1andMAT1-2 idiomorphs occurring in different isolates. PCR amplification using mating type primers developed in this study, showed that the MAT1-1-1andMAT1-2-1 genes were also present in different isolates of B. rouxiae. Pairing of isolates representing the two mating types of both species,using a variety of techniques failed to produce sexual structures. Although we have found no direct evidence that they reproduce sexually, these fungi are clearly heterothallic with both mating types occurring in some countries suggesting that a cryptic sexual cycle could exist for them.

IMA Genome-F 9: Draft genome sequence of Annulohypoxylon stygium, Aspergillus mulundensis, Berkeleyomyces basicola (syn. Thielaviopsis basicola), Ceratocystis smalleyi, two Cercospora beticola strains, Coleophoma cylindrospora, Fusarium fracticaudum, Phialophora cf. hyalina, and Morchella septimelata

Draft genomes of the species Annulohypoxylon stygium, Aspergillus mulundensis, Berkeleyomyces basicola (syn. Thielaviopsis basicola), Ceratocystis smalleyi, two Cercospora beticola strains, Coleophoma cylindrospora, Fusarium fracticaudum, Phialophora cf. hyalina and Morchella septimelata are presented. Both mating types (MAT1-1 and MAT1-2) of Cercospora beticola are included. Two strains of Coleophoma cylindrospora that produce sulfated homotyrosine echinocandin variants, FR209602, FR220897 and FR220899 are presented. The sequencing of Aspergillus mulundensis, Coleophoma cylindrospora and Phialophora cf. hyalina has enabled mapping of the gene clusters encoding the chemical diversity from the echinocandin pathways, providing data that reveals the complexity of secondary metabolism in these different species. Overall these genomes provide a valuable resource for understanding the molecular processes underlying pathogenicity (in some cases), biology and toxin production of these economically important fungi.

Diversity and evolution of polyketide biosynthesis gene clusters in the Ceratocystidaceae

Polyketides are secondary metabolites with diverse biological activities. Polyketide synthases (PKS) are often encoded from genes clustered in the same genomic region. Functional analyses and genomic studies show that most fungi are capable of producing a repertoire of polyketides. We considered the potential of Ceratocystidaceae for producing polyketides using a comparative genomics approach. Our aims were to identify the putative polyketide biosynthesis gene clusters, to characterize them and predict the types of polyketide compounds they might produce. We used sequences from nineteen species in the genera, Ceratocystis, Endoconidiophora, Davidsoniella, Huntiella, Thielaviopsis and Bretziella, to identify and characterize PKS gene clusters, by employing a range of bioinformatics and phylogenetic tools. We showed that the genomes contained putative clusters containing a non-reducing type I PKS and a type III PKS. Phylogenetic analyses suggested that these genes were already present in the ancestor of the Ceratocystidaceae. By contrast, the various reducing type I PKS-containing clusters identified in these genomes appeared to have distinct evolutionary origins. Although one of the identified clusters potentially allows for the production of melanin, their functional characterization will undoubtedly reveal many novel and important compounds implicated in the biology of the Ceratocystidaceae.

Unexpected placement of the MAT1-1-2 gene in the MAT1-2 idiomorph of Thielaviopsis

Sexual reproduction in the Ascomycota is controlled by genes encoded at the mating-type or MAT1 locus. The two allelic versions of this locus in heterothallic species, referred to as idiomorphs, are defined by the MAT1-1-1 (for the MAT1-1 idiomorph) and MAT1-2-1 (for the MAT1-2 idiomorph) genes. Both idiomorphs can contain additional genes, although the contents of each is typically specific to and conserved within particular Pezizomycotina lineages. Using full genome sequences, complemented with conventional PCR and Sanger sequencing, we compared the mating-type idiomorphs in heterothallic species of Thielaviopsis (Ceratocystidaceae). The analyses showed that the MAT1-1 idiomorph of T. punctulata, T. paradoxa, T. euricoi, T. ethacetica and T. musarum harboured only the expected MAT1-1-1 gene. In contrast, the MAT1-2 idiomorph of T. punctulata, T. paradoxa and T. euricoi encoded the MAT1-2-1, MAT1-2-7 and MAT1-1-2 genes. Of these, MAT1-2-1 and MAT1-2-7 are genes previously reported in this idiomorph, while MAT1-1-2 is known only in the MAT1-1 idiomorph. Phylogenetic analysis showed that the Thielaviopsis MAT1-1-2 groups with the known homologues of this gene in other Microascales, thus confirming its annotation. Previous work suggests that MAT1-1-2 is involved in fruiting body development, a role that would be unaffected by its idiomorphic position. This notion is supported by our findings for the MAT1 locus structure in Thielaviopsis species. This also serves as the first example of a MAT1-1-specific gene restricted to only the MAT1-2 idiomorph.

IMA Genome-F 8: Draft genome of Cercospora zeina, Fusarium pininemorale, Hawksworthiomyces lignivorus, Huntiella decipiens and Ophiostoma ips

The genomes of Cercospora zeina, Fusarium pininemorale, Hawksworthiomyces lignivorus, Huntiella decipiens, and Ophiostoma ips are presented in this genome announcement. Three of these genomes are from plant pathogens and otherwise economically important fungal species. Fusarium pininemorale and H. decipiens are not known to cause significant disease but are closely related to species of economic importance. The genome sizes range from 25.99 Mb in the case of O. ips to 4.82 Mb for H. lignivorus. These genomes include the first reports of a genome from the genus Hawksworthiomyces. The availability of these genome data will allow the resolution of longstanding questions regarding the taxonomy of these species. In addition these genome sequences through comparative studies with closely related organisms will increase our understanding of how these species or close relatives cause disease.

Genome-Wide Analysis of Secondary Metabolite Gene Clusters in Ophiostoma ulmi and Ophiostoma novo-ulmi Reveals a Fujikurin-Like Gene Cluster with a Putative Role in Infection

The emergence of new microbial pathogens can result in destructive outbreaks, since their hosts have limited resistance and pathogens may be excessively aggressive. Described as the major ecological incident of the twentieth century, Dutch elm disease, caused by ascomycete fungi from the Ophiostoma genus, has caused a significant decline in elm tree populations (Ulmus sp.) in North America and Europe. Genome sequencing of the two main causative agents of Dutch elm disease (Ophiostoma ulmi and Ophiostoma novo-ulmi), along with closely related species with different lifestyles, allows for unique comparisons to be made to identify how pathogens and virulence determinants have emerged. Among several established virulence determinants, secondary metabolites (SMs) have been suggested to play significant roles during phytopathogen infection. Interestingly, the secondary metabolism of Dutch elm pathogens remains almost unexplored, and little is known about how SM biosynthetic genes are organized in these species. To better understand the metabolic potential of O. ulmi and O. novo-ulmi, we performed a deep survey and description of SM biosynthetic gene clusters (BGCs) in these species and assessed their conservation among eight species from the Ophiostomataceae family. Among 19 identified BGCs, a fujikurin-like gene cluster (OpPKS8) was unique to Dutch elm pathogens. Phylogenetic analysis revealed that orthologs for this gene cluster are widespread among phytopathogens and plant-associated fungi, suggesting that OpPKS8 may have been horizontally acquired by the Ophiostoma genus. Moreover, the detailed identification of several BGCs paves the way for future in-depth research and supports the potential impact of secondary metabolism on Ophiostoma genus’ lifestyle.

Genera of phytopathogenic fungi: GOPHY 1

Genera of Phytopathogenic Fungi (GOPHY) is introduced as a new series of publications in order to provide a stable platform for the taxonomy of phytopathogenic fungi. This first paper focuses on 21 genera of phytopathogenic fungi: Bipolaris, Boeremia, Calonectria, Ceratocystis, Cladosporium, Colletotrichum, Coniella, Curvularia, Monilinia, Neofabraea, Neofusicoccum, Pilidium, Pleiochaeta, Plenodomus, Protostegia, Pseudopyricularia, Puccinia, Saccharata, Thyrostroma, Venturia and Wilsonomyces. For each genus, a morphological description and information about its pathology, distribution, hosts and disease symptoms are provided. In addition, this information is linked to primary and secondary DNA barcodes of the presently accepted species, and relevant literature. Moreover, several novelties are introduced, i.e. new genera, species and combinations, and neo-, lecto- and epitypes designated to provide a stable taxonomy. This first paper includes one new genus, 26 new species, ten new combinations, and four typifications of older names.

Resolved phylogeny and biogeography of the root pathogen Armillaria and its gasteroid relative, Guyanagaster

Background

Armillaria is a globally distributed mushroom-forming genus composed primarily of plant pathogens. Species in this genus are prolific producers of rhizomorphs, or vegetative structures, which, when found, are often associated with infection. Because of their importance as plant pathogens, understanding the evolutionary origins of this genus and how it gained a worldwide distribution is of interest. The first gasteroid fungus with close affinities to Armillaria—Guyanagaster necrorhizus—was described from the Neotropical rainforests of Guyana. In this study, we conducted phylogenetic analyses to fully resolve the relationship of G. necrorhizus with Armillaria. Data sets containing Guyanagaster from two collecting localities, along with a global sampling of 21 Armillaria species—including newly collected specimens from Guyana and Africa—at six loci (28S, EF1α, RPB2, TUB, actin-1 and gpd) were used. Three loci—28S, EF1α and RPB2—were analyzed in a partitioned nucleotide data set to infer divergence dates and ancestral range estimations for well-supported, monophyletic lineages.

Results

The six-locus phylogenetic analysis resolves Guyanagaster as the earliest diverging lineage in the armillarioid clade. The next lineage to diverge is that composed of species in Armillaria subgenus Desarmillaria. This subgenus is elevated to genus level to accommodate the exannulate mushroom-forming armillarioid species. The final lineage to diverge is that composed of annulate mushroom-forming armillarioid species, in what is now Armillaria sensu stricto. The molecular clock analysis and ancestral range estimation suggest the most recent common ancestor to the armillarioid lineage arose 51 million years ago in Eurasia. A new species, Guyanagaster lucianii sp. nov. from Guyana, is described.

Conclusions

The armillarioid lineage evolved in Eurasia during the height of tropical rainforest expansion about 51 million years ago, a time marked by a warm and wet global climate. Species of Guyanagaster and Desarmillaria represent extant taxa of these early diverging lineages. Desarmillaria represents an armillarioid lineage that was likely much more widespread in the past. Guyanagaster likely evolved from a gilled mushroom ancestor and could represent a highly specialized endemic in the Guiana Shield. Armillaria species represent those that evolved after the shift in climate from warm and tropical to cool and arid during the late Eocene. No species in either Desarmillaria or Guyanagaster are known to produce melanized rhizomorphs in nature, whereas almost all Armillaria species are known to produce them. The production of rhizomorphs is an adaptation to harsh environments, and could be a driver of diversification in Armillaria by conferring a competitive advantage to the species that produce them.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0877-3) contains supplementary material, which is available to authorized users.

IMA Genome-F 7: Draft genome sequences for Ceratocystis fagacearum, C. harringtonii, Grosmannia penicillata, and Huntiella bhutanensis

Draft genomes for the fungi Ceratocystis fagacearum, C. harringtonii, Grosmannia penicillata, and Huntiella bhutanensis are presented. Ceratocystis fagacearum is a major causal agent of vascular wilt of oaks and other trees in the family Fagaceae. Ceratocystis harringtonii, previously known as C. populicola, causes disease in Populus species in the USA and Canada. Grosmannia penicillata is the causal agent of bluestain of sapwood on various conifers, including Picea spp. and Pinus spp. in Europe. Huntiella bhutanensis is a fungus in Ceratocystidaceae and known only in association with the bark beetle Ips schmutzenhorferi that infests Picea spinulosa in Bhutan. The availability of these genomes will facilitate further studies on these fungi.