Random interbreeding between cryptic lineages of the Common Raven: evidence for speciation in reverse

DNA sequence studies frequently reveal evidence of cryptic lineages in morphologically uniform species, many of which turn out to be evolutionarily distinct species. The Common Raven (Corvus corax) includes two deeply divergent mtDNA lineages: one lineage seems restricted to western North America and the other is Holarctic in distribution. These deep clades hint of the possibility of cryptic species in the western United States. We tested this hypothesis in a population consisting of an equal proportion of both mtDNA clades, by quantifying mating patterns and associated fitness consequences with respect to mtDNA. We also tested for morphological, behavioural and ecological correlates of sex and mtDNA clade membership. Mate pairings were random with respect to mtDNA clades, and there were no differences in reproductive success between assortatively and nonassortatively mated pairs. We found no differences in survival or resource use between clades. There were no differences in morphological or behavioural characters between mtDNA clades, except one clade trended towards greater mobility. These results suggest there are no barriers to gene flow between mtDNA clades and argue that the mtDNA clades have remerged in this population, likely due to a lack of ecological or signal differentiation between individuals in each lineage. Hence, in Common Ravens, phylogeographic structure in mtDNA is a reflection of likely past isolation rather than currently differentiated species.     

Georgian and kurd mtDNA sequence analysis shows a lack of correlation between languages and female genetic lineages

Mitochondrial DNA sequences from Georgians and Kurds were analyzed in order to test the possible correlation between female lineages and languages in these two neighboring West Eurasian groups. Mitochondrial sequence pools in both populations are very similar despite their different linguistic and prehistoric backgrounds. Both populations present mtDNA lineages that clearly belong to the European gene pool, as shown by 1) similar nucleotide and sequence diversities; 2) a large number of sequences shared with the rest of European samples; 3) nonsignificant genetic distances; and 4) classification of the present lineages into the major European mtDNA haplogroups already described. The outlier position of the populations from the Caucasus according to classical genetic markers is not recognized in the present Georgian mtDNA sequence pool. This result suggests that the differentiation of mtDNA sequences in West Eurasia and the outlier features of Caucasian populations should be attributed to different processes. Moreover, the putative linguistic relationship between Caucasian groups and the Basques, another outlier population within Europe for classical genetic markers, is not detected by the analysis of mtDNA sequences.     

Nuclear DNA recapitulates the cryptic mitochondrial lineages of Lumbricus rubellus and suggests the existence of cryptic species in an ecotoxological soil sentinel

Mitochondrial DNA analysis has revealed two distinct phylogenetic lineages within the ecotoxological sentinel earthworm model Lumbricus rubellus Hoffmeister, 1843. The existence of these lineages could complicate ecotoxicological studies that use the species as a sentinel for soil contamination testing, as they may respond differently to contamination; however, as mitochondrial haplotypes are not always expected to segregate in the same way as chromosomal DNA in natural populations, we further investigated this issue by using nuclear DNA markers (microsatellites) to measure genetic diversity, differentiation, and gene flow in sympatric populations of the two L. rubellus lineages at two sites in South Wales. Our results show that sympatric populations of the two lineages are more genetically differentiated than geographically distant populations of the same lineage, and Bayesian clustering analysis revealed no evidence of gene flow between the lineages at either site. Additionally, DNA sequencing of these microsatellite loci uncovered substantial differentiation between lineages at homologous flanking regions. Overall our findings indicate a high degree of nuclear genetic differentiation between the two lineages of L. rubellus, implying reproductive isolation at the two study sites and therefore the potential existence of cryptic species. The existence of two cryptic taxa has major implications for the application of L. rubellus as an ecotoxicological sentinel. It may therefore be necessary to consider the lineages as separate taxa during future ecotoxicological studies. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 780–795.     

From glacial refugia to modern populations: new assemblages of organelle genomes generated by differential cytoplasmic gene flow in transcontinental black spruce

Assessing species' range-wide cytoplasmic diversity provides valuable insights regarding their dispersal and adaptive potential in a changing environment. Transcontinental chloroplast (cpDNA) and mitochondrial DNA (mtDNA) population structures were compared to identify putative ancestral and new cytoplasmic genome assemblages in black spruce (Picea mariana), a North American boreal conifer. Mean within-population diversity and allelic richness for cpSSR markers were 0.80 and 4.21, respectively, and diminished westward. Population differentiation based on G(ST) was lower for cpDNA than for mtDNA (G(ST) =0.104 and 0.645, respectively) but appeared comparable when estimated using Jost differentiation index (D=0.459 and 0.537, respectively). Further analyses resulted in the delineation of at least three genetically distinct cpDNA lineages partially congruent with those inferred from mtDNA data, which roughly corresponded to western, central and eastern Canada. Additionally, the patterns of variation in Alaska for both cpDNA and mtDNA markers suggested that black spruce survived the last glacial maximum in this northern region. The range-wide comparison of the geographic extent of cytoplasmic DNA lineages revealed that extensive pollen gene flow between ancestral lineages occurred preferentially from west to east during the postglacial expansion of the species, while seed-mediated gene flow remained geographically restricted. This differential gene flow promoted intraspecific cytoplasmic capture that generated new assemblages of cpDNA and mtDNA genomes during the Holocene. Hence, black spruce postglacial colonization unexpectedly resulted in an increase in genetic diversity with possible adaptive consequences.     

Mitochondrial DNA diversity in South America and the genetic history of Andean highlanders

We analyzed mtDNA sequence variation in 590 individuals from 18 south Amerindian populations. The spatial pattern of mtDNA diversity in these populations fits well the model proposed on the basis of Y-chromosome data. We found evidence of a differential action of genetic drift and gene flow in western and eastern populations, which has led to genetic divergence in the latter but not in the former. Although it is not possible to identify a pattern of genetic variation common to all South America, when western and eastern populations are analyzed separately, the mtDNA diversity in both regions fits the isolation-by-distance model, suggesting independent evolutionary dynamics. Maximum-likelihood estimates of divergence times between central and south Amerindian populations fall between 13,000 and 19,000 years, which is consistent with a Pleistocenic peopling of South America. Moreover, comparison of among-population variability of mtDNA and Y-chromosome DNA seems to indicate that South America is the only continent where the levels of differentiation are similar for maternal and paternal lineages.     

Population analysis using the nuclear white gene detects Pliocene/Pleistocene lineage divergence within Anopheles nuneztovari in South America

Anopheles (Nyssorhynchus) nuneztovari Gabaldón (Diptera: Culicidae), a locally important malaria vector in some regions of South America, has been hypothesized to consist of at least two cryptic incipient species. We investigated its phylogeographic structure in several South American localities to determine the number of lineages and levels of divergence using the nuclear white gene, a marker that detected two recently diverged genotypes in the primary Neotropical malaria vector Anopheles darlingi Root. In An. nuneztovari, five distinct lineages (1-5) were elucidated: (1) populations from northeastern and central Amazonia; (2) populations from Venezuela east and west of the Andes; (3) populations from Colombia and Venezuela west of the Andes; (4) southeastern and western Amazonian Brazil populations, and (5) southeastern and western Amazonian Brazil and Bolivian populations. There was a large amount of genetic differentiation among these lineages. The deepest and earliest divergence was found between lineage 3 and lineages 1, 2 and 4, which probably accounts for the detection of lineage 3 in some earlier studies. The multiple lineages within Amazonia are partially congruent with previous mtDNA and ITS2 data, but were undetected in many earlier studies, probably because of their recent (Pleistocene) divergence and the differential mutation rates of the markers. The estimates for the five lineages, interpreted as recently evolved or incipient species, date to the Pleistocene and Pliocene. We hypothesize that the diversification in An. nuneztovari is the result of an interaction between the Miocene/Pliocene marine incursion and Pleistocene climatic changes leading to refugial isolation. The identification of cryptic lineages in An. nuneztovari could have a significant impact on local vector control measures.     

Genetic population structure of the paper wasp Polistes olivaceus (Hymenoptera: Vespidae) in Bangladesh

Dispersal triggers gene flow, which in turn strongly affects the ensuing genetic population structure of a species. Using nuclear microsatellite loci and mitochondrial DNA (mtDNA), we estimated the genetic population structure of the wasp Polistes olivaceus throughout Bangladesh. The level of population differentiation using nuclear markers (F _ST) appeared to be much lower than that estimated using mtDNA haplotype sequences (Ф_ST), even after correcting for effective population size differences between the two markers. These results suggest a philopatric tendency, in which gynes disperse less than males. We observed no isolation by distance among the study populations at either the nuclear or mtDNA level, suggesting nonequilibrium between gene flow and drift as a result of very frequent interpopulation movement. For the nuclear markers, an individual assignment test showed no genetically and geographically distinct groups. Instead, phylogenetic analyses as well as a minimum spanning network using mtDNA haplotypes consistently revealed two distinct lineages. The distribution of haplotypes indicated western populations with a single lineage and offered clear evidence for restricted gene flow across the Jamuna–Padma–Upper Meghna river system. Mismatch distributions exhibited a unimodal distribution, which along with a starlike haplotype network, suggested a population expansion in lineage I but not in lineage II. Overall, these results suggest that gene flow among populations of P. olivaceus was affected by both female philopatry and a major river system across Bangladesh.     

Inferring the degree of incipient speciation in secondary contact zones of closely related lineages of Palearctic green toads (Bufo viridis subgroup)

Reproductive isolation between lineages is expected to accumulate with divergence time, but the time taken to speciate may strongly vary between different groups of organisms. In anuran amphibians, laboratory crosses can still produce viable hybrid offspring >20 My after separation, but the speed of speciation in closely related anuran lineages under natural conditions is poorly studied. Palearctic green toads (Bufo viridis subgroup) offer an excellent system to address this question, comprising several lineages that arose at different times and form secondary contact zones. Using mitochondrial and nuclear markers, we previously demonstrated that in Sicily, B. siculus and B. balearicus developed advanced reproductive isolation after Plio-Pleistocene divergence (2.6 My, 3.3–1.9), with limited historic mtDNA introgression, scarce nuclear admixture, but low, if any, current gene flow. Here, we study genetic interactions between younger lineages of early Pleistocene divergence (1.9 My, 2.5–1.3) in northeastern Italy (B. balearicus, B. viridis). We find significantly more, asymmetric nuclear and wider, differential mtDNA introgression. The population structure seems to be molded by geographic distance and barriers (rivers), much more than by intrinsic genomic incompatibilities. These differences of hybridization between zones may be partly explained by differences in the duration of previous isolation. Scattered research on other anurans suggests that wide hybrid zones with strong introgression may develop when secondary contacts occur <2 My after divergence, whereas narrower zones with restricted gene flow form when divergence exceeds 3 My. Our study strengthens support for this rule of thumb by comparing lineages with different divergence times within the same radiation.     

Phylogenetic Relationships Among Fragmented Asian Black Bear (Ursus Thibetanus) Populations in Western Japan

The number of Asian black bears (Ursus thibetanus) in Japan has been reduced and their habitats fragmented and isolated because of human activities. Our previous study examining microsatellite DNA loci revealed significant genetic differentiation among four local populations in the western part of Honshu. Here, an approximate 700-bp nucleotide sequence of mitochondrial DNA (mtDNA) control region was analysed in 119 bears to infer the evolutionary history of these populations. Thirteen variable sites and variation in the number of Ts at a T-repeat site were observed among the analysed sequences, which defined 20 mtDNA haplotypes with the average sequence divergence of 0.0051 (SD = 0.00001). The observed haplotype frequencies differed significantly among the four populations. Phylogeographic analysis of the haplotypes suggested that black bears in this region have gone through two different colonisation histories, since the observed haplotypes belonged to two major monophyletic lineages and the lineages were distributed with an apparent border. The spatial genetic structure revealed by using mtDNA was different from that observed using microsatellite DNA markers, probably due to female philopatry and male-biased dispersal. Since nuclear genetic diversity will be lost in the three western populations because of the small population size and genetic isolation, their habitats need to be preserved, and these four populations should be linked to each other by corridors to promote gene flow from the easternmost population with higher nuclear genetic diversity.     

Incongruence of mitochondrial and nuclear gene trees in the Carabid beetles Ohomopterus

We studied the molecular phylogeny of the carabid subgenus Ohomopterus (genus Carabus), using two mitochondrial (mt) DNA regions (16SrRNA and NADH dehydrogenase subunit 5) and three nuclear DNA regions (wingless, phosphoenolpyruvate carboxykinase, and an anonymous locus). We revisited the previously reported incongruence between the distribution of mtDNA markers and morphologically defined species (Su et al., 1996; J. Mol. Evol. 43:662-671), which those authors attributed to "type switching", a concerted change in many morphological characters that results in the repeated evolution of a particular morphological type. Our mtDNA gene tree obtained from 44 individuals representing all 15 currently recognized species of Ohomopterus revealed that haplotypes isolated from individuals of a single "species" were frequently separated into distant clades, confirming the previous report. The three nuclear markers generally conformed better-with the morphologically defined species than did the mitochondrial markers. The phylogenetic signal in mtDNA and nuclear DNA data differed strongly, and these two partitions were significantly incongruent with each other according to the incongruence length difference test of Farris et al. (1994; Cladistics 10:315-320), although the three nuclear partitions were not homogeneous either. Our results did not support the type-switching hypothesis that had been proposed to fit the morphological data to the mitochondrial gene tree: The incongruence of the mtDNA tree with other nuclear markers indicates that the mtDNA-based tree does not reflect species history any better than the morphological data do. Incongruence of gene trees in Ohomopterus may have been promoted by the complex processes of geographic isolation and hybridization in the Japanese Archipelago that have led to occasional gene flow and recombination between separated entities. The occurrence of reticulate patterns in this group is intriguing, because species of Ohomopterus exhibit extremely divergent genitalic structures that represent a highly efficient reproductive isolation mechanism.     

Phylogeographic sympatry and isolation of the Eurasian badgers (Meles, Mustelidae, Carnivora): Implications for an alternative analysis using maternally as well as paternally inherited genes

In the present study, to further understand the phylogenetic relationships among the Eurasian badgers (Meles, Mustelidae, Carnivora), which are distributed widely in the Palearctic, partial sequences of the mitochondrial DNA (mtDNA) control region (539-545 base-pairs) as a maternal genetic marker, and the sex-determining region on the Y-chromosome gene (SRY: 1052-1058 base-pairs), as a paternal genetic marker, were examined. The present study revealed ten SRY haplotypes from 47 males of 112 individuals of the Eurasian Continent and Japan. In addition, 39 mtDNA haplotypes were identified from those animals. From the phylogeography of both the uniparentally inherited genes, four lineages were recognized as Japanese, eastern Eurasian, Caucasian, and western Eurasian. The distribution patterns of the mtDNA lineages showed the existence of a sympatric zone between the eastern and western Eurasian lineages around the Volga River in western Russia. Furthermore, the present study suggested that in the Japanese badgers, the larger genetic differentiation of the Shikoku population was attributable to geographic history in the Japanese islands.     

Reproductive isolation among cryptic species in the ectomycorrhizal genus Strobilomyces: population-level CAPS marker-based genetic analysis

Among the higher fungi, reproductively isolated cryptic species exist that are morphologically difficult to distinguish owing to a lack of taxonomically useful morphological characters. Mating tests are helpful for detecting reproductive isolation between cryptic species, but are often difficult to perform for higher fungi, especially ectomycorrhizal fungi. In order to identify cryptic species of the ectomycorrhizal genus Strobilomyces more efficiently, lineages were defined based on the nucleotide sequence of two mitochondrial genes. Then the gene flow among lineages was measured using cleaved amplified polymorphic sequences (CAPS) markers designed for single copy nuclear genes. No heterozygosity was observed between different lineages, but within the same lineage heterozygosity was present at the ratio expected given Hardy Weinberg equilibrium. These results show that the mtDNA lineages are separate Mendelian populations, possibly cryptic species that are reproductively isolated from each other.     

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.     

Less Pollen-Mediated Gene Flow for More Signatures of Glacial Lineages: Congruent Evidence from Balsam Fir cpDNA and mtDNA for Multiple Refugia in Eastern and Central North America

The phylogeographic structure and postglacial history of balsam fir (Abies balsamea), a transcontinental North American boreal conifer, was inferred using mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) markers. Genetic structure among 107 populations (mtDNA data) and 75 populations (cpDNA data) was analyzed using Bayesian and genetic distance approaches. Population differentiation was high for mtDNA (dispersed by seeds only), but also for cpDNA (dispersed by seeds and pollen), indicating that pollen gene flow is more restricted in balsam fir than in other boreal conifers. Low cpDNA gene flow in balsam fir may relate to low pollen production due to the inherent biology of the species and populations being decimated by recurrent spruce budworm epidemics, and/or to low dispersal of pollen grains due to their peculiar structural properties. Accordingly, a phylogeographic structure was detected using both mtDNA and cpDNA markers and population structure analyses supported the existence of at least five genetically distinct glacial lineages in central and eastern North America. Four of these would originate from glacial refugia located south of the Laurentide ice sheet, while the last one would have persisted in the northern Labrador region. As expected due to reduced pollen-mediated gene flow, congruence between the geographic distribution of mtDNA and cpDNA lineages was higher than in other North American conifers. However, concordance was not complete, reflecting that restricted but nonetheless detectable cpDNA gene flow among glacial lineages occurred during the Holocene. As a result, new cpDNA and mtDNA genome combinations indicative of cytoplasmic genome capture were observed.     

Multiple glacial refugia and contemporary dispersal shape the genetic structure of an endemic amphibian from the Pyrenees

Historical factors (colonization scenarios, demographic oscillations) and contemporary processes (population connectivity, current population size) largely contribute to shaping species’ present‐day genetic diversity and structure. In this study, we use a combination of mitochondrial and nuclear DNA markers to understand the role of Quaternary climatic oscillations and present‐day gene flow dynamics in determining the genetic diversity and structure of the newt Calotriton asper (Al. Dugès, 1852), endemic to the Pyrenees. Mitochondrial DNA did not show a clear phylogeographic pattern and presented low levels of variation. In contrast, microsatellites revealed five major genetic lineages with admixture patterns at their boundaries. Approximate Bayesian computation analyses and linear models indicated that the five lineages likely underwent separate evolutionary histories and can be tracked back to distinct glacial refugia. Lineage differentiation started around the Last Glacial Maximum at three focal areas (western, central and eastern Pyrenees) and extended through the end of the Last Glacial Period in the central Pyrenees, where it led to the formation of two more lineages. Our data revealed no evidence of recent dispersal between lineages, whereas borders likely represent zones of secondary contact following expansion from multiple refugia. Finally, we did not find genetic evidence of sex‐biased dispersal. This work highlights the importance of integrating past evolutionary processes and present‐day gene flow and dispersal dynamics, together with multilocus approaches, to gain insights into what shaped the current genetic attributes of amphibians living in montane habitats.     

Genetic differentiation and hybridization in two naturally occurring sympatric trout Salmo spp. forms from a small karstic lake

In this study, multiple molecular markers [genotyping of 12 nuclear microsatellite loci and the protein‐coding gene ldh‐c1* plus sequencing of the mitochondrial DNA (mtDNA) control region] were employed to investigate the genetic structure of the two trout forms, Salmo cettii and Salmo fibreni, inhabiting Lake Posta Fibreno, central Italy. The two forms were found to share a unique mtDNA haplotype, belonging to a widespread Mediterranean haplogroup (AD). Bayesian clustering analyses showed that these two forms correspond to well‐defined autochthonous gene pools. Genetic introgression between the two gene pools, however, was observed, whose frequency appears to correlate with the environmental features of the spawning sites. The interplay of selection for the spawning sites, philopatry and natural selection can be argued to maintain genetic differentiation despite the lack of complete reproductive isolation.     

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y‐chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least‐divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi‐absence of Y‐chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non‐neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.     

Evolutionary and demographic processes shaping geographic patterns of genetic diversity in a keystone species,the African forest elephant (Loxodonta cyclotis)

The past processes that have shaped geographic patterns of genetic diversity may be difficult to infer from current patterns. However, in species with sex differences in dispersal, differing phylogeographic patterns between mitochondrial (mt) and nuclear (nu) DNA may provide contrasting insights into past events. Forest elephants (Loxodonta cyclotis) were impacted by climate and habitat change during the Pleistocene, which likely shaped phylogeographic patterns in mitochondrial (mt) DNA that have persisted due to limited female dispersal. By contrast, the nuclear (nu) DNA phylogeography of forest elephants in Central Africa has not been determined. We therefore examined the population structure of Central African forest elephants by genotyping 94 individuals from six localities at 21 microsatellite loci. Between forest elephants in western and eastern Congolian forests, there was only modest genetic differentiation, a pattern highly discordant with that of mtDNA. Nuclear genetic patterns are consistent with isolation by distance. Alternatively, male‐mediated gene flow may have reduced the previous regional differentiation in Central Africa suggested by mtDNA patterns, which likely reflect forest fragmentation during the Pleistocene. In species like elephants, male‐mediated gene flow erases the nuclear genetic signatures of past climate and habitat changes, but these continue to persist as patterns in mtDNA because females do not disperse. Conservation implications of these results are discussed.     

Phylogeography and population history of Atlantic mackerel (Scomber scombrus L.): a genealogical approach reveals genetic structuring among the eastern Atlantic stocks

Despite the resolving power of DNA markers, pelagic and migratory marine fish species generally show very little geographical population structuring. In mackerel (Scomber scombrus L.) population differentiation has been detected only at a transatlantic scale. By applying two regions in mitochondrial DNA (mtDNA) (D-loop and cytochrome b (cytb)) in combination with genealogical and frequency-based statistical approaches, our data suggest population differentiation among eastern Atlantic spawning stocks. In contrast, and indicative of homing behaviour, no genetic structuring was observed among shoals of individuals outside the spawning season. Among spawning stocks, mtDNA D-loop sequences detected differentiation within the eastern Atlantic, while the cytb gene detected transatlantic differentiation. The impact of recurrent events (e.g. gene flow restricted by isolation by distance) and historic events (e.g. population range expansions) among spawning stocks was investigated applying a nested cladistic analysis of geographical distribution of cytb haplotype lineages. In the eastern Atlantic, historical population range expansion, presumably in connection with recolonization of northern areas after the last glaciation, is suggested to be the main factor determining mtDNA lineage distribution. This was supported by estimates of mtDNA nucleotide diversity, where the highest diversity was observed for the stock spawning in the Bay of Biscay, for which the size estimate is only 15% of the largest stock (Celtic Sea). In addition to revealing population differentiation, our data demonstrate the importance of sampling strategy and the power of applying statistical methods addressing both ongoing and historical population processes.     

Historical introgression and the role of selective vs. neutral processes in structuring nuclear genetic variation (AFLP) in a circumpolar marine fish, the capelin (Mallotus villosus)

The capelin (Mallotus villosus) is a widespread marine fish species for which previous work has identified geographically distinct mtDNA clades, the frontiers of which are well within adult and larval dispersal capabilities. Here, we use AFLPs to test for the presence of nuclear gene flow among clades. In addition, we evaluate genetic structuring within one clade, the Northwest Atlantic (NWA). We found that each of the mtDNA clades corresponds with a unique nuclear DNA genetic cluster. Within the NWA clade, we detected individuals with small but significant amounts of genetic ancestry from other clades, likely due to historical introgression. Further support for historical introgression comes from analyses of variance in locus-specific differentiation, which support introgression between some clades and divergence without gene flow between others. Within the NWA, we identified two genetic clusters that correspond to sites in geographically adjacent areas. However, these clusters differ primarily at 'outlier' loci, and a genetic subdivision (K=2) was not supported by genetic clustering programs using neutral loci. Significant neutral F(ST) differentiation was found only between sites that otherwise differed at outlier loci. Thus, these populations may be in the initial stages of 'isolation by adaptation'. These results suggest strong between-clade reproductive isolation despite opportunities for gene flow and support the hypothesis that selection can contribute to divergence in otherwise 'open' systems.