What an Intron May Tell: Several Sexual Biospecies Coexist in Meriderma spp. (Myxomycetes)

A new species of Diacheopsis (Myxomycetes) and a new habitat for myxomycetes

We describe a new species, Diacheopsis resinae (Myxomycetes), collected from a microhabitat new for myxomycetes: stem wounds of coniferous trees (Norway spruce) where the resin is overgrown with a community of resinicolous fungi. The 80 known collections come from the Vosges (France), the Black Forest (Germany), Swabian Alp (Germany), and several localities in Denmark and Norway. Observations, but as well as metabarcoding of substrate samples with fungal (ITS [internal transcribed spacer]), bacterial (16S rDNA), and myxomycete (18S nuc rDNA) primers from eight trunks, revealed the new myxomycete to co-occur with resin-degrading ascomycetes (Infundichalara microchona, Lophium arboricola, Zythia resinae). The gram-negative bacterial genera Endobacter and Sphingomonas were found to be abundant in the substrate and may be a food source for the myxomycete. Fruit bodies were found mostly during the more humid winter season, with a peak in January/February. Partial sequences of two independent molecular markers (18S nuc rDNA, EF1α [elongation factor 1-alpha] gene) were obtained for 41 accessions, which form a monophyletic cluster in a two-gene phylogeny of Stemonititidales but do not group with other species of Diacheopsis, thus rendering this genus paraphyletic. The new species, although exclusively developing sessile sporocarps and morphologically undoubtedly falling into the genus Diacheopsis, is most closely related to species of Lamproderma, especially L. album, L. zonatum, and L. zonatopulchellum. Within D. resinae, three groups can be differentiated, which show nearly complete reproductive isolation, as judged from a recombination analysis of the two unlinked markers and the allelic combinations of the EF1α gene.

The Systematics and Phylogeny of Myxomycetes: Yesterday, Today, and Tomorrow

Myxomycetes are amoeboid fungus-like organisms (Amoebozoa) with a unique life cycle characterized by a great morphological diversity of fruiting bodies. Due to the similarity of these structures to the fruiting bodies of some representatives of Ascomycota and Basidiomycota, myxomycetes have been classified as fungi since the first known scientific description in 1654. Only in the 19th century, when their life cycle was studied, did the difference of this group from fungi become clear. During the same period, microscopic structures of fruiting bodies, as well as ornamentation of the spore surface, began to be considered as diagnostic features. Due to this, in the period from the end of the 19th to the middle of the 20th century, a rather stable system was formed. However, as further studies have shown, both macro- and micromorphological characters are often quite variable, depend on environmental conditions, and often result from convergent evolution, which causes difficulties in defining species and taxonomic units of higher ranks. Since the first decade of the 21st century, due to the development of molecular genetic methods and the accumulation of data on nucleotide sequences of marker genes together with the improvement of microscopic studies, it has been possible to obtain data on the evolutionary relationships of different groups of myxomycetes. A milestone in this process was the publication of the first phylogenetic system of myxomycetes in 2019. This work was the starting point for a number of studies on the relationships between different groups of myxomycetes at a lower taxonomic level. Thus, there has been a surge in the number of studies that bring us closer to constructing a natural system.

A new species of the Trichia botrytis complex from the cloud forest in Ecuador

The paper describes a new myxomycete species, Trichia tuberculata, from the Ecuadorian cloud forest. The phylogeny constructed with nuclear 18S rDNA and mitochondrial 17S rDNA sequences indicates that the taxon is closely related to recently described species T. acetocorticola, T. nubila, and T. pinicola, all within the T. botrytis complex. The morphology of the inner peridium surface and capillitium also aligns T. tuberculata with these species. However, the distinctive tuberculate peridium, short elater tips, and perhaps the distribution range (so far found in the equatorial zone only) distinguish T. tuberculata from related taxa.

DNA barcodes reliably differentiate between nivicolous species of Diderma (Myxomycetes, Amoebozoa) and reveal regional differences within Eurasia

The nivicolous species of the genus Diderma are challenging to identify, and there are several competing views on their delimitation. We analyzed 102 accessions of nivicolous Diderma spp. that were sequenced for two or three unlinked genes to determine which of the current taxonomic treatments is better supported by molecular species delimitation methods. The results of a haplotype web analysis, Bayesian species delimitation under a multispecies coalescent model, and phylogenetic analyses on concatenated alignments support a splitting approach that distinguishes six taxa: Diderma alpinum, D. europaeum, D. kamchaticum, D. meyerae, D. microcarpum and D. niveum. The first two approaches also support the separation of Diderma alpinum into two species with allopatric distribution. An extended dataset of 800 specimens (mainly from Europe) that were barcoded with 18S rDNA revealed only barcode variants similar to those in the species characterized by the first data set, and showed an uneven distribution of these species in the Northern Hemisphere: Diderma microcarpum and D. alpinum were the only species found in all seven intensively sampled mountain regions. Partial 18S rDNA sequences serving as DNA barcodes provided clear signatures that allowed for unambiguous identification of the nivicolous Diderma spp., including two putative species in D. alpinum.

Tasmaniomyxa umbilicata, a new genus and new species of myxomycete from Tasmania

A new genus and species of myxomycete, Tasmaniomyxa umbilicata, is described based on numerous observations in Tasmania and additional records from southeastern Australia and New Zealand. The new taxon is characterized by an unusual combination of characters from two families: Lamprodermataceae and Didymiaceae. With Lamprodermataceae the species shares limeless sporocarps, a shining membranous peridium, an epihypothallic stalk, and a cylindrical columella. Like Didymiaceae, it has a soft, flaccid, sparsely branched capillitium, with rough tubular threads that contain fusiform nodes and are firmly connected to the peridium. Other characters of T. umbilicata that also occur in many Didymiaceae are the peridium dehiscing into petaloid lobes, the yellow, motile plasmodium, and the spores ornamented with larger, grouped and smaller, scattered warts. The transitional position of the new taxon is reflected by a three-gene phylogeny, which places T. umbilicata at the base of the branch of all lime-containing Physarales, thus justifying its description as a monotypic genus.

High Environmentally Induced Plasticity in Spore Size and Numbers of Nuclei per Spore in Physarum albescens (Myxomycetes)

Spore size enables dispersal in plasmodial slime molds (Myxomycetes) and is an important taxonomic character. We recorded size and the number of nuclei per spore for 39 specimens (colonies of 50-1000 sporocarps) of the nivicolous myxomycete Physarum albescens, a morphologically defined taxon with several biological species. For each colony, three sporocarps were analyzed from the same spore mount under brightfield and DAPI-fluorescence, recording ca. 14,000 spores per item. Diagrams for spore size distribution showed narrow peaks of mostly uninucleate spores. Size was highly variable within morphospecies (10.6-13.5 µm, 11-13%), biospecies (3-13%), even within spatially separated colonies of one clone (ca. 8%); but fairly constant for a colony (mean variation 0.4 µm, ca. 1.5%). ANOVA explains most of this variation by the factor locality (within all colonies: 32.7%; within a region: 21.4%), less by biospecies (13.5%), whereas the contribution of intra-colony variation was negligible (<0.1%). Two rare aberrations occur: 1) multinucleate spores and 2) oversized spores with a double or triple volume of normal spores. Both are not related to each other or limited to certain biospecies. Spore size shows high phenotypic plasticity, but the low variation within a colony points to a strong genetic background.

Phenotypic plasticity in plasmodial slime molds and molecular phylogeny of terrestrial vs. aquatic species

Fifty years ago, the enigmatic Brazilian myxomycete-species Didymium aquatile was described and analyzed with respect to the structure of the plasmodium and its spores. In this study, we compare this rare plasmodial slime mold with another, temporarily aquatic taxon from Europe, Didymium nigripes. Phenotypic plasticity of D. nigripes was investigated under various environmental conditions. Large changes in the morphology of the plasmodia were observed. For species identification, characteristics of the fruiting bodies are key features. However, Didymium aquatile was only characterized by its “abnormal” plasmodia, but no molecular data were available. Here, we analyzed DNA-sequences of 22 species of the genera Didymium and Diderma with a focus on this South American taxon via molecular genetics. A comparison of 18S-rDNA-sequences from D. aquatile and 21 other Didymium (and Diderma)-species indicates that D. aquatile is a reproductively isolated morpho-species. Phenotypic plasticity of D. nigripes is documented with respect to plasmodium morphology and the formation of fruiting bodies, as an example of an adaptation of a terrestrial species to aquatic environments.

A workflow for low-cost automated image analysis of myxomycete spore numbers, size and shape

Measuring spore size is a standard method for the description of fungal taxa, but in manual microscopic analyses the number of spores that can be measured and information on their morphological traits are typically limited. To overcome this weakness we present a method to analyze the size and shape of large numbers of spherical bodies, such as spores or pollen, by using inexpensive equipment. A spore suspension mounted on a slide is treated with a low-cost, high-vibration device to distribute spores uniformly in a single layer without overlap. Subsequently, 10,000 to 50,000 objects per slide are measured by automated image analysis. The workflow involves (1) slide preparation, (2) automated image acquisition by light microscopy, (3) filtering to separate high-density clusters, (4) image segmentation by applying a machine learning software, Waikato Environment for Knowledge Analysis (WEKA), and (5) statistical evaluation of the results. The technique produced consistent results and compared favorably with manual measurements in terms of precision. Moreover, measuring spore size distribution yields information not obtained by manual microscopic analyses, as shown for the myxomycete Physarum albescens. The exact size distribution of spores revealed irregularities in spore formation resulting from the influence of environmental conditions on spore maturation. A comparison of the spore size distribution within and between sporocarp colonies showed large environmental and likely genetic variation. In addition, the comparison identified specimens with spores roughly twice the normal size. The successful implementation of the presented method for analyzing myxomycete spores also suggests potential for other applications.

Genetic structure of the protist Physarum albescens (Amoebozoa) revealed by multiple markers and genotyping by sequencing

Myxomycetes are terrestrial protists with many presumably cosmopolitan species dispersing via airborne spores. A truly cosmopolitan species would suffer from outbreeding depression hampering local adaptation, while locally adapted species with limited distribution would be at a higher risk of extinction in changing environments. Here, we investigate intraspecific genetic diversity and phylogeography of Physarum albescens over the entire Northern Hemisphere. We sequenced 324 field collections of fruit bodies for 1-3 genetic markers (SSU, EF1A, COI) and analysed 98 specimens with genotyping by sequencing. The structure of the three-gene phylogeny, SNP-based phylogeny, phylogenetic networks, and the observed recombination pattern of three independently inherited gene markers can be best explained by the presence of at least 18 reproductively isolated groups, which can be seen as cryptic species. In all intensively sampled regions and in many localities, members of several phylogroups coexisted. Some phylogroups were found to be abundant in only one region and completely absent in other well-studied regions, and thus may represent regional endemics. Our results demonstrate that the widely distributed myxomycete species Ph. albescens represents a complex of at least 18 cryptic species, and some of these seem to have a limited geographical distribution. In addition, the presence of groups of presumably clonal specimens suggests that sexual and asexual reproduction coexist in natural populations of myxomycetes.

Didymium pseudonivicola: A new myxomycete from the austral Andes emerges from broad-scale morphological and molecular analyses of D. nivicola collections

A new nivicolous myxomycete is described as a result of a comprehensive study of Didymium nivicola collections from the entire range of its occurrence. Statistical analysis of 12 morphological characters, phylogenetic analyses of nuc 18S rDNA and elongation factor 1-alpha gene (EF1A), and a delimitation method (automatic barcode gap diversity) have been applied to corroborate the identity of the new species. A preliminary morphological analysis of D. nivicola revealed high variability of South American populations where four types of spore ornamentation were noted. However, results of molecular study and statistical analysis of morphological characters did not support recognition of these four forms but the distinction of two morphotypes. Consequently, two species have been recognized: D. nivicola and the newly proposed D. pseudonivicola. The new species can be distinguished from D. nivicola by distinctly larger and mostly plasmodiocarpic sporophores, which are scattered to gregarious, paler spores, and by the paler, more delicate and more elastic capillitium. Spore ornamentation of D. pseudonivicola is uniform and can be described as distinctly spiny (pilate under scanning electron microscope [SEM]), whereas those of D. nivicola is more variable, where spines (pilae under SEM) are delicate, distinct, or conspicuous. Additionally, whereas D. nivicola is a species distributed worldwide, D. pseudonivicola occurs only in the austral Andes of Argentina and Chile.

Design of potentially universal SSU primers in myxomycetes using next-generation sequencing

Unlike fungi, which have a universally accepted barcode marker, universal primers still lack in myxomycetes. Typically, DNA barcode primers were designed based on comparing existing myxomycetes sequences and targeting the conserved regions. However, the extreme genetic diversity within major myxomycetes groups and the frequent occurrence of group I introns have made the development of universal DNA barcode a severe challenge. The emergence of next-generation sequencing provides an opportunity to address this problem. We sequenced the mixed genomic DNA of 81 myxomycetes and extracted the SSU gene's reads using next-generation sequencing. After alignment and assembly, we designed a set of SSU primers that matched all potential SNPs, avoided all known group I intron insertion sites, and were highly conserved between major myxomycetes orders. This set of SSU primers has the potential to become one of the universal primer combinations. Due to the high genetic divergence caused by long and complicated evolutionary histories, the lack of universal barcode primers is common in protists. Our research provides a new method to solve this problem.

Range-wide Phylogeography of a Nivicolous Protist Didymium nivicola Meyl. (Myxomycetes, Amoebozoa): Striking Contrasts Between the Northern and the Southern Hemisphere

Soil protists play a crucial role in terrestrial ecosystems and often show immense taxonomic diversity. However, for many groups, distribution patterns remain largely unknown. We investigated range-wide intraspecific diversity of a specialized airborne protist (Didymium nivicola Meyl.) that occupies a narrow ecological niche associated with long-lasting snow cover. We sampled 122 collections covering all areas where the species was recorded worldwide. We obtained 105 and 41 sequences of small ribosomal subunit rDNA (SSU) and elongation factor 1-alpha (EF1A), respectively. While the species is very diverse in the austral Andes, Southern Hemisphere (SH; 17 SSU ribotypes and 12 EF1A genotypes identified), its populations are genetically uniform across three continents of the Northern Hemisphere (NH; single ribotype, single genotype). Our results indicate the austral Andes as a possible diversification centre for D. nivicola where populations seem to reproduce sexually. Two main parts of the range display highly contrasting genetic patterns, thus biogeographical history and dynamics. Current distribution of D. nivicola in the NH is likely a result of a dispersal event from the SH and subsequent long-distance dispersal (LDD) that might be associated with a shift to asexual mode of reproduction.

New protocol for successful isolation and amplification of DNA from exiguous fractions of specimens: a tool to overcome the basic obstacle in molecular analyses of myxomycetes

Herbarium collections provide an essential basis for a wide array of biological research and, with development of DNA-based methods, they have become an invaluable material for genetic analyses. Yet, the use of such material is hindered by technical limitations related to DNA degradation and to quantity of biological material. The latter is inherent for some biological groups, as best exemplified by myxomycetes which form minute sporophores. It is estimated that ca. two-thirds of myxomycete taxa are represented by extremely scanty material. As DNA isolation methods applied so far in myxomycete studies require destructive sampling of many sporophores, a large part of described diversity of the group remains unavailable for phylogenetic studies or barcoding. Here, we tested several procedures of DNA isolation and amplification to seek for an efficient and possibly non-destructive method of sampling. Tests were based on herbarium specimens of 19 species representing different taxonomic orders. We assayed several variants of isolation based on silica gel membrane columns, and a newly designed procedure using highly reduced amount of biological material (small portion of spores), based on fine disruption of spores and direct PCR. While the most frequently used column-based method led to PCR success in 89.5% of samples when a large amount of material was used, its performance dropped to 52% when based on single sporophores. Single sporophores provided amplicons in 89.5% of samples when using a kit dedicated to low-amount DNA samples. Our new procedure appeared the most effective (94.7%) while it used only a small fraction of spores, being nearly non-destructive; it was also the most cost-effective. We thus demonstrate that combination of adequate handling of spore micro-disruption coupled with application of direct PCR can be an efficient way to circumvent technical limitations for genetic studies in myxomycetes and thus can substantially improve taxon sampling for phylogeny and barcoding. Additionally, this approach gives a unique possibility to apply both molecular and morphological assays to the same structure (sporophore), which then can be further stored as documentation.

Different Degrees of Niche Differentiation for Bacteria, Fungi, and Myxomycetes Within an Elevational Transect in the German Alps

We used direct DNA amplification from soil extracts to analyze microbial communities from an elevational transect in the German Alps by parallel metabarcoding of bacteria (16S rRNA), fungi (ITS2), and myxomycetes (18S rRNA). For the three microbial groups, 5710, 6133, and 261 operational taxonomic units (OTU) were found. For the latter group, we can relate OTUs to barcodes from fruit bodies sampled over a 4-year period. The alpha diversity of myxomycetes was positively correlated with that of bacteria. Vegetation type was found to be the main explanatory parameter for the community composition of all three groups and a substantial species turnover with elevation was observed. Bacteria and fungi display similar community responses, driven by symbiont species and plant substrate quality. Myxamoebae show a more patchy distribution, though still clearly stratified between taxa, which seems to be a response to both structural properties of the habitat and interaction with specific bacterial and fungal taxa. Finally, we report a high number of myxomycete OTUs not represented in a reference database from fructifications, which might represent novel species.

A four year survey reveals a coherent pattern between occurrence of fruit bodies and soil amoebae populations for nivicolous myxomycetes

Among soil-inhabiting protists, myxomycetes stand out by their macroscopic fructifications which have allowed studies on their ecology and distribution for more than two hundred years. One of the most distinct ecological guilds in myxomycetes are the nivicolous or "snowbank" myxomycete species, which produce fruit bodies at the edge of melting snowbanks in spring. Relationship between the occurrence of fructifications and myxamoebae remain unknown. In this study we used modern molecular techniques, by direct DNA amplification from soil extracts (NGS metabarcoding) to compare the distribution of soil-inhabiting myxamoebae found in 2016 with fructifications from the same sites collected over the course of four years (2013, 2015-17) along an elevational transect in the northern German Alps. A coherent community composition between fructification and soil myxamoebae, though with species-specific differences in relative abundance, was revealed. Although patterns varied among species, myxamoebae were found at both low and high elevations, whereas fruit bodies were mainly found at higher elevations, likely explained by the presence of a stable and long-lasting snow cover. In addition, a year to year comparison of fructification records support the hypothesis that the abundance of fructifications strongly depends on the onset of snowfall in the previous autumn and the soil temperature regime throughout the winter.

Genetic barcoding of dark-spored myxomycetes (Amoebozoa)-Identification, evaluation and application of a sequence similarity threshold for species differentiation in NGS studies

Unicellular, eukaryotic organisms (protists) play a key role in soil food webs as major predators of microorganisms. However, due to the polyphyletic nature of protists, no single universal barcode can be established for this group, and the structure of many protistean communities remains unresolved. Plasmodial slime moulds (Myxogastria or Myxomycetes) stand out among protists by their formation of fruit bodies, which allow for a morphological species concept. By Sanger sequencing of a large collection of morphospecies, this study presents the largest database to date of dark-spored myxomycetes and evaluate a partial 18S SSU gene marker for species annotation. We identify and discuss the use of an intraspecific sequence similarity threshold of 99.1% for species differentiation (OTU picking) in environmental PCR studies (ePCR) and estimate a hidden diversity of putative species, exceeding those of described morphospecies by 99%. When applying the identified threshold to an ePCR data set (including sequences from both NGS and cloning), we find 64 OTUs of which 21.9% had a direct match (>99.1% similarity) to the database and the remaining had on average 90.2 ± 0.8% similarity to their best match, thus thought to represent undiscovered diversity of dark-spored myxomycetes.

Transcriptome reprogramming during developmental switching in Physarum polycephalum involves extensive remodeling of intracellular signaling networks

Activation of a phytochrome photoreceptor triggers a program of Physarum polycephalum plasmodial cell differentiation through which a mitotic multinucleate protoplasmic mass synchronously develops into haploid spores formed by meiosis and rearrangement of cellular components. We have performed a transcriptome-wide RNAseq study of cellular reprogramming and developmental switching. RNAseq analysis revealed extensive remodeling of intracellular signaling and regulation in switching the expression of sets of genes encoding transcription factors, kinases, phosphatases, signal transduction proteins, RNA-binding proteins, ubiquitin ligases, regulators of the mitotic and meiotic cell cycle etc. in conjunction with the regulation of genes encoding metabolic enzymes and cytoskeletal proteins. About 15% of the differentially expressed genes shared similarity with members of the evolutionary conserved set of core developmental genes of social amoebae. Differential expression of genes encoding regulators that act at the transcriptional, translational, and post-translational level indicates the establishment of a new state of cellular function and reveals evolutionary deeply conserved molecular changes involved in cellular reprogramming and differentiation in a prototypical eukaryote.

Speciation in progress? A phylogeographic study among populations of Hemitrichia serpula (Myxomycetes)

Myxomycetes (plasmodial slime molds, Amoebozoa) are often perceived as widely distributed, confounding to the “everything is everywhere” hypothesis. To test if gene flow within these spore-dispersed protists is restricted by geographical barriers, we chose the widespread but morphologically unmistakable species Hemitrichia serpula for a phylogeographic study. Partial sequences from nuclear ribosomal RNA genes (SSU) revealed 40 ribotypes among 135 specimens, belonging to three major clades. Each clade is dominated by specimens from a certain region and by one of two morphological varieties which can be differentiated by SEM micrographs. Partial sequences of the protein elongation factor 1 alpha (EF1A) showed each clade to possess a unique combination of SSU and EF1A genotypes. This pattern is best explained assuming the existence of several putative biospecies dominating in a particular geographical region. However, occasional mismatches between molecular data and morphological characters, but as well heterogeneous SSU and heterozygous EF1A sequences, point to ongoing speciation. Environmental niche models suggest that the putative biospecies are rather restricted by geographical barriers than by macroecological conditions. Like other protists, myxomycetes seem to follow the moderate endemicity hypothesis and are in active speciation, which is most likely shaped by limited gene flow and reproductive isolation.