Continuous Gene Flow Research Articles

Speciation with gene flow in a narrow endemic West Virginia cave salamander (Gyrinophilus subterraneus)

Due to their limited geographic distributions and specialized ecologies, cave species are often highly endemic and can be especially vulnerable to habitat degradation within and surrounding the cave systems they inhabit. We investigated the evolutionary history of the West Virginia Spring Salamander (Gyrinophilus subterraneus), estimated the population trend from historic and current survey data, and assessed the current potential for water quality threats to the cave habitat. Our genomic data (mtDNA sequence and ddRADseq-derived SNPs) reveal two, distinct evolutionary lineages within General Davis Cave corresponding to G. subterraneus and its widely distributed sister species, Gyrinophilus porphyriticus, that are also differentiable based on morphological traits. Genomic models of evolutionary history strongly support asymmetric and continuous gene flow between the two lineages, and hybrid classification analyses identify only parental and first generation cross (F1) progeny. Collectively, these results point to a rare case of sympatric speciation occurring within the cave, leading to strong support for continuing to recognize G. subterraneus as a distinct and unique species. Due to its specialized habitat requirements, the complete distribution of G. subterraneus is unresolved, but using survey data in its type locality (and currently the only known occupied site), we find that the population within General Davis Cave has possibly declined over the last 45 years. Finally, our measures of cave and surface stream water quality did not reveal evidence of water quality impairment and provide important baselines for future monitoring. In addition, our unexpected finding of a hybrid zone and partial reproductive isolation between G. subterraneus and G. porphyriticus warrants further attention to better understand the evolutionary and conservation implications of occasional hybridization between the species.

Divergence With Gene Flow and Contrasting Population Size Blur the Species Boundary in Cycas Sect. Asiorientales, as Inferred From Morphology and RAD-Seq Data.

The divergence process of incipient species is fascinating but elusive by incomplete lineage sorting or gene flow. Species delimitation is also challenging among those morphologically similar allopatric species, especially when lacking comprehensive data. Cycas sect. Asiorientales, comprised of C. taitungensis and C. revoluta in the Ryukyu Archipelago and Taiwan, diverged recently with continuous gene flow, resulting in a reciprocal paraphyletic relationship. Their previous evolutionary inferences are questioned from few genetic markers, incomplete sampling, and incomprehensive morphological comparison by a long-term taxonomic misconception. By whole range sampling, this study tests the geographic mode of speciation in the two species of Asiorientales by approximate Bayesian computation (ABC) using genome-wide single nucleotide polymorphisms (SNPs). The individual tree was reconstructed to delimit the species and track the gene-flow trajectory. With the comparison of diagnostic morphological traits and genetic data, the allopatric speciation was rejected. Alternatively, continuous but spatially heterogeneous gene flow driven by transoceanic vegetative dispersal and pollen flow with contrasting population sizes blurred their species boundary. On the basis of morphological, genetic, and evolutionary evidence, we synonymized these two Cycas species. This study highlights not only the importance of the Kuroshio Current to species evolution but also the disadvantage of using species with geographically structured genealogies as conservation units.

Whole-genome analysis of multiple wood ant population pairs supports similar speciation histories, but different degrees of gene flow, across their European ranges.

The application of demographic history modelling and inference to the study of divergence between species has become a cornerstone of speciation genomics. Speciation histories are usually reconstructed by analysing single populations from each species, assuming that the inferred population history represents the actual speciation history. However, this assumption may not be met when species diverge with gene flow, for example, when secondary contact may be confined to specific geographic regions. Here, we tested whether divergence histories inferred from heterospecific populations may vary depending on their geographic locations, using the two wood ant species Formica polyctena and F. aquilonia. We performed whole‐genome resequencing of 20 individuals sampled in multiple locations across the European ranges of both species. Then, we reconstructed the histories of distinct heterospecific population pairs using a coalescent‐based approach. Our analyses always supported a scenario of divergence with gene flow, suggesting that divergence started in the Pleistocene (c. 500 kya) and occurred with continuous asymmetrical gene flow from F. aquilonia to F. polyctena until a recent time, when migration became negligible (2–19 kya). However, we found support for contemporary gene flow in a sympatric pair from Finland, where the species hybridise, but no signature of recent bidirectional gene flow elsewhere. Overall, our results suggest that divergence histories reconstructed from a few individuals may be applicable at the species level. Nonetheless, the geographical context of populations chosen to represent their species should be taken into account, as it may affect estimates of migration rates between species when gene flow is spatially heterogeneous.

Wild sunflower goes viral: Citizen science and comparative genomics allow tracking the origin and establishment of invasive sunflower in the Levant.

Globalization and intensified volume of trade and transport around the world are accelerating the rate of biological invasions. It is therefore increasingly important to understand the processes through which invasive species colonize new habitats, often to the detriment of native flora. The initial steps of an invasion are particularly critical, as the introduced species relies on limited genetic diversity to adapt to a new environment. However, our understanding of this critical stage of the invasion is currently limited. We used a citizen science approach and social media to survey the distribution of invasive sunflower in Israel. We then sampled and sequenced a representative collection and compared it with available genomic data sets of North American wild sunflower, landraces and cultivars. We show that invasive wild sunflower is rapidly establishing throughout Israel, probably from a single, recent introduction from Texas, while maintaining high genetic diversity through ongoing gene flow. Since its introduction, invasive sunflower has spread quickly to most regions, and differentiation was detected despite extensive gene flow between clusters. Our findings suggest that rapid spread followed by continuous gene flow between diverging populations can serve as an efficient mechanism for maintaining sufficient genetic diversity at the early stages of invasion, promoting rapid adaptation and establishment in the new territory.

Microsatellite Markers from Peronospora tabacina, the Cause of Blue Mold of Tobacco, Reveal Species Origin, Population Structure, and High Gene Flow.

Peronospora tabacina is an obligate parasite that causes blue mold of tobacco. The pathogen reproduces primarily by sporangia, whereas the sexual oospores are rarely observed. A collection of 122 isolates of P. tabacina was genotyped using nine microsatellites to assess the population structure of individuals from subpopulations collected from central, southern, and western Europe; the Middle East; Central America; North America; and Australia. Genetic variations among the six subpopulations accounted for ∼8% of the total variation, including moderate levels of genetic differentiation, high gene flow among these subpopulations, and a positive correlation between geographic and genetic distance (r = 0.225; P < 0.001). Evidence of linkage disequilibrium (P < 0.001) showed that populations contained partially clonal subpopulations but that subpopulations from Australia and Mediterranean Europe did not. High genetic variation and population structure among samples could be explained by continuous gene flow across continents via infected transplant exchange and/or long-distance dispersal of sporangia via wind currents. This study analyzed the most numerous P. tabacina collection and allowed conclusions regarding the migration, mutation, and evolutionary history of this obligate biotrophic oomycete. The evidence pointed to the species origin in Australia and identified intracontinental and intercontinental migration patterns of this important pathogen.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Invasion genetics of the Asian hornet Vespa velutina nigrithorax in Southern Europe

In 2004, Vespa velutina was first seen in France. Since then, this fierce honey bee predator spread across many countries, giving rise to one of the most phenomenal insect invasions in Europe. An early study in France showed a genetically depauperate population, originating from a single multi-mated queen introduced from China. Here, we further unveil V. velutina invasion genetics in Europe by surveying the Iberian and Italian peninsulas using cytonuclear markers. Our results show that the French population acted as the colonists’ source in Spain, Portugal and Italy, leading to rejecting the hypothesis of multiple introductions from the native range. While Spain and Italy were colonized predominantly by leading-edge expansions from the French core population, in Portugal the invasion started from long-distance jump. Both processes were accompanied by a significant reduction in genetic diversity, with stronger losses for Portugal (Ar = 17.4%; uHe = 42.3%) than for Spain (Ar = 9.0%; uHe = 20.6%) or Italy (Ar = 16.3%; uHe = 26.8%). Signatures of differentiation and population structure, associated to the founding event in Portugal, enabled detection of secondary contact between the front derived from the primary propagule introduced in France and the front derived from the secondary propagule introduced in Portugal. Detection of first-generation migrants in the three countries suggests continuous gene flow that is bringing in new alleles, and this effect is stronger in Portugal, as reflected by a 20.3% increase in allelic richness. Overall, this study provides further insights into the invasion genetics of V. velutina in Europe, which can aid developing strategies to manage this major threat to beekeeping.

Plant diversity and ecology on the Qinghai–Tibet Plateau

Plant diversity and ecology on the Qinghai–Tibet Plateau

Genome evolution of the psammophyte Pugionium for desert adaptation and further speciation

Deserts exert strong selection pressures on plants, but the underlying genomic drivers of ecological adaptation and subsequent speciation remain largely unknown. Here, we generated de novo genome assemblies and conducted population genomic analyses of the psammophytic genus Pugionium (Brassicaceae). Our results indicated that this bispecific genus had undergone an allopolyploid event, and the two parental genomes were derived from two ancestral lineages with different chromosome numbers and structures. The postpolyploid expansion of gene families related to abiotic stress responses and lignin biosynthesis facilitated environmental adaptations of the genus to desert habitats. Population genomic analyses of both species further revealed their recent divergence with continuous gene flow, and the most divergent regions were found to be centered on three highly structurally reshuffled chromosomes. Genes under selection in these regions, which were mainly located in one of the two subgenomes, contributed greatly to the interspecific divergence in microhabitat adaptation.

Honey Bee Habitat Sharing Enhances Gene Flow of the Parasite Nosema ceranae.

Host-parasite co-evolution is a process of reciprocal, adaptive genetic change. In natural conditions, parasites can shift to other host species, given both host and parasite genotypes allow this. Even though host-parasite co-evolution has been extensively studied both theoretically and empirically, few studies have focused on parasite gene flow between native and novel hosts. Nosema ceranae is a native parasite of the Asian honey bee Apis cerana, which infects epithelial cells of mid-guts. This parasite successfully switched to the European honey bee Apis mellifera, where high virulence has been reported. In this study, we used the parasite N. ceranae and both honey bee species as model organisms to study the impacts of two-host habitat sharing on parasite diversity and virulence. SNVs (Single Nucleotide Variants) were identified from parasites isolated from native and novel hosts from sympatric populations, as well as novel hosts from a parapatric population. Parasites isolated from native hosts showed the highest levels of polymorphism. By comparing the parasites isolated from novel hosts between sympatric and parapatric populations, habitat sharing with the native host significantly enhanced parasite diversity, suggesting there is continuing gene flow of parasites between the two host species in sympatric populations.

Genetic structure of Fusarium verticillioides populations from maize in Iran

Fusarium verticillioides is one of the most important fungal pathogens of maize since it causes severe yield losses and produces the mycotoxins fumonisins that represent a major concern for human and animal health. Information about genetic diversity and population structure of fungal pathogens is essential for developing disease management strategies. The aim of this research was to investigate the genetic structure of F. verticillioides isolated from different provinces of Iran through determination of mating type idiomorphs, phylogenetic analyses based on translation elongation factor-1 alpha (EF-1α), RNA Polymerase II Subunit (RPB2), beta-tubulin (tub2) and Calmodulin (cmdA) genes and genetic diversity analyses based on 6 simple-sequence repeats (SSRs). Both mating types were detected in Iranian populations of F. verticillioides, particularly in Qazvin and Khuzestan, with equal frequency, which highlighted that sexual reproduction is favorable under field conditions. However, the linkage disequilibrium indices did not support the hypothesis of random mating in Khuzestan and Fars. Although assessment of nucleotide diversity based on housekeeping genes showed low level of variation among strains, genotype diversity based on SSRs revealed a high level of genetic diversity within Iranian populations. AMOVA analysis highlighted that the genetic variation of F. verticillioides in Iran was mainly distributed within population of a single area (97%), while a small proportion of genetic variation (3%) resided among populations. These patterns of variation are likely explained by the continuous gene flow among populations isolated from different areas. On the other hand, principal coordinate analysis indicated that the distribution of genetic variation among populations could be explained by the geographical distances. Consequently, to reduce pathogen gene flow among regions, the quarantine processes in Iran should be intensified.

Compensation of Wild Plants Weakens the Effects of Crop-Wild Gene Flow on Wild Rice Populations.

Crop-wild gene flow may alter the fitness of the recipient i.e., crop-wild hybrids, then potentially impact wild populations, especially for the gene flow carrying selective advantageous crop alleles, such as transgenes conferring insect resistance. Given the continuous crop-wild gene flow since crop domestication and the occasionally stressful environments, the extant wild populations of most crops are still “wild.” One interpretation for this phenomenon is that wild populations have the mechanism buffered for the effects of crop alleles. However, solid evidence for this has been scarce. We used wild rice (Oryza rufipogon) and transgenic (Bt/CpTI) rice (O. sativa) as a crop-wild gene flow model and established cultivated, wild, and F7 hybrid rice populations under four levels of insect (Chilo suppressalis) pressure. Then, we measured the trait performance of the plants and estimated fitness to test the compensatory response of relatively high fitness compared to the level of insect damage. The performance of all plants varied with the insect pressure level; wild plants had higher insect-tolerance that was expressed as over- or equal-compensatory responses to insect damage, whereas crop and hybrids exhibited under-compensatory responses. The higher compensation resulted in a better performance of wild rice under insect pressure where transgenes conferring insect resistance had a somewhat beneficial effect. Remarkable hybrid vigour and the benefit effect of transgenes increased the fitness of hybrids together, but this joint effect was weakened by the compensation of wild plants. These results suggest that compensation to environmental stress may reduce the potential impacts of crop alleles on wild plants, thereby it is a mechanism maintaining the “wild” characteristics of wild populations under the scenario of continuous crop-wild gene flow.

Phylogeography of Prunus armeniaca L. revealed by chloroplast DNA and nuclear ribosomal sequences

To clarify the phytogeography of Prunus armeniaca L., two chloroplast DNA fragments (trnL-trnF and ycf1) and the nuclear ribosomal DNA internal transcribed spacer (ITS) were employed to assess genetic variation across 12 P. armeniaca populations. The results of cpDNA and ITS sequence data analysis showed a high the level of genetic diversity (cpDNA: HT = 0.499; ITS: HT = 0.876) and a low level of genetic differentiation (cpDNA: FST = 0.1628; ITS: FST = 0.0297) in P. armeniaca. Analysis of molecular variance (AMOVA) revealed that most of the genetic variation in P. armeniaca occurred among individuals within populations. The value of interpopulation differentiation (NST) was significantly higher than the number of substitution types (GST), indicating genealogical structure in P. armeniaca. P. armeniaca shared genotypes with related species and may be associated with them through continuous and extensive gene flow. The haplotypes/genotypes of cultivated apricot populations in Xinjiang, North China, and foreign apricot populations were mixed with large numbers of haplotypes/genotypes of wild apricot populations from the Ili River Valley. The wild apricot populations in the Ili River Valley contained the ancestral haplotypes/genotypes with the highest genetic diversity and were located in an area considered a potential glacial refugium for P. armeniaca. Since population expansion occurred 16.53 kyr ago, the area has provided a suitable climate for the population and protected the genetic diversity of P. armeniaca.

Different population size change and migration histories created genetic diversity of three oaks in Tokai region, central Japan.

To understand genetic diversity in focal species, it is important to consider the possibility of speciation with gene flow, especially in species with porous genomes such as oaks. We studied genetic diversity and structure in three oaks, Quercus mongolica var. mongolicoides (QM), Q. mongolica var. crispula (QC) and Q. serrata (QS), growing in the Tokai region, central Japan. QM is semi-endemic to the region while the others are common taxa. We also conducted demographic modeling to infer their population size change and migration histories using an approximate Bayesian computation (ABC) approach. The three taxa showed distinct genetic structures but there was genetic admixture among the taxa, especially between QM and QC. ABC analysis of population size change revealed that the population size of QM was stable during and after the last glacial period, while QC and QS showed population expansion after the last glacial maximum. ABC analysis of population divergence and migration revealed that continuous gene flow between QM and QC after their divergence was supported, while between QM and QS, and between QC and QS, secondary contact after sufficient isolation was supported. These historical migration patterns among the three taxa indicate that QM and QC are currently in the early stage or gray zone of speciation, whereas speciation of the other two taxon pairs is considered to have almost been established. Observed gene flow patterns and strength between QM and QC, and between QM and QS, were explained by both flowering patterns and historical distributions, but those between QC and QS were not.

Low Persistence of Genetic Rescue Across Generations in the Arctic Fox (Vulpes lagopus).

Genetic rescue can facilitate the recovery of small and isolated populations suffering from inbreeding depression. Long-term effects are however complex, and examples spanning over multiple generations under natural conditions are scarce. The aim of this study was to test for long-term effects of natural genetic rescue in a small population of Scandinavian Arctic foxes (Vulpes lagopus). By combining a genetically verified pedigree covering almost 20 years with a long-term dataset on individual fitness (n = 837 individuals), we found no evidence for elevated fitness in immigrant F2 and F3 compared to native inbred foxes. Population inbreeding levels showed a fluctuating increasing trend and emergence of inbreeding within immigrant lineages shortly after immigration. Between 0–5 and 6–9 years post immigration, the average number of breeding adults decreased by almost 22% and the average proportion of immigrant ancestry rose from 14% to 27%. Y chromosome analysis revealed that 2 out of 3 native male lineages were lost from the gene pool, but all founders represented at the time of immigration were still contributing to the population at the end of the study period through female descendants. The results highlight the complexity of genetic rescue and suggest that beneficial effects can be brief. Continuous gene flow may be needed for small and threatened populations to recover and persist in a longer time perspective.

Inbreeding depression contributes to the maintenance of habitat segregation between closely related monkeyflower species.

Incompletely reproductively isolated species often segregate into different microhabitats, even when they are able to survive and reproduce in both habitats. Longer term evolutionary factors may contribute to this lack of cross-habitat persistence. When reproductive interference reduces immigrant fitness, assortative mating, including self-fertilization, increases immigrants' fitness in a single generation, but longer term, inbreeding depression may reduce the chance of population persistence. Two California monkeyflower species repeatedly segregate into drier and wetter areas in their zone of sympatry. To test whether inbreeding depression may contribute to the maintenance of this segregation pattern, we transplanted outbred and successively inbred Mimulus guttatus and Mimulus nudatus into their native habitats and heterospecific habitats. We measured germination, survival, and seed set and found that recurrent selfing reduced all aspects of fitness in both species, most strongly in foreign habitats. A simulation model, parameterized from the transplant experiment, found that inbreeding reduced fitness to such an extent that sequentially inbred populations of either species would be unable to persist in heterospecific-occupied habitats in the absence of continued gene flow. These results demonstrate that individual immigrants are unlikely to form persistent populations and thus, inbreeding depression contributes to the absence of fine-scale coexistence in this species pair.

No severe genetic bottleneck in a rapidly range-expanding bumblebee pollinator.

Genetic bottlenecks can limit the success of populations colonizing new ranges. However, successful colonizations can occur despite bottlenecks, a phenomenon known as the genetic paradox of invasion. Eusocial Hymenoptera such as bumblebees (Bombus spp.) should be particularly vulnerable to genetic bottlenecks, since homozygosity at the sex-determining locus leads to costly diploid male production (DMP). The Tree Bumblebee (Bombus hypnorum) has rapidly colonized the UK since 2001 and has been highlighted as exemplifying the genetic paradox of invasion. Using microsatellite genotyping, combined with the first genetic estimates of DMP in UK B. hypnorum, we tested two alternative genetic hypotheses (‘bottleneck’ and ‘gene flow’ hypotheses) for B. hypnorum's colonization of the UK. We found that the UK population has not undergone a recent severe genetic bottleneck and exhibits levels of genetic diversity falling between those of widespread and range-restricted Bombus species. Diploid males occurred in 15.4% of reared colonies, leading to an estimate of 21.5 alleles at the sex-determining locus. Overall, the findings show that this population is not bottlenecked, instead suggesting that it is experiencing continued gene flow from the continental European source population with only moderate loss of genetic diversity, and does not exemplify the genetic paradox of invasion.

Inconsistency Between Socio-Spatial and Genetic Structure in a Coastal Dolphin Population

Identifying population structure and boundaries among communities of wildlife exposed to anthropogenic threats is key to successful conservation management. Previous studies on the demography, social and spatial structure of Indo-Pacific bottlenose dolphins (Tursiops aduncus) suggested four nearly discrete behavioral communities in Perth metropolitan waters, Western Australia. We investigated the genetic structure of these four communities using highly polymorphic microsatellite markers and part of the hypervariable segment of the mitochondrial control region. Overall, there was no evidence of spatial genetic structure. We found significant, yet very small genetic differentiation between some communities, most likely due to the presence of highly related individuals within these communities. Our findings of high levels of contemporary migration and highly related individuals among communities point toward a panmictic genetic population with continuous gene flow among each of the communities. In species with slow life histories and fission-fusion dynamics, such asTursiopsspp., genetic and socio-spatial structures may reflect different timescales. Thus, despite genetic similarity, each social community should be considered as a distinct ecological unit to be conserved because they are exposed to different anthropogenic threats and occur in different ecological habitats, social structure being as important as genetic information for immediate conservation management. The estuarine community, in particular, is highly vulnerable and appropriate conservation measures are needed in order to maintain its connectivity with the adjacent, semi-enclosed coastal communities.

Isolation by environment and recurrent gene flow shaped the evolutionary history of a continentally distributed Neotropical treefrog

AbstractAimPhylogeographic studies show how historical and current changes in landscapes shape the geographic distribution of genetic diversity in species of animals and plants. In particular, for the species of the Diagonal of Open Formations (DOF), the compartmentalization of the Central Brazilian Plateau (CBP) during the Tertiary and climatic oscillations during the Quaternary have often been invoked to explain the origin and current patterns of biodiversity. We investigated how landscape changes and climatic oscillations shaped the distribution and diversification history of a widespread South American treefrog.LocationSouth American Diagonal of Open Formations (DOF) including Caatinga, Cerrado, and Chaco biomes.TaxonTreefrog Boana raniceps.MethodsWe used a multi‐locus dataset from 288 individual frogs collected at 115 localities throughout most of the species’ distribution. We used population assignment analysis, species distribution models, historical demography models, approximate Bayesian computation and landscape genetic analyses to test alternative hypotheses of diversification.ResultsWe found two genetic lineages that diverged during the mid‐Pleistocene with continued gene flow. Approximate Bayesian computation supported a scenario of isolation with migration until the Last Glacial Maximum, followed by more recent population expansion in north‐eastern Brazil and stability at the southwest in South America. Isolation by environment was the best predictor of genetic distance between populations, which is in accordance with their different environmental niches. As Boana raniceps is a lowland species, steep slopes in the CBP likely restrained gene flow enough to sustain population divergence. We found evidence for major range contraction during the Last Glacial Maximum, raising the possibility of synergic action of climate change and the CBP compartmentalization in regulating migration.Main conclusionsOur findings highlight how landscape and climatic changes can shape the diversification of DOF biota. Past climatic fluctuations and environmental resistance due to topography acted in concert, forming a semipermeable barrier to gene flow, promoting intraspecific differentiation in a continentally distributed species.

Population history of North Africa based on modern and ancient genomes.

Compared with the rest of the African continent, North Africa has provided limited genomic data. Nonetheless, the genetic data available show a complex demographic scenario characterized by extensive admixture and drift. Despite the continuous gene flow from the Middle East, Europe and sub-Saharan Africa, an autochthonous genetic component that dates back to pre-Holocene times is still present in North African groups. The comparison of ancient and modern genomes has evidenced a genetic continuity in the region since Epipaleolithic times. Later population movements, especially the gene flow from the Middle East associated with the Neolithic, have diluted the genetic autochthonous component, creating an east to west gradient. Recent historical movements, such as the Arabization, have also contributed to the genetic landscape observed currently in North Africa and have culturally transformed the region. Genome analyses have not shown evidence of a clear correlation between cultural and genetic diversity in North Africa, as there is no genetic pattern of differentiation between Tamazight (i.e. Berber) and Arab speakers as a whole. Besides the gene flow received from neighboring areas, the analysis of North African genomes has shown that the region has also acted as a source of gene flow since ancient times. As a result of the genetic uniqueness of North African groups and the lack of available data, there is an urgent need for the study of genetic variation in the region and its implications in health and disease.

Gene flow relates to evolutionary divergence among populations at the range margin.

BackgroundMorphological differentiation between populations resulting from local adaptations to environmental conditions is likely to be more pronounced in populations with increasing genetic isolation. In a previous study a positive clinal variation in body size was observed in isolated Roesel’s bush-cricket, Metrioptera roeselii, populations, but were absent from populations within a continuous distribution at the same latitudinal range. This observational study inferred that there was a phenotypic effect of gene flow on climate-induced selection in this species.MethodsTo disentangle genetic versus environmental drivers of population differences in morphology, we measured the size of four different body traits in wild-caught individuals from the two most distinct latitudinally-matched pairs of populations occurring at about 60°N latitude in northern Europe, characterised by either restricted or continuous gene flow, and corresponding individuals raised under laboratory conditions.ResultsIndividuals that originated from the genetically isolated populations were always bigger (femur, pronotum and genital appendages) when compared to individuals from latitudinally-matched areas characterised by continuous gene flow between populations. The magnitude of this effect was similar for wild-caught and laboratory-reared individuals. We found that previously observed size cline variation in both male and female crickets was likely to be the result of local genetic adaptation rather than phenotypic plasticity.ConclusionsThis strongly suggests that restricted gene flow is of major importance for frequencies of alleles that participate in climate-induced selection acting to favour larger phenotypes in isolated populations towards colder latitudes.