Repeated sampling detects gene flow in a flightless ground beetle in a fragmented landscape

Secondary clines level down in the course of time if the gene flow is not interrupted. Temporally repeated sampling of populations in a cline allows the investigation not only of its occurrence but also of the estimation of the amount of ongoing gene flow. We reinvestigated an allozyme gradient in Carabus auronitens populations in the Westphalian Lowlands (northwestern Germany) 15 to 20 years after it had originally been recorded. A total of 977 individuals of this flightless woodland species from 29 sample sites were genotyped at the diallelic Est-1 locus in 2005-2006 and compared to former findings, collected in 1985-1994 from the same populations. Both data sets showed clinal variation. Pairwise differences between the samples of both data sets indicated significant decrease in the steepness of the cline during the past 15 to 20 years. The estimated average gene flow per generation is 0.6% of each beetle population. Ongoing gene flow in the flightless ground beetle C. auronitens led to a less pronounced cline despite a stable degree of fragmentation (and connectivity) of the landscape. Migration and gene flow were obviously enabled by the numerous hedgerows. The corridors are seen to be a prerequisite for migration between populations and for possible future range shifts of forest insect species.     

The influence of gene flow on latitudinal clines of photoperiodic adult diapause in the Drosophila auraria species-complex

Temperate species belonging to the Drosophila auraria species complex, D. auraria Peng, D. biauraria Bock & Wheeler, D. triauraria Bock & Wheeler and D. subauraria Kimura, enter reproductive diapause to pass the winter in response to short daylengths. These species from Japan showed latitudinal clines in critical daylength which is longer in populations from higher latitudes. The slopes of these clines coincided well with that of the cline which is approximately predicted from climatic data, suggesting that these clines result from adaptation of the species to the latitudinal gradient of climatic conditions. Between the mainlands and the surrounding islands of Japan, the slopes of clines did not differ significantly, but the deviation from the regression line was usually smaller in mainland populations. It is assumed that gene flow reduces the genetic variation among mainland populations and results in the development of smooth clines. In the plain of east China, D. triauraria did not show clinal variation in critical daylength, although the development of the cline is expected from climatic data. Extensive gene flow among Chinese populations is considered to prevent the development of a cline.     

Adaptation to a seasonally varying environment: a strong latitudinal cline in reproductive diapause combined with high gene flow in Drosophila montana

Adaptation to seasonal changes in the northern hemisphere includes an ability to predict the forthcoming cold season from gradual changes in environmental cues early enough to prepare for the harsh winter conditions. The magnitude and speed of changes in these cues vary between the latitudes, which induces strong selection pressures for local adaptation.We studied adaptation to seasonal changes in Drosophila montana, a northern maltfly, by defining the photoperiodic conditions leading to adult reproductive diapause along a latitudinal cline in Finland and by measuring genetic differentiation and the amount of gene flow between the sampling sites with microsatellites. Our data revealed a clear correlation between the latitude and the critical day length (CDL), in which half of the females of different cline populations enter photoperiodic reproductive diapause. There was no sign of limited gene flow between the cline populations, even though these populations showed isolation by distance. Our results show that local adaptation may occur even in the presence of high gene flow, when selection for locally adaptive life-history traits is strong. A wide range of variation in the CDLs of the fly strains within and between the cline populations may be partly due to gene flow and partly due to the opposing selection pressures for fly reproduction and overwinter survival. This variation in the timing of diapause will enhance populations' survival over the years that differ in the severity of the winter and in the length of the warm period and may also help them respond to long-term changes in environmental conditions.     

Differential fertilization success between two populations of eastern oyster, Crassostrea virginica

Identification of mechanisms promoting prezygotic reproductive isolation and their prevalence are key goals in evolutionary biology because of their potential role in speciation. In marine broadcast-spawning species, molecular interactions between gamete surface proteins are more important than mating behavior for determining reproductive compatibility. Evidence for differential fertilization capacity has been reported from experiments utilizing competing sperm from two males sampled within populations and between species, but to our knowledge conspecific populations that might have diverged in allopatry have never been tested on the basis of sperm competition. In the present study, the gametic compatibility and embryo survivorship from matings between two allopatric populations of Crassostrea virginica, the eastern oyster, on either side of a genetic step cline were investigated. Fertilization success, embryo survival, and paternity data all indicated an absence of strong reproductive barriers between the two oyster populations, implicating other mechanisms for maintenance of the cline step. Sperm from northern male oysters showed a tendency to produce more larvae than expected when competing with sperm from southern male oysters. Although the northern male advantage was not strong, the trend implies that long-distance dispersal across the step cline might more successfully result in north-to-south gene flow than the reverse, providing a mechanistic hypothesis explaining the asymmetric cline shape.     

Secondary contact zones of closely‐related Erebia butterflies overlap with narrow phenotypic and parasitic clines

Zones of secondary contact between closely related taxa are a common legacy of the Quaternary ice ages. Despite their abundance, the factors that keep species apart and prevent hybridization are often unknown. Here, we study a very narrow contact zone between three closely related butterfly species of the Erebia tyndarus species complex. Using genomic data, we first determined whether gene flow occurs and then assessed whether it might be hampered by differences in chromosome number between some species. We found interspecific gene flow between sibling species that differ in karyotype by one chromosome. Conversely, only F1 hybrids occurred between two species that have the same karyotype, forming a steep genomic cline. In a second step, we fitted clines to phenotypic, ecological and parasitic data to identify the factors associated with the genetic cline. We found clines for phenotypic data and the prevalence of the endosymbiont parasite Wolbachia to overlap with the genetic cline, suggesting that they might be drivers for separating the two species. Overall, our results highlight that some gene flow is possible between closely related species despite different chromosome numbers, but that other barriers restrict such gene flow.     

Strong genetic clines and geographical variation in gene flow in the rocky intertidal barnacle Balanus glandula

A long-standing issue in marine biology is identifying spatial scales at which populations of sessile adults are connected by planktonic offspring. We examined the genetic continuity of the acorn barnacle Balanus glandula, an abundant member of rocky intertidal communities of the northeastern Pacific Ocean, and compared these genetic patterns to the nearshore oceanography described by trajectories of surface drifters. Consistent with its broad dispersal potential, barnacle populations are genetically similar at both mitochondrial (cytochrome oxidase I) and nuclear (elongation factor 1-alpha) loci across broad swaths of the species' range. In central California, however, there is a striking genetic cline across 475 km of coastline between northern and southern populations. These patterns indicate that gene flow within central California is far more restricted spatially than among other populations. Possible reasons for the steep cline include the slow secondary introgression of historically separated populations, a balance between diversifying selection and dispersal, or some mix of both. Geographic trajectories of oceanic drifters closely parallel geographical patterns of gene flow. Drifters placed to the north (Oregon; approximately 44 degrees N) and south (Santa Barbara, California; approximately 34 degrees N) of the cline disperse hundreds of kilometers within 40 days, yet over the long-term their trajectories never overlapped. The lack of communication between waters originating in Oregon and southern California probably helps to maintain strong genetic differentiation between these regions. More broadly, the geographical variation in gene flow implies that focusing on species-level averages of gene flow can mask biologically important variance within species which reflects local environmental conditions and historical events.     

Cline coupling and uncoupling in a stickleback hybrid zone

Strong ecological selection on a genetic locus can maintain allele frequency differences between populations in different environments, even in the face of hybridization. When alleles at divergent loci come into tight linkage disequilibrium, selection acts on them as a unit and can significantly reduce gene flow. For populations interbreeding across a hybrid zone, linkage disequilibria between loci can force clines to share the same slopes and centers. However, strong ecological selection on a locus can also pull its cline away from the others, reducing linkage disequilibrium and weakening the barrier to gene flow. We looked for this “cline uncoupling” effect in a hybrid zone between stream resident and anadromous sticklebacks at two genes known to be under divergent natural selection (Eda and ATP1a1) and five morphological traits that repeatedly evolve in freshwater stickleback. These clines were all steep and located together at the top of the estuary, such that we found no evidence for cline uncoupling. However, we did not observe the stepped shape normally associated with steep concordant clines. It thus remains possible that these clines cluster together because their individual selection regimes are identical, but this would be very surprising given their diverse roles in osmoregulation, body armor, and swimming performance.     

Distinguishing between primary and secondary intergradation among morphologically differentiated populations of Anolis marmoratus

Distinguishing between primary and secondary intergradation among differentiated populations, and the relative importance of drift and selection, are persistent problems in evolutionary biology. An historical perspective on population interactions can provide insight into the nature of contacts, and thus help resolve these questions. Continuously distributed populations of Anolis marmoratus from the island of Basse Terre in the Guadeloupean archipelago of the Lesser Antilles show a striking degree of geographic variation in morphology. Initial surveys of mtDNA variation from throughout the Guadeloupean Archipelago revealed one case where levels of sequence difference and phylogenetic relationships of alleles from morphologically differentiated populations from the east coast of Basse Terre were consistent with primary intergradation. In this paper, I examine the genetic population structure of a series of populations spanning this north-south cline in morphological variation to test the hypothesis of primary intergradation. Sequences of the mitochondrial cytochrome-b gene from 50 individuals representing five populations spanning the cline were obtained and fourteen unique haplotypes (differing by 2% or less) were detected. Patterns of nucleotide substitution among haplotypes do not deviate from neutral expectation indicating no effect of selection at the level of mtDNA sequences. Estimates of population structure and gene flow were made using both summary statistics for nucleotide diversity (N_at) and cladistic methods. The results are sensitive to the choice of gene flow model, and this is discussed in detail. Mitochondrial variation in the northern populations may not be at equilibrium, and the phylogeny of alleles is consistent with a recent increase in effective population size. Estimates of nucleotide diversity, gene flow, and the phylogenetic relationships of haplotypes indicate that the southern-most population (representing the extreme of morphological variation along this cline) has been relatively isolated from populations to the north and has experienced a reduced effective population size. The apparent clinal variation between the southern population and the others may therefore reflect secondary contact and introgression rather than primary intergradation.     

Natural selection along an environmental gradient: a classic cline in mouse pigmentation

We revisited a classic study of morphological variation in the oldfield mouse ( Peromyscus polionotus ) to estimate the strength of selection acting on pigmentation patterns and to identify the underlying genes. We measured 215 specimens collected by Francis Sumner in the 1920s from eight populations across a 155-km, environmentally variable transect from the white sands of Florida's Gulf coast to the dark, loamy soil of southeastern Alabama. Like Sumner, we found significant variation among populations: mice inhabiting coastal sand dunes had larger feet, longer tails, and lighter pigmentation than inland populations. Most striking, all seven pigmentation traits examined showed a sharp decrease in reflectance about 55 km from the coast, with most of the phenotypic change occurring over less than 10 km. The largest change in soil reflectance occurred just south of this break in pigmentation. Geographic analysis of microsatellite markers shows little interpopulation differentiation, so the abrupt change in pigmentation is not associated with recent secondary contact or reduced gene flow between adjacent populations. Using these genetic data, we estimated that the strength of selection needed to maintain the observed distribution of pigment traits ranged from 0.0004 to 21%, depending on the trait and model used. We also examined changes in allele frequency of SNPs in two pigmentation genes, Mc1r and Agouti , and show that mutations in the cis -regulatory region of Agouti may contribute to this cline in pigmentation. The concordance between environmental variation and pigmentation in the face of high levels of interpopulation gene flow strongly implies that natural selection is maintaining a steep cline in pigmentation and the genes underlying it.     

The deviation from linkage equilibrium with multiple loci varying in a stepping-stone cline

The balance between the creation of associations between alleles at different loci by immigration and the convergence to linkage equilibrium due to the recombination process is studied in a theoretical model. The geographical structure of the model is a stepping-stone chain of populations linking two genetically constant source populations. The model assumes an arbitrary number of autosomal loci and considers genetic variation (two alleles at each locus) that is not subject to natural selection. The gene frequencies at each locus will then show a linear cline through the stepping-stone chain of populations. The deviation from linkage equilibrium through the stepping-stone cline is characterized by an equation for linear measures that provide the linkage disequilibrium measures for a given set of loci in terms of the gene frequencies and the linkage disequilibria in the source populations and in terms of the linkage disequilibrium measures through the cline for lower numbers of loci. Numerical examples of this iterative solution are given, and it is shown that the build-up of the higher order Bennett-disequilibria through the cline is considerably more pronounced than the build-up of two-locus disequilibria.     

Molecular Analysis of an Allozyme Cline: Alcohol Dehydrogenase in Drosophila Melanogaster on the East Coast of North America

Clines may either be selectively maintained or be the by-product of nonadaptive processes related to population structure and history. Drosophila melanogaster populations on the east coast of North America show a latitudinal cline in the frequencies of two common electrophoretically distinguishable alleles at the alcohol dehydrogenase locus (Adh), designated Adh-S and Adh-F. This cline may either be adaptive or an artifact of a possible recent dual founding of North American D. melanogaster populations in which frequencies of Adh alleles differed between founder populations. By means of a high resolution restriction-mapping technique, we studied the distribution of 113 haplotypes derived from 44 polymorphic DNA markers within the Adh region in 1533 individuals from 25 populations throughout the cline. We found significant clinal differentiation at the polymorphisms determining the mobility-difference causing amino acid replacement between Adh-F and Adh-S alleles. Hitchhiking was limited, despite extensive linkage disequilibrium, and other sites did not vary clinally. Such a pattern of differentiation implies that selection is responsible for the cline. To investigate whether selection acts only on the Adh-F/S site, we performed a ``selective equivalence'''' test under the assumption that all variability within the specified allelic class is selectively neutral. This revealed selective equivalence among Adh-S-bearing haplotypes, whose frequencies showed no differentiation throughout the cline, implying high levels of frequency-homogenizing gene flow. Geographical heterogeneity among Adh-F-bearing haplotypes implied the action of selection on one or more additional variants in linkage disequilibrium with Adh-F. In a further study of a subset of the data (n = 1076 from 18 populations), we found a combined insertion/deletion polymorphism, designated &1, located in the 5'' adult intron and in linkage disequilibrium with Adh-F, to show more marked clinal variation than Adh-F/S. Although the unequivocal identification of the precise target(s) of selection requires further study, we suggest that clinal selection may be acting epistatically on the Adh-F/S and &1 polymorphisms.     

中国板栗居群间等位酶基因频率的空间分布

中国板栗21个自然群间等位酶遗传变异的空间自相关分析及F－统计分析结果表明：其多数等位基因频率在居群间呈随机分布模式,缺乏一定的空间结构;而部分等位基因表现为渐变或双向渐变的非随机分布模式,又具有特定空间结构。中国板栗遗传变异空间结构模式的形成可能是长距离基因流、自然气候、人类活动、地理距离隔离等诸因素综合作用的结果。文中还就居群等位基因分布格局的成因进行了讨论：在第四纪冰川后,中国板栗以长江流域中下游的孑遗中心为起点,等位基因分别沿着向北和向南的不同方向迁移形成现在的居群结构;季风气候和人类活动干扰是削弱居群分化的主要因素,而基于环境梯度的选择,是形成由北向南渐变分布的原因。     

Northern gene flow into southeastern East Asians inferred from genome-wide array genotyping

The population history of Southeast (SE) China remains poorly understood due to the sparse sampling of present-day populations and limited modeling with ancient genomic data. We report genome-wide genotyping data from 207 present-day Han Chinese and Hmong-Mien (HM)-speaking She people from Fujian and Taiwan Island, SE China. We coanalyzed 66 Early Neolithic to Iron Age ancient Fujian and Taiwan Island individuals obtained from previously published works to explore the genetic continuity and admixture based on patterns of genetic variations of the high-resolution time transect. We found the genetic differentiation between northern and southern East Asians was defined by a north–south East Asian genetic cline and our studied southern East Asians were clustered in the southern end of this cline. The southeastern coastal modern East Asians are genetically similar to other southern indigenous groups as well as geographically close to Neolithic-to-Iron Age populations, but they also shared excess alleles with post-Neolithic Yellow River ancients, which suggested a southward gene flow on the modern southern coastal gene pool. In addition, we identified one new HM genetic cline in East Asia with the coastal Fujian HM-speaking She localizing at the intersection between HM and Han clines. She people show stronger genetic affinity with southern East Asian indigenous populations, with the main ancestry deriving from groups related to southeastern ancient indigenous rice farmers. The southeastern Han Chinese could be modeled with the primary ancestry deriving from the group related to the Yellow River Basin millet farmers and the remaining from groups related to rice farmers, which was consistent with the northern China origin of modern southeastern Han Chinese and in line with the historically and archaeologically attested southward migrations of Han people and their ancestors. Our estimated north–south admixture time ranges based on the decay of the linkage disequilibrium spanned from the Bronze Age to historic periods, suggesting the recent large-scale population migrations and subsequent admixture participated in the formation of modern Han in SE Asia.     

Admixture of genetic lineages of different glacial origin: a case study of <Emphasis Type="Italic">Abies alba</Emphasis> Mill. in the Carpathians

Genetic exchange between divergent lineages of silver fir (Abies alba Mill.) was studied in the Ukrainian Carpathians where two expanding populations originating from different glacial refugia meet. The study included 33 silver fir populations from Ukraine, Romania and Slovakia which were investigated using the maternally inherited mitochondrial nad5-4 marker and biparentally inherited nuclear microsatellites. The boundary between mitochondrial lineages is very sharp; only six populations containing a mixture of different haplotypes were found. Bayesian analysis of population structure based on seven nSSR loci revealed the existence of two clusters which coincided fairly well with mitochondrial lineages. Both haplotype frequencies and proportions of clusters identified by the Bayesian analysis exhibited a clinal transition over the contact zone, with cline widths of 17.6 km for mitochondrial haplotype frequencies (reflecting gene flow by seeds) and 119.6 km for Bayesian clusters based on nSSR (reflecting gene flow by pollen). Allelic richness and gene diversity differ significantly between mitochondrial lineages, the Balkan group being more variable, but an increase in gene diversity towards the boundary between lineages was observed only within the Balkan lineage. The observed patterns are suggested to reflect the postglacial colonization and historical gene flow. They demonstrated that in silver fir as a wind-dispersed species, the colonization front is quite continuous, and the survival of migrant seeds in already established populations is low. The relative contribution of pollen-mediated gene flow to genetic exchange between divergent lineages associated with glacial refugia is much higher than in the case of seeds, but pollen dispersal distance is lower than suggested by earlier studies.     

Genetic structure of natural populations of Castanea sativa in Turkey: evidence of a hybrid zone

This study points out the evidence of a hybrid zone between two groups of genetically differentiated populations of chestnut (Castanea sativa Mill.) in Turkey. Genetic structure, gene flow and introgression levels, based on 16 allozyme loci, were investigated on 34 population samples spanning the entire C. sativa distribution area in this country. The occurrence of the hybrid zone, located in the Bithynian region, was inferred in a group of populations showing the following genetic characteristics: (i) enhanced genetic variability and intermediate allelic frequencies between those of the Western and Eastern groups of populations; (ii) sharp and concordant changes in allele frequencies; (iii) decreased gene flow with the Western and Eastern populations. Starting from the cline width estimated to be 324 km, strength of selection was evaluated from the gene flow distance, as indicated from the degree of genetic structuring outside the hybrid zone. Evolutionary processes shaping the observed genetic differentiation and introgression are discussed on the basis of palynological data, palaeoclimatic events and evidence of hybridization found in other plant and animal species in the same region.     

Clinal variation in morphology and assortative mating in the soldier beetle, Chauliognathus pennsylvanicus (Goleoptera: Cantharidae)

In north Georgia populations of the soldier beetle, Chauliognathus pennsylvanicus , the length of the elytral spot varies clinally. At the southern end of a 200 km cline the distribution of spot length is unimodal with longer spot lengths predominating while at the northern end of the cline the distribution is bimodal but with shorter spot lengths being more frequent. North of the cline only short elytral spot lengths are observed, while the converse is true south of the cline. The strength of assortative mating on the basis of elytral spot length increases from south to north along the cline, resulting in complete pre-mating isolation between short and long spot length morphs at the north end of the cline. Laboratory mate choice tests indicate that assortative mating in the field is not the result of differential timing of activity or microhabitat choice but rather that it represents a real behavioural preference. Individuals from monomorphic populations on either side of the cline do not mate assortatively in the laboratory, indicating that reproductive isolation has evolved on the cline.     

中国栗疫病菌群体遗传结构的空间自相关性分析

应用空间自相关分析方法对中国栗疫病菌17个居群RAPD遗传变异的空间结构进行研究,以探讨栗疫病菌居群遗传变异的分布特征及其形成机制。结果表明：中国栗疫病菌居群缺乏空间结构,绝大多数RAPD位点变异为随机分布的空间模式,但部分位点表现出渐变、斑块和双向渐变的非随机分布模式,又显示了一定的空间结构。推测其形成原因可能是长距离的基因流、人类活动、地理隔离以及栗疫病菌本身的繁殖特性综合作用的结果,并依据部分位点呈单向渐变的模式推测西南地区为中国栗疫病菌的起源中心。     

Assessment of gene flow across a hybrid zone in red-tailed chipmunks (Tamias ruficaudus)

The role of hybridization in animal speciation is controversial and recent research has challenged the long-standing criterion of complete reproductive isolation to attain species status. The speciation-with-gene-flow model posits that the genome is semi-permeable and hybridization may be a phase in the process of divergence. Here, we apply these concepts to a previously identified zone of mtDNA introgression between the two strongly morphologically differentiated subspecies of red-tailed chipmunk ( Tamias ruficaudus ) in the US Inland Northwest. Using multilocus genotype data from the southern, older contact zone, we demonstrate that neutral gene flow is unusually low between the subspecies across the Lochsa River. This is geographically congruent with the discontinuity in bacular morphology, indicating that the cline of mitochondrial DNA (mtDNA) haplotypes is displaced. Furthermore, we elucidate the evolutionary forces responsible by testing hypotheses of lineage sorting and hybridization. We determined that introgressive hybridization is the cause of mtDNA/morphology incongruence because there are non-zero levels of migration and gene flow. Although our estimate of the age of the hybrid zone has wide credibility intervals, the hybridization events occurred in the Late Pleistocene and the divergence occurred in the Middle Pleistocene. Finally, we assessed substructure within and adjacent to the hybrid zone and found that the hybrid zone constitutes a set of populations that are genetically differentiated from parental sets of populations; therefore, hybridization in this system is not likely an evolutionary sink, but has generated novel combinations of genotypes.     

Testing limits to adaptation along altitudinal gradients in rainforest Drosophila

Given that evolution can generate rapid and dramatic shifts in the ecological tolerance of a species, what prevents populations adapting to expand into new habitat at the edge of their distributions? Recent population genetic models have focused on the relative costs and benefits of migration between populations. On the one hand, migration may limit adaptive divergence by preventing local populations from matching their local selective optima. On the other hand, migration may also contribute to the genetic variance necessary to allow populations to track these changing optima. Empirical evidence for these contrasting effects of gene flow in natural situations are lacking, largely because it remains difficult to acquire. Here, we develop a way to explore theoretical models by estimating genetic divergence in traits that confer stress resistance along similar ecological gradients in rainforest Drosophila. This approach allows testing for the coupling of clinal divergence with local density, and the effects of genetic variance and the rate of change of the optimum on the response to selection. In support of a swamping effect of migration on phenotypic divergence, our data show no evidence for a cline in stress-related traits where the altitudinal gradient is steep, but significant clinal divergence where it is shallow. However, where clinal divergence is detected, sites showing trait means closer to the presumed local optimum have more genetic variation than sites with trait means distant from their local optimum. This pattern suggests that gene flow also aids a sustained response to selection.