FEMS Yeast Research, 16, 2016, fov103 doi: 10.1093/femsyr/fov103 Advance Access Publication Date: 13 November 2015 Research Article RESEARCH ARTICLE Local climatic conditions constrain soil yeast diversity patterns in Mediterranean forests, woodlands and scrub biome Andrey M. Yurkov1,∗ , Oliver Röhl2,†,‡ , Ana Pontes2,‡ , Cláudia Carvalho2 , Cristina Maldonado2 and José Paulo Sampaio2 1 Leibniz Institute DSMZ-German Collection of Microorganisms and Cell Cultures, Inhoffenstraße 7B, 38124 Brunswick, Germany and 2 UCIBIO-REQUIMTE, Departamento de Ciências da Vida, Faculdade de Ciências e Tecnologia, Universidade Nova de Lisboa, 2829-516 Caparica, Portugal ∗ Corresponding author: Leibniz Institute DSMZ-German Collection of Microorganisms and Cell Cultures, Inhoffenstraße 7B, 38124 Brunswick, Germany. Tel: +493512616239; Fax: +493512616418; E-mail: [email protected] † Present address: Geobotany, Faculty for Biology and Biotechnology, Ruhr-Universität Bochum, 44780 Bochum, Germany. ‡ These authors contributed equally to this work. One sentence summary: The authors provide the first comprehensive inventory of soil yeasts in a Mediterranean ecosystem and evaluate the effect of microclimate factors on yeast community in a transect from humid broadleaf forests to dry maquis. Editor: Cletus Kurtzman ABSTRACT Soil yeasts represent a poorly known fraction of the soil microbiome due to limited ecological surveys. Here, we provide the first comprehensive inventory of cultivable soil yeasts in a Mediterranean ecosystem, which is the leading biodiversity hotspot for vascular plants and vertebrates in Europe. We isolated and identified soil yeasts from forested sites of Serra da Arrábida Natural Park (Portugal), representing the Mediterranean forests, woodlands and scrub biome. Both cultivation experiments and the subsequent species richness estimations suggest the highest species richness values reported to date, resulting in a total of 57 and 80 yeast taxa, respectively. These values far exceed those reported for other forest soils in Europe. Furthermore, we assessed the response of yeast diversity to microclimatic environmental factors in biotopes composed of the same plant species but showing a gradual change from humid broadleaf forests to dry maquis. We observed that forest properties constrained by precipitation level had strong impact on yeast diversity and on community structure and lower precipitation resulted in an increased number of rare species and decreased evenness values. In conclusion, the structure of soil yeast communities mirrors the environmental factors that affect aboveground phytocenoses, aboveground biomass and plant projective cover. Keywords: yeasts; soils; Mediterranean basin; biodiversity hotspot; Basidiomycetes; microclimate INTRODUCTION Yeasts inhabit soils worldwide (Botha 2011; Yurkov, Kemler and Begerow 2012). Even though first isolations of yeasts from soils have been made more than a century ago (reviewed by Guilliermond 1920), our knowledge of the determinants of soil yeast diversity and density is restricted to a limited number of studies focused on biotopes like arable fields (e.g. Starkey and Henrici 1927; Bouthilet 1951; Sláviková and Vadkertiová 2003; Lynch and Thorn 2006) and a few major forest types (e.g. di Received: 27 September 2015; Accepted: 10 November 2015 C FEMS 2015. All rights reserved. For permissions, please e-mail: [email protected] 1 2 FEMS Yeast Research, 2016, Vol. 16, No. 1 Menna 1965; Sláviková and Vadkertiová 2000; Maksimova and Chernov 2004; Mestre et al. 2011; Yurkov, Kemler and Begerow 2012). As a consequence, many biotopes and soil types still lack an inventory of the correspondent yeast community, a situation that is further complicated by the difficulty in tracing species in older publications due to changing species concepts and different nomenclatures. In the last few decades, the use of ribosomal DNA sequences greatly improved yeast identification and increased the depth of ecological studies. Moreover, species recognitions and comparisons across different studies become possible through sequence data deposited in public repositories such as NCBI GenBank and MycoBank. In addition to their contribution to biodiversity, soil yeasts are gaining increasing importance in biotechnology as enzyme producers and as sources of oils (e.g. Buzzini and Margesin 2014; Schulze et al. 2014; Tanimura et al. 2014). However, undersampling strongly limits the discovery of new soil yeast species (e.g. Yurkov, Kemler and Begerow 2011) and consequently yeast diversity in most soil types remains largely unknown. Forest biotopes harbour diverse above- and belowground cultivable yeast populations predominantly consisting of basidiomycetous yeasts (Sláviková and Vadkertiová 2000; Maksimova and Chernov 2004; Fonseca and Inácio 2006; Yurkov, Kemler and Begerow 2012; Yurkov et al. 2015a). According to current knowledge, the type of forest determines the size, diversity and structure of the yeast community (reviewed in Yurkov et al. 2015a). However, different types of forests occur in distinct climatic zones and this makes the assessment of yeast biodiversity biased by the effects of local climatic parameters, both aboveground and in soils or soil litter. Chernov (2005) and Vishniac (2006) investigated the structure of the community of soil yeasts along vast latitudinal gradients and reported that species composition changed significantly with rainfall values. However, it is important to document that the most pronounced differences were observed between sampling sites in regions that were distantly separated, i.e. in a different natural zone and under a different type of vegetation. To our knowledge, there was no study designed to specifically evaluate the effects of either vegetation properties or climatic conditions on soil yeasts. However, there is some indirect evidence that the two factors are indeed important constraints of the yeast community. Yeast communities isolated in a spruce forest (zonal vegetation type) for example differed from those observed in birch and alder forests (intrazonal vegetation type) (Maksimova and Chernov 2004). Similarly, the yeast community of Sphagnum mosses present in forest ground cover showed a gradual change along the temperature-moisture ecocline between boreal spruce forest (zonal) and swamp (intrazonal) biotopes (Kachalkin and Yurkov 2012). In agreement with those observations, soil yeast communities isolated from clear- cut planted spruce forests differed substantially from the yeast populations found in near-natural unmanaged beech forests (Yurkov, Kemler and Begerow 2012). It was also observed that total yeast numbers and diversity did not depend on basic soil properties (Birkhofer et al. 2012), thus suggesting the importance of aboveground biotope properties for the yeast community. Here we report on a study designed to assess the effects of local microclimatic conditions on basic parameters of soil yeast communities. The goal of this investigation was to study abundance, diversity (observed and estimated), composition and structure (evenness, number of singletons, proportion of rare species) of the yeast communities in soils characterized by a highly similar vegetation cover but strong differences in precipitation levels. The study was performed in Serra da Arrábida Natural Park (Portugal) on sites previously studied for phylloplane yeasts (Inácio et al. 2002). This study also provides the first comprehensive inventory of cultivable soil yeasts in a Mediterranean ecosystem and includes in the same sampling strategy the three most relevant Mediterranean biomes—forests, woodlands and scrubland vegetation. According to the classification of World Wide Fund for Nature, the Mediterranean forests, woodlands and scrub biome (PA1214 ecoregion), together with four other scrubland types, support more than 25 000 species of plants, of which more than half are endemic. Similarly to other countries, in Portugal this Eco-region is threatened by urbanization, agriculture, frequent fires and invasive species (Olson and Dinerstein 1998). This is the first study providing a comprehensive analysis of soil biodiversity in this type of ecosystem. MATERIALS AND METHODS Study sites and sampling The Serra da Arrábida is a small chain of limestone outcrops with a maximum elevation of 500 m running in an east-west direction, parallel to the Southern edge of the Setúbal peninsula (Fig. 1) and falling off as steep cliffs into the sea. It encompasses different areas with specific formations of typical Mediterranean vegetation of great interest, which has led to the creation of the Arrábida Natural Park in 1976 (http://www.icnf.pt/portal/ap/p-nat/pnar). The region is generally characterized by an Atlantic-Mediterranean climate and consists of different microclimatic areas determined by the varying orientation of the landscape and orography. The area is characterized by a Mediterranean climate with 14–16◦ C average temperatures and 600–700 mm annual precipitation (Guerreiro 2008). Two rather contrasting habitats can be found on the northern and southern slopes of the Arrábida mountain range, a humid site with more pronounced Atlantic influence Palmela Lisbon Setúbal Arrábida Natural Park Azeitao 10 km rea da mi Hu y Dr to id um 50 km bh Su Sesimbra 1 km N1: Humid, forest Figure 1. Location of study sites in the Serra da Arrábida Natural Park. S1: Subhumid, forest S2: Dry, shrubland Yurkov et al. located on the northern slope and a subhumid to semiarid site on the southern slope. Soils are red Mediterranean soils on calcareous or dolomitic bedrock, specifically Rendzic Leptosols associated with Rendzic Cambisols on the northern slope, and Calcic Cambisols with rock outcrops on the southern slope (Schrittenlocher 1997). According to the information available (Schrittenlocher 1997; Guerreiro 2008), topsoils under humid and semihumid forest are moist (Topographic Wetness Index (Beven and Kirkby 1979), TWI = 13–17) and are characterized by high values of available water capacity (140–200 mm) at root depth. In contrast, soils under scrubland vegetation are dry (TWI 3–5) and have poor available water capacity of 50–100 mm (Schrittenlocher 1997; Guerreiro 2008). Soils in the areas studied are slightly acid to neutral, pH = 6–7 (Schrittenlocher 1997). Other climatic parameters and basic soil properties of the study area have been published elsewhere (Schrittenlocher 1997; Guerreiro 2008). Differences in precipitation have a strong effect on vegetation cover resulting in a gradual change of vegetation from humid forest (northern slope) to semihumid chaparral forest (southern slope, valleys) and dry patchy maquis scrubland (southern slope, open areas). The predominant plant cover consists of forests and shrub formations largely composed of the same plant species, namely oaks (Quercus coccifera, Q. faginea, Q. ilex, Q. suber), Montpellier maple (Acer monspessulanum), juniper (Juniperus navicularis), wild olive (Olea europaea), turpentine tree (Pistacia lentiscus), tree heath (Erica arborea), strawberry tree (Arbutus unedo), cistus (Cystus albidus) and phillyrea (Phillyrea latifolia). These plant species provide a good representation of the ancient Mediterranean vegetation and may probably represent remnants of the primordial forest formations that developed under the specific bioclimatic conditions of Arrábida. In September 2013, soils were collected from three different areas of the Serra da Arrábida, in humid forest (northern slope, plot N1), in semihumid chaparral forest (Mata do Solitário, plot S1) and in maquis scrubland (Alto do Jaspe, plot S2) (Fig. 1). On each 20 × 20 m plot, four 10-m long transects were made from the centre of the plot to the four cardinal directions (north, east, south and west). Soils were collected in equal proportions every 2 m along a transect resulting in five soil samples combined in a composite sample (100–200 g in total) and making a total of four composite soil samples per plot. Coarse woody debris, roots and stones (>5 mm) were aseptically removed in the field. Soil samples were placed in sterile paper bags, transferred to the laboratory and kept at –80◦ C before analysis. Prior to the analysis, soils were slowly thawed at 4◦ C and sieved sequentially through 5, 2 and 1 mm meshes. Roots, stones and woody particles were thereby removed in every step. Five subsamples were randomly taken (20 subsamples per plot), from the four composite samples collected on a single plot, suspended in sterile water and plated on three plates each, 180 plates in total. The use of composite soil samples gave reliable results in previous surveys (e.g. Yurkov, Kemler and Begerow 2011) and was also tested on a few samples for Serra da Arrábida soils in the present study (data not shown). Isolation of cultures Frozen soil samples were thawed at 4◦ C overnight and placed in 50 ml plastic tubes, suspended 1:5, 1:10 and 1:20 (w/v) in sterile water and shaken on an orbital shaker at 200 rpm for 1 hour. An aliquot of 0.15 ml was distributed on the surface of glucose-yeast extract-peptone agar acidified with lactate to 3 pH 4.5 (Yurkov, Kemler and Begerow 2011). Plates were incubated at room temperature for 2 days and then at lower temperature (4◦ C) to prevent fast development of moulds. Plates were checked after 7, 14 and 21 days of incubation. Colonies were differentiated into macromorphological types using dissection microscopy, counted and one to two representatives of every colony type per plate were purified. Identification of cultures Detailed protocols describing DNA extraction, amplification, purification and sequencing are given by Yurkov, Kemler and Begerow (2011) and Yurkov, Kemler and Begerow (2012). In short, PCR fingerprinting with minisatellite-specific oligonucleotides derived from the core sequence of bacteriophage M13 with the sequence given by Sampaio et al. (2001) or microsatellite-specific oligonucleotides (GTG)5 , (ATG)5 and (GAC)5 as single PCR primer (Gadanho and Sampaio 2002) were used to group pure cultures. Strains showing identical electrophoretic profiles were considered as conspecific and only one to two representatives of them were chosen for further identification by sequencing of rRNA gene regions. Yeast cultures were identified using nucleotide sequences of the D1/D2 domains of the large subunit (26S/28S or LSU) rRNA gene and in a few cases the internal transcribed spacer (ITS) region was also used. The nucleotide sequences were compared with sequences deposited in the NCBI (www.ncbi.nih.gov) and CBS (www.cbs.knaw.nl) databases, respectively. Nucleotide sequences were deposited in GenBank under the accession numbers given in Table 1. Statistical data analyses Yeast quantity and community structure were determined for each subsample. Yeast quantity was calculated as colony forming units (CFU) per gram of soil at natural humidity. Frequency of occurrence was calculated as the number of samples where a species was observed as a proportion of the total number of samples. Relative abundance was calculated as proportion of a particular species in the sample and is based on colony counts. Species diversity was assessed using Shannon diversity index and species evenness was assessed using Pielou community evenness index. Additional details on these methods are given by Yurkov, Kemler and Begerow (2012). Cultivation results from a total of 428 plates were included in the final analysis. Results from other subsamples and culture plates were not included in the final dataset, as they yielded no yeast cultures either due to low fungal quantity in a particular replicate or due to fast development of moulds, which made isolation and appropriate quantification of yeasts difficult. Species accumulation curves were calculated with EstimateS 9.1.0 using 100 randomizations and sampling without replacement, and the upper limit for infrequent species was set to five samples (Colwell et al. 2012). As distinct yeast species could form colonies with similar morphology and thus make the separation into different types doubtful, we used only presence/absence (incidence data) in our community matrix. The latter did not depend on the morphological differentiation and relied solely on molecular species identification. Four estimators of species richness were used: Chao 2 richness estimator, ICE incidencebased coverage estimator, first-order Jackknife richness estimator (Jackknife 1) and Bootstrap richness estimator. Other details on the use of species richness estimators are given by Yurkov, Kemler and Begerow (2011). Cryptococcus aerius Cryptococcus aerius 2 Cryptococcus terreus Cryptococcus sp. 8 Cryptococcus uzbekistanensis Cryptococcus sp. 8 Cryptococcus arrabidensis Cryptococcus sp. 1 Cryptococcus sp. 2 Cryptococcus sp. 9 Cryptococcus magnus Cryptococcus oeirensis Cryptococcus sp. 7 Cryptococcus laurentii Cryptococcus terrestris Cryptococcus carnescens Cryptococcus heimaeyensis Cryptococcus sp. 3 Cryptococcus sp. 4 Cryptococcus cuniculi Cryptococcus podzolicus Cryptotrichosporon sp. 1 Trichosporon moniliiforme Lineage Filobasidiales Filobasidiales Filobasidiales Filobasidiales Filobasidiales Filobasidiales Filobasidiales Filobasidiales Filobasidiales Filobasidiales Filobasidiales Filobasidiales Tremellales Tremellales Tremellales Tremellales Tremellales Tremellales Tremellales Tremellales Tremellales Trichosporonales Trichosporonales Pezizomycotina Pezizomycotina Pezizomycotina Saccharomycotina Saccharomycotina Saccharomycotina Saccharomycotina Saccharomycotina Saccharomycotina Saccharomycotina Saccharomycotina Saccharomycotina Saccharomycotina Saccharomycotina Basidiomycota, Agaricomycotina Aureobasidium namibiae Aureobasidium pullulans Dothidea sp. Barnettozyma californica Candida magnoliae Candida parapsilosis Candida quercitrusa Candida railenensis Debaryomyces hansenii Kuraishia molischiana Lachancea thermotolerans Meyerozyma guilliermondii Schwanniomyces polymorphus Torulaspora delbrueckii Ascomycota Indentification results Aerius Aerius Aerius Aerius Albidus Albidus Arrabidensis Arrabidensis Arrabidensis Arrabidensis Floriforme Floriforme Amylolyticus Bulleromyces Bulleromyces Dimennae Dimennae Filobasidiella Filobasidiella Kwoniella Podzolicus Cryptotrichosporon Cutaneum Barnettozyma Starmerella Lodderomyces Kurtzmaniella Kurtzmaniella Debaryomyces Kurashia Lachancea Meyerozyma Schwanniomyces Torulaspora Clade OR 177 OR 727 OR 263 OR 721 OR 812 OR 67 OR 255 OR 78 OR 402 OR 262 OR 759 OR 826 OR 975 OR 333 OR 780 OR 129 OR 313 OR 918 OR 849 OR 338 OR 981 OR 395 OR 562 OR 431 OR 756 OR 1009 OR 575 OR 46 OR 435 OR 191 OR 687 OR 547 OR 554 OR 967 OR 506 OR 363 OR 446 Strain ID PYCC 6815 PYCC 6809 PYCC 6483 PYCC 6814 PYCC 6805 PYCC 6816 PYCC 6486 PYCC 6482 PYCC 6481 PYCC 6817 PYCC 6806 PYCC 6810 PYCC 6480 PYCC 6919 PYCC 6812 PYCC 6477 PYCC 6813 PYCC 6811 Collection ID Taxon example KT225566 KT225567 KT304207 KT314198 KT253546 KT253545 KT253538 KT314191 KT314192 KT253539 KT304206 KT253541 KT314197 KT304205 KT253546 KT253540 KT304204 KT314193 KT314194 KT304202 KT253542 KT314199 KT304216 KT304200 LN906592 KT314203 KT225561 KT304199 KT225562 KT225563 KT225565 KT304211 KT253550 KT253551 KT261409 KT253548 KT261411 GenBank KC433824 FR716531 HF558655 HF558655 JX188136 JX188127 AF181535 EU002810 FN428974 FN428974 JX188126 AM160646 KM555196 AY315663 FN428902 HG421433 HG421432 JN939462 JN939461 EU002791 FN428938 KP020115 FN868265 FN824492 FJ150942 GQ153222 FR772338 KJ630501 HF545671 KF830174 JX188108 HF545662 DQ128167 CU928180 JX188191 FN545820 HE799671 Closest match 99%, AF075486 99%, AF075486 99%, AF075479 99%, AF075479 99%, AF181508 97%, FJ473376 100%, AF181535 97%, AF181535 98%, AF181535 98%, AF181535 99%, AF181851 100%, AF181519 94%, AB261011 99%, FN689393 100%, LK023756 100%, AB035054 100%, DQ000317 96%, AF075484 99%, AF105391 99%, DQ333885 99%, AF075481 97%, AY550002 100%, HE616749 99%, FJ150937 100%, FJ150942 97%, AY930108 100%, U75957 99%, U45722 100%, U45754 99%, U45831 99%, U45800 100%, JN940504 100%, U70178 99%, JQ689022 99%, U45709 100%, AB054994 100%, AF105391 Closest type GenBank accession numbers 1.04 – 18.57 1.56 1.56 s.i. – 2.08 – – – – – 1.04 s.i. – 1.04 – – s.i. 1.56 – s.i. 5.73 – – 1.04 – s.i. 4.69 – 1.56 s.i. 3.65 18.23 1.56 s.i. Plot N1 4.85 1.21 24.24 1.21 – – 2.42 – – s.i. – – – – – – – – – – – – – 1.21 – – – – – 6.06 1.21 – – s.i. 13.33 s.i. – Plot S1 Frequency, % s.i. – 26.27 4.31 8.24 – 1.57 s.i. – 1.96 s.i. – s.i. s.i. s.i. – 3.92 s.i. – s.i. s.i. – 1.57 3.14 s.i. – 0.78 s.i. – – – – 3.14 s.i. – – Plot S2 Table 1. Yeast species inventory and frequency of occurrence of soil yeasts isolated three forest types. Total quantity, species richness, diversity and community evenness indices are provided below. 4 FEMS Yeast Research, 2016, Vol. 16, No. 1 Rhodotorula slooffiae Curvibasidium cygneicollum Curvibasidium pallidicorallinum Curvibasidium sp.1 Leucosporidium scottii Rhodosporidium babjevae Rhodotorula mucilaginosa Sporobolomyces roseus Rhodotorula sp. 1 Rhodotorula sp. 2 Species richness, N Shannon diversity index, H Pielou evenness index, J Frequent species (>10%), % Rare species (<1%), % Singletons, N Basidiomycota, Pucciniomycotina Cryptococcus ramirezgomezianus Cryptococcus sp. 5 Cryptococcus sp. 6 Holtermanniella takashimae Cystofilobasidium capitatum Cystofilobasidium sp. 1 Cystofilobasidium sp. 2 Cystofilobasidium sp. 3 Cystofilobasidium sp. 4 Guehomyces pullulans Cryptococcus huempii Indentification results Table 1. (Continued). Cystobasidiomycetes Mycrobotryomycetes Mycrobotryomycetes Mycrobotryomycetes Mycrobotryomycetes Mycrobotryomycetes Mycrobotryomycetes Mycrobotryomycetes Mycrobotryomycetes Mycrobotryomycetes Trichosporonales Trichosporonales Trichosporonales Holtermanniales Cystofilobasidiales Cystofilobasidiales Cystofilobasidiales Cystofilobasidiales Cystofilobasidiales Cystofilobasidiales Cystofilobasidiales Lineage Sporidiobolales Cystobasidiales Curvibasidium Curvibasidium Curvibasidium Leucosporidiales Sporidiobolales Sporidiobolales Sporidiobolales Cystofilobasidium Cystofilobasidium Cystofilobasidium Cystofilobasidium Cystofilobasidium Guehomyces Huempii Humicola Humicola Humicola Clade OR 1010 OR 442 OR 860 OR 28 OR 548 OR 337 OR 228 OR 157 OR 209 OR 742 OR 335 OR 1016 OR 701 OR 241 OR 124 OR 392 OR 240 OR 238 OR 878 OR 523 OR 30 Strain ID PYCC 6808 PYCC 6818 PYCC 6807 PYCC 6856 PYCC 6804 PYCC 6484 Collection ID Taxon example KT261410 KT304209 KT304210 KT314200 KT261408 KT304213 KT304214 KT304215 KT314204 KT261412 KT253543 KT314195 KT314196 KT304212 KT253547 KM409623 KM409621 KT314201 KT314202 KT253549 KT304217 GenBank AB052212 JX188149 JX188150 JX188149 AF070419 JX188219 KF411535 AM160644 KC006703 KF411559 HF558654 KP020115 FN824491 JX188168 JF501392 FN667853 EU002815 KC006882 KC006882 FN868259 AF189844 Closest match 99%, AF189965 99%, AF189928 100%, AF444736 99%, AF189950 100%, AF070419 99%, AF070420 100%, AF070432 100%, AF070441 91%, AY212993 95%, AF070432 100%, AB126584 97%, AY550002 97%, AB126586 100%, FM242574 99%, AF075465 98%, AF075465 98%, EU085532 91%, EU085532 91%, AF075465 100%, AF105394 99%, AF189844 Closest type GenBank accession numbers – 2.6 7.81 1.56 1.56 s.i. 8.33 – – – 31 3.08 0.87 6 26 9 5.73 – – – – – – – – 2.6 s.i. Plot N1 – 1.82 3.03 s.i. – – 9.09 s.i. s.i. – 25 2.81 0.86 8 32 8 – – 4.24 s.i. 1.21 6.67 5.45 s.i. – 7.88 – Plot S1 Frequency, % s.i. – 9.02 9.8 – s.i. 13.33 s.i. – 0.78 34 2.92 0.8 6 59 17 s.i. s.i. – – 1.57 2.35 – – 1.18 – 0.78 Plot S2 Yurkov et al. 5 6 FEMS Yeast Research, 2016, Vol. 16, No. 1 Candida parapsilosis Cryptococcus uzbekistanensis Curvibasidium cygneicollum Cryptococcus terrestris Schwanniomyces polymorphus Cryptococcus ramirezgomezianus Guehomyces pullulans Cryptococcus huempii N1 Rhodosporidium babjevae Cryptococcus sp. 1 Candida quercitrusa Cryptococcus sp. 8 31 species Aureobasidium namibiae Cryptococcus aerius Cryptococcus terreus S2 Curvibasidium pallidocorallinum S1 Curvibasidium sp. 1 Lachancea thermotolerans Meyerozyma guilliermondii Rhodotorula mucilaginosa 34 species 25 species Incidence 15 - 11 % Cystofilobasidium capitatum 10 - 5 % Cystofilobasidium sp. 1 4-2% Sporobolomyces roseus 1-0% Figure 2. Total species richness, dominating and shared yeasts species isolated from the three plots. Grey shades show average frequency of occurrence of common yeast species. Colour codes follow those of Fig. 1. RESULTS All analysed soils yielded yeasts. Highest total yeast counts values varied from 1.3 × 105 (plot S1) to 1.3 × 106 (plot S2) CFU per gram of soil at natural humidity. A total of 613 cultures were isolated and identified to 57 yeast species, 13 belonging to the Ascomycetes and 44 to the Basidiomycetes (Table 1). Only eight species (14%) occurred in all three sampling sites. The highest number of species recovered in a single plot was 34 (59% of the total number of species found). These values appear to reflect contrasting environmental conditions in the study area. The most species-rich community (34 species) was that of soil under dry maquis scrubland (plot S2) followed by humid forest community (31 species, plot N1) and that of semihumid chaparral forest (25 species, plot S1). In all three plots, basidiomycetous yeasts prevailed over ascomycetous species (Fig. 2). The most frequent basidiomycetous species were Cryptococcus terreus (plot N1 and S1), Rhodotorula mucilaginosa (plot S1) and Cr. aerius (plot S2). Among ascomycetous yeasts, Candida quercitrusa (plots N1 and S1), Meyerozyma guilliermondii (plots N1 and S1), Aureobasidium namibiae (plots N1 and S2) and Lachancea thermotolerans (plots N1 and S2) were the most frequent species (Fig. 2). The Pielou evenness index that compares the proportion of the different species in the community ranged between 0.87 and 0.80, which indicates that structure of the yeast communities in the three plots was not even (Table 1). The class of most frequent species (those with >10% of isolates) had the same number of representatives in all three plots: two species or 6–8% of the total richness. While the majority of species in forested sites (N1 and S1) were found at frequencies ranging from 1% to 10%, 21 species (68% of species) in plot N1 and 15 species (60% of species) in plot S1, the majority of species (20 species, 59%) isolated from dry scrubland (plot S2) were rare (often singletons) and were found with a frequency of less than 1%. Shannon diversity index reflects both species richness and community structure. The highest diversity value was recorded for humid forest soil (H = 3.08, plot N1) followed by dry maquis scrubland (H = 2.92, plot S2) and semihumid chaparral forest (H = 2.81, plot S1). Despite the substantial difference in the observed species richness values between the two plots located on the southern slope (25 and 35 species in plots S1 and S2, respectively), the Shannon diversity indices calculated for these plots were very close and reflected little contribution of rare species to soil yeast communities. Diversity estimations using Chao 2, ICE, first-order Jackknife and Bootstrap estimators predict up to 80 yeast species to be present from all three investigated plots (Fig. 3). Species richness curves obtained with Chao 2 and ICE estimators were close to saturation after the given sampling effort. Chao 2 estimator resulted in the highest diversity values. The observed value of 60 species laid within 75%–90% of the predicted species richness values (depending on the estimator used) showing that our study provides a reliable basis for yeast species richness assessment in the studied soils. Estimations performed for each sampling site suggest that the cultivation approach used s results in species rarefaction curves that are close to saturation for soils in the two forested areas (plots N1 and S1) but not for soil Yurkov et al. 50 All Plots Chao 2 Jack 1 ICE 80 Bootstrap 60 40 Number of species, N Number of species, N 100 20 0 Plot N1 Jack 1 ICE Chao 2 Bootstrap 40 30 20 10 0 50 100 150 200 0 250 Number of sub-samples 40 0 70 Plot S1 25 50 75 Number of sub-samples 100 Plot S2 Chao 2 Jack 1 ICE Chao 2 Bootstrap 30 20 10 Number of species, N Number of species, N 7 60 ICE Jack 1 50 Bootstrap 40 30 20 10 0 0 25 50 75 Number of sub-samples 100 0 0 25 50 75 Number of sub-samples 100 Figure 3. Randomized species accumulation curves for the total soil yeast community of the Serra da Arrábida soils and for each of the three plots (N1, S1 and S2) obtained with incidence-based coverage (ICE), Chao 2, first-order Jackknife (Jack 1) and bootstrap richness estimators. Shadowed areas correspond to the 75% expected diversity cut-off. Colour codes follow those of Fig. 1. associated with scrubland vegetation (plot S2) (Fig. 3). Species richness estimators calculated similar values for both forest types (Fig. 3), i.e. 40–43 species in plot N1 and 29–34 species in plot S1, respectively. In soils under scrubland vegetation, estimated species richness values varied largely from 45 to 64 species reflecting a contradictory contribution of rare species to the total diversity. Overall, species richness estimators performed differently depending on the soil yeast community structure. Increasing number of rare species resulted in higher species richness values predicted by Chao 2 and ICE estimators, which distinguish between rare and frequent species. A total of 57 yeast species were isolated and identified during this study. They belong to four lineages of Fungi, Pezizomycotina (2 species), Saccharomycotina (11 species), Agaricomycotina (34 species) and Pucciniomycotina (10 species). According to the genetic distances obtained with sequences deposited in NCBI GenBank and MycoBank databases, 20 yeast taxa represented potential new species (Table 1). Detailed analysis and formal description of these new species will be the aim of future studies. While many putative new species were isolated in low frequencies, a few species (e.g. Cystofilobasidium sp. 1 in S1 and S2 plots, and Cystofilobasidium sp. 2 in S1 plot, respectively) were rather frequent and were isolated from different plots. DISCUSSION Current estimations suggest that there are 1.5–5.1 million fungal species, from which only a small part (about 100 000 species) has been described (Blackwell 2011; Taylor et al. 2014). Furthermore, in some substrates such as plants and soils the number of expected fungi may exceed the number of plants by 10:1 (Blackwell 2011). Presently, culture-independent surveys of soil fungi receive much attention as the availability of high-throughput sequencing techniques is constantly increasing. Similarly to other groups of soil microorganisms, diversity values of fungal taxa (expressed as molecular OTUs) are high (e.g. Bueé et al. 2009; Wubet et al. 2012; Vořı́šková et al. 2014). Phylogenetic lineages represented by fungi with yeast stages are regularly present in sequence (or clone) libraries but despite their high occurrence they are rarely identified to the species level and instead are assigned to large polyphyletic genera such as Cryptococcus and Rhodotorula. In other cases, environmental sequences matching those of the tremellomycetous yeasts Cr. aerius, Cr. terreus, Cr. terricola, Cr. podzolicus, Trichosporon dulcitum, T. porosum and Guehomyces pullulans are among the most frequent soil OTUs revealed in culture-independent surveys (e.g. Lynch and Thorn 2006; Bueé et al. 2009; Taylor et al. 2014; Vořı́šková et al. 2014). All these widespread pedobionts were also isolated during numerous culture-dependent studies (e.g. Wuczkowski and Prillinger 2004; Vishniac 2006; Mestre et al. 2011; Yurkov, Kemler and Begerow 2011, 2012; Yurkov et al. 2012). Therefore, in spite of the growing number of soil biodiversity assessments using high-throughput sequencing techniques, none of them has yet reported a substantial uncultured diversity within the different phylogenetic lineages that include yeasts. 8 FEMS Yeast Research, 2016, Vol. 16, No. 1 This study provides the first soil yeast species inventory in the Mediterranean forests, woodlands and scrub biome, which is one of the major Mediterranean biomes. Our study yielded yeasts that are frequently reported from soils worldwide, namely Cr. aerius, Cr. terreus, Cr. podzolicus, G. pullulans and T. moniliiforme, all members of the Basidiomycota, and Barnettozyma californica, Debaryomyces hansenii, M. guilliermondii and Schwanniomyces polymorphus from the Ascomycota (e.g. Botha 2006; Vishniac 2006; Yurkov, Chernov and Tiunov 2008; Mestre et al. 2011; Yurkov, Kemler and Begerow 2012). A few species isolated in the present study provide a link between above- and belowground communities. The following yeasts were previously found on plant material collected in Serra da Arrábida: Cr. carnescens, Cr. flavescens, Cr. magnus, Cr. oierensis, Curvibasidium cygneicollum, Cystobasidium slooffiae, Rhodosporidium babjevae, R. mucilaginosa and Sporobolomyces roseus (Inácio 2003; Sampaio et al. 2004). These yeasts were isolated from aboveground in high numbers and represent in our opinion typical transient species that are introduced in the soil habitat with falling leafs or other dying plant material. The two ascomycetous yeasts C. quercitrusa and C. railenensis can reflect the type of forest as they were previously reported in association with oaks (e.g. Isaeva et al. 2009; Lachance et al. 2011), and L. thermotolerans was repeatedly recovered from oak bark using a selective medium (Sampaio and Gonçalves 2008). A few yeasts found in Serra da Arrábida soils were previously isolated by other scientists in Portugal and are deposited in PYCC, e.g. Cr. aerius, Cr. podzolicus, Cr. terreus, Cu. cygneicollum, L. thermotolerans, M. guilliermondii and S. occidentalis. Analysis of nucleotide sequences demonstrated that out of the 57 yeast species that we found, 20 (35%) showed low identity to any known species and may thus represent novel taxa. These results are in agreement with the previous observation by Yurkov, Kemler and Begerow (2012) who showed beech forest soil in Germany to harbour many undescribed species. Similarly to other studies, yeasts of the genus Cryptococcus accounted for the highest diversity and were members of the following clades and lineages: Albidus, Magnus (Floriforme), Aerius, Arrabidensis in Filobasidiales, and Aureus, Filobasidiella, Kwoniella and Victoriae in Tremellales (Table 1). Remarkably, this study yielded two new species, represented by strains OR849 and OR918, showing close affinity to Cr. sensu stricto, i.e., the Filobasidiella clade, the group comprising the human pathogens Cr. neoformans and Cr. gattii. Nucleotide sequences of strains of two typical soil-related species Cr. aerius and Cr. terreus also showed some deviation from those of type strains. However, a reliable analysis of relationships within these two species groups is difficult as ribosomal gene fragments show little interspecific divergence (Scorzetti et al. 2002; Scholten 2011; A Yurkov unpublished data). According to our results, strains of Cr. aerius and Cr. terreus show as many as 12 and 10 nucleotide substitutions in LSU (D1/D2 domains of 26S) rRNA, respectively. The first common DNA-Barcode for Fungi, the ITS region, is less variable than LSU as shown by Schoch et al. (2012) and our own results, and shows four and three nucleotide substitutions in Cr. aerius and Cr. terreus, respectively. Thus, a combination of the two genetic markers traditionally used in yeast biodiversity assessments may be insufficient to resolve yeast species presently comprised within Cr. aerius and Cr. terreus. This assumption is supported to some extent by results of the recent multilocus sequence analysis of Cr. flavescens –Cr. terrestris species complex (Yurkov et al. 2015b). Analysis of five loci including genes involved in sexual reproduction showed the limited utility of ribosomal ITS and LSU sequences for species discrimination and resulted in the discovery of two new cryptic species. Studies by Fonseca, Scorzetti and Fell (2000) and Scorzetti et al. (2002) made a solid basis for sequence-based species delimitation in Filobasidiales, including species in the Aerius clade. Our observations suggest that relationships within this group need to be resolved with a multigene approach in the future. Another interesting observation is the unexpected high occurrence and diversity of yeasts of the genus Cystofilobasidium. Besides Cy. capitatum, a species frequently reported from forest litter and dried plant material (e.g. Babeva and Chernov 1995), the present study yielded several potential new species (Fig. 4). Some of these yeasts were observed in high numbers, thus questioning their allochtonous origin in soil, as previously suggested for Cy. capitatum (Babeva and Chernov 1995). Specifically Cystofilobasidium sp. 1 and 2 belong to the group of frequent species isolated from soils under S1 and S2 vegetation. Furthermore, queries in the PYCC database (http://pycc.bio-aware.com/) yielded a strain conspecific to Cystofilobasidium sp. 2 that was also isolated in Serra de Arrábida by Joao Inácio (PYCC database, GenBank EU002815). This suggests that yeasts of the genus Cystofilobasidium may be regular inhabitants of Mediterranean forests, woodlands and scrub biome in Portugal. The yeast community of Mediterranean soils is hyperdiverse The term biodiversity hotspot was introduced to distinguish biogeographic regions with a significant reservoir of biodiversity, total richness and proportion of endemic species (e.g. Myers et al. 2000). In Europe, the Mediterranean basin has been recognized as the primary biodiversity hotspot with 45% of endemic species of plants and vertebrates (Olson and Dinerstein 1998; Myers et al. 2000). Unlike plants and animals, microorganisms inhabiting natural biotopes in this region received little attention. The present study reports the first diversity values for Mediterranean forests, woodlands and scrub biome. Yeasts are found in soils worldwide and data available from previous studies suggest that soil yeast communities are species poor in a single location (low alpha-diversity) but very dissimilar on a broader spatial scale (high beta-diversity). To date, analysis of forest soils yielded 18–26 yeast species isolated from a single biotope, as discussed in Yurkov, Kemler and Begerow (2011). Variation among these values can be due to the different methods used for yeast identification but cryptic species are not considered to greatly influence these numbers (Yurkov, Kemler and Begerow 2011; Yurkov et al. 2015a). In the present study, a few Cryptococcus spp. members of Arrabidensis, Humicola and Kwoniella clades could be considered as cryptic species but these yeasts were rarely observed in the same plots (Table 1). On average, boreal forest soils are characterized by lower yeast species richness values than broadleaf forest soils (Maksimova and Chernov 2004; Wuczkowski and Prillinger 2004; Yurkov, Kemler and Begerow 2012). In the present study, we found 57 yeast species, a value that is 2- to 3-fold higher than the values reported before, namely 20–25 species in Russia (Maksimova and Chernov 2004; Golubtsova, Glushakova and Chernov 2007; Glushakova, Kachalkin and Chernov 2015), 18 species in Slovakia (Sláviková and Vadkertiová 2000), 19–24 species in Austria (Wuczkowski and Prillinger 2004), 10–26 species in Germany (Yurkov, Kemler and Begerow 2011, 2012; Yurkov et al. 2012) and 22 species in Argentina (Mestre et al. 2011). Because Yurkov et al. 9 OR158 OR273 OR242 OR185 OR392 69 OR407 Cystofilobasidium sp. 1 OR694 OR394 99 OR223 OR428 Cystofilobasidium “capitatum” MOM_884 (HG421444) 95 OR398 Cystofilobasidium capitatum CBS 6358T (AF075465) 87 97 90 99 74 OR124 OR408 OR411 Cystofilobasidium feraegula CBS 7201T (AF075487) Cystofilobasidium lacus-mascardii CBS 10642T (AY15864) Cystofilobasidium macerans CBS 2206T (AF189848) Cystofilobasidium infirmo-miniatum CBS 323T (AF075505) 68 Cystofilobasidium bisporidii CBS 6346T (EU085532) OR934 100 OR878 99 100 Cystofilobasidium sp. 4 Cystofilobasidium “capitatum” HAI-Y-699 (KC006882) OR238 98 Cystofilobasidium sp. 3 Cystofilobasidium sp. 8C5 (KM275210) Cystofilobasidium sp. EX-F1547 (DQ640493) 0.02 100 Cystofilobasidium sp. K146b (HF934016) OR661 OR233 OR235 Cystofilobasidium sp. 2 OR240 OR282 Cystofilobasidium sp. PYCC 6093 (EU002815) “Mrakia” curviuscula CBS 9136T (EF118826) Figure 4. Maximum likelihood analysis (RAxML, GTRCAT option) of an alignment of the LSU (D1/D2 domains) rRNA gene for Cystofilobasidium spp. strains isolated during this study. The numbers given on branches are frequencies (>50%) with which a given branch appeared in 100 bootstrap replications. The scale indicates the number of expected substitutions accumulated per site. The tree is rooted with Mrakia curviuscula CBS 9136T . a direct comparison between studies is difficult due to different species recognition approaches and sampling effort, the most reliable comparisons of yeast diversity can be achieved using a restricted number of recent studies (e.g. Wuczkowski and Prillinger 2004; Vishniac 2006; Mestre et al. 2011; Yurkov, Kemler and Begerow 2011, 2012; Yurkov et al. 2012; Glushakova, Kachalkin and Chernov 2015). The only study to date utilizing species richness estimators reported 14–17 yeast species to be expected from the analysis of soil in a single beech forest (Yurkov, Kemler and Begerow 2011). Another study that used the same sampling approach that we used here reported the isolation of 24 species from beech forests in three regions of Germany (Yurkov, Kemler and Begerow 2012). Accordingly, species richness estimations suggest substantially higher yeast species richness in the Mediterranean soils that we analysed than in soils under broadleaf vegetation in Central Europe. The present study performed in a single biotope representing Mediterranean forests, woodlands and scrub biome resulted in about 3-fold higher species richness values than those reported for boreal and temperate forests. Tropical forests are well recognized for their extraordinary diversity of plants and animals (e.g. Myers et al. 2000); however, our knowledge on soil yeasts in these habitats is scarce. The only recent study reporting moderate yeast diversity values, namely up to 9–11 species isolated from soils underneath four cloud and rain forests in Costa Rica (Vishniac 2006). 10 FEMS Yeast Research, 2016, Vol. 16, No. 1 The influence of climatic conditions on soil yeast diversity Near-natural and unmanaged forests are often characterized by high heterogeneity of forest and soil cover. Forest alteration through management activities cause substantial changes also in belowground microbial communities of broadleaf forests and affects both prokaryotes and fungi, including yeasts (Birkhofer et al. 2012; Yurkov et al. 2012; Yurkov, Kemler and Begerow 2012). Species richness estimations performed in the present study showed strong differences between plots and suggested a total of 80 yeast species to be expected from soils in Serra da Arrábida Natural Park. Based on species richness estimations, maquis scrublands are expected to have the highest diversity, followed by humid and subhumid forests. The analysis of soil yeast community structure showed that the occurrence of rare species accounted for higher species richness estimations for scrubland soil. Specifically, species-rank distribution and distribution of shared species showed that the number of frequent species (>10% of total isolates) was equal among plots and that the frequent species corresponded to 6–8% of the total richness (Table 1). At the same time the number of species occurring with 1–10% frequency declined from humid forest to chaparral forest and scrubland. Rare species (<1% of total isolates) were prominent (59% of the total richness) in soils under dry maquis scrubland. Out of 20 rare species detected in scrubland soils, 17 were singletons and this is twice the number of singletons found in forest plots (Table 1). These differences in species occurrence also resulted in decreasing community evenness values from humid forest (J = 0.87) to scrubland vegetation (J = 0.8) (Table 1). As mentioned before, the three plots were selected to depict a gradual change of environmental conditions, especially the level of precipitation. Although the studied plots are characterized by highly similar plant communities, strong differences in water availability shaped these phytocenoses, and so they differed substantially in properties such as aboveground biomass and plant projective cover. The humid forest on the northern slope (plot N1) was very dense with well-developed understory and decaying wood logs. The semihumid chaparral forest on the southern slope (plot S1) was located in a valley that on average provided plants with more water than the neighbouring scrubland areas. Although large trees dominated semihumid forest, it was more spaced with little understory. In scrublands, the plants grew in patches and were separated by open land and rocks. Plants were often suppressed and trees did not develop their full size although they were still able to grow under such conditions, e.g. oaks, strawberry trees and wild olives. A few studies have addressed the importance of forest type (dominating tree species) and forest properties (management history, presence of decaying wood, etc.) on soil yeast communities. In the present study, we addressed the influence of local climatic parameters that do not alter plant species composition but change vegetation properties such as biomass and plant projective cover. The latter can be also considered as a measure of substrate discontinuity for microfungi as they do not form abundant hyphae and, thus, are likely to be restricted to a local pool of nutrients in a microhabitat. Our results suggest that the discontinuity of the vegetation cover may affect soil microorganisms, exemplified by yeasts, shifting the structure of the community from the dominance of frequent species to the prevalence of minor and, possibly, transient species. In the present study, we provide evidence that supports the hypothesis that local climatic conditions can lead to a substantial alteration of soil yeast communities even on a small spatial scale. CONCLUSIONS Soil-related basidiomycetous yeasts have been recently recognized as promising candidates for biotechnological applications, e.g. as enzyme and oil producers. Despite this relevance, limited sampling hampers a proper diversity assessment and consequently the utilization of novel species and strains. The sampling technique that we have optimized in previous studies provides a solid basis for the efficient recovery of soil yeasts in different biotopes. The analysis of cultivable soil yeasts in the Mediterranean forest, woodland and scrub biome revealed remarkably high yeast diversity values, a finding that is in line with the current view of the Mediterranean basin as a biodiversity hotspot. Results obtained in the present study challenge the present knowledge of soil yeast communities as species-poor assemblages and suggest that diversity patterns follow forest cover properties and habitat fragmentation. ACKNOWLEDGEMENTS Authors are grateful to Álvaro Fonseca for his assistance in plot selection and sampling and Brian J. Tindall for correcting the English text. FUNDING This work was financially supported by the Fundação para a Ciência e a Tecnologia, Portugal grant numbers PTDC/BIABIC/4585/2012 and UID/Multi/04378/2013. Conflict of interest. None declared. REFERENCES Babeva IP, Chernov IY. Geographic aspects of yeast ecology. Physiol Gen Biol Rev 1995;9:1–54. Beven KJ, Kirkby MJ. A physically based, variable contributing area model of basin hydrology. Hydrolog Sci J 1979;24: 43–69. Birkhofer K, Schöning I, Alt F, et al. General relationships between abiotic soil properties and soil biota across spatial scales and different land-use types. PLoS One 2012;7:e43292. Blackwell M. The fungi: 1, 2, 3 . . . 5.1 million species? Am J Bot 2011;98:426–38. Botha A. Yeasts in soil. In: Rosa CA, Peter G (eds). Biodiversity and Ecophysiology of Yeasts. The Yeast Handbook. Berlin-Heidelberg, Germany: Springer, 2006, 221–40. Botha A. The importance and ecology of yeasts in soil. Soil Biol Biochem 2011;43:1–8. Bouthilet RJ. A taxonomic study of soil yeasts. Mycopath Mycol Appl 1951;6:79–85. Bueé M, Reich M, Murat C, et al. 454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity. New Phytol 2009;184:449–56. Buzzini P, Margesin R. Cold-Adapted Yeasts. Berlin-Heidelberg: Springer, 2014. Chernov IY. The latitude-zonal and spatial-successional trends in the distribution of yeasts. Zh Obshch Biol 2005;66:123–35. Colwell RK, Chao A, Gotelli NJ, et al. Models and estimators linking individual-based and sample-based rarefaction, extrapolation and comparison of assemblages. J Plant Ecol 2012;5:3–21. Di Menna ME. Yeasts in New Zealand soils. New Zeal J Bot 1965;3:194–203. Yurkov et al. Fonseca Á, Inácio J. Phylloplane yeasts. In: Rosa CA, Peter G (eds). Biodiversity and Ecophysiology of Yeasts. The Yeast Handbook. Berlin-Heidelberg, Germany: Springer, 2006, 263–301. Fonseca Á, Scorzetti G, Fell JW. Diversity in the yeast Cryptococcus albidus and related species as revealed by ribosomal DNA sequence analysis. Can J Microbiol 2000;46:7–27. Gadanho M, Sampaio JP. Polyphasic taxonomy of the basidiomycetous yeast genus Rhodotorula: Rh. glutinis sensu stricto and Rh. dairenensis comb. nov. FEMS Yeast Res 2002;2:47–58. Glushakova AM, Kachalkin AV, Chernov IY. Soil yeast communities under the aggressive invasion of Sosnowsky’s hogweed (Heracleum sosnowskyi). Eurasian Soil Sci+ 2015;48:201–7. Golubtsova YV, Glushakova AM, Chernov IY. The seasonal dynamics of yeast communities in the rhizosphere of soddypodzolic soils. Eurasian Soil Sci+ 2007;40:875–9. Guerreiro SR. Contributo para a caracterização e gestão da vegetação da serra da Arrábida. Lisbon: Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa, 2008. Guilliermond A. The Yeasts. New York: John Wiley and Sons, 1920. Inácio J. Yeast Occurrence and Diversity on the Phylloplane of Selected Plants From the Arrábida Natural Park. Lisbon: Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa, 2003. Inácio J, Pereira P, Carvalho DM, et al. Estimation and diversity of phylloplane mycobiota on selected plants in a Mediterranean-type ecosystem in Portugal. Microbial Ecol 2002;44:344–53. Isaeva OV, Glushakova AM, Yurkov AM, et al. The yeast Candida railenensis in the fruits of English oak (Quercus robur L.). Microbiology+ 2009;78:355–9. Kachalkin AV, Yurkov AM. Yeast communities in Sphagnum phyllosphere along the temperature-moisture ecocline in the boreal forest-swamp ecosystem and description of Candida sphagnicola sp. nov. Anton Leeuw 2012;102:29–43. Lachance MA, Boekhout T, Scorzetti G, et al. Candida Berkhout (1923). In: Kurtzman CP, Fell JW, Boekhout T (eds). The Yeasts, A Taxonomic Study, 5th edn. Amsterdam: Elsevier, 2011, 987– 1278. Lynch MD, Thorn RG. Diversity of basidiomycetes in Michigan agricultural soils. Appl Environ Microb 2006;72:7050–6. Maksimova IA, Chernov IY. Community structure of yeast fungi in forest biogeocenoses. Microbiology+ 2004;73:474–81. Mestre MC, Rosa CA, Safar SV, et al. Yeast communities associated with the bulk-soil, rhizosphere and ectomycorrhizosphere of a Nothofagus pumilio forest in northwestern Patagonia, Argentina. FEMS Microbiol Ecol 2011;78:531–41. Myers N, Mittermeier RA, Mittermeier CG, et al. Biodiversity hotspots for conservation priorities. Nature 2000;403:853–8. Olson DM, Dinerstein E. The Global 200: A representation approach to conserving the Earth’s most biologically valuable ecoregions. Conserv Biol 1998;12:502–15. Sampaio JP, Gadanho M, Santos S, et al. Polyphasic taxonomy of the basidiomycetous yeast genus Rhodosporidium: Rhodosporidium kratochvilovae and related anamorphic species. Int J Syst Evol Micr 2001;51:687–97. Sampaio JP, Golubev WI, Fell JW, et al. Curvibasidium cygneicollum gen. nov., sp. nov. and Curvibasidium pallidicorallinum sp. nov., novel taxa in the Microbotryomycetidae (Urediniomycetes), and their relationship with Rhodotorula fujisanensis and Rhodotorula nothofagi. Int J Syst Evol Micr 2004;54: 1401–7. Sampaio JP, Gonçalves P. Natural populations of Saccharomyces kudriavzevii in Portugal are associated with oak bark and are sympatric with S. cerevisiae and S. paradoxus. Appl Environ Microb 2008;74:2144–52. 11 Schoch CL, Seifert KA, Huhndorf S, et al. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. P Natl Acad Sci USA 2012;109:6241–6. Scholten N. DNA-Barcoding und Molekular-Phylogenetische Untersuchungen der Polyhyletischen Bodenhefe Cryptococcus. Bochum: Ruhr-Universität Bochum, 2011. Schrittenlocher R. The soils of the Parque Natural da Arrábida. Finisterra 1997;32:15–43. Schulze I, Hansen S, Großhans S, et al. Characterization of newly isolated oleaginous yeasts - Cryptococcus podzolicus, Trichosporon porosum and Pichia segobiensis. AMB Express 2014;4:1–11. Scorzetti G, Fell JW, Fonseca Á, et al. Systematics of basidiomycetous yeasts: a comparison of large subunit D1/D2 and internal transcribed spacer rDNA regions. FEMS Yeast Res 2002;2:495–517. Sláviková E, Vadkertiová R. The occurrence of yeasts in the forest soils. J Basic Microbiol 2000;40:207–12. Sláviková E, Vadkertiová R. The diversity of yeasts in the agricultural soil. J Basic Microbiol 2003;43:430–6. Starkey RL, Henrici AT. The occurrence of yeasts in soils. Soil Sci 1927;23:33–46. Tanimura A, Takashima M, Sugita T, et al. Cryptococcus terricola is a promising oleaginous yeast for biodiesel production from starch through consolidated bioprocessing. Sci Rep 2014;4:4776. Taylor DL, Hollingsworth TN, McFarland JW, et al. A first comprehensive census of fungi in soil reveals both hyperdiversity and fine-scale niche partitioning. Ecol Monogr 2014;84: 3–20. Vishniac HS. A multivariate analysis of soil yeasts isolated from a latitudinal gradient. Microb Ecol 2006;52:90–103. Vořı́šková J, Brabcová V, Cajthaml T, et al. Seasonal dynamics of fungal communities in a temperate oak forest soil. New Phytol 2014;201:269–78. Wubet T, Christ S, Schöning I, et al. Differences in soil fungal communities between European beech (Fagus sylvatica L.) dominated forests are related to soil and understory vegetation. PLoS One 2012;7:e47500. Wuczkowski M, Prillinger H. Molecular identification of yeasts from soils of the alluvial forest national park along the river Danube downstream of Vienna, Austria (“Nationalpark Donauauen”). Microbiol Res 2004;159:263–75. Yurkov AM, Chernov IY, Tiunov AV. Influence of Lumbricus terrestris earthworms on the structure of the yeast community of forest litter. Microbiology+ 2008;71:107–11. Yurkov AM, Kemler M, Begerow D. Species accumulation curves and incidence-based species richness estimators to appraise the diversity of cultivable yeasts from beech forest soils. PLoS One 2011;6:e23671. Yurkov AM, Kemler M, Begerow D. Assessment of yeast diversity in soils under different management regimes. Fungal Ecol 2012;5:24–35. Yurkov A, Wehde T, Kahl T, et al. Aboveground deadwood deposition supports development of soil yeasts. Diversity 2012;4:453–74. Yurkov A, Inácio J, Chernov IY, et al. Yeast biogeography and the effects of species recognition approaches: the case study of widespread basidiomycetous species from birch forests in Russia. Curr Microbiol 2015a;70:587–601. Yurkov A, Guerreiro MA, Sharma L, et al. Correction: Multigene assessment of the species boundaries and sexual status of the basidiomycetous yeasts Cryptococcus flavescens and C. terrestris (Tremellales). PLoS One 2015b;10:e0126996.
© Copyright 2026 Paperzz