Increased genetic divergence between two closely related fir species in areas of range overlap

Because of introgressive hybridization, closely related species can be more similar to each other in areas of range overlap (parapatry or sympatry) than in areas where they are geographically isolated from each other (allopatry). Here, we report the reverse situation based on nuclear genetic divergence between two fir species, Abies chensiensis and Abies fargesii, in China, at sites where they are parapatric relative to where they are allopatric. We examined genetic divergence across 126 amplified fragment length polymorphism (AFLP) markers in a set of 172 individuals sampled from both allopatric and parapatric populations of the two species. Our analyses demonstrated that AFLP divergence was much greater between the species when comparisons were made between parapatric populations than between allopatric populations. We suggest that selection in parapatry may have largely contributed to this increased divergence.     

The genetics and ecology of sympatric speciation: A case study

Mating occurs on the larval host plant in allRhagoletis species (Diptera: Tephritidae). We show how this attribute, when coupled with certain differences in other biological traits, strongly influences the mode of speciation. In species of thesuavis species group, host shifts have never occurred during speciation, and larvae feed in the husks of any walnut species(Juglans spp.), which are highly toxic. Taxa are allopatric or parapatric and exhibit deep phylogenetic nodes suggesting relatively ancient speciation events. Traits responsible for species and mate recognition, particularly in parapatric species, are morphologically distinct and strongly sexually dimorphic. All aspects of their biology, genetics and distribution are consistent with a slow rate of allopatric speciation followed by morphological divergence in secondary contact. In contrast, speciation in thepomonella species group has always involved a shift to a new, usually unrelated, non-toxic host, and all taxa within these groups are sympatric, monophagous and morphologically indistinguishable from one another. Phylogenetic nodes are very shallow, indicating recent sympatric speciation. Sympatric divergence is promoted by genetic variation which allows a portion of the original species to shift to a new habitat or host. Evidence suggests that changes in a few key loci responsible for host selection and fitness on a new host may initiate host shifts. By exploiting different habitats, competition for resources between diverging populations is reduced or avoided. We provide evidence that in phytophagous and parasitic insects sufficient intrinsic barriers to gene flow can evolve between sister populations as they adapt to different habitats or hosts to allow each population to establish independent evolutionary lineages in sympatry.     

Parapatric divergence of sympatric morphs in a salamander: incipient speciation on Long Island?

Speciation is often categorized based on geographic modes (allopatric, parapatric or sympatric). Although it is widely accepted that species can arise in allopatry and then later become sympatrically or parapatrically distributed, patterns in the opposite direction are also theoretically possible (e.g. sympatric lineages or ecotypes becoming parapatric), but such patterns have not been shown at a macrogeographic scale. Here, we analyse genetic, climatic, ecological and morphological data and show that two typically sympatric colour morphs of the salamander Plethodon cinereus (redback and leadback) appear to have become parapatrically distributed on Long Island, New York, with pure‐redback populations in the west and pure‐leadback populations in the east (and polymorphic populations in between and on the mainland). In addition, the pure‐leadback populations in eastern Long Island are genetically, ecologically and morphologically divergent from both mainland and other Long Island populations, suggesting the possibility of incipient speciation. This parapatric separation seems to be related to the different ecological preferences of the two morphs, preferences which are present on the mainland and across Long Island. These results potentially support the idea that spatial segregation of sympatric ecotypes may sometimes play an important part in parapatric speciation.     

Along the speciation continuum in sticklebacks

Speciation can be viewed as a continuum, potentially divisible into several states: (1) continuous variation within panmictic populations, (2) partially discontinuous variation with minor reproductive isolation, (3) strongly discontinuous variation with strong but reversible reproductive isolation and (4) complete and irreversible reproductive isolation. Research on sticklebacks (Gasterosteidae) reveals factors that influence progress back and forth along this continuum, as well as transitions between the states. Most populations exist in state 1, even though some of these show evidence of disruptive selection and positive assortative mating. Transitions to state 2 seem to usually involve strong divergent selection coupled with at least a bit of geographic separation, such as parapatry (e.g. lake and stream pairs and mud and lava pairs) or allopatry (e.g. different lakes). Transitions to state 3 can occur when allopatric or parapatric populations that evolved under strong divergent selection come into secondary contact (most obviously the sympatric benthic and limnetic pairs), but might also occur between populations that remained in parapatry or allopatry. Transitions to state 4 might be decoupled from these selective processes, because the known situations of complete, or nearly complete, reproductive isolation (Japan Sea and Pacific Ocean pair and the recognized gasterosteid species) are always associated with chromosomal rearrangements and environment‐independent genetic incompatibilities. Research on sticklebacks has thus revealed complex and shifting interactions between selection, adaptation, mutation and geography during the course of speciation.     

ASYMMETRIC RESPONSE PATTERNS TO SUBSPECIES‐SPECIFIC SONG DIFFERENCES IN ALLOPATRY AND PARAPATRY IN THE GRAY‐BREASTED WOOD‐WREN

Song divergence between closely related taxa may play a critical role in the evolutionary processes of speciation and hybridization. We explored song variation between two Ecuadorian subspecies of the gray‐breasted wood‐wren (Henicorhina leucophrys) and tested the impact of song divergence on response behaviors. Songs were significantly different between the two subspecies, even between two parapatric populations 10 km apart. Playback experiments revealed an asymmetric response pattern to these divergent subspecies specific songs; one subspecies responded more to songs of its own subspecies than to the other subspecies’ songs, whereas the second responded equally strongly to songs of both subspecies. While song parameters revealed a mixed pattern of divergence between allopatric and parapatric populations, the majority of spectral characteristics showed increased divergence in parapatry, suggestive of character displacement. This increased song divergence in parapatry appeared to affect behavioral responses to playback as discriminating responses were most prominent in parapatry and against parapatric songs. The clear behavioral impact of subspecies‐specific song differences supports a potential role for song as an acoustic barrier to gene flow. The asymmetric nature of the responses suggests that song divergence could affect the direction of gene flow and the position of the subspecies‐specific transition.     

The evolution of host preference in allopatric vs. parapatric populations of Timema cristinae walking-sticks

Divergent habitat preferences can contribute to speciation, as has been observed for host-plant preferences in phytophagous insects. Geographic variation in host preference can provide insight into the causes of preference evolution. For example, selection against maladaptive host-switching occurs only when multiple hosts are available in the local environment and can result in greater divergence in regions with multiple vs. a single host. Conversely, costs of finding a suitable host can select for preference even in populations using a single host. Some populations of Timema cristinae occur in regions with only one host-plant species present (in allopatry, surrounded by unsuitable hosts) whereas others occur in regions with two host-plant species adjacent to one another (in parapatry). Here, we use host choice and reciprocal-rearing experiments to document genetic divergence in host preference among 33 populations of T. cristinae. Populations feeding on Ceanothus exhibited a stronger preference for Ceanothus than did populations feeding on Adenostoma. Both allopatric and parapatric pairs of populations using the different hosts exhibited divergent host preferences, but the degree of divergence tended to be greater between allopatric pairs. Thus, gene flow between parapatric populations apparently constrains divergence. Host preferences led to levels of premating isolation between populations using alternate hosts that were comparable in magnitude to previously documented premating isolation caused by natural and sexual selection against migrants between hosts. Our findings demonstrate how gene flow and different forms of selection interact to determine the magnitude of reproductive isolation observed in nature.     

Testing the role of climate in speciation: New methods and applications to squamate reptiles (lizards and snakes)

Climate may play important roles in speciation, such as causing the range fragmentation that underlies allopatric speciation (through niche conservatism) or driving divergence of parapatric populations along climatic gradients (through niche divergence). Here, we developed new methods to test the frequency of climate niche conservatism and divergence in speciation, and applied it to species pairs of squamate reptiles (lizards and snakes). We used a large‐scale phylogeny to identify 242 sister species pairs for analysis. From these, we selected all terrestrial allopatric pairs with sufficient occurrence records (n = 49 pairs) and inferred whether each originated via climatic niche conservatism or climatic niche divergence. Among the 242 pairs, allopatric pairs were most common (41.3%), rather than parapatric (19.4%), partially sympatric (17.7%), or fully sympatric species pairs (21.5%). Among the 49 selected allopatric pairs, most appeared to have originated via climatic niche divergence (61–76%, depending on the details of the methods). Surprisingly, we found greater climatic niche divergence between allopatric sister species than between parapatric pairs, even after correcting for geographic distance. We also found that niche divergence did not increase with time, further implicating niche divergence in driving lineage splitting. Overall, our results suggest that climatic niche divergence may often play an important role in allopatric speciation, and the methodology developed here can be used to address the generality of these findings in other organisms.     

Hawthorn-infesting populations of Rhagoletis pomonella in Mexico and speciation mode plurality

Categorizing speciation into dichotomous allopatric versus nonallopatric modes may not always adequately describe the geographic context of divergence for taxa. If some of the genetic changes generating inherent barriers to gene flow between populations evolved in geographic isolation, whereas others arose in sympatry, then the mode of divergence would be mixed. The apple maggot fly, Rhagoletis pomonella, has contributed to this emerging concept of a mixed speciation mode "plurality." Genetic studies have implied that a source of diapause life-history variation associated with inversions and contributing to sympatric host race formation and speciation for R. pomonella in the United States may have introgressed from the Eje Volcanico Trans Mexicano (EVTM; a.k.a. the Altiplano) in the past. A critical unresolved issue concerning the introgression hypothesis is how past gene flow occurred given the current 1200-km disjunction in the ranges of hawthorn-infesting flies in the EVTM region of Mexico and the southern extreme of the U.S. population in Texas. Here, we report the discovery of a hawthorn-infesting population of R. pomonella in the Sierra Madre Oriental Mountains (SMO) of Mexico. Sequence data from 15 nuclear loci and mitochondrial DNA imply that the SMO flies are related to, but still different from, U.S. and EVTM flies. The host affiliations, diapause characteristics, and phylogeography of the SMO population are consistent with it having served as a conduit for gene flow between Mexico and the United States. We also present evidence suggesting greater permeability of collinear versus rearranged regions of the genome to introgression, in accord with recent models of chromosomal speciation. We discuss the implications of the results in the context of speciation mode plurality. We do not argue for abandoning the terms sympatry or allopatry, but caution that categorizing divergence into either/or geographic modes may not describe the genetic origins of all species. For R. pomonella in the United States, the proximate selection pressures triggering race formation and speciation stem from sympatric host shifts. However, some of the phenological variation contributing to host-related ecological adaptation and reproductive isolation in sympatry at the present time appears to have an older history, having originated and become packaged into inversion polymorphism in allopatry.     

ALLOPATRIC ORIGIN OF SYMPATRIC POPULATIONS OF LAKE WHITEFISH (COREGONUS CLUPEAFORMIS) AS REVEALED BY MITOCHONDRIAL-DNA RESTRICTION ANALYSIS

In the paper, restriction-fragment length polymorphisms in mitochondrial DNA (mtDNA) were studied to test the hypothesis that sympatric populations of lake whitefish in the Allegash basin have recently diverged through sympatric speciation. Thirteen restriction enzymes were used to analyze mtDNA of 156 specimens representing 13 populations from eastern Canada and northern Maine where normal and dwarf phenotypes of whitefish exist in sympatry and allopatry. Two monophyletic assemblages of populations that exhibit different geographic distributions were identified. One showed an eastern distribution that expands from Cape Breton to the Allegash basin and the other exhibits a more western distribution. The Allegash basin was the only area of overlap. The western assemblage exhibited the normal size phenotype in all cases, whereas the eastern assemblage exhibited the normal size phenotype in allopatric conditions and the dwarf size phenotype in sympatry. The existence of sympatric pairs in the Allegash basin result from the secondary contact of two monophyletic groups of whitefish that evolved allopatrically in separate refugia during the last glaciation events. The weak mtDNA difference of sympatric pairs suggests that speciation of lake whitefish in eastern North America was accompanied by only minor alterations of the ancestral gene pool.     

THE RATE TEST OF SPECIATION: ESTIMATING THE LIKELIHOOD OF NON‐ALLOPATRIC SPECIATION FROM REPRODUCTIVE ISOLATION RATES IN DROSOPHILA

Among the most debated subjects in speciation is the question of its mode. Although allopatric (geographical) speciation is assumed the null model, the importance of parapatric and sympatric speciation is extremely difficult to assess and remains controversial. Here I develop a novel approach to distinguish these modes of speciation by studying the evolution of reproductive isolation (RI) among taxa. I focus on the Drosophila genus, for which measures of RI are known. First, I incorporate RI into age‐range correlations. Plots show that almost all cases of weak RI are between allopatric taxa whereas sympatric taxa have strong RI. This either implies that most reproductive isolation (RI) was initiated in allopatry or that RI evolves too rapidly in sympatry to be captured at incipient stages. To distinguish between these explanations, I develop a new “rate test of speciation” that estimates the likelihood of non‐allopatric speciation given the distribution of RI rates in allopatry versus sympatry. Most sympatric taxa were found to have likely initiated RI in allopatry. However, two putative candidate species pairs for non‐allopatric speciation were identified (5% of known Drosophila). In total, this study shows how using RI measures can greatly inform us about the geographical mode of speciation in nature.     

Genetic differentiation among populations of a Hispaniolan trunk anole that exhibit geographical variation in dewlap colour

Long neglected by classic island biogeographical theory, speciation within and among islands is increasingly recognized as a major contributor to insular diversity. Although the factors responsible for island speciation remain poorly understood, this process appears critically dependent on geographical variation and speciation in allopatry or parapatry. Here, we investigate geographical variation and speciation in a complex of Hispaniolan trunk anoles (Anolis distichus), where populations with strikingly distinct dewlap colours and patterns correspond with deeply divergent mtDNA structure. Using a multilocus, population‐level analysis, we investigate whether these phenotypically and mitochondrially distinct populations exhibit the type of nuclear differentiation expected among species or incipient species. Along a transect that extends across a recently recessed marine barrier, our results are consistent with the persistence of an abrupt phenotypic and mitochondrial transition following secondary contact, in spite of little or no evidence for a reduction in nuclear gene flow. Along a second transect extending across a steep environmental gradient, our phenotypic and microsatellite data suggest a sharp genetic break with little or no admixture, whereas mtDNA recovers a signature of extensive unidirectional introgression. Together, these results are consistent with previous studies of Lesser Antillean anoles, suggesting that allopatric divergence alone is insufficient for speciation, whereas reduced gene flow and partial reproductive isolation may accumulate in the presence of ecological gradients.     

Phylogeographic patterns of host-race evolution in Tephritis conura (Diptera: Tephritidae)

Host-race evolution is a prime candidate for sympatric speciation because host shifts must take place in the presence of both hosts. However, the geographic context in which the shift takes place may have strong allopatric or peripatric components if the primary host within a localized area is scarce or even goes extinct. Inference of the relative importance of the geographic mode of speciation may be gained from phylogeographic imprints. Here, we investigate the phylogeography of host races of the tephritid fly Tephritis conura from sympatric, parapatric and allopatric populations of Cirsium heterophyllum and Cirsium oleraceum (Asteraceae) in Europe, for addressing the age and direction, and the geographic context of host-race formation. Haplotype networks of the host races differed significantly in molecular diversity and topology. However, host-race haplotypes were nested within the same network, with a central haplotype H1 being the most common haplotype in both host races. C. heterophyllum flies were genetically more diverse and substructured than flies from C. oleraceum, suggesting that the latter is the derived race. The phylogeographic imprint indicates either that C. heterophyllum flies colonized C. oleraceum during range expansion or that heterophyllum flies went extinct in an area where oleraceum flies persisted (followed by re-immigration). These imprints are in concordance with peripatric diversification, probably in the European Alps and related to the last ice age, where the host-race diversification was largely completed before the postglacial range expansion on C. oleraceum to current areas of sympatry and parapatry with C. heterophyllum.     

Quantitative traits and mode of speciation in Martinique anoles

We investigate extensive quantitative trait variation (dewlap hue, colour pattern, dorsum hue, body proportions and scalation) in the Martinique anole across eight transects representing nascent parapatric ecological speciation, nascent allopatric speciation and allopatric divergence without sufficient genetic structure to suggest speciation. Quantitative trait divergence can be extremely large between adjacent sets of populations, but with one exception that this is associated with difference in habitat rather than past allopatry. Nascent ecological speciation shows the greatest level of quantitative trait divergence across all character sets including those implicated in natural, as well as sexual selection. The sole example of nascent allopatric speciation is associated with fairly strong quantitative trait divergence among most character sets, but not the set most implicated in natural (rather than sexual) selection. The role of sexual selection in ecological speciation is discussed, both in terms of female choice with assortative mating and male–male competition with condition‐dependant sexual signals.     

Proximate sources of marine biodiversity

When temperature and other kinds of barrier divide formerly continuous populations and confine them to more restricted geographical areas, there is an evolutionary reaction that will, over time, result in the formation of endemic species. In such cases, an allopatric speciation process is considered to have taken place because reproductive isolation was caused by physical means instead of by natural selection. In contrast, when populations exist in a very high-diversity area and remain undivided by physical events, they exhibit a tendency to speciate by means of sympatry (or parapatry). This process, sometimes called competitive or ecological speciation, does involve reproductive isolation by means of natural selection. Populations that exist in geographical provinces bounded by physical barriers add to the overall diversity through the production of endemic species. This increase by species packing is relatively slow due to the very gradual tempo of the allopatric speciation process. Populations existing in centres of origin add to the general diversity through the production of species that are dominant in terms of their ability to spread over large parts of the world. It is proposed that such species are usually formed by sympatric speciation, a process that can be c. 20 times faster than species formation by allopatry. It is not suggested that sympatry is exclusive to centres of origin, nor that allopatry is confined to peripheral provinces. Both processes are widespread, but there do appear to be distinctive geographical concentrations. Considering that numbers of widespread species produced by centres of origin may eventually become subdivided by barriers, and thus give rise to descendants by allopatry, it is difficult to say how much of our present species diversity has come from one source or the other. Both speciation by sympatry from centres of origin and speciation by allopatry in peripheral provinces appear to be important sources of marine biodiversity.     

Grapes, galls, and geography: the distribution of nuclear and mitochondrial DNA variation across host-plant species and regions in a specialist herbivore

Studies of patterns of molecular variation in natural populations can provide important insights into a number of evolutionary problems. Among these, the question of whether geographic factors are more important than ecological factors in promoting population differentiation and ultimately speciation has been an important and contentious area in evolutionary biology. Systems involving herbivorous insects have played a leading role in this discussion. This study examined the distribution of molecular variation in a highly specialized gall-forming insect, grape phylloxera (Daktulosphaira vitifoliae Fitch), that is found on both sympatric and allopatric host-plant species of the genus Vitis. In addition, the relationship of insects in the introduced range in the United States to ancestral populations in the native range was examined. Evidence for differentiation along host-plant lines from both nuclear (RAPD) and mitochondrial (COI) DNA was confounded with the effect of geography. Differentiation was found where hosts were allopatric or parapatric, but no evidence was found for such differentiation on two hosts, V. vulpina and V. aestivalis, that are broadly sympatric. The question of population differentiation onto these sympatric hosts can be considered to be resolved--it has not occurred in spite of a long history of association. Evidence was equivocal, but suggestive of a period of divergence in allopatry prior to reestablishment of contact, for insects associated with another host plant species, V. cinerea, found in both sympatric and parapatric populations. A low level of diversity and placement of samples collected from the grape species V. riparia at the tip of a phylogenetic tree supports the hypothesis that this host has been recently colonized from populations from the Mississippi Valley. A polyphyletic origin for biotype B grape phylloxera was supported: Although most samples collected from vineyards in the introduced range in California had similar haplotypes, they were closely related to natives on V. vulpina from the Atlantic Coast-Piedmont region. All samples collected from vineyards in Oregon and Washington were closely related to natives on V. riparia in the northern United States.     

HISTORICAL DIVERGENCE AND GENE FLOW: COALESCENT ANALYSES OF MITOCHONDRIAL,AUTOSOMAL AND SEX‐LINKED LOCI IN PASSERINA BUNTINGS

Quantifying the role of gene flow during the divergence of closely related species is crucial to understanding the process of speciation. We collected DNA sequence data from 20 loci (one mitochondrial, 13 autosomal, and six sex‐linked) for population samples of Lazuli Buntings (Passerina amoena) and Indigo Buntings (Passerina cyanea) (Aves: Cardinalidae) to test explicitly between a strict allopatric speciation model and a model in which divergence occurred despite postdivergence gene flow. Likelihood ratio tests of coalescent‐based population genetic parameter estimates indicated a strong signal of postdivergence gene flow and a strict allopatric speciation model was rejected. Analyses of partitioned datasets (mitochondrial, autosomal, and sex‐linked) suggest the overall gene flow patterns are driven primarily by autosomal gene flow, as there is no evidence of mitochondrial gene flow and we were unable to reject an allopatric speciation model for the sex‐linked data. This pattern is consistent with either a parapatric divergence model or repeated periods of allopatry with gene flow occurring via secondary contact. These results are consistent with the low fitness of female avian hybrids under Haldane's rule and demonstrate that sex‐linked loci likely are important in the initial generation of reproductive isolation, not just its maintenance.     

Phylogeography of speciation: allopatric divergence and secondary contact between outcrossing and selfing Clarkia

The origins of hybrid zones between parapatric taxa have been of particular interest for understanding the evolution of reproductive isolation and the geographic context of species divergence. One challenge has been to distinguish between allopatric divergence (followed by secondary contact) versus primary intergradation (parapatric speciation) as alternative divergence histories. Here, we use complementary phylogeographic and population genetic analyses to investigate the recent divergence of two subspecies of Clarkia xantiana and the formation of a hybrid zone within the narrow region of sympatry. We tested alternative phylogeographic models of divergence using approximate Bayesian computation (ABC) and found strong support for a secondary contact model and little support for a model allowing for gene flow throughout the divergence process (i.e. primary intergradation). Two independent methods for inferring the ancestral geography of each subspecies, one based on probabilistic character state reconstructions and the other on palaeo-distribution modelling, also support a model of divergence in allopatry and range expansion leading to secondary contact. The membership of individuals to genetic clusters suggests geographic substructure within each taxon where allopatric and sympatric samples are primarily found in separate clusters. We also observed coincidence and concordance of genetic clines across three types of molecular markers, which suggests that there is a strong barrier to gene flow. Taken together, our results provide evidence for allopatric divergence followed by range expansion leading to secondary contact. The location of refugial populations and the directionality of range expansion are consistent with expectations based on climate change since the last glacial maximum. Our approach also illustrates the utility of combining phylogeographic hypothesis testing with species distribution modelling and fine-scale population genetic analyses for inferring the geography of the divergence process.     

Cryptic barriers to dispersal within a lake allow genetic differentiation of Eurasian perch

Gene flow between coexisting or nearby populations normally prevents genetic divergence and local adaptation. Despite this, there are an increasing number of reports of sympatric sister taxa, indicating potential divergence and speciation in the face of gene flow. A large number of such reported cases involve lake-dwelling fish, which are expected to run into few physical barriers to dispersal within their aquatic habitat. However, such cases may not necessarily reflect sympatric speciation if cryptic dispersal barriers are common in lakes and other aquatic systems. In this study, we examined genetic differentiation in perch (Perca fluviatilis L.) from nine locations in a single, small lake (24 km(2)), using microsatellites. We detected significant genetic differentiation in all but two pairwise comparisons. These patterns were not consistent with divergence by distance or the existence of kin groups. Instead, they suggest that cryptic barriers to dispersal exist within the lake, allowing small-scale genetic divergence. Such an observation suggests that allopatric (or parapatric) divergence may be possible, even in small, apparently homogenous environments such as lakes. This has important consequences for how we currently view evidence from nature for sympatric speciation.     

Sympatric,parapatric or allopatric: the most important way to classify speciation?

The most common classification of modes of speciation begins with the spatial context in which divergence occurs: sympatric, parapatric or allopatric. This classification is unsatisfactory because it divides a continuum into discrete categories, concentrating attention on the extremes, and it subordinates other dimensions on which speciation processes vary, such as the forces driving differentiation and the genetic basis of reproductive isolation. It also ignores the fact that speciation is a prolonged process that commonly has phases in different spatial contexts. We use the example of local adaptation and partial reproductive isolation in the intertidal gastropod Littorina saxatilis to illustrate the inadequacy of the spatial classification of speciation modes. Parallel divergence in shell form in response to similar environmental gradients in England, Spain and Sweden makes this an excellent model system. However, attempts to demonstrate 'incipient' and 'sympatric' speciation involve speculation about the future and the past. We suggest that it is more productive to study the current balance between local adaptation and gene flow, the interaction between components of reproductive isolation and the genetic basis of differentiation.     

GENETIC EVIDENCE FOR REPRODUCTIVE ISOLATION AND MULTIPLE ORIGINS OF SYMPATRIC TROPHIC ECOTYPES OF WHITEFISH (COREGONUS)

We assessed variation in mitochondrial DNA (mtDNA) by restriction fragment length polymorphism (RFLP) analysis and in nuclear genes by allozyme analysis among sympatric pairs of limnetic and benthic ecotypes of whitefish (Coregonus) coexisting in three lakes of southern Yukon to address three evolutionary questions regarding their origins. Are sympatric low and high gill-raker count ecotypes genetically differentiated? Are they issued from monophyletic or polyphyletic evolutionary events? If they are polyphyletic in origins, did they originate from multiple allopatric speciation events or intralacustrine radiation? Our results corroborated previous genetic and ecological studies of these ecotypes, indicating that they represent genetically distinct reproductive units, and therefore refuting the hypothesis of phenotypic polymorphism within a single population. However, the amount of gene flow between ecotypes varied among lakes, correlating with the extent of morphological differentiation and the potential for premating reproductive isolation. The results indicated a polyphyletic origin of ecotypes whereby each of them have been expressed independently more than once. In the two lakes of Squanga Creek drainage, the existence of sympatric pairs was best explained by the secondary contact of two monophyletic whitefish groups that evolved in allopatry during the last glaciation events. In Dezadeash L. of Alsek R. drainage, our results could not verify either sympatric or allopatric (or microallopatric) origin of ecotypes. Regardless of the mode of speciation involved in their origins, these sympatric whitefish populations provided further evidence that Pleistocene glaciation events created conditions favoring rapid divergence and phenotypic differentiation among northern freshwater fishes.