Large-scale fungal strain sequencing unravels the molecular diversity in mating loci maintained by long-term balancing selection.

Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of strains. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci define mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty strains of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. The sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. 17,550 mating types were predicted based on the allelic classes. In vitro crosses allowed us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. We conclude that the genetic diversity of mating loci in Trichaptum is due to long-term balancing selection, with limited recombination and duplication activity. The large number of sequenced strains highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Global field collection data confirm an affinity of brown rot fungi for coniferous habitats and substrates.

Unlike 'white rot' (WR) wood-decomposing fungi that remove lignin to access cellulosic sugars, 'brown rot' (BR) fungi selectively extract sugars and leave lignin behind. The relative frequency and distribution of these fungal types (decay modes) have not been thoroughly assessed at a global scale; thus, the fate of one-third of Earth's aboveground carbon, wood lignin, remains unclear. Using c. 1.5 million fungal sporocarp and c. 30 million tree records from publicly accessible databases, we mapped and compared decay mode and tree type (conifer vs angiosperm) distributions. Additionally, we mined fungal record metadata to assess substrate specificity per decay mode. The global average for BR fungi proportion (BR/(BR + WR records)) was 13% and geographic variation was positively correlated (R2 = 0.45) with conifer trees proportion (conifer/(conifer + angiosperm records)). Most BR species (61%) were conifer, rather than angiosperm (22%), specialists. The reverse was true for WR (conifer: 19%; angiosperm: 62%). Global BR proportion patterns were predicted with greater accuracy using the relative distributions of individual tree species (R2 = 0.82), rather than tree type. Fungal decay mode distributions can be explained by tree type and, more importantly, tree species distributions, which our data suggest is due to strong substrate specificities.

Introgression between highly divergent fungal sister species.

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

The Indoor Mycobiomes of Daycare Centers Are Affected by Occupancy and Climate.

Many children spend considerable time in daycare centers and may be influenced by the indoor microorganisms there, including fungi. In this study, we investigate the indoor mycobiomes of 125 daycare centers distributed along strong environmental gradients throughout Norway. Dust samples were collected from doorframes outside and inside buildings using a community science sampling approach. Fungal communities in the dust samples were analyzed using DNA metabarcoding of the internal transcribed spacer 2 (ITS2) region. We observed a marked difference between the outdoor and indoor mycobiomes. The indoor mycobiomes included considerably more yeasts and molds than the outdoor samples, with Saccharomyces, Mucor, Malassezia, and Penicillium being among the most dominant fungal genera. Changes in the indoor fungal richness and composition correlated with numerous variables related to both outdoor and indoor conditions; there was a clear geographic structure in the indoor mycobiome composition that mirrored the outdoor climate, ranging from humid areas in western Norway to drier and colder areas in eastern Norway. Moreover, the number of children in the daycare centers, as well as various building features, influenced the indoor mycobiome composition. We conclude that the indoor mycobiomes in Norwegian daycare centers are structured by multiple factors and are dominated by yeasts and molds. This study exemplifies how community science sampling enables DNA-based analyses of a high number of samples covering wide geographic areas. IMPORTANCE With an alarming increase in chronic diseases like childhood asthma and allergies, there is an increased focus on the exposure of young children to indoor biological and chemical air pollutants. Our study of 125 daycares throughout Norway demonstrates that the indoor mycobiome not only reflects cooccurring outdoor fungi but also includes a high abundance of yeast and mold fungi with an affinity for indoor environments. A multitude of factors influence the indoor mycobiomes in daycares, including the building type, inhabitants, as well as the outdoor environment. Many of the detected yeasts and molds are likely associated with the human body, where some have been coupled with allergies and respiratory problems. Our results call for further studies investigating the potential impact of the identified daycare-associated mycobiomes on children's health.

Population genomics of a forest fungus reveals high gene flow and climate adaptation signatures.

Genome sequencing of spatially distributed individuals sheds light on how evolution structures genetic variation. Populations of Phellopilus nigrolimitatus, a red-listed wood-inhabiting fungus associated with old-growth coniferous forests, have decreased in size over the last century due to a loss of suitable habitats. We assessed the population genetic structure and investigated local adaptation in P. nigrolimitatus, by establishing a reference genome and genotyping 327 individuals sampled from 24 locations in Northern Europe by RAD sequencing. We revealed a shallow population genetic structure, indicating large historical population sizes and high levels of gene flow. Despite this weak substructuring, two genetic groups were recognized; a western group distributed mostly in Norway and an eastern group covering most of Finland, Poland and Russia. This substructuring may reflect coimmigration with the main host, Norway spruce (Picea abies), into Northern Europe after the last ice age. We found evidence of low levels of genetic diversity in southwestern Finland, which has a long history of intensive forestry and urbanization. Numerous loci were significantly associated with one or more environmental factors, indicating adaptation to specific environments. These loci clustered into two groups with different associations with temperature and precipitation. Overall, our findings indicate that the current population genetic structure of P. nigrolimitatus results from a combination of gene flow, genetic drift and selection. The acquisition of similar knowledge especially over broad geographic scales, linking signatures of adaptive genetic variation to evolutionary processes and environmental variation, for other fungal species will undoubtedly be useful for assessment of the combined effects of habitat fragmentation and climate change on fungi strongly bound to old-growth forests.

Analysing indoor mycobiomes through a large-scale citizen science study in Norway.

In the built environment, fungi can cause important deterioration of building materials and have adverse health effects on occupants. Increased knowledge about indoor mycobiomes from different regions of the world, and their main environmental determinants, will enable improved indoor air quality management and identification of health risks. This is the first citizen science study of indoor mycobiomes at a large geographical scale in Europe, including 271 houses from Norway and 807 dust samples from three house compartments: outside of the building, living room and bathroom. The fungal community composition determined by DNA metabarcoding was clearly different between indoor and outdoor samples, but there were no significant differences between the two indoor compartments. The 32 selected variables, related to the outdoor environment, building features and occupant characteristics, accounted for 15% of the overall variation in community composition, with the house compartment as the key factor (7.6%). Next, climate was the main driver of the dust mycobiomes (4.2%), while building and occupant variables had significant but minor influences (1.4% and 1.1%, respectively). The house-dust mycobiomes were dominated by ascomycetes (⁓70%) with Capnodiales and Eurotiales as the most abundant orders. Compared to the outdoor samples, the indoor mycobiomes showed higher species richness, which is probably due to the mixture of fungi from outdoor and indoor sources. The main indoor indicator fungi belonged to two ecological groups with allergenic potential: xerophilic moulds and skin-associated yeasts. Our results suggest that citizen science is a successful approach for unravelling the built microbiome at large geographical scales.

Contrasting demographic histories revealed in two invasive populations of the dry rot fungus Serpula lacrymans.

Globalization and international trade have impacted organisms around the world leading to a considerable number of species establishing in new geographic areas. Many organisms have taken advantage of human-made environments, including buildings. One such species is the dry rot fungus Serpula lacrymans, which is the most aggressive wood-decay fungus in indoor environments in temperate regions. Using population genomic analyses of 36 full genome sequenced isolates, we demonstrated that European and Japanese isolates are highly divergent and the populations split 3000-19,000 generations ago, probably predating human influence. Approximately 250 generations ago, the European population went through a tight bottleneck, probably corresponding to the fungus colonization of the built environment in Europe. The demographic history of these populations, probably lead to low adaptive potential. Only two loci under selection were identified using a Fst outlier approach, and selective sweep analyses identified three loci with extended haplotype homozygosity. The selective sweep analyses found signals in genes possibly related to decay of various substrates in Japan and in genes involved DNA replication and protein modification in Europe. Our results suggest that the dry rot fungus independently established in indoor environments in Europe and Japan and that invasive species can potentially establish large populations in new habitats based on a few colonizing individuals.

In vitro evidence of root colonization suggests ecological versatility in the genus Mycena.

The root-associated habit has evolved on numerous occasions in different fungal lineages, suggesting a strong evolutionary pressure for saprotrophic fungi to switch to symbiotic associations with plants. Species within the ubiquitous, saprotrophic genus Mycena are frequently major components in molecular studies of root-associated fungal communities, suggesting that an evaluation of their trophic status is warranted. Here, we report on interactions between a range of Mycena species and the plant Betula pendula. In all, 17 Mycena species were inoculated onto B. pendula seedlings. Physical interactions between hyphae and fine roots were examined using differential staining and fluorescence microscopy. Physiological interactions were investigated using 14 C and 32 P to show potential transfer between symbionts. All Mycena species associated closely with fine roots, showing hyphal penetration into the roots, which in some cases were intracellular. Seven species formed mantle-like structures around root tips, but none formed a Hartig net. Mycena pura and Mycena galopus both enhanced seedling growth, with M. pura showing significant transfer of 32 P to the seedlings. Our results support the view that several Mycena species can associate closely with plant roots and some may potentially occupy a transitional state between saprotrophy and biotrophy.

Rapid Divergence of Genome Architectures Following the Origin of an Ectomycorrhizal Symbiosis in the Genus Amanita.

Fungi are evolutionary shape shifters and adapt quickly to new environments. Ectomycorrhizal (EM) symbioses are mutualistic associations between fungi and plants and have evolved repeatedly and independently across the fungal tree of life, suggesting lineages frequently reconfigure genome content to take advantage of open ecological niches. To date analyses of genomic mechanisms facilitating EM symbioses have involved comparisons of distantly related species, but here, we use the genomes of three EM and two asymbiotic (AS) fungi from the genus Amanita as well as an AS outgroup to study genome evolution following a single origin of symbiosis. Our aim was to identify the defining features of EM genomes, but our analyses suggest no clear differentiation of genome size, gene repertoire size, or transposable element content between EM and AS species. Phylogenetic inference of gene gains and losses suggests the transition to symbiosis was dominated by the loss of plant cell wall decomposition genes, a confirmation of previous findings. However, the same dynamic defines the AS species A. inopinata, suggesting loss is not strictly associated with origin of symbiosis. Gene expansions in the common ancestor of EM Amanita were modest, but lineage specific and large gene family expansions are found in two of the three EM extant species. Even closely related EM genomes appear to share few common features. The genetic toolkit required for symbiosis appears already encoded in the genomes of saprotrophic species, and this dynamic may explain the pervasive, recurrent evolution of ectomycorrhizal associations.

Wood Modification by Furfuryl Alcohol Resulted in a Delayed Decomposition Response in Rhodonia (Postia) placenta.

The aim of this study was to investigate differential expression profiles of the brown rot fungus Rhodonia placenta (previously Postia placenta) harvested at several time points when grown on radiata pine (Pinus radiata) and radiata pine with three different levels of modification by furfuryl alcohol, an environmentally benign commercial wood protection system. The entire gene expression pattern of a decay fungus was followed in untreated and modified wood from initial to advanced stages of decay. The results support the current model of a two-step decay mechanism, with the expression of genes related to initial oxidative depolymerization, followed by an accumulation of transcripts of genes related to the hydrolysis of cell wall polysaccharides. When the wood decay process is finished, the fungus goes into starvation mode after five weeks when grown on unmodified radiata pine wood. The pattern of repression of oxidative processes and oxalic acid synthesis found in radiata pine at later stages of decay is not mirrored for the high-furfurylation treatment. The high treatment level provided a more unpredictable expression pattern throughout the incubation period. Furfurylation does not seem to directly influence the expression of core plant cell wall-hydrolyzing enzymes, as a delayed and prolonged, but similar, pattern was observed in the radiata pine and the modified experiments. This indicates that the fungus starts a common decay process in the modified wood but proceeds at a slower pace as access to the plant cell wall polysaccharides is restricted. This is further supported by the downregulation of hydrolytic enzymes for the high treatment level at the last harvest point (mass loss, 14%). Moreover, the mass loss does not increase during the last weeks. Collectively, this indicates a potential threshold for lower mass loss for the high-furfurylation treatment.IMPORTANCE Fungi are important decomposers of woody biomass in natural habitats. Investigation of the mechanisms employed by decay fungi in their attempt to degrade wood is important for both the basic scientific understanding of ecology and carbon cycling in nature and for applied uses of woody materials. For wooden building materials, long service life and carbon storage are essential, but decay fungi are responsible for massive losses of wood in service. Thus, the optimization of durable wood products for the future is of major importance. In this study, we have investigated the fungal genetic response to furfurylated wood, a commercial environmentally benign wood modification approach that improves the service life of wood in outdoor applications. Our results show that there is a delayed wood decay by the fungus as a response to furfurylated wood, and new knowledge about the mechanisms behind the delay is provided.

A single ectomycorrhizal plant root system includes a diverse and spatially structured fungal community.

Although only a relatively small proportion of plant species form ectomycorrhizae with fungi, it is crucial for growth and survival for a number of widespread woody plant species. Few studies have attempted to investigate the fine scale spatial structure of entire root systems of adult ectomycorrhizal (EcM) plants. Here, we use the herbaceous perennial Bistorta vivipara to map the entire root system of an adult EcM plant and investigate the spatial structure of its root-associated fungi. All EcM root tips were sampled, mapped and identified using a direct PCR approach and Sanger sequencing of the internal transcribed spacer region. A total of 32.1% of all sampled root tips (739 of 2302) were successfully sequenced and clustered into 41 operational taxonomic units (OTUs). We observed a clear spatial structuring of the root-associated fungi within the root system. Clusters of individual OTUs were observed in the younger parts of the root system, consistent with observations of priority effects in previous studies, but were absent from the older parts of the root system. This may suggest a succession and fragmentation of the root-associated fungi even at a very fine scale, where competition likely comes into play at different successional stages within the root system.

Population genomic analyses reveal possible drivers of population divergence.

Recent advances in sequencing technology and efficiency enable new and improved methods to investigate how populations diverge and species evolve. Fungi have relatively small and simple genomes and can often be cultured in the laboratory. Fungal populations can thus be sequenced for a relatively low cost, which makes them ideal for population genomic analyses. In several recent population genomic studies, wild populations of fungal model organisms and human pathogens have been analysed, for example Neurospora crassa (Ellison et al. ), Saccharomyces uvarum (Almeida et al. ), Coccidioides spp. (Neafsey et al. ) and Cryptococcus gatti (Engelthaler et al. ). In this issue of Molecular Ecology, Branco et al. () apply population genomic tools to understand population divergence and adaptation in a symbiotic (mycorrhizal) fungus. This study exemplifies the possibilities of diving deeper into the genomic features involved in population divergence and speciation, also for nonmodel organisms, and how molecular and analytical tools will improve our understanding of the patterns and mechanisms that underlie adaptation to habitats, population divergence and dispersal limitation of fungi.

Population structure of Serpula lacrymans in Europe with an outlook to the French population.

In this study the genetic variation and population structure in a French population of the dry rot fungus S. lacrymans was investigated using 14 microsatellites markers and compared to the rest of Europe. In that comparison the French population possessed the same allelic diversity as rest of Europe. A weak geographic structuring of the genetic variation was observed across Europe, where the French isolates to some extent separated from the rest of Europe, indicating that weak barriers to gene flow exists. Eighty percent of the isolates had unique multilocus microsatellite genotypes, which corresponds to high recombination and dispersal by sexual spores. Deviations from Hardy-Weinberg expectations were observed in multiple loci. In most loci there was an excess of heterozygotes, which could be due to either non-random mating, presence of more than two nuclei in the secondary mycelia or another unrecognized process. A total of six vegetative compatibility (VC) groups were present in Europe, out of which four were sampled in France. One VC group was over-represented in France while two others were underrepresented, as compared to the rest of Europe.

The Ribosomal Operon Database: A Full-Length rDNA Operon Database Derived From Genome Assemblies.

Current rDNA reference sequence databases are tailored towards shorter DNA markers, such as parts of the 16/18S marker or the internally transcribed spacer (ITS) region. However, due to advances in long-read DNA sequencing technologies, longer stretches of the rDNA operon are increasingly used in environmental sequencing studies to increase the phylogenetic resolution. There is, therefore, a growing need for longer rDNA reference sequences. Here, we present the ribosomal operon database (ROD), which includes eukaryotic full-length rDNA operons fished from publicly available genome assemblies. Full-length operons were detected in 34.1% of the 34,701 examined eukaryotic genome assemblies from NCBI. In most cases (53.1%), more than one operon variant was detected, which can be due to intragenomic operon copy variability, allelic variation in non-haploid genomes, or technical errors from the sequencing and assembly process. The highest copy number found was 5947 in Zea mays. In total, 453,697 unique operons were detected, with 69,480 operon variant clusters remaining after intragenomic clustering at 99% sequence identity. The operon length varied extensively across eukaryotes, ranging from 4136 to 16,463 bp, which will lead to considerable polymerase chain reaction (PCR) bias during amplification of the entire operon. Clustering the full-length operons revealed that the different parts (i.e., 18S, 28S, and the hypervariable regions V4 and V9 of 18S) provide divergent taxonomic resolution, with 18S, the V4 and V9 regions being the most conserved. The ROD will be updated regularly to provide an increasing number of full-length rDNA operons to the scientific community.

Reticulate evolution and rapid development of reproductive barriers upon secondary contact in a forest fungus.

Reproductive barriers between sister species of the mushroom-forming fungi tend to be stronger in sympatry, leading to speculation on whether they are being reinforced by selection against hybrids. We have used population genomic analyses together with in vitro crosses of a global sample of the wood decay fungus Trichaptum abietinum to investigate reproductive barriers within this species complex and the processes that have shaped them. Our phylogeographic analyses show that T. abietinum is delimited into six major genetic groups: one in Asia, two in Europe, and three in North America. The groups present in Europe are interfertile and admixed, whereas our crosses show that the North American groups are reproductively isolated. In Asia, a more complex pattern appears, with partial intersterility between subgroups that likely originated independently and more recently than the reproductive barriers in North America. We found pre-mating barriers in T. abietinum to be moderately correlated with genomic divergence, whereas mean growth reduction of the mated hybrids showed a strong correlation with increasing genomic divergence. Genome-wide association analyses identified candidate genes with programmed cell death annotations, which are known to be involved in intersterility in distantly related fungi, although their link here remains unproven. Our demographic modeling and phylogenetic network analyses fit a scenario where reproductive barriers in Trichaptum abietinum could have been reinforced upon secondary contact between groups that diverged in allopatry during the Pleistocene glacial cycles. Our combination of experimental and genomic approaches demonstrates how T. abietinum is a tractable system for studying speciation mechanisms.

Extreme overall mushroom genome expansion in Mycena s.s. irrespective of plant hosts or substrate specializations.

Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

Cytochrome P450 Complement May Contribute to Niche Adaptation in Serpula Wood-Decay Fungi.

Serpula wood-decay fungi occupy a diverse range of natural and man-made ecological niches. Serpula himantioides is a forest-floor generalist with global coverage and strong antagonistic ability, while closely related species Serpula lacrymans contains specialist sister strains with widely differing ecologies. Serpula lacrymans var. shastensis is a forest-floor specialist in terms of resource preference and geographic coverage, while Serpula lacrymans var. lacrymans has successfully invaded the built environment and occupies a building-timber niche. To increase understanding of the cellular machinery required for niche adaptation, a detailed study of the P450 complement of these three strains was undertaken. Cytochrome P450 monooxygenases are present in all fungi and typically seen in high numbers in wood decay species, with putative roles in breakdown of plant extractives and lignocellulose metabolism. Investigating the genomes of these related yet ecologically diverse fungi revealed a high level of concordance in P450 complement, but with key differences in P450 family representation and expression during growth on wood, suggesting P450 proteins may play a role in niche adaptation. Gene expansion of certain key P450 families was noted, further supporting an important role for these proteins during wood decay. The generalist species S. himantioides was found to have the most P450 genes with the greatest family diversity and the highest number of P450 protein families expressed during wood decay.

Evolutionary history of Serpulaceae (Basidiomycota): molecular phylogeny, historical biogeography and evidence for a single transition of nutritional mode.

BACKGROUND: The fungal genus Serpula (Serpulaceae, Boletales) comprises several saprotrophic (brown rot) taxa, including the aggressive house-infecting dry rot fungus Serpula lacrymans. Recent phylogenetic analyses have indicated that the ectomycorrhiza forming genera Austropaxillus and Gymnopaxillus cluster within Serpula. In this study we use DNA sequence data to investigate phylogenetic relationships, historical biogeography of, and nutritional mode transitions in Serpulaceae. RESULTS: Our results corroborate that the two ectomycorrhiza-forming genera, Austropaxillus and Gymnopaxillus, form a monophyletic group nested within the saprotrophic genus Serpula, and that the Serpula species S. lacrymans and S. himantioides constitute the sister group to the Austropaxillus-Gymnopaxillus clade. We found that both vicariance (Beringian) and long distance dispersal events are needed to explain the phylogeny and current distributions of taxa within Serpulaceae. Our results also show that the transition from brown rot to mycorrhiza has happened only once in a monophyletic Serpulaceae, probably between 50 and 22 million years before present. CONCLUSIONS: This study supports the growing understanding that the same geographical barriers that limit plant- and animal dispersal also limit the spread of fungi, as a combination of vicariance and long distance dispersal events are needed to explain the present patterns of distribution in Serpulaceae. Our results verify the transition from brown rot to ECM within Serpulaceae between 50 and 22 MyBP.

Contrasting genetic structuring in the closely related basidiomycetes Trichaptum abietinum and Trichaptum fuscoviolaceum (Hymenochaetales).

Trichaptum abietinum and Trichaptum fuscoviolaceum (Hymenochaetales, Basidiomycota) are closely related saprotrophic fungi, widely distributed on coniferous wood in temperate regions worldwide. Three intersterility groups have previously been detected in T. abietinum, while no prezygotic barriers have been proven within T. fuscoviolaceum. The aim of this study was to reveal the phylogeography and genetic relationship between these two closely related species and to explore whether the previously observed intersterility groups in T. abietinum are reflected in the genetic data. We assembled worldwide fruit body collections of both species (N = 314) and generated DNA sequences from three nuclear (ITS2, LSU, IGS) and one mitochondrial rDNA region (mtLSU). The two species are genetically well separated in all analyses. In correspondence with observations from earlier mating studies, our results revealed that T. fuscoviolaceum is genetically more uniform than T. abietinum. Multiple genetic sub-groups exist in T. abietinum that may correspond to the previously observed intersterility groups. However, there is low consistency across the investigated loci in delimiting the different sub-groups, except for a consistent North American group. As for many other widespread fungi, a complex phylogeographic pattern is found in T. abietinum which may have been formed by geographic, as well as multiple genetic intersterility barriers.

Spatiotemporal variation of the indoor mycobiome in daycare centers.

BACKGROUND: Children spend considerable time in daycare centers in parts of the world and are exposed to the indoor micro- and mycobiomes of these facilities. The level of exposure to microorganisms varies within and between buildings, depending on occupancy, climate, and season. In order to evaluate indoor air quality, and the effect of usage and seasonality, we investigated the spatiotemporal variation in the indoor mycobiomes of two daycare centers. We collected dust samples from different rooms throughout a year and analyzed their mycobiomes using DNA metabarcoding. RESULTS: The fungal community composition in rooms with limited occupancy (auxiliary rooms) was similar to the outdoor samples, and clearly different from the rooms with higher occupancy (main rooms). The main rooms had higher abundance of Ascomycota, while the auxiliary rooms contained comparably more Basidiomycota. We observed a strong seasonal pattern in the mycobiome composition, mainly structured by the outdoor climate. Most markedly, basidiomycetes of the orders Agaricales and Polyporales, mainly reflecting typical outdoor fungi, were more abundant during summer and fall. In contrast, ascomycetes of the orders Saccharomycetales and Capnodiales were dominant during winter and spring. CONCLUSIONS: Our findings provide clear evidences that the indoor mycobiomes in daycare centers are structured by occupancy as well as outdoor seasonality. We conclude that the temporal variability should be accounted for in indoor mycobiome studies and in the evaluation of indoor air quality of buildings. Video abstract.