Species delimitation of common reef corals in the genus Pocillopora using nucleotide sequence phylogenies, population genetics and symbiosis ecology

Stony corals in the genus Pocillopora are among the most common and widely distributed of Indo-Pacific corals and, as such, are often the subject of physiological and ecological research. In the far Tropical Eastern Pacific (TEP), they are major constituents of shallow coral communities, exhibiting considerable variability in colony shape and branch morphology and marked differences in response to thermal stress. Numerous intermediates occur between morphospecies that may relate to extensive hybridization. The diversity of the Pocillopora genus in the TEP was analysed genetically using nuclear ribosomal (ITS2) and mitochondrial (ORF) sequences, and population genetic markers (seven microsatellite loci). The resident dinoflagellate endosymbiont (Symbiodinium sp.) in each sample was also characterized using sequences of the internal transcribed spacer 2 (ITS2) rDNA and the noncoding region of the chloroplast psbA minicircle. From these analyses, three symbiotically distinct, reproductively isolated, nonhybridizing, evolutionarily divergent animal lineages were identified. Designated types 1, 2 and 3, these groupings were incongruent with traditional morphospecies classification. Type 1 was abundant and widespread throughout the TEP; type 2 was restricted to the Clipperton Atoll; and type 3 was found only in Panama and the Galapagos Islands. Each type harboured a different Symbiodinium'species lineage' in Clade C, and only type 1 associated with the 'stress-tolerant'Symbiodinium glynni (D1). The accurate delineation of species and implementation of a proper taxonomy may profoundly improve our assessment of Pocillopora's reproductive biology, biogeographic distributions, and resilience to climate warming, information that must be considered when planning for the conservation of reef corals.     

Population genetics of an ecosystem-defining reef coral Pocillopora damicornis in the Tropical Eastern Pacific

Background

Coral reefs in the Tropical Eastern Pacific (TEP) are amongst the most peripheral and geographically isolated in the world. This isolation has shaped the biology of TEP organisms and lead to the formation of numerous endemic species. For example, the coral Pocillopora damicornis is a minor reef-builder elsewhere in the Indo-West Pacific, but is the dominant reef-building coral in the TEP, where it forms large, mono-specific stands, covering many hectares of reef. Moreover, TEP P. damicornis reproduces by broadcast spawning, while it broods mostly parthenogenetic larvae throughout the rest of the Indo-West Pacific. Population genetic surveys for P. damicornis from across its Indo-Pacific range indicate that gene flow (i.e. larval dispersal) is generally limited over hundreds of kilometers or less. Little is known about the population genetic structure and the dispersal potential of P. damicornis in the TEP.

Methodology

Using multilocus microsatellite data, we analyzed the population structure of TEP P. damicornis among and within nine reefs and test for significant genetic structure across three geographically and ecologically distinct regions in Panama.

Principal Findings/Conclusions

We detected significant levels of population genetic structure (global R_ST = 0.162), indicating restricted gene flow (i.e. larvae dispersal), both among the three regions (R_RT = 0.081) as well as within regions (R_SR = 0.089). Limited gene flow across a distinct environmental cline, like the regional upwelling gradient in Panama, indicates a significant potential for differential adaptation and population differentiation. Individual reefs were characterized by unexpectedly high genet diversity (avg. 94%), relatively high inbreeding coefficients (global F_IS = 0.183), and localized spatial genetic structure among individuals (i.e. unique genets) over 10 m intervals. These findings suggest that gene flow is limited in TEP P. damicornis populations, particularly among regions, but even over meter scales within populations.     

Tectonic survey of the Central Pacific coast of Costa Rica

Three fossil coral samples (Pocillopora elegans and Pocillopora damicornis) were collected from a beachrock in Puerto Escondido, central Pacific coast of Costa Rica. The age of the samples ranged from 1272+/-60 to 2010+/-60 years BP. They represent some of the few Holocene fossil corals known from the Pacific coast of the Americas. The uplift rate for the central Pacific coast of Costa Rica was calculated in 0.4 to 2.1 mm year(-1).     

Coral-associated bacterial communities on Ningaloo Reef, Western Australia

Coral-associated microbial communities from three coral species (Pocillopora damicornis, Acropora tenuis and Favites abdita) were examined every 3 months (January, March, June, October) over a period of 1 year on Ningaloo Reef, Western Australia. Tissue from corals was collected throughout the year and additional sampling of coral mucus and seawater samples was performed in January. Tissue samples were also obtained in October from P. damicornis coral colonies on Rottnest Island off Perth, 1200 km south of Ningaloo Reef, to provide comparisons between coral-microbial associates in different locations. The community structures of the coral-associated microorganisms were analysed using phylogenetic analysis of 16S rRNA gene clone libraries, which demonstrated highly diverse microbial profiles among all the coral species sampled. Principal component analysis revealed that samples grouped according to time and not species, indicating that coral-microbial associations may be a result of environmental drivers such as oceanographic characteristics, benthic community structure and temperature. Tissue samples from P. damicornis at Rottnest Island revealed similarities in bacteria to the samples at Ningaloo Reef. This study highlights that coral-associated microbial communities are highly diverse; however, the complex interactions that determine the stability of these associations are not necessarily dependent on coral host specificity.     

Spermatozoan ultrastructure of scleractinian corals from the eastern Pacific: Pocilloporidae and Agariciidae

Spermatozoa ofPocillopora damicornis, Pocillopora elegans (Astrocoeniina, Pocilloporidae) andPavona gigantea (Fungiina, Agariciidae) from the eastern Pacific (Isla del Caño, Costa Rica) were examined using transmission electron microscopy. The hermaphroditic pocilloporidsP. damicornis andP. elegans are spermiomorphologically very similar to hermaphroditic acroporids, being characterized by bullet-shaped nuclei and elongated mitochondria. Such traits have not been found in other families. Thus, the suborder Astrocoeniina, including pocilloporids and acroporids, can clearly be distinguished from other scleractinian suborders. This separation underlines the isolated position of the Astrocoeniina within the order Scleractinia following the evolutionary scheme of Wells. A conical sperm type, known from gonochoric species (sexes separate) from all families except Acroporidae and Pocilloporidae, was found in the gonochoric agariciidPavona gigantea. This supports previous findings that gonochoric corals share a unique and common sperm structure regardless of which family they belong to. However, no gonochoric Astrocoeniina have ever been examined. Hence, the question whether the sperm type common to gonochorists is also represented in Astrocoeniina, which would undermine the ultrastructural distinction of Astrocoeniina and other suborders seen among hermaphrodites, as well as the systematic value of sperm structure within scleractinian corals, remains open.     

Temperature-regulated bleaching and lysis of the coral Pocillopora damicornis by the novel pathogen Vibrio coralliilyticus

Coral bleaching is the disruption of symbioses between coral animals and their photosynthetic microalgal endosymbionts (zooxanthellae). It has been suggested that large-scale bleaching episodes are linked to global warming. The data presented here demonstrate that Vibrio coralliilyticus is an etiological agent of bleaching of the coral Pocillopora damicornis. This bacterium was present at high levels in bleached P. damicornis but absent from healthy corals. The bacterium was isolated in pure culture, characterized microbiologically, and shown to cause bleaching when it was inoculated onto healthy corals at 25 degrees C. The pathogen was reisolated from the diseased tissues of the infected corals. The zooxanthella concentration in the bacterium-bleached corals was less than 12% of the zooxanthella concentration in healthy corals. When P. damicornis was infected with V. coralliilyticus at higher temperatures (27 and 29 degrees C), the corals lysed within 2 weeks, indicating that the seawater temperature is a critical environmental parameter in determining the outcome of infection. A large increase in the level of the extracellular protease activity of V. coralliilyticus occurred at the same temperature range (24 to 28 degrees C) as the transition from bleaching to lysis of the corals. We suggest that bleaching of P. damicornis results from an attack on the algae, whereas bacterium-induced lysis and death are promoted by bacterial extracellular proteases. The data presented here support the bacterial hypothesis of coral bleaching.     

Ten polymorphic STR loci in the cosmopolitan reef coral, Pocillopora damicornis

We report the development of 10 polymorphic molecular markers containing short tandem repeats in the cosmopolitan reef-building coral, Pocillopora damicornis, an important model species for coral health, physiology, ecology, and genetics. The availability of polymorphic DNA markers in P. damicornis can act as impetus for investigations into inheritance and population genetics, as well as novel investigations into host-symbiont ecology and evolution. Coral bleaching and gene flow studies performed with these markers can have direct conservation implications.     

Restricted gene flow and local adaptation highlight the vulnerability of high‐latitude reefs to rapid environmental change

Global climate change poses a serious threat to the future health of coral reef ecosystems. This calls for management strategies that are focused on maximizing the evolutionary potential of coral reefs. Fundamental to this is an accurate understanding of the spatial genetic structure in dominant reef‐building coral species. In this study, we apply a genotyping‐by‐sequencing approach to investigate genome‐wide patterns of genetic diversity, gene flow, and local adaptation in a reef‐building coral, Pocillopora damicornis, across 10 degrees of latitude and a transition from temperate to tropical waters. We identified strong patterns of differentiation and reduced genetic diversity in high‐latitude populations. In addition, genome‐wide scans for selection identified a number of outlier loci putatively under directional selection with homology to proteins previously known to be involved in heat tolerance in corals and associated with processes such as photoprotection, protein degradation, and immunity. This study provides genomic evidence for both restricted gene flow and local adaptation in a widely distributed coral species, and highlights the potential vulnerability of leading‐edge populations to rapid environmental change as they are locally adapted, reproductively isolated, and have reduced levels of genetic diversity.     

Elevated Temperature Alters the Lunar Timing of Planulation in the Brooding Coral Pocillopora damicornis

Reproductive timing in corals is associated with environmental variables including temperature, lunar periodicity, and seasonality. Although it is clear that these variables are interrelated, it remains unknown if one variable in particular acts as the proximate signaler for gamete and or larval release. Furthermore, in an era of global warming, the degree to which increases in ocean temperatures will disrupt normal reproductive patterns in corals remains unknown. Pocillopora damicornis, a brooding coral widely distributed in the Indo-Pacific, has been the subject of multiple reproductive ecology studies that show correlations between temperature, lunar periodicity, and reproductive timing. However, to date, no study has empirically measured changes in reproductive timing associated with increased seawater temperature. In this study, the effect of increased seawater temperature on the timing of planula release was examined during the lunar cycles of March and June 2012. Twelve brooding corals were removed from Hobihu reef in Nanwan Bay, southern Taiwan and placed in 23 and 28°C controlled temperature treatment tanks. For both seasons, the timing of planulation was found to be plastic, with the high temperature treatment resulting in significantly earlier peaks of planula release compared to the low temperature treatment. This suggests that temperature alone can influence the timing of larval release in Pocillopora damicornis in Nanwan Bay. Therefore, it is expected that continued increases in ocean temperature will result in earlier timing of reproductive events in corals, which may lead to either variations in reproductive success or phenotypic acclimatization.     

Diversity and distribution of symbiodinium associated with seven common coral species in the Chagos Archipelago, central Indian Ocean

The Chagos Archipelago designated as a no-take marine protected area in 2010, lying about 500 km south of the Maldives in the Indian Ocean, has a high conservation priority, particularly because of its fast recovery from the ocean-wide massive coral mortality following the 1998 coral bleaching event. The aims of this study were to examine Symbiodinium diversity and distribution associated with scleractinian corals in five atolls of the Chagos Archipelago, spread over 10,000 km(2). Symbiodinium clade diversity in 262 samples of seven common coral species, Acropora muricata, Isopora palifera, Pocillopora damicornis, P. verrucosa, P. eydouxi, Seriatopora hystrix, and Stylophora pistillata were determined using PCR-SSCP of the ribosomal internal transcribed spacer 1 (ITS1), PCR-DDGE of ITS2, and phylogenetic analyses. The results indicated that Symbiodinium in clade C were the dominant symbiont group in the seven coral species. Our analysis revealed types of Symbiodinium clade C specific to coral species. Types C1 and C3 (with C3z and C3i variants) were dominant in Acroporidae and C1 and C1c were the dominant types in Pocilloporidae. We also found 2 novel ITS2 types in S. hystrix and 1 novel ITS2 type of Symbiodinium in A. muricata. Some colonies of A. muricata and I. palifera were also associated with Symbiodinium A1. These results suggest that corals in the Chagos Archipelago host different assemblages of Symbiodinium types then their conspecifics from other locations in the Indian Ocean; and that future research will show whether these patterns in Symbiodinium genotypes may be due to local adaptation to specific conditions in the Chagos.     

Coral-bacterial communities before and after a coral mass spawning event on Ningaloo Reef

Bacteria associated with three coral species, Acropora tenuis, Pocillopora damicornis and Tubastrea faulkneri, were assessed before and after coral mass spawning on Ningaloo Reef in Western Australia. Two colonies of each species were sampled before and after the mass spawning event and two additional samples were collected for P. damicornis after planulation. A variable 470 bp region of the 16 S rRNA gene was selected for pyrosequencing to provide an understanding of potential variations in coral-associated bacterial diversity and community structure. Bacterial diversity increased for all coral species after spawning as assessed by Chao1 diversity indicators. Minimal changes in community structure were observed at the class level and data at the taxonomical level of genus incorporated into a PCA analysis indicated that despite bacterial diversity increasing after spawning, coral-associated community structure did not shift greatly with samples grouped according to species. However, interesting changes could be detected from the dataset; for example, α-Proteobacteria increased in relative abundance after coral spawning and particularly the Roseobacter clade was found to be prominent in all coral species, indicating that this group may be important in coral reproduction.     

Reproductive character displacement by inversion of coiling in clausiliid snails (Gastropoda, Pulmonata)

In land snails, a change in the direction of coiling, being associated with a shift in the position of the genital apparatus, may act as a barrier against hybridization between sympatric species. Putative reproductive character displacement by an inversion in chirality has been reported in only a few land snails, based on observations in the field and interbreeding experiments. In this study, we present a new case of possible reproductive character displacement in the direction of coiling, in the clausiliid snail Isabellaria dextrorsa . This species is dextral, in contrast with its nearest relatives, including I. torifera and I. lophauchena , which share plesiomorphic sinistral coiling. Whereas I. dextrorsa occurs in sympatry and even syntopically with I. lophauchena throughout most of its range, the sinistral species have a mosaic distribution. Phylogenetic analyses of mitochondrial cytochrome c oxidase subunit I (COI) sequences demonstrated that I. dextrorsa constitutes a clade with I. torifera . In this clade, a shift in coiling direction occurred at least twice, maybe triggered by the presence of a sympatric congeneric sinistral species. The analyses separated the sequences of all I. dextrorsa samples from those of sympatric and syntopic I. lophauchena samples. The failure to demonstrate gene flow between these species is consistent with the hypothesis of genetic isolation by reproductive character displacement.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 88 , 155–164.     

Natural hybridization and reproductive isolation between two Primula species

Natural hybridization frequently occurs in plants and can facilitate gene flow between species,possibly resulting in species refusion. However,various reproductive barriers block the formation of hybrids and maintain species integrity. Here,we conducted a field survey to examine natural hybridization and reproductive isolation(RI) between sympatric populations of Primula secundiflora and P. poissonii using ten nuclear simple sequence repeat(SSR) loci. Although introgressive hybridization occurred,species boundaries between P. secundiflora and P. poissonii were maintained through nearly complete reproductive isolation. These interfertile species provide an excellent model for studying the RI mechanisms and evolutionary forces that maintain species boundaries.     

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell’s wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.Subject terms: Population genetics, Speciation     

Fatty acid variations in symbiotic dinoflagellates from Okinawan corals

The fatty acid composition of polar lipids and triacylglycerols was determined in different morphophysiological types of symbiotic dinoflagellates (SD) isolated from the hydrocoral Millepora intricata and the scleractinian corals Pocillopora damicornis, Seriatopora caliendrum, Seriatopora hystrix and Stylophora pistillata from a fringing reef of Sesoko Island, Okinawa, Japan. The distribution of the fatty acids among the morphophysiologically distinct types of SD reported in these corals makes it possible to readily distinguish one type of SD from the other. Moreover, differences were found both in polar lipids and triacylglycerols. The polar lipids of SD from M. intricata showed a very distinctive fatty acid profile. A combination of large proportions of 18:4 (n-3), 18:5 (n-3), 22:5 (n-6), and 22:6 (n-3) and negligible amounts of 20:4 (n-6), and 20:5 (n-3) in SD from M. intricata was particularly noteworthy. The fatty acid profiles of SD from P. damicornis and SD isolated from S. caliendrum and S. hystrix differed in the proportion of 18:4 (n-3) and 22:6 (n-3). It is suggested that fatty acids might provide useful information on possible taxonomic differences among symbiotic dinoflagellates. It is assumed that biochemical differences can reflect the genetic diversity of the morphophysiological types of SD associated with several species of hermatypic corals from this region.     

History rather than hybridization determines population structure and adaptation in Populus balsamifera

Hybridization between species is known to greatly affect their genetic diversity and, therefore, their evolution. Also, within species, there may be genetic clusters between which gene flow is limited, which may impact natural selection. However, few studies have looked simultaneously at the influence of among‐species and within‐species gene flow. Here, we study the influence of hybridization between Populus balsamifera and Populus trichocarpa on population structure and adaptation in P. balsamifera. We did this by sampling a total of 1517 individuals from across the ranges of these two species, and by genotyping them using a combination of 93 nuclear and 17 cpDNA SNPs. We found that hybridization is mostly limited to the contact zone where the species’ distributions overlap. Within P. balsamifera, we found multiple levels of population structure. Interestingly, the border between the Eastern and Central clusters is very sharp, whereas the border between the Central and Western clusters is diffuse. Outlier analysis revealed that three loci associated with the sharp border were also associated with climate. We hypothesize that the observed clusters derive from three refugia during the Pleistocene ice ages. Between the Central and Western clusters, post‐glacial long‐distance gene flow has led to the diffusion of their border. In the Eastern cluster, we hypothesize that endogenous genomic barriers have developed, leading to the sharp border and a spurious climate association. We conclude that the large‐scale genetic structure of P. balsamifera is mostly shaped by historical factors and the influence of interspecific hybridization is limited.     

A sea water barrier to coral gene flow

Land is not the only barrier to dispersal encountered by marine organisms. For sedentary shallow water species, there is an additional, marine barrier, 5000 km of uninterrupted deep‐water stretch between the central and the eastern Pacific. This expanse of water, known as the ‘Eastern Pacific Barrier’, has been separating faunas of the two oceanic regions since the beginning of the Cenozoic. Species with larvae that cannot stay in the plankton for the time it takes to cross between the two sides have been evolving independently. That the eastern Pacific does not share species with the rest of the Pacific was obvious to naturalists two centuries ago (Darwin 1860). Yet, this rule has exceptions. A small minority of species are known to straddle the Eastern Pacific Barrier. One such exception is the scleractinian coral Porites lobata (Fig.  1 ). This species is spread widely throughout the Indo‐Pacific, where it is one of the major reef‐builders, but it is also encountered in the eastern Pacific. Are eastern and central Pacific populations of this coral connected by gene flow? In this issue of Molecular Ecology, Baums et al. (2012) use microsatellite data to answer this question. They show that P. lobata populations in the eastern Pacific are cut off from genetic influx from the rest of the Pacific. Populations within each of the two oceanic regions are genetically connected (though those in the Hawaiian islands are also isolated). Significantly, the population in the Clipperton Atoll, the westernmost island in the eastern Pacific, genetically groups with populations from the central Pacific, suggesting that crossing the Eastern Pacific Barrier by P. lobata propagules does occasionally occur.     

Divergence population genetic analysis of hybridization between rhesus and cynomolgus macaques

The geographic ranges of rhesus ( Macaca mulatta ) and cynomolgus ( M. fascicularis ) macaques adjoin in Indochina where they appear to hybridize. We used published and newly generated DNA sequences from 19 loci spanning ~20 kb to test whether introgression has occurred between these macaque species. We studied introgression at the level of nuclear DNA and distinguished between incomplete lineage sorting of ancestral polymorphisms or interspecific gene flow. We implemented a divergence population genetics approach by fitting our data to an isolation model implemented in the software IMa. The model that posits no gene flow from the rhesus into the cynomolgus macaque was rejected ( P  = 1.99 × 10⁻⁸). Gene flow in this direction was estimated as 2 Nm ~1.2, while gene flow in the reverse direction was nonsignificantly different from zero ( P  = 0.16). The divergence time between species was estimated as ~1.3 million years. Balancing selection, a special case of incomplete sorting, was taken into consideration, as well as potential crossbreeding in captivity. Parameter estimates varied between analyses of subsets of data, although we still rejected isolation models. Geographic sampling of the data, where samples of cynomolgus macaques derived from Indochina were excluded, revealed a lost signature of gene flow, indicating that interspecific gene flow is restricted to mainland Indochina. Our results, in conjunction with those by others, justify future detailed analyses into the genetics of reproductive barriers and reticulate evolution in these two genome-enabled primates. Future studies of the natural hybridization between rhesus and cynomolgus macaques would expand the repertoire of systems available for speciation studies in primates.     

Gene flow of Acanthaster planci (L.) in relation to ocean currents revealed by microsatellite analysis

Population outbreaks of the coral-eating starfish, Acanthaster planci , are hypothesized to spread to many localities in the Indo-Pacific Ocean through dispersal of planktonic larvae. To elucidate the gene flow of A. planci across the Indo-Pacific in relation to ocean currents and to test the larval dispersal hypothesis, the genetic structure among 23 samples over the Indo-Pacific was analysed using seven highly polymorphic microsatellite loci. The F -statistics and genetic admixture analysis detected genetically distinct groups in accordance with ocean current systems, that is, the Southeast African group (Kenya and Mayotte), the Northwestern Pacific group (the Philippines and Japan), Palau, the North Central Pacific group (Majuro and Pohnpei), the Great Barrier Reef, Fiji, and French Polynesia, with a large genetic break between the Indian and Pacific Oceans. A pattern of significant isolation by distance was observed among all samples ( P =  0.001, r  = 0.88, n  = 253, Mantel test), indicating restricted gene flow among the samples in accordance with geographical distances. The data also indicated strong gene flow within the Southeast African, Northwestern Pacific, and Great Barrier Reef groups. These results suggest that the western boundary currents have strong influence on gene flow of this species and may trigger secondary outbreaks.     

Genomic footprints of an old affair: Single nucleotide polymorphism data reveal historical hybridization and the subsequent evolution of reproductive barriers in two recently diverged grasshoppers with partly overlapping distributions

Secondary contact in close relatives can result in hybridization and the admixture of previously isolated gene pools. However, after an initial period of hybridization, reproductive isolation can evolve through different processes and lead to the interruption of gene flow and the completion of the speciation process. Omocestus minutissimus and O. uhagonii are two closely related grasshoppers with partially overlapping distributions in the Central System mountains of the Iberian Peninsula. To analyse spatial patterns of historical and/or contemporary hybridization between these two taxa and understand how species boundaries are maintained in the region of secondary contact, we sampled sympatric and allopatric populations of the two species and obtained genome‐wide single nucleotide polymorphism data using a restriction site‐associated DNA sequencing approach. We used Bayesian clustering analyses to test the hypothesis of contemporary hybridization in sympatric populations and employed a suite of phylogenomic approaches and a coalescent‐based simulation framework to evaluate alternative hypothetical scenarios of interspecific gene flow. Our analyses rejected the hypothesis of contemporary hybridization but revealed past introgression in the area where the distributions of the two species overlap. Overall, these results point to a scenario of historical gene flow after secondary contact followed by the evolution of reproductive isolation that currently prevents hybridization among sympatric populations.