Massive transoceanic gene flow in a freshwater turtle (Testudines: Geoemydidae: Mauremys rivulata)

The freshwater turtle Mauremys rivulata ranges from the Adriatic coast of the Balkan Peninsula through the Aegean region and coastal western and southern Turkey southwards to Israel. In addition, it occurs on several Aegean islands, Crete and Cyprus. Previous investigations using mtDNA sequences found virtually no genetic differentiation across its distribution range, despite some major biogeographical barriers for terrestrial and freshwater biota. Thus, the absence of any phylogeographical differentiation would be unexpected. To re‐examine genetic differentiation within M. rivulata, here we use a comprehensive rangewide sampling and information of 13 unlinked polymorphic microsatellite loci and compare these data against mtDNA variation. Our microsatellite analyses reveal a weak population structuring which conflicts, however, with most biogeographical barriers. We conclude that the genetic structure in the vast majority of the species' range has been shaped by massive transoceanic gene flow. This explanation is unlikely for the northernmost populations, which seem rather to be genetically impacted by intentionally released foreign turtles.     

Population structure and gene flow of the syntopic turtles Emys and Mauremys from coastal and inland regions of Anatolia (Turkey): results from mitochondrial and microsatellite data

Revealing the genetic basis of the existence of different species living together in different geographic regions provides clarification of this phylogeographic differentiation. In this study, we investigated the population genetics and evaluated the level of genetic variation of inland and coastal populations of Mauremys and Emys in Turkey. Tissue samples of 196 terrapins were studied which were collected from syntopic coastal (Gölbent-Söke/Ayd?n; M. rivulata and E. orbicularis) and inland populations (Bahçesaray/Aksaray; M. caspica and E. orbicularis). DNA was isolated using the InnuPREP DNA Mini Kit. Mitochondrial DNA sequences and allelic variation at 13 microsatellite loci for Mauremys and 12 microsatellite loci for Emys were examined.  Three haplotypes were found for Emys orbicularis (Im, Ip and Iw) collected from the coastal region and two haplotypes for Emys orbicularis (Ig and Im) collected from inland. Two haplotypes were identified for M. caspica (Cmt8 and Cmt9) and three haplotypes were identified for M. rivulata (Rmt3, Rmt24 and Rmt26). Using microsatellites and the software STRUCTURE the most probable value for K was revealed as two 2 for both species. The F_ST value between M. rivulata and M. caspica was 0.39, and between the coastal and inland populations of E. orbicularis 0.09. It can be concluded that Emys populations tend to evolve by somehow preserving the allelic richness they have and Mauremys populations continue to differentiate so that new species emerge in the evolutionary process to reach the ideal allelic structure.

     

Human disturbance causes the formation of a hybrid swarm between two naturally sympatric fish species

Most evidence for hybrid swarm formation stemming from anthropogenic habitat disturbance comes from the breakdown of reproductive isolation between incipient species, or introgression between allopatric species following secondary contact. Human impacts on hybridization between divergent species that naturally occur in sympatry have received considerably less attention. Theory predicts that reinforcement should act to preserve reproductive isolation under such circumstances, potentially making reproductive barriers resistant to human habitat alteration. Using 15 microsatellites, we examined hybridization between sympatric populations of alewife (Alosa pseudoharengus) and blueback herring (A. aestivalis) to test whether the frequency of hybridization and pattern of introgression have been impacted by the construction of a dam that isolated formerly anadromous populations of both species in a landlocked freshwater reservoir. The frequency of hybridization and pattern of introgression differed markedly between anadromous and landlocked populations. The rangewide frequency of hybridization among anadromous populations was generally 0–8%, whereas all landlocked individuals were hybrids. Although neutral introgression was observed among anadromous hybrids, directional introgression leading to increased prevalence of alewife genotypes was detected among landlocked hybrids. We demonstrate that habitat alteration can lead to hybrid swarm formation between divergent species that naturally occur sympatrically, and provide empirical evidence that reinforcement does not always sustain reproductive isolation under such circumstances.     

Introgression between invasive and native blue mussels (genus Mytilus) in the central California hybrid zone

The ecological and genetic factors determining the extent of introgression between species in secondary contact zones remain poorly understood. Here, we investigate the relative importance of isolating barriers and the demographic expansion of invasive Mytilus galloprovincialis on the magnitude and the direction of introgression with the native Mytilus trossulus in a hybrid zone in central California. We use double‐digest restriction‐site‐associated DNA sequencing (ddRADseq) to genotype 1337 randomly selected single nucleotide polymorphisms and accurately distinguish early and advanced generation hybrids for the first time in the central California Mytilus spp. hybrid zone. Weak levels of introgression were observed in both directions but were slightly more prevalent from the native M. trossulus into the invasive M. galloprovincialis. Few early and advanced backcrossed individuals were observed across the hybrid zone confirming the presence of strong barriers to interbreeding. Heterogeneous patterns of admixture across the zone of contact were consistent with the colonization history of M. galloprovincialis with more extensive introgression in northern localities furthest away from the putative site of introduction in southern California. These observations reinforce the importance of dynamic spatial and demographic expansions in determining patterns of introgression between close congeners, even in those with high dispersal potential and well‐developed reproductive barriers. Our results suggest that the threat posed by invasive M. galloprovincialis is more ecological than genetic as it has displaced, and continues to displace the native M. trossulus from much of central and southern California.     

Combining hydrodynamic modelling with genetics: can passive larval drift shape the genetic structure of Baltic Mytilus populations?

While secondary contact between Mytilus edulis and Mytilus trossulus in North America results in mosaic hybrid zone formation, both species form a hybrid swarm in the Baltic. Despite pervasive gene flow, Baltic Mytilus species maintain substantial genetic and phenotypic differentiation. Exploring mechanisms underlying the contrasting genetic composition in Baltic Mytilus species will allow insights into processes such as speciation or adaptation to extremely low salinity. Previous studies in the Baltic indicated that only weak interspecific reproductive barriers exist and discussed the putative role of adaptation to environmental conditions. Using a combination of hydrodynamic modelling and multilocus genotyping, we investigate how oceanographic conditions influence passive larval dispersal and hybrid swarm formation in the Baltic. By combining our analyses with previous knowledge, we show a genetic transition of Baltic Mytilus species along longitude 12°‐13°E, that is a virtual line between Malmö (Sweden) and Stralsund (Germany). Although larval transport only occurs over short distances (10–30 km), limited larval dispersal could not explain the position of this genetic transition zone. Instead, the genetic transition zone is located at the area of maximum salinity change (15–10 psu). Thus, we argue that selection results in weak reproductive barriers and local adaptation. This scenario could maintain genetic and phenotypic differences between Baltic Mytilus species despite pervasive introgressive hybridization.     

Bioinvasion threatens the genetic integrity of native diversity and a natural hybrid zone: smooth‐shelled blue mussels (Mytilus spp.) in the Strait of Magellan

Smooth‐shelled blue mussels of the Mytilus edulis species complex are widely distributed bivalve molluscs whose introductions threaten native marine biodiversity (non‐indigenous species – NIS). The aim of the present study was to identify the species and hybrids of Mytilus present in the Magellan Region (southern Chile). Results indicate that three mussel species of the Mytilus edulis complex are found in the region – M. edulis, M. chilensis (or the Southern Hemisphere lineage of Mytilus galloprovincialis), and M. galloprovincialis of Northern Hemisphere origin. For the first time, alleles of the introduced M. trossulus are reported from the Southern Hemisphere. In the Strait of Magellan the native Pacific blue mussel, Mytilus chilensis and the native Atlantic blue mussel, Mytilus edulis, meet and mix at a natural hybrid zone (about 125 km in length). This is the first record of a natural Mytilus hybrid zone in the Southern Hemisphere and is also the first record of the co‐occurrence of genes from all four Mytilus species in any one region. These results contribute to the knowledge of the biodiversity and delimitation of mussel species in southern South America, and highlight how introduced species may threaten the genetic integrity of native species through hybridization and introgression.     

Genomic divergence and lack of introgressive hybridization between two 13‐year periodical cicadas support life cycle switching in the face of climate change

Life history evolution spurred by post‐Pleistocene climatic change is hypothesized to be responsible for the present diversity in periodical cicadas (Magicicada), but the mechanism of life cycle change has been controversial. To understand the divergence process of 13‐year and 17‐year cicada life cycles, we studied genetic relationships between two synchronously emerging, parapatric 13‐year periodical cicada species in the Decim group, Magicicada tredecim and M. neotredecim. The latter was hypothesized to be of hybrid origin or to have switched from a 17‐year cycle via developmental plasticity. Phylogenetic analysis using restriction‐site‐associated DNA sequences for all Decim species and broods revealed that the 13‐year M. tredecim lineage is genomically distinct from 17‐year Magicicada septendecim but that 13‐year M. neotredecim is not. We detected no significant introgression between M. tredecim and M. neotredecim/M. septendecim thus refuting the hypothesis that M. neotredecim are products of hybridization between M. tredecim and M. septendecim. Further, we found that introgressive hybridization is very rare or absent in the contact zone between the two 13‐year species evidenced by segregation patterns in single nucleotide polymorphisms, mitochondrial lineage identity and head width and abdominal sternite colour phenotypes. Our study demonstrates that the two 13‐year Decim species are of independent origin and nearly completely reproductively isolated. Combining our data with increasing observations of occasional life cycle change in part of a cohort (e.g. 4‐year acceleration of emergence in 17‐year species), we suggest a pivotal role for developmental plasticity in Magicicada life cycle evolution.     

Preliminary analysis of a hybrid zone between two subspecies of Zootermopsis nevadensis

For species largely allopatric in distribution, zones of contact provide an opportunity for hybridization, testing grounds for species boundaries, and may result in the formation of a new species. Thus, hybrid zones have the potential to provide important insights into speciation. In this study, we performed a preliminary analysis of a hybrid zone between two subspecies of the dampwood termite, Zootermopsis nevadensis (Z. n. nuttingi and Z. n. nevadensis) near Bartle, CA, using 12 microsatellite loci and a mitochondrial gene. Fifty-seven colonies collected in 36 locations were analyzed. The analysis of genetic markers revealed a large hybrid zone approximately 104 km in width. Although stepped clines best explained the data, we are unable to rule out the existence of a mosaic hybrid zone. We inferred weak selection (s* < 0.05%) against hybrids, but the data also suggested the existence of a barrier to gene flow from Z. n. nevadensis to Z. n. nuttingi, but not in the other direction. Given the large zone of contact, extensive sampling is needed to obtain a more complete characterization of the hybrid zone. However, the results of this study suggest that despite the accumulation of phenotypic and genetic differences and time since divergence (~2 million years), Z. n. nuttingi and Z. n. nevadensis are capable of extensive hybridization.     

Molecular patterns of introgression in a classic hybrid zone between the Australian tree frogs,Litoria ewingii and L. paraewingi: evidence of a tension zone

Hybrid zones provide a rare opportunity to explore the processes involved in reproductive isolation and speciation. The southern hybrid zone between the southeastern Australian tree frogs Litoria ewingii and L. paraewingi has been comprehensively studied over the last 40 years, primarily using reproductive compatibility experiments and male advertisement calls. We used mitochondrial DNA (mtDNA) and eight nuclear microsatellite markers to characterize this hybrid zone along a historically studied transect and to test various dispersal‐dependent and dispersal‐independent hybrid zone models. The species are genetically distinct and the level of hybridization within the contact zone is low, with the majority of admixed individuals representing later‐generation hybrids. Based on previous experimental genetic compatibility studies, we predicted that hybrids with L. paraewingi mtDNA would be more frequent than hybrids with L. ewingii mtDNA. Surprisingly, a greater proportion of the identified hybrids had L. ewingii mtDNA. Geographical cline analyses showed a sharp transition in allele frequencies across the transect, and both the mtDNA and microsatellite data showed concordant cline centres, but were best supported by a model that allowed width to vary. Overall, the L. ewingii–L. paraewingi hybrid zone is best characterized as a tension zone, due to the narrow cline width, concordant genetic clines and low levels of hybridization.     

Influence of niche similarity on hybridization between Myriophyllum sibiricum and M. spicatum

The impact of ecological factors on natural hybridization is of widespread interest. Here, we asked whether climate niche influences hybridization between the two closely related plant species Myriophyllum sibiricum and M. spicatum. Eight microsatellite loci and two chloroplast fragments were used to investigate the occurrence of hybridization between these two species in two co‐occurring regions: north‐east China (NEC) and the Qinghai‐Tibetan Plateau (QTP). The climate niches of the species were quantified by principal component analysis with bioclimatic data, and niche comparisons were performed between the two species in each region. Reciprocal hybridization was observed, and M. sibiricum was favoured as the maternal species. Furthermore, hybrids were rare in NEC but common in the QTP. Accordingly, in NEC, the two species were climatically distinct, and hybrids only occurred in the narrow geographical or ecological transition zone, whereas in the QTP, obvious niche overlaps were found for the two species, and hybrids occurred in multiple contact zones. This association between hybridization pattern and climate niche similarity suggests that the level of hybridization was promoted by niche overlap. Compared with the parental species, similar climate niches were found for the hybrid populations in the QTP, indicating that other environmental factors rather than climate were important for hybrid persistence. Our findings highlight the significance of climate niche with respect to hybridization patterns in plants.     

Secondary contact between Lycaeides idas and L. melissa in the Rocky Mountains: extensive admixture and a patchy hybrid zone

Studies of hybridization have increased our understanding of the nature of species boundaries, the process of speciation, and the effects of hybridization on the evolution of populations and species. In the present study we use genetic and morphological data to determine the outcome and consequences of secondary contact and hybridization between the butterfly species Lycaeides idas and L. melissa in the Rocky Mountains. Admixture proportions estimated from structure and geographical cline analysis indicate L. idas and L. melissa have hybridized extensively in the Rocky Mountains and that reproductive isolation was insufficient to prevent introgression for much of the genome. Geographical patterns of admixture suggest that hybridization between L. idas and L. melissa has led to the formation of a hybrid zone. The hybrid zone is relatively wide, given estimates of dispersal for Lycaeides butterflies, and does not show strong evidence of cline concordance among characters. We believe the structure of the Lycaeides hybrid zone might be best explained by the patchy distribution of Lycaeides, local extinction and colonization of habitat patches, environmental variation and weak overall selection against hybrids. We found no evidence that hybridization in the Rocky Mountains has resulted in the formation of independent hybrid species, in contrast to the outcome of hybridization between L. idas and L. melissa in the Sierra Nevada. Finally, our results suggest that differences in male morphology between L. idas and L. melissa might contribute to isolation, or perhaps even that selection has favoured the spread of L. melissa male genitalia alleles.     

Multiple origins of parthenogenesis,and a revised species phylogeny for the Southeast Asian butterfly lizards,Leiolepis

Hybridization often occurs in areas of secondary contact between closely related species. In some cases these hybridization events can create hybrid offspring that are reproductively viable as new parthenogenetic species. The genus Leiolepis contains nine species that collectively range throughout continental Southeast Asia. Of these, four are unisexual (some diploid and some triploid). We analyzed a multi‐locus dataset within a multi‐lineage coalescent framework to infer the origins of these parthenogenetic hybrid species. Our results provide evidence that repeated hybridization events between L. reevesii and L. guttata have led to the formation of all four distinct parthenogenetic species. Our data further suggest there have been low levels of mitochondrial introgression between L. belliana and L. reevesii at their contact zone in southern Cambodia. This work addresses contentious species boundaries and provides the first taxon‐complete hypothesis of relationships for the butterfly lizards. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 1080–1093.     

Tracing species replacement in Iberian marbled newts

Secondary contact between closely related species can lead to the formation of hybrid zones, allowing for interspecific gene flow. Hybrid zone movement can take place if one of the species possesses a competitive advantage over the other, ultimately resulting in species replacement. Such hybrid zone displacement is predicted to leave a genomic footprint across the landscape in the form of asymmetric gene flow (or introgression) of selectively neutral alleles from the displaced to the advancing species. Hybrid zone movement has been suggested for marbled newts in the Iberian Peninsula, supported by asymmetric gene flow and a distribution relict (i.e., an enclave) of Triturus marmoratus in the range of T. pygmaeus. We developed a panel of nuclear and mitochondrial SNP markers to test for the presence of a T. marmoratus genomic footprint in the Lisbon peninsula, south of the enclave. We found no additional populations of T. marmoratus. Analysis with the software Structure showed no genetic traces of T. marmoratus in T. pygmaeus. A principal component analysis showed some variation within the local T. pygmaeus, but it is unclear if this represents introgression from T. marmoratus. The results may be explained by (a) species replacement without introgressive hybridization and (b) displacement with hybridization followed by the near‐complete erosion of the footprint by purifying selection. We predict that testing for a genomic footprint north of the reported enclave would confirm that species replacement in these marbled newts occurred with hybridization.     

Temporal and spatial dynamics of competitive parapatry in chewing lice

We synthesize observations from 1979 to 2016 of a contact zone involving two subspecies of pocket gophers (Thomomys bottae connectens and T. b. opulentus) and their respective chewing lice (Geomydoecus aurei and G. centralis) along the Río Grande Valley in New Mexico, U.S.A., to test predictions about the dynamics of the zone. Historically, the natural flood cycle of the Rio Grande prevented contact between the two subspecies of pocket gophers. Flood control measures completed in the 1930s permitted contact, thus establishing the hybrid zone between the pocket gophers and the contact zone between their lice (without hybridization). Since that time, the pocket gopher hybrid zone has stabilized, whereas the northern chewing louse species has replaced the southern louse species at a consistent rate of ~150 m/year. The 0.2–0.8 width of the replacement zone has remained constant, reflecting the constant rate of chewing louse species turnover on a single gopher and within a local pocket gopher population. In contrast, the full width of the replacement zone (northernmost G. centralis to southernmost G. aurei) has increased annually. By employing a variety of metrics of the species replacement zone, we are better able to understand the dynamics of interactions between and among the chewing lice and their pocket gopher hosts. This research provides an opportunity to observe active species replacement and resulting distributional shifts in a parasitic organism in its natural setting.     

Genetic and morphological evidence of a geographically widespread hybrid zone between two crocodile species,Crocodylus acutus and Crocodylus moreletii

Hybrid zones represent natural laboratories to study gene flow, divergence and the nature of species boundaries between closely related taxa. We evaluated the level and extent of hybridization between Crocodylus moreletii and Crocodylus acutus using genetic and morphological data on 300 crocodiles from 65 localities. To our knowledge, this is the first genetic study that includes the entire historic range and sympatric zone of the two species. Contrary to expectations, Bayesian admixture proportions and maximum‐likelihood estimates of hybrid indexes revealed that most sampled crocodiles were admixed and that the hybrid zone is geographically extensive, extending well beyond their historical region of sympatry. We identified a few geographically isolated, nonadmixed populations of both parental species. Hybrids do not appear to be F₁s or recent backcrosses, but rather are more likely later‐generation hybrids, suggesting that hybridization has been going on for several to many generations and is mostly the result of natural processes. Crocodylus moreletii is not the sister species of C. acutus, suggesting that the hybrid zone formed from secondary contact rather than primary divergence. Nonadmixed individuals from the two species were distinguishable based on morphological characters, whereas hybrids had a complex mosaic of morphological characters that hinders identification in the wild. Very few nonadmixed C. acutus and C. moreletii populations exist in the wild. Consequently, the last nonadmixed C. moreletii populations have become critically endangered. Indeed, not only the parental species but also the naturally occurring hybrids should be considered for their potential conservation value.     

Females drive asymmetrical introgression from rare to common species in Darwin's tree finches

The consequences of hybridization for biodiversity depend on the specific ecological and evolutionary context in which it occurs. Understanding patterns of gene flow among hybridizing species is crucial for determining the evolutionary trajectories of species assemblages. The recently discovered hybridization between two species of Darwin's tree finches (Camarhynchus parvulus and C. pauper) on Floreana Island, Galápagos, presents an exciting opportunity to investigate the mechanisms causing hybridization and its potential evolutionary consequences under conditions of recent habitat disturbance and the introduction of invasive pathogens. In this study, we combine morphological and genetic analysis with pairing observations to explore the extent, direction and drivers of hybridization and to test whether hybridization patterns are a result of asymmetrical pairing preference driven by females of the rarer species (C. pauper). We found asymmetrical introgression from the critically endangered, larger‐bodied C. pauper to the common, smaller‐bodied C. parvulus, which was associated with a lack of selection against heterospecific males by C. pauper females. Examination of pairing data showed that C. parvulus females paired assortatively, whereas C. pauper females showed no such pattern. This study shows how sex‐specific drivers can determine the direction of gene flow in hybridizing species. Furthermore, our results suggest the existence of a hybrid swarm comprised of C. parvulus and hybrid birds. We discuss the influence of interspecific abundance differences and susceptibility to the invasive parasite Philornis downsi on the observed hybridization and recommend that the conservation of this iconic species group should be managed jointly rather than species‐specific.     

Genetic variation in Plethodon cinereus and Plethodon hubrichti from in and around a contact zone

Climate change poses several challenges to biological communities including changes in the frequency of encounters between closely related congeners as a result of range shifts. When climate change leads to increased hybridization, hybrid dysfunction or genetic swamping may increase extinction risk—particularly in range‐restricted species with low vagility. The Peaks of Otter Salamander, Plethodon hubrichti, is a fully terrestrial woodland salamander that is restricted to ~18 km of ridgeline in the mountains of southwestern Virginia, and its range is surrounded by the abundant and widespread Eastern Red‐backed Salamander, Plethodon cinereus. In order to determine whether these two species are hybridizing and how their range limits may be shifting, we assessed variation at eight microsatellite loci and a 1,008 bp region of Cytochrome B in both species at allopatric reference sites and within a contact zone. Our results show that hybridization between P. hubrichti and P. cinereus either does not occur or is very rare. However, we find that diversity and differentiation are substantially higher in the mountaintop endemic P. hubrichti than in the widespread P. cinereus, despite similar movement ability for the two species as assessed by a homing experiment. Furthermore, estimation of divergence times between reference and contact zone populations via approximate Bayesian computation is consistent with the idea that P. cinereus has expanded into the range of P. hubrichti. Given the apparent recent colonization of the contact zone by P. cinereus, future monitoring of P. cinereus range limits should be a priority for the management of P. hubrichti populations.     

Phylogeography of the Lacerta viridis complex: mitochondrial and nuclear markers provide taxonomic insights

Based on broad, nearly rangewide sampling, we reanalysed the phylogeography of the Lacerta viridis complex using the mitochondrial cytochrome b gene and the intron 7 of the nuclear β‐fibrinogen gene. Using the mitochondrial marker, we identified in phylogenetic analyses 10 terminal clades clustering in four deeply divergent main lineages whose relationships are weakly resolved. These lineages correspond to Lacerta bilineata, L. viridis, the previously identified Adriatic or West Balkan lineage and a newly discovered fourth lineage from the Anatolian Black Sea coast and the south‐eastern Balkan Peninsula. Except for the latter lineage, there is considerable phylogeographic structuring in each lineage, with higher diversity in the south of the distribution ranges. This pattern indicates the existence of two distinct microrefugia in the Italian Peninsula and Sicily and of up to seven microrefugia in the Balkan Peninsula, but of only one refugium along the Black Sea coast of Anatolia. We identified secondary contact zones of the main lineages and of terminal clades within these lineages. However, most of the formerly described putative contact zone of L. bilineata and L. viridis turned out to be a contact zone between the Adriatic lineage and L. viridis, but L. bilineata seems to be involved only marginally. Our nuclear marker could not unambiguously resolve whether there is gene flow in contact zones. Thus, further research is necessary to decide whether the four main lineages are conspecific or whether they represent distinct biological species. We restrict the name L. v. meridionalis to the newly identified genetic lineage from Turkey and south‐eastern Europe, synonymize some previously recognized taxa and suggest a tentative nomenclature for the L. viridis complex.     

Does genetic distance between parental species influence outcomes of hybridization among coral reef butterflyfishes?

Christmas Island is located at the overlap of the Indian and Pacific Ocean marine provinces and is a hot spot for marine hybridization. Here, we evaluate the ecological framework and genetic consequences of hybridization between butterflyfishes Chaetodon guttatissimus and Chaetodon punctatofasciatus. Further, we compare our current findings to those from a previous study of hybridization between Chaetodon trifasciatus and Chaetodon lunulatus. For both species groups, habitat and dietary overlap between parental species facilitate frequent heterospecific encounters. Low abundance of potential mates promotes heterospecific pair formation and the breakdown of assortative mating. Despite similarities in ecological frameworks, the population genetic signatures of hybridization differ between the species groups. Mitochondrial and nuclear data from C. guttatissimus × C. punctatofasciatus (1% divergence at cyt b) show bidirectional maternal contributions and relatively high levels of introgression, both inside and outside the Christmas Island hybrid zone. In contrast, C. trifasciatus × C. lunulatus (5% cyt b divergence) exhibit unidirectional mitochondrial inheritance and almost no introgression. Back‐crossing of hybrid C. guttatissimus × C. punctatofasciatus and parental genotypes may eventually confound species‐specific signals within the hybrid zone. In contrast, hybrids of C. trifasciatus and C. lunulatus may coexist with and remain genetically distinct from the parents. Our results, and comparisons with hybridization studies in other reef fish families, indicate that genetic distance between hybridizing species may be a factor influencing outcomes of hybridization in reef fish, which is consistent with predictions from terrestrially derived hybridization theory.     

Comparative rangewide phylogeography of four endemic Taiwanese bat species

Phylogeographic reconstructions of codistributed taxa can help reveal the interplay between abiotic factors, such as altitude and climate, and species‐specific attributes, in shaping patterns of population genetic structure. Recent studies also demonstrate the value of both rangewide sampling and species distribution modelling (SDM) in comparative phylogeography. Here, we combine these approaches to study the population histories of four phylogenetically related forest‐dependent bat species. All are endemic to the mountainous island of Taiwan but show differences in their tolerance to altitude, with Murina gracilis considered to be a high‐altitude specialist, M. recondita and Kerivoula sp. low‐altitude specialists, and M. puta an altitudinal generalist. We tested the prediction that contrasting habitat preferences would impact on patterns of past and contemporary gene flow and found broad concordance between the results of population genetic analyses and species distribution models based on the Model for Interdisciplinary Research on Climate. Both lowland species showed evidence of genetic divergence between the east and west of the island, consistent with SDMs that indicated the Central Mountain Range (CMR) has presented a long‐term and continuous barrier to gene flow since before the Last Glacial Maximum. In contrast, Murina gracilis and M. puta showed lower degrees of historical isolation and genetic differentiation associated with the CMR, reflecting greater gene flow, possibly coupled with past population growth in M. puta. Together our results highlight the usefulness of combining distribution models with phylogeographic analyses to understand the drivers of genetic structure.