A Population Genomic Assessment of Three Decades of Evolution in a Natural Drosophila Population

Population genetics seeks to illuminate the forces shaping genetic variation, often based on a single snapshot of genomic variation. However, utilizing multiple sampling times to study changes in allele frequencies can help clarify the relative roles of neutral and non-neutral forces on short time scales. This study compares whole-genome sequence variation of recently collected natural population samples of Drosophila melanogaster against a collection made approximately 35 years prior from the same locality—encompassing roughly 500 generations of evolution. The allele frequency changes between these time points would suggest a relatively small local effective population size on the order of 10,000, significantly smaller than the global effective population size of the species. Some loci display stronger allele frequency changes than would be expected anywhere in the genome under neutrality—most notably the tandem paralogs Cyp6a17 and Cyp6a23, which are impacted by structural variation associated with resistance to pyrethroid insecticides. We find a genome-wide excess of outliers for high genetic differentiation between old and new samples, but a larger number of adaptation targets may have affected SNP-level differentiation versus window differentiation. We also find evidence for strengthening latitudinal allele frequency clines: northern-associated alleles have increased in frequency by an average of nearly 2.5% at SNPs previously identified as clinal outliers, but no such pattern is observed at random SNPs. This project underscores the scientific potential of using multiple sampling time points to investigate how evolution operates in natural populations, by quantifying how genetic variation has changed over ecologically relevant timescales.     

Empirical Distributions of F ST from Large-Scale Human Polymorphism Data

Studies of the apportionment of human genetic variation have long established that most human variation is within population groups and that the additional variation between population groups is small but greatest when comparing different continental populations. These studies often used Wright’s F _ST that apportions the standardized variance in allele frequencies within and between population groups. Because local adaptations increase population differentiation, high-F _ST may be found at closely linked loci under selection and used to identify genes undergoing directional or heterotic selection. We re-examined these processes using HapMap data. We analyzed 3 million SNPs on 602 samples from eight worldwide populations and a consensus subset of 1 million SNPs found in all populations. We identified four major features of the data: First, a hierarchically F _ST analysis showed that only a paucity (12%) of the total genetic variation is distributed between continental populations and even a lesser genetic variation (1%) is found between intra-continental populations. Second, the global F _ST distribution closely follows an exponential distribution. Third, although the overall F _ST distribution is similarly shaped (inverse J), F _ST distributions varies markedly by allele frequency when divided into non-overlapping groups by allele frequency range. Because the mean allele frequency is a crude indicator of allele age, these distributions mark the time-dependent change in genetic differentiation. Finally, the change in mean-F _ST of these groups is linear in allele frequency. These results suggest that investigating the extremes of the F _ST distribution for each allele frequency group is more efficient for detecting selection. Consequently, we demonstrate that such extreme SNPs are more clustered along the chromosomes than expected from linkage disequilibrium for each allele frequency group. These genomic regions are therefore likely candidates for natural selection.     

Low genetic diversity but local genetic differentiation in endemic Minasia (Asteraceae) species from Brazil

Espinhaço Range is a Brazilian chain of mountains, extremely rich in endemic species. Minasia has six species, all perennial herbs endemic to this range. Twenty-two populations were sampled from all Minasia species and assayed for allozyme variation. The species showed low genetic variation, compatible with the expectation for endemics. Genetic identities in Serra do Cipó and Diamantina were high, even between populations of different species. On the other hand M. cabralensis, from Serra do Cabral, was more differentiated from the other species. We also observed a marked genetic differentiation within M. cabralensis. Most strikingly, two close populations of this species showed very different allele frequencies. Our findings highlight the importance of local differentiation in campos rupestres. The observed genetic structure indicates that substantial portions of genetic diversity could be lost with the extinction of only one population, which is especially threatening considering the already low genetic diversity.     

Genetic Differentiation between Geographically Distant Populations of DROSOPHILA MELANOGASTER

We have studied allozyme variation at 26 gene loci in nine populations of Drosophila melanogaster originating on five different continents. The distant populations show significant genetic differentiation. However, only half of the loci studied have contributed to this differentiation; the other half show identical patterns in all populations. The genetic differentiation in North American, European and African populations is correlated with the major climatic differences between north and south. These differences arise mainly from seven loci that show gene-frequency patterns suggestive of latitudinal clines in allele frequencies. The clinal variation is such that subtropical populations are more heterozygous than temperate populations. These results are discussed in relation to the selectionist and neutralist hypotheses of genetic variation in natural populations.     

Temporal Stability of Genetic Variability and Differentiation in the Three-Spined Stickleback (Gasterosteus aculeatus)

Temporal variation in allele frequencies, whether caused by deterministic or stochastic forces, can inform us about interesting demographic and evolutionary phenomena occurring in wild populations. In spite of the continued surge of interest in the genetics of three-spined stickleback (Gasterosteus aculeatus) populations, little attention has been paid towards the temporal stability of allele frequency distributions, and whether there are consistent differences in effective size (N_e) of local populations. We investigated temporal stability of genetic variability and differentiation in 15 microsatellite loci within and among eight collection sites of varying habitat type, surveyed twice over a six-year time period. In addition, N_es were estimated with the expectation that they would be lowest in isolated ponds, intermediate in larger lakes and largest in open marine sites. In spite of the marked differences in genetic variability and differentiation among the study sites, the temporal differences in allele frequencies, as well as measures of genetic diversity and differentiation, were negligible. Accordingly, the N_e estimates were temporally stable, but tended to be lower in ponds than in lake or marine habitats. Hence, we conclude that allele frequencies in putatively neutral markers in three-spined sticklebacks seem to be temporally stable – at least over periods of few generations – across a wide range of habitat types differing markedly in levels of genetic variability, effective population size and gene flow.     

Genetic polymorphisms of blood proteins in the troops of Japanese macaques,Macaca fuscata: III. Erythrocyte carbonic anhydrase polymorphism inMacaca fuscata

In order to clarify the genetic relationships between troops of Japanese macaques, the authors have been looking for the genetic variation of blood proteins by electrophoretical technique. In this work the genetic variants of erythrocyte carbonic anhydrase isozymes of 1273 blood samples from 34 troops of this species were examined. The genetic polymorphisms were observed in several troops. The variable troops of CA I were Fukushima, Takasakiyama, Kohchi, Shimane, Shodoshima I, K, and T, Kashima, Kawara, Takasakiyama A, B, and C, and Tomogashima. The variant allele found in these troops is onlyd ₂ allele, which was probably identical with that reported byTashian et al. (1971). The frequencies ofd ₂ allele in the Shodoshima and Kashima troops were very high, and this phenomenon was interpreted as being a result of the founder effect and the random fluctuation of gene frequencies in these troops.     

Validation of SNP Allele Frequencies Determined by Pooled Next-Generation Sequencing in Natural Populations of a Non-Model Plant Species

Sequencing of pooled samples (Pool-Seq) using next-generation sequencing technologies has become increasingly popular, because it represents a rapid and cost-effective method to determine allele frequencies for single nucleotide polymorphisms (SNPs) in population pools. Validation of allele frequencies determined by Pool-Seq has been attempted using an individual genotyping approach, but these studies tend to use samples from existing model organism databases or DNA stores, and do not validate a realistic setup for sampling natural populations. Here we used pyrosequencing to validate allele frequencies determined by Pool-Seq in three natural populations of Arabidopsis halleri (Brassicaceae). The allele frequency estimates of the pooled population samples (consisting of 20 individual plant DNA samples) were determined after mapping Illumina reads to (i) the publicly available, high-quality reference genome of a closely related species (Arabidopsis thaliana) and (ii) our own de novo draft genome assembly of A. halleri. We then pyrosequenced nine selected SNPs using the same individuals from each population, resulting in a total of 540 samples. Our results show a highly significant and accurate relationship between pooled and individually determined allele frequencies, irrespective of the reference genome used. Allele frequencies differed on average by less than 4%. There was no tendency that either the Pool-Seq or the individual-based approach resulted in higher or lower estimates of allele frequencies. Moreover, the rather high coverage in the mapping to the two reference genomes, ranging from 55 to 284x, had no significant effect on the accuracy of the Pool-Seq. A resampling analysis showed that only very low coverage values (below 10-20x) would substantially reduce the precision of the method. We therefore conclude that a pooled re-sequencing approach is well suited for analyses of genetic variation in natural populations.     

Latitudinal countergradient variation in the common frog (Rana temporaria) development rates--evidence for local adaptation

Adaptive genetic differentiation along a climatic gradient as a response to natural selection is not necessarily expressed at phenotypic level if environmental effects on population mean phenotypes oppose the genotypic effects. This form of cryptic evolution--called countergradient variation--has seldom been explicitly demonstrated for terrestrial vertebrates. We investigated the patterns of phenotypic and genotypic differentiation in developmental rates of common frogs (Rana temporaria) along a ca. 1600 km latitudinal gradient across Scandinavia. Developmental rates in the field were not latitudinally ordered, but displayed large variation even among different ponds within a given latitudinal area. In contrast, development rates assessed in the laboratory increased strongly and linearly with increasing latitude, suggesting a genetic capacity for faster development in the northern than the southern larvae. Experiments further revealed that environmental effects (temperature and food) could easily override the genetic effects on developmental rates, providing a possible mechanistic explanation as to why the genetic differentiation was not seen in the samples collected from the wild. Our results suggest that the higher developmental rates of the northern larvae are likely to be related to selection stemming from seasonal time constrains, rather than from selection dictated by low ambient temperatures per se. All in all, the results provide a demonstration of environmental effects concealing substantial latitudinally ordered genetic differentiation understandable in terms of adaptation to clinal variation in time constrains.     

Genetic studies on serum transferrins in Atlantic salmon

Published and unpublished data on genetic variation at the transferrin locus ( TF *) in Atlantic salmon from rivers in eastern North America sampled from 1968 to 1970 were reanalysed and compared with data for samples collected in 1998 from nine of the same rivers. Genetic differentiation among rivers was highly significant as was spatial differentiation among tributary samples within the Miramichi River system, the largest rivers studied. Comparison of allele frequencies in rivers also sampled in 1998 show no overall evidence of significant genetic change after 30 years, spanning 9 generations. The results strongly support the stability of the patterns of spatial genetic differentiation and support the occurrence in Atlantic salmon of reproductive isolation among rivers and among tributaries within large river systems.     

Pantophysin (Pan I) locus divergence between inshore v. offshore and northern v. southern populations of Atlantic cod in the north‐east Atlantic

More than 6000 cod Gadus morhua , sampled in coastal and offshore waters stretching from the Barents Sea down to the North Sea, were analysed for frequencies of alleles at the scnDNA pantophysin locus ( Pan I)[formerly called synaptophysin ( Syp I)]. The significant allele frequency difference between the two major stocks of cod in Norway, north‐east Arctic cod (NEAC) and Norwegian coastal cod (NCC), was upheld in all years of the investigation (1993 to 2001), and applied both to larval cod and post‐juveniles of various ages. On a north‐south axis, the appearance of a latitudinal cline of post‐juvenile (≥1 year) allele frequencies was exposed. The intermediate allele frequencies in coastal areas of northern Norway, seem to a large extent to be caused by intermingling of the two stocks, although the existence of populations of coastal cod with alternative Pan I frequencies could not be ruled out. The role of selection is yet unresolved. Depth of the sampling location seemed to have an effect on the allele frequencies and their temporal stability, while there was no indication of seasonal variation in the frequencies. Breeding structure was the most likely cause for upholding the extreme divergence in Pan I frequencies between NEAC and NCC.     

Quantification and reduction of bias from sampling larvae to infer population and landscape genetic structure

A well-designed sampling scheme is critical for obtaining accurate results from population genetic studies. Larval samples contain only the genetic material of successful breeders, often of a single year, and may be biased towards particular families. To quantify the bias of using larval samples to infer population and landscape genetic structure and explore how this bias may be reduced using sibship analysis, we analysed eight microsatellite loci from 484 tissue samples of larvae and adults of Columbia spotted frogs (Rana luteiventris) and long-toed salamanders (Ambystoma macrodactylum) at nine breeding sites in north Idaho. Differences in allele frequencies between adult and larval samples were not detected after full-siblings were removed from the larval data set for Columbia spotted frogs; for long-toed salamanders, these differences remained at two out of four ponds. Data from Columbia spotted frog larvae indicated higher levels of differentiation among populations (median difference in F_ST = 0.020, P < 0.01), as predicted by population genetic theory, whereas data from larval samples of long-toed salamanders showed some evidence of lower levels of differentiation among populations (median difference in F_ST = 0.012, P = 0.06). For both species, removing all but one individual from each full-sibling family led to parameter estimates that were closer to those calculated from adult samples for both population and landscape genetic measures. Removal of full-siblings is likely to improve estimates of population genetic parameters; however, knowledge of the species’ breeding system is essential for understanding additional sources of bias when inferring population genetic structure from larval samples.     

Evidence for Two Highly Differentiated Herring Groups at Goose Bank in the Barents Sea and the Genetic Relationship to Pacific Herring, Clupea pallasi

Genetic studies on Atlantic herring, Clupea harengus, have generally revealed a low level of genetic variation over large geographic areas. Genetically distinct herring populations in some of the Norwegian fjords are exceptions, and juvenile herring from the large oceanic herring, Norwegian Spring Spawners (NSS), are often found in mixture with local fjord populations as well as widely distributed in the Barents Sea. Research surveys in the eastern Barents Sea (Goose Bank) in 1993, 1994 and 2001 included collection of herring samples for allozyme analyses. As expected the results identified juveniles from NSS stock, but an additional unique group of herring (low vertebrae number), being almost fixed for alternative alleles at several allozyme loci, was detected. In some cases, the two groups of herring were taken in the same trawl catches as documented by highly significant departure from Hardy—Weinberg expectation with large excess of homozygotes providing evidence for population mixing. Large genetic differences (Nei's genetic distance = 1.53; F_ST = 0.754) were detected in pairwise comparisons based on five allozyme loci. The two herring groups were also compared with reference samples of Pacific herring, Clupea pallasi, including one sample from Japan Sea and three Alaskan samples. UPGMA dendrogram based on five allozyme loci revealed a close genetic relationship between the low vertebrae herring in the Barents Sea and the group of samples of Pacific herring. Although significant different in allele frequencies, one of the herring samples clustered together with the reference sample from Bering Sea with genetic distance of 0.008 and F_ST value of 0.032. The close genetic relationship found in this paper, suggest a re-evaluation of the taxonomic status of the Barents Sea herring populations investigated.     

Genetic Analysis of a Population of Tribolium. VI. Polymorphism and Demographic Equilibrium

Demographic and genetic data of continuously growing populations of Tribolium castaneum initiated with identical age structures but with different frequencies of the unsaturated fatty acid-sensitive allele were collected for 68 weeks.—The life history data provided the following insights: Genotypic differences for total number of offspring were due primarily to larval viability. The total lifetime offspring production of the genotypes predicted a stable polymorphic genetic equilibrium. The genotypic reproductive functions forecast that a stable age structure was not a prerequisite for genetic equilibrium.—Those cultures initially segregating for the unsaturated fatty acid-sensitive allele converged to an equilibrium allele frequency of 0.25 and a genotypic array composed of equal numbers of +/+ and +/cos individuals.—The numbers of larvae, pupae and adults during the first six weeks were quadratic functions of the initial frequency of the sensitive allele. Qualitative age structure changes that followed were similar in all cultures and demographic equilibrium was realized at week 50. The overall demographic pattern during the 68-week study was interpreted in terms of the interactions among the numbers of small larvae, large larvae plus pupae, and adults.     

Confidence Intervals for Population Allele Frequencies: The General Case of Sampling from a Finite Diploid Population of Any Size

The estimation of population allele frequencies using sample data forms a central component of studies in population genetics. These estimates can be used to test hypotheses on the evolutionary processes governing changes in genetic variation among populations. However, existing studies frequently do not account for sampling uncertainty in these estimates, thus compromising their utility. Incorporation of this uncertainty has been hindered by the lack of a method for constructing confidence intervals containing the population allele frequencies, for the general case of sampling from a finite diploid population of any size. In this study, we address this important knowledge gap by presenting a rigorous mathematical method to construct such confidence intervals. For a range of scenarios, the method is used to demonstrate that for a particular allele, in order to obtain accurate estimates within 0.05 of the population allele frequency with high probability (%), a sample size of is often required. This analysis is augmented by an application of the method to empirical sample allele frequency data for two populations of the checkerspot butterfly (Melitaea cinxia L.), occupying meadows in Finland. For each population, the method is used to derive % confidence intervals for the population frequencies of three alleles. These intervals are then used to construct two joint % confidence regions, one for the set of three frequencies for each population. These regions are then used to derive a % confidence interval for Jost''s D, a measure of genetic differentiation between the two populations. Overall, the results demonstrate the practical utility of the method with respect to informing sampling design and accounting for sampling uncertainty in studies of population genetics, important for scientific hypothesis-testing and also for risk-based natural resource management.     

Lack of congruence between morphometric evolution and genetic differentiation suggests a recent dispersal and local habitat adaptation of the Madeiran lizard Lacerta dugesii

Genetic differentiation among nine populations of the endemic lizard Lacerta dugesii Milne-Edwards 1829 (Lacertidae) from four groups of islands constituting the Archipelago of Madeira, was investigated by protein electrophoresis at 23 enzyme loci. Among twenty polymorphic loci, the total genetic diversity was due primarily to intra-population variation. The allele and genotypic frequencies among populations showed some heterogeneity, allowing the species to present a structuring pattern compatible with their geographical clustering. Some evidence suggests that selection acting on some loci in different ecological conditions may be responsible for the clustering of the populations studied. There was no apparent isolation effect expected under an "island" model of population divergence, and no correlation was found between genetic and geographic distances among populations. Morphological variation of the proposed three L. dugesii subspecies is not congruent with the allozyme analysis. This most probably suggests a rapid colonization of the islands followed by a strong effect of selection operating over the morphological characters used to define the subspecies.     

Geographic variation and temporal population differentiation of chum salmon <Emphasis Type="Italic">Oncorhynchus keta</Emphasis> from some regions of the Russian Far East

Allozyme variation of chum salmon Oncorhynchus keta was examined in the populations from nine rivers of the Russian Far East. A total of 12 loci were tested, and eight of these were shown to be polymorphic. The greatest contribution to the samples differentiation was made by the EST-D1* and sIDHP-2* loci. Most of the allele diversity was distributed among the individuals within the samples. The among-sample differences within races were statistically significant and accounted for slightly more than 7% of total variation. The differences between the races constituted 0.1% and were statistically nonsignificant. Genetic similarity of seasonal races within single locality allowed suggestion on their secondary formation after settling of founder population into a certain region.     

Population genetic structure of mussels from the Baltic Sea

In a macrogeographic survey, the population genetic structure of mussels from various regions of the Baltic Sea, a large semi-enclosed brackish-water basin, was examined with reference toMytilus edulis andM. galloprovincialis samples from the North Sea, Irish coast and southern Portugal. Electrophoretically detectable variation was analysed at 6 polymorphic enzyme loci (Ap, Est-D, Lap-2, Odh, Pgi andPgm). Evidence was provided of a remarkably large amount of biochemical genetic differentiation among ecologically and morphologically divergent mussel populations in the Baltic. Patterns of allele frequencies in low-salinity populations from the area of the Baltic Proper were demonstrated to be widely homogeneous but contrast strongly with those of the western Baltic, the latter resembling populations from marine habitats of the North Sea. Associated with a pronounced salinity gradient, the spatial heterogeneity in gene-pool structure is indicated by steep clines of allele frequency changes in the area of the eastern Danish isles. The adaptive significance of the observed allozymic variation is suggested. From genetic distance estimates, the subdivision of population structure is discussed in relation to the significant amount of differentiation detected withinMytilus populations to date and to the evolutionary time required for the divergence of Baltic mussel populations. The allozymic data provide evidence for the genetic distinctiveness of mussels from the low-salinity areas of the Baltic. Their position at the specific or subspecific level of classification requires further consideration.     

Population variation of elytral ridge occurrence in ladybirds <Emphasis Type="Italic">Harmonia axyridis</Emphasis> Pallas

Intra- and interpopulation variation of occurrence of elytral ridge (neutral morphological character with the known genetic determination) in some populations of Harmonia axyridis from the Russian part of the species area was studied for the first time. Comparative analysis of the frequencies of the recessive allele r of a corresponding gene and the frequencies of recessive homozygotes in 32 samples from 16 localities has shown the lack of temporal and microgeographic variations of this trait. The character of its geographic variation con-firms the hypothesis advanced earlier on the basis of analysis of elytral pattern variation about the existence of two subspecies of H. axyridis and suggests the existence of a zone of secondary hybridization between them in the central part of the area.     

Pattern of genetic differentiation in the Maculinea alcon species group (Lepidoptera, Lycaenidae) in Central Europe

The taxonomic status of Alcon Blues living in Central and Western Europe (conventionally: Maculinea alcon and Maculinea rebeli ) is quite confused. Some authors distinguish them as separate species of the M. alcon species group, while others consider them as subspecies or simply ecological forms. Our aim was to study the pattern of genetic differentiation among several populations of these taxa with reference to their taxonomic position. Imagoes were collected from 27 localities in Central Europe between 1999 and 2003. Samples originated from four regions: northern Hungary, central Hungary, Romania: Transylvania and Slovenia. Enzyme polymorphism was analysed using polyacrylamide gel electrophoresis. In all samples 16 enzyme loci were studied. In the analysis of the data, F -statistics was computed and the total genetic variation was partitioned into within and between population components. Nei's genetic distances were calculated and UPGMA dendrogram was constructed on the basis of the distance matrix. Hierarchical F -statistics and amova was computed to study the pattern of genetic differentiation among the samples. Principal component analysis (PCA) was also carried out using the allele frequencies of the samples. The results of all analyses indicated a high level of differentiation among the samples. The differences among the samples collected in different years from the same population accounted for a considerable amount of the between sample variation. The distribution of the between sample variation did not exhibit a species pattern. The results of amova and PCA indicated that the geographic pattern was slightly more expressed in the between sample variation compared with the species pattern of it. These results seem to contradict with the conventional taxonomical subdivision.     

Population genetic differentiation of white-spotted char Salvelinus leucomaenis (Pallas) in Russian Far East

The population genetic structure of white-spotted char Salvelinus leucomaenis (family Salmonidae) was determined based on variations in ten microsatellite DNA loci in samples from different parts of the species range in the Russian Far East. In a number of samples, variations in 21 allozyme loci were examined, of which five loci were found to be polymorphic. The overall diversity level at the examined markers was comparable to that observed in a closely related char species, Dolly Varden. The estimates of interpopulation genetic differentiation were highly statistically significant in most of the pairwise comparisons among the samples. The overall evaluation of the spatial genetic differentiation in white-spotted char constituted F _ST = 0.203 and R _ST = 0.202. Cluster analysis and multidimensional scaling based on microsatellite allele frequencies indicated the possible subdivision of the examined samples into two main groups, i.e., northern (represented by the regions of the north of Khabarovsk krai, Kamchatka, Yama Bay) and southern (with regions including Sakhalin, Primorye, and the Kuril Islands). The allozome data demonstrated a similar pattern of differentiation. The level of intra- and interpopulation genetic diversity in the southern group was higher than in the northern group. The isolation-by-distance test did not identify a significant correlation between genetic and geographic distances among the samples. The data obtained enabled the suggestion that the genetic structure of the populations of white-spotted char was shaped by the influence of historical geological climatic rearrangements of its range and the genetic drift because of relatively low population number and limited in the extent migration activity of its anadromous form.