Parasitic flatworms infecting thorny skate,Amblyraja radiata: Infection by the monogeneans Acanthocotyle verrilli and Rajonchocotyle emarginata in Svalbard

Parasite diversity above the Arctic circle remains understudied even for commercially valuable host taxa. Thorny skate, Amblyraja radiata, is a common bycatch species with a growing commercial value. Its natural range covers both sides of the North Atlantic including the Arctic zone. Svalbard is a Norwegian archipelago located on the northwest corner of the Barents Shelf which sustains a spectacular species diversity. So far, several monogenean species have been reported infecting thorny skate across the Atlantic Ocean. In the present study, we intend to fill in the knowledge gap on monogenean parasites infecting thorny skate in the northern part of its range and thus indirectly assess the connectivity between the thorny skate populations off the Svalbard coast and from previously studied locations. 46 monogenean individuals were recovered from 11 specimens of thorny skate. Following morphological and molecular assessment, two species of monogeneans, Acanthocotyle verrilli and Rajonchocotyle emarginata, were identified. The results serve as the northernmost record for both parasite genera and the first record of monogenean species off Svalbard. Detailed morphometric evaluation revealed a relatively high level of morphological variation in A. verrilli compared to its congeners. Phylogenetic reconstruction placed A. verrilli in a well-supported clade with A. imo. Our study also suggests high diagnostic significance of sclerotised structures in the identification of Rajonchocotyle. Even though the occurrence of two directly transmitted parasite species supports the previously suggested long-distance migration of A. radiata, future studies employing highly variable genetic markers are needed to assess the ongoing and historical migration patterns.     

RAD sequencing reveals genomewide divergence between independent invasions of the European green crab (Carcinus maenas) in the Northwest Atlantic

Genomic studies of invasive species can reveal both invasive pathways and functional differences underpinning patterns of colonization success. The European green crab (Carcinus maenas) was initially introduced to eastern North America nearly 200 years ago where it expanded northwards to eastern Nova Scotia. A subsequent invasion to Nova Scotia from a northern European source allowed further range expansion, providing a unique opportunity to study the invasion genomics of a species with multiple invasions. Here, we use restriction‐site‐associated DNA sequencing‐derived SNPs to explore fine‐scale genomewide differentiation between these two invasions. We identified 9137 loci from green crab sampled from 11 locations along eastern North America and compared spatial variation to mitochondrial COI sequence variation used previously to characterize these invasions. Overall spatial divergence among invasions was high (pairwise F_ST ~0.001 to 0.15) and spread across many loci, with a mean F_ST ~0.052 and 52% of loci examined characterized by F_ST values >0.05. The majority of the most divergent loci (i.e., outliers, ~1.2%) displayed latitudinal clines in allele frequency highlighting extensive genomic divergence among the invasions. Discriminant analysis of principal components (both neutral and outlier loci) clearly resolved the two invasions spatially and was highly correlated with mitochondrial divergence. Our results reveal extensive cryptic intraspecific genomic diversity associated with differing patterns of colonization success and demonstrates clear utility for genomic approaches to delineating the distribution and colonization success of aquatic invasive species.     

Genetic discontinuity among regional populations of <Emphasis Type="Italic">Lophelia pertusa</Emphasis> in the North Atlantic Ocean

Knowledge of the degree to which populations are connected through larval dispersal is imperative to effective management, yet little is known about larval dispersal ability or population connectivity in Lophelia pertusa, the dominant framework-forming coral on the continental slope in the North Atlantic Ocean. Using nine microsatellite DNA markers, we assessed the spatial scale and pattern of genetic connectivity across a large portion of the range of L. pertusa in the North Atlantic Ocean. A Bayesian modeling approach found four distinct genetic groupings corresponding to ocean regions: Gulf of Mexico, coastal southeastern U.S., New England Seamounts, and eastern North Atlantic Ocean. An isolation-by-distance pattern was supported across the study area. Estimates of pairwise population differentiation were greatest with the deepest populations, the New England Seamounts (average F _ST = 0.156). Differentiation was intermediate with the eastern North Atlantic populations (F _ST = 0.085), and smallest between southeastern U.S. and Gulf of Mexico populations (F _ST = 0.019), with evidence of admixture off the southeastern Florida peninsula. Connectivity across larger geographic distances within regions suggests that some larvae are broadly dispersed. Heterozygote deficiencies were detected within the majority of localities suggesting deviation from random mating. Gene flow between ocean regions appears restricted, thus, the most effective management scheme for L. pertusa involves regional reserve networks.     

Phylogeographic analysis reveals two cryptic species of the endangered fern <Emphasis Type="Italic">Ceratopteris</Emphasis><Emphasis Type="Italic">thalictroides</Emphasis> (L.) Brongn. (Parkeriaceae) in China

Ceratopteris thalictroides (L.) Brongn. (Parkeriaceae) is a difficult fern species to taxonomically classify. Three cryptic species were revealed in the previous studies, referred to as the north type, the south type, and the third type. Because much of the distribution range of C. thalictroides in China was not included in the sampling of the previous studies, the taxonomic complexity of C. thalictroides in China remained uncertain. In order to identify the uncharacterized cryptic species, we examined four chloroplast DNA (cpDNA) non-coding regions and compared sequence variation within this species complex. Sequence data were obtained from 143 individuals in 24 populations throughout the natural distribution of the species in China. Nineteen haplotypes were identified. Molecular systematic and phylogeographical analyses revealed two genetically distinct clusters of cpDNA haplotypes in China. One cluster included haplotypes associated with the north type, and another with the south type cryptic species. The N _ST value was significantly higher than the G _ST value (N _ST = 0.768 > G _ST = 0.434, P < 0.05), indicating the presence of a significant phylogeographical structure of C. thalictroides in China. The results of AMOVA analysis showed a significant inter-group differentiation (F _ST = 0.918; P < 0.001). Analyses based on different, but complementary methods suggest that in China, C. thalictroides contains only two of the cryptic species (the north and south types). Two haplotypes, H8 and H17, of the interior node in the minimum-spanning network (MSN) of cpDNA haplotypes are widespread. The origin of the widespread haplotypes in China may have resulted from long-distance dispersal to China.     

Population genetic history of the dreissenid mussel invasions: expansion patterns across North America

This study tests population genetic patterns across the Eurasian dreissenid mussel invasions of North America—encompassing the zebra mussel Dreissena polymorpha (1986 detection) and the quagga mussel D. rostriformis bugensis (detected in 1990, which now has largely displaced the former in the Great Lakes). We evaluate their source-spread relationships and invasion genetics using 9–11 nuclear microsatellite loci for 583 zebra mussels (21 sites) and 269 quagga mussels (12 sites) from Eurasian and North American range locations, with the latter including the Great Lakes, Mississippi River basin, Atlantic coastal waterways, Colorado River system, and California reservoirs. Additionally, mtDNA cytochrome b gene sequences are used to verify species identity. Our results indicate that North American zebra mussels originate from multiple non-native northern European populations, whereas North American quagga mussels trace to native estuaries in the Southern Bug and Dnieper Rivers. Invasive populations of both species show considerable genetic diversity and structure (zebra F _ST = 0.006–0.263, quagga F _ST = 0.008–0.267), without founder effects. Most newer zebra mussel populations have appreciable genetic diversity, whereas quagga mussel populations from the Colorado River and California show some founder effects. The population genetic composition of both species changed over time at given sites; with some adding alleles from adjacent populations, some losing them, and all retaining closest similarity to their original composition. Zebra mussels from Kansas and California appear genetically similar and assign to a possible origin from the St. Lawrence River, whereas quagga mussels from Nevada and California assign to a possible origin from Lake Ontario. These assignments suggest that overland colonization pathways via recreational boats do not necessarily reflect the most proximate connections. In conclusion, our microsatellite results comprise a valuable baseline for resolving present and future dreissenid mussel invasion pathways.     

Phylogeography and population genetics of the cryptic bonnethead shark Sphyrna aff. tiburo in Brazil and the Caribbean inferred from mtDNA markers

Resolving the identity, phylogeny and distribution of cryptic species within species complexes is an essential precursor to management. The bonnethead shark, Sphyrna tiburo, is a small coastal shark distributed in the Western Atlantic from North Carolina (U.S.A.) to southern Brazil. Genetic analyses based on mitochondrial markers revealed that bonnethead sharks comprise a species complex with at least two lineages in the Northwestern Atlantic and the Caribbean (S. tiburo and Sphyrna aff. tiburo, respectively). The phylogeographic and phylogenetic analysis of two mitochondrial markers [control region (mtCR) and cytochrome oxidase I (COI)] showed that bonnethead sharks from southeastern Brazil correspond to S. aff. tiburo, extending the distribution of this cryptic species >5000 km. Bonnethead shark populations are only managed in the U.S.A. and in the 2000s were considered to be regionally extinct or collapsed in southeast Brazil. The results indicate that there is significant genetic differentiation between S. aff. tiburo from Brazil and other populations from the Caribbean (Φ_ST = 0.9053, P < 0.000), which means that collapsed populations in the former are unlikely to be replenished from Caribbean immigration. The species identity of bonnethead sharks in the Southwest Atlantic and their relationship to North Atlantic and Caribbean populations still remains unresolved. Taxonomic revision and further sampling are required to reevaluate the status of the bonnethead shark complex through its distribution range.     

GENETIC POPULATION STRUCTURE OF PSEUDO‐NITZSCHIA PUNGENS (BACILLARIOPHYCEAE) FROM THE PACIFIC NORTHWEST AND THE NORTH SEA1

Several species of the diatom Pseudo‐nitzschia produce the neurotoxin domoic acid (DA). Consumption of fish and shellfish that have accumulated this potent excitotoxin has resulted in severe illness and even death in humans, marine mammals, and seabirds. Pseudo‐nitzschia pungens (Grunow ex Cleve) Hasle is a cosmopolitan diatom commonly occurring in the waters of the Pacific Northwest (PNW) and the eastern North Atlantic, including the North Sea. However, genetic and physiological relationships among populations throughout this large geographic distribution have not been assessed. Population genetic parameters (e.g., Hardy–Weinberg equilibrium, linkage equilibrium, F_ST) calculated for P. pungens collected from the Juan de Fuca eddy region in the PNW indicated the presence of two distinct groups that were more divergent from each other than either was from a P. pungens sample from the North Sea. Geographic heterogeneity was also detected within each of the two PNW groups. These results suggested that the populations of P. pungens recently mixed in the Juan de Fuca eddy region (a seasonally retentive feature off the coasts of Washington State, USA, and Vancouver Island, Canada) but did not exchange genetic material by sexual reproduction. Alternatively, these two groups may be cryptic (morphologically identical, but reproductively isolated) species. Identifying cryptic diversity in Pseudo‐nitzschia is important for bloom prediction and aiding the identification of molecular markers that can be used for rapid detection assay development.     

Historical influences dominate the population genetic structure of a sedentary marine fish,Atlantic wolffish (Anarhichas lupus), across the North Atlantic Ocean

Genetic variation was assessed in Atlantic wolffish, Anarhichas lupus, across the North Atlantic Ocean using microsatellite and amplified fragment length polymorphism (AFLP) markers. Despite unusual life history attributes such as large benthic eggs, large larvae, a limited pelagic stage and relatively sedentary adults, which suggest potential for strong population structure, range‐wide F_ST values were comparable to other marine fishes (≤0.035). Nevertheless, both significant genetic differentiation among regions and isolation by distance were observed, suggesting limited dispersal in this species. AFLP loci, evaluated on a subset of samples, revealed slightly higher F_ST values, but similar patterns of differentiation and isolation‐by‐distance estimates, compared to microsatellites. The genetic structure of Atlantic wolffish has likely been shaped by its post‐glacial history of recolonization, North Atlantic current patterns and continuity of habitat on continental shelves.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Coast to coast: High genomic connectivity in North American scoters

Dispersal shapes demographic processes and therefore is fundamental to understanding biological, ecological, and evolutionary processes acting within populations. However, assessing population connectivity in scoters (Melanitta sp.) is challenging as these species have large spatial distributions that span remote landscapes, have varying nesting distributions (disjunct vs. continuous), exhibit unknown levels of dispersal, and vary in the timing of the formation of pair bonds (winter vs. fall/spring migration) that may influence the distribution of genetic diversity. Here, we used double‐digest restriction‐associated DNA sequence (ddRAD) and microsatellite genotype data to assess population structure within the three North American species of scoter (black scoter, M. americana; white‐winged scoter, M. deglandi; surf scoter, M. perspicillata), and between their European congeners (common scoter, M. nigra; velvet scoter, M. fusca). We uncovered no or weak genomic structure (ddRAD Φ_ST < 0.019; microsatellite F_ST < 0.004) within North America but high levels of structure among European congeners (ddRAD Φ_ST > 0.155, microsatellite F_ST > 0.086). The pattern of limited genomic structure within North America is shared with other sea duck species and is often attributed to male‐biased dispersal. Further, migratory tendencies (east vs. west) of female surf and white‐winged scoters in central Canada are known to vary across years, providing additional opportunities for intracontinental dispersal and a mechanism for the maintenance of genomic connectivity across North America. In contrast, the black scoter had relatively elevated levels of divergence between Alaska and Atlantic sites and a second genetic cluster found in Alaska at ddRAD loci was concordant with its disjunct breeding distribution suggestive of a dispersal barrier (behavioral or physical). Although scoter populations appear to be connected through a dispersal network, a small percentage (<4%) of ddRAD loci had elevated divergence which may be useful in linking areas (nesting, molting, staging, and wintering) throughout the annual cycle.     

Population structure of the dusky pipefish (Syngnathus floridae) from the Atlantic and Gulf of Mexico,as revealed by mitochondrial DNA and microsatellite analyses

Aim To elucidate the historical phylogeography of the dusky pipefish (Syngnathus floridae) in the North American Atlantic and Gulf of Mexico ocean basins. Location Southern Atlantic Ocean and northern Gulf of Mexico within the continental United States. Methods A 394‐bp fragment of the mitochondrial cytochrome b gene and a 235‐bp fragment of the mitochondrial control region were analysed from individuals from 10 locations. Phylogenetic reconstruction, haplotype network, mismatch distributions and analysis of molecular variance were used to infer population structure between ocean basins and time from population expansion within ocean basins. Six microsatellite loci were also analysed to estimate population structure and gene flow among five populations using genetic distance methods (F_ST, Nei’s genetic distance), isolation by distance (Mantel’s test), coalescent‐based estimates of genetic diversity and migration patterns, Bayesian cluster analysis and bottleneck simulations. Results Mitochondrial analyses revealed significant structuring between ocean basins in both cytochrome b (Φ_ST = 0.361, P < 0.0001; Φ_CT = 0.312, P < 0.02) and control region (Φ_ST = 0.166, P < 0.0001; Φ_CT = 0.128, P < 0.03) sequences. However, phylogenetic reconstructions failed to show reciprocal monophyly in populations between ocean basins. Microsatellite analyses revealed significant population substructuring between all locations sampled except for the two locations that were in closest proximity to each other (global F_ST value = 0.026). Bayesian analysis of microsatellite data also revealed significant population structuring between ocean basins. Coalescent‐based analyses of microsatellite data revealed low migration rates among all sites. Mismatch distribution analysis of mitochondrial loci supports a sudden population expansion in both ocean basins in the late Pleistocene, with the expansion of Atlantic populations occurring more recently. Main conclusions Present‐day populations of S. floridae do not bear the mitochondrial DNA signature of the strong phylogenetic discontinuity between the Atlantic and Gulf coasts of North America commonly observed in other species. Rather, our results suggest that Atlantic and Gulf of Mexico populations of S. floridae are closely related but nevertheless exhibit local and regional population structure. We conclude that the present‐day phylogeographic pattern is the result of a recent population expansion into the Atlantic in the late Pleistocene, and that life‐history traits and ecology may play a pivotal role in shaping the realized geographical distribution pattern of this species.     

QST < FST As a signature of canalization

A key aim of evolutionary biology – inferring the action of natural selection on wild species – can be achieved by comparing neutral genetic differentiation between populations (F_ST) with quantitative genetic variation (Q_ST). Each of the three possible outcomes of comparisons of Q_ST and F_ST (Q_ST > F_ST, Q_ST = F_ST, Q_ST < F_ST) is associated with an inference (diversifying selection, genetic drift, uniform selection, respectively). However, published empirical and theoretical studies have focused on the Q_ST > F_ST outcome. We believe that this reflects the absence of a straightforward biological interpretation of the Q_ST < F_ST pattern. We here report recent evidence of this neglected evolutionary pattern, provide guidelines to its interpretation as either a canalization phenomenon or a consequence of uniform selection and discuss the significant importance this issue will have for the area of evolutionary biology.     

Population structure and cryptic speciation in bonnethead sharks Sphyrna tiburo in the south‐eastern U.S.A. and Caribbean

Population structure and lineage diversification within a small, non‐dispersive hammerhead shark species, the bonnethead shark Sphyrna tiburo, was assessed. Sphyrna tiburo is currently described as one continuously distributed species along the Atlantic continental margins of North, Central and South America, but recent genetic analysis of an insular population (Trinidad) suggests the possibility of cryptic speciation. To address this issue S. tiburo were sampled at six sites along c. 6200 km of continuous, continental coastline and from one island location (Grand Bahama) across a discontinuity in their distribution (the Straits of Florida), in order to test if they constitute a single lineage over this distribution. A total of 1030 bp of the mitochondrial control region (CR) was obtained for 239 S. tiburo, revealing 73 distinct haplotypes, high nucleotide diversity (0·01089) and a pair of highly divergent lineages estimated to have separated 3·61–5·62 million years ago. Mitochondrial cytochrome oxidase I and nuclear internal transcribed spacer loci show the same pattern. Divergence is similar within S. tiburo to that observed between established elasmobranch sister species, providing further evidence of cryptic speciation. A global AMOVA based on CR confirms that genetic diversity is primarily partitioned among populations (Φ_ST = 0·828, P < 0·001) because the divergent lineages are almost perfectly segregated between Belize and North America–The Bahamas. An AMOVA consisting solely of the North American and Bahamian samples is also significantly different from zero (Φ_ST = 0·088, P < 0·001) and pairwise F_ST is significantly different between all sites. These findings suggest that S. tiburo comprises a species complex and supports previous research indicating fine population structure, which has implications for fisheries management and biodiversity conservation.     

Limited genetic variation and structure in softshell clams (Mya arenaria) across their native and introduced range

To offset declines in commercial landings of the softshell clam, Mya arenaria, resource managers are engaged in extensive stocking of seed clams throughout its range in the northwest Atlantic. Because a mixture of native and introduced stocks can disrupt locally adapted genotypes, we investigated genetic structure in M. arenaria populations across its current distribution to test for patterns of regional differentiation. We sequenced mitochondrial cytochrome oxidase I for a total of 212 individuals from 12 sites in the northwest Atlantic (NW Atlantic), as well as two introduced sites, the northeast Pacific (NE Pacific), and the North Sea Europe (NS Europe). Populations exhibited extremely low genetic variation, with one haplotype dominating (65–100%) at all sites sampled. Despite being introduced in the last 150–400 years, both NE Pacific and NS Europe populations had higher diversity measures than those in the NW Atlantic and both contained private haplotypes at frequencies of 10–27% consistent with their geographic isolation. While significant genetic structure (F _ST = 0.159, P < 0.001) was observed between NW Atlantic and NS Europe, there was no evidence for genetic structure across the pronounced environmental clines of the NW Atlantic. Reduced genetic diversity in mtDNA combined with previous studies reporting reduced genetic diversity in nuclear markers strongly suggests a recent population expansion in the NW Atlantic, a pattern that may result from the retreat of ice sheets during Pleistocene glacial periods. Lack of genetic diversity and regional genetic differentiation suggests that present management strategies for the commercially important softshell clam are unlikely to have a significant impact on the regional distribution of genetic variation, although the possibility of disrupting locally adapted stocks cannot be excluded.     

Diversity and genetic connectivity among populations of a threatened tree (<Emphasis Type="Italic">Dalbergia nigra</Emphasis>) in a recently fragmented landscape of the Brazilian Atlantic Forest

In this study we evaluated the influence of recent landscape fragmentation on the dynamics of remnant fragments from the Brazilian Atlantic Forest. This biome is one of the richest in the world and has been extensively deforested and fragmented. We sampled five populations of the threatened Dalbergia nigra, a tree endemic to the Brazilian Atlantic Forest, two located in a large reserve of continuous forest and three in fragments of different sizes and levels of disturbance. In order to assess historical changes, considering the longevity of the analyzed species, 119 adults and 116 saplings were genotyped for six microsatellite loci. Lower levels of genetic diversity were found in the most impacted fragments when compared to the most preserved population located inside the reserve, and there was significant genetic structure among the populations studied (pairwise F _ST = 0.031–0.152; pairwise D _EST = 0.039–0.301). However, genetic structure among saplings (F _ST = 0.056; D _EST = 0.231) was significantly lower than among adults (F _ST = 0.088; D _EST = 0.275). Estimates of contemporary gene flow based on assignment tests corroborated this result, suggesting that fragmentation led to an increase in gene flow. This connectivity among remnant fragments could mitigate the loss of genetic diversity through a metapopulation dynamic, but the high rate of habitat loss and the unknown long-term genetic effects add uncertainty. These results, taken together with the presence of private alleles in disturbed populations, highlight the importance of preserving the extant fragments.     

Chromosomal fusion and life history‐associated genomic variation contribute to within‐river local adaptation of Atlantic salmon

Chromosomal inversions have been implicated in facilitating adaptation in the face of high levels of gene flow, but whether chromosomal fusions also have similar potential remains poorly understood. Atlantic salmon are usually characterized by population structure at multiple spatial scales; however, this is not the case for tributaries of the Miramichi River in North America. To resolve genetic relationships between populations in this system and the potential for known chromosomal fusions to contribute to adaptation, we genotyped 728 juvenile salmon using a 50 K SNP array. Consistent with previous work, we report extremely weak overall population structuring (Global F_ST = 0.004) and failed to support hierarchical structure between the river's two main branches. We provide the first genomic characterization of a previously described polymorphic fusion between chromosomes 8 and 29. Fusion genomic characteristics included high LD, reduced heterozygosity in the fused homokaryotes, and strong divergence between the fused and the unfused rearrangement. Population structure based on fusion karyotype was five times stronger than neutral variation (F_ST = 0.019), and the frequency of the fusion was associated with summer precipitation supporting a hypothesis that this rearrangement may contribute local adaptation despite weak neutral differentiation. Additionally, both outlier variation among populations and a polygenic framework for characterizing adaptive variation in relation to climate identified a 250‐Kb region of chromosome 9, including the gene six6 that has previously been linked to age‐at‐maturity and run‐timing for this species. Overall, our results indicate that adaptive processes, independent of major river branching, are more important than neutral processes for structuring these populations.     

Genetically distinct populations of northern shrimp,Pandalus borealis,in the North Atlantic: adaptation to different temperatures as an isolation factor

The large‐scale population genetic structure of northern shrimp, Pandalus borealis, was investigated over the species’ range in the North Atlantic, identifying multiple genetically distinct groups. Genetic divergence among sample localities varied among 10 microsatellite loci (range: F_ST = ?0.0002 to 0.0475) with a highly significant average (F_ST = 0.0149; P < 0.0001). In contrast, little or no genetic differences were observed among temporal replicates from the same localities (F_ST = 0.0004; P = 0.33). Spatial genetic patterns were compared to geographic distances, patterns of larval drift obtained through oceanographic modelling, and temperature differences, within a multiple linear regression framework. The best‐fit model included all three factors and explained approximately 29% of all spatial genetic divergence. However, geographic distance and larval drift alone had only minor effects (2.5–4.7%) on large‐scale genetic differentiation patterns, whereas bottom temperature differences explained most (26%). Larval drift was found to promote genetic homogeneity in parts of the study area with strong currents, but appeared ineffective across large temperature gradients. These findings highlight the breakdown of gene flow in a species with a long pelagic larval phase (up to 3 months) and indicate a role for local adaptation to temperature conditions in promoting evolutionary diversification and speciation in the marine environment.     

Post‐glacial colonization routes coincide with a life‐history breakpoint along a latitudinal gradient

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post‐glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical Q_ST–F_ST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The Q_ST estimates between the colonization routes were four times higher than F_ST, indicating a prominent role for natural selection. Q_ST within the colonization routes did not generally differ from F_ST, but we found temperature‐dependent adaptive divergence close to the contact zone. These results indicate that lineage‐specific variation can account for much of the adaptive divergence along a latitudinal gradient.     

Size differentiation in Finnish house sparrows follows Bergmann's rule with evidence of local adaptation

Bergmann's rule predicts that individuals are larger in more poleward populations and that this size gradient has an adaptive basis. Hence, phenotypic divergence in size traits between populations (P_ST) is expected to exceed the level of divergence by drift alone (F_ST). We measured 16 skeletal traits, body mass and wing length in 409 male and 296 female house sparrows Passer domesticus sampled in 12 populations throughout Finland, where the species has its northernmost European distributional margin. Morphometric differentiation across populations (P_ST) was compared with differentiation in 13 microsatellites (F_ST). We find that twelve traits phenotypically diverged more than F_ST in both sexes, and an additional two traits diverged in males. The phenotypic divergence exceeded F_ST in several traits to such a degree that findings were robust also to strong between‐population environmental effects. Divergence was particularly strong in dimensions of the bill, making it a strong candidate for the study of adaptive molecular genetic divergence. Divergent traits increased in size in more northern populations. We conclude that house sparrows show evidence of an adaptive latitudinal size gradient consistent with Bergmann's rule on the modest spatial scale of ca. 600 km.