Interspecies introgressive hybridization in spiny frogs Quasipaa (Family Dicroglossidae) revealed by analyses on multiple mitochondrial and nuclear genes

Introgression may lead to discordant patterns of variation among loci and traits. For example, previous phylogeographic studies on the genus Quasipaa detected signs of genetic introgression from genetically and morphologically divergent Quasipaa shini or Quasipaa spinosa. In this study, we used mitochondrial and nuclear DNA sequence data to verify the widespread introgressive hybridization in the closely related species of the genus Quasipaa, evaluate the level of genetic diversity, and reveal the formation mechanism of introgressive hybridization. In Longsheng, Guangxi Province, signs of asymmetrical nuclear introgression were detected between Quasipaa boulengeri and Q. shini. Unidirectional mitochondrial introgression was revealed from Q. spinosa to Q. shini. By contrast, bidirectional mitochondrial gene introgression was detected between Q. spinosa and Q. shini in Lushan, Jiangxi Province. Our study also detected ancient hybridizations between a female Q. spinosa and a male Q. jiulongensis in Zhejiang Province. Analyses on mitochondrial and nuclear genes verified three candidate cryptic species in Q. spinosa, and a cryptic species may also exist in Q. boulengeri. However, no evidence of introgressive hybridization was found between Q. spinosa and Q. boulengeri. Quasipaa exilispinosa from all the sampling localities appeared to be deeply divergent from other communities. Our results suggest widespread introgressive hybridization in closely related species of Quasipaa and provide a fundamental basis for illumination of the forming mechanism of introgressive hybridization, classification of species, and biodiversity assessment in Quasipaa.     

Testing hypotheses for genealogical discordance in a rainforest lizard

Genealogical discordance, or when different genes tell distinct stories although they evolved under a shared history, often emerges from either coalescent stochasticity or introgression. In this study, we present a strong case of mito‐nuclear genealogical discordance in the Australian rainforest lizard species complex of Saproscincus basiliscus and S. lewisi. One of the lineages that comprises this complex, the Southern S. basiliscus lineage, is deeply divergent at the mitochondrial genome but shows markedly less divergence at the nuclear genome. By placing our results in a comparative context and reconstructing the lineages' demography via multilocus and coalescent‐based approximate Bayesian computation methods, we test hypotheses for how coalescent variance and introgression contribute to this pattern. These analyses suggest that the observed genealogical discordance likely results from introgression. Further, to generate such strong discordance, introgression probably acted in concert with other factors promoting asymmetric gene flow between the mitochondrial and nuclear genomes, such as selection or sex‐biased dispersal. This study offers a framework for testing sources of genealogical discordance and suggests that historical introgression can be an important force shaping the genetic diversity of species and their populations.     

Cytonuclear discordance among Southeast Asian black rats (Rattus rattus complex)

Black rats are major invasive vertebrate pests with severe ecological, economic and health impacts. Remarkably, their evolutionary history has received little attention, and there is no firm agreement on how many species should be recognized within the black rat complex. This species complex is native to India and Southeast Asia. According to current taxonomic classification, there are three taxa living in sympatry in several parts of Thailand, Cambodia and Lao People's Democratic Republic, where this study was conducted: two accepted species (Rattus tanezumi, Rattus sakeratensis) and an additional mitochondrial lineage of unclear taxonomic status referred to here as ‘Rattus R3’. We used extensive sampling, morphological data and diverse genetic markers differing in rates of evolution and parental inheritance (two mitochondrial DNA genes, one nuclear gene and eight microsatellite loci) to assess the reproductive isolation of these three taxa. Two close Asian relatives, Rattus argentiventer and Rattus exulans, were also included in the genetic analyses. Genetic analyses revealed discordance between the mitochondrial and nuclear data. Mitochondrial phylogeny studies identified three reciprocally monophyletic clades in the black rat complex. However, studies of the phylogeny of the nuclear exon interphotoreceptor retinoid‐binding protein gene and clustering and assignation analyses with eight microsatellites failed to separate R. tanezumi and R3. Morphometric analyses were consistent with nuclear data. The incongruence between mitochondrial and nuclear (and morphological) data rendered R. tanezumi/R3 paraphyletic for mitochondrial lineages with respect to R. sakeratensis. Various evolutionary processes, such as shared ancestral polymorphism and incomplete lineage sorting or hybridization with massive mitochondrial introgression between species, may account for this unusual genetic pattern in mammals.     

Phylogeography and lineage‐specific patterns of genetic diversity and molecular evolution in a group of North American skinks

Geography influences the evolutionary trajectory of species by mediating opportunities for hybridization, gene flow, demographic shifts and adaptation. We sought to understand how geography and introgression can generate species‐specific patterns of genetic diversity by examining phylogeographical relationships in the North American skink species Plestiodon multivirgatus and P. tetragrammus (Squamata: Scincidae). Using a multilocus dataset (three mitochondrial genes, four nuclear genes; a total of 3455 bp) we discovered mito‐nuclear discordance, consistent with mtDNA introgression. We further tested for evidence of species‐wide mtDNA introgression by using comparisons of genetic diversity, selection tests and extended Bayesian skyline analyses. Our findings suggest that P. multivirgatus acquired its mitochondrial genome from P. tetragrammus after their initial divergence. This putative species‐wide mitochondrial capture was further evidenced by statistically indistinguishable substitution rates between mtDNA and nDNA in P. multivirgatus. This rate discrepancy was observed in P. multivirgatus but not P. tetragrammus, which has important implications for studies that combine mtDNA and nDNA sequences when inferring time since divergence between taxa. Our findings suggest that by facilitating opportunities for interspecific introgression, geography can alter the course of molecular evolution between recently diverged lineages.     

Mito‐nuclear discordance with evidence of shared ancestral polymorphism and selection in cactophilic species of Drosophila

The Drosophila serido haplogroup is a monophyletic group composed of the following four cryptic and cactophilic species that are endemic to eastern Brazil: D. borborema, D. gouveai, D. seriema and D. serido. Here, we investigate the mito‐nuclear discordance in these species found among the cytochrome c oxidase subunit I (COI) mitochondrial gene, the autosomal alpha‐Esterase‐5 (α‐Est5) and the X‐linked period gene (per). Our analysis indicates that shared polymorphisms in these three molecular markers may be explained by the maintenance of ancestral polymorphisms rather than introgressive hybridization. The primary structures of COI, per and α‐Est5 genes evolve primarily under purifying selection, but we detected some sites that evolved under positive selection in α‐Est5. Considering the high variability of cacti species in eastern Brazil and the role attributed to Drosophila esterases in digestion metabolism and/or the detoxification of several compounds found in cactus tissues, we conjecture about the role of natural selection triggered by host shifts as an important factor in the intraspecific diversification of the D. serido haplogroup.     

Phylogeographic history of Patagonian lizards of the Liolaemus elongatus complex (Iguania: Liolaemini) based on mitochondrial and nuclear DNA sequences

In this study, we present a phylogeographic analysis of a group of lizards distributed in north‐western Patagonia, the Liolaemus elongatus complex. We sequenced 581 individuals for one mitochondrial gene (cytochrome‐b), and for a subset, we sequenced another mitochondrial gene (12S rRNA) and two nuclear genes: kinesin family member 24 (KIF24) and the anonymous nuclear locus LDAB1D. We estimated gene trees, mitochondrial and nuclear haploytpe networks, standard molecular diversity indices, genetic distances between lineages and Bayesian skyline plots. Our results provide evidence for recognition of seven species previously described within the L. elongatus complex: Liolaemus antumalguen, Liolaemus chillanensis, Liolaemus carlosgarini, Liolaemus burmeisteri, Liolaemus smaug, Liolaemus elongatus and Liolaemus crandalli, but we did not find sufficient evidence to support Liolaemus choique, Liolaemus shitan or Liolaemus sp. 6 as distinct species. We identified four candidate species (Liolaemus sp. 1, Liolaemus sp. 2, Liolaemus sp. 3 and Liolaemus sp. 7), and we discuss evolutionary processes that may have contributed to the origin of these lineages and their taxonomic and conservation implications.     

Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.)

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species‐specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri‐specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure ), bidirectional estimates of gene flow (migrate ‐n and IMa2) and phylogeny‐based methods (*beast , buck y) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species‐specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri‐like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.     

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito‐genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito‐genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex‐biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt‐b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt‐b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial–nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt‐b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito‐nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex‐biased dispersal. A weak pattern of female‐biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex‐biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris‐type mito‐genome in thermal adaptation may also contribute to the introgression between the two species.     

Phylogeography of the sergeants Abudefduf sexfasciatus and A. vaigiensis reveals complex introgression patterns between two widespread and sympatric Indo‐West Pacific reef fishes

On evolutionary timescales, sea level oscillations lead to recurrent spatio‐temporal variation in species distribution and population connectivity. In this situation, applying classical concepts of biogeography is challenging yet necessary to understand the mechanisms underlying biodiversity in highly diverse marine ecosystems such as coral reefs. We aimed at studying the outcomes of such complex biogeographic dynamics on reproductive isolation by sampling populations across a wide spatial range of a species‐rich fish genus: the sergeants (Pomacentridae: Abudefduf). We generated a mutlilocus data set that included ten morpho‐species from 32 Indo‐West Pacific localities. We observed a pattern of mito‐nuclear discordance in two common and widely distributed species: Abudefduf sexfasciatus and Abudefduf vaigiensis. The results showed three regional sublineages (Indian Ocean, Coral Triangle region, western Pacific) in A. sexfasciatus (0.6–1.5% divergence at cytb). The other species, A. vaigiensis, is polyphyletic and consists of three distinct genetic lineages (A, B and C) (9% divergence at cytb) whose geographic ranges overlap. Although A. vaigiensis A and A. sexfasciatus were found to be distinct based on nuclear information, A. vaigiensis A was found to be nested within A. sexfasciatus in the mitochondrial gene tree. A. sexfasciatus from the Coral Triangle region and A. vaigiensis A were not differentiated from each other at the mitochondrial locus. We then used coalescent‐based simulation to characterize a spatially widespread but weak gene flow between the two species. We showed that these fishes are good candidates to investigate the evolutionary complexity of the discrepancies between phenotypic and genetic similarity in closely related species.     

From hybridization to introgression between two closely related sympatric ant species

Interspecific hybridization is becoming more frequent worldwide due to increasing global changes and translocations of organisms. For individual organisms, the most significant negative consequences are sterility or inviability of hybrid offspring. However, hybridization sometimes leads to fertile offspring, promoting introgression from one species into another. In such situations, hybridization can play a key role in evolution and speciation. Combining hypervariable DNA (microsatellites) and mitochondrial DNA markers with the use of several modeling methods allow an efficient detection of hybridization processes. The present study therefore investigates hybridization between two ant species, Tetramorium immigrans and T. caespitum, using multiple methods, and systematically comparing results with simulated data to ensure accurate identification of hybrids. Introgression was revealed both by backcross detection based on 14 nuclear microsatellite loci and by mitochondrial‐nuclear discordance based on comparison with mitochondrial DNA (cytochrome c oxidase subunit I). Results were spatially consistent, with hybrids located at latitudes where parental species are sympatric. The causes and consequences of hybridization and introgression between T. caespitum and T. immigrans remain to be further investigated, especially because T. immigrans could be an invasive species in France.     

Diversification and asymmetrical gene flow across time and space: lineage sorting and hybridization in polytypic barking frogs

Young species complexes that are widespread across ecologically disparate regions offer important insights into the process of speciation because of their relevance to how local adaptation and gene flow influence diversification. We used mitochondrial DNA and up to 28 152 genomewide single nucleotide polymorphisms from polytypic barking frogs (Craugastor augusti complex) to infer phylogenetic relationships and test for the signature of introgressive hybridization among diverging lineages. Our phylogenetic reconstructions suggest (i) a rapid Pliocene–Pleistocene radiation that produced at least nine distinct lineages and (ii) that geographic features of the arid Central Mexican Plateau contributed to two independent northward expansions. Despite clear lineage differentiation (many private alleles and high between‐lineage F_ST scores), D‐statistic tests, which differentiate introgression from ancestral polymorphism, allowed us to identify two putative instances of reticulate gene flow. Partitioned D‐statistics provided evidence that these events occurred in the same direction between clades but at different points in time. After correcting for geographic distance, we found that lineages involved in hybrid gene flow interactions had higher levels of genetic variation than independently evolving lineages. These findings suggest that the nature of hybrid compatibility can be conserved overlong periods of evolutionary time and that hybridization between diverging lineages may contribute to standing levels of genetic variation.     

Mitonuclear discordance as a confounding factor in the DNA taxonomy of monogonont rotifers

Discordance between mitochondrial and nuclear phylogenies is being increasingly recognized in animals and may confound DNA‐based taxonomy. This is especially relevant for taxa whose microscopic size often challenges any effort to distinguish between cryptic species without the assistance of molecular data. Regarding mitonuclear discordance, two strikingly contrasting scenarios have been recently demonstrated in the monogonont rotifers of the genus Brachionus. While strict mitonuclear concordance was observed in the marine B. plicatilis species complex, widespread hybridization‐driven mitonuclear discordance was revealed in the freshwater B. calyciflorus species complex. Here, we investigated the frequency of occurrence and the potential drivers of mitonuclear discordance in three additional freshwater monogonont rotifer taxa, and assessed its potential impact on the reliability of DNA taxonomy results based on commonly used single markers. We studied the cryptic species complexes of Keratella cochlearis, Polyarthra dolichoptera and Synchaeta pectinata. Phylogenetic reconstructions were based on the mitochondrial barcoding marker cytochrome c oxidase subunit I gene and the nuclear internal transcribed spacer 1 locus, which currently represent the two most typical genetic markers used in rotifer DNA taxonomy. Species were delimited according to each marker separately using a combination of tree‐based coalescent, distance‐based and allele‐sharing‐based approaches. Mitonuclear discordance was observed in all species complexes with incomplete lineage sorting and unresolved phylogenetic reconstructions recognized as the likely drivers. Evidence from additional sources, such as morphology and ecology, is thus advisable for deciding between often contrasting mitochondrial and nuclear species scenarios in these organisms.     

Phylogeny and species delineation in the marine diatom Pseudo‐nitzschia (Bacillariophyta) using cox1, LSU,and ITS2 rRNA genes: A perspective in character evolution

Analyses of the mitochondrial cox1, the nuclear‐encoded large subunit (LSU), and the internal transcribed spacer 2 (ITS2) RNA coding region of Pseudo‐nitzschia revealed that the P. pseudodelicatissima complex can be phylogenetically grouped into three distinct clades (Groups I–III), while the P. delicatissima complex forms another distinct clade (Group IV) in both the LSU and ITS2 phylogenetic trees. It was elucidated that comprehensive taxon sampling (sampling of sequences), selection of appropriate target genes and outgroup, and alignment strategies influenced the phylogenetic accuracy. Based on the genetic divergence, ITS2 resulted in the most resolved trees, followed by cox1 and LSU. The morphological characters available for Pseudo‐nitzschia, although limited in number, were overall in agreement with the phylogenies when mapped onto the ITS2 tree. Information on the presence/absence of a central nodule, number of rows of poroids in each stria, and of sectors dividing the poroids mapped onto the ITS2 tree revealed the evolution of the recently diverged species. The morphologically based species complexes showed evolutionary relevance in agreement with molecular phylogeny inferred from ITS2 sequence–structure data. The data set of the hypervariable region of ITS2 improved the phylogenetic inference compared to the cox1 and LSU data sets. The taxonomic status of P. cuspidata and P. pseudodelicatissima requires further elucidation.     

Parallel plumage colour evolution and introgressive hybridization in wheatears

Genetic and phenotypic mosaics, in which various phenotypes and different genomic regions show discordant patterns of species or population divergence, offer unique opportunities to study the role of ancestral and introgressed genetic variation in phenotypic evolution. Here, we investigated the evolution of discordant phenotypic and genetic divergence in a monophyletic clade of four songbird taxa—pied wheatear (O. pleschanka), Cyprus wheatear (Oenanthe cypriaca), and western and eastern subspecies of black‐eared wheatear (O. h. hispanica and O. h. melanoleuca). Phenotypically, black back and neck sides distinguish pied and Cyprus wheatears from the white‐backed/necked black‐eared wheatears. Meanwhile, mitochondrial variation only distinguishes western black‐eared wheatear. In the absence of nuclear genetic data, and given frequent hybridization among eastern black‐eared and pied wheatear, it remains unclear whether introgression is responsible for discordance between mitochondrial divergence patterns and phenotypic similarities, or whether plumage coloration evolved in parallel. Multispecies coalescent analyses of about 20,000 SNPs obtained from RAD data mapped to a draft genome assembly resolve the species tree, provide evidence for the parallel evolution of colour phenotypes and establish western and eastern black‐eared wheatears as independent taxa that should be recognized as full species. The presence of the entire admixture spectrum in the Iranian hybrid zone and the detection of footprints of introgression from pied into eastern black‐eared wheatear beyond the hybrid zone despite strong geographic structure of ancestry proportions furthermore suggest a potential role for introgression in parallel plumage colour evolution. Our results support the importance of standing heterospecific and/or ancestral variation in phenotypic evolution.     

Multilocus phylogeny of the widely distributed South American lizard clade Eulaemus (Liolaemini,Liolaemus)

The lizard genus Liolaemus and different clades within it have been the focus of several recent phylogenetic studies mainly based on morphology and mtDNA. Although there is general consensus for recognizing two clades (subgenera) within the genus, [Liolaemus (sensu stricto) and Eulaemus], phylogenetic relationships within each subgenus remain difficult to elucidate, given incomplete taxonomic sampling and large discordance between published studies. Here, new phylogenetic relationships for the Eulaemus subgenus are proposed based on the largest molecular data set ever used for this clade, which includes 188 individuals and 14 loci representing different parts of the genome (mtDNA, anonymous nuclear loci and nuclear protein‐coding loci). This data set was analysed using two species tree approaches (*beast and MDC). Levels of discordance among methods were found, and with previously published studies, but results are robust enough to propose new phylogenetic hypotheses for the Eulaemus clade. Specifically well‐resolved and well‐supported novel hypotheses are provided within the lineomaculatus section, and we formally recognize the zullyae clade, the sarmientoi clade and the hatcheri group. We also resolve species relationships within the montanus section, and particularly within the melanops series. We found discordance between mitochondrial and nuclear trees and discussed alternative hypotheses for the lineomaculatus and montanus sections, as well as the challenge in resolving phylogenetic relationships for large clades in general.     

Estimating the temporal and spatial extent of gene flow among sympatric lizard populations (genus Sceloporus) in the southern Mexican highlands

Interspecific gene flow is pervasive throughout the tree of life. Although detecting gene flow between populations has been facilitated by new analytical approaches, determining the timing and geography of hybridization has remained difficult, particularly for historical gene flow. A geographically explicit phylogenetic approach is needed to determine the overlap of ancestral populations. In this study, we performed population genetic analyses, species delimitation, simulations and a recently developed approach of species tree diffusion to infer the phylogeographic history, timing and geographic extent of gene flow in lizards of the Sceloporus spinosus group. The two species in this group, S. spinosus and S. horridus, are distributed in eastern and western portions of Mexico, respectively, but populations of these species are sympatric in the southern Mexican highlands. We generated data consisting of three mitochondrial genes and eight nuclear loci for 148 and 68 individuals, respectively. We delimited six lineages in this group, but found strong evidence of mito‐nuclear discordance in sympatric populations of S. spinosus and S. horridus owing to mitochondrial introgression. We used coalescent simulations to differentiate ancestral gene flow from secondary contact, but found mixed support for these two models. Bayesian phylogeography indicated more than 60% range overlap between ancestral S. spinosus and S. horridus populations since the time of their divergence. Isolation–migration analyses, however, revealed near‐zero levels of gene flow between these ancestral populations. Interpreting results from both simulations and empirical data indicate that despite a long history of sympatry among these two species, gene flow in this group has only recently occurred.     

Range instability leads to cytonuclear discordance in a morphologically cryptic ground squirrel species complex

The processes responsible for cytonuclear discordance frequently remain unclear. Here, we employed an exon capture data set and demographic methods to test hypotheses generated by species distribution models to examine how contrasting histories of range stability vs. fluctuation have caused cytonuclear concordance and discordance in ground squirrel lineages from the Otospermophilus beecheyi species complex. Previous studies in O. beecheyi revealed three morphologically cryptic and highly divergent mitochondrial DNA lineages (named the Northern, Central and Southern lineages based on geography) with only the Northern lineage exhibiting concordant divergence for nuclear genes. Here, we showed that these mtDNA lineages likely formed in allopatry during the Pleistocene, but responded differentially to climatic changes that occurred since the last interglacial (~120,000 years ago). We find that the Northern lineage maintained a stable range throughout this period, correlating with genetic distinctiveness among all genetic markers and low migration rates with the other lineages. In contrast, our results suggested that the Southern lineage expanded from Baja California Sur during the Late Pleistocene to overlap and potentially swamp a contracting Central lineage. High rates of intraspecific gene flow between Southern lineage individuals among expansion origin and expansion edge populations largely eroded Central ancestry from autosomal markers. However, male‐biased dispersal in this system preserved signals of this past hybridization and introgression event in matrilineal‐biased X‐chromosome and mtDNA markers. Our results highlight the importance of range stability in maintaining the persistence of phylogeographic lineages, whereas unstable range dynamics can increase the tendency for lineages to merge upon secondary contact.     

Speciation processes in putative island endemic sister bat species: false impressions from mitochondrial DNA and microsatellite data

Cases of geographically restricted co‐occurring sister taxa are rare and may point to potential divergence with gene flow. The two bat species Murina gracilis and Murina recondita are both endemic to Taiwan and are putative sister species. To test for nonallopatric divergence and gene flow in these taxa, we generated sequences using Sanger and next‐generation sequencing, and combined these with microsatellite data for coalescent‐based analyses. MtDNA phylogenies supported the reciprocally monophyletic sister relationship between M. gracilis and M. recondita; however, clustering of microsatellite genotypes revealed several cases of species admixture suggesting possible introgression. Sequencing of microsatellite flanking regions revealed that admixture signatures stemmed from microsatellite allele homoplasy rather than recent introgressive hybridization, and also uncovered an unexpected sister relationship between M. recondita and the continental species Murina eleryi, to the exclusion of M. gracilis. To dissect the basis of these conflicts between ncDNA and mtDNA, we analysed sequences from 10 anonymous ncDNA loci with *beast and isolation‐with‐migration and found two distinct clades of M. eleryi, one of which was sister to M. recondita. We conclude that Taiwan was colonized by the ancestor of M. gracilis first, followed by the ancestor of M. recondita after a period of allopatric divergence. After colonization, the mitochondrial genome of M. recondita was replaced by that of the resident M. gracilis. This study illustrates how apparent signatures of sympatric divergence can arise from complex histories of allopatric divergence, colonization and hybridization, thus highlighting the need for rigorous analyses to distinguish between such scenarios.     

Insight into the population structure of hardhead silverside,Atherinomorus stipes (Teleostei: Atherinidae), in Belize and the Florida Keys using nd2

Little is known about the natural history, biology, and population genetic structure of the Hardhead Silverside, Atherinomorus stipes, a small schooling fish found around islands throughout the Caribbean. Our field observations of A. stipes in the cays of Belize and the Florida Keys found that populations tend to be in close association with the shoreline in mangrove habitats. Due to this potential island‐based population structuring, A. stipes represents an ideal system to examine questions about gene flow and isolation by distance at different geographic scales. For this study, the mitochondrial gene nd2 was amplified from 394 individuals collected from seven different Belizean Cays (N = 175) and eight different Floridian Keys (N = 219). Results show surprisingly high haplotype diversity both within and between island‐groups, as well as a high prevalence of unique haplotypes within each island population. The results are consistent with models that require gene flow among populations as well as in situ evolution of rare haplotypes. There was no evidence for an isolation by distance model. The nd2 gene tree consists of two well‐supported monophyletic groups: a Belizean‐type clade and a Floridian‐type clade, indicating potential species‐level differentiation.