A new Rhizopogon species associated with Pinus amamiana in Japan

Pinus amamiana is an endangered Pinus species found only on Yakushima Island and Tanegashima Island, Japan. We surveyed remaining P. amamiana forests and found some sporocarps of Rhizopogon (Boletales), many species of which exhibit strict host specificity to a narrow range of Pinaceae trees and play critical roles in host establishment. Based on morphological characteristics and molecular phylogeny, here we describe Rhizopogon yakushimensis sp. nov. This new species belongs to a new clade, phylogenetically related to the subgenera Versicolores and Rhizopogon. We also confirmed its ectomycorrhizal association with P. amamiana by comparing rDNA ITS sequences between the sporocarps and ectomycorrhizal root tips.     

Discovery of Rhizopogon associated with Larix from northeastern Siberia: Insights into host shift of ectomycorrhizal fungi

Rhizopogon (Boletales) is an ectomycorrhizal fungal genus that exhibits a strong specificity to Pinaceae. This strict association occurs almost exclusively with Pinus and Pseudotsuga, while associations with other genera in Pinaceae are inconclusive. Here, we describe Rhizopogon laricinus sp. nov. associated with Larix cajanderi distributed in northeastern Siberia, where forest fires are frequent. We confirmed the host identity by comparing rDNA internal transcribed spacer (ITS) sequences obtained from basidiomata and ectomycorrhizal root tips collected at the same sites. Morphological characteristics and molecular identification revealed that R. laricinus is a new species associated with Larix unequivocally. The molecular phylogeny based on ITS sequences placed this species sister to the subgenus Roseoli, which is specific to Pinus, and not to the Pseudotsuga-specific subgenus Villosuli. Thus, R. laricinus evolution does not correspond to host phylogeny as Larix and Pseudotsuga form a monophyletic clade. Instead, ecological traits of Rhizopogon for adapting to disturbed habitats may have driven the host shift to Larix under high-frequency fire ecosystems.     

Revisiting the host effect on ectomycorrhizal fungal communities: implications from host–fungal associations in relict Pseudotsuga japonica forests

Host identity is among the most important factors in structuring ectomycorrhizal (ECM) fungal communities. Both host–fungal coevolution and host shifts can account for the observed host effect, but their relative significance in ECM fungal communities is not well understood. To investigate these two host-related mechanisms, we used relict forests of Pseudotsuga japonica, which is an endangered endemic species in Japan. As with other Asian Pseudotsuga species, P. japonica has been isolated from North American Pseudotsuga spp. since the Oligocene and has evolved independently as a warm-temperate species. We collected 100 soil samples from four major localities in which P. japonica was mixed with other conifers and broadleaf trees. ECM tips in the soil samples were subjected to molecular analyses to identify both ECM fungi and host species. While 136 ECM fungal species were identified in total, their communities were significantly different between host groups, confirming the existence of the host effect on ECM fungal communities. None of the 68 ECM fungal species found on P. japonica belonged to Pseudotsuga-specific lineages (e.g., Rhizopogon and Suillus subgroups) that are common in North America. Most of ECM fungi on P. japonica were shared with other host fungi or phylogenetically close to known ECM fungi on other hosts in Asia. These results suggest that after migrating, Pseudotsuga-specific fungal lineages may have become extinct in small isolated populations in Japan. Instead, most of the ECM fungal symbionts on P. japonica likely originated from host shifts in the region.     

Ectomycorrhizal fungi diversity in a white sand forest in western Amazonia

The genera Dicymbe and Aldina (Fabaceae) host ectomycorrhizal fungi (EcM) and are common in white sand forests (WSFs), a highly specialized habitat with a high level of plant endemism compared with terra-firme forests. In this study, we visited four times a 1-ha permanent plot established in a small patch of a WSF in the south of Colombia Amazonia. Forty-eight species of EcM fungi were recovered from sporocarps and 15 ITS species-level were detected from root tips. Seventeen species were new reports to Colombia and seven corresponded to undescribed species. These results confirm that this WSF supports a significant EcM fungal diversity. Most of the species found in this study have been previously reported to be associated with other legume and/or dipterocarp species from geographically distant forests. The long-distance occurrence combined with low host specificity, suggest the possibility of gene flow between geographically distant populations of EcM fungi in neotropical lowland rainforests.     

DNA-based fungal diversity in Madagascar and arrival of the ectomycorrhizal fungi to the island

Madagascar is known for its high diversity and endemism of fauna and flora. Fungi, however, have been largely overlooked in diversity and evolution studies on the island, and whether fungi exhibit the same patterns as animals and plants has yet to be further examined. We collected fungal sporocarps and ectomycorrhizal (EcM) roots during three opportunistic surveys in five forests in Madagascar and generated a dataset of fungal Internal Transcribed Spacer (ITS) DNA sequences. We analyzed them together with all publicly available fungal ITS DNA sequences and identified 620 Operational Taxonomic Units (OTUs) from Madagascar, 10% of which contained only sequences from our surveys. Two hundred and ninety-two OTUs belonged to EcM species with /russula-lactarius, /boletus, /tomentella-telephora, /cortinarius, and /amanita being the most abundant EcM lineages. Overall, 60% of all fungi and 81% of the EcM species from Madagascar appear to be endemic. We conducted a phylogenetic analysis using all the OTUs in Amanitaceae, Boletaceae, and Russulaceae families to elucidate their relative timing of arrival in Madagascar. We found that most EcM species from Madagascar in the three families diverged less than 22 million years ago (mya), long after the separation of India and Madagascar (88 mya), which is consistent with a dispersal mediated process of arrival to the island. Our study provides the first comprehensive view of the overall DNA-based fungal diversity in Madagascar and the current state of knowledge of EcM fungi based on DNA sequences, useful for further ecological and evolutionary studies. Abstract in Malagasy is available with online material.     

Pleistocene origin and colonization history of Lobelia columnaris Hook. f. (Campanulaceae: Lobelioideae) across sky islands of West Central Africa

We aimed to infer ancestral area and historical colonization of Lobelia columnaris in the sky islands of Bioko and Cameroon through dated phylogeny using chloroplast genomes. Specifically, we aim to answer the following questions: (1) What are the phylogenetic relationships among Bioko Island and Cameroon populations? (2) Are the older populations found in the older sky islands? We assembled novel plastomes from 20 individuals of L. columnaris from 5 mountain systems. The plastome data were explored with phylogenetic analyses using Maximum Likelihood and Bayesian Inference. The populations of L. columnaris have a monophyletic origin, subdivided into three plastomes‐geographic clades. The plastid phylogenomic results and age of the sky islands indicate that L. columnaris colonized first along with the Cameroon Volcanic Line''s young sky islands of Bioko. The crown group (1.54 Ma) split the population in Bioko and mainland Cameroon. It is possible that Bioko was the ancestral area and likely isolated during cold and dry conditions in forest refugia. Presumably, the colonization history occurred during the middle‐late Pleistocene from South Bioko''s young sky island to North Bioko and the northern old sky islands in Cameroon. Furthermore, the central depression with lowland forest between North and South Bioko is a current geographic barrier that keeps separating the populations of Bioko from each other. Also, the shallow sea channel keeps isolated the populations of Bioko and the mainland populations. The Pleistocene climatic oscillations led to the divergence of the Cameroon and Bioko populations into three clades. L. columnaris colonized the older sky islands in mainland Cameroon after establishment in Bioko''s younger sky islands. Contrary to expectations, the biogeography history was an inverse progression with respect to the age of the Afromontane sky islands.     

High diversity and widespread occurrence of mitotic spore mats in ectomycorrhizal Pezizales

Fungal mitospores may function as dispersal units and/ or spermatia and thus play a role in distribution and/or mating of species that produce them. Mitospore production in ectomycorrhizal (EcM) Pezizales is rarely reported, but here we document mitospore production by a high diversity of EcM Pezizales on three continents, in both hemispheres. We sequenced the internal transcribed spacer (ITS) and partial large subunit (LSU) nuclear rDNA from 292 spore mats (visible mitospore clumps) collected in Argentina, Chile, China, Mexico and the USA between 2009 and 2012. We collated spore mat ITS sequences with 105 fruit body and 47 EcM root sequences to generate operational taxonomic units (OTUs). Phylogenetic inferences were made through analyses of both molecular data sets. A total of 48 OTUs from spore mats represented six independent EcM Pezizales lineages and included truffles and cup fungi. Three clades of seven OTUs have no known meiospore stage. Mitospores failed to germinate on sterile media, or form ectomycorrhizas on Quercus, Pinus and Populus seedlings, consistent with a hypothesized role of spermatia. The broad geographic range, high frequency and phylogenetic diversity of spore mats produced by EcM Pezizales suggests that a mitospore stage is important for many species in this group in terms of mating, reproduction and/or dispersal.     

Ancient geographical gaps and paleo-climate shape the phylogeography of an endemic bird in the sky islands of southern India

Background

Sky islands, formed by the highest reaches of mountain tracts physically isolated from one another, represent one of the biodiversity-rich regions of the world. Comparative studies of geographically isolated populations on such islands can provide valuable insights into the biogeography and evolution of species on these islands. The Western Ghats mountains of southern India form a sky island system, where the relationship between the island structure and the evolution of its species remains virtually unknown despite a few population genetic studies.

Methods and Principal Findings

We investigated how ancient geographic gaps and glacial cycles have partitioned genetic variation in modern populations of a threatened endemic bird, the White-bellied Shortwing Brachypteryx major, across the montane Shola forests on these islands and also inferred its evolutionary history. We used Bayesian and maximum likelihood-based phylogenetic and population-genetic analyses on data from three mitochondrial markers and one nuclear marker (totally 2594 bp) obtained from 33 White-bellied Shortwing individuals across five islands. Genetic differentiation between populations of the species correlated with the locations of deep valleys in the Western Ghats but not with geographical distance between these populations. All populations revealed demographic histories consistent with population founding and expansion during the Last Glacial Maximum. Given the level of genetic differentiation north and south of the Palghat Gap, we suggest that these populations be considered two different taxonomic species.

Conclusions and Significance

Our results show that the physiography and paleo-climate of this region historically resulted in multiple glacial refugia that may have subsequently driven the evolutionary history and current population structure of this bird. The first avian genetic study from this biodiversity hotspot, our results provide insights into processes that may have impacted the speciation and evolution of the endemic fauna of this region.     

Isolation by environment in White‐breasted Nuthatches (Sitta carolinensis) of the Madrean Archipelago sky islands: a landscape genomics approach

Understanding landscape processes driving patterns of population genetic differentiation and diversity has been a long‐standing focus of ecology and evolutionary biology. Gene flow may be reduced by historical, ecological or geographic factors, resulting in patterns of isolation by distance (IBD) or isolation by environment (IBE). Although IBE has been found in many natural systems, most studies investigating patterns of IBD and IBE in nature have used anonymous neutral genetic markers, precluding inference of selection mechanisms or identification of genes potentially under selection. Using landscape genomics, the simultaneous study of genomic and ecological landscapes, we investigated the processes driving population genetic patterns of White‐breasted Nuthatches (Sitta carolinensis) in sky islands (montane forest habitat islands) of the Madrean Archipelago. Using more than 4000 single nucleotide polymorphisms and multiple tests to investigate the relationship between genetic differentiation and geographic or ecological distance, we identified IBE, and a lack of IBD, among sky island populations of S. carolinensis. Using three tests to identify selection, we found 79 loci putatively under selection; of these, seven matched CDS regions in the Zebra Finch. The loci under selection were highly associated with climate extremes (maximum temperature of warmest month and minimum precipitation of driest month). These results provide evidence for IBE – disentangled from IBD – in sky island vertebrates and identify potential adaptive genetic variation.     

Soil spore banks of ectomycorrhizal fungi in endangered Japanese Douglas-fir forests

Interactions between trees and ectomycorrhizal fungi are critical to the growth and survival of both partners. However, ectomycorrhizal symbiosis has barely been explored in endangered trees, and no information is available regarding soil spore banks of ectomycorrhizal fungi from forests of threatened trees. Here, we evaluated soil spore banks of ectomycorrhizal fungi from endangered Japanese Douglas-fir (Pseudotsuga japonica) forests using bioassay approaches with congeneric P. menziesii and Pinus densiflora seedlings in combination with molecular identification techniques. Rhizopogon togasawariana was predominant in soil propagule banks and was found in all remaining P. japonica forests when assayed with P. menziesii, while no colonization of this fungus was observed on Pinus seedlings. Given the observed specificity of R. togasawariana for P. menziesii and its phylogenetic position within the Pseudotsuga-specific Rhizopogon lineage, its geographical distribution is likely restricted to the remaining Japanese Douglas-fir forests, indicating a high extinction risk for this fungus as well as its endangered host. Spore banks of R. togasawariana remained highly infective after preservation for 1 year or heat treatment at 70 °C, suggesting an ecological strategy of establishing ectomycorrhizal associations on regenerating Japanese Douglas-fir seedlings after disturbance, as observed in other Rhizopogon–Pinaceae combinations. Therefore, the regeneration of Japanese Douglas-fir seedlings may depend largely on the soil spore banks dominated by R. togasawariana, which has co-evolved with the Japanese Douglas-fir for over 30 million years. More attention must be paid to underground ectomycorrhizal fungi for the conservation of endangered tree species, especially in the era of human-induced mass extinction.     

Comparison of root-associated communities of native and non-native ectomycorrhizal hosts in an urban landscape

Non-native tree species are often used as ornamentals in urban landscapes. However, their root-associated fungal communities remain yet to be examined in detail. Here, we compared richness, diversity and community composition of ectomycorrhizosphere fungi in general and ectomycorrhizal (EcM) fungi in particular between a non-native Pinus nigra and a native Quercus macrocarpa across a growing season in urban parks using 454-pyrosequencing. Our data show that, while the ectomycorrhizosphere community richness and diversity did not differ between the two host, the EcM communities associated with the native host were often more species rich and included more exclusive members than those of the non-native hosts. In contrast, the ectomycorrhizosphere communities of the two hosts were compositionally clearly distinct in nonmetric multidimensional ordination analyses, whereas the EcM communities were only marginally so. Taken together, our data suggest EcM communities with broad host compatibilities and with a limited numbers of taxa with preference to the non-native host. Furthermore, many common fungi in the non-native Pinus were not EcM taxa, suggesting that the fungal communities of the non-native host may be enriched in non-mycorrhizal fungi at the cost of the EcM taxa. Finally, while our colonization estimates did not suggest a shortage in EcM inoculum for either host in urban parks, the differences in the fungi associated with the two hosts emphasize the importance of using native hosts in urban environments as a tool to conserve endemic fungal diversity and richness in man-made systems.     

Host phylogeny is the primary determinant of ectomycorrhizal fungal community composition in the permafrost ecosystem of eastern Siberia at a regional scale

Ectomycorrhizal (EM) fungal communities have been studied worldwide; however, those in the very cold and dry continental climate zone of northern Eurasia remain understudied. We investigated EM fungal community structure on plant roots and its determinants in eastern Siberia. We identified 291 EM fungal taxa belonging to 37 fungal genera from nine sites spanning 2100 km. In a variation partitioning analysis, host plant phylogeny was the primary factor that explained variation in fungal community composition, followed by spatial distance, soil, and climate. Host specificity and preference were attributed to differences in EM fungal community composition among host plants. The EM fungal community on Larix cajanderi, the dominant canopy tree in the region, was characterized by a high proportion of Suillus and Rhizopogon species. This implies that these specialist fungal symbionts have a close ecological relationship with pioneer Larix trees to adapt to the harsh continental climate of Siberia.     

A complex history of introgression and vicariance in a threatened montane skink (<Emphasis Type="Italic">Pseudemoia cryodroma</Emphasis>) across an Australian sky island system

Species endemic to sky island systems are isolated to mountain peaks and high elevation plateaux both geographically and ecologically, making them particularly vulnerable to the effects of climate change. Pressures associated with climate change have already been linked to local extinctions of montane species, emphasizing the importance of understanding the genetic diversity and population connectivity within sky islands systems for the conservation management of remaining populations. Our study focuses on the endangered alpine skink Pseudemoia cryodroma, which is endemic to the Victorian Alps in south-eastern Australia, and has a disjunct distribution in montane habitats above 1100 m a.s.l. Using mitochondrial DNA (mtDNA) and microsatellite loci, we investigated species delimitation, genetic connectivity and population genetic structure across the geographic range of this species. We found discordance between genetic markers, indicating historical mtDNA introgression at one of the study sites between P. cryodroma and the closely related, syntopic P. entrecasteauxii. Molecular diversity was positively associated with site elevation and extent of suitable habitat, with inbreeding detected in three of the five populations. These results demonstrate the complex interaction between geography and habitat in shaping the population structure and genetic diversity of P. cryodroma, and highlight the importance of minimising future habitat loss and fragmentation for the long-term persistence of this species.     

Flying between Sky Islands: The Effect of Naturally Fragmented Habitat on Butterfly Population Structure

High elevation montane areas are called “sky islands” when they occur as a series of high mountains separated by lowland valleys. Different climatic conditions at high elevations makes sky islands a specialized type of habitat, rendering them naturally fragmented compared to more continuous habitat at lower elevations. Species in sky islands face unsuitable climate in the intervening valleys when moving from one montane area to another. The high elevation shola-grassland mosaic in the Western Ghats of southern India form one such sky island complex. The fragmented patches make this area ideal to study the effect of the spatial orientation of suitable habitat patches on population genetic structure of species found in these areas. Past studies have suggested that sky islands tend to have genetically structured populations, possibly due to reduced gene flow between montane areas. To test this hypothesis, we adopted the comparative approach. Using Amplified Fragment Length Polymorphisms, we compared population genetic structures of two closely related, similar sized butterfly species: Heteropsis oculus, a high elevation shola-grassland specialist restricted to the southern Western Ghats, and Mycalesis patnia, found more continuously distributed in lower elevations. In all analyses, as per expectation the sky island specialist H. oculus exhibited a greater degree of population genetic structure than M. patnia, implying a difference in geneflow. This difference in geneflow in turn appears to be due to the natural fragmentation of the sky island complexes. Detailed analysis of a subset of H. oculus samples from one sky island complex (the Anamalais) showed a surprising genetic break. A possible reason for this break could be unsuitable conditions of higher temperature and lower rainfall in the intervening valley region. Thus, sky island species are not only restricted by lack of habitat continuity between montane areas, but also by the nature of the intervening habitat.     

Regional and local controls on postglacial vegetation and fire in the Siskiyou Mountains, northern California, USA

The Siskiyou Mountains of northwestern California and southwestern Oregon are a floristic hotspot, and the high diversity of conifers there likely results from a combination of geological, ecological, climatological and historical factors. To evaluate how past climate variability has influenced the composition, structure and fire regime of the Siskiyou forests, pollen, charcoal, and lithological evidence was examined from two lakes along a moisture gradient to reconstruct the vegetation, fire and climate history. The late-glacial period was characterized by subalpine parkland and infrequent fire at both sites. During the late-glacial/Early Holocene transition period, subalpine parkland was replaced by a closed forest of Pinus, Cupressaceae, Abies and Pseudotsuga and more frequent fires a 1000 years earlier at the wetter site, and it is likely that reduced Pacific Ocean upwelling created warmer drier conditions at the coast. In the Early Holocene, Pinus, Cupressaceae were less abundant and fire less frequent at the coastal site during a period of increased coastal upwelling and fog production. In the Late Holocene, Abies, Pseudotsuga, Pinus, and Quercus vaccinifolia increased in the forest at both sites suggesting a widespread response to cooling. Fewer fires at the wetter site may account for the abundance of Picea breweriana within the last 1000 years. The comparison of the two records implies that large-scale controls in climate during the last 14,000 cal yr BP have resulted in major changes in vegetation and fire regime. Asynchrony in the ecosystem response of wetter and drier sites arises from small-scale spatial variations in effective moisture and temperature resulting from topographically-influenced microclimates and coastal-to-inland climate gradients.     

Variability of ecological and autotrophic host specificity in a mycoheterotrophic system: Pterospora andromedea and associated fungal and conifer hosts

Few studies of tripartite mycoheterotrophic systems have examined ecological specificity across broad geographic ranges or addressed autotrophic host specificity. Pterospora andromedea was selected as an ideal candidate to examine ecological specificity of a mycoheterotrophic system as it is widely distributed, has been shown to have high levels of symbiont specificity with Rhizopogon subgenus Amylopogon, and is found with several autotrophic hosts. Pairs of P. andromedea + Rhizopogon spp. samples were co-collected across North America and were sequenced using trnL and ITS, respectively. Bayesian phylogenetic reconstructions between the co-collected taxa were used to examine ecological specificity, and for subsequent tests for autotroph specificity. P. andromedea lineages exhibited both high specificity and relaxed specificity for fungal symbionts and autotrophic hosts across the geographic landscape under allopatric and sympatric conditions. This strong evidence for geographic mosaics of specificity in mycoheterotrophic systems is an important future consideration in determining the evolutionary ecology of mycoheterotrophs.     

Genomic insight into “sky island” species diversification in a mountainous biodiversity hotspot

“Sky island” species diversification contributes greatly to mountainous biodiversity. However, the underlying genomic divergence and the inferred drivers remain largely unknown. In this study, we examined the diversification history of five diploid species with three exclusively endemic to the sky islands (mountains) of the Himalaya–Hengduan Mountains biodiversity hotspot. All of them together comprise a clade of the genus Eutrema (Brassicaceae). We resequenced genomes of multiple individuals of the found populations for each species. We recovered the inconsistent phylogenetic relationships between plastome and nuclear‐genome trees for one species. Based on nuclear population genomic data, we detected high genetic divergence between five species with limited gene flow. Four species seemed to diverge mainly through geographical isolation, whereas one arose through hybrid origin. The origins of the sampled five species were dated to within the late Miocene when mountains were uplifted and climates oscillated. All species decreased their population sizes since the inferred origin of each species initially, but only two of them expanded after the Last Glacial Maximum. Together, these findings suggest that geographic isolation plays an important role in driving the sky island species diversification of the sampled species in addition to the occasional gene flow that might have led to the hybrid origin of some sky island species, similar to the species diversification of sea islands.     

The population biology of Bromus tectorum in forests: distinguishing the opportunity for dispersal from environmental restriction

Summary With increasing elevation and corresponding changes in the macroclimate, forest zones in the Intermountain Region of western North America are often dominated in turn by Pinus ponderosa, Pseudotsuga menziesii, Abies grandis, an Thuja plicata. Bromus tectorum (cheatgrass), and introduced annual grass now abundant in the Region's steppe, is uncommon in mature stands representative of these forest zones. In order to determine whether B. tectorum is largely excluded from these forests by insufficient seed dispersal or environmental restriction(s), the grass's demography was compared in each of four years among populations experimentally-introduced into mature forests. The number of recruits did not differ among the Pinus, Pseudotsuga, and Abies sites; recruitment was however significantly lower on the coolest site dominated by Thuja. Emergence in both the low elevation Pinus and Pseudotsuga sites was about the same in autumn, winter, and spring. In the cooler, moister Abies and Thuja sites, emergence was limited to autumn and early winter. Survival in these forest sites ranged between 0 and 87%. The percentage of the total population to survive until harvest was highest in the Pseudotsuga site, intermediate in the Pinus and Abies sites, and lowest in the Thuja site. Compared with B. tectorum in the steppe, the surviving plants were small, and few produced seeds. All parents were members of either the autumn or winter cohort, and most parents produced only one seed. No seeds were produced at the Thuja site. Although phenotypic plasticity apparently contributes to the wide ecological amplitude of this grass, its growing season on these sites in most years is too short for it to reproduce. Consequently, these forest zones broadly define the current environmental limits to the distribution of cheatgrass in this portion of its new geographic range.     

Co-invasive exotic pines and their ectomycorrhizal symbionts show capabilities for wide distance and altitudinal range expansion

We asked if exotic Pinus elliotti seedlings can survive and form ectomycorrhizas at higher elevations and long distances from their current range, and which ECM partners disperse to these soils. We selected three plots at four grassland sites along an altitudinal gradient (900, 1600, 2200, and 2700 m asl) established at c. 110, 3000, 6000, and 9000 m from the closest pine plantation, respectively. We combined field experiments with glasshouse assays to assess survival and ECM fungi in roots and soils. A pine plantation close to the lowest site was also selected for DNA metabarcoding of soils. Pine seedlings survived at all altitudes but not all formed mycorrhizas. They formed mycorrhizas with Suillus granulatus at 900, 1600, and 2200 m asl (i.e. up to 6000 m from the closest pine plantation), and with Rhizopogon pseudoroseolus and Thelephora terrestris at lower altitudes and distances. Twelve ECM fungal OTUs were found in grasslands and 34 were detected in the pine plantation. Although richness and abundance of ECM fungi decreased with increasing distance from the pine plantation, there was at least one non-native ECM fungal species present in each sampling site, even at 2700 masl and 9000 m distance from the closest plantation. This study provides evidence that the availability of suitable fungal symbionts might constrain but not hinder the expansion of a pine species over wide distances and altitudinal zones even in areas with no native ECM fungi.