A genome-informed higher rank classification of the biotechnologically important fungal subphylum Saccharomycotina.

The subphylum Saccharomycotina is a lineage in the fungal phylum Ascomycota that exhibits levels of genomic diversity similar to those of plants and animals. The Saccharomycotina consist of more than 1 200 known species currently divided into 16 families, one order, and one class. Species in this subphylum are ecologically and metabolically diverse and include important opportunistic human pathogens, as well as species important in biotechnological applications. Many traits of biotechnological interest are found in closely related species and often restricted to single phylogenetic clades. However, the biotechnological potential of most yeast species remains unexplored. Although the subphylum Saccharomycotina has much higher rates of genome sequence evolution than its sister subphylum, Pezizomycotina, it contains only one class compared to the 16 classes in Pezizomycotina. The third subphylum of Ascomycota, the Taphrinomycotina, consists of six classes and has approximately 10 times fewer species than the Saccharomycotina. These data indicate that the current classification of all these yeasts into a single class and a single order is an underappreciation of their diversity. Our previous genome-scale phylogenetic analyses showed that the Saccharomycotina contains 12 major and robustly supported phylogenetic clades; seven of these are current families (Lipomycetaceae, Trigonopsidaceae, Alloascoideaceae, Pichiaceae, Phaffomycetaceae, Saccharomycodaceae, and Saccharomycetaceae), one comprises two current families (Dipodascaceae and Trichomonascaceae), one represents the genus Sporopachydermia, and three represent lineages that differ in their translation of the CUG codon (CUG-Ala, CUG-Ser1, and CUG-Ser2). Using these analyses in combination with relative evolutionary divergence and genome content analyses, we propose an updated classification for the Saccharomycotina, including seven classes and 12 orders that can be diagnosed by genome content. This updated classification is consistent with the high levels of genomic diversity within this subphylum and is necessary to make the higher rank classification of the Saccharomycotina more comparable to that of other fungi, as well as to communicate efficiently on lineages that are not yet formally named. Taxonomic novelties: New classes: Alloascoideomycetes M. Groenew., Hittinger, Opulente & A. Rokas, Dipodascomycetes M. Groenew., Hittinger, Opulente & A. Rokas, Lipomycetes M. Groenew., Hittinger, Opulente, A. Rokas, Pichiomycetes M. Groenew., Hittinger, Opulente & A. Rokas, Sporopachydermiomycetes M. Groenew., Hittinger, Opulente & A. Rokas, Trigonopsidomycetes M. Groenew., Hittinger, Opulente & A. Rokas. New orders: Alloascoideomycetes: Alloascoideales M. Groenew., Hittinger, Opulente & A. Rokas; Dipodascomycetes: Dipodascales M. Groenew., Hittinger, Opulente & A. Rokas; Lipomycetes: Lipomycetales M. Groenew., Hittinger, Opulente & A. Rokas; Pichiomycetes: Alaninales M. Groenew., Hittinger, Opulente & A. Rokas, Pichiales M. Groenew., Hittinger, Opulente & A. Rokas, Serinales M. Groenew., Hittinger, Opulente & A. Rokas; Saccharomycetes: Phaffomycetales M. Groenew., Hittinger, Opulente & A. Rokas, Saccharomycodales M. Groenew., Hittinger, Opulente & A. Rokas; Sporopachydermiomycetes: Sporopachydermiales M. Groenew., Hittinger, Opulente & A. Rokas; Trigonopsidomycetes: Trigonopsidales M. Groenew., Hittinger, Opulente & A. Rokas. New families: Alaninales: Pachysolenaceae M. Groenew., Hittinger, Opulente & A. Rokas; Pichiales: Pichiaceae M. Groenew., Hittinger, Opulente & A. Rokas; Sporopachydermiales: Sporopachydermiaceae M. Groenew., Hittinger, Opulente & A. Rokas. Citation: Groenewald M, Hittinger CT, Bensch K, Opulente DA, Shen X-X, Li Y, Liu C, LaBella AL, Zhou X, Limtong S, Jindamorakot S, Gonçalves P, Robert V, Wolfe KH, Rosa CA, Boekhout T, Cadez N, Péter G, Sampaio JP, Lachance M-A, Yurkov AM, Daniel H-M, Takashima M, Boundy-Mills K, Libkind D, Aoki K, Sugita T, Rokas A (2023). A genome-informed higher rank classification of the biotechnologically important fungal subphylum Saccharomycotina. Studies in Mycology 105: 1-22. doi: 10.3114/sim.2023.105.01 This study is dedicated to the memory of Cletus P. Kurtzman (1938-2017), a pioneer of yeast taxonomy.

Into the wild: new yeast genomes from natural environments and new tools for their analysis.

Genomic studies of yeasts from the wild have increased considerably in the past few years. This revolution has been fueled by advances in high-throughput sequencing technologies and a better understanding of yeast ecology and phylogeography, especially for biotechnologically important species. The present review aims to first introduce new bioinformatic tools available for the generation and analysis of yeast genomes. We also assess the accumulated genomic data of wild isolates of industrially relevant species, such as Saccharomyces spp., which provide unique opportunities to further investigate the domestication processes associated with the fermentation industry and opportunistic pathogenesis. The availability of genome sequences of other less conventional yeasts obtained from the wild has also increased substantially, including representatives of the phyla Ascomycota (e.g. Hanseniaspora) and Basidiomycota (e.g. Phaffia). Here, we review salient examples of both fundamental and applied research that demonstrate the importance of continuing to sequence and analyze genomes of wild yeasts.

Diversity of Secondary Metabolism in Aspergillus nidulans Clinical Isolates.

The filamentous fungus Aspergillus nidulans has been a primary workhorse used to understand fungal genetics. Much of this work has focused on elucidating the genetics of biosynthetic gene clusters (BGCs) and the secondary metabolites (SMs) they produce. SMs are both niche defining in fungi and of great economic importance to humans. Despite the focus on A. nidulans, very little is known about the natural diversity in secondary metabolism within this species. We determined the BGC content and looked for evolutionary patterns in BGCs from whole-genome sequences of two clinical isolates and the A4 reference genome of A. nidulans Differences in BGC content were used to explain SM profiles determined using liquid chromatography-high-resolution mass spectrometry. We found that in addition to genetic variation of BGCs contained by all isolates, nine BGCs varied by presence/absence. We discovered the viridicatumtoxin BGC in A. nidulans and suggest that this BGC has undergone a horizontal gene transfer from the Aspergillus section Nigri lineage into Penicillium sometime after the sections Nigri and Nidulantes diverged. We identified the production of viridicatumtoxin and several other compounds previously not known to be produced by A. nidulans One isolate showed a lack of sterigmatocystin production even though it contained an apparently intact sterigmatocystin BGC, raising questions about other genes and processes known to regulate this BGC. Altogether, our work uncovers a large degree of intraspecies diversity in BGC and SM production in this genetic model species and offers new avenues to understand the evolution and regulation of secondary metabolism.IMPORTANCE Much of what we know about the genetics underlying secondary metabolite (SM) production and the function of SMs in the model fungus Aspergillus nidulans comes from a single reference genome. A growing body of research indicates the importance of biosynthetic gene cluster (BGC) and SM diversity within a species. However, there is no information about the natural diversity of secondary metabolism in A. nidulans We discovered six novel clusters that contribute to the considerable variation in both BGC content and SM production within A. nidulans We characterize a diverse set of mutations and emphasize how findings of single nucleotide polymorphisms (SNPs), deletions, and differences in evolutionary history encompass much of the variation observed in nonmodel systems. Our results emphasize that A. nidulans may also be a strong model to use within-species diversity to elucidate regulatory cross talk, fungal ecology, and drug discovery systems.

Evidence for a high frequency of simultaneous double-nucleotide substitutions.

Point mutations are generally assumed to involve changes of single nucleotides. Nevertheless, the nature and known mechanisms of mutation do not exclude the possibility that several adjacent nucleotides may change simultaneously in a single mutational event. Two independent approaches are used here to estimate the frequency of simultaneous double-nucleotide substitutions. The first examines switches between TCN and AGY (where N is any nucleotide and Y is a pyrimidine) codons encoding absolutely conserved serine residues in a number of proteins from diverse organisms. The second reveals double-nucleotide substitutions in primate noncoding sequences. These two complementary approaches provide similar high estimates for the rate of doublet substitutions, on the order of 0.1 per site per billion years.

What can comparative genomics tell us about species concepts in the genus Aspergillus?

Understanding the nature of species" boundaries is a fundamental question in evolutionary biology. The availability of genomes from several species of the genus Aspergillus allows us for the first time to examine the demarcation of fungal species at the whole-genome level. Here, we examine four case studies, two of which involve intraspecific comparisons, whereas the other two deal with interspecific genomic comparisons between closely related species. These four comparisons reveal significant variation in the nature of species boundaries across Aspergillus. For example, comparisons between A. fumigatus and Neosartorya fischeri (the teleomorph of A. fischerianus) and between A. oryzae and A. flavus suggest that measures of sequence similarity and species-specific genes are significantly higher for the A. fumigatus - N. fischeri pair. Importantly, the values obtained from the comparison between A. oryzae and A. flavus are remarkably similar to those obtained from an intra-specific comparison of A. fumigatus strains, giving support to the proposal that A. oryzae represents a distinct ecotype of A. flavus and not a distinct species. We argue that genomic data can aid Aspergillus taxonomy by serving as a source of novel and unprecedented amounts of comparative data, as a resource for the development of additional diagnostic tools, and finally as a knowledge database about the biological differences between strains and species.

Intron insertion as a phylogenetic character: the engrailed homeobox of Strepsiptera does not indicate affinity with Diptera.

The phylogenetic relationships of the order Strepsiptera are unclear. Affiliation to Coleoptera has been proposed, however this implies that dipteran halteres and strep-sipteran haltere-like organs evolved convergently. An alternative is a sister group relationship with Diptera. In this case, halteres could be homologous but a radical homeotic mutation may have switched their position to the Strepsipteran mesothorax. Ribosomal DNA sequence analysis has been used to support Dipteran affiliation, although this is controversial. Here we investigate the potential of an intron insertion site as a phylogenetic character. We find that the en homeobox gene of the strepsipteran Stichotrema dallatorreanum lacks a derived intron insertion shared by representatives of Diptera and Lepidoptera. We argue against a close affiliation between Strepsiptera and Diptera.

Rare genomic changes as a tool for phylogenetics.

DNA sequence data have offered valuable insights into the relationships between living organisms. However, most phylogenetic analyses of DNA sequences rely primarily on single nucleotide substitutions, which might not be perfect phylogenetic markers. Rare genomic changes (RGCs), such as intron indels, retroposon integrations, signature sequences, mitochondrial and chloroplast gene order changes, gene duplications and genetic code changes, provide a suite of complementary markers with enormous potential for molecular systematics. Recent exploitation of RGCs has already started to yield exciting phylogenetic information.

Understanding patterns of genetic diversity in the oak gallwasp Biorhiza pallida: demographic history or a Wolbachia selective sweep?

The endosymbiont Wolbachia can be responsible for selective sweeps on mitochondrial DNA variability within species. Similar signals can also result from demographic processes, although crucially the latter affect nuclear as well as mitochondrial loci. Here we present data on Wolbachia infection status and phylogeographic patterning for a widely distributed insect host, the oak gallwasp Biorhiza pallida (Hymenoptera: Cynipidae). Two hundred and eighteen females from eight European countries were screened for Wolbachia. All individuals from Hungary, Italy, France, U.K., Ireland, Switzerland, Sweden, and northern and southern Spain were infected with a single group A strain of Wolbachia, while populations in central Spain were not infected. A mitochondrial marker (cytochrome b) shows low variation and departure from neutrality in infected populations, but greater variation and no deviation from neutrality in Wolbachia-free populations. This pattern is compatible with a Wolbachia-induced selective sweep. However, we also find parallel differences between infected and uninfected populations for nuclear markers (sequence data for ITS1 and ITS2). All markers support the existence of a deep split between populations in Spain (some free of Wolbachia), and those in the rest of Europe (all infected). Allelic variation for five allozyme loci is also consistent with the Spain-rest of Europe split. Concordant patterns for nuclear and mitochondrial markers suggest that differences in the nature and extent of genetic diversity between these two regions are best explained by differing demographic histories (perhaps associated with range expansion from Pleistocene glacial refugia), rather than a Wolbachia-associated selective sweep.

Differential success in northwards range expansion between ecotypes of the marble gallwasp Andricus kollari: a tale of two lifecycles.

The Marble gallwasp Andricus kollari has a native range divided into two geographically separated lifecycles. In Eastern Europe and Turkey, the lifecycle involves a sexual generation on Turkey oak, Quercus cerris, while in Iberia and North Africa the sexual generation host is cork oak, Q. suber. Over the last 500 years, A. kollari has expanded its range into northern Europe, following human planting of Q. cerris from Italy and the Balkans. We ask: (i) what is the genetic relationship between eastern and western distributions of Andricus kollari? Can we determine which lifecycle is ancestral, and how long ago they diverged? (ii) To what extent have eastern and western native ranges contributed to northwards range expansion? (iii) Is there any evidence for hybridization between the two life cycle types? We present analyses of allozyme data for 13 polymorphic loci and of sequence variation for a 433 bp fragment of the mitochondrial cytochrome b gene. These show: (i) that four haplotype lineages (one in Spain, two in Hungary/Italy and one in Turkey) diverged more or less simultaneously between 1 and 2 million years ago, suggesting the existence of at least four refuges through recent ice age cycles. Our data cannot resolve which lifecycle type is ancestral. (ii) Populations north of putative refuges are divided into two sets. Populations in south-west France are allied to Spain, while all remaining populations in northern Europe have been colonized from Italy and the Balkans. (iii) The transition from one race to another in south-west France is marked by abrupt transitions in the frequency of refuge-specific private alleles and corresponds closely to the northern limit of the distribution of cork oak. Although hybrids were detected in north-west France, none were detected where the two lifecycles meet in south-western France. The biology of oak gallwasps predicts that any hybrid zone will be narrow, and limited to regions where Q. cerris and Q. suber meet. Our data suggest that eastern and western A. kollari are effectively separate species.

The incidence and diversity of Wolbachia in gallwasps (Hymenoptera; Cynipidae) on oak.

Wolbachia bacteria infect approximately 20% of all insect species, and cause a range of alterations to host reproduction, including imposition of thelytoky. The incidence and phenotypic impact of Wolbachia remains to be established in many insect taxa, and considerable research effort is currently focused on its association with particular reproductive modes and the relative importance of the various pathways via which infection occurs. Gallwasps represent an attractive system for addressing these issues for two reasons. First, they show a diversity of reproductive modes (including arrhenotoky, thelytoky and cyclical parthenogenesis) in which the impact of Wolbachia infection can be examined. Second, they occupy two intimately linked trophic niches (gall-inducers and inquilines) between which there is potential for the horizontal exchange of Wolbachia infection. In the arrhenotokous gallwasp lineages screened to date (the herb-galling 'Aylacini' and the rose-galling Diplolepidini), Wolbachia infection always induces thelytoky. The impact of Wolbachia in other arrhenotokous clades, and in the cyclically parthenogenetic clades remains unknown. Here we use polymerase chain reaction (PCR) screening and sequence data for two Wolbachia genes (wsp and ftsZ) to examine the prevalence and incidence of Wolbachia infection in 64 species (a total of 609 individuals) in two further tribes: the arrhenotokous inquilines (tribe Synergini), and the cyclically parthenogenetic oak gallwasps (tribe Cynipini). We ask: (i) whether Wolbachia infection has any apparent impact on host reproduction in the two tribes and (ii) whether there is any correlation between Wolbachia infection and the apparent lack of an arrhenotokous generation in many oak gallwasp life cycles. We show: (i) that Wolbachia infection is rare in the Cynipini. Infected species show no deviation from cyclical parthenogenesis, and infection is no more common in species known only from a thelytokous generation; (ii) that there is a higher incidence of infection within the arrhenotokous inquilines, and generally in gallwasp tribes without cyclical parthenogensis; (iii) all Wolbachia-positive inquiline species are known to possess males, implying either that Wolbachia infection does not result in loss of sex in this tribe or, more probably, that (as for some rose gallwasps) Wolbachia infection leads to loss of sex in specific populations; and (iv) although we find some inquilines and gall inducers to be infected with Wolbachia having the same wsp sequence, these hosts are not members of the same gall communities, arguing against frequent horizontal transmission between these two trophic groups. We suggest that exchange may be mediated by the generalist parasitoids common in oak galls.