{'fraction_non_alphanumeric': 0.03805376471284083, 'fraction_numerical': 0.018320813864079966, 'mean_word_length': 5.013025695741087, 'pattern_counts': {'":': 0, '<': 11, '<?xml version=': 0, '>': 10, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Phylogeny and character evolution in the Dacrymycetes, and systematics of Unilacrymaceae and Dacryonaemataceae fam. nov.

We present a multilocus phylogeny of the class Dacrymycetes, based on data from the 18S, ITS, 28S, RPB1, RPB2, TEF-1α, 12S, and ATP6 DNA regions, with c. 90 species including the types of most currently accepted genera. A variety of methodological approaches was used to infer phylogenetic relationships among the Dacrymycetes, from a supermatrix strategy using maximum likelihood and Bayesian inference on a concatenated dataset, to coalescence-based calculations, such as quartet-based summary methods of independent single-locus trees, and Bayesian integration of single-locus trees into a species tree under the multispecies coalescent. We evaluate for the first time the taxonomic usefulness of some cytological phenotypic characters, i.e., vacuolar contents (vacuolar bodies and lipid bodies), number of nuclei of recently discharged basidiospores, and pigments, with especial emphasis on carotenoids. These characters, along with several others traditionally used for the taxonomy of this group (basidium shape, presence and morphology of clamp connections, morphology of the terminal cells of cortical/marginal hyphae, presence and degree of ramification of the hyphidia), are mapped on the resulting phylogenies and their evolution through the class Dacrymycetes discussed. Our analyses reveal five lineages that putatively represent five different families, four of which are accepted and named. Three out of these four lineages correspond to previously circumscribed and published families (Cerinomycetaceae, Dacrymycetaceae, and Unilacrymaceae), and one is proposed as the new family Dacryonaemataceae. Provisionally, only a single order, Dacrymycetales, is accepted within the class. Furthermore, the systematics of the two smallest families, Dacryonaemataceae and Unilacrymaceae, are investigated to the species level, using coalescence-based species delimitation on multilocus DNA data, and a detailed morphological study including morphometric analyses of the basidiospores. Three species are accepted in Dacryonaema, the type, Da. rufum, the newly combined Da. macnabbii (basionym Dacrymyces macnabbii), and a new species named Da. macrosporum. Two species are accepted in Unilacryma, the new U. bispora, and the type, U. unispora, the latter treated in a broad sense pending improved sampling across the Holarctic.


INTRODUCTION
The Dacrymycetes (Basidiomycota) is one of the three widely accepted classes in the Agaricomycotina, sister to the Agarico mycetes. It is a monophyletic group comprised of widespread, saprotrophic fungi, mainly characterised by the strikingly long bisterigmate (forked) holobasidia (Hibbett et al. 2007), Unilacry ma unispora being the only known exception (Shirouzu et al. 2013). The class is further characterised at the ultrastructural level by the presence of continuous parenthosomes covering the dolipore-type septa (Oberwinkler 1994, Wells 1994). In addition, basidiomata are gelatinous to waxy or fleshy-cartilaginous, displaying a wide variation in size and shape, and very frequently yellowish to orange due to the pre sence of carotenoids (Gill & Steglich 1987), the main pigments of this group (Fig. 1).
The Dacrymycetes are typically brown-rotters (Worrall et al. 1997), although some have been reported to degrade considerable amounts of lignin (Seifert 1983). Gene contents from two Calocera genomes also suggests a brown-rotting lifestyle , and Dacrymycetes are the clade with the oldest stem age among Basidiomycota wood decomposers. As a result, their ancestors were likely among the first basidio-mycetous fungi able to decay wood (Oberwinkler 1994, Kirk et al. 2008, Floudas et al. 2012. Floudas et al. (2012), based on genome-scale data (26 DNA regions) in which Dacrymycetes was represented by a single species, estimated the stem age of the class to around 350 million years ago (Ma). A similar result, based on phylogenetic analysis of up to six DNA regions, was reached by Zhao et al. (2017), who reported a mean stem age of 358 Ma. Garnica et al. (2016) estimated the stem age to around 300 Ma using four DNA regions, but the Dacrymycetes clade was represented by only two species in the Dacrymy cetaceae. Recently, Varga et al. (2019) performed a variety of dating analyses on a phylogenomic dataset including three Dacrymycetaceae species and obtained mean stem ages for the Dacrymycetes around 278-504 Ma, although mostly in the 360 -385 Ma range. In view of this, the ancestors of this class and the Agaricomycetes seem to have split during the Upper Devonian or the Carboniferous.
The Dacrymycetes was traditionally comprised of a single order, Dacry mycetales, with two families, Dacrymycetaceae and Ceri nomycetaceae, but a second order called Unilacrymales was erected to accommodate the single known species with simple basidia, Unilacryma unispora (Shirouzu et al. 2013). The two orders were accepted in a recent study dealing with high-level taxonomy of the Basidiomycota (Zhao et al. 2017), although no species of Cerinomycetaceae were included. 2008). The shape of the fruitbodies, basidioma growth direction, presence or absence of a stipe, wall thickness of internal hyphae, and cell differentiation of cortical/marginal hyphae (partially related to the presence/absence of a sterile cortex) are among the most important taxonomic characters traditionally used to distinguish the genera ( Fig. 1-2). In a revision based mainly on these characters, Oberwinkler (2014) accepted the following 13 genera: Arrhyti dia, Calocera, Cerinomyces, Da crymyces, Dacryomitra, Dacryonaema, Dacryopinax, Dacryo scyphus, Ditiola, Femsjonia, Guepiniopsis, Heterotextus, and Unilacryma. However, some of the mentioned characters, e.g., the shape of the fruitbodies, were demonstrated to be homo-plastic by Shirouzu et al. (2013), and several common genera are polyphyletic (e.g., Calocera, Cerinomyces, Dacrymyces). In the same study, it was shown that other phenotypic traits, such as the presence or absence of clamp connections, not previously used to distinguish genera, are much better correlated with the phylogeny based on DNA data. For instance, species with and without clamp connections in Calocera, or in Cerinomyces, are not closely related.
Carotenoid pigments typical of the Dacrymycetes are synthetised by other fungal groups, but they are relatively unusual among non-photosynthetic eukaryotes (Echavarri-Erasun & Johnson 2002). They are prominent only in a few fungal taxa,  and as a result they have been used as taxonomic characters both in Ascomycota (e.g., Eckblad 1968, Schumacher 1990, Baral 2004 and Basidiomycota (e.g., Arpin & Fiasson 1971, Olariaga et al. 2015. Their acquisition and secondary loss in different phylogenetic lineages are well-known (Hansen et al. 2013, Olariaga et al. 2015, but remains to be investigated in the Dacrymycetes. Species delimitations are challenging in the Dacrymycetes. According to Hibbett et al. (2007), this is the smallest class of the subphylum Agaricomycotina, with 110-120 species known so far (Kirk et al. 2008, Shirouzu et al. 2009. However, judging from recent studies and our experience with fresh and herbarium samples, the real number of species is certainly larger, probably over 200. Basidiospore dimensions and septation, and the morphology of hyphidia are important characters for morphological species delimitation (McNabb 1964(McNabb , 1965a(McNabb -e, 1966(McNabb , 1973, but their variation may be considerable even within a single fruitbody. Shirouzu et al. (2017) provided the most complete phylogeny in terms of included vouchers with both DNA data and morphological data from fruitbodies. Some of the tentatively identified species turned out to be polyphyletic (e.g., Calocera cf. guepinioides, Cerinomyces crustulinus, Ce. pallidus, Dacrymyces microsporus, D. minor, D. novaezelan diae, D. stillatus) or appeared on exceptionally long branches within their clade, indicating there is considerable sequence variation (e.g., Calocera cornea, Ca. lutea, Dacrymyces cap itatus, Heterotextus miltinus). These groups, among others, may contain so far unrecognised species.
The class was subject to a worldwide monograph by Robert McNabb, who described new taxa from a variety of countries and ecozones in a series of studies (McNabb 1964(McNabb , 1965a(McNabb -e, 1966(McNabb , 1973 based on numerous other previous treatments, including those of, e.g., Brasfield (1938) and Kennedy (1959a, b) from the USA, or Kobayasi (1939a, b) from East Asia. During the last decade, significant contributions to the understanding of this group, including several studies based on morphological, molecular and ecological data, mainly on the Japanese and Australasian species, have been published by, e.g., Shirouzu et al. (2009Shirouzu et al. ( , 2013Shirouzu et al. ( , 2016Shirouzu et al. ( , 2017. In Europe, this group sparked the interest of Neuhoff (1936), and was later addressed in considerable detail by Donk (1966) and Raitviir (1967, including a substantial part of northern Asia). Reid (1974) presented a monograph focused on the British taxa, including 20 species, two of which were proposed as new. Since then, however, there have been no thorough revision of the European taxa with modern techniques. In the Nordic countries, Nannfeldt's (1947) monotypic genus Dacryonaema stands out as infrequently recorded, mainly from Northern Europe (Torkelsen 1997, http:// www.gbif.org/species/2512811), although it has also been reported from North America (Brough & Bandoni 1975, Redhead 1997. No DNA data were available for Dacryonaema at the beginning of this study, and the phylogenetic relationships of this morphologically deviant genus with the remaining species of Dacrymycetes are largely unknown (Oberwinkler 2014). On the other hand, the monotypic genus Unilacryma, described from North America and recorded from Japan, has recently been found also in Europe (Dam & Dam 2014). Shirouzu et al. (2013) commented on the morphological characteristics of the only known species, U. unispora, and compared it with another morphologically close species, Dacrymyces ovisporus, based on basidia and basidiospore characters. Dacrymyces ovisporus was described from Germany and has a wide distribution, including, e.g., northern Europe (Tor kelsen 1997) as well as North America (Bandoni 1963). No DNA sequence data of this species are available yet, and its taxonomic position remains unclear.
The objectives of this paper are: (i) to establish phylo genetic relationships among genera and species in the Dacrymycetes (including the genus Dacry onaema) based on multilocus DNA data; (ii) to identify phenotypic characters that allow the recognition of the phylogenetic lineages at family level, paying particular attention to those previously overlooked; and (iii) to revise the systematics and nomenclature of Dacryo naemataceae and Unilacrymaceae based on an integrative taxonomy approach, with special emphasis on the European taxa.
The third aim was taxonomically and geographically restricted because a wider scope would have been hampered by the extensive diversity and the large number of published names. The systematics of the larger Cerinomycetaceae and Dacry mycetaceae will be investigated in collaboration with other researchers in forthcoming studies.

Taxon and marker sampling
We sampled representatives for an overall phylogeny based on their morphology and phylogenetic position in previous studies. Unilacryma unispora s.lat. in the Unilacrymaceae, and Da cryonaema rufum with an unknown placement. Additionally, we included Dacrymyces macnabbii and D. ovisporus, which have not been included in any phylogeny yet, and two unidentified species that, according to morphological and preliminary molecular data, may belong in Dacryonaema or Unilacryma. We paid particular attention to selecting only visibly clean and non-parasitised fruitbodies. In addition, we added all available Dacrymycetes samples in GenBank being represented by at least three of the eight DNA regions selected (see below), for this purpose considering the ITS (ITS1 + 5.8S + ITS2) as a single region (Table 1). This includes mostly samples from the AFTOL project (Hibbett et al. 2007), the studies on the Dacrymycetes by Shirouzu et al. (2007Shirouzu et al. ( , 2009Shirouzu et al. ( , 2013Shirouzu et al. ( , 2016, and four Dacrymycetes species in MycoCosm (Grigoriev et al. 2014). In order to better assess the consistency of some phenotypic characters at the family level, we expanded the taxonomic sampling to include most available species with both a reasonable amount of molecular and morphological data. For samples with at least two DNA regions available in GenBank, we included data of up to two individuals of each taxon identified to species level when it was either not previously included or represented by a single individuum. When selecting among multiple available sequences, preference was given first to type specimens and second to the longest sequences; most of these additional sequences were generated by Shirouzu et al. (2017). Finally, and as an exception, a partial 28S sequence of Dacryoscyphus chrysochilus (AY604567, from the holotype specimen) was included because this was the only available DNA sequence of the type of Dacryoscyphus.  Table 1 GenBank accession numbers for Dacrymycetes sequences used in the phylogenetic analyses (except those from the four genomes indicated in Table 2). Taxon names refer to our post-hoc knowledge and are frequently provisional pending specific studies on certain groups. Environmental nucleic acid sequences without vouchers are indicated as 'ENAS', and are tentatively assigned to a species name based on sequence similarity only. New sequences generated for this study are marked in bold. *Holotype specimen. Unpub. = Unpublished DNA sequence from the AFTOL database.   This sampling comprised c. 90 ingroup species (not all of them identified) and included, at least in a broad sense, the types of all genera accepted by McNabb (1964McNabb ( , 1965aMcNabb ( -e, 1966McNabb ( , 1973, and all genera accepted by Oberwinkler (2014) except Arrhytidia. The type of Arrhytidia (A. flava) was considered a species of Dacrymyces (D. marginatus), placed close to the D. chrysocomus group by McNabb (1973). No recent collections identified as such are known.
We selected eight DNA regions for phylogenetic analy ses: three linked, multicopy regions in the nuclear ribosomal DNA, viz., 1) the nuclear ribosomal small subunit (18S); 2) the internal transcribed spacer (ITS) region (including the ITS1, 5.8S, and ITS2); and 3) the D1-D3 domains of the nuclear ribosomal large subunit (28S). We also selected three unlinked, low-copy, protein-coding, nuclear regions, viz., 4) a c. 900 -1 100 bp long fragment from the largest subunit of the DNA-directed RNA polymerase II gene between domains A and C (RPB1); 5) a c. 1 500 -1 700 bp long fragment from the second largest subunit of the DNA-directed RNA polymerase II gene between domains 6 and 11 (RPB2); and 6) a c. 900-1 000 bp long frag ment at the end of the translation elongation factor 1 alpha (TEF-1α) gene. Finally, we included two mitochondrial regions, viz., 7) the mitochondrial ribosomal small subunit (12S); and 8) the ATP synthase subunit 6 (ATP6).
An outgroup was selected by retrieving from MycoCosm (Grigoriev et al. 2014) representatives from 20 of the c. 26 orders currently accepted of Tremellomycetes (three species) and Agari comycetes (17 species). Samples were selected (Table 2) to minimise the amount of missing data in the outgroup. For each sample, nrDNA and mitochondrial regions were retrieved through BLAST searches using already generated sequences (e.g., from GenBank) of the same species or close relatives, keeping the best scoring hit. RPB1, RPB2, and TEF-1α sequences were retrieved directly from the annotated contigs through the MycoCosm search option.
For morphological and cytological studies, the following fresh samples from Sweden of species also molecularly analysed were used: Dacrymyces chrysocomus s.lat. (Uppland,Uppsala,Lunsen,N59 We also performed an independent and denser sampling to investigate species delimitations in Dacryonaema and Unila cryma. In total, we included DNA data of eight individuals of Da. rufum from Canada, Finland, Norway, and Sweden, 12 of Da. macnabbii from the UK (holotype), Finland, and Sweden, six of an unidentified and potentially new Dacryonaema spe-  cies from Finland and Norway, seven of Unilacryma unispora s.lat. from Finland, Japan, and Sweden, and six of a bisporic Unilacryma from Finland and Sweden. In addition, all available GenBank sequences that had a high similarity in BLAST searches (≥ 99 %) to any of the studied DNA regions were added (Table 1). Several additional specimens of Da. macnab bii, Da. rufum, U. unispora s.lat., and the bisporic Unilacryma were collected and morphologically studied to gather a better knowledge of their intraspecific variation (see taxonomic part).

Morphological study
Whenever possible, samples were studied morphologically and photographed when fresh, either immediately after being collected or otherwise kept in a refrigerated (4 °C) humid chamber for a maximum of 48 h. Alternatively, when samples were found in a dry state in the field, they were either studied directly, if the structures were well-developed, or moistened with tap water by submerging the entire substrate a few minutes, and then placed in a refrigerated humid chamber for one to seven days before study. General terminology largely follows Oberwinkler (1994), and for clamp connections we follow Stalpers (1978). In particular, 'hyphidia' is preferred over 'dikaryophyses', as the last term is not consistently applied for dikaryotic hyphidia in other groups of the Agaricomycotina, particularly in aphyllophoroid fungi. Also, we prefer the term 'young basidium' instead of 'probasidium' to avoid confusion with the thick-walled structures in Pucciniomycotina and Ustilaginomycotina, from which the basidium forms (Clemençon 2004).
Macromorphology was studied under a Leica MZ 75 dissecting microscope, and macroscopic photographs were taken with a coupled Infinity 1 macro camera. Micromorphological study was performed by mounting thin hand-cut sections of the fruitbodies in water or a 5 -10 % w/v aqueous solution of KOH, and gently squeezing them under the cover slip. All measurements were taken at 630× magnification in KOH solution. The hilar appendix was not included in basidiospore measurements. Microscopic photographs were taken in a Zeiss Axio Imager A1 microscope, fitted with differential interference contrast (DIC) and equipped with an AxioCam ICc3 digital camera. For some of the most voluminous microscopic structures, a series of photographs were stacked using Combine ZP.
An aqueous solution of 1 % w/v potassium triiodide (KI 3 , Lugol's solution) was used to check the reaction of the carotenoid contents in some species. Nuclear staining was done by using toluidine blue. The best results were achieved when fresh samples were pre-treated with 5 -10 % KOH, quickly rinsed with tap water, and stained with a recently prepared alkaline solution (pH ≈ 10.5 -11) of 0.1 % toluidine blue.
We performed analyses of variance (ANOVA) on the basidiospore length, basidiospore width, and the length/width ratio (Q) of Dacryonaema spp. and Unilacryma spp. samples, in order to detect significant differences (p < 0.001) among the species in each genus. Diagnostic plots for residuals (residuals vs fitted values, normal Q-Q, scale-location, and residuals vs leverage) were checked to ensure that no model assumptions were severely violated. Cohen's f was used as an estimator of effect size. Since Dacryonaema included three species, when any significant difference was found in the ANOVA, the Tukey's honest significant difference (Tukey's HSD) post hoc test was applied to identify pairs of means different at p < 0.001. All analyses were done using the basic package of R (R Core Team 2018) and sjstats v. 0.17.4 (Lüdecke 2019  After study, specimens were dried at room temperature and deposited in UPS. Additional herbarium specimens from F, K, NCU, O, S, TENN, UBC, and UPS were studied (herbarium acronyms following Thiers 2019). Geographical coordinates are provided in the World Geodetic System 1984 (WGS84).

Lab work and primer design
Genomic DNA was extracted from a single fruitbody or part of a fruitbody (c. 1 mm 3 ), rarely (species with very small fruitbodies like Dacryonaema rufum) from 2 -5 adjacent fruitbodies. The samples were disrupted in 1.5 ml microcentrifuge tubes with a plastic micropestle either after being frozen in dried state with liquid nitrogen, or directly in the lysis buffer. We used the DNeasy Plant Mini Kit (Qiagen, Germany), following the manufacturers' instructions, but increased the lysis buffer incubation time to up to two hours. For some samples, including the type specimens of Dacrymyces macnabbii (K(M) 81678) and Unilacryma unispora (NCU F-0026842), we followed the Chelex extraction protocol of Ferencova et al. (2017), in order to minimise the amount of material used.
Forty-six different primers were used for amplifying the DNA regions mentioned above (Table 3), comprising both general primers used in other studies as well as newly designed primers for 12S, RPB1, RPB2, and fragments of ITS1 and ITS2 (specific primers for these last two loci designed only for Dacryonaema and Unilacryma). ITS4B (Gardes & Bruns 1993) could not be used because the primer site has four mismatches in the genus Unilacryma, one of them affecting the last base at the 3'-end, and three mismatches in the genus Dacryonaema, in both cases resulting in unsuccessful PCR amplifications. Primer design was carried out by identifying conserved regions in Dacrymycetes alignments containing already generated sequences of those loci from MycoCosm, GenBank and own sequences obtained using general primers. We tried to design the primers so that their CG content was as close as possible to 50 %, with a length in the range of 17-25 bp, melting temperatures differing less than 5 °C within a pair, and avoiding hairpin and duplex formation. Calculations were done using OligoAnalyzer v. 3.1 (Integrated DNA Technologies®).
PCR reactions were performed in a final volume of 20 -25 μl using MasterMix (Qiagen, Germany). Numerous primer combinations and PCR conditions were tried for each DNA region, so we only report the most successful PCR conditions used for amplifying the majority of samples. The PCR programmes always consisted on an initial denaturation of 2 min at 94 °C and a final extension step of 7 min at 72 °C. The central cycles and primer pairs for the different loci greatly varied and are summarised in  PCR products were purified using ExoSAP-IT TM (Applied Biosystems) or Illustra TM ExoProStar TM (GE Healthcare), following the manufacturer's instructions but adding 1/4 of the recommended volume of enzyme mix and extending the reaction time to 30 min. Purified PCR products were sequenced by Macrogen with the same primers used for amplification.

Sequence edition and alignment
Raw chromatograms were assembled and edited using Sequencher v. 4.1.4 (Gene Codes, USA). Even if most DNA extractions were generated from a single fruitbody, chromatograms with scattered double peaks were found in several samples. Positions with a double signal were often similar or identical among different samples (thus, from independent extractions and PCRs) of a single species, more abundant in introns and, when affecting the exons of protein coding genes, they were almost always restricted to the third codon position, resulting in alternative codons for the same amino acid (synonymous mutations). This observation suggests that more than one allele was present for a particular DNA region in some samples. For the purpose of phylogenetic reconstruction, these positions were coded using ambiguity codes. In case there was a polymorphism affecting the length of the alleles (e.g., in ITS, and often only a single base indel), the shortest copy was retained.
Two datasets were assembled and analysed independently, one for higher level taxa (overall phylogeny) and one for the species delimitation analyses of the genera Dacryonaema and Unilacryma. The ITS1 region was considered to start immediately after the highly conserved motif 'AGGACATTA' at the end of the 18S. The 5' limit of the 5.8S region was defined by comparing the aligned sequences to annotated sequences from ITSoneBD (Santamaria et al. 2017). The 3'-end of the 5.8S and the 5' start of the 28S were defined using the annotate option of the ITS2 database (Keller et al. 2009). The ITS1 and ITS2 outgroup sequences were excluded for being largely unalignable among the distantly related taxa included in the dataset for the overall phylogeny. Likewise, repetitions at the end of the 28S (in Cerinomyces cf. crustulinus, C. lagerheimii, Dacrymyces flabelliformis, D. subantarcticensis, and Dacryopi nax indacocheae) were removed, as well as two intron regions above 600 bp each in the 12S (present in Vuilleminia comedens and Cerinomyces ceraceus genomic data, and in the newly generated sequences of Calocera cornea, Cerinomyces tor tus s.lat. 1, and Dacrymyces cf. minutus s.lat. 2). Additionally, we excluded a highly variable and gappy region in the 12S appearing between two easily alignable regions (the one in the 5' ending with the motif 'AGGGTTYGTAG', and the one in 3' starting with the motif 'TAGAATCW'). For the protein-coding DNA regions, intron positions were identified using open reading frames from the JGI and annotated sequences from GenBank as references (with few exceptions), and also excluded from the matrices. All sequences were trimmed to fit to the longest sequences generated by us and to remove poor read ends.
For the overall phylogeny, all alignments were estimated using PASTA v. 1.7.8 (Mirarab et al. 2015). We performed three independent runs with 10 iterations, maximum subset size being half of the sequences included in the alignment, the mask option deactivated, and keeping the best likelihood score alignment. MAFFT (L-INS-i option) was used as the aligner, OPAL as the merger, and RAxML as the tree estimator. For the ribosomal DNA regions we used the GTR + Γ model, and for the exon regions of the protein coding genes, the WAG + I + Γ model for the datasets temporarily translated into amino acids. The alignment with the best score after the three runs was manually ATP6-3/ATP6-2 or ATP6-3/ATP6-4 Programme according to Kretzer & Bruns (1999)  corrected (minor modifications) and submitted to a last run in PASTA, using the parameters indicated above. Since PASTA already optimises the likelihood under a phylogenetic framework, and the resulting alignments seem to be phylogenetically accurate (Mirarab et al. 2015), neither further manual refinement was done nor parts of the alignments were excluded. All the matrices independently aligned in PASTA contained less than 15 % of missing data at the ends of the sequences, except for the RPB2 dataset where c. 1/3 of the sequences (all from GenBank) were only represented by a fragment of 843 bp covering an easily alignable exon region (a single indel due to three additional codons in the Tremellomycetes sequences).
For the Dacryonaema + Unilacryma dataset, the 5.8S, TEF-1α, and ATP6 alignments contained no indels. The RPB2 alignment contained a single, unambiguous indel of three bases in length, corresponding to a codon position. The 18S and 28S alignments contained only three and two unambiguous single-base indels, respectively, and were manually aligned. The Dacryonaema and Unilacryma RPB1 datasets contained a single, unambiguous codon indel and were aligned and analysed separately because the incongruence detected in the overall phylogeny (see 'Congruence among datasets' below) suggests that they may be non-orthologous. The ITS1 and ITS2 regions were independently aligned using BAli-Phy v. 3.3 (Suchard & Redelings 2006). Since BAli-Phy is computationally intensive, we reduced the datasets by selecting two or three sequences within the same species, keeping the maximum variation but avoiding those sequences with more missing data or ambiguous base calls. A GTR + I + Γ model was used for the substitutions, whereas an rs07 model was used for inserations and deletions. For each dataset, we executed 6 simultaneous runs with 5 × 10 4 iterations per run, checked convergence, and computed the alignments discarding the first 25 % of the runs as burn-in. The average standard deviation of splits in the posterior sample turned out to be smaller than 0.005, and all parameters analysed in Tracer v. 1.7 (Rambaut et al. 2018) had effective sample sizes exceeding 1 000 for individual runs (≥ 7 000 when combined). The resulting alignments were used as backbones for manually re-introducing the remaining sequences. Indels were coded as a separate binary dataset (gap present = 1, gap absent = 0) in FastGap v. 1.2 (Borchsenius 2009), using the simple indel coding method (Simmons & Ochoterena 2000). All final alignments were deposited in TreeBASE (http://purl.org/phylo/treebase/ phylows/study/TB2:S24147).

Congruence among datasets
Dataset congruence was assessed for both the overall phylogeny and the Dacryonaema + Unilacryma datasets through maximum likelihood (ML) analyses using IQ-TREE v. 1.6.8 . We performed standard, non-parametric bootstrap analyses with 200 replicates for the nrDNA dataset, each of the three nuclear protein-coding genes, and the mitochondrial dataset, assuming no recombination within each dataset. Each of the protein coding genes was partitioned by codon position, the nrDNA dataset was originally divided into five partitions (18S, ITS1, 5.8S, ITS2, and 28S), and the 12S and ATP6 genes were considered as potentially different partitions. Additionally, we also tested for incongruence among partitions of the nrDNA dataset. Incongruence was considered significant when a strongly supported clade from one phylogeny was contradicted by another strongly supported clade in another phylogeny (Mason-Gamer & Kellogg 1996). A bootstrap support (BS) threshold of 75 % was used, at or above which a clade was considered strongly supported. A significant case of incongruence (Dacryonaema in the RPB1 dataset) in the overall phylogeny analyses was inspected for alternative explanations.
Possible deviations from nucleotide composition heterogeneity among taxa were checked with the IQ-TREE composition test.
The amount of recombination, if any, was tested with GARD (Kosakovsky Pond et al. 2006) and SlidingBayes v. 0.94 (Paraskevis et al. 2005). Finally, we checked for heterotachy by performing analyses under the edge-unlinked partition model (Lopez et al. 2002) and the GHOST model (Crotty et al. 2017) in IQ-TREE. In the end, these processes as well as incomplete lineage sorting turned out to be unlikely candidates for the observed topological conflict, and consequently the Dacryonaema RPB1 sequences were excluded from downstream analyses.

Phylogenetic analyses for the overall phylogeny
Two approaches were used for calculating phylogenies: (i) tree estimation based on the concatenated dataset of all loci; and (ii) coalescent-based species tree estimation joining independent single-locus tree phylogenies. In this set of analyses, gaps were treated as missing data. For the concatenated analyses, two strategies were used for phylogenetic reconstructions: maximum likelihood (ML) and Bayesian Inference (BI). ML analyses were performed using IQ-TREE v. 1.6.8 , and the model selection and partitioning scheme for the mitochondrial and nuclear regions were calculated by using the Bayesian Information Criterion in the integrated version of ModelFinder (Kalyaanamoorthy et al. 2017). Fourteen potential partitions were defined for the nuclear regions: 18S, ITS1, 5.8S, ITS2, 28S, and one for each codon position of RBP1, RPB2, and TEF-1α, and four potential partitions were defined for the mitochondrial regions: 12S and each codon position of ATP6. Branch lengths were linked and proportional across partitions. The analysis was repeated twice, starting from random trees. Branch support was assessed by standard non-parametric bootstrap, performing 500 replicates, and through the Shimodaira-Hasegawa-like approximate likelihood ratio test (SH-aLRT, Guindon et al. 2010), performing 1 000 replicates. The purpose of using the SH-aLRT was to identify very short branches (with no signal from the data) that may have received spuriously high BS support (Guindon et al. 2010). Significant support is considered when BS ≥ 75 % and SH-aLRT ≥ 85 %.
Bayesian inference on the concatenated dataset was carried out using MrBayes v. 3.2.6 (Ronquist et al. 2012). We applied the same partitioning scheme obtained in the previous analysis with IQ-TREE, unlinking model parameters across partitions, and using model jumping to integrate over models and their parameter values in each partition (Huelsenbeck et al. 2004), allowing a gamma distributed rate heterogeneity across sites (approximated by four categories) and a proportion of invariant sites. We used a (1, 1, 1, 1, 1, 1) Dirichlet prior for the substitution rates, a (1, 1, 1, 1) Dirichlet prior on the state frequencies, and a uniform (0, 1) prior for the proportion of invariable sites. Branch lengths were linked and proportional across partitions. We used a compound Dirichlet prior on branch lengths in the overall phylogeny (1, 0.0588, 1, 1), based on tree length estimates from the best tree of the ML analysis. We set the temperature increment parameter to 0.1 after having observed swap rates between the cold and hot chains in preliminary runs. Four parallel runs starting with random trees were executed, each with four chains, for a maximum of 1 × 10 8 generations.
The analysis was automatically stopped before reaching the maximum number of generations if the average standard deviation of split frequencies dropped below 0.01. Sampling was done every 1 000th tree. The first 50 % of the analysis was eliminated as burn-in, and Tracer v. 1.7 (Rambaut et al. 2018) was used to ensure that effective sample size (ESS) for each parameter exceeded 200. RWTY (Warren et al. 2017) was used to check the convergence of tree topologies. The 50 % majorityrule tree with posterior probabilities (PP MrBayes , consi dered signi-ficant when ≥ 0.95) and average branch lengths was calculated from the post-burn-in trees, and visualised in FigTree v. 1.4 (Rambaut 2016).
ML and Bayesian approaches were also used for inferring a coalescent-based species tree from the overall phylogeny data. An ML species tree was estimated with ASTRAL III (Zhang et al. 2018, Rabiee et al. 2019) v. 5.6.3, using the independent ML phylogenetic trees obtained for each locus in IQ-TREE.
Since ASTRAL uses a summary method for computing a species tree based on single-locus trees, limited accuracy is expected when the number of loci is low. Considering this limitation, we conducted two analyses, first by considering the effectively linked parts of the nrDNA as a single locus (analysis assuming five loci), and secondly by splitting the nrDNA into its units: 18S, ITS1, 5.8S, ITS2, and 28S (analysis of nine loci).
The second analysis has the advantage of accounting for any unsupported discordances among those nrDNA regions, but has the shortcoming of giving too much weight to the nrDNA compared to other loci if recombination is absent among the mentioned nrDNA regions. Branches receiving ≤ 10 % BS support in each single-locus tree were collapsed using Newick Utilities v. 1.6 (Junier & Zdobnov 2010), as recommended by Zhang et al. (2018). Branch support was calculated as local posterior probabilities based on quartet support (QS, Sayyari & Mirarab 2016), and was considered significant when QS ≥ 0.95. Species assignments of individuals mostly followed our identifications and those from the literature (e.g., Shirouzu et al. 2009Shirouzu et al. , 2017, and tended to be conservative, i.e., in some cases very closely related and potentially cryptic species might have been merged: Calocera viscosa s.lat., Dacrymyces chrysospermus s.lat., and Unilacryma unispora s.lat. In these three cases, all samples in each group formed a fully supported clade in the previous concatenated ML and Bayesian analyses. Exceptions were made when we detected either non-monophyly or a striking divergence among individuals in the ML analysis, e.g., the three clades of Calocera cf. guepinioides were considered as three independent species, and Dacrymyces cf. minutus s.lat. 1 and D. aff. stillatus, from Shirouzu et al. (2009), were treated as different species to our D. cf. minutus s.lat. 2 and D. stillatus, respectively.
StarBEAST 2, implemented in BEAST 2.5.1 (Bouckaert et al. 2014, Ogilvie et al. 2017, was used to infer a dated species chronogram for the overall phylogeny, based on the multi species coalescent. In the absence of fossils for the included taxa, we used a secondary calibration for the Dacrymycetes, setting the stem age of the group by using a normal prior distribution with values for µ = 350.0 Ma and σ = 50.0, according to the results from Floudas et al. (2012) (also largely agreeing with Zhao et al. 2017), under an uncorrelated lognormal relaxed clock model (Drummond et al. 2006). We used the uninformative 1/X prior on clock rates to allow the calibration to guide those clock rates, and a calibrated Yule tree model. Clock and tree models were unlinked among the following loci: nrDNA, RPB1, RPB2, TEF 1α, and the mitochondrial regions. The original partitions (18S, ITS1, 5.8S, ITS2, 28S, 12S, and codon positions of each protein coding gene) were used since it is biologically implausible to combine partitions from different loci under the coalescent (e.g., codon positions of unlinked protein coding genes). We used bModelTest  for coestimating site models during the runs, averaging among all extended named models. The result of model averaging was summarised using BModelAnalyser (part of the BEAST 2 package, Bouckaert et al. 2014). Given the general life cycle in Agaricomycotina, we considered a biparental inheritance for the mitochondrion, and we used the analytical population size integration as the population model, since the number of individuals per species was low, and we were not specifically interested in per-species population sizes. Individuals were assigned to species as indicated above for the ML analysis in ASTRAL III. Initially, we executed six parallel runs starting from random trees for 5 × 10 8 generations, sampling every 10 4th tree, then checked the traces and ESSs in Tracer v. 1.7, and discarded three of them due to failure to reach convergence for some parameters (different parameters in each run and none affecting the statistics of the species tree, see Results and Discussion). The other three runs were resumed until completing 1 × 10 9 generations. The first 50 % of these three runs was eliminated as burn-in.
Tracer v. 1.7 (Rambaut et al. 2018) was used to check effective sample size (ESS) of model parameters (ideally > 200). RWTY (Warren et al. 2017) was used to check the convergence of the species tree topologies. The maximum clade credibility species tree with posterior probabilities (considered significant when ≥ 0.95) and node ages was calculated using TreeAnnotator v. 2.5.1 (part of the BEAST 2 package, Bouckaert et al. 2014) and visualised in FigTree v. 1.4.

Character evolution analyses
We inferred the history of selected morphological characters, including transformation counts and ancestral states at the nodes corresponding to the Dacrymycetes and the four major families Cerinomycetaceae, Dacrymycetaceae, Dacryonaema taceae, and Unilacrymaceae. We avoided further analyses of characters for which convincing inferences were already provided by Shirouzu et al. (2013), and not expected to be altered by our sampling, e.g., basidiospore septation pattern, fruitbody morphology, and presence or absence of a stalk. The following traits were analysed: (i) basidium morphology (0 = 2-spored, 1 = 1-spored); (ii) clamp connections presence and morphology (0 = clamps absent, 1 = short and rounded clamps present at most septa, 2 = long and angulose clamps mainly restricted to the bases of basidia and hyphidia, and subhymenial hyphae); (iii) hyphidia (0 = conspicuous and branched, 1 = unbranched or indistinct); (iv) terminal cells of cortical/marginal hyphae (0 = filamentose and simple, 1 = clearly dilated and simple, 2 = filamentose and branched); (v) carotenoids (0 = absent or indistinct, 1 = conspicuous at least in some cells); (vi) cytoplasmic lipid content (0 = indistinct and with no visible lipid drops in the cytoplasm, 1 = lipid drops visible when the lipid bodies coalesce); (vii) number of nuclei in recently discharged basidiospores before the formation of the first septum (0 = uninucleate, 1 = binucleate). Taxa with unknown states were coded as having prior probabilities equally divided among all known states. For simplicity, the character state matrix is shown next to the species tree obtained in the StarBEAST 2 analysis (Fig.  9). Our hypotheses, based on prior knowledge of the group, were as follows: (i) the two-sterigmate basidia, presence of clamps, branched hyphidia, filamentose and simple terminal cells of cortical/marginal hyphae, absence of conspicuous carotenoids, high cytoplasmic lipid content, and uninucleate young basidiospores are plesiomorphic in the Dacrymycetes; (ii) the long, angulose clamp connections at subhymenial hyphae and the bases of basidia and hyphidia, and the filamentose and branched terminal cells of cortical/marginal hyphae are plesiomorphic in the Dacryonaemataceae; (iii) binucleate young basidiospores and low cytoplasmic lipid content are acquired states plesiomorphic in the Cerinomycetaceae; (iv) the presence of carotenoids and unbranched hyphidia are acquired states plesiomorphic in the Dacrymycetaceae; (v) the 1-spored basi dium, the dilated terminal cells of cortical/marginal hyphae, and the clampless septa are apomorphic conditions within different families in the Dacrymycetes.
Character transformation counts and ancestral state reconstructions were performed by stochastic mapping with SIM-MAP (Bollback 2006), using the makeSimmap() function of the R package phytools (Revell 2012). We performed the analysis on a subsample of 1 000 trees from both the posterior distribution of the MrBayes phylograms and the StarBEAST 2 chronograms to explore whether there is any difference between assuming that the phenotypic change is proportional to genetic change or to time, following Kistenich et al. (2018). The MrBayes phylograms were pruned to retain only one sample per species (the one on the shortest branch) and to exclude the outgroup. Outgroup taxa were also excluded from the StarBEAST 2 chronograms. Since the analytical approaches used for obtaining the phylogenetic trees were also different (supermatrix vs supertrees), the resulting topologies may differ as well. Under these circumstances, it may be difficult to know whether any difference in results are caused by the topology or by the scaling of branch lengths. Therefore, in case a contrasting result was found for a given character, a third group of analyses was performed using chronograms generated from the MrBayes phylograms. These chronograms were calculated using penalised likelihoods under a correlated model (Paradis 2013) with the chronos() function of the R package ape v. 5.0 (Paradis & Schliep 2019). In all cases, ten simulations were performed for each tree. Otherwise, Kistenich et al. (2018) was followed, including the use of an instantaneous rate matrix, Q, estimated from the empirical data, and a prior distribution of states on the root node estimated from the stationary distribution of Q. However, we summarised the inferred ancestral states at the nodes of interest from all mapped trees, instead of later downsampling them, since the number of initial simulations was smaller. Reconstructions of ancestral states were focused on five nodes in the backbone of the Dacrymycetes, i.e., the crown node for the class Dacrymycetes and the four crown nodes for each named family. We found a quite obvious artefact in the reconstruction of the number of basidiospore nuclei at the root node (Dacrymycetes), which was inferred with 100 % probability as binucleate with the settings indicated above. This problem was caused by an unusual calculation of the Q matrix and was solved by using a symmetric model for that particular character. The problem can also be fixed by sampling Q from the posterior probability distribution of Q using MCMC, but this is much more computationally costly and impractical for large datasets.
Following Kistenich et al. (2018), we also checked whether ancestral state reconstructions at or near the root node were influenced by the assumptions made about the distribution of states at the root. The rayDISC() function of the corHMM v. 1.22 R package (Beaulieu et al. 2013) was used for this purpose.
Reconstructions were made on the pruned majority-rule consensus tree from the MrBayes analyses and on the maximum clade credibility tree from the StarBEAST 2 analysis. The same coding and state characters used for the stochastic mapping were applied. The three different assumptions about the weight given to the states at the root node allowed in rayDISC were used for calculating marginal reconstructions: (i) equal weighting among all states; (ii) based on the estimated transition rates (Yang 2006); and (iii) estimated from the data and the model (Maddison et al. 2007, FitzJohn et al. 2009).

Phylogenetic and species delimitation analyses in Dacryonaema and Unilacryma
The phylogeny of Dacryonaema and Unilacryma species based on the concatenated dataset was inferred following the same methods indicated for the overall phylogeny, i.e., ML analyses in IQ-TREE v. 1.6.7 and BI in MrBayes v. 3.2.6. The only differences were that we used the compound Dirichlet prior on branch lengths (1, 1.3333, 1, 1) for the Bayesian analysis, based on the different tree length estimates from the best tree of the ML analysis, and the temperature increment parameter was set to 0.2, resulting in good mixing.
For the purpose of delimiting species in Dacryonaema and Uni lacryma, we used the multispecies coalescent approach implemented in STACEY v. 1.2.4 (Jones 2017). Clock and tree model parameters were estimated independently for the nrDNA and each protein coding locus. We used an uncorrelated lognormal relaxed clock model (Drummond et al. 2006) and a Yule tree model. The nrDNA was divided into three partitions, (i) 18S + 5.8S + 28S, (ii) ITS1 + ITS2 and (iii) coded indels, and each protein coding locus was divided in (i) 1 st + 2 nd codon positions, (ii) the 3 rd codon position, and (iii) coded indels for each locus, if present. With this simplified scheme we tried to avoid overpartitioning, since the amount of variation in the Dacryonaema + Unilacryma dataset was relatively low, particularly the first and second codon positions of the protein coding genes, and parameter estimates when treating those partitions independently were sometimes poor in preliminary analyses. Parameters of the site models were independently estimated for each DNA partition using bModelTest , taking into consideration all transition/transversion split models, while the Lewis MK model with two states was used for the binary partitions. The most visited models were summarised using BModelAnalyser (part of the BEAST 2 package, Bouckaert et al. 2014). We used the CIPRES Science Gateway (Miller et al. 2010) to run ten MCMC parallel analyses for 1 × 10 9 generations, sampling every 5 × 10 4th tree, and with the collapse height parameter ε = 10 -4 . Considering the amount of missing data (most individuals from GenBank were represented by a single locus) and the taxonomic information from both morphology and the phylogenetic analyses on the concatenated dataset, we used an informative prior for the collapse weight parameter (ω). Since the most likely number of species according to the tree topology and information from morphology was around five or six, and impossible below two (that would imply considering all samples of Dacryonaema and Unilacryma not only as a single genus, but as a single species), we set ω = 0.908 (highest probability around κ = 5-6), with 0.5 (κ = 25) and 0.98 (κ = 2) as the lower and upper values of the interval, respectively. The first half of each run was discarded as burn-in, and the remaining posterior sample was used to calculate the most likely number of clusters (i.e., putative species), using SpeciesDelimitationAnalyzer (

Vacuolar cell contents and its variation with mounting media
The cytoplasm of living cells in all studied species of Dacrymy cetes was found to be filled with conspicuous, highly refractive vacuolar bodies (VBs) that are variable in size and shape within specimens but not between species (Fig. 3a, c, g, i, m, o, s, u). The VBs disappear when treated with alkaline lethal media, such as KOH, or when cells die by other causes, e.g., desiccation (Fig. 3b, d -f, h, j-l, n, p-r, t, v-x).
Lipid bodies (LBs) have been observed in nearly all studied samples but, unlike VBs, they show considerable variation. Conspicuous LBs are always absent in living cells, being reduced to hardly discernible vacuoles of granular appearance (≤ 0.2 μm diam, unaffected by KOH treatment) due to the presence of the much more conspicuous VBs. When cells die, the LBs become more visible and tend to gradually fuse, forming larger oil drops when the amount of LBs allows so. The difference in the amount of visible cytoplasmatic lipid content among groups is more pronounced in basidiospores, less so in basidia, and even less in somatic cells.
In Cerinomycetaceae, the cytoplasm of samples treated with KOH appears either as undifferentiated or, at most, as only slightly granular, with an inconspicuous lipid content (Fig. 3b, d -e). Big oil drops (> 2 μm diam) are rarely seen in dead cells of samples of this family (Fig. 3f), neither in rapidly desiccated samples nor in samples exposed to reagents that favour the coalescence of LBs, such as ethanol or HCl. However, exceptions occur, and D. tortus s.lat. 1 displayed a rather high lipid content, comparable to some species in the Dacrymycetaceae.
In Dacryonaemataceae, there is a large variation among the three species (Fig. 3h, j -l). Dacryonaema rufum has a comparatively low cytoplasmic lipid content, and LBs are not visible or hardly discernible in dead cells, even after forcing their coalescence (Fig. 3h, j). Dacryonaema macnabbii and especially Da. macrosporum may show a considerable amount of LBs, and some big lipid drops may form when the LBs coalesce ( Fig. 3k -l).
In all studied species of Dacrymycetaceae and Unilacryma ceae the cytoplasm of various cells (particularly in basidia and basidio spores) appears clearly granular when samples are mounted in KOH (Fig. 3n, p, t, v). In addition, a few big oil drops (> 2 μm diam) are commonly observed in at least some of the dead cells ( Fig. 3q -r, w-x).
Coalesced LBs are often visible in well-preserved herbarium samples, even decades after being collected. However, the cytoplasm tends to become rather uniformly refractive with age in some old herbarium specimens, in which case internal details are lost.

Other variations with mounting media
Cell size and shape displayed little variation among living or dead cells mounted in water and dead cells mounted in KOH (Fig. 3). On the other hand, due to the commonly large size of many structures and the elasticity of the cell walls, dimensions are easily altered by simply applying pressure on the cover slip, or by premature mounting of insufficiently rehydrated specimens. The terminal cells of cortical/marginal hyphae in Dacryonaema and a few other genera (e.g., Heterotextus) are an exception, since walls become distinctly gelatinised and expanded in KOH.
Hyphae with somewhat roughened walls (e.g., some internal hyphae in various species, Fig. 2d) are often more distinctly roughened when the sections are mounted in KOH, perhaps due to the precipitation of some substances of the intercellular gelatinous matrix.

Nuclei
Nuclei are masked by the VBs in the living state, and are best studied in recently dead cells. In species with highly pigmented cytoplasms, nuclei appear as 'empty' areas due to the absence of pigments, and are relatively easy to observe in fresh samples mounted in KOH (Fig. 4c, e-g, k-m). They are also observable under the same conditions in species with a very low amount of lipid contents (Fig. 4j, n, q). Nuclei can sometimes be seen in slides prepared with Lugol's solution as slightly paler areas of the cytoplasm when the cells are still alive or have recently died (Fig. 4t). Sometimes, they appear as slightly darker areas (due to staining) when the cells are dead. The best contrast is often achieved by using alkaline toluidine blue solutions, where the nuclei (or at least the nucleoli) stain deeply blue, while the cytoplasm stains less intensely and often more purplish, depending on the species and small variations in pH ( Fig. 4a - Vegetative cells in all fresh samples studied are dikaryotic, and mitospores formed directly from those hyphae, e.g., the arthrospores of Dacrymyces stillatus are dikaryotic as well . Basidia under formation are originally dikaryotic, but the two haploid nuclei quickly fuse into a single, large, diploid nucleus ( Fig. 4e -f, h). At a certain state of development, often when the sterigmata are already formed, meiosis occurs and sometimes the four postmeiotic nuclei are revealed (Fig. 4g, i), the position of which varies within the mature basidium.
In Cerinomycetaceae, young basidiospores constantly have two nuclei, even when still attached to the basidium (Fig. 4j, n -p, observed also in two additional species not included in phylogenetic analyses).
In Dacrymycetaceae and Dacryonaemataceae, a single nucleus passes into each basidiospore and is not immediately followed by another nuclear division, resulting in young basidiospores that are always uninucleate (Fig. 4k, q -s, u, observed also in c. 15 additional Dacrymycetaceae species not included in the phylogenetic analyses). After discharge, they remain uninucleate until mitosis occurs (including septation), and septate basidiospores contain a single nucleus in each daughter cell v). This has been clearly observed in basidiospores with up to eight cells. When the number of cells is higher, nuclei are obscured by the narrowness of each individual cell. Scattered binucleate, aseptate basidiospores have been seen in samples where several other basidiospores were already septate, probably representing a developmental stage immediately prior to septation.
In Unilacrymaceae, nuclei are very difficult to observe in the basidiospores due to the combination of large volume, high content of lipid bodies, and absence of a colourful cytoplasm (see 'pigments' below). In addition, the differential staining of nuclei with toluidine blue is not as effective as in samples of the other families. Nevertheless, our observations indicate that basidiospores are very likely uninucleate prior to septation (Fig 4t).
Nuclear details are often lost some time after desiccation, except typically for xerotolerant species (e.g., in Dacryonaema), where nuclei can be visible at least in the basidiospores for longer periods after being collected.

Pigments
There is considerable colour variation in the fresh fruit bodies of the four families of the Dacrymycetes. Fruitbodies in the Dacrymycetaceae are often yellow or orange (Fig. 5d), although other colours, e.g., brownish or whitish, occur. In the Cerino mycetaceae, Dacryonaemataceae, and Unilacryma ceae duller colours, e.g., pale cream, ochre, brown or grey, are the rule (Fig. 5a -c), although young fruitbodies of some species may be rather orange. The spore print is often yellow to orange in many species of Dacrymycetaceae with also yellow to orange fruitbodies (Fig. 5e), pale cream or whitish in paler species and in all samples of Cerinomycetaceae, Dacryo naemataceae and Unilacrymaceae from which we were able to obtain a spore print.
Yellow-orange carotenoid pigments were observed mainly in the lipid bodies of the cytoplasm (Fig. 5f, i -l, p -q, s). In some species, these pigments are only clearly observable after having been concentrated by the coalescence of the lipid bodies and are easily overlooked in living or recently dead cells . A distinct carotenoid content was observed in most fresh specimens of Dacrymycetaceae and also in some herbarium samples. In Calocera cornea, the yellowish colour may be particularly inconspicuous and is best seen in the basidium cytoplasm. In Unilacrymaceae and Da. macrosporum, carotenoids are typically not visible (Fig. 5r), but an inconspicuous yellowish cream tinge (never bright yellow or bright orange) after the fusion of lipid bodies may be seen in fruitbodies growing in very bright places and receiving more sun-light. In Cerinomyce taceae, Da. rufum, and Da. macnabbii, yellow or orange colours were not discernible within the cytoplasm and, as noted above, no big oil drops that could concentrate these pigments are normally formed when the cells die ( Fig. 5n-o). These observations do not imply that Cerinomycetaceae, Dacryonaemataceae, and Unilacrymaceae are necessarily devoid of carotenoids, but indicate that the amount of carotenoids, if present, is distinctly lower than in most species of Dacrymycetaceae, or of a different composition. The carotenoids progressively degrade after the samples are collected. For this reason, the yellow or orange colours tend to slowly fade in herbarium specimens and may disappear after a few decades.
The carotenoid contents in the studied species of Dacrymyce taceae show a distinct greenish to greenish grey reaction with Lugol's solution (Fig. 5g, m), although in Calocera cornea it is often inconspicuous (it may be seen on spore prints). This reaction has not been observed in any species of Cerinomyce taceae, Dacryonaemataceae, or Unilacrymaceae (Fig. 5h),  not even when a pale yellowish cream tinge was present in the cytoplasm.
Additional pigments, both intracellular and parietal, are sometimes seen in various species of all four families. An intracellular brownish pigment, seen as small granules in water but partially soluble in KOH, has been noted in samples of all three species of Dacryonaema, and it is particularly conspicuous in some marginal cells of Da. macrosporum (Fig. 5t). Parietal diffuse pigments were observed in Dacrymyces tortus s.lat. and all species of Dacryonaema and Unilacryma, in Dacryopinax ele gans, and in a few Dacrymyces s.lat. species not included in the present phylogenetic analyses. They are more conspicuous in the sterile marginal cells and may vary from clearly brownish to yellowish brown, brownish orange, or reddish brown. In Dacryonaema, they are sometimes better distinguished from the intracellular pigments when mounting samples in KOH (Fig. 5u). These pigments seem to provide at least part of the dark colour in species like Da. rufum, and are more persistent than the caro tenoids in herbarium specimens. All pigments tend to be more abundant in fruitbodies growing in sun-lit places.

Molecular data
Counting each of the three parts of the nrDNA as independent sequences, 402 new sequences were generated in this study ( No introns were detected in the nrDNA or in the ATP6 sequences of any of the included species. Five introns were identified in the amplified RPB1 fragment. All species in Dacrymycetes lack the third RPB1 intron (spliceosomal intron 2 in Matheny et al. 2002), which is consistently present in all species in the outgroup. All species included in the study have a single intron in RPB2 at domain 7 (Arabidopsis intron 16).
In TEF1α, 12 introns were identified in the alignment, seven of them exclusively present in the Tremellomycetes. The second closest intron to the 5'-end is present in some lineages of the

Dacryonaema macnabbii UPS F-940952
Filobasidium wieringae JGI    Dacrymycetes (Fig. 7). We did not find any evidence of heterogeneous nucleotide composition, recombination, or heterotachy that could explain the topologi cal differences. One possibility is that an RPB1 paralogue was sequenced in Dacryonaema. The RPB1 sequence of one sample of Dacryonaema macnabbii (UPS F-940952), amplified and sequenced twice, contained at its end substitutions that resulted in codon positions coding for unique amino acids, and a sequence from another sample (UPS F-940954) had a single base indel towards the end (both strands sequenced and no double signal noted), affecting the open reading frame. These problems were not noted in any of the other Dacrymycetes RPB1 sequences. Excluding suspicious parts of the alignment did not change the topological incongruence of the phylogeny. We consider it unlikely that the incongruence was caused by incomplete lineage sorting, since it affected very deep and comparatively long branches in the backbone, while relationships among closely related species were congruent. No other cases of incongruence were detected.
The final normalised quartet score in ASTRAL III was 0.983 in the five-locus analysis and 0.955 in the nine-locus analysis, indicating very low discordance among the single-locus trees. The trees computed from each of those two analyses were nearly identical. Only three runs reached convergence in the StarBEAST 2 analysis. Nearly every parameter had an ESS > 200 (often > 1 000), including species tree parameters (posterior, likelihood, tree height, tree length), except the cySpe-ciationRate, which did not fully converge (100 < ESS < 200, but see Discussion). RWTY plots did not show convergence problems among the species tree topologies. The ASTRAL III and StarBEAST 2 species trees were very similar, with minor differences only in unsupported groups. Thus, only the species tree from StarBEAST 2 with QS values from ASTRAL and PP values from StarBEAST analyses is shown in Fig. 9.
Several internal subclades can be identified in the Dacrymyceta ceae, and we highlight eight of them with taxonomic or nomenclatural relevance: D1-D8. All are well-supported except D8 that received high support only from SH-aLRT (100 %), PP MrBayes (1.00), and PP BEAST (1.00), and D3 and D5 that received high support in all analyses except ASTRAL III. For details on the included species and morphological characterization of these subclades, see Discussion.

Character evolution
Posterior probabilities of character states (pp) from the SIM-MAP simulations at selected nodes were similar between the MrBayes phylograms and the StarBEAST 2 chronograms, except for the hyphidia ( RayDISC results (Table 6) largely agreed with those from SIMMAP, including the different inferences for the hyphidia when a different tree was used. However, state reconstructions for the basidium type and the amount of lipid content varied considerably at the root node (ancestor of the Dacrymycetes) depending on the assumed distribution of states at that node. A 1-spored basidium was inferred as the ancestral state with 100 % pp assuming an equal distribution of states in MrBayes phylograms, while a 2-spored basidium was inferred with 100 % pp both when the distribution of states is calculated using the estimated transition rates (both trees) and when it is estimated from the data and the model in MrBayes phylograms. The amount of cytoplasmic lipid content is inferred as high when the distribution of states is calculated using the estimated transition rates (96 -98 % pp), while it is inferred as low when it is estimated from the data and the model (100 % pp in both trees) and also assuming an equal distribution of states for the BEAST chronograms (100 % pp).
The median transformation counts (Table 7) show that the 1-spored basidium seems to be a rarely acquired apomorphic state, probably as a reduction of the 2-spored basidium, although there is a substantial uncertainty in the counts from the chronograms. Clamp connections were lost only in a couple of internal groups in the Dacrymycetaceae (remarkably in subclade D8, Fig. 9), and reversions to the state with clamp connections were uncommon. The hyphidium branching seems to have increased in most clades, although the acquisition of branched hyphidia was rare in Dacrymycetaceae. The terminal cells of cortical/ marginal hyphae showed some flexibility, as the transformation counts show numerous transitions from simple and not swollen cells into more or less conspicuously dilated cells in the Dacrymycetaceae. Some reversions to the ancestral state were also inferred. The branched and not swollen cells most probably appeared only once (in the Dacryonaemataceae).
The carotenoids were most likely acquired only once, by the ancestor of the Dacrymycetaceae, with a few secondary losses. The visible lipid content of the cytoplasm in basidia and basidiospores was reduced a few times (e.g., Dacryonaema rufum and some species in the Cerinomycetaceae), although increases in visibility also occurred. The binucleate young basidiospore state is an acquired state in the Cerinomycetaceae.

Dacryonaema and Unilacryma phylogeny and species delimitation analyses
No incongruence among loci were detected. Both ML and Bayesian analyses based on the concatenated dataset resulted in trees with a similar topology. The best partitioning scheme and models for each partition in the ML analysis were: (i) 18S + 5.8S + TEF-1α 2 nd , K2P + I; (ii) ITS1 + RPB1 ( The MrBayes analysis was automatically halted after 8 × 10 6 generations, when the ASDSF dropped below 0.01. At this point, all ESSs after discarding the burn-in were > 1 000, and all tree and model parameters had PSRFs = 1.000 ± 0.001. The STACEY analysis yielded a chronogram with a similar topology to the phylograms obtained using the concatenated dataset. Thus, only the STACEY chronogram with support values from all analyses is shown in Fig. 10. According to SpeciesDe-limitationAnalyzer, the most probable delimitation involves six species (40.7 % posterior probability). All other delimitations had a posterior probability of < 5 %.
Dacryonaemataceae and Unilacrymaceae were resolved as monophyletic, highly supported groups. Within the Dacryonae mataceae, three fully or nearly fully supported main clades are resolved (BS = 100 %, SH-aLRT = 99-100 %, PP MrBayes = 1.00, PP STACEY = 1.00, and the same clades are inferred as putative species (Da. macnabbii, Da. macrosporum, and Da. rufum) in the STACEY species delimitation analysis. The posterior probability among specimens to belong to the same inferred group was high according to the similarity matrix (pp > 0.88). The posterior probability that any of the samples in one of the three putative species may belong to any of the other two putative species was zero.
Within Unilacrymaceae, two fully supported clades were resolved, representing U. bispora and U. unispora s.lat. The posterior probability of all specimens of U. bispora belonging to the same cluster inferred by STACEY was high (pp > 0.99), and the probability that any of the samples belong to any of the U. unispora s.lat. groups was zero. The clade formed by U. uni spora s.lat. samples is further divided into two subclades, one including all North European samples (BS = 71 %, SH-aLRT = 96 %, PP MrBayes = 1.00, PP STACEY = 0.98) and the other including all East Asian samples (BS = 77 %, SH-aLRT = 98 %, PP MrBayes = 0.58, PP STACEY = 0.95). The STACEY analysis inferred these subclades as two putative species. Within the North European group, all specimens had high posterior probability of belonging to the same group (pp > 0.99). Within the East Asian group, the posterior probability of each sample belonging to the same group was highly variable (0.67 < pp < 0.94). The posterior probability of North European and Japanese samples belonging to the same cluster was low (pp < 0.02).

Phenotypic characters
Cytology as a source of taxonomic characters Vacuolar cell contents have been rarely noted in Dacrymycetes by other researchers, and there are scarce records of them in the literature. For instance, they were probably drawn by Nannfeldt (1947) for some basidiospores of Dacryonaema rufum (e.g., his f. 2 showing vacuolated ungerminated basidiospores after 5 days in water). Also, the photograph for explaining the ontogeny of basidia in Dacrymyces stillatus included by Oberwinkler (2014: f. 3 of that study) showed conspicuous VBs in some basidia. Unfortunately, we did not see any taxonomically relevant variation of the VBs in the Dacrymycetes.
The amount of lipid contents, routinely used as a taxonomic character in several groups of Ascomycota (Baral 1992), has been poorly explored in the Basidiomycota. The higher lipid content in the cytoplasm of Dacrymycetaceae and Unilacrymaceae compared to Cerinomycetaceae seems useful for identifying the latter family, although D. tortus s.lat. 1 is a clear exception in which conspicuous lipid drops are easily observed. This species might be confused with others in the Dacrymycetaceae, especially as taxa in the D. tortus complex also have dacrymycetoid fruitbodies. They rarely coalesce and may have obvious orange tinges when young. The basidia are sometimes comparatively long, in which case their dimensions and shape overlap with those of some Dacrymycetaceae.
There are not many studies focused on the nuclear stages in the Dacrymycetes. The first mentions appeared in Dangeard (1895), Istvanffi (1895), Juel (1898), and Maire (1902, who noted that a single nucleus passes into each young basidiospore when attached to the basidium. Bodman (1938) carefully explained the karyological states of basidia and basi diospores of Dacryopinax spathularia, and suggested that of the four postmeiotic nuclei, two migrate to each basidiospore (thus, basidiospores are haploid and uninucleate) and two degenerate in the basidium. Lange (1955), for the Dacrymycetaceae (as 'Dacryomycetace ae'), suggested an ontogeny of the basidia in agreement with Bodman (1938), and depicted dikaryotic vege tative cells and 'oidia-like' mitospores, that might correspond to a particular interpretation of the arthrosporic state of, e.g., D. stillatus. Oberwinkler (2014) described nuclear stages in agreement with Bodman (1938) and Lange (1955) when detailing the general life cycle of Dacrymycetes, exemplified by Dacrymy ces stillatus. Duncan & Chisholm (1982) discovered a third nuclear division (pattern B of Duncan & Galbraith 1972) in the basidia of Dacrymyces stillatus, but basidiospores were still uninucleate (six nuclei remaining in the old basidium). McLaughlin et al. (2016) studied the ontogeny and cytology of a newly described species, Dacryopinax primogenitus, and observed the dikaryotic nature of the vegetative cells and conidia, as well as a similar pattern of post-meiotic nuclear division noted by Duncan & Chisholm (1982).

Most of this information closely agrees with our observations in
Dacrymycetaceae. Although the observation of some nuclear stages was possible without stains, a clear observation of the nuclei often needs specific staining. Most nuclear staining techniques involve a long process (sometimes days) that re quires the fixation of the samples prior to staining, and uses some chemicals unusual in rutine taxonomic studies (Bodman 1938). As an alternative, we tried a fast (< 1 min) staining with alkaline toluidine blue, not requiring previous fixation of the samples (see Material and Methods). This gave reasonably good results, allowing us to confirm several of our observations, despite the fact that only a low number of cells were often properly stained. However, we have not been able to observe the pattern B of the third nuclear division reported by Duncan & Chisholm (1982) in any species.
Our observations of the uninucleate state in basidiospores of Dacrymycetaceae agree with nearly all previous information about Dacrymycetaceae (see above), except for Vandendries (1938), who claimed that basidiospores in Calocera cornea are dikaryotic. Nannfeldt (1947) likely illustrated uninucleate basidiospores in Dacryonaema rufum when still attached to the basidium (f. 2 on p. 328 in that study), although it was not explicitly indicated whether or not the central, rounded organelle depicted was the nucleus. Although further research is needed to confirm the number of nuclei in Unilacrymaceae basidiospores, young (= non-septate) basidiospores in both species seem to be uninucleate (e.g., Fig. 4t). The ancestral state reconstruction infers the uninucleate state of the basidiospores as the most probable one for the ancestors of the Dacrymycetes and also of Dacrymycetaceae, Unilacrymaceae, and the genus Dacryonaema (Table 5 -6, Fig. 9).
Unlike the other families, binucleate young basidiospores seem to be the rule in the Cerinomycetaceae. Maekawa (1987) already indicated that recently discharged basidiospores of Ce. aculeatus were mostly binucleate. In addition to this species, we have consistently observed binucleate young basidiospores in five out of six species of Cerinomycetaceae that were studied based on living specimens (three of them included in the present phylogenetic analyses). Binucleate basidiospores in the Basidiomycota are not uncommon in some taxonomic groups, as noted by Duncan & Galbraith (1972). The binucleate condition is typically a product of postmeiotic nuclear division pattern D, i.e., the original haploid post-meiotic nucleus in the young basidiospore further divides inside the basidiospore and both daughter nuclei remain. Dacrymyces tortus s.lat. 1, in which observations of nuclei are difficult because of high amounts of lipid bodies in the basidiospore cytoplasm, showed both uninucleate and binucleate basidiospores when young and aseptate. This would still agree with pattern D, if we consider the third nuclear division more delayed than in other species, and then very recently discharged basidiospores may be uninucleate for a short period of time. Should further stu dies on a wider number of species confirm our observations, a comparatively long binucleate state of the basidiospores (i.e., not starting only immediately prior to septation) may represent a new synapo-morphy and distinguishing character for the Cerinomycetaceae. With the present data, ancestral state reconstruction analyses indeed infer the binucleate state of the basidiospores as plesiomorphic but unique to the Cerinomycetaceae.

Lineagespecific presence of carotenoids
It is generally accepted that the yellow/orange colours in the Da crymycetes are caused by the presence of various carotenoid pigments (Gill & Steglich 1987, Oberwinkler 2014). The precise composition has been investigated in Calocera, Dacrymyces, and Dacryopinax (e.g., Goodwin 1953, Hanna & Bulat 1953, Fiasson 1973, Czeczuga 1980, and β-carotene has been found to be among the most abundant substances. Carotenoids are lipidophilous and thus typically found in lipid globules or in membranes (Echavarri-Erasun & Johnson 2002). We clearly observed them only in the lipid bodies of the cytoplasm, as reported for many other fungi (e.g., Cederberg & Neujahr 1970, Riley 1978, Griffin 1994). However, Vail & Lilly (1968) reported the highest amount of carotenoids in the 'cell wall fraction' of Dacryopinax spathularia. The colour of fresh fruitbodies in the Dacrymycetes varies considerably, especially at the family level, and species outside the Dacrymycetaceae are commonly creamish, brownish, or greyish, instead of bright yellow or orange ( Fig. 5a-d). Already Nannfeldt (1947) noted that neither the basidia nor the basidiospores of Dacryonaema rufum contain yellow lipid drops, in contrast with their common occurrence in many species of Dacrymycetaceae. Oberwinkler (2014) reported that some taxa, such as Ditiola haasii and Cerinomy ces species were exceptions to the rule that the basidiocarp pigmentation is caused by yellowish to orange carotenoids in the Dacrymycetes. The spore print is often reported as yellow to orange, and this is the case in many species within the Da crymycetaceae (e.g., Calocera viscosa, Fig. 5e, Dacrymyces aureosporus, D. chrysospermus), but there are also reports of species with pale or whitish spore prints (e.g., Cerinomyces albosporus, Boidin & Gilles 1986). Other published information may not be as accurate. For in stance, some taxa with clearly yellowish spore prints have been incorrectly reported as having a whitish spore print (e.g., the generic descriptions of Dacrymyces, Ditiola, Femsjonia and Guepiniopsis in Laessøe & Petersen 2014).
Our results indicate that there has been a sharp increase in the amount of carotenoids in the Dacrymycetaceae, whereas these pigments are not clearly observable, perhaps even absent, in Cerinomycetaceae, Unilacrymaceae, and in Dacryonaema. Also, the colour of the fruitbodies and that of colonies grown in pure culture seem to be related. Shirouzu et al. (2009) only reported yellow to orange colours for species in the Dacrymy cetaceae (pale or whitish colonies also occurred), while the colonies in the clades corresponding to Cerinomycetaceae and Unilacrymaceae were constantly white.
Carotenoids have been demonstrated to play a role in the cell protection against oxidative stress and photoreactions (Will et al. 1984, Will III & Scovel 1989, Moliné et al. 2010, Avalos & Limón 2015. This function was demonstrated in pure cultures of Dacryopinax spathularia (Goldstrohm & Lilly 1965). However, several of the species most frequently found on exposed wood, receiving large amounts of sun-light, possess small or no amounts of carotenoids, e.g., Dacryonaema, Unilacryma, and at least one Cerinomyces species (D. tortus s.lat. 3). Therefore, a putative ecological advantage of these pigments may also be related to other functions of carotenoids (e.g., as precursors of some active compounds, Avalos & Limón 2015).
Iodine solutions are known to react with the carotenoid contents in different groups of fungi, typically exemplified by the greening of the paraphyses of several ascomycetes with Lugol's or Melzer solutions (Malenҫon & Bertault 1971, Schumacher 1990, other coloured reactions also occur in this group). The carotenoid contents in the Dacrymycetaceae display a greenish to greenish grey reaction in Lugol's solution, while this reaction has not been observed in any species of Cerinomycetaceae, Dacryonaema, or Unilacryma. The combination of direct observation of the fruitbodies and cell cytoplasm, plus their greenish reaction with iodine, may allow researchers to quickly and objectively assess the presence and relative amount of carotenoids.
There are, however, some species in the Dacrymycetaceae with colours other than the common yellow or orange, such as the whitish or very pale coloured Calocera cavarae, Ca. pallido spathulata (Reid 1974), Ca. fusca, Ditiola haasii (Oberwinkler 1989), and Dacrymyces cylindricus (Shirouzu et al. 2017). For these taxa, it is reasonable to suggest a secondary loss of the pigmentation, since DNA sequence data have demonstrated that they belong to at least three different internal clades of the Dacrymycetaceae dominated by species with yellow-orange pigments (Shirouzu et al. 2013(Shirouzu et al. , 2017. This agrees with both the present phylogenies and our ancestral state reconstructions including Ca. fusca and D. cylindricus (Fig. 9). The loss of carotenoids in certain species or specimens has already been reported in other fungi (e.g., Baral 2004 in Sarcoscy pha, Hansen et al. 2013 in some clades of Pyronemataceae, Olariaga et al. 2015 in Cantharellus). In addition, fruitbodies of species that normally produce carotenoids often produce more pigment when exposed to sunlight. In contrast, when they grow in absence of light they are often paler or even whitish, a phenomenon well known in some other Basidiomycota (Pilz et al. 2003, Olariaga et al. 2015. The Dacrymycetes are not an exception. Such uncommon cases do not invalidate the usefulness of carotenoid content for taxonomic purposes. With the current data, any species with a conspicuous and high carotenoid content would be tentatively best placed in the Dacrymycetaceae, since the presence of distinct carotenoids is inferred as being acquired in the ancestor of the family, and few reversions to the plesiomorphic carotenoid-deficient condition occurred (Table 5 -7). We recommend that future studies on Dacrymycetes include careful observations of fresh, unstained samples and the subsequent reporting of approximate carotenoid content.

Other morphological characters
There is only one previous study including ancestral reconstruction analyses of multiple morphological characters in the Dacrymycetes (Shirouzu et al. 2013). They proposed that the macromorphology of the fruitbodies is highly homoplastic and rarely a synapomorphy for any higher clade, and the same was true for the wall thickness of marginal hyphae and, to some extent, also basidiospore septation. By contrast, other characters, such as the presence or absence of clamp connections, turned out to characterise some phylogenetic lineages.
In agreement with Shirouzu et al. (2013), our phylogeny shows that clamped septa seem to be the plesiomorphic state in the Dacrymycetes, while simple septa appear in some derived clades only within the Dacrymycetaceae. In addition, we found that clamp connection morphology and location may be informative for high-level taxa delimitation. Loop-like, long and angulose clamp connections (i.e., long medallions), almost restricted to subhymenial hyphae and the bases of basidia and hyphidia, are characteristic of the three known species of Dacryonaema. This specific morphology and distribution do not consistently appear in any other group, as was inadvertently indicated by Reid (1974) when he commented on the uniqueness of the Da. macnabbii (as Dacrymyces macnabbii) clamps. Loop-like, but normally shorter and more rounded clamp connections are present in some Dacrymycetes, as noted by, e.g., Martin (1949). We have observed them in some species, for instance in the Dacrymyces tortus complex and in Unilacryma, and infrequent transitions to Dacryonaema-like, longer clamp connections may occur. It is interesting to note that Unilacrymaceae, the sister group to Dacryonaema, also lacks clamp connections at some septa, but this feature may occur in some species of Cerinomy cetaceae and Dacrymycetaceae with typically clamped septa, due to the formation of secondary simple septa.
The morphology of cortical/marginal hyphae was used as a key character by Oberwinkler (1994Oberwinkler ( , 2014 to circumscribe some genera, and it is commonly found in the literature as important to characterise, e.g., Guepiniopsis and Heterotextus (McNabb 1965c-d). Our ancestral state reconstruction analyses show that the presence of branched and anastomosing terminal cells with walls that become gelatinised in KOH is plesiomorphic and synapomorphic for Dacryonaema but absent in other clades, where terminal cells of cortical/marginal hyphae tend to be simple and anastomoses rare. Indeed, the ancestral condition predicted for the Dacrymycetes, Cerinomycetaceae, Dacry mycetaceae, and Unilacrymaceae is the presence of simple and not swollen terminal cells, although some uncertainty exists depending on the assumptions made at the root node (Table 5 -6). These terminal cells are conspicuously dilated and morphologically diverse in some groups (e.g., Femsjonia, Heterotextus), suggesting that this character should be carefully addressed in future studies, particularly in the Dacrymyceta ceae.
The type of basidium, i.e., forked or simple, cannot any longer be considered a good character to distinguish Unilacryma from other clades, as Unilacryma bispora constantly has bisterigmate basidia. Simple, 1-spored basidia are highly unusual in the Basidiomycota, and in most cases species characterised by 1-spored basidia belong to genera with a variable number of spores produced per basidium, e.g., Heteroacanthella, Lactarius, Oliveonia, Russula, Septobasidium, and Tremella s.lat. (Zamora et al. 2014, 2016, Vidal et al. 2019. This result, together with the phylogenetic position of Dacryonaema, allowed inferring the bisterigmate basidium as the most probable ancestral state for the Dacrymycetes in most analyses. However, the rayDisc analysis assuming equal distribution of states at the root node unexpectedly reconstructed the 1-spored basidium as the ancestral state for the Dacrymycetes. This is the least parsimonious option and also the equally-distributed states assumption is anyway unrealistic. McNabb (1973) emphasised the degree of branching of the hyphidia for delimiting species within Dacrymyces s.lat., and this is still used to characterise some species, e.g., Dacry myces dendrocalami (Oberwinkler & Tschen 1989, Shirouzu et al. 2009). In our analyses, the morphology of the hyphidia has been simplified to only two extremes (absent or simple vs conspicuous and branched) to be able to make use of imprecise descriptions from the literature. For instance, McNabb (1973) commented on the presence (sometimes occasional) of simple hyphidia in the majority of Dacrymycetes, in which branched hyphidia are not found, while Oberwinkler (2014) reported that the hymenium is composed of only basidia, a few species excepted. This character should be revisited in the future, since there are some clades, e.g., the one including D. chrysocomus and D. ovisporus, where hyphidia are often distinct, yet simple. In others, e.g., in Calocera s.str., and the D. chrysospermus and D. stillatus clades, true hyphidia are probably absent and the observed hyphae tend to develop into new basidia and new gene rative hyphae. While the plesiomorphic state in Dacryona ema and Unilacryma seems to be the branched hyphidia, and the ancestral state for Dacrymycetaceae the simple or absent hyphi dia, important differences are found between the ancestors of the Dacrymycetes and Cerinomycetaceae. These differences are related to both the type of phylogenetic trees and the method (Table 5 -6), a relatively unusual result (Litsios & Salamin 2012, Cusimano & Renner 2014) that should be considered before drawing definitive conclusions. Consequently, given the current data, we consider the ancestral state of the hyphidia in the Dacrymycetes and Cerinomycetaceae to be uncertain.
There are additional characters that should be investigated in the future. Reid (1974) found that the basidium in Calocera furcata usually has an apical protuberance that is normally absent in C. cornea. This character was also noted by Van de Put (2014). The morphology of mitotic spores, whether dikaryotic arthrospores, other conidia formed from vegetative cells, or microconidia formed from the basidiospores, seems to be useful for species delimitation. Shirouzu et al. (2009) reported clear differences in mitospore morphology in fruitbodies of D. pinacearum, D. stillatus, and D. subarcticus, and Kirschner & Yang (2006) described morphologically characteristic stauroconidia in Dacryoscyphus chrysochilus. Reid (1974) described Calocera pallidospathulata by emphasizing the deviant cylindrical microconidia produced by mature basidiospores when compared with other species of Calocera, which normally produce rounded microconidia. These traits and others (e.g., characteristics in culture) are promising for delinea ting taxonomic groups, especially within the Dacrymycetaceae.

Differences among phylogenetic analyses
Coalescence-based phylogenetic inference often outperforms concatenation when discordance among single-locus phylogenies exists, assuming that the conflicts are a consequence of incomplete lineage sorting (Pamilo & Nei 1988, Kubatko & Degnan 2007, Degnan & Rosenberg 2009, Fernández-Mazuecos et al. 2018. Precise estimates of the effective population sizes, through inclusion of several individuals per species, benefit the species tree inference (Heled & Drummond 2010), while the superior accuracy of coalescence vs concatenation is often only relevant for relatively short branches (Ogilvie et al. 2016). With our final dataset containing no incongruence after the exclusion of Dacryonaema RPB1 sequences, less than ten independent loci, and some species being represented by a single individuum, both analytical approaches supported approximately the same groups. The few alternative relationships among some species were unsupported and affected very short branches, particularly in the Dacrymycetaceae. On the other hand, despite the improvements in performance in StarBEAST 2, species tree inferences often require very long runs for the MCMC to converge (Ogilvie et al. 2017). Our RWTY analyses indicate that species tree topologies did indeed reach convergence, and we only noted some convergence problems with the speciation rate parameter (cySpeciationRate, with 100 < ESS < 200), an observation made previously (Dorsey et al. 2018).
Phylogram and chronogram topologies obtained in the present study are generally congruent with the most recent phylogenies based only on data from the nrDNA (e.g., Shirouzu et al. 2016Shirouzu et al. , 2017, at least concerning the supported groups, and slightly deviate from previous studies including fewer taxa (Shirouzu et al. 2013). Our improved resolution was probably achieved by the combination of more taxa and DNA loci.

Taxonomic boundaries at ordinal and familiar ranks
Nowadays, there are two main currents for describing higher level taxa. One is the most classical, focused on the diagnosability of the resulting named groups (Christenhusz et al. 2015) considering their evolutionary relationships (thus, in agreement with phylogenies based on molecular data), as followed by, e.g., the consortium of experts in The Angiosperm Phylogeny Group (2016). The alternative method is 'temporal banding', which considers estimated stem ages of the clades to define taxonomic ranks (Avise & Johns 1999, Holt & Jønsson 2014. However, adhering too tightly to the temporal banding method may have undesirable taxonomic and nomenclatural consequences. For instance, it may force the joining, under the same higher-ranked taxon name, of lower-ranked taxa that are strikingly different based on phenotypic characters, thus resulting in heterogeneous assemblages that are difficult to recognise. Conversely, it could force the splitting of morphologically homogenous groups for little gain (Lücking 2019).We think that trying to reconcile the systems may be the best solution, i.e., proposing only diagnosable groups while considering their stem ages with some flexibility. As guides for our study, the mean stem age estimates proposed by Kraichak et al. (2018) is in the range 176 -194 Ma for orders and 111-135 Ma for families in the Lecanoro mycetidae and Ostropomycetidae (Ascomycota). Zhao et al. (2017), using a wide sampling in the Basidiomycota, reported stem ages of most orders around 120 -290 Ma. Although Zhao et al. (2017) were not focused on the family level, the oldest crown ages estimated for families were c. 150 -160 Ma (e.g., in Moniliellaceae and Phleogenaceae), meaning that their stem ages must be even older and broadly overlapping with the order ages. Considering these premises, we can circumscribe at least five groups within the Dacrymycetes that deserve a taxonomic rank above the generic level, agreeing with the current use of ranks in the Basidiomycota. Four of them are well characterised by phenotypic characters and highly supported monophyletic groups in all analyses, three of which are already named as families, viz., Cerinomycetaceae, Dacrymycetaceae, and Unila crymaceae. The specimen NBRC 110597 may represent another family judging from its phylogenetic position, but since we lack morphological data and its phylogenetic placement is based only on nrDNA data from a single sample, we consider it appropriate to leave it unnamed for the time being.
The clade containing the genus Dacryonaema includes taxa with long medallion clamps concentrated in the subhymenium and the bases of basidia and hyphidia, with septa normally lacking clamp connections elsewhere. The cortical/marginal hyphae are variable but some terminal cells are always abundantly branched and anastomosed, becoming gelatinised in KOH. Hyphidia are branched and carotenoids are inconspicuous or absent. Basidiospores are ± cylindric-allantoid and may remain aseptate or become transversely 1-septate (but rarely up to 3-septate) at maturity. In addition, all TEF-1α sequences have an intron (the sixth in the analysed fragment) that is absent in other samples among the Dacrymycetes. This clade is proposed here as a new family, Dacryonaemataceae, since its morphological characteristics do not allow the placement of Dacryonaema in any of the previously circumscribed families. This decision is also in agreement with the relatively old stem age of the clade, as our results indicate that the ancestors of Dacryonaemataceae and Unilacrymaceae may have split around 169 Ma. This estimate was calculated setting an age for the Dacrymycetes within the ranges recently obtained by He et al. (2019) and Varga et al. (2019).
The clade representing the family Unilacrymaceae, which only contains the genus Unilacryma, is comprised of species with short clamp connections at most septa and often simple and not swollen terminal cells of marginal hyphae that are not highly gelatinised in KOH. Hyphidia are branched and there is a high lipid content in the cytoplasm of the basidia and basidiospores. Carotenoids are inconspicuous or absent. The basidiospores are ovoid to ellipsoid, presumably uninucleate when young, and become muriform at maturity. The discovery of the new species U. bispora demonstrates that simple basidia, one of the main morphological characters used to recognise the family by Shirouzu et al. (2013), is not a synapomorphy of this group. On the other hand, the presence of ovoid-ellipsoid basidiospores that become muriform, another important character noted in Shirouzu et al. (2013), is also not enough to distinguish the family, since this morphology is present in Dacrymyces ovispo rus (Fig. 5s). This species, however, is a typical member of the Dacrymycetaceae and is characterised by the conspicuously high carotenoid cell content, and is additionally distinguished from Unilacryma by the simple to (rarely) sparingly branched hyphidia (Fig. 2l). In addition, D. chrysocomus, a close relative to D. ovisporus in our phylogenetic analyses, has predominantly cylindric-allantoid basidiospores that frequently become muriform at maturity. Scattered shorter, ellipsoid basidiospores may be present, linking the morphology of the cylindric-allantoid basidiospores with the ovoid ones.
The clade representing the family Cerinomycetaceae contains species with short clamp connections at most or all septa and often simple and not swollen terminal cells of cortical/marginal hyphae that are not clearly gelatinised in KOH. The hyphidia are often abundantly to sparingly branched, but sometimes simple or absent. The amount of lipid bodies in the cytoplasm of basidia and basidiospores is normally low (with one exception), and carotenoids are inconspicuous or absent. The basidiospores are ± cylindric-allantoid, binucleate during some time before the formation of the first septum, and may remain aseptate or become transversely 1-3-septate at maturity. Basi dia tend to be relatively short compared to those of the other clades, especially to Dacryonaema and Unilacryma. The fruitbody macromorphology varies from resupinate, coalescing to dacrymycetoid, and the consistency is ± dry, fleshy-cartilaginous or gelatinous.
The clade representing the family Dacrymycetaceae includes several genera, e.g., Calocera, Dacrymyces, Dacryopinax, Fem sjonia, and Heterotextus. Some of these genera are non-monophyletic, which has also been demonstrated in previous studies (Shirouzu et al. 2007(Shirouzu et al. , 2009(Shirouzu et al. , 2013. This clade is very species-rich and morphologically diverse, containing species with both clampless and clamped septa, variously shaped terminal cells of cortical/marginal hyphae. Basidiospore shape and septation are highly variable, although basidiospores tend to be cylindric-allantoid in most cases. All species studied possess basidia and basidiospore cytoplasm with high lipid contents and young basidiospores with a single nucleus before the formation of the first septum. The presence of conspicuous carotenoids and the absence of highly branched hyphidia are characteristic of most species, but exceptions occur. All TEF-1α sequences obtained for this family lacked the third intron in the analysed fragment. Taxa at the order level are not well characterised, and some doubts arise when assessing the placement of Dacryonae mataceae. A conservative classification that would agree with the tree topology and partly with the estimated node ages, and in addition does not require any nomenclatural changes, would be to adopt the two orders proposed by Shirouzu et al. (2013) and include Dacryonaemataceae in the order Unilacry males. However, since the 1-spored basidium is no longer a synapomorphy of the Unilacrymales, and since basidiospore shape would also not be diagnostic with the putative inclusion of Dacryonaemataceae, the phenotypic characterization, i.e., diagnosability, of Unilacrymales becomes problematic. Some characters are shared by the Dacryonaemataceae and Uni lacrymaceae, e.g., the clamped hyphae, branched hyphidia, absence of conspi cuous carotenoids, and the similar ecology, but those traits are to some extent shared with members of the Cerinomyceta ceae. The amount of observable lipids is not diagnostic, as no lipid drops are seen in the cytoplasm of Da. rufum. As noted above, basidia are often much larger in Dacryonaemataceae and Unilacrymaceae when compared with Cerinomycetaceae, but in the Dacrymycetaceae some species also have relatively large basidia. The differences between Ce rinomycetaceae and Dacrymycetaceae are also considerable, making it a challenge to identify diagnostic characters that can include all members of both families while at the same time excluding all taxa in Dacryonaemataceae and Unilacrymaceae. Consequently, sorting the four named families into two orders seems unpractical.
There are two alternatives instead of accepting the two currently recognised orders in the Dacrymycetes. One would be to create a new order for Dacryonaemataceae, and the second option to accept the four families (putatively five) within a single order Dacrymycetales. Erecting an order for Dacryonaemataceae might be justifiable from both a morphological point of view and by using the divergence times from Zhao et al. (2017). However, it would then be logical to also include the Cerinomycetaceae in a separate order, resulting in the recognition of four orders, with a single family in each one. We refrain from taking this action in the present study, as it would currently not provide any substantial benefit. The recognition of the four described families in a single order Dacrymycetales is the other conservative approach that does not require any nomenclatural change, and seems reasonable based on morphological and molecular data. Considering the ample overlap among stem and crown ages of orders and families in Zhao et al. (2017), and the uncertainty of our age estimates, we believe that this provisional decision does not disagree with the temporal banding method. Further studies based on genome-scale data could provide a more definite classification for the higher ranked taxa, as well as more accurate age estimates.

Taxonomy below the family level and future directions
Relationships among species groups in the Dacrymycetaceae are not well resolved, and generic delimitations and species circumscriptions need substantial efforts to reach taxonomic and nomenclatural stability. Subclade D1 includes D. chrysocomus s.lat., D. ovisporus, and D. variisporus, and is characterised by the presence of clamp connections, thick-walled terminal cells of cortical/marginal hyphae, which may form a cortex concolorous with the hymenium, and by the large basidiospores that are either transversely multiseptate and/or muriform at maturity. Subclade D2 corresponds to Heterotextus (H. miltinus s.lat., generic type) and has clamp connections, a cortex of thickwalled dilated terminal cells, concolorous with the hymenium, and multiseptate but not muriform, basidiospores. Subclade D3 is formed by a group of sessile to stalked Dacrymyces s.lat. with clamp connections and 0 -3-septate, thin-walled basidiospores (D. ancyleus, D. citrinus, D. cylindricus, and D. cyrtosporus). Subclade D4 contains a group of stalked species of Calocera s.lat., Dacrymyces s.lat. and Dacryopinax s.lat. with clamp connections and thin-walled, 0 -3-septate basidiospores (Calocera guepinioides s.lat., Ca. fusca, D. fla belliformis, and Dacryopinax sphenocarpa). Subclade D5 is a morphologically heterogeneous group with sessile and stalked species, with and without clamp connections, with or without branched hyphidia and with thin-to slightly thick-walled, 0 -3-septate basidiospores, including Ca. arborea, D. den drocalami, and D. cf. adpressus. Subclade D6 corresponds to the genus Femsjonia (F. monospora and F. peziziformis s.lat., generic type), which has clamp connections, a contrasting pale cortex of thin-to thick-walled, cylindrical or ± dilated terminal cells, and 1-multiseptate, thin-walled basidiospores. Subclade D7 includes specimens provisionally identified as D. minutus, with clamp connections, a cortex of thick-walled and ± dilated terminal cells, concolorous with the hymenium, and 3-septate basidiospores. Finally, subclade D8, the 'core' Dacrymyceta ceae, includes most species without clamp connections (with the exception of D. pachysporus), and the generic types of Calocera (Ca. viscosa), Dacrymyces (D. stillatus), Dacryomitra (Dacryomitra pusilla), Dacryopinax (Dacryopinax elegans), Dacryoscyphus (Dacryoscyphus chrysochilus), Ditiola (Di. radi cata), and Guepi niopsis (G. buccina). The macromorphology in this large last group is highly variable, with fruitbodies that may be resupinate, clavarioid, cochleariform, flabelliform, pulvinate, etc., and with thin-to thick-walled basidiospores displaying various septation patterns. While some groups in subclade D8 are morphologically homogeneous, e.g., Cerinomyces without clamp connections, or the D. chrysospermus group with multiseptate basidiospores, there are also exceptions. For example, D. novaezelandiae s.lat. and D. sanaugustinii, which also have multiseptate basidiospores, are phylogenetically related to Dacryopinax elegans, a species with 3-septate, thick-walled basidiospores and a very different macromorphology. The lack of support for some of these groups is likely worsened by the limited amount of available DNA data. With the primers and PCR programmes provided in the present study, it should be unproblematic to obtain at least three or four unlinked DNA loci from recently collected (≤ 20 years old) samples.
The inclusion of several new samples in our Dacrymycetes phylogeny has led us to detect some misidentifications that have taxonomic and nomenclatural relevance. For example, D. stil latus and Di. radicata, the types of two accepted genera, have been variously interpreted in the literature. Dacrymyces stillatus is likely a species complex, and an improved sampling of taxa and markers is needed to clarify its identity. The specimens provisionally identified here as D. stillatus s.str. agree well with the protologue and with the current European concept of the species, which is a taxon very commonly found in the Nordic countries. This species has clear molecular differences with the 'D. stillatus' sample used by Shirouzu et al. (2013), and also a somewhat different morphology according to the published description and iconography. Ditiola radicata is a common species in Europe, and the Swedish and Finnish specimens included in our phylogeny agree very well with the general concept of this species. However, this understanding of the species is notably different from 'Di. radicata' sensu Shirouzu et al. (2013), obtained from culture HHB-12010, which nests in the same clade as specimens previously identified as D. capi tatus and D. lacrymalis. This clade is probably another species complex, within which D. capitatus s.lat. displays substantial morphological variation and misidentifications are likely. The sequences retrieved from the JGI as 'Calocera cornea' were also substantially different from those of our samples of Ca. cornea and agree well with preliminary data we have obtained from our own specimens of Ca. furcata. The present study is not aimed at resolving the taxonomy of these taxa, but the disagreements mentioned here serve to caution against uncritical taxonomic rearrangements and nomenclatural novelties based on previously generated data.

Species delimitations in Dacryonaema and Unilacryma
According to our results, the five species of Dacryonaema and Unilacryma found in the Nordic countries seem to be well differentiated, and distinguishing them by morphological or molecular data (or both) is unproblematic. The three species of Dacryo naema accepted here are highly supported by the STACEY species delimitation analysis, receive high support in all other phylogenetic analyses (Fig. 10), and are distinguished by several morphological characters (see Taxonomic part). Therefore, we do not consider them mere infraspecific divergent populations (Sukumaran & Knowles 2017). In addition, distribution ranges overlap at least in northern Europe, yet the three species retain their morphological distinguishing features in sympatry.
Unilacryma bispora is an easily circumscribed species, with little taxonomic problems except for its potential confusion with Dacrymyces ovisporus (see Taxonomic part). The STACEY analysis supported this as a distinct species, and it appears monophyletic in all phylogenetic analyses. Furthermore, basidia and basidiospore morphology (Fig. 6, 14 -15) clearly separate it from its closest relative, U. unispora. Species delimitations in the Unilacryma unispora s.lat. clade are much less clear. While, as a group, it is distinguished from U. bispora, it is not entirely clear whether more than one species should be recognised here. As noted in the results, the U. unispora s.lat. clade is formed by two partially supported subclades, one joining all North European specimen and the other comprising the two Japanese specimens and additional environmental sequences from eastern Asia. It should be noted that the only molecular difference between the two clades shared by all included samples are two transitions (C/T) in the 28S dataset. Additional molecular differences are present between the 18S, ITS, and RPB2 sequences of the Japanese specimens and those from the Nordic countries. Even if the two clades are distinguished by STACEY as two putative species, this pattern can easily be obtained when analysing few samples of isolated populations in species with a broad distribution (Sukumaran & Knowles 2017, Chambers & Hillis 2019. Also, the few differences detected might just be part of their intraspecific variation. Furthermore, the clade formed by the Japanese samples is poorly supported in the concatenated Bayesian analysis, although it is well-supported in the ML analysis. Since we have no information about other populations of U. unispora s.lat. in Eurasia and North America, and the amount of missing data from the East Asian clade is high due to the single region obtained from environmental studies, keeping all of them temporarily as a single species is the only justifiable option with the present data. A phylogeographic study and a careful morphological analysis involving numerous vouchered samples from throughout the northern hemisphere is needed to decide whether U. unispora is a single species or a complex of closely related species.

TAXONOMY
We present the following provisional key to families in the class Dacrymycetes: Fruitbodies ± firmly gelatinous when fresh, sometimes cartilaginous, xerotolerant, sessile or stalked, with the fertile parts pulvinate to depressed or subglobose. Hymenium amphigenous, or confined to the upper part of the fruitbody when fruitbodies are discoid, then with a distinct sterile cortex. Cortical/marginal hyphae variable, frequently clampless, but sparse loop-like clamp connections may occur, some with terminal cells clearly branched and anastomosing, strongly gelatinised in KOH; unbranched cylindrical, thick-walled and slightly dilated terminal cells may also occur, then normally tapering at the apex. Internal hyphae thick-walled, clampless or with very sparse clamps. Subhymenium composed by more or less thin-walled hyphae, with long and angulose, medallion clamps. Basidia 2-spored, bearing a loop-like basal clamp connection; apex frequently U-shaped, sometimes Y-shaped, and more rarely with an apical protuberance. Hyphidia abundantly branched and often anastomosing, with a loop-like basal clamp connection. Recently discharged and still aseptate basidiospores uninucleate. Mature basidiospores 0 -1-septate (exceptionally 3-septate), thin-walled, hyaline, cylindric-allantoid. Spore print whitish. Microconidia ellipsoid. Cell cytoplasm without conspicuous carotenoids under the light microscope. Brownish intracellular and diffuse parietal pigments often visible at least in some cortical/marginal hyphae. Basidia and basidiospores with or without conspicuous lipid drops when old.

Key to species
1. Fruitbodies for a long time appearing as a conical, dark brown to dark reddish brown (rarely amber-coloured), sterile stalk, which finally develops an apical, often rounded or irregularly shaped, but not centrally depressed, pruinose fertile head, ≤ 0.7 mm diam when fresh and mature; dead basidia and basidiospores with a rather even cytoplasm, with invisible or hardly discernible lipid bodies; basidiospores (2.9 -)3.2 -4.3(-4.9) μm wide . . . . . . . . . . . . . . Da. rufum 1. Fruitbodies very variable in colour, from yellowish or cream coloured to brownish or dark brown, pulvinate to pustulate at first, later slightly convex or becoming flattened or centrally depressed (pezizoid), sessile or sometimes shortly stalked, often ≥ 0.7 mm diam when fresh and mature; dead basidia and basidiospores with a slightly to conspicuously granular cytoplasm, often with some discernible lipid bodies; basidiospores (3.5 -) 4.3 -6.6(-7. Basidiocarps ± firmly-gelatinous, from pale cream coloured to yellowish, ochraceous, or brownish, sometimes with an orange tinge, darker when dried, xerotolerant, at first erumpent or pulvinate, soon discoid, flattened to centrally depressed (pezizoid), sessile or shortly stalked, attached to the substrate by a central point and sometimes rooting, (0.5 -)0.7-1.5(-1.8) mm diam when mature. Cortical/marginal hyphae conspicuous, variable, at least some with terminal cells abundantly branched and partially anastomosing, with walls that become gelatinised in KOH; some slightly to clearly thick-walled and somewhat dilated cells may also appear. Internal hyphae (2 -) 2.5 -7 µm wide, slightly to clearly thick-walled, often roughened, embedded in a gelatinous matrix, clampless or, rarely, with scattered and obscure long clamps. Hymenium usually confined to the upper surface of the fruitbodies, rarely becoming amphigenous with age, forming numerous basidia intermixed with hyphidia. Subhymenium with thin-to slightly thick-walled hyphae, with long and ± angulose medallions, including at the bases of basidia and hyphidia. Hyphidia welldeveloped, profusely branched, normally abundant, 2-4(-5) μm thick in the upper half, terminal branches 1-2(-2.5) μm thick, often with simple septa throughout their length. Young basidia cylindrical to very narrowly clavate, sometimes with a tortuous base. Mature basidia 2-spored (very rarely 3-spored basidia have been seen), (42 -)45 -100(-125) × (4 -)4.5 -7(-10) μm, sterigmata 11-20 μm long, often forming an obtuse angle or arising from slightly lateral position from the apex, basidia sometimes with a small apical papilla between the sterigmata. Basidiospores (8.7-)10.4-13.9(-15.4) × (3.5-)4.3-5.5(-6.2), Q = (1.6-)2 -3.2(-3.7), ellipsoid to cylindric-allantoid. Basidia and basidiospore cytoplasm hyaline, often slightly granular in KOH due to the presence of some lipid bodies, which may form small lipid drops after coalescence. Microconidia from germinating basidiospores rarely seen, hyaline, narrowly ellipsoid to cylindric, about 3.5 -4.5 × 1-1.5 μm. Ecology & Distribution -Apparently common on exposed coniferous wood, particularly Pinus sylvestris branches still hanging from the tree, but also found on fallen wood lying on the ground. Its distribution is uncertain because it has been confused with Dacrymyces tortus and other species, but published records and studied specimens demonstrate that it is present in the temperate and boreal biomes of north-western Europe. We have studied specimens of Da. macnabbii from Finland, Norway, Sweden, and the United Kingdom. It has also been reported from Denmark by Mathiesen (1991), whose descriptions and iconography agree well with our observations. The only two Swedish specimens previously identified as Da. macnabbii (S 250593 and UPS F-561899) were misidentified and belong to the D. tortus group.
Nomenclature & Taxonomy -The yellowish to brownish colour of the fruitbodies, as well as the abundantly branched hyphidia, the long loop-like clamp connections present at the bases of basidia, and the basidiospore dimensions, as already reported in the protologue (Reid 1974), easily distinguish this species. The holotype specimen is particularly well distinguished from the morphologically close Da. macrosporum thanks to its clearly ochraceous colour and the short basidiospores (10.5 -13 μm as indicated by D.A. Reid). In addition, all sequenced samples accepted here as Da. macnabbii form a well-supported clade including Reid's type collection, clearly independent from the clade formed by all Da. macrosporum samples. Some species in the D. tortus group may be macromorphologically very similar to both Da. macnabbii and Da. macro sporum, and a microscopical study is often needed to separate them. The cortical/marginal hyphae of D. tortus s.lat. are not finely and abundantly branched, their walls barely gelatinise in KOH (they may show a somewhat roughened surface), the hyphidia are less abundantly branched, basidia are smaller (particularly shorter), and clamp connections are much more abundant in all hyphae and not of the form of long or angulose medallions (short medallions may be present particularly on internal hyphae). Diagnosis. Morphologically similar to Dacryonaema macnabbii, but fruitbodies frequently darker (brown to dark reddish brown or blackish), cell cytoplasm with more distinct lipid bodies (higher lipid content), and basidiospores distinctly larger, (12.9 -)14.7-18.6(-23) × (4.4 -)5.2 -6.6(-7.1) μm (on average > 14.5 × 5.3 μm).
Ecology & Distribution -Insufficiently known; it seems to be much rarer than Da. macnabbii and Da. rufum. All studied specimens were growing on rather weathered and decomposed Pinus sylvestris wood, on branches lying on the ground, where it may grow intermixed with species in the D. tortus group (e.g., in Fig. 11a, the slightly paler fruitbody partially shown at the bottom of the picture is D. tortus s.lat. 3).
Nomenclature & Taxonomy -This species is morphologically distinguished from Da. macnabbii by its darker colour and, in particular, by the larger basidiospores. However, since the des cription above is based on only c. 250 fruitbodies (c. 150 from the holotype and isotype) from six individuals, its intraspecific variation is insufficiently known, so the traits selected as diagnostic characters should be taken with some caution. Lightcoloured specimens of Da. macrosporum and dark coloured specimens of Da. macnabbii may be very difficult to separate, unless they are fertile and a certain number of well-developed basidiospores could be measured. Dacryonaema fruitbodies are comparatively long-lived and environmental conditions may affect basidiospore size, particularly if the maturation is aborted when the fruitbodies suddenly dry up. In addition, basidia in Da. macrosporum tend to be slightly wider and the amount of visible lipids in the basidia and basidiospore cytoplasm is slightly higher. Unfortunately, those characters show a considerable overlap and, rarely, specimens of Da. macnabbii with unusually wide basidia and a considerable amount of lipid bodies occur. In any case, the two species are clearly separated by DNA sequence data. For the distinction between Da. macrosporum and species in the D. tortus group, see taxonomic comments under Da. macnabbii. The two isotypes are part of the same individual, while the holotype is a different individual, collected on a different piece of wood. However, the three specimens are to be considered as duplicates because they are part of the same gathering (ICN, Art. 8.2, Turland et al. 2018 Basidiocarps firm-gelatinous to cartilaginous, xerotolerant, often consisting on a well-developed (but sometimes very short), reddish brown to dark reddish brown, narrowly conical stalk, 0.3 -2(-2.5) mm long, 0.3 -0.7(-0.8) mm wide at the base and tapering to around 0.2 -0.4 mm wide near the apex. Tardily, an apical, concolorous or slightly paler, rounded, fertile head is developed, becoming pruinose and thus looking distinctly paler when fully mature, 0.3 -0.8 mm diam. Fruitbodies sometimes orangish, amber-coloured or with a pinkish hue when growing in the dark or when developing and maturing quickly, distinctly shorter in the last case. External hyphae of the stalk typically thick-walled, strongly glued to gether and difficult to dissect; cortical/marginal hyphae variable, at least some with thin-walled, anastomosing and branched terminal cells, gelatinising in KOH, others thin-to thick-walled and mostly unbranched, sometimes a few slightly dilated thickwalled terminal cells can be seen; cytoplasmic contents hyaline or with some brownish pigments; walls hyaline to slightly yellowish brown or brownish. Internal hyphae 3 -10 µm wide, thick-walled, clampless or with very scattered and obscure long clamp connections, slightly embedded in a cartilaginous matrix and tightly adhering. Hymenium amphigenous, producing numerous basidia intermixed with somewhat scattered hyphidia. Subhymenium formed by thin-to slightly thickwalled (rarely clearly thick-walled) hyphae, bearing long and ± angulose medallion clamp connections, including at the bases of basidia and hyphidia. Hyphidia distinct but sometimes sparse, branched, 2 -4.5 μm thick in the upper half, terminal branches 1-2.5 μm thick, often with simple septa throughout their length. Young basidia cylindric to narrowly clavate. Mature basidia 2-spored, 35 -80(-85) × (2.5 -)3.5 -5.5(-6.5) μm, sterigmata 9 -23 μm long, often forming an obtuse angle or arising from slightly lateral position from the apex, rarely with an apical papilla. Basidiospores (8.7-)10.1-13.4(-14.8) × (2.9 -) 3.2 -4.3(-4.9), Q = (2.2 -)2.5 -3.9(-4.9), cylindric-allantoid. Basidia and basidiospore cytoplasm hyaline, normally without visible lipid contents, rarely with scattered and minute lipid drops. Microconidia from germinating basidiospores hyaline, narrowly ellipsoid to cylindric, about 3.5 -4.5(-5) × 1-1.5 μm.
Ecology & Distribution -Rather common and widespread, known from the temperate and boreal biomes of North America and Europe but not recorded from Eastern Asia (Shirouzu et al. 2009(Shirouzu et al. , 2016. It is apparently restricted to coniferous wood, mainly Pinus sylvestris in Europe, but also occurs on other genera and species (e.g., Pseudotsuga) in North America. It very frequently occurs on anthropic wood, and is rather widespread on fences, roofs, benches, etc. The ecology of this species was thoroughly detailed by Nannfeldt (1947).
Nomenclature & Taxonomy -No specimens are cited, nor are illustrations or locality data present in the protologue of Sphaeronema rufum (Fries 1818). However, in the later Systema Mycologicum 2, Fries (1823) cited the exiccata 'Scleromyceti Sueciae n. 27'. According to the current Art. F.3.9 of the ICN, and since the name is sanctioned, those specimens are available for lectotype designation by being elements from the context of the sanctioning work. Two duplicates of Scleromyceti Sueciae n. 27 with the original label are kept in UPS. One is mounted on an envelope joined with other samples of the Scleromyceti Sueciae (first 100 numbers, UPS F-117472), and the other is mounted on a free envelope (UPS F-117078). The two speci mens are conspecific, well-developed and typical, despite all fruitbodies being young and sterile. We chose UPS F-117078 as the lectotype in order to facilitate future revisions. Since the selected lectotype is taxonomically unambiguous, the designation of an epitype is unjustifiable (ICN, Art. 9.9), regardless of, e.g., the lack of molecular data from this sample.
We have studied two further specimens collected by E. Fries. One (S F250639) was collected in Sweden without other locotypic indication, and is labelled as 'Original!' by a later hand than that of Fries. Both 'Sphaeronemella rufa (Fr.) Sacc.' and 'Sphaeria rufa Fr.' were written with the same ink and by the same hand, so the label was likely written after Saccardo (1884) made the combination in Sphaeronemella. The other specimen (UPS F-117079) was collected in Femsjö, but the label has no date and Femsjö is not cited as a collecting locality neither in the protologue nor in the sanctioning work. As a result, none of these specimens can be considered as original material.
The initially conical and typically dark reddish brown fruitbodies of Da. rufum have a unique morphology among all Dacrymy cetes, being easily characterised in the field whether in its sterile or fertile state. Indeed, this species is not likely to be confused with any other Dacrymycetes, but some ascomycetes producing dark peritheca or pycnidia with long necks might appear similar to the naked eye. In fact, Fries (1818Fries ( , 1823 placed various ascomycetes in the same genus and described the fruitbodies as perithecia. When fertile, Da. rufum is more difficult to confuse with any other species, although perhaps some confusion with calicioid fungi could arise. A careful examination under the microscope or even with a hand-lens would quickly dispel any doubt. Clamp connections in this species have been largely overlooked, even in the detailed study of Nannfeldt (1947), perhaps because they are normally restricted to the sub- hymenium and the bases of basidia and hyphidia, and because it is difficult to dissect Da. rufum specimens due to their hard texture. However, they were noted and accurately illustrated by Brough & Bandoni (1975), but this was later considered as an exception in some North American specimens (e.g., Laessøe & Petersen 2014). All fertile specimens of Da. rufum studied by us have conspicuous clamp connections at the bases of basidia, regardless of their geographical origin. The branched hyphidia were also not noted by Nannfeldt (1947), who stated: "no sterile elements have been seen in the hymenium". However, although they may be scarce and difficult to find in some samples, and several slides may have to be prepared to locate them, they are present in all fertile specimens studied by us. Brough & Bandoni (1975) also found and illustrated branched hyphidia in the hymenium of this species.
Additional specimens examined (due to the high number of specimens studied, for European material we only cite those specimens used for morphological and molecular analyses, and typification, but c. 50 additional specimens, mainly from the Nordic Countries and preserved in O, S, and UPS, were revised as well). Basidiocarps moderately to clearly gelatinous, but not of soft consistency, xerotolerant, pulvinate to depressed pulvinate when mature, with an almost amphigenous hymenium. Cortical/ marginal hyphae normally unbranched, ± thick-walled, with terminal cells cylindrical to slightly tapering, not conspicuously swollen, with frequent clamp connections but also with some simple septa; anastomoses sparse. Context hyphae slightly to clearly thick-walled, with ± rounded and short clamp connections visible at several septa, but some loop-like clamp connections and simple septa also often present. Subhymenium composed by thin to slightly thick-walled, clamped hyphae. Basidia 1-or 2-spored, bearing a basal clamp connection; when 2-spored apex often U-shaped, but basidia with an apical protuberance or with sterigmata forming an acute angle also occur. Hyphidia sparsely to abundantly branched, rarely anastomosing among them, with a basal clamp connection. Recently discharged and still aseptate basidiospores presumably uninucleate. Basidiospores subglobose to narrowly ovoid or ellipsoid, soon with one transverse septum and later becoming muriform due to the formation of additional longitudinal, oblique and sometimes transverse septa, thin-walled, hyaline or with a very subtle yellowish tinge. Spore print whitish. Microconidia ellipsoid. Cell cytoplasm without conspicuous carotenoids under the light microscope, at most, some basidia and basidiospores with a very subtle creamish tinge. Parietal diffuse yellowish brown pigments sometimes seen at some cortical/marginal hyphae. Basidia and basidiospores with conspicuous lipid bodies, becoming ≥ 2 μm lipid drops at least in some old cells.
Ecology & Distribution -Unilacryma bispora is a very common species that has only been observed fruiting on decorticated, old and frequently exposed wood of Pinus sylvestris, especially on branches still attached to living trees. It occurs in both temperate and boreal forests, but it seems to be particularly common in the boreal biome, reaching the northern boreal zone. It is currently known from Finland, Norway, Russia, and Sweden. In a study on the fungal communities inhabiting Picea abies stumps (Kubart et al. 2016), environmental sequences probably related to this species were found in three Swedish nature reserves (Altarliden, Västerbotten; Fänstjärnskogen, Värmland; Ugglehöjden, Närke). The representative ITS sequence deposited in GenBank (KM494136) is identical to some of the ITS sequences of U. bispora generated in this study. Based on the results from the species delimitation analyses and considering the observed distribution in the field, we consider it very likely that those environmental sequences indicate that the species has a wider ecology than we have observed.
Nomenclature & Taxonomy -Unilacryma bispora is readily distinguished by the combination of small, pulvinate and often amber-coloured basidiocarps, 2-spored basidia, ovoid to ellipsoid basidiospores that become muriform when mature, branched hyphidia, and absence of obvious carotenoids under the light microscope. It is similar to U. unispora except for the slightly broader, bisterigmate basidia, the slightly narrower and more elongate basidiospores, and the more northern distribution. The basidiospore shape is nevertheless somewhat plastic, and some specimens have predominantly ovoid basidiospores (e.g., UPS F-941254), while others show mostly ellipsoid ones (e.g., UPS F-941255). There may be little overlap among the length/width ratio of the basidiospores in certain specimens, but since most of them show intermediate shapes, and we have not detected any other differences, this should be regarded as part of the intraspecific variation. The holotype and the isotypes were collected on different, unconnected pine branches and are different individuals of the same gathering (ICN, Art. 8.2).
Following the literature and published identification keys (e.g., Torkelsen 1997), it might be confused with Dacrymyces ovi sporus, a species with somewhat similar basidia and basidiospores, but with clearly larger, yellow-orange basidiocarps (Fig.  1f), and with very obvious carotenoid contents in the cytoplasm of all cells, including the basidiospores (Fig. 5s). The comments by Oberwinkler (2014) on 'Dacrymyces ovisporus', indicating a macromorphology identical to that of Unilacryma unispora, and a very close micromorphology, might mean that the specimens studied by Oberwinkler actually belong to U. bispora. In its most commonly accepted sense, Dacrymyces ovisporus is additionally distinguished by often staining the substrate yellow (Torkelsen 1997), by hyphidia that are simple or rarely with 1-2 short branches (Fig. 2l), and by both the basidiospore shape and septation, which is more similar to that of U. unispora (thus, more rounded and becoming muriform earlier), this last character being well represented in the original plate of the protologue (Brefeld 1888). Further details on D. ovisporus will be provided in a forthcoming study. The two species are distantly related according to our phylogeny, belonging to two different families (Fig. 8 -9). Basidiocarps ± gelatinous, amber-coloured to brownish, xerotolerant, isolated, gregarious or crowded, initially erumpent, pulvinate, tardily becoming depressed-pulvinate, without a conspicuous margin, 0.3 -1.0(-1.2) mm diam. Cortical/mar ginal hyphae slightly to clearly thick-walled, with cylindrical or attenuate terminal cells, 4 -7 µm wide, unbranched or rarely sparingly branched, with sparse anastomoses and both clamped and clampless septa. Internal hyphae 3 -7 µm wide, slightly to clearly thick-walled, embedded in a gelatinous matrix, with clamp connections but also frequently with some simple septa. Hymenium almost amphigenous, progressively forming numerous basidia intermixed with some hyphidia. Subhyme nium with thin to slightly thick-walled, clamped hyphae. Hyphidia hyaline, sparsely to heavily branched, (2.5 -)3 -4.5 μm thick in the upper half, terminal branches (1-)1.5-3 μm thick, often with clamped septa throughout their length. Young basidia initials cylindric to narrowly claviform, 5 -7.5 μm thick when close to maturity. Mature basidia developing a single thick sterigma (very rarely two irregularly formed sterigma-like structures may be seen), about 18 -32 μm long, bearing a single terminal and asymmetrically attached basidiospore. Total basidium length reaching 70 -130 μm. Basidiospores (11.9 -)14 -17.4(-19.5) × (10.2-)11.2-13.5(-15.7) μm, Q = (1.1-)1.2-1.4(-1.6), ovoid to broadly ellipsoid, rarely globose or ellipsoid, with a conspicuous 1.0 -1.5 μm long hilar appendix, becoming transversely 1-septate and developing additional oblique, longitudinal and transverse septa, then becoming muriform, with a single hardly visible nucleus when still aseptate. Basidium and basidiospore cytoplasm hyaline or with a very pale yellowish tinge, with numerous small lipid bodies that slowly coalesce in dead state to form conspicuous lipid drops. Basidiospores germinating by the production of hyphae or microconidia. Microconidia rarely Fig. 15 Unilacryma unispora s.lat. a -d. Fruitbodies in different degrees of development, all in fresh condition except a; e -m. Micromorphology; e. marginal hyphae; f-g. branched hyphidia; h. terminal cells of marginal hyphae, note the presence of both clampless (black arrows) and clamped (white arrow) septa; i. young basidia; j -k. mature basidia; l. basidiospores; m. two young basidia and one mature basidium bearing a single sterigma with a terminal basidiospore. a. From the holotype; b, e -g, j -k. from UPS F-941279; c. from UPS F-941286; d. from UPS F-941284; h -i, m. from UPS F-941283; l. from UPS F-941280. -Scale bars: a -d = 1 mm, e -m = 10 µm. observed, narrowly ellipsoid to cylindric, about 3 -4 × 1.5 μm, hyaline.
Ecology & Distribution -Very common, frequently growing on decorticated, old and rather exposed wood of Juniperus spp. and Pinus spp., and particularly common on dead branches still attached to trees. In areas where the distributions of both Unilacryma species overlap, they may grow together on the same branch and even intermixed. In Europe, U. unispora has been previously reported only from The Netherlands by Dam & Dam (2014). In a wide sense, it is broadly distributed in the Holarctic temperate biome (North America, East Asia, and Europe), and extends into the southern boreal zone, but it has not been found in the middle and northern boreal forests of Sweden and Finland where U. bispora is common.
Nomenclature & Taxonomy -When describing Platygloea unispora, Olive (1944) failed to include the required Latin diagnosis, resulting in an invalid name (ICN, Art. 39.1). The taxon was validated by Olive (1958), who provided the required diagno sis as well as a holotype. Olive (1944) is often incorrectly cited as the place of publication of the basionym, which renders the intended new combinations in Achroomyces and Dacrymy ces invalid, none of the requirements of Art. 41.8 being fulfilled.
In a broad taxonomic sense, U. unispora is well distinguished by 1-spored basidia and subglobose to ovoid basidiospores that become transversely septate and finally muriform with age. The macromorphology is nearly identical to that of U. bispora, so a micromorphological study of the samples is always needed to identify them. The European specimens agree morphologically very well with the description of Olive (1944) and additional descriptions and iconography of North American specimens (e.g., Wells 1994), as well as with the North American specimens we studied, including the holotype. Unfortunately, we were not able to obtain DNA sequence data of the studied North American specimens, not even after using specific primers to amplify short fragments of the ITS. The description of the Japanese specimens by Shirouzu et al. (2009) also agrees relatively well, although they reported basidia and basidiospores to be on average slightly smaller, marginal hyphae as thin-walled, and branched hyphidia as absent. In addition, the samples from North Europe and East Asia form two different clades (Fig. 10). Nevertheless, with the current morphological and genetic data, it is unclear whether these clades represent two distinct species. Therefore, we assume here a conservative approach and treat the North American, East Asian, and European specimens as different populations of a single species.