In orchid species that have populations occurring in strongly contrasting habitats, mycorrhizal divergence and other habitat-specific adaptations may lead to the formation of reproductively isolated taxa and ultimately to species formation. However, little is known about the mycorrhizal communities associated with recently diverged sister taxa that occupy different habitats.In this study, 454 amplicon pyrosequencing was used to investigate mycorrhizal communities associating with Epipactis helleborine in its typical forest habitat and with its presumed sister species E. neerlandica that almost exclusively occurs in coastal dune habitats. Samples of the phylogenetically more distant E. palustris, which co-occurred with E. neerlandica, were also included to investigate the role of habitat-specific conditions on mycorrhizal communities.A total of 105 operational taxonomic units (OTUs) of putative orchid mycorrhizal fungi were observed in the three studied species. The majority of these fungi were endophytic fungi of Helotiales and ectomycorrhizal fungi belonging to Thelephoraceae, Sebacinaceae and Inocybaceae. In addition, a large number of other ectomycorrhizal taxa were detected, including Cortinarius, Cenococcum, Tuber, Geopora, Wilcoxina, Meliniomyces, Hebeloma, Tricholoma, Russula and Peziza Mycorrhizal communities differed significantly between the three species, but differences were most pronounced between the forest species (E. helleborine) and the two dune slack species (E. neerlandica and E. palustris).The results clearly showed that recently diverged orchid species that occupy different habitats were characterized by significantly different mycorrhizal communities and call for more detailed experiments that aim at elucidating the contribution of habitat-specific adaptations in general and mycorrhizal divergence in particular to the process of speciation in orchids.