The Genetic Structure of Natural Populations of Drosophila melanogaster. Xx. Comparison of Genotype-Environment Interaction in Viability between a Northern and a Southern Population

In order to examine the operation of diversifying selection as the maintenance mechanism of excessive additive genetic variance for viability in southern populations in comparison with northern populations of Drosophila melanogaster, two sets of experiments were conducted using second chromosomes extracted from the Ogasawara population (a southern population in Japan) and from the Aomori population (a northern population in Japan). Chromosomal homozygote and heterozygote viabilities were estimated in eight kinds of artificially produced breeding environments. The main findings in the present investigation are as follows: (1) Significant genotype-environment interaction was observed using chromosomes extracted from the Ogasawara population. Indeed, the estimate of the genotype-environment interaction variance for heterozygotes was significantly larger than that of the genotypic variance. On the other hand, when chromosomes sampled from the Aomori population were examined, that interaction variance was significant only for homozygotes and its value was no more than one quarter of that for the chromosomes from the Ogasawara population. (2) The average genetic correlation between any two viabilities of the same lines estimated in the eight kinds of breeding environments for the chromosomes sampled from the Ogasawara population was smaller than that for the chromosomes from the Aomori population both in homozygotes and in heterozygotes, especially in the latter. (3) The stability of heterozygotes over homozygotes against fluctuations of environmental conditions was seen in the chromosomes from the Ogasawara population, but not from the Aomori population. (4) From the excessive genotype-environment interaction variance compared with the genotypic variance in heterozygotes, it was suggested for the chromosomes from the Ogasawara population that the reversal of viability order between homozygotes took place in some environments at the locus level. On the basis of these findings, it is strongly suggested that diversifying selection is operating in a southern population of D. melanogaster on some of the viability polygenes which are probably located outside the structural loci, and the excessive additive genetic variance of viability in southern populations is maintained by this type of selection.     

The Genetic Structure of Natural Populations of DROSOPHILA MELANOGASTER. Xix. Genotype-Environment Interaction in Viability

To investigate whether or not an excess of additive genetic variance for viability detected in southern natural populations of Drosophila melanogaster was created by diversifying selection, genotype-environment interaction was tested as follows. (1) Two karyotype chromosomes were used: 61 second chromosomes with the standard karyotype and 63 second chromosomes carrying In(2L)t. Their homozygote viabilities were larger than 50% of the average viability of random heterozygotes. (2) The effects of two factors (culture media and yeasts) were examined at three levels (the culture media: tomato, corn and banana; and the yeasts: sake, brewer's and baker's). The results of 16 three by three factorial experiments by the Cy method in the same karyotype groups for relative viabilities of homozygotes and heterozygotes elucidated the following findings: (1) there was no significant difference between the two karyotype groups, (2) the variance components of genotype-environment interaction were highly significant, (3) the variance component of heterozygotes was significantly smaller than that of homozygotes. From the experimental findings and previous results, diversifying selection in natural populations acting on viability polygenes to increase the additive genetic variance was suggested. The relation of the present result to protein polymorphism is also discussed.     

The Genetic Variance for Viability and Its Components in a Local Population of DROSOPHILA MELANOGASTER

Two hundred and ninety second chromosomes extracted from a natural population of Drosophila melanogaster were analyzed to estimate the genetic variance of viability and its components by means of a partial diallel cross (Design II of Comstock and Robinson 1952). The additive and dominance variances are estimated to be 0.009 and 0.0012. Using the dominance variance and the inbreeding depression, the effective number of overdominant loci contributing to the variance in viability is estimated to be very small, a dozen or less. Either the actual number of loci is small, or the distribution of viabilities is strongly skewed with a large majority of very weakly selected loci. The additive variance in viability appears to be too large to be accounted for by recurrent harmful mutants or by overdominant loci at equilibrium with various genetic parameters estimated independently. The excess might be due to frequency-dependent selection, to negative correlations between viability and fertility, or possibly to the presence of a mutator. The selection for viability and fertility, or possibly to the presence of a mutator. The selection for viability at the average polymorphic locus must be very slight, of the order of 10(-3) or less.     

Pleiotropic Models of Polygenic Variation, Stabilizing Selection, and Epistasis

We show that in polymorphic populations many polygenic traits pleiotropically related to fitness are expected to be under apparent ``stabilizing selection' independently of the real selection acting on the population. This occurs, for example, if the genetic system is at a stable polymorphic equilibrium determined by selection and the nonadditive contributions of the loci to the trait value either are absent, or are random and independent of those to fitness. Stabilizing selection is also observed if the polygenic system is at an equilibrium determined by a balance between selection and mutation (or migration) when both additive and nonadditive contributions of the loci to the trait value are random and independent of those to fitness. We also compare different viability models that can maintain genetic variability at many loci with respect to their ability to account for the strong stabilizing selection on an additive trait. Let V(m) be the genetic variance supplied by mutation (or migration) each generation, V(g) be the genotypic variance maintained in the population, and n be the number of the loci influencing fitness. We demonstrate that in mutation (migration)-selection balance models the strength of apparent stabilizing selection is order V(m)/V(g). In the overdominant model and in the symmetric viability model the strength of apparent stabilizing selection is approximately 1/(2n) that of total selection on the whole phenotype. We show that a selection system that involves pairwise additive by additive epistasis in maintaining variability can lead to a lower genetic load and genetic variance in fitness (approximately 1/(2n) times) than an equivalent selection system that involves overdominance. We show that, in the epistatic model, the apparent stabilizing selection on an additive trait can be as strong as the total selection on the whole phenotype.     

Statistical Studies on Protein Polymorphism in Natural Populations II. Gene Differentiation between Populations

With the aim of testing the validity of the mutation-drift hypothesis, we examined the pattern of genetic differentiation between populations by using data from Drosophila, fishes, reptiles, and mammals. The observed relationship between genetic identity and correlation of heterozygosities of different populations or species was generally in good agreement with the theoretical expectations from the mutation-drift theory, when the variation in mutation rate among loci was taken into account. In some species of Drosophila, however, the correlation was unduly high. The relationship between the mean and variance of genetic distance was also in good agreement with the theoretical prediction in almost all organisms. We noted that both the distribution of heterozygosity within species and the pattern of genetic differentiation between species can be explained by the same set of genetic parameters in each group of organisms. Alternative hypotheses for explaining these observations are discussed.     

Molecular and Quantitative Genetic Differentiation in Sitobion avenae Populations from Both Sides of the Qinling Mountains

Quantitative trait differences are often assumed to be correlated with molecular variation, but the relationship is not certain, and empirical evidence is still scarce. To address this issue, we sampled six populations of the cereal aphid Sitobion avenae from areas north and south of the Qinling Mountains, and characterized their molecular variation at seven microsatellite loci and quantitative variation at nine life-history traits. Our results demonstrated that southern populations had slightly longer developmental times of nymphs but much higher lifetime fecundity, compared to northern populations. Of the nine tested quantitative characters, eight differed significantly among populations within regions, as well as between northern and southern regions. Genetic differentiation in neutral markers was likely to have been caused by founder events and drift. Increased subdivision for quantitative characters was found in northern populations, but reduced in southern populations. This phenomenon was not found for molecular characters, suggesting the decoupling between molecular and quantitative variation. The pattern of relationships between F_ST and Q_ST indicated divergent selection and suggested that local adaptation play a role in the differentiation of life-history traits in tested S. avenae populations, particularly in those traits closely related to reproduction. The main role of natural selection over genetic drift was also supported by strong structural differences in G-matrices among S. avenae populations. However, cluster analyses did not result in two groups corresponding to northern and southern regions. Genetic differentiation between northern and southern populations in neutral markers was low, indicating considerable gene flow between them. The relationship between molecular and quantitative variation, as well as its implications for differentiation and evolution of S. avenae populations, was discussed.     

Epistasis and the Genetic Divergence of Photoperiodism between Populations of the Pitcher-Plant Mosquito, Wyeomyia Smithii

Parallel crosses between each of two southern (ancestral) and one northern (derived) population of the pitcher-plant mosquito, Wyeomyia smithii, were made to determine the genetic components of population divergence in critical photoperiod, a phenological trait that measures adaptation to seasonality along a climatic gradient. Joint scaling tests were used to analyze means and variances of first- and second-generation hybrids in order to determine whether nonadditive genetic variance, especially epistatic variance, contributed to divergence in critical photoperiod. In both crosses, digenic epistatic effects were highly significant, indicating that genetic divergence cannot have resulted solely from differences in additively acting loci. For one cross that could be tested directly for such effects, higher order epistasis and/or linkage did not contribute to the divergence of critical photoperiod between the constituent populations.     

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.     

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl

Gradients of variation—or clines—have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by‐product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.     

Selection in Finite Populations with Multiple Alleles. III. Genetic Divergence with Centripetal Selection and Mutation

Natural selection for an intermediate level of gene or enzyme activity has been shown to lead to a high frequency of heterotic polymorphisms in populations subject to mutation and random genetic drift. The model assumes a symmetrical spectrum of mutational variation, with the majority of variants having only minor effects on the probability of survival. Each mutational event produces a variant which is novel to the population. Allelic effects are assumed to be additive on the scale of enzyme activity, heterosis arising whenever a heterozygote has a mean level of activity closer to optimal than that of other genotypes in the population.-A new measure of genetic divergence between populations is proposed, which is readily interpreted genetically, and increases approximately linearly with time under centripetal selection, drift and mutation. The parameter is closely related to the rate of accumulation of mutational changes in a cistron over an evolutionary time span.-A survey of published data concerning polymorphic loci in man and Drosophila suggests than an alternative model, based on the superiority of hybrid molecules, is not of general importance. Thirteen loci giving rise to hybrid zones on electrophoresis have a mean heterozygote frequency of 0.22 +/-.06, compared with a value of 0.23 +/-.04 for 16 loci classified as producing no hybrid enzyme.     

Genetic diversity and structure of natural fragmented Chamaecyparis obtusa populations as revealed by microsatellite markers

The genetic diversity and population structure of hinoki (Chamaecyparis obtusa) in Japan were investigated by examining the distribution of alleles at 13 microsatellite loci in 25 natural populations from Iwaki in northern Japan to Yakushima Island in southern Japan. On average, 26.9 alleles per locus were identified across all populations and 4.0% of the genetic variation was retained among populations (G _ST = 0.040). According to linkage disequilibrium analysis, estimates of effective population size and detected evidence of bottleneck events, the genetic diversity of some populations may have declined as a result of fragmentation and/or over-exploitation. The central populations located in the Chubu district appear to have relatively large effective population sizes, while marginal populations, such as the Yakushima, Kobayashi and Iwaki populations, have smaller effective population sizes and are isolated from the other populations. Microsatellite analysis revealed the genetic uniqueness of the Yakushima population. Although genetic differentiation between populations was low, we detected a gradual cline in the genetic structure and found that locus Cos2619 may be non-neutral with respect to natural selection.     

Genetic differentiation of pink salmon oncorhynchus gorbuscha Walbaum in the Asian part of the range

Genetic variation at 19 allozyme (including 11 polymorphic) and 10 microsatellite loci was examined in the population samples of odd- and even-broodline pink salmon from the southern part of Sakhalin Island, Southern Kuril Islands, and the northern coast of the Sea of Okhotsk. The estimates of relative interpopulation component of genetic variation over the allozyme loci, per broodline, were on average 0.43% (GST), while over the microsatellite loci it was 0.26% (the theta(ST) coefficient, F-statistics based on the allele frequency variance), and 0.90% (the rho(ST) coefficient, R-statistics based on the allele size variance). The values of interlinear component constituted 2.34, 0.31, and 1.05% of the total variation, respectively. Using the allozyme loci, statistically significant intralinear heterogeneity was demonstrated among the regions, as well as among the populations of Southern Sakhalin Island. Multivariate scaling based on the allozyme data demonstrated regional clustering of the sample groups, representing certain populations during the spawning run or in different years. Most of the microsatellite loci examined were found to be highly polymorphic (mean heterozygosity > 0.880). The estimates of interlinear, interregional, and interpopulation variation over these loci in terms of theta(ST) values were substantially lower than in terms of rho(ST) values. Regional genetic differentiation, mostly expressed at the allozyme loci among the populations from the northern and southern parts of the Sea of Okhotsk (i.e., between the Sakhalin and Kuril populations), was less expressed at the microsatellite loci. The differentiation between these regions observed can be considered as the evidence in favor of a large-scale isolation by distance characterizing Asian pink salmon. It is suggested that in pink salmon, low genetic differentiation at neutral microsatellite loci can be explained by extremely high heterozygosity,of the loci themselves, as well as by the migration gene exchange among the populations (the estimate of the genetic migration coefficient inferred from the "private" allele data constituted 2.6 to 3.4%), specifically, by the ancient migration exchange, which occurred during postglacial colonization and colonization of the range.     

Neutral additive genetic variance in a metapopulation

For neutral, additive quantitative characters, the amount of additive genetic variance within and among populations is predictable from Wright's FST, the effective population size and the mutational variance. The structure of quantitative genetic variance in a subdivided metapopulation can be predicted from results from coalescent theory, thereby allowing single-locus results to predict quantitative genetic processes. The expected total amount of additive genetic variance in a metapopulation of diploid individual is given by 2Ne sigma m2 (1 + FST), where FST is Wright's among-population fixation index, Ne is the eigenvalue effective size of the metapopulation, and sigma m2 is the mutational variance. The expected additive genetic variance within populations is given by 2Ne sigma e2(1-FST), and the variance among demes is given by 4FSTNe sigma m2. These results are general with respect to the types of population structure involved. Furthermore, the dimensionless measure of the quantitative genetic variance among populations, QST, is shown to be generally equal to FST for the neutral additive model. Thus, for all population structures, a value of QST greater than FST for neutral loci is evidence for spatially divergent evolution by natural selection.     

Geographic patterns in the reproductive ecology of Agave lechuguilla (Agavaceae) in the Chihuahuan desert. II. Genetic variation, differentiation, and inbreeding estimates

Plants with natural variation in their floral traits and reproductive ecology are ideal subjects for analyzing the effects of natural selection and other evolutionary forces on genetic structure of natural populations. Agave lechuguilla shows latitudinal changes in floral morphology, color, and nectar production along its distribution through north-central Mexico. Both the type and abundance of its pollinators also change with latitude. Using starch electrophoresis, we examined the levels and patterns of variation of 13 polymorphic allozyme loci in 11 populations of A. lechuguilla. The overall level of genetic variability was high (H(e) = 0.394), but the levels of genetic variation had no geographic pattern. However, the southern populations exhibited an excess of heterozygotes in relation to expectations for Hardy-Weinberg equilibrium, whereas the northern populations had an excess of homozygotes. Total differentiation among populations was low (θ = 0.083), although gene flow estimates (Nm) varied among groups of populations: southern populations had the lowest levels of genetic differentiation, suggesting high levels of gene flow; northern populations had greater levels of genetic differentiation (θ = 0.115), suggesting low gene flow among them. The patterns and inferences of the genetic structure of the population at the molecular level is consistent with variation in floral traits and pollinator visitation rates across the range of the species.     

Variance of Neutral Genetic Variances within and between Populations for a Quantitative Character

The variances of genetic variances within and between finite populations were systematically studied using a general multiple allele model with mutation in terms of identity by descent measures. We partitioned the genetic variances into components corresponding to genetic variances and covariances within and between loci. We also analyzed the sampling variance. Both transient and equilibrium results were derived exactly and the results can be used in diverse applications. For the genetic variance within populations, sigma 2 omega, the coefficient of variation can be very well approximated as [formula: see text] for a normal distribution of allelic effects, ignoring recurrent mutation in the absence of linkage, where m is the number of loci, N is the effective population size, theta 1(0) is the initial identity by descent measure of two genes within populations and t is the generation number. The first term is due to genic variance, the second due to linkage disequilibrium, and third due to sampling. In the short term, the variation is predominantly due to linkage disequilibrium and sampling; but in the long term it can be largely due to genic variance. At equilibrium with mutation [formula: see text] where u is the mutation rate. The genetic variance between populations is a parameter. Variance arises only among sample estimates due to finite sampling of populations and individuals. The coefficient of variation for sample gentic variance between populations, sigma 2b, can be generally approximated as [formula: see text] when the number of loci is large where S is the number of sampling populations.     

Reproductive Strategies Explain Genetic Diversity in Atlantic Salmon, Salmo salar

Synopsis We investigated the relationship between conservation status and genetic variability in European and North American Atlantic salmon, Salmo salar, populations, many of which have suffered severe bottlenecks. A negative north--south cline exists for the status of population conservation in this species. A literature review of genetic variability and demographic parameters of wild Atlantic salmon populations resulted in no statistical associations between population conservation status and genetic variation at enzyme or VNTR loci. We found however, a negative relationship between male parr maturation rates and geographical latitude for both American and European populations. The increase in effective population size due to participation by mature male parr and the increased proportions of these males in smaller (southern) populations could explain the lack of expected relationship between genetic variation and conservation status.     

Contrasting Levels of Variation in Neutral and Quantitative Genetic Loci on Island Populations of Moor Frogs (<Emphasis Type="Italic">Rana arvalis</Emphasis>)

Reduced levels of genetic variability and a prominent differentiation in both neutral marker genes and phenotypic traits are typical for many island populations as compared to their mainland conspecifics. However, whether genetic diversity in neutral marker genes reflects genetic variability in quantitative traits, and thus, their evolutionary potential, remains typically unclear. Moreover, the phenotypic differentiation on islands could be attributable to phenotypic plasticity, selection or drift; something which seldom has been tested. Using eight polymorphic microsatellite loci and quantitative genetic breeding experiments we conducted a detailed comparison on genetic variability and differentiation between Nordic islands (viz. Gotland, Öland and Læsø) and neighbouring mainland populations of moor frogs (Rana arvalis). As expected, the neutral variation was generally lower in island than in mainland populations. But as opposed to this, higher levels of additive genetic variation (V _A) in body size and tibia length were found on the island of Gotland as compared to the mainland population. When comparing the differentiation seen in neutral marker genes (F _ST) with the differentiation in genes coding quantitative traits (Q _ST) two different evolutionary scenarios were found: while selection might explain a smaller size of moor frogs on Gotland, the differentiation seen in tibia length could be explained by genetic drift. These results highlight the limited utility of microsatellite loci alone in inferring the causes behind an observed phenotypic differentiation, or in predicting the amount of genetic variation in ecologically important quantitative traits.     

Genetic profile of cosmopolitan populations: effects of hidden subdivision

Natural populations of many organisms exhibit excess of rare alleles in comparison with the predictions of the neutral mutation hypothesis. It has been shown before that either a population bottleneck or the presence of slightly deleterious mutations can explain this phenomenon. A third explanation is presented in this work, showing that hidden subdivision within a population can also lead to an excess of rare alleles in the total population when the expectations of the neutral model are based on the allele frequency profile of the entire population data. With two examples (mitochondrial DNA-morph distribution and isozyme allele frequency distributions), it is shown that most cosmopolitan human populations exhibit excess of rare as well as total allele counts, when these are compared with the expectations of the neutral mutation hypothesis. The mitochondrial data demonstrate that such excesses can be detected from genetic variation at a single locus as well, and this is not due to stochastic error of allele frequency distributions. Contrast of the present observations with the allele frequency profiles in agglomerated tribal populations from South and Central America shows that even when the neutral expectations hold for individual subpopulations, if all subpopulations are grouped into a single population, the pooled data exhibit an excess of total number of alleles that is mainly due to the excess of rare alleles. Therefore, a primary cause of the excess number of rare alleles could be the hidden subdivision, and the magnitude of the excess indicates the extent of substructuring. The two components of hidden subdivision are: 1) Number of subpopulations, and 2) the average genetic distance among them. The implications of this observation in estimating mutation rate are discussed indicating the difficulties of comparing mutation rates from different population surveys.