Search

Morchella spongiola is a highly prized mushroom for its delicious flavor and medical value and is one of the most flourishing, representative, and dominant macrofungi in the Qilian Mountains of the Qinghai-Tibet Plateau subkingdoms (QTPs). However, the understanding of M. spongiola remains largely unknown, and its taxonomy is ambiguous. In this study, we redescribed a unique species of M. spongiola, i.e., micromorphology, molecular data, genomics, and comparative genomics, and the historical biogeography of M. spongiola were estimated for 182 single-copy homologous genes. A high-quality chromosome-level reference genome of M. spongiola M12-10 was obtained by combining PacBio HiFi data and Illumina sequencing technologies; it was approximately 57.1 Mb (contig N50 of 18.14 Mb) and contained 9775 protein-coding genes. Comparative genome analysis revealed considerable conservation and unique characteristics between M. spongiola M12-10 and 32 other Morchella species. Molecular phylogenetic analysis indicated that M. spongiola M12-10 is similar to the M. prava/Mes-7 present in sandy soil near rivers, differentiating from black morels ~ 43.06 Mya (million years ago), and diverged from M. parva/Mes-7 at approximately 12.85 Mya (in the Miocene epoch), which is closely related to the geological activities in the QTPs (in the Neogene). Therefore, M. spongiola is a unique species rather than a synonym of M. vulgaris/Mes-5, which has a distinctive grey-brown sponge-like ascomata. This genome of M. spongiola M12-10 is the first published genome sequence of the species in the genus Morchella from the QTPs, which could aid future studies on functional gene identification, germplasm resource management, and molecular breeding efforts, as well as evolutionary studies on the Morchella taxon in the QTPs., (© 2024. The Author(s).)

Dynamic transposition of transposable elements (TEs) in fungal pathogens has significant impact on genome stability, gene expression, and virulence to the host. In Magnaporthe oryzae , genome plasticity resulting from TE insertion is a major driving force leading to the rapid evolution and diversification of this fungus. Despite their importance in M. oryzae population evolution and divergence, our understanding of TEs in this context remains limited. Here, we conducted a genome-wide analysis of TE transposition dynamics in the 11 most abundant TE families in M. oryzae populations. Our results show that these TEs have specifically expanded in recently isolated M. oryzae rice populations, with the presence/absence polymorphism of TE insertions highly concordant with population divergence on Geng/ Japonica and Xian/ Indica rice cultivars. Notably, the genes targeted by clade-specific TEs showed clade-specific expression patterns and are involved in the pathogenic process, suggesting a transcriptional regulation of TEs on targeted genes. Our study provides a comprehensive analysis of TEs in M. oryzae populations and demonstrates a crucial role of recent TE bursts in adaptive evolution and diversification of the M. oryzae rice-infecting lineage., Importance: Magnaporthe oryzae is the causal agent of the destructive blast disease, which caused massive loss of yield annually worldwide. The fungus diverged into distinct clades during adaptation toward the two rice subspecies, Xian/ Indica and Geng/ Japonica . Although the role of TEs in the adaptive evolution was well established, mechanisms underlying how TEs promote the population divergence of M. oryzae remain largely unknown. In this study, we reported that TEs shape the population divergence of M. oryzae by differentially regulating gene expression between Xian/ Indica -infecting and Geng/ Japonica -infecting populations. Our results revealed a TE insertion-mediated gene expression adaption that led to the divergence of M. oryzae population infecting different rice subspecies., Competing Interests: The authors declare no conflict of interest.

The genetic architecture of mating-type loci in lichen-forming fungi has been characterized in very few taxa. Despite the limited data, and in contrast to all other major fungal lineages, arrangements that have both mating-type alleles in a single haploid genome have been hypothesized to be absent from the largest lineage of lichen-forming fungi, the Lecanoromycetes. We report the discovery of both mating-type alleles from the haploid genomes of three species within this group. Our results demonstrate that Lecanoromycetes are not an outlier among Ascomycetes., (Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution 2024.)

Verticillium dahliae is a widespread soilborne fungus that causes Verticillium wilt on numerous economically important plant species. In tomato, until now, three races have been characterized based on the response of differential cultivars to V. dahliae, but the genetic basis of race divergence in V. dahliae remains undetermined. To investigate the genetic basis of race divergence, we sequenced the genomes of two race 2 strains and four race 3 strains for comparative analyses with two known race 1 genomes. The genetic basis of race divergence was described by the pathogenicity-related genes among the three races, orthologue analyses, and genomic structural variations. Global comparative genomics showed that chromosomal rearrangements are not the only source of race divergence and that race 3 should be split into two genotypes based on orthologue clustering. Lineage-specific regions (LSRs), frequently observed between genomes of the three races, encode several predicted secreted proteins that potentially function as suppressors of immunity triggered by known effectors. These likely contribute to the virulence of the three races. Two genes in particular that can act as markers for race 2 and race 3 ( VdR2e and VdR3e , respectively) contribute to virulence on tomato, and the latter acts as an avirulence factor of race 3. We elucidated the genetic basis of race divergence through global comparative genomics and identified secreted proteins in LSRs that could potentially play critical roles in the differential virulence among the races in V. dahliae. IMPORTANCE Deciphering the gene-for-gene relationships during host-pathogen interactions is the basis of modern plant resistance breeding. In the Verticillium dahliae-tomato pathosystem, two races (races 1 and 2) and their corresponding avirulence ( Avr ) genes have been identified, but strains that lack these two Avr genes exist in nature. In this system, race 3 has been described, but the corresponding Avr gene has not been identified. We de novo -sequenced genomes of six strains and identified secreted proteins within the lineage-specific regions (LSRs) distributed among the genomes of the three races that could potentially function as manipulators of host immunity. One of the LSR genes, VdR3e , was confirmed as the Avr gene for race 3. The results indicate that differences in transcriptional regulation may contribute to race differentiation. This is the first study to describe these differences and elucidate roles of secreted proteins in LSRs that play roles in race differentiation.

The lichen, to which the name Lecidea lichenicola is found to have been misapplied, was first described from England and is an extreme specialist of chalk pebbles. It has long been known that it is not closely related to Lecidea in the strict sense, but its true evolutionary relationships have been unknown. Here we use metagenome-assembled genome data to place this fungus in a six-locus phylogeny of Ascomycota, and find strong support for its placement in the class Lichinomycetes. Multiple gene trees using existing data from Lichinomycetes support its further placement within the family Lichinaceae. Based on a revision of types and original descriptions, we conclude that the earliest name for this species is Lecidea obsoleta (syn. Thrombium cretaceum). We neotypify that name by a modern collection and accommodate it in the new genus Watsoniomyces. Type and other original material of L. lichenicola (syn. Discocera lichenicola) was re-examined and found not to be on chalk and to represent a different lichen, Trapelia glebulosa. Watsoniomyces is the first described member of Lichinomycetes with an endolithic thallus., (Copyright © 2021 British Mycological Society. Published by Elsevier Ltd. All rights reserved.)

A key factor to take actions against phytosanitary problems is the accurate and rapid detection of the causal agent. Here, we develop a molecular diagnostics system based on comparative genomics to easily identify fusariosis and specific pathogenic species as the Fusarium kuroshium, the symbiont of the ambrosia beetle Euwallaceae kuroshio Gomez and Hulcr which is responsible for Fusarium dieback disease in San Diego CA, USA. We performed a pan-genome analysis using sixty-three ascomycetes fungi species including phytopathogens and fungi associated with the ambrosia beetles. Pan-genome analysis revealed that 2,631 orthologue genes are only shared by Fusarium spp., and on average 3,941 (SD ± 1,418.6) are species-specific genes. These genes were used for PCR primer design and tested on DNA isolated from i) different strains of ascomycete species, ii) artificially infected avocado stems and iii) plant tissue of field-collected samples presumably infected. Our results let us propose a useful set of primers to either identify any species from Fusarium genus or, in a specific manner, species such as F. kuroshium, F. oxysporum, and F. graminearum. The results suggest that the molecular strategy employed in this study can be expanded to design primers against different types of pathogens responsible for provoking critical plant diseases., Competing Interests: The authors have declared that no competing interests exist.

Monilinia laxa is an important fungal plant pathogen causing brown rot on many stone and pome fruits worldwide. Mitochondrial genome (mitogenome) plays a critical role in evolutionary biology of the organisms. This study aimed to characterize the complete mitogenome of M. laxa by using next-generation sequencing and approaches of de novo assembly and annotation. The total length of the mitogenome of M. laxa was 178,357 bp, and its structure was circular. GC content of the mitogenome was 30.1%. Annotation of the mitogenome presented 2 ribosomal RNA (rRNA) genes, 32 transfer RNA genes (tRNA), 1 gene encoding mitochondrial ribosomal protein S3, 14 protein-coding genes and 15 open reading frame encoding hypothetical proteins. Moreover, the group I mobile introns encoding homing endonucleases including LAGLIDADG and GIY-YIG families were found both within coding regions (genic) and intergenic regions of the mitogenome, indicating an enlarged size and a dynamic structure of the mitogenome. Furthermore, a comparative mitogenomic analysis was performed between M. laxa and the three closely related fungal phytopathogen species (Botryotinia fuckeliana, Sclerotinia sclerotiorum and, S. borealis). Due to the number and distribution of introns, the large extent of structural rearrangements and diverse mitogenome sizes were detected among the species investigated. Monilinia laxa presented the highest number of homing endonucleases among the fungal species considered in the analyses. This study is the first to report a detailed annotation of the mitogenome of an isolate of M. laxa, providing a solid basis for further investigations of mitogenome variations for the other Monilinia pathogens causing brown rot disease.

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.

The symbiont endophytic fungi in tobacco are highly diverse and difficult to classify. Here, we sequenced the genomes of Curvularia trifolii and Leptosphaerulina chartarum isolated from tobacco plants. Finally, 41.68 Mb and 37.95 Mb nuclear genomes were sequenced for C. trifolii and L. chartarum with the scaffold N50, accounting for 638.94 Kb and 284.12 Kb, respectively. Meanwhile, we obtained 68,926 bp and 59,100 bp for their mitochondrial genomes. To more accurately classify C. trifolii and L. chartarum , we extracted seven nuclear genes and 12 mitochondrial genes from these two genomes and their closely related species. The genes were then used for calculation of evolutionary rates and for phylogenetic analysis. Results showed that it was difficult to achieve consistent results using a single gene due to their different evolutionary rates, while the phylogenetic trees obtained by combining datasets showed stable topologies. It is, therefore, more accurate to construct phylogenetic relationships for endophytic fungi based on multi-gene datasets. This study provides new insights into the distribution and characteristics of endophytic fungi in tobacco.

The basal ascomycetes in genus Taphrina have strict host specificity and coevolution with their host plants, making them appealing models for studying the genomic basis of ecological divergence and host adaption. We therefore performed genome sequencing and comparative genomics of different Taphrina species with distinct host ranges to reveal their evolution. We identified frequent chromosomal rearrangements and highly dynamic lineage-specific (LS) genomic regions in Taphrina genomes. The LS regions occur at the flanking regions of chromosomal breakpoints, and are greatly enriched for DNA repeats, non-core genes, and in planta up-regulated genes. Furthermore, we identified hundreds of candidate secreted effector proteins (CSEPs) that are commonly organized in gene clusters that form distinct AT-rich isochore-like regions. Nearly half of the CSEPs constitute two novel superfamilies with modular structures unique to Taphrina. These CSEPs are commonly up-regulated during infection, enriched in the LS regions, evolved faster, and underwent extensive gene gain and loss in different species. In addition to displaying signatures of positive selection, functional characterization of selected CSEP genes confirmed their roles in suppression of plant defence responses. Overall, our results showed that extensive chromosomal rearrangements and rapidly evolving CSEP superfamilies play important roles in speciation and host adaptation in the early-branching ascomycetous fungi., (© 2020 The Authors. Molecular Plant Pathology published by British Society for Plant Pathology and John Wiley & Sons Ltd.)

Shiraia bambusicola is a rare medicinal fungus found in China that causes bamboo plants to decay and die with severe infection. Hypocrellin, its main active ingredient, is widely used in several fields, such as medicine, agriculture, and food industry. In this study, to clarify the genomic components, taxonomic status, pathogenic genes, secondary metabolite synthesis pathways, and regulatory mechanisms of S. bambusicola , whole-genome sequencing, assembly, and functional annotation were performed using high-throughput sequencing and bioinformatics approaches. It was observed that S. bambusicola has 33 Mb genome size, 48.89% GC content, 333 scaffolds, 2590 contigs, 10,703 genes, 82 tRNAs, and 21 rRNAs. The total length of the repeat sequence is 2,151,640 bp. The annotation of 5945 proteins was obtained from InterProScan hits based on the Gene Ontology database. Phylogenetic analysis showed that S. bambusicola belongs to Shiraiaceae, a new family of Pleosporales. It was speculated that there are more than two species or genus in Shiraiaceae. According to the annotation, 777 secreted proteins were associated with virulence or detoxification, including 777 predicted by the PHI database, 776 by the CAZY and Fungal Cytochrome P 450 database, and 441 by the Proteases database. The 252 genes associated with the secondary metabolism of S. bambusicola were screened and enriched into 28 pathways, among which the terpenoids, staurosporine, aflatoxin, and folate synthesis pathways have not been reported in S. bambusicola The T1PKS was the main gene cluster among the 28 secondary metabolite synthesis gene clusters in S. bambusicola The analysis of the T3PKS gene cluster related to the synthesis of hypocrellin showed that there was some similarity between S. bambusicola and 10 other species of fungi; however, the similarity was very low wherein the highest similarity was 17%. The genomic information of S. bambusicola obtained in this study was valuable to understand its genetic function and pathogenicity. The genomic information revealed that several enzyme genes and secreted proteins might be related to their host interactions and pathogenicity. The annotation and analysis of its secondary metabolite synthesis genes and gene clusters will be an important reference for future studies on the biosynthesis and regulation mechanism of the secondary metabolites, contributing to the discovery of new metabolites and accelerating drug development and application., (Copyright © 2020 Ren et al.)

The Venturia genus comprises fungal species that are pathogens on Rosaceae host plants, including V. inaequalis and V. asperata on apple, V. aucupariae on sorbus and V. pirina on pear. Although the genetic structure of V. inaequalis populations has been investigated in detail, genomic features underlying these subdivisions remain poorly understood. Here, we report whole genome sequencing of 87 Venturia strains that represent each species and each population within V. inaequalis We present a PacBio genome assembly for the V. inaequalis EU-B04 reference isolate. The size of selected genomes was determined by flow cytometry, and varied from 45 to 93 Mb. Genome assemblies of V. inaequalis and V. aucupariae contain a high content of transposable elements (TEs), most of which belong to the Gypsy or Copia LTR superfamilies and have been inactivated by Repeat-Induced Point mutations. The reference assembly of V. inaequalis presents a mosaic structure of GC-equilibrated regions that mainly contain predicted genes and AT-rich regions, mainly composed of TEs. Six pairs of strains were identified as clones. Single-Nucleotide Polymorphism (SNP) analysis between these clones revealed a high number of SNPs that are mostly located in AT-rich regions due to misalignments and allowed determining a false discovery rate. The availability of these genome sequences is expected to stimulate genetics and population genomics research of Venturia pathogens. Especially, it will help understanding the evolutionary history of Venturia species that are pathogenic on different hosts, a history that has probably been substantially influenced by TEs., (Copyright © 2019 Le Cam et al.)

Purpureocillium lilacinum has been widely used as a commercial biocontrol agent for the control of plant parasitic nematodes. Whole genome analysis promotes the identification of functional genes and the exploration of their molecular mechanisms. The Cyclophilin ( CYP ) gene family belongs to the immunophillin superfamily, and has a conserved cyclophilin-like domain (CLD). CYPs are widely identified in prokaryotes and eukaryotes, and can be divided into single- and multi-domain proteins. In the present study, 10 CYP genes possessing the CLD, named PlCYP1-P10 , were identified from the genome of P. lilacinum strain 36-1. Those 10 PlCYPs were predicted to have different cellular localizations in P. lilacinum . Phylogenetic and gene structure analysis revealed the evolutionary differentiation of CYPs between Ascomycotina and Saccharomycotina fungi, but conservation within the Ascomycotina fungi. Motif and gene structure distributions further support the result of phylogenetic analysis. Each PlCYP gene had a specific expression pattern in different development stages of P. lilacinum and its parasitism stage on eggs of Meloidogyne incognita . In addition, the 10 PlCYP genes exhibited different expression abundances in response to abiotic stresses, among which PlCYP4 was highly expressed at a high temperature (35 °C), while PlCYP6 was up-regulated under 5 mM of H 2 O 2 stress. Furthermore, the heterologous expression of PlCYP4 and PlCYP6 in Escherichia coli enhanced the cellular tolerance against a high temperature and H 2 O 2 . In summary, our study indicates the potential functions of PlCYPs in virulence and the stress response, and also provides a frame for further analysis of the CYP gene family in Ascomycotina fungi.

Synchytrium endobioticum is an obligate biotrophic soilborne Chytridiomycota (chytrid) species that causes potato wart disease, and represents the most basal lineage among the fungal plant pathogens. We have chosen a functional genomics approach exploiting knowledge acquired from other fungal taxa and compared this to several saprobic and pathogenic chytrid species. Observations linked to obligate biotrophy, genome plasticity and pathogenicity are reported. Essential purine pathway genes were found uniquely absent in S. endobioticum, suggesting that it relies on scavenging guanine from its host for survival. The small gene-dense and intron-rich chytrid genomes were not protected for genome duplications by repeat-induced point mutation. Both pathogenic chytrids Batrachochytrium dendrobatidis and S. endobioticum contained the largest amounts of repeats, and we identified S. endobioticum specific candidate effectors that are associated with repeat-rich regions. These candidate effectors share a highly conserved motif, and show isolate specific duplications. A reduced set of cell wall degrading enzymes, and LysM protein expansions were found in S. endobioticum, which may prevent triggering plant defense responses. Our study underlines the high diversity in chytrids compared to the well-studied Ascomycota and Basidiomycota, reflects characteristic biological differences between the phyla, and shows commonalities in genomic features among pathogenic fungi.

Background: Many small peptides regulate eukaryotic cell biology. In fungi, some of these peptides are produced after KEX2 protease activity on proteins displaying repetitions of identical or nearly identical motifs. Following this endoprotease activity, peptides are released in the extracellular space. This type of protein maturation is involved in the production of the α-type sexual pheromone in Ascomycota. In other cases, this processing allows the production of secreted peptides regulating fungal cell wall structure or acting as mycotoxins. In this work, we report for the first time a genome-wide search of KEX2-processed repeat proteins that we call KEPs. We screened the secreted proteins of 250 fungal species to compare their KEP repertoires with regard to their lifestyle, morphology or lineage., Results: Our analysis points out that nearly all fungi display putative KEPs, suggesting an ancestral origin common to all opisthokonts. As expected, our pipeline identifies mycotoxins but also α-type sexual pheromones in Ascomycota that have not been explored so far, and unravels KEP-derived secreted peptides of unknown functions. Some species display an expansion of this class of proteins. Interestingly, we identified conserved KEPs in pathogenic fungi, suggesting a role in virulence. We also identified KEPs in Basidiomycota with striking similarities to Ascomycota α-type sexual pheromones, suggesting they may also play alternative roles in unknown signalling processes., Conclusions: We identified putative, new, unexpected secreted peptides that fall into different functional categories: mycotoxins, hormones, sexual pheromones, or effectors that promote colonization during host-microbe interactions. This wide survey will open new avenues in the field of small-secreted peptides in fungi that are critical regulators of their intimate biology and modulators of their interaction with the environment.

Sodiomyces alkalinus is one of the very few alkalophilic fungi, adapted to grow optimally at high pH. It is widely distributed at the plant-deprived edges of extremely alkaline lakes and locally abundant. We sequenced the genome of S. alkalinus and reconstructed evolution of catabolic enzymes, using a phylogenomic comparison. We found that the genome of S. alkalinus is larger, but its predicted proteome is smaller and heavily depleted of both plant-degrading enzymes and proteinases, when compared to its closest plant-pathogenic relatives. Interestingly, despite overall losses, S. alkalinus has retained many proteinases families and acquired bacterial cell wall-degrading enzymes, some of them via horizontal gene transfer from bacteria. This fungus has very potent proteolytic activity at high pH values, but slowly induced low activity of cellulases and hemicellulases. Our experimental and in silico data suggest that plant biomass, a common food source for most fungi, is not a preferred substrate for S. alkalinus in its natural environment. We conclude that the fungus has abandoned the ancestral plant-based diet and has become specialized in a more protein-rich food, abundantly available in soda lakes in the form of prokaryotes and small crustaceans., (2018 The Authors. Molecular Ecology Published by John Wiley & Sons Ltd.)

Background: Powdery mildew (PM) is one of the most important and widespread plant diseases caused by biotrophic fungi. Notably, while monocot (grass) PM fungi exhibit high-level of host-specialization, many dicot PM fungi display a broad host range. To understand such distinct modes of host-adaptation, we sequenced the genomes of four dicot PM biotypes belonging to Golovinomyces cichoracearum or Oidium neolycopersici., Results: We compared genomes of the four dicot PM together with those of Blumeria graminis f.sp. hordei (both DH14 and RACE1 isolates), B. graminis f.sp. tritici, and Erysiphe necator infectious on barley, wheat and grapevine, respectively. We found that despite having a similar gene number (6620-6961), the PM genomes vary from 120 to 222 Mb in size. This high-level of genome size variation is indicative of highly differential transposon activities in the PM genomes. While the total number of genes in any given PM genome is only about half of that in the genomes of closely related ascomycete fungi, most (~ 93%) of the ascomycete core genes (ACGs) can be found in the PM genomes. Yet, 186 ACGs were found absent in at least two of the eight PM genomes, of which 35 are missing in some dicot PM biotypes, but present in the three monocot PM genomes, indicating remarkable, independent and perhaps ongoing gene loss in different PM lineages. Consistent with this, we found that only 4192 (3819 singleton) genes are shared by all the eight PM genomes, the remaining genes are lineage- or biotype-specific. Strikingly, whereas the three monocot PM genomes possess up to 661 genes encoding candidate secreted effector proteins (CSEPs) with families containing up to 38 members, all the five dicot PM fungi have only 116-175 genes encoding CSEPs with limited gene amplification., Conclusions: Compared to monocot (grass) PM fungi, dicot PM fungi have a much smaller effectorome. This is consistent with their contrasting modes of host-adaption: while the monocot PM fungi show a high-level of host specialization, which may reflect an advanced host-pathogen arms race, the dicot PM fungi tend to practice polyphagy, which might have lessened selective pressure for escalating an with a particular host.

Morchella is a popular edible fungus worldwide due to its rich nutrition and unique flavor. Many research efforts were made on the domestication and cultivation of Morchella all over the world. In recent years, the cultivation of Morchella was successfully commercialized in China. However, the biology is not well understood, which restricts the further development of the morel fungus cultivation industry. In this paper, we performed de novo sequencing and assembly of the genomes of two monospores with a different mating type (M04M24 and M04M26) isolated from the commercially cultivated strain M04. Gene annotation and comparative genome analysis were performed to study differences in CAZyme (Carbohydrate-active enzyme) enzyme content, transcription factors, duplicated sequences, structure of mating type sites, and differences at the gene and functional levels between the two monospore strains of M. importuna . Results showed that the de novo assembled haploid M04M24 and M04M26 genomes were 48.98 and 51.07 Mb, respectively. A complete fine physical map of M. importuna was obtained from genome coverage and gene completeness evaluation. A total of 10,852 and 10,902 common genes and 667 and 868 endemic genes were identified from the two monospore strains, respectively. The Gene Ontology (GO) and KAAS (KEGG Automatic Annotation Serve) enrichment analyses showed that the endemic genes performed different functions. The two monospore strains had 99.22% collinearity with each other, accompanied with certain position and rearrangement events. Analysis of complete mating-type loci revealed that the two monospore M. importuna strains contained an independent mating-type structure and remained conserved in sequence and location. The phylogenetic and divergence time of M. importuna was analyzed at the whole-genome level for the first time. The bifurcation time of morel and tuber was estimated to be 201.14 million years ago (Mya); the two monospore strains with a different mating type represented the evolution of different nuclei, and the single copy homologous genes between them were also different due to a genetic differentiation distance about 0.65 Mya. Compared with truffles, M. importuna had an extension of 28 clusters of orthologous genes (COGs) and a contraction of two COGs. The two different polar nuclei with different degrees of contraction and expansion suggested that they might have undergone different evolutionary processes. The different mating-type structures, together with the functional clustering and enrichment analysis results of the endemic genes of the two different polar nuclei, imply that M. importuna might be a heterothallic fungus and the interaction between the endemic genes may be necessary for its complete life history. Studies on the genome of M. importuna facilitate a better understanding of morel biology and evolution.

This study describes a novel ericoid mycorrhizal fungus (ErMF), Gamarada debralockiae Midgley and Tran-Dinh gen. nov. sp. nov. Additionally, catabolism was explored from a genomic perspective. The nuclear and mitochondrial genomes of G. debralockiae were sequenced. Morphological characteristics were assessed on various media. Catabolic genes of G. debralockiae were explored using SignalP and dbCAN. Phylogenetic comparisons were undertaken using Phylogeny.fr. The 58.5-Mbp draft genome of G. debralockiae contained 17,075 putative genes. The complete mitochondrial genome was 28,168 bp in length. In culture, G. debralockiae produces slow-growing non-sporulating colonies. Gamarada debralockiae has many putative secreted catabolic enzymes. Phylogeny indicated G. debralockiae was distinct from known ascomycetous ErMF: Pezoloma ericae, Meliniomyces spp., Oidiodendron spp., and Cairneyella variabilis. It is closely related to many undescribed plant root-associated fungi and its nearest described relative is Hyphodiscus brevicollaris. Gamarada debralockiae has been recovered from virtually all Australian ericoid mycorrhizal studies and biogeographic data suggests the taxon is widespread in Australia. Gamarada debralockiae has similar catabolic potential to C. variabilis and co-occurs with C. variabilis at Australian sites. Plants that host multiple ErMF may benefit from subtle differences in catabolism that improve access to nitrogen and phosphorus from within recalcitrant organic matter.

Yeasts of the Meyerozyma guilliermondii species complex are widespread in nature and can be isolated from a variety of sources, from the environment to arthropods to hospital patients. To date, the species complex comprises the thoroughly studied and versatile M. guilliermondii , the hard to distinguish M. caribbica , and Candida carpophila Here we report the whole genome sequencing and de novo assembly of four M. caribbica isolates, identified with the most recent molecular techniques, derived from four Diptera species. The four novel assemblies present reduced fragmentation and comparable metrics (genome size, gene content) to the available genomes belonging to the species complex. We performed a phylogenomic analysis comprising all known members of the species complex, to investigate evolutionary relationships within this clade. Our results show a compact phylogenetic structure for the complex and indicate the presence of a sizable core set of genes. Furthermore, M. caribbica , despite a broad literature on the difficulties of discerning it from M. guilliermondii , seems to be more closely related to C. carpophila Finally, we believe that there is evidence for considering these four genomes to be the first published for the species M. caribbica Raw reads and assembled contigs have been made public to further the study of these organisms., (Copyright © 2018 De Marco et al.)

Parastagonospora nodorum , the causal agent of Septoria nodorum blotch in wheat, has emerged as a model necrotrophic fungal organism for the study of host-microbe interactions. To date, three necrotrophic effectors have been identified and characterized from this pathogen, including SnToxA, SnTox1, and SnTox3. Necrotrophic effector identification was greatly aided by the development of a draft genome of Australian isolate SN15 via Sanger sequencing, yet it remained largely fragmented. This research presents the development of nearly finished genomes of P. nodorum isolates Sn4, Sn2000, and Sn79-1087 using long-read sequencing technology. RNAseq analysis of isolate Sn4, consisting of eight time points covering various developmental and infection stages, mediated the annotation of 13,379 genes. Analysis of these genomes revealed large-scale polymorphism between the three isolates, including the complete absence of contig 23 from isolate Sn79-1087, and a region of genome expansion on contig 10 in isolates Sn4 and Sn2000. Additionally, these genomes exhibit the hallmark characteristics of a "two-speed" genome, being partitioned into two distinct GC-equilibrated and AT-rich compartments. Interestingly, isolate Sn79-1087 contains a lower proportion of AT-rich segments, indicating a potential lack of evolutionary hotspots. These newly sequenced genomes, consisting of telomere-to-telomere assemblies of nearly all 23 P. nodorum chromosomes, provide a robust foundation for the further examination of effector biology and genome evolution., (Copyright © 2018 Richards et al.)

Synteny can be maintained for certain genomic regions across broad phylogenetic groups. In these homologous genomic regions, sites that are under relaxed purifying selection, such as intergenic regions, could be used broadly as markers for population genetic and phylogenetic studies on species complexes. To explore the potential of this approach, we found 125 Collinear Orthologous Regions (COR) ranging from 1 to >10kb across nine genomes representing the Lecanoromycetes and Eurotiomycetes (Pezizomycotina, Ascomycota). Twenty-six of these COR were found in all 24 eurotiomycete genomes surveyed for this study. Given the high abundance and availability of fungal genomes we believe this approach could be adopted for other large groups of fungi outside the Pezizomycotina. Asa proof of concept, we selected three Collinear Orthologous Regions (COR1b, COR3, and COR16), based on synteny analyses of several genomes representing three classes of Ascomycota: Eurotiomycetes, Lecanoromycetes, and Lichinomycetes. COR16, for example, was found across these three classes of fungi. Here we compare the resolving power of these three new markers with five loci commonly used in phylogenetic studies of fungi, using section Polydactylon of the cyanolichen-forming genus Peltigera (Lecanoromycetes) - a clade with several challenging species complexes. Sequence data were subjected to three species discovery and two validating methods. COR markers substantially increased phylogenetic resolution and confidence, and highly contributed to species delimitation. The level of phylogenetic signal provided by each of the COR markers was higher than the commonly used fungal barcode ITS. High cryptic diversity was revealed by all methods. As redefined here, most species represent lineages that have relatively narrower, and more homogeneous biogeographical ranges than previously understood. The scabrosoid clade consists of ten species, seven of which are new. For the dolichorhizoid clade, twenty-two new species were discovered for a total of twenty-nine species in this clade., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Spliceosomal introns can occupy nearby rather than identical positions in orthologous genes (intron sliding or shifting). Stwintrons are complex intervening sequences, where an 'internal' intron interrupts one of the sequences essential for splicing, generating after its excision, a newly formed canonical intron defined as 'external'. In one experimentally demonstrated configuration, two alternatively excised internal introns, overlapping by one G, disrupt respectively the donor and the acceptor sequence of an external intron, leading to mRNAs encoding identical proteins. In a gene encoding a DHA1 antiporter in Pezizomycotina, we find a variety of predicted intron configurations interrupting the DNA stretch encoding a conserved peptidic sequence. Some sport a stwintron where the internal intron interrupts the donor of the external intron (experimentally confirmed for Aspergillus nidulans). In others, we found and demonstrate (for Trichoderma reesei) alternative, overlapping internal introns. Discordant canonical introns, one nt apart, are present in yet other species, exactly as predicted by the alternative loss of either of the internal introns at the DNA level from an alternatively spliced stwintron. An evolutionary pathway of 1 nt intron shift, involving an alternatively spliced stwintron intermediate is proposed on the basis of the experimental and genomic data presented., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Background: G-protein coupled receptors (GPCRs) are the largest family of transmembrane receptors in fungi, where they play important roles in signal transduction. Among them, the Pth11-related GPCRs form a large and divergent protein family, and are only found in fungi in Pezizomycotina. However, the evolutionary process and potential functions of Pth11-related GPCRs remain largely unknown., Results: Twenty genomes of fungi in Pezizomycotina covering different nutritional strategies were mined for putative Pth11-related GPCRs. Phytopathogens encode much more putative Pth11-related GPCRs than symbionts, saprophytes, or entomopathogens. Based on the phylogenetic tree, these GPCRs can be divided into nine clades, with each clade containing fungi in different taxonomic orders. Instead of fungi from the same order, those fungi with similar nutritional strategies were inclined to share orthologs of putative Pth11-related GPCRs. Most of the CFEM domain-containing Pth11-related GPCRs, which were only included in two clades, were detected in phytopathogens. Furthermore, many putative Pth11-related GPCR genes of phytopathogens were upregulated during invasive plant infection, but downregulated under biotic stress. The expressions of putative Pth11-related GPCR genes of saprophytes and entomopathogens could be affected by nutrient conditions, especially the carbon source. The gene expressions revealed that Pth11-related GPCRs could respond to biotic/abiotic stress and invasive plant infection with different expression patterns., Conclusion: Our results indicated that the Pth11-related GPCRs existed before the diversification of Pezizomycotina and have been gained and/or lost several times during the evolutionary process. Tandem duplications and trophic variations have been important factors in this evolution.

Despite the recent surge in mitochondrial (mt) genome sequencing, Kingdom Fungi remains underrepresented with respect to mtDNA. We describe the mt genome of the conifer needle endophyte, Phialocephala scopiformis DAOMC 229536 (Helotiales, Ascomycota). This strain is of interest to the Canadian forestry industry as it produces the anti-insectan compound rugulosin. Sequence was obtained from whole genome shotgun sequencing. Comparison to the only other published Phialocephala mt genome, Phialocephala subalpina, indicates that the suite of common mt genes - cox1-3, cob, nad1-6 and 4L, atp6, 8 and 9, as well as rrnL and rrnS - has retained an identical order. Nad4L remains one of the most conserved mitochondrial genes within Phialocephala. Members of the closely related Phialocephala fortinii s.l. - Acephala appalanata species complex (PAC) share too much sequence similarity to properly resolve lineages using ITS barcoding alone. Using P. scopiformis sequence as an outgroup, we determined ancestral gene states that help confirm clades within Phialocephala. Our results show: (1) the complete mt genome of P. scopiformis, representing the 10th complete mt genome for the order Helotiales (containing >3800 species), and (2) how large-scale genomic patterns, such as mitochondrial gene order, can be used to confirm lineages within fungal species complexes., (Copyright © 2016 British Mycological Society. All rights reserved.)

Background: Renewable biopolymers, such as cellulose, starch and chitin are highly resistance to enzymatic degradation. Therefore, there is a need to upgrade current degradation processes by including novel enzymes. Lytic polysaccharide mono-oxygenases (LPMOs) can disrupt recalcitrant biopolymers, thereby enhancing hydrolysis by conventional enzymes. However, novel LPMO families are difficult to identify using existing methods. Therefore, we developed a novel profile Hidden Markov model (HMM) and used it to mine genomes of ascomycetous fungi for novel LPMOs., Results: We constructed a structural alignment and verified that the alignment was correct. In the alignment we identified several known conserved features, such as the histidine brace and the N/Q/E-X-F/Y motif and previously unidentified conserved proline and glycine residues. These residues are distal from the active site, suggesting a role in structure rather than activity. The multiple protein alignment was subsequently used to build a profile Hidden Markov model. This model was initially tested on manually curated datasets and proved to be both sensitive (no false negatives) and specific (no false positives). In some of the genomes analyzed we identified a yet unknown LPMO family. This new family is mostly confined to the phyla of Ascomycota and Basidiomycota and the class of Oomycota. Genomic clustering indicated that at least some members might be involved in the degradation of β-glucans, while transcriptomic data suggested that others are possibly involved in the degradation of pectin., Conclusions: The newly developed profile hidden Markov Model was successfully used to mine fungal genomes for a novel family of LPMOs. However, the model is not limited to bacterial and fungal genomes. This is illustrated by the fact that the model was also able to identify another new LPMO family in Drosophila melanogaster. Furthermore, the Hidden Markov model was used to verify the more distant blast hits from the new fungal family of LPMOs, which belong to the Bivalves, Stony corals and Sea anemones. So this Hidden Markov model (Additional file 3) will help the broader scientific community in identifying other yet unknown LPMOs.

The range of hosts that parasites can infect is a key determinant of the emergence and spread of disease. Yet, the impact of host range variation on the evolution of parasite genomes remains unknown. Here, we show that codon optimization underlies genome adaptation in broad host range parasites. We found that the longer proteins encoded by broad host range fungi likely increase natural selection on codon optimization in these species. Accordingly, codon optimization correlates with host range across the fungal kingdom. At the species level, biased patterns of synonymous substitutions underpin increased codon optimization in a generalist but not a specialist fungal pathogen. Virulence genes were consistently enriched in highly codon-optimized genes of generalist but not specialist species. We conclude that codon optimization is related to the capacity of parasites to colonize multiple hosts. Our results link genome evolution and translational regulation to the long-term persistence of generalist parasitism.

Objective: To explore the genome sequence of Aureobasidium pullulans CCTCC M2012223, analyze the key genes related to the biosynthesis of important metabolites, and provide genetic background for metabolic engineering., Methods: Complete genome of A. pullulans CCTCC M2012223 was sequenced by Illumina HiSeq high throughput sequencing platform. Then, fragment assembly, gene prediction, functional annotation, and GO/COG cluster were analyzed in comparison with those of other five A. pullulans varieties., Results: The complete genome sequence of A. pullulans CCTCC M2012223 was 30756831 bp with an average GC content of 47.49%, and 9452 genes were successfully predicted. Genome-wide analysis showed that A. pullulans CCTCC M2012223 had the biggest genome assembly size. Protein sequences involved in the pullulan and polymalic acid pathway were highly conservative in all of six A. pullulans varieties. Although both A. pullulans CCTCC M2012223 and A. pullulans var. melanogenum have a close affinity, some point mutation and inserts were occurred in protein sequences involved in melanin biosynthesis., Conclusion: Genome information of A. pullulans CCTCC M2012223 was annotated and genes involved in melanin, pullulan and polymalic acid pathway were compared, which would provide a theoretical basis for genetic modification of metabolic pathway in A. pullulans.

Sporothrix globosa is a thermo-dimorphic fungus belonging to a pathogenic clade that also includes Sporothrix schenckii, which causes human and animal sporotrichosis. Here, we present the first genome assemblies of two S. globosa strains providing data for future comparative genomic studies in pathogenic Sporothrix species., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Understanding the phylogenetic relationships among the yeasts of the subphylum Saccharomycotina is a prerequisite for understanding the evolution of their metabolisms and ecological lifestyles. In the last two decades, the use of rDNA and multilocus data sets has greatly advanced our understanding of the yeast phylogeny, but many deep relationships remain unsupported. In contrast, phylogenomic analyses have involved relatively few taxa and lineages that were often selected with limited considerations for covering the breadth of yeast biodiversity. Here we used genome sequence data from 86 publicly available yeast genomes representing nine of the 11 known major lineages and 10 nonyeast fungal outgroups to generate a 1233-gene, 96-taxon data matrix. Species phylogenies reconstructed using two different methods (concatenation and coalescence) and two data matrices (amino acids or the first two codon positions) yielded identical and highly supported relationships between the nine major lineages. Aside from the lineage comprised by the family Pichiaceae, all other lineages were monophyletic. Most interrelationships among yeast species were robust across the two methods and data matrices. However, eight of the 93 internodes conflicted between analyses or data sets, including the placements of: the clade defined by species that have reassigned the CUG codon to encode serine, instead of leucine; the clade defined by a whole genome duplication; and the species Ascoidea rubescens These phylogenomic analyses provide a robust roadmap for future comparative work across the yeast subphylum in the disciplines of taxonomy, molecular genetics, evolutionary biology, ecology, and biotechnology. To further this end, we have also provided a BLAST server to query the 86 Saccharomycotina genomes, which can be found at http://y1000plus.org/blast., (Copyright © 2016 Shen et al.)

Background: The Rhynchosporium species complex consists of hemibiotrophic fungal pathogens specialized to different sweet grass species including the cereal crops barley and rye. A sexual stage has not been described, but several lines of evidence suggest the occurrence of sexual reproduction. Therefore, a comparative genomics approach was carried out to disclose the evolutionary relationship of the species and to identify genes demonstrating the potential for a sexual cycle. Furthermore, due to the evolutionary very young age of the five species currently known, this genus appears to be well-suited to address the question at the molecular level of how pathogenic fungi adapt to their hosts., Results: The genomes of the different Rhynchosporium species were sequenced, assembled and annotated using ab initio gene predictors trained on several fungal genomes as well as on Rhynchosporium expressed sequence tags. Structures of the rDNA regions and genome-wide single nucleotide polymorphisms provided a hypothesis for intra-genus evolution. Homology screening detected core meiotic genes along with most genes crucial for sexual recombination in ascomycete fungi. In addition, a large number of cell wall-degrading enzymes that is characteristic for hemibiotrophic and necrotrophic fungi infecting monocotyledonous hosts were found. Furthermore, the Rhynchosporium genomes carry a repertoire of genes coding for polyketide synthases and non-ribosomal peptide synthetases. Several of these genes are missing from the genome of the closest sequenced relative, the poplar pathogen Marssonina brunnea, and are possibly involved in adaptation to the grass hosts. Most importantly, six species-specific genes coding for protein effectors were identified in R. commune. Their deletion yielded mutants that grew more vigorously in planta than the wild type., Conclusion: Both cryptic sexuality and secondary metabolites may have contributed to host adaptation. Most importantly, however, the growth-retarding activity of the species-specific effectors suggests that host adaptation of R. commune aims at extending the biotrophic stage at the expense of the necrotrophic stage of pathogenesis. Like other apoplastic fungi Rhynchosporium colonizes the intercellular matrix of host leaves relatively slowly without causing symptoms, reminiscent of the development of endophytic fungi. Rhynchosporium may therefore become an object for studying the mutualism-parasitism transition.

Ascomycete yeasts are metabolically diverse, with great potential for biotechnology. Here, we report the comparative genome analysis of 29 taxonomically and biotechnologically important yeasts, including 16 newly sequenced. We identify a genetic code change, CUG-Ala, in Pachysolen tannophilus in the clade sister to the known CUG-Ser clade. Our well-resolved yeast phylogeny shows that some traits, such as methylotrophy, are restricted to single clades, whereas others, such as l-rhamnose utilization, have patchy phylogenetic distributions. Gene clusters, with variable organization and distribution, encode many pathways of interest. Genomics can predict some biochemical traits precisely, but the genomic basis of others, such as xylose utilization, remains unresolved. Our data also provide insight into early evolution of ascomycetes. We document the loss of H3K9me2/3 heterochromatin, the origin of ascomycete mating-type switching, and panascomycete synteny at the MAT locus. These data and analyses will facilitate the engineering of efficient biosynthetic and degradative pathways and gateways for genomic manipulation., Competing Interests: C.H.C. and T.W.J. are employees of Xylome Corporation, which is developing nonconventional yeasts for biotechnological applications.

The fungal pathogen, Ceratocystis manginecans, has caused serious canker and wilt disease on Mangifera indica (mango), on legume tree species in Oman and Pakistan and on Acacia spp. in Indonesia. A Ceratocystis species, with similar morphology to C. manginecans, has recently been reported in Vietnam, causing severe disease of Acacia trees. Previous population genetic studies on isolates from M. indica in Oman and Pakistan have shown that the pathogen represents a single clonal haplotype, indicative of an introduced pathogen. The aim of this study was to investigate the genetic diversity and population structure of 160 C. manginecans isolates, from four host-associated populations from Oman, Pakistan, Indonesia and Vietnam. This was done by applying a combination of 14 previously developed microsatellite markers and a new set, designed in this study from two different C. manginecans genomes. Sequence data confirmed that the isolates in Vietnam are the same species as those in Indonesia and were thus identified as C. manginecans. Unlike the populations in Oman and Pakistan, relatively high levels of genetic variation were found for the isolates from Indonesia and Vietnam. The Vietnam population was significantly differentiated from the other populations and isolates from this area had the highest level of genetic diversity thus far encountered for the pathogen., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Background: Aschersonia badia [(Ab) Teleomorph: Hypocrella siamensis] is an entomopathogenic fungus that specifically infects scale insects and whiteflies. We present the whole genome sequence of Ab and its comparison with two clavicipitaceous fungi Metarhizium robertsii (MR: generalist entomopathogen) and M. acridum (MAC: acridid-specific entomopathogen) that exhibit variable host preferences. Here, through comparative analysis of pathogen-host interacting genes, carbohydrate active enzymes, secondary metabolite biosynthesis genes, and sexuality genes, we explore the proteins with possible virulence functions in clavicipitaceous fungi. Comprehensive overview of GH18 family chitinases has been provided to decipher the role of chitinases in claviceptaceous fungi that are either host specific or generalists., Results: We report the 28.8 Mb draft genome of Ab and its comparative genome analysis with MR and MAC. The comparative analyses suggests expansion in pathogen-host interacting gene families and carbohydrate active enzyme families in MR, whilst their contraction in Ab and MAC genomes. The multi-modular NRPS gene (dtxS1) responsible for biosynthesis of the secondary metabolite destruxin in MR is not conserved in Ab, similar to the specialist pathogen MAC. An additional siderophore biosynthetic gene responsible for acquisition of iron was identified in MR. Further, the domain survey of chitinases suggest that the CBM50 (LysM) domains, which participate in chitin-binding functions, were not observed in MAC, but were present in Ab and MR. However, apparent differences in frequency of CBM50 domains associated with chitinases of Ab and MR was identified, where MR chitinases displayed a higher proportion of associated CBM50 domains than Ab chitinases., Conclusions: This study suggests differences in distribution of dtxS1 and chitinases in specialists (Ab and MAC) and generalists (MR) fungi. Our analysis also suggests the presence of a siderophore biosynthetic gene in the MR genome which perhaps aids in enhanced virulence potential and host range. The variation in association of CBMs, being higher in generalists (MR) and lower in specialists (Ab and MAC) fungi may further be responsible for the differences in host affiliation.

This paper describes a novel species of ericoid mycorrhizal fungus from Australia, Cairneyella variabilis, Midgley and Tran-Dinh, gen. nov. sp. nov. The genome of C. variabilis was sequenced and a draft genome assembled. The draft genome of C. variabilis is 52.4 Mbp in length, and to our knowledge, this is the first study to present a genome of an ericoid mycorrhizal fungus from the southern hemisphere. Using the SignalP and dbCAN bioinformatic pipelines, a study of the catabolic potential of C. variabilis was undertaken and showed genes for an array of degradative enzymes, most of which appear to be secreted from the hyphae, to access a suite of different carbon sources. Isolates of C. variabilis have been previously shown to utilise cellulose, carboxymethyl cellulose (CMC), cellobiose, xylan, pectin, starch and tannic acid for growth, and in the current study, putative enzymes for these processes were revealed. These enzymes likely play key roles in nutrient cycling and other edaphic processes in heathland environments. ITS phylogenetic analyses showed C. variabilis to be distinct from the fungi of the "Hymenoscyphus ericae aggregate".

Background: Magnaporthaceae, a family of ascomycetes, includes three fungi of great economic importance that cause disease in cereal and turf grasses: Magnaporthe oryzae (rice blast), Gaeumannomyces graminis var. tritici (take-all disease), and Magnaporthe poae (summer patch disease). Recently, the sequenced and assembled genomes for these three fungi were reported. Here, the genomes were compared for orthologous genes in order to identified genes that are unique to the Magnaporthaceae family of fungi. In addition, ortholog clustering was used to identify a core proteome for the Magnaporthaceae, which was examined for diversifying and purifying selection and evidence of two-speed genome evolution., Results: A genome-scale comparative study was conducted across 74 fungal genomes to identify clusters of orthologous genes unique to the three Magnaporthaceae species as well as species specific genes. We found 1149 clusters that were unique to the Magnaporthaceae family of fungi with 295 of those containing genes from all three species. Gene clusters involved in metabolic and enzymatic activities were highly represented in the Magnaporthaceae specific clusters. Also highly represented in the Magnaporthaceae specific clusters as well as in the species specific genes were transcriptional regulators. In addition, we examined the relationship between gene evolution and distance to repetitive elements found in the genome. No correlations between diversifying or purifying selection and distance to repetitive elements or an increased rate of evolution in secreted and small secreted proteins were observed., Conclusions: Taken together, these data show that at the genome level, there is no evidence to suggest multi-speed genome evolution or that proximity to repetitive elements play a role in diversification of genes.

Background: Ochroconis mirabilis, a recently introduced water-borne dematiaceous fungus, is occasionally isolated from human skin lesions and nails. We identified an isolate of O. mirabilis from a skin scraping with morphological and molecular studies. Its genome was then sequenced and analysed for genetic features related to classification and biological characteristics., Results: UM 578 was identified as O. mirabilis based on morphology and internal transcribed spacer (ITS)-based phylogeny. The 34.61 Mb assembled genome with 13,435 predicted genes showed less efficiency of this isolate in plant cell wall degradation. Results from the peptidase comparison analysis with reported keratin-degrading peptidases from dermatophytes suggest that UM 578 is very unlikely to be utilising these peptidases to survive in the host. Nevertheless, we have identified peptidases from M10A, M12A and S33 families that may allow UM 578 to invade its host via extracellular matrix and collagen degradation. Furthermore, the lipases in UM 578 may have a role in supporting the fungus in host invasion. This fungus has the potential ability to synthesise melanin via the 1,8-dihydroxynaphthalene (DHN)-melanin pathway and to produce mycotoxins. The mating ability of this fungus was also inspected in this study and a mating type gene containing alpha domain was identified. This fungus is likely to produce taurine that is required in osmoregulation. The expanded gene family encoding the taurine catabolism dioxygenase TauD/TdfA domain suggests the utilisation of taurine under sulfate starvation. The expanded glutathione-S-transferase domains and RTA1-like protein families indicate the selection of genes in UM 578 towards adaptation in hostile environments., Conclusions: The genomic analysis of O. mirabilis UM 578 provides a better understanding of fungal survival tactics in different habitats.

Xylona heveae has only been isolated as an endophyte of rubber trees. In an effort to understand the genetic basis of endophytism, we compared the genome contents of X. heveae and 36 other Ascomycota with diverse lifestyles and nutritional modes. We focused on genes that are known to be important in the host-fungus interaction interface and that presumably have a role in determining the lifestyle of a fungus. We used phylogenomic data to infer the higher-level phylogenetic position of the Xylonomycetes, and mined ITS sequences to explore its taxonomic and ecological diversity. The X. heveae genome contains a low number of enzymes needed for plant cell wall degradation, suggesting that Xylona is a highly adapted specialist and likely dependent on its host for survival. The reduced repertoire of carbohydrate active enzymes could reflect an adaptation to intercellulary growth and to the avoidance of the host's immune system, suggesting that Xylona has a strictly endophytic lifestyle. Phylogenomic data resolved the position of Xylonomycetes as sister to Lecanoromycetes and Eurotiomycetes and placed the beetle-endosymbiont Symbiotaphrina as a member of this class. ITS data revealed that Trinosporium is also part of the Xylonomycetes, extending the taxonomic and ecological diversity of this group., (Copyright © 2015 The British Mycological Society. All rights reserved.)

The Ascomycete Onygenales order embraces a diverse group of mammalian pathogens, including the yeast-forming dimorphic fungal pathogens Histoplasma capsulatum, Paracoccidioides spp. and Blastomyces dermatitidis, the dermatophytes Microsporum spp. and Trichopyton spp., the spherule-forming dimorphic fungal pathogens in the genus Coccidioides, and many nonpathogens. Although genomes for all of the aforementioned pathogenic species are available, only one nonpathogen had been sequenced. Here, we enhance comparative phylogenomics in Onygenales by adding genomes for Amauroascus mutatus, Amauroascus niger, Byssoonygena ceratinophila, and Chrysosporium queenslandicum--four nonpathogenic Onygenales species, all of which are more closely related to Coccidioides spp. than any other known Onygenales species. Phylogenomic detection of gene family expansion and contraction can provide clues to fungal function but is sensitive to taxon sampling. By adding additional nonpathogens, we show that LysM domain-containing proteins, previously thought to be expanding in some Onygenales, are contracting in the Coccidioides-Uncinocarpus clade, as are the self-nonself recognition Het loci. The denser genome sampling presented here highlights nearly 800 genes unique to Coccidiodes, which have significantly fewer known protein domains and show increased expression in the endosporulating spherule, the parasitic phase unique to Coccidioides spp. These genomes provide insight to gene family expansion/contraction and patterns of individual gene gain/loss in this diverse order--both major drivers of evolutionary change. Our results suggest that gene family expansion/contraction can lead to adaptive radiations that create taxonomic orders, while individual gene gain/loss likely plays a more significant role in branch-specific phenotypic changes that lead to adaptation for species or genera., (Copyright © 2016 Whiston and Taylor.)

Magnaporthaceae is a family of ascomycetes that includes three fungi of great economic importance: Magnaporthe oryzae, Gaeumannomyces graminis var. tritici, and Magnaporthe poae. These three fungi cause widespread disease and loss in cereal and grass crops, including rice blast disease (M. oryzae), take-all disease in wheat and other grasses (G. graminis), and summer patch disease in turf grasses (M. poae). Here, we present the finished genome sequence for M. oryzae and draft sequences for M. poae and G. graminis var. tritici. We used multiple technologies to sequence and annotate the genomes of M. oryzae, M. poae, and G. graminis var. tritici. The M. oryzae genome is now finished to seven chromosomes whereas M. poae and G. graminis var. tritici are sequenced to 40.0× and 25.0× coverage respectively. Gene models were developed by the use of multiple computational techniques and further supported by RNAseq data. In addition, we performed preliminary analysis of genome architecture and repetitive element DNA., (Copyright © 2015 Okagaki et al.)

As part of a genome sequencing project for Ophiocordyceps sinensis, strain 1229, a complete mitochondrial (mt) genome was assembled as a single circular dsDNA of 157,510 bp, one of the largest reported for fungi. Conserved genes including the large and small rRNA subunits, 27 tRNA and 15 protein-coding genes, were identified. In addition, 58 non-conserved open reading frames (ncORFs) in the intergenic and intronic regions were also identified. Transcription analyses using RNA-Seq validated the expression of most conserved genes and ncORFs. Fifty-two introns (groups I and II) were found within conserved genes, accounting for 68.5% of the genome. Thirty-two homing endonucleases (HEs) with motif patterns LAGLIDADG (21) and GIY-YIG (11) were identified in group I introns. The ncORFs found in group II introns mostly encoded reverse transcriptases (RTs). As in other hypocrealean fungi, gene contents and order were found to be conserved in the mt genome of O. sinensis, but the genome size was enlarged by longer intergenic regions and numerous introns. Intergenic and intronic regions were composed of abundant repetitive sequences usually associated with mobile elements. It is likely that intronic ncORFs, which encode RTs and HEs, may have contributed to the enlarged mt genome of O. sinensis.

A number of sap-sucking insects harbor endosymbionts, which are thought to play an important role in the development of their hosts. One of the most important rice pests, the brown planthopper (BPH), Nilaparvata lugens (Stål), harbors an obligatory yeast-like symbiont (YLS) that cannot be cultured in vitro. Genomic information on this YLS would be useful to better understand its evolution. In this study, we performed genome sequencing of the YLS using both 454 and Illumina approaches, generating a draft genome that shows a slightly smaller genome size and relatively higher GC content than most ascomycete fungi. A phylogenomic analysis of the YLS supported its close relationship with insect pathogens. We analyzed YLS-specific genes and the categories of genes that are likely to have changed in the YLS during its evolution. The loss of mating type locus demonstrated in the YLS sheds light on the evolution of eukaryotic symbionts. This information about the YLS genome provides a helpful guide for further understanding endosymbiotic associations in hemiptera and the symbiotic replacement of ancient bacteria with a multifunctional YLS seems to have been a successful change., (© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Small non-coding RNAs constitute a large family of regulatory molecules with diverse function. Knowledge about small RNAs in fungi remains limited. Here, a RNA-seq-based analysis of small RNA transcriptome from the wheat pathogenic fungus Zymoseptoria tritici under nutrient-stress and nutrient-rich growth conditions is presented. Small RNAs mapped to Z. tritici genomic features including sense and antisense strands of genes, repetitive elements, rRNAs, snRNA, snoRNA and tRNA loci throughout the chromosomes. The majority of small RNAs from rich condition mapped to all the chromosomes, primarily the sense strand of chromosome seven, but under stress condition, small RNAs were enriched in accessory genome. Inspection of the predominant small RNAs matching to protein-coding genes, repeats and rRNAs revealed that intron and LTP retrotransposon-derived small RNAs accumulated upon stress, whereas rRNA-derived small RNAs were enriched for nutrient-rich condition. U enrichment at position 1 of small RNAs mapping to repeats and antisense-strand genes was only observed from stress library. Additionally, target analysis of one predicted milRNA suggests its potential role in fungal pathogenesis and manipulation of host immunity. Taken together, the findings provide a new insight into the function of small RNAs in regulating growth and development, stress adaptation, and possibly pathogenesis in Z. tritici., (Copyright © 2015 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.)

Background: Pseudogymnoascus spp. is a wide group of fungi lineages in the family Pseudorotiaceae including an aggressive pathogen of bats P. destructans. Although several lineages of P. spp. were shown to produce ascospores in culture, the vast majority of P. spp. demonstrates no evidence of sexual reproduction. P. spp. can tolerate a wide range of different temperatures and salinities and can survive even in permafrost layer. Adaptability of P. spp. to different environments is accompanied by extremely variable morphology and physiology., Results: We sequenced genotypes of 14 strains of P. spp., 5 of which were extracted from permafrost, 1 from a cryopeg, a layer of unfrozen ground in permafrost, and 8 from temperate surface environments. All sequenced genotypes are haploid. Nucleotide diversity among these genomes is very high, with a typical evolutionary distance at synonymous sites dS ≈ 0.5, suggesting that the last common ancestor of these strains lived >50 Mya. The strains extracted from permafrost do not form a separate clade. Instead, each permafrost strain has close relatives from temperate environments. We observed a strictly clonal population structure with no conflicting topologies for ~99% of genome sequences. However, there is a number of short (~100-10,000 nt) genomic segments with the total length of 67.6 Kb which possess phylogenetic patterns strikingly different from the rest of the genome. The most remarkable case is a MAT-locus, which has 2 distinct alleles interspersed along the whole-genome phylogenetic tree., Conclusions: Predominantly clonal structure of genome sequences is consistent with the observations that sexual reproduction is rare in P. spp. Small number of regions with noncanonical phylogenies seem to arise due to some recombination events between derived lineages of P. spp., with MAT-locus being transferred on multiple occasions. All sequenced strains have heterothallic configuration of MAT-locus.

Exserohilum rostratum was the cause of most cases of fungal meningitis and other infections associated with the injection of contaminated methylprednisolone acetate produced by the New England Compounding Center (NECC). Until this outbreak, very few human cases of Exserohilum infection had been reported, and very little was known about this dematiaceous fungus, which usually infects plants. Here, we report using whole-genome sequencing (WGS) for the detection of single nucleotide polymorphisms (SNPs) and phylogenetic analysis to investigate the molecular origin of the outbreak using 22 isolates of E. rostratum retrieved from 19 case patients with meningitis or epidural/spinal abscesses, 6 isolates from contaminated NECC vials, and 7 isolates unrelated to the outbreak. Our analysis indicates that all 28 isolates associated with the outbreak had nearly identical genomes of 33.8 Mb. A total of 8 SNPs were detected among the outbreak genomes, with no more than 2 SNPs separating any 2 of the 28 genomes. The outbreak genomes were separated from the next most closely related control strain by ∼136,000 SNPs. We also observed significant genomic variability among strains unrelated to the outbreak, which may suggest the possibility of cryptic speciation in E. rostratum., (Copyright © 2014, American Society for Microbiology. All Rights Reserved.)