Nonallopatric and parallel origin of local reproductive barriers between two snail ecotypes

Theory suggests that speciation is possible without physical isolation of populations (hereafter, nonallopatric speciation), but recent nonallopatric models need the support of irrefutable empirical examples. We collected snails (Littorina saxatilis) from three areas on the NW coast of Spain to investigate the population genetic structure of two ecotypes. Earlier studies suggest that these ecotypes may represent incipient species: a large, thick-shelled 'RB' ecotype living among the barnacles in the upper intertidal zone and a small, thin-shelled 'SU' ecotype living among the mussels in the lower intertidal zone only 10-30 m away. The two ecotypes overlap and hybridize in a midshore zone only 1-3 m wide. Three different types of molecular markers [allozymes, mitochondrial DNA (mtDNA) and microsatellites] consistently indicated partial reproductive isolation between the RB and the SU ecotypes at a particular site. However, each ecotype was related more closely to the other ecotype from the same site than to the same ecotype from another site further along the Galician coast (25-77 km away). These findings supported earlier results based solely on allozyme variation and we could now reject the possibility that selection produced these patterns. The patterns of genetic variation supported a nonallopatric model in which the ecotypes are formed independently at each site by parallel evolution and where the reproductive barriers are a byproduct of divergent selection for body size. We argue that neither our laboratory hybridization experiments nor our molecular data are compatible with a model based on allopatric ecotype formation, secondary overlap and introgression.     

Lack of early laboratory postzygotic reproductive isolation between two ecotypes of Littorina saxatilis (Mollusca, Gastropoda) showing strong premating sexual isolation

Two sympatric and divergent adaptive ecotypes of Littorina saxatilis (RB and SU) are known to hybridize showing partial premating isolation in the wild. Previous studies have revealed that morphological intermediate forms (presumably hybrids) present fitness (viability, sexual selection and fecundity) similar to that from pure ecotypes at the mid-shore. However, the absence of postzygotic isolation due to genetic incompatibility cannot be ruled out unless it is measured directly on true F ₁ hybrids. In this study, we overcome this problem and present data on 56 individual crosses including the four possible mating combinations (RB/RB, RB/SU, SU/RB and SU/SU) to compare fertilization and fecundity rates (including young progeny viability) between the four type crosses. Pooled RB female crosses showed apparently larger fertility and fecundity than pooled SU female crosses, probably because of differences in fecundity and laboratory survivorship between ecotypes. However, similar fertilization and fecundity rates were found for both RB and SU females when mated with different male types, supporting the idea that genetic-incompatibility-based postzygotic isolation can be ignored as a major determinant of this polymorphism in nature.     

Selection on hybrids of ecologically divergent ecotypes of a marine snail: the relative importance of exogenous and endogenous barriers

Unravelling the form of selection acting on hybrids of ecotypes undergoing ecological speciation is essential to understand the mechanisms behind the evolution of reproductive isolation in the face of gene flow. Shell phenotype is known to be affected by natural selection and is involved in the fitness of the marine snail Littorina saxatilis. Here, we studied the association between shell traits and fitness in hybrids in order to determine the relative role of exogenous and endogenous selection in this hybrid zone of L. saxatilis. We show that directional selection is the predominant mode of selection among hybrids. We also show its heterogeneity, affecting different shell traits, within populations at the level of the microhabitat. Therefore, endogenous selection mechanisms are most probably lacking in this hybrid zone and exogenous barriers (pre‐ and post‐zygotic) are possibly one of the main forces behind the evolution of barriers to gene flow between these ecologically divergent ecotypes. This study shows how this barrier might represent an important type of reproductive isolation within ecological speciation, and this should be taken into account in future studies of speciation in hybrid zones. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 391–400.     

Genetic differentiation and estimation of effective population size and migration rates in two sympatric ecotypes of the marine snail Littorina saxatilis

On exposed rocky shores in Galicia (northwest Spain), a striking polymorphism exists between two ecotypes (RB and SU) of Littorina saxatilis that occupy different levels of the intertidal zone and exhibit an incomplete reproductive isolation. The setting has been suggested to represent ongoing sympatric speciation by ecological adaptation of the two ecotypes to their respective habitats. In this article we address whether or not the ecotypes have developed their own population structures in response to the rigors of their corresponding environments and life histories. We analyzed four to five allozymic loci from three surveys of the same sites, spanning a 14-year period. An experimental design including three localities with two transects per locality and three shore levels allowed studying temporal and spatial population structure and estimation of effective population sizes (N(e)), neighborhood sizes (N(n)), and migration rates (m). Genetic differentiation was significantly lower in RB populations (theta(ST) = 0.067) than in SU ones (theta(ST) = 0.124). Mean estimates of N(e), N(n), and m did not differ significantly between ecotypes, but local ecotype differences in migration between the two closest localities (larger migration rates in RB than in SU populations) could explain the pattern in population differentiation.     

Habitat preference and flowering‐time variation contribute to reproductive isolation between diploid and autotetraploid Anacamptis pyramidalis

Tetraploid lineages are typically reproductively isolated from their diploid ancestors by post‐zygotic isolation via triploid sterility. Nevertheless, polyploids often also exhibit ecological divergence that could contribute to reproductive isolation from diploid ancestors. In this study, we disentangled the contribution of different forms of reproductive isolation between sympatric diploid and autotetraploid individuals of the food‐deceptive orchid Anacamptis pyramidalis by quantifying the strength of seven reproductive barriers: three prepollination, one post‐pollination prezygotic and three post‐zygotic. The overall reproductive isolation between the two cytotypes was found very high, with a preponderant contribution of two prepollination barriers, that is phenological and microhabitat differences. Although the contribution of post‐zygotic isolation (triploid sterility) is confirmed in our study, these results highlight that prepollination isolation, not necessarily involving pollinator preference, can represent a strong component of reproductive isolation between different cytotypes. Thus, in the context of polyploidy as quantum speciation, that generates reproductive isolation via triploid sterility, ecological divergence can strengthen the reproductive isolation between cytotypes, reducing the waste of gametes in low fitness interploidy crosses and thus favouring the initial establishment of the polyploid lineage. Under this light, speciation by polyploidy involves ecological processes and should not be strictly considered as a nonecological form of speciation.     

Case studies and mathematical models of ecological speciation. 3: Ecotype formation in a Swedish snail

Formation of partially reproductively isolated ecotypes in the rough periwinkle, Littorina saxatilis , may be a case of incipient nonallopatric ecological speciation. To better understand the dynamics of ecotype formation, its timescale, driving forces and evolutionary consequences, we developed a spatially explicit, individual-based model incorporating relevant ecological, spatial and mate selection data for Swedish L. saxatilis . We explore the impact of bounded hybrid superiority, ecological scenarios and mate selection systems on ecotype formation, gene flow and the evolution of prezygotic isolation. Our model shows that ecotypes are expected to form rapidly in parapatry under conditions applicable to Swedish L. saxatilis and may proceed to speciation. However, evolution of nonrandom mating had complex behaviour. Ecotype evolution was inhibited by pre-existing mating preferences, but facilitated by the evolution of novel preferences. While in many scenarios positive assortative mating reduced gene flow between ecotypes, in others negative assortative mating arose, preferences were lost after ecotype formation, preferences were confined to one ecotype or the ancestral ecotype became extinct through sexual selection. Bounded hybrid superiority (as observed in nature) enhanced ecotype formation but increased gene flow. Our results highlight that ecotype formation and speciation are distinct processes: factors that contribute to ecotype formation can be detrimental to speciation and vice versa. The complex interactions observed between local adaptation and nonrandom mating imply that generalization from data is unreliable without quantitative theory for speciation.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

Incipient speciation in the Anastrepha fraterculus cryptic species complex: reproductive compatibility between A. sp.1 aff. fraterculus and A. sp.3 aff. fraterculus

Reproductive isolation among biological strains can be detected by analyzing reproductive barriers between populations. The Anastrepha fraterculus Wiedemann (Diptera: Tephritidae) cryptic species complex comprises eight morphotypes defined by morphometric characteristics. The present study analyzed the existence of pre‐ and post‐zygotic isolation among populations of A. sp.1 aff. fraterculus (A. sp.1) and A. sp.3 aff. fraterculus (A. sp.3) determined by tests of sexual compatibility, analyses of the temporal pattern of calling behavior, and the reproductive performance of the offspring of homotypic and heterotypic crosses. Pre‐zygotic tests indicated asymmetric matings with preferences for homotypic crosses. The A. sp.3 population demonstrated distinct characteristics, including low copulation percentages and a high proportion of homotypic copulations under crowded conditions. The females of the A. sp.1 population demonstrated lower copulation latency. The analyzed morphotypes demonstrated two peaks in pheromone emission activity, at the beginning and at the end of the morning and were differentiated in the abundance of males. The calling behavior of the hybrids between A. sp.1 male and A. sp.3 female crosses had an intermediate pattern compared to either parental type. Post‐zygotic compatibility tests indicated low viability of the cross between A. sp.3 males and A. sp.1 females. Pronounced asymmetries were found in the sex ratios of the offspring of heterotypic crosses, and only hybrids from A. sp.1 male vs. A. sp.3 female crosses produced descendants. The observed pre‐ and post‐zygotic incompatibilities reflect the probable occurrence of incipient speciation between A. sp.1 and A. sp.3. A more adequate taxonomic classification of this species complex that considers the distinct characteristics of each morphotype will be needed to improve environmentally wise control methods against this insect pest.     

Proteomic evidence of a paedomorphic evolutionary process within a marine snail species: a strategy for adapting to extreme ecological conditions?

The exposed and sheltered ecotypes of the marine snail Littorina saxatilis from European rocky shores are considered a key model system to study adaptation and ecological speciation. Previous studies showed that two ecotypes (RB and SU) of this species in NW Spain have adapted differently to different shore levels and microhabitats. In order to understand how this divergent adaptive process has been accomplished, we followed a quantitative proteomic approach to investigate the proteome variation in a number of different biological factors, that is, ecotype, ontogeny and their interactions. This approach allowed testing the hypothesis that one of the ecotypes has evolved by paedomorphosis, and also whether or not the molecular mechanisms related to ecotype differentiation are set up in early developmental stages. Additionally, the identification of some candidate proteins using mass spectrometry provides some functional insights into these evolutionary processes. Results from this study provided evidence of higher ontogenetic differentiation at proteome level in the RB (metamorphic) than in SU (paedomorphic) ecotype that point to the possibility of juvenile stage retention in this latter ecotype. The level of protein expression (proteome) differences between ecotypes maintained nearly constant from late embryonic stages to adulthood, although some proteins involved in these changes considerably differed in embryonic compared to other ontogenetic stages. Paedomorphosis may be the evolutionary response of the SU ecotype of solving the trade‐off during sexually immaturity that is caused by the evolution of small size arising from adaptation to the wave‐exposed habitat. Some potential candidate genes of adaptation related to energetic metabolism have been identified, providing a promising baseline for future functional analyses.     

Insights into the role of differential gene expression on the ecological adaptation of the snail <Emphasis Type="Italic">Littorina saxatilis</Emphasis>

Background  

In the past 40 years, there has been increasing acceptance that variation in levels of gene expression represents a major source of evolutionary novelty. Gene expression divergence is therefore likely to be involved in the emergence of incipient species, namely, in a context of adaptive radiation. In this study, a genome-wide expression profiling approach (cDNA-AFLP), validated by quantitative real-time polymerase chain reaction (qPCR) were used to get insights into the role of differential gene expression on the ecological adaptation of the marine snail Littorina saxatilis. This gastropod displays two sympatric ecotypes (RB and SU) which are becoming one of the best studied systems for ecological speciation.     

Proteomic comparison between two marine snail ecotypes reveals details about the biochemistry of adaptation

The proteomic changes occurring during speciation are fundamental to understand this process, though they have been rarely addressed until present. Therefore, we compared the proteome of two ecotypes (RB and SU) of the marine snail Littorina saxatilis, a case of sympatric incomplete speciation, originated as a byproduct of adaptation to distinct habitats. Thus, the RB ecotype is able to resist stresses of desiccation and temperature on the upper shore, whereas the SU ecotype defies strong physical disturbances due to wave action. Qualitative analyses of 2-DE gels demonstrated 21 proteins differentially expressed (1.4% of the proteome, 1.2% after considering type-I errors), while quantitative changes accounted for differences in 22 spots (16% of the proteome, 11% after considering type-I errors). These results suggest that adaptative phenotypic plasticity, natural selection, or both maintain these ecotypes in sympatry. Among the proteins identified by MS, we found that fructose-bisphosphate aldolase and arginine kinase were up-regulated in the SU ecotype, suggesting an enhancement of the level of energy available as ATP, in order to withstand its wave-exposed habitat.     

Immigrant and extrinsic hybrid seed inviability contribute to reproductive isolation between forest and dune ecotypes of Epipactis helleborine (Orchidaceae)

Reproductive isolation caused by divergent natural selection arising from differences between ecological environments or ecological interactions represents a major mechanism contributing to speciation, but its relative importance is poorly known. In this study, controlled reciprocal crossings and seed germination experiments were combined with genetic and morphometric analyses to test the hypothesis that previously described differences in mycorrhizal communities between forest and dune ecotypes of Epipactis helleborine were sufficiently strong to create complete reproductive isolation between the two ecotypes. Molecular analyses using 770 SNP markers showed that the two ecotypes were genetically distinct and that populations of the dune ecotype were genetically impoverished compared to populations of the forest ecotype. Morphologically, the two ecotypes were also significantly different, with plants of the dune ecotype generally being smaller than plants of the forest ecotype. The results further showed that immigrant seeds had a significantly lower probability of protocorm formation than native seeds, indicating strong immigrant inviability. Although both ecotypes were able to cross easily and to produce a large number of viable seeds, hybrid seeds showed significantly lower protocorm formation than pure seeds, further contributing to reproductive isolation. Overall, these results indicate that interfertile populations of a widespread orchid adapting to contrasting environments diverge as a consequence of concurrent selection acting against immigrants and hybrids and suggest that mycorrhizal fungi can play a role in the early stages of plant speciation.     

Intergenomic epistasis causes asynchronous hatch times in whitefish hybrids,but only when parental ecotypes differ

Support for the theory of ecological speciation requires evidence for ecological divergence between species which directly or indirectly causes reproductive isolation. This study investigates effects of ecological vs. genetic disparity of parental species on the presence of endogenous selection (deformation and mortality rates) and potential sources of exogenous selection (growth rates and hatch timing) on hybrids. Hybrid embryonic development is analysed in a common‐garden full‐sib cross of three species belonging to two different ecotypes within the Coregonus lavaretus species flock in the central Alpine region of Europe. Although hatch timing was similar across the three species, embryonic growth rates and egg sizes differed between ecotypes. This led to a mismatch between embryonic growth rate and egg size in hybrid crosses that reveals epistasis between the maternal and embryonic genomes and transgressive hatch times that were asynchronous with control crosses. A strong constraint of egg size to embryo size at late development was also evident. We argue that this demonstrates potential for coadaptation of a maternal trait (egg size) with offspring growth rate to be an important source of selection against hybridization between ecotypes with different egg sizes. Implications for the measurement and quantification of early life‐history traits affected by this additive relationship, such as hatch day and larval size, are also discussed.     

WIDESPREAD HOST‐DEPENDENT HYBRID UNFITNESS IN THE PEA APHID SPECIES COMPLEX

Linking adaptive divergence to hybrid unfitness is necessary to understand the ecological factors contributing to reproductive isolation and speciation. To date, this link has been demonstrated in few model systems, most of which encompass ecotypes that occupy relatively early stages in the speciation process. Here we extend these studies by assessing how host‐plant adaptation conditions hybrid fitness in the pea aphid, Acyrthosiphon pisum. We made crosses between and within five pea aphid biotypes adapted to different host plants and representing various stages of divergence within the complex. Performance of F1 hybrids and nonhybrids was assessed on a “universal” host that is favorable to all pea aphid biotypes in laboratory conditions. Although hybrids performed equally well as nonhybrids on the universal host, their performance was much lower than nonhybrids on the natural hosts of their parental populations. Hence, hybrids, rather than being intrinsically deficient, are maladapted to their parents’ hosts. Interestingly, the impact of this maladaptation was stronger in certain hybrids from crosses involving the most divergent biotype, suggesting that host‐dependent postzygotic isolation has continued to evolve late in divergence. Even though host‐independent deficiencies are not excluded, hybrid maladaptation to parental hosts supports the hypothesis of ecological speciation in this complex.     

Mate search and aggregation behaviour in the Galician hybrid zone of Littorina saxatilis

In Galician rocky shores two ecotypes of the snail L. saxatilis can be found in sympatry. A ridged and banded ecotype (RB-morph) and a smooth and unbanded ecotype (SU-morph) overlap in midshore with the production of some hybrids. The distinct morphs mate assortatively and there is evidence of a partial reproductive barrier between them. This sexual isolation is caused by a nonrandom microdistribution and mate choice behaviour. Mucus trail-following, movement rate and aggregation behaviour were studied to determine their roles in the mating behaviour and sexual isolation of this species. Morph-specific mucus trail-following could not, in our experiments, explain either of these two processes. The reasons for the aggregation of morphs were investigated by Monte Carlo simulations of data from natural populations, which showed that size aggregation (refuge sizes fit different sized morphs differently) could explain only about 36% of the morph aggregation in adult snails. In the laboratory, morph aggregation was still present, and simulations suggested that size aggregation was the possible explanation. Thus, morph aggregation in Galician L. saxatilis has to be explained also by other causes in addition to size aggregation. These may be a combination of contrasting preferences for barnacle and mussel patches in the two morphs, and possibly longer copulation and pair formation time with similar sized snails of the same morph. Thus aggregation behaviour, but not trail-following, contributes to incipient reproductive isolation and perhaps sympatric speciation in Galician L. saxatilis populations.     

Post‐zygotic hybrid viability in sympatric species pairs: a case study from European lampreys

Ecological speciation mechanisms are widely assumed to play an important role in the early stages of divergence between incipient species, and this especially true of fishes. In the present study, we tested explicitly for post‐zygotic barriers to gene flow between a sympatric, recently diverged lamprey species pair that likely arose through ecological divergence. Experimental in vitro hybridization between anadromous parasitic Lampetra fluviatilis and resident nonparasitic Lampetra planeri resulted in a high proportion of embryos capable of attaining the burrowing pro‐larval stage, strongly indicating no post‐zygotic barriers to gene flow between these species. A sympatric, locally‐adapted resident parasitic form of L. fluviatilis was also found to successfully hybridize with both members of this species pair. The consequences of these findings are discussed in the context of lamprey speciation. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 108 , 378–383.     

The evolution of intrinsic reproductive isolation in the genus Cakile (Brassicaceae)

In theory, adaptive divergence can increase intrinsic post‐zygotic reproductive isolation (RI), either directly via selection on loci associated with RI, or indirectly via linkage of incompatibility loci with loci under selection. To test this hypothesis, we measured RI at five intrinsic post‐zygotic reproductive barriers between 18 taxa from the genera Cakile and Erucaria (Brassicaceae). Using a comparative framework, we tested whether the magnitude of RI was associated with genetic distance, geographic distance, ecological divergence and parental mating system. Early stages of post‐zygotic RI related to F₁ viability (i.e. initial seed set) tended to be stronger than later stages related to F₁ fecundity (i.e. flower number, fruit number). Mating system significantly influenced early stages of RI, such that RI was lowest when the mother was selfing and father was outcrossing, consistent with an imbalance between sink strength and resistance to provisioning. We found little evidence that adaptive divergence accelerates the evolution of intrinsic post‐zygotic RI, consistent with a nonecological model of evolution that predicts the nonlinear accumulation of RI and RI asymmetry with time (i.e. genetic distance), irrespective of adaptive divergence. Thus, although certain aspects of ecological divergence do not appear to have contributed strongly to the evolution of RI in this system, divergence in mating system actually reduced RI, suggesting that mating system evolution may play a significant role in speciation dynamics.     

Is embryo abortion a post‐zygotic barrier to gene flow between Littorina ecotypes?

Genetic incompatibilities contribute to reproductive isolation between many diverging populations, but it is still unclear to what extent they play a role if divergence happens with gene flow. In contact zones between the "Crab" and "Wave" ecotypes of the snail Littorina saxatilis, divergent selection forms strong barriers to gene flow, while the role of post‐zygotic barriers due to selection against hybrids remains unclear. High embryo abortion rates in this species could indicate the presence of such barriers. Post‐zygotic barriers might include genetic incompatibilities (e.g. Dobzhansky–Muller incompatibilities) but also maladaptation, both expected to be most pronounced in contact zones. In addition, embryo abortion might reflect physiological stress on females and embryos independent of any genetic stress. We examined all embryos of >500 females sampled outside and inside contact zones of three populations in Sweden. Females' clutch size ranged from 0 to 1,011 embryos (mean 130 ± 123), and abortion rates varied between 0% and 100% (mean 12%). We described female genotypes by using a hybrid index based on hundreds of SNPs differentiated between ecotypes with which we characterized female genotypes. We also calculated female SNP heterozygosity and inversion karyotype. Clutch size did not vary with female hybrid index, and abortion rates were only weakly related to hybrid index in two sites but not at all in a third site. No additional variation in abortion rate was explained by female SNP heterozygosity, but increased female inversion heterozygosity added slightly to increased abortion. Our results show only weak and probably biologically insignificant post‐zygotic barriers contributing to ecotype divergence, and the high and variable abortion rates were marginally, if at all, explained by hybrid index of females.     

Postglacial climate changes and rise of three ecotypes of harbour porpoises,Phocoena phocoena,in western Palearctic waters

Despite no obvious barriers to gene flow in the marine realm, environmental variation and ecological specializations can lead to genetic differentiation in highly mobile predators. Here, we investigated the genetic structure of the harbour porpoise over the entire species distribution range in western Palearctic waters. Combined analyses of 10 microsatellite loci and a 5085 base‐pair portion of the mitochondrial genome revealed the existence of three ecotypes, equally divergent at the mitochondrial genome, distributed in the Black Sea (BS), the European continental shelf waters, and a previously overlooked ecotype in the upwelling zones of Iberia and Mauritania. Historical demographic inferences using approximate Bayesian computation (ABC) suggest that these ecotypes diverged during the last glacial maximum (c. 23–19 kilo‐years ago, kyrbp ). ABC supports the hypothesis that the BS and upwelling ecotypes share a more recent common ancestor (c. 14 kyrbp ) than either does with the European continental shelf ecotype (c. 28 kyrbp ), suggesting they probably descended from the extinct populations that once inhabited the Mediterranean during the glacial and post‐glacial period. We showed that the two Atlantic ecotypes established a narrow admixture zone in the Bay of Biscay during the last millennium, with highly asymmetric gene flow. This study highlights the impacts that climate change may have on the distribution and speciation process in pelagic predators and shows that allopatric divergence can occur in these highly mobile species and be a source of genetic diversity.     

Ongoing speciation within the Anastrepha fraterculus cryptic species complex: the case of the Andean morphotype

The Anastrepha fraterculus (Wiedemann) (Diptera: Tephritidae) cryptic species complex is currently composed of seven taxonomically recognized morphotypes. Both, pre‐ and post‐zygotic isolation has been documented among four of these morphotypes, revealing that in fact they appear to be distinct biological entities. In order to progress in the full delimitation of species within the complex, we examined reproductive isolation between a Colombian population of the Andean morphotype and populations belonging to four other morphotypes spanning from Mexico to Argentina. Flies from the Andean morphotype exhibited strong pre‐zygotic mating isolation through temporal partitioning of mating activity. Post‐zygotic isolation was observed for crosses of males of all morphotypes and Andean morphotype females, yet most of the F1 hybrid ♂ × F1 hybrid ♀ self‐crosses showed normal levels of fertility, a finding suggesting a nuclear–cytoplasmic interaction according to previous studies. Overall, the Andean morphotype within the complex also appears to be a distinct biological entity. We discuss the implications of these findings for the understanding of speciation mechanisms in the Neotropical genus Anastrepha.