Hybrid speciation driven by multilocus introgression of ecological traits.

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

Contrasting signatures of genomic divergence during sympatric speciation.

The transition from 'well-marked varieties' of a single species into 'well-defined species'-especially in the absence of geographic barriers to gene flow (sympatric speciation)-has puzzled evolutionary biologists ever since Darwin1,2. Gene flow counteracts the buildup of genome-wide differentiation, which is a hallmark of speciation and increases the likelihood of the evolution of irreversible reproductive barriers (incompatibilities) that complete the speciation process3. Theory predicts that the genetic architecture of divergently selected traits can influence whether sympatric speciation occurs4, but empirical tests of this theory are scant because comprehensive data are difficult to collect and synthesize across species, owing to their unique biologies and evolutionary histories5. Here, within a young species complex of neotropical cichlid fishes (Amphilophus spp.), we analysed genomic divergence among populations and species. By generating a new genome assembly and re-sequencing 453 genomes, we uncovered the genetic architecture of traits that have been suggested to be important for divergence. Species that differ in monogenic or oligogenic traits that affect ecological performance and/or mate choice show remarkably localized genomic differentiation. By contrast, differentiation among species that have diverged in polygenic traits is genomically widespread and much higher overall, consistent with the evolution of effective and stable genome-wide barriers to gene flow. Thus, we conclude that simple trait architectures are not always as conducive to speciation with gene flow as previously suggested, whereas polygenic architectures can promote rapid and stable speciation in sympatry.

Convergence in sympatric swallowtail butterflies reveals ecological interactions as a key driver of worldwide trait diversification.

Ecological interactions can promote phenotypic diversification in sympatric species. While competition can enhance trait divergence, other ecological interactions may promote convergence in sympatric species. Within butterflies, evolutionary convergences in wing color patterns have been reported between distantly related species, especially in females of palatable species, where mimetic color patterns are promoted by predator communities shared with defended species living in sympatry. Wing color patterns are also often involved in species recognition in butterflies, and divergence in this trait has been reported in closely related species living in sympatry as a result of reproductive character displacement. Here, we investigate the effect of sympatry between species on the convergence vs. divergence of their wing color patterns in relation to phylogenetic distance, focusing on the iconic swallowtail butterflies (family Papilionidae). We developed an unsupervised machine learning-based method to estimate phenotypic distances among wing color patterns of 337 species, enabling us to finely quantify morphological diversity at the global scale among species and allowing us to compute pairwise phenotypic distances between sympatric and allopatric species pairs. We found phenotypic convergence in sympatry, stronger among distantly related species, while divergence was weaker and restricted to closely related males. The convergence was stronger among females than males, suggesting that differential selective pressures acting on the two sexes drove sexual dimorphism. Our results highlight the significant effect of ecological interactions driven by predation pressures on trait diversification in Papilionidae and provide evidence for the interaction between phylogenetic proximity and ecological interactions in sympatry, acting on macroevolutionary patterns of phenotypic diversification.

Sympatric speciation of the spiny mouse from Evolution Canyon in Israel substantiated genomically and methylomically.

SignificanceWhether sympatric speciation (SS) is rare or common is still debated. Two populations of the spiny mouse, Acomys cahirinus, from Evolution Canyon I (EC I) in Israel have been depicted earlier as speciating sympatrically by molecular markers and transcriptome. Here, we investigated SS both genomically and methylomically, demonstrating that the opposite populations of spiny mice are sister taxa and split from the common ancestor around 20,000 years ago without an allopatric history. Mate choice, olfactory receptors, and speciation genes contributed to prezygotic/postzygotic reproductive isolation. The two populations showed different methylation patterns, facilitating adaptation to their local environment. They cope with abiotic and biotic stresses, due to high solar interslope radiation differences. We conclude that our new genomic and methylomic data substantiated SS.

A genetic switch for male UV iridescence in an incipient species pair of sulphur butterflies.

Mating cues evolve rapidly and can contribute to species formation and maintenance. However, little is known about how sexual signals diverge and how this variation integrates with other barrier loci to shape the genomic landscape of reproductive isolation. Here, we elucidate the genetic basis of ultraviolet (UV) iridescence, a courtship signal that differentiates the males of Colias eurytheme butterflies from a sister species, allowing females to avoid costly heterospecific matings. Anthropogenic range expansion of the two incipient species established a large zone of secondary contact across the eastern United States with strong signatures of genomic admixtures spanning all autosomes. In contrast, Z chromosomes are highly differentiated between the two species, supporting a disproportionate role of sex chromosomes in speciation known as the large-X (or large-Z) effect. Within this chromosome-wide reproductive barrier, linkage mapping indicates that cis-regulatory variation of bric a brac (bab) underlies the male UV-iridescence polymorphism between the two species. Bab is expressed in all non-UV scales, and butterflies of either species or sex acquire widespread ectopic iridescence following its CRISPR knockout, demonstrating that Bab functions as a suppressor of UV-scale differentiation that potentiates mating cue divergence. These results highlight how a genetic switch can regulate a premating signal and integrate with other reproductive barriers during intermediate phases of speciation.

The Evolution of Local Co-occurrence in Birds in Relation to Latitude, Degree of Sympatry, and Range Symmetry.

AbstractRecent speciation rates and the degree of range-wide sympatry are usually higher farther from the equator. Is there also a higher degree of secondary syntopy (coexistence in local assemblages in sympatry) at higher latitudes and, subsequently, an increase in local species richness? We studied the evolution of syntopy in passerine birds using worldwide species distribution data. We chose recently diverged species pairs from subclades not older than 5 or 7 million years, range-wide degree of sympatry not lower than 5% or 25%, and three definitions of the breeding season. We related their syntopy to latitude, the degree of sympatry (breeding range overlap), range symmetry, and the age of split. Syntopy was positively related to latitude, but it did not differ between tropical and temperate regions, instead increasing from the Southern to the Northern Hemisphere. Syntopy was also higher in species pairs with a higher degree of sympatry and more symmetric ranges, but it did not predict local species richness. Following speciation, species in the Northern Hemisphere presumably achieve positive local co-occurrence faster than elsewhere, which could facilitate their higher speciation rates. However, this does not seem to be linked to local species richness, which is probably governed by other processes.

Competitive Displacement and Agonistic Character Displacement, or the Ghost of Interference Competition.

AbstractInterference competition can drive species apart in habitat use through competitive displacement in ecological time and agonistic character displacement (ACD) over evolutionary time. As predicted by ACD theory, sympatric species of rubyspot damselflies (Hetaerina spp.) that respond more aggressively to each other in staged encounters differ more in microhabitat use. However, the same pattern could arise from competitive displacement if dominant species actively exclude subordinate species from preferred microhabitats. The degree to which habitat partitioning is caused by competitive displacement can be assessed with removal experiments. We carried out removal experiments with three species pairs of rubyspot damselflies. With competitive displacement, removing dominant species should allow subordinate species to shift into the dominant species' microhabitat. Instead, we found that species-specific microhabitat use persisted after the experimental removals. Thus, the previously documented association between heterospecific aggression and microhabitat partitioning in this genus is most likely a product of divergence in habitat preferences caused by interference competition in the evolutionary past.

Diet in phenotypically divergent sympatric species of African weakly electric fish (genus: Campylomormyrus)-A hybrid capture/HTS metabarcoding approach.

Ecological speciation within the mormyrid genus Campylomormyrus resulted in sympatric species exhibiting divergence in their feeding apparatus and electric organ discharge (EOD). This study documents the overall diet of the genus Campylomormyrus and examines the hypothesis that the Campylomormyrus radiation is caused by adaptation to different food sources. We performed diet assessment of five sympatric Campylomormyrus species (C. alces, C. compressirostris, C. curvirostris, C. tshokwe, C. numenius) and their sister taxon Gnathonemus petersii with markedly different snout morphologies and EODs using hybrid capture/HTS DNA metabarcoding of their stomach contents. Our approach allowed for high taxonomic resolution of prey items, including benthic invertebrates, allochthonous invertebrates and vegetation. Comparisons of the diet compositions using quantitative measures and diet overlap indices revealed that all species are able to exploit multiple food niches in their habitats, that is fauna at the bottom, the water surface and the water column. A major part of the diet is larvae of aquatic insects, such as dipterans, coleopterans and trichopterans, known to occur in holes and interstitial spaces of the substrate. The results indicate that different snout morphologies and the associated divergence in the EOD could translate into different prey spectra. This suggests that the diversification in EOD and/or morphology of the feeding apparatus could be under functional adaptation.

Recurrent gene flow events occurred during the diversification of clownfishes of the skunk complex.

Clownfish (subfamily Amphiprioninae) are an iconic group of coral reef fish that evolved a mutualistic interaction with sea anemones, which triggered the adaptive radiation of the clade. Within clownfishes, the "skunk complex" is particularly interesting. Besides ecological speciation, interspecific gene flow and hybrid speciation are thought to have shaped the evolution of the group. We investigated the mechanisms characterizing the diversification of this complex. By taking advantage of their disjunct geographical distribution, we obtained whole-genome data of sympatric and allopatric populations of the three main species of the complex (Amphiprion akallopisos, A. perideraion and A. sandaracinos). We examined population structure, genomic divergence and introgression signals and performed demographic modelling to identify the most realistic diversification scenario. We excluded scenarios of strict isolation or hybrid origin of A. sandaracinos. We discovered moderate gene flow from A. perideraion to the ancestor of A. akallopisos + A. sandaracinos and weak gene flow between the species in the Indo-Australian Archipelago throughout the diversification of the group. We identified introgressed regions in A. sandaracinos and detected in A. perideraion two large regions of high divergence from the two other species. While we found that gene flow has occurred throughout the species' diversification, we also observed that recent admixture was less pervasive than initially thought, suggesting a role of host repartition or behavioural barriers in maintaining the genetic identity of the species in sympatry.

Phylogenomics of the narrowly endemic Eurycheilichthys (Siluriformes: Loricariidae): Sympatric species with non-sister relationships suggest mainly allopatric speciation.

Armored catfishes of the genus Eurycheilichthys are endemic to Southern Brazil and Misiones (Argentina) comprising nine species of small size, with a high degree of sympatry and species diversity distributed in two river basins. Here we use new genome-wide data to infer a species phylogeny and test species boundaries for this poorly known group. We estimate 1) the phylogenetic relationships of the species of Eurycheilichthys based on 29,350 loci in 65 individuals of nine species plus outgroups, and 2) the population structure and differentiation based on 43,712 loci and 62 individuals to estimate how geography may have acted on speciation and formation of the sympatric species groups. Analyses support the monophyly of the genus and suggest two species-inclusive clades (East and West) with high support and very recently diverged species. Western clade contains E. limulus (from upper Jacuí River basin) that is sister to Western species of the Taquari-Antas basin plus E. paucidens. The Eastern clade contains E. pantherinus (from Uruguay River basin) sister to the Eastern species of the Taquari-Antas basin E. coryphaenus, plus the central-distributed species E. planus and E. vacariensis, and the more widely-distributed species E. luisae. Eurycheilichthys luisae is not monophyletic and may contain one or more cryptic species or hybrid individuals. A stronger diversity on structure of lineages on the Taquari-Antas, when compared to upper Uruguay and Jacuí River basins, and the fact that most of the sympatrically distributed taxa have non-sister relationships suggest a scenario of mainly allopatric speciation and may indicate a more dynamic landscape with headwater capture events among these tributaries.

A consequential one-night stand: Episodic historical hybridization leads to mitochondrial takeover in sympatric desert ant-eating spiders.

Disentangling the genomic intricacies underlying speciation and the causes of discordance between sources of evidence can offer remarkable insights into evolutionary dynamics. The ant-eating spider Zodarion nitidum, found across the Middle East and Egypt, displays yellowish and blackish morphs that co-occur sympatrically. These morphs additionally differ in behavioral and physiological features and show complete pre-mating reproductive isolation. In contrast, they possess similar sexual features and lack distinct differences in their mitochondrial DNA. We analyzed both Z. nitidum morphs and outgroups using genome-wide and additional mitochondrial DNA data. The genomic evidence indicated that Yellow and Black are reciprocally independent lineages without signs of recent admixture. Interestingly, the sister group of Yellow is not Black but Z. luctuosum, a morphologically distinct species. Genomic gene flow analyses pinpointed an asymmetric nuclear introgression event, with Yellow contributing nearly 5 % of its genome to Black roughly 320,000 years ago, intriguingly aligning with the independently estimated origin of the mitochondrial DNA of Black. We conclude that the blackish and yellowish morphs of Z. nitidum are long-diverged distinct species, and that the ancient and modest genomic introgression event registered resulted in a complete mitochondrial takeover of Black by Yellow. This investigation underscores the profound long-term effects that even modest hybridization events can have on the genome of organisms. It also exemplifies the utility of phylogenetic networks for estimating historical events and how integrating independent lines of evidence can increase the reliability of such estimations.

Divergent warning patterns influence male and female mating behaviours in a tropical butterfly.

Traits under divergent ecological selection that also function during mating can be important in maintaining species boundaries. Few studies have considered mutual mate choice, where both males and females base mating decisions on the same trait. Wing colouration in Heliconius butterflies evolved as a warning signal but also functions as a mating cue. We investigated the contribution of visual preference to assortative mating in an aposematic butterfly Heliconius cydno in the context of reproductive isolation with its sympatric, visually distinct relative Heliconius melpomene. Heliconius cydno have conspicuous white bands on their forewings, whereas those of H. melpomene are red in colour. We predicted that both sexes of H. cydno contributed to assortative mating by exhibiting visual preference towards conspecific wing colouration. We analysed published and new data from preference experiments, in which males were presented with conspecific and H. melpomene females. We also recorded female responses and mating outcomes in choice experiments, involving conspecific males with either the original white or artificially painted red forewing bands. Both sexes of H. cydno responded more positively towards the conspecific colouration, and males strongly preferred females of its own colours. In contrast, male colouration did not predict mating outcomes in female choice experiments. As courtships are initiated by males in butterflies, our findings suggest that female visual preference might be of secondary importance in H. cydno. Our data also suggest that the contribution of visual preference to reproductive isolation might be unequal between H. cydno and its sympatric relative H. melpomene.

Sympatry and parapatry among rocky reef cichlids of Lake Victoria explained by female mating preferences.

Work on the Lake Victoria cichlids Pundamilia nyererei (red dorsum males, deeper water), Pundamilia pundamilia (blue males, shallower water) and related species pairs has provided insights into processes of speciation. Here, we investigate the female mating behaviour of 5 Pundamilia species and 4 of their F1 hybrids through mate choice trials and paternity testing. Complete assortative mating was observed among all sympatric species. Parapatric species with similar depth habitat distributions interbred whereas other parapatric and allopatric species showed complete assortative mating. F1 hybrids mated exclusively with species accepted by females of the parental species. The existence of complete assortative mating among some currently allopatric species suggests that pre-existing mating barriers could be sufficient to explain current patterns of co-existence, although, of course, many other factors may be involved. Regardless of the mechanism, mating preferences may influence species distribution in potentially hybridizing taxa, such as in the adaptive radiation of cichlid fish. We suggest that this at least partly explains why some species fail to establish breeding populations in locations where they are occasionally recorded. Our results support the notion that the mating preferences of potentially cross-breeding species ought to be included in coexistence theory.

Tempo and mode of morphological evolution are decoupled from latitude in birds.

The latitudinal diversity gradient is one of the most striking patterns in nature, yet its implications for morphological evolution are poorly understood. In particular, it has been proposed that an increased intensity of species interactions in tropical biota may either promote or constrain trait evolution, but which of these outcomes predominates remains uncertain. Here, we develop tools for fitting phylogenetic models of phenotypic evolution in which the impact of species interactions-namely, competition-can vary across lineages. Deploying these models on a global avian trait dataset to explore differences in trait divergence between tropical and temperate lineages, we find that the effect of latitude on the mode and tempo of morphological evolution is weak and clade- or trait dependent. Our results indicate that species interactions do not disproportionately impact morphological evolution in tropical bird families and question the validity of previously reported patterns of slower trait evolution in the tropics.

Is self-incompatibility a reproductive barrier for hybridization in a sympatric species?

PREMISE: Barriers at different reproductive stages contribute to reproductive isolation. Self-incompatibility (SI) systems that prevent self-pollination could also act to control interspecific pollination and contribute to reproductive isolation, preventing hybridization. Here we evaluated whether SI contributes to reproductive isolation among four co-occurring Opuntia species that flower at similar times and may hybridize with each other. METHODS: We assessed whether Opuntia cantabrigiensis, O. robusta, O. streptacantha, and O. tomentosa, were self-compatible and formed hybrid seeds in five manipulation treatments to achieve self-pollination, intraspecific cross-pollination, open pollination (control), interspecific crosses or apomixis, then recorded flowering phenology and synchrony. RESULTS: All species flowered in the spring with a degree of synchrony, so that two pairs of species were predisposed to interspecific pollination (O. cantabrigiensis with O. robusta, O. streptacantha with O. tomentosa). All species had distinct reproductive systems: Opuntia cantabrigiensis is self-incompatible and did not produce hybrid seeds as an interspecific pollen recipient; O. robusta is a dioecious species, which formed a low proportion of hybrid seeds; O. streptacantha and O. tomentosa are self-compatible and produced hybrid seeds. CONCLUSIONS: Opuntia cantabrigiensis had a strong pollen-pistil barrier, likely due to its self-incompatibility. Opuntia robusta, the dioecious species, is an obligate outcrosser and probably partially lost its ability to prevent interspecific pollen germination. Given that the self-compatible species can set hybrid seeds, we conclude that pollen-pistil interaction and high flowering synchrony represent weak barriers; whether reproductive isolation occurs later in their life cycle (e.g., germination or seedling survival) needs to be determined.

Central African dwarf crocodiles found in syntopy are comparably divergent to South American dwarf caimans.

Recent molecular taxonomic advancements have expanded our understanding of crocodylian diversity, revealing the existence of previously overlooked species, including the Congo dwarf crocodile (Osteolaemus osborni) in the central Congo Basin rainforests. This study explores the genomic divergence between O. osborni and its better-known relative, the true dwarf crocodile (Osteolaemus tetraspis), shedding light on their evolutionary history. Field research conducted in the northwestern Republic of the Congo uncovered a locality where both species coexist in sympatry/syntopy. Genomic analysis of sympatric individuals reveals a level of divergence comparable to that between ecologically similar South American dwarf caimans (Paleosuchus palpebrosus and Paleosuchus trigonatus), suggesting parallel speciation in the Afrotropics and Neotropics during the Middle to Late Miocene, 10-12 Ma. Comparison of the sympatric and allopatric dwarf crocodiles indicates no gene flow between the analysed sympatric individuals of O. osborni and O. tetraspis. However, a larger sample will be required to answer the question of whether or to what extent these species hybridize. This study emphasizes the need for further research on the biology and conservation status of the Congo dwarf crocodile, highlighting its significance in the unique biodiversity of the Congolian rainforests and thus its potential as a flagship species.

Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe.

BACKGROUND: Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry. METHODS AND RESULTS: MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank. CONCLUSIONS: Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

Niche partitioning and individual specialisation in resources and space use of sympatric fur seals at their range margin.

Ecological theory predicts niche partitioning between high-level predators living in sympatry as a mechanism to minimise the selective pressure of competition. Accordingly, male Australian fur seals Arctocephalus pusillus doriferus and New Zealand fur seals A. forsteri that live in sympatry should exhibit partitioning in their broad niches (in habitat and trophic dimensions) in order to coexist. However, at the northern end of their distributions in Australia, both are recolonising their historic range after a long absence due to over-exploitation, and their small population sizes suggest competition should be weak and may allow overlap in niche space. We found some niche overlap, yet clear partitioning in diet trophic level (δ15N values from vibrissae), spatial niche space (horizontal and vertical telemetry data) and circadian activity patterns (timing of dives) between males of each species, suggesting competition may remain an active driver of niche partitioning amongst individuals even in small, peripheral populations. Consistent with individual specialisation theory, broad niches of populations were associated with high levels of individual specialisation for both species, despite putative low competition. Specialists in isotopic space were not necessarily specialists in spatial niche space, further emphasising their diverse individual strategies for niche partitioning. Males of each species displayed distinct foraging modes, with Australian fur seals primarily benthic and New Zealand fur seals primarily epipelagic, though unexpectedly high individual specialisation for New Zealand fur seals might suggest marginal populations provide exceptions to the pattern generally observed amongst other fur seals.

Character displacement drives trait divergence in a continental fauna.

Coexisting (sympatric) pairs of closely related species are often characterized by exaggerated trait differences. This widespread pattern is consistent with adaptation for reduced similarity due to costly interactions (i.e., "character displacement")-a classic hypothesis in evolutionary theory. But it is equally consistent with a community assembly bias in which lineages with greater trait differences are more likely to establish overlapping ranges in the first place (i.e., "species sorting"), as well as with null expectations of trait divergence through time. Few comparative analyses have explicitly modeled these alternatives, and it remains unclear whether trait divergence is a general prerequisite for sympatry or a consequence of interactions between sympatric species. Here, we develop statistical models that allow us to distinguish the signature of these processes based on patterns of trait divergence in closely related lineage pairs. We compare support for each model using a dataset of bill shape differences in 207 pairs of New World terrestrial birds representing 30 avian families. We find that character displacement models are overwhelmingly supported over species sorting and null expectations, indicating that exaggerated bill shape differences in sympatric pairs result from enhanced divergent selection in sympatry. We additionally detect a latitudinal gradient in character displacement, which appears strongest in the tropics. Our analysis implicates costly species interactions as powerful drivers of trait divergence in a major vertebrate fauna. These results help substantiate a long-standing but equivocally supported linchpin of evolutionary theory.

Exploring mycorrhizal diversity in sympatric mycoheterotrophic plants: a comparative study of Monotropastrum humile var. humile and M. humile var. glaberrimum.

Mycoheterotrophic plants (MHPs) rely on their mycorrhizal fungus for carbon and nutrient supply, thus a shift in mycobionts may play a crucial role in speciation. This study aims to explore the mycorrhizal diversity of two closely related and sympatric fully MHPs, Monotropastrum humile var. humile (Mhh) and M. humile var. glaberrimum (Mhg), and determine their mycorrhizal associations. A total of 1,108,710 and 1,119,071 ectomycorrhizal fungal reads were obtained from 31 Mhh and 31 Mhg, and these were finally assigned to 227 and 202 operational taxonomic units, respectively. Results show that sympatric Mhh and Mhg are predominantly associated with different fungal genera in Russulaceae. Mhh is consistently associated with members of Russula, whereas Mhg is associated with members of Lactarius. Associating with different mycobionts and limited sharing of fungal partners might reduce the competition and contribute to their coexistence. The ectomycorrhizal fungal communities are significantly different among the five forests in both Mhh and Mhg. The distinct mycorrhizal specificity between Mhh and Mhg suggests the possibility of different mycobionts triggered ecological speciation between sympatric species.