Evolutionary inference from QST

Q(ST) is a commonly used metric of the degree of genetic differentiation among populations displayed by quantitative traits. Typically, Q(ST) is compared to F(ST) measured on putatively neutral loci; if Q(ST)=F(ST), this is taken as evidence of spatially heterogeneous and diversifying selection. This paper reviews the uses, assumptions and statistics of Q(ST) and F(ST) comparisons. Unfortunately, Q(ST)/F(ST) comparisons are statistically challenging. For a single trait, Q(ST) must be compared not to the mean F(ST) but to the distribution of F(ST) values. The sources of biases and sampling error for Q(ST) are reviewed, and a new method for comparing Q(ST) and F(ST) is suggested. Simulation results suggest that the distribution of neutral F(ST) and Q(ST) values are little affected by various deviations from the island model. Consequently, the distributions of Q(ST) and F(ST) are well approximated by the Lewontin-Krakauer prediction, even with realistic deviations from the island-model assumptions.     

Population differentiation and nuclear gene flow in the Dominican anole (Anolis oculatus)

Allele frequency data from nuclear microsatellite loci were used to investigate patterns of nuclear gene flow and population structure in the morphologically variable Dominican anole (Anolis oculatus). All six loci used proved to be highly polymorphic, with an average of 18.8 alleles per locus. Test for Hardy-Weinberg equilibrium revealed small numbers of heterozygote deficiencies at single loci in single populations and consistent patterns of increasingly significant heterozygote deficiency in global tests across populations and loci. No significant relationship between FST and patristic distances estimated from mitochondrial DNA sequences was detected and estimates of FIS were significantly higher in females than in males, indicating that gene flow may be sex-biased and mediated mainly by male migration. A highly significant correlation between linearized FST and loge (geographical distance) indicates that geographical proximity is a significant factor in the genetic structure of A. oculatus populations. However, levels of gene flow between morphologically differentiated parapatric populations are frequently seen to be relatively high. This supports the hypothesis of natural selection being the driving force behind the development and maintenance of morphological variation and shows that adaptive differentiation may be maintained despite the homogenizing influence of gene flow. Generally, the morphologically variable populations of A. oculatus are seen to be poor candidates for in situ speciation, but an exceptional case on the west coast of Dominica indicates that isolation resulting from vicariant events may lead to rapid differentiation at both mitochondrial and nuclear loci. This provides a possible mechanism for anole speciation on other Caribbean islands.     

Macrogeographic population structure of the tsetse fly, Glossina pallidipes (Diptera: Glossinidae)

Tsetse flies are confined to sub-Saharan Africa where they occupy discontinuous habitats. In anticipation of area-wide control programmes, estimates of gene flow among tsetse populations are necessary. Genetic diversities were partitioned at eight microsatellite loci and five mitochondrial loci in 21 Glossina pallidipes Austin populations. At microsatellite loci, Nei's unbiased gene diversity averaged over loci was 0.659 and the total number of alleles was 214, only four of which were shared among all populations. The mean number of alleles per locus was 26.8. Random mating was observed within but not among populations (fixation index FST=0.18) and 81% of the genetic variance was within populations. Thirty-nine mitochondrial variants were detected. Mitochondrial diversities in populations varied from 0 to 0.85 and averaged 0.42, and FST=0.51. High levels of genetic differentiation were characteristic, extending even to subpopulations separated by tens and hundreds of kilometres, and indicating low rates of gene flow.     

Comparing relative rates of pollen and seed gene flow in the island model using nuclear and organelle measures of population structure

We describe a method for comparing nuclear and organelle population differentiation (F(ST)) in seed plants to test the hypothesis that pollen and seed gene flow rates are equal. Wright's infinite island model is used, with arbitrary levels of self-fertilization and biparental organelle inheritance. The comparison can also be applied to gene flow in animals. Since effective population sizes are smaller for organelle genomes than for nuclear genomes and organelles are often uniparentally inherited, organelle F(ST) is expected to be higher at equilibrium than nuclear F(ST) even if pollen and seed gene flow rates are equal. To reject the null hypothesis of equal seed and pollen gene flow rates, nuclear and organelle F(ST)'s must differ significantly from their expected values under this hypothesis. Finite island model simulations indicate that infinite island model expectations are not greatly biased by finite numbers of populations (>/=100 subpopulations). The power to distinguish dissimilar rates of pollen and seed gene flow depends on confidence intervals for fixation index estimates, which shrink as more subpopulations and loci are sampled. Using data from the tropical tree Corythophora alta, we rejected the null hypothesis that seed and pollen gene flow rates are equal but cannot reject the alternative hypothesis that pollen gene flow is 200 times greater than seed gene flow.     

Population genetic analysis of white shrimp,Litopenaeus setiferus,using microsatellite genetic markers

The white shrimp (Litopenaeus setiferus) is a commercially and recreationally valuable species, yet little is known of its population structure or genetic diversity. White shrimp are distributed along the Atlantic coast of the United States and from the west coast of Florida to the Bay of Campeche, Mexico. In this study, shrimp were collected from North Carolina, South Carolina (four separate collections were taken from 1995 to 1999), Georgia, the Atlantic and Gulf coasts of Florida, Louisiana, Texas and Mexico. DNA was isolated from these individuals, and genetic variation was assessed at six microsatellite loci. These loci were, for the most part, highly polymorphic with an average expected heterozygosity of 0.68. Deviations from Hardy-Weinberg proportions were observed over all samples, but experimental results suggested the presence of null alleles, which confounded a biological interpretation of this result. Pairwise tests of the similarity of allele frequency distributions and distance measure analyses showed broad-scale genetic homogeneity superimposed over occasional indications of random geographical and temporal differentiation. FST and RST estimates over all loci and samples were 0.002 or less and indicated little population structure. Weak but significant genetic differentiation was evident only between pooled western Atlantic and pooled Gulf of Mexico samples. Within the Gulf of Mexico or within the western Atlantic, the large-scale genetic homogeneity observed may be a consequence of genetic mixing resulting from pelagic larvae and adult migrations, while the random local genetic differentiation may be a result of genetic sampling or experimental sampling error. The weak differentiation between shrimp from the Gulf of Mexico and the western Atlantic can be explained by a relatively recent separation of these two populations and/or small amounts of ongoing gene flow.     

Isolation by distance in the Atlantic cod, Gadus morhua, at large and small geographic scales

Genetic isolation by distance (IBD) has rarely been described in marine species with high potential for dispersal at both the larval and adult life-history stages. Here, we report significant relationships between inferred levels of gene flow and geographic distance in the Atlantic cod, Gadus morhua, at 10 nuclear restriction-fragment-length-polymorphism (RFLP) loci at small regional scales in the western north Atlantic region (< 1,600 km) that mirror those previously detected over its entire geographic range (up to 7,300 km). Highly significant allele frequency differences were observed among eight northwestern Atlantic populations, although the mean FST for all 10 loci was only 0.014. Despite this weak population structuring, the distance separating populations explained between 54% and 62% of the variation in gene flow depending on whether nine or 10 loci were used to estimate Nm. Across the species' entire geographic range, highly significant differences were observed among six regional populations at nine of the 10 loci (mean FST = 0.068) and seven loci exhibited significant negative relationships between gene flow and distance. At this large geographic scale, natural selection acting in the vicinity of one RFLP locus (GM798) had a significant effect on the correlation between gene flow and distance, and eliminating it from the analysis caused the coefficient of determination to increase from 17% to 62%. The role of vicariance was assessed by sequentially removing populations from the analysis and was found to play a minor role in contributing to the relationship between gene flow and distance at either geographic scale. The correlation between gene flow and distance detected in G. morhua at small and large spatial scales suggests that dispersal distances and effective population sizes are much smaller than predicted for the species and that the recent age of populations, rather than extensive gene flow, may be responsible for its weak population structure. Our results suggest that interpreting limited genetic differences among populations as reflecting high levels of ongoing gene flow should be made with caution.     

Estimation of Levels of Gene Flow from DNA Sequence Data

We compare the utility of two methods for estimating the average levels of gene flow from DNA sequence data. One method is based on estimating FST from frequencies at polymorphic sites, treating each site as a separate locus. The other method is based on computing the minimum number of migration events consistent with the gene tree inferred from their sequences. We compared the performance of these two methods on data that were generated by a computer simulation program that assumed the infinite sites model of mutation and that assumed an island model of migration. We found that in general when there is no recombination, the cladistic method performed better than FST while the reverse was true for rates of recombination similar to those found in eukaryotic nuclear genes, although FST performed better for all recombination rates for very low levels of migration (Nm = 0.1).     

Techniques for estimating genetic admixture and applications to the problem of the origin of the Icelanders and the Ashkenazi Jews

Summary A method is introduced for simultaneously using multiple loci to estimate admixutre and test goodness of fit of the model of admixture. Deviation of observed frequencies from expectation caused by sources of error such as sampling and/or drift is allowed for all loci in all populations. This allows investigation of the effects of different assumptions about sources of error on the estimates. Admixture is then investigated for Icelanders and Ashkenazi Jews. Results indicate that the Icelanders have a large Norse contribution, and that the Jews may have a small to moderate contribution from the European gene pool. There are some indications that ABO and G6PD give abnormal estimates of admixture compared to other loci, and that the Jewish gene pool may be derived from additional populations in addition to the populations considered.     

Reconciling patterns of inter-ocean molecular variance from four classes of molecular markers in blue marlin (Makaira nigricans)

Different classes of molecular markers occasionally yield discordant views of population structure within a species. Here, we examine the distribution of molecular variance from 14 polymorphic loci comprising four classes of molecular markers within approximately 400 blue marlin individuals (Makaira nigricans). Samples were collected from the Atlantic and Pacific Oceans over 5 years. Data from five hypervariable tetranucleotide microsatellite loci and restriction fragment length polymorphism (RFLP) analysis of whole molecule mitochondrial DNA (mtDNA) were reported and compared with previous analyses of allozyme and single-copy nuclear DNA (scnDNA) loci. Temporal variance in allele frequencies was nonsignificant in nearly all cases. Mitochondrial and microsatellite loci revealed striking phylogeographic partitioning among Atlantic and Pacific Ocean samples. A large cluster of alleles was present almost exclusively in Atlantic individuals at one microsatellite locus and for mtDNA, suggesting that, if gene flow occurs, it is likely to be unidirectional from Pacific to Atlantic oceans. Mitochondrial DNA inter-ocean divergence (FST) was almost four times greater than microsatellite or combined nuclear divergences including allozyme and scnDNA markers. Estimates of Neu varied by five orders of magnitude among marker classes. Using mathematical and computer simulation approaches, we show that substantially different distributions of FST are expected from marker classes that differ in mode of inheritance and rate of mutation, without influence of natural selection or sex-biased dispersal. Furthermore, divergent FST values can be reconciled by quantifying the balance between genetic drift, mutation and migration. These results illustrate the usefulness of a mitochondrial analysis of population history, and relative precision of nuclear estimates of gene flow based on a mean of several loci.     

Isolation and characterization of microsatellite loci in western North American Anodonta species

We have developed and characterized 13 microsatellite loci from a group of Anodonta species in western North America, and demonstrated their utility in populations representing two major clades in this genus. Allelic diversity and polymorphic information content were high for all loci, although these characteristics varied across populations. Deviations from Hardy-Weinberg genotypic ratios were not detected, although the estimated frequency of null alleles was high in one population for one locus. This is the first set of microsatellite loci to be developed for freshwater mussels in western North America, and will be useful for describing gene flow patterns among populations.     

Isolation and characterization of microsatellites in the bird-pollinated, autohexaploid, Eremophila glabra ssp. glabra (R.Br. (Ostenf.)) (Myoporaceae), an Australian endemic plant

Thirty-eight microsatellite loci were developed for the bird pollinated, autohexaploid, Eremophila glabra ssp. glabra. A genomic library was screened with dinucleotide and trinucleotide sequence repeats. Polymorphism ranged from one to 21 alleles per locus. Twenty-four loci exhibited null alleles, based on patterns of inheritance between maternal and progeny phenotypes. Cross-species amplification of nine Eremophila species was successful for most primers, indicating wide applicability across the genus. These microsatellites will be used to study the gene flow patterns of fragmented populations of E. glabra ssp. glabra.     

Analysis of microsatellite variation in Pinus radiata reveals effects of genetic drift but no recent bottlenecks

Most conifer species occur in large continuous populations, but radiata pine, Pinus radiata, occurs only in five disjunctive natural populations in California and Mexico. The Mexican island populations were presumably colonized from the mainland millions of years ago. According to Axelrod (1981), the mainland populations are relicts of an earlier much wider distribution, reduced some 8,000 years ago, whereas according to Millar (1997, 2000), the patchy metapopulation-like structure is typical of the long-term population demography of the species. We used 19 highly polymorphic microsatellite loci to describe population structure and to search for signs of the dynamics of population demography over space and time. Frequencies of null alleles at microsatellite loci were estimated using an approach based on the probability of identity by descent. Microsatellite genetic diversities were high in all populations [expected heterozygosity (H(e)) = 0.68-0.77], but the island populations had significantly lower estimates. Variation between loci in genetic differentiation (F(ST)) was high, but no locus deviated statistically significantly from the rest at an experiment wide level of 0.05. Thus, all loci were included in subsequent analysis. The average differentiation was measured as F(ST) = 0.14 (SD 0.012), comparable with earlier allozyme results. The island populations were more diverged from the other populations and from an inferred common ancestral gene pool than the mainland ones. All populations showed a deficiency of expected heterozygosity given the number of alleles, the mainland populations more so than the island ones. The results thus do not support a recent important contraction in the mainland range of radiata pine.     

Isolation and characterization of microsatellites in the woody shrub, Banksia sphaerocarpa var. caesia (Proteaceae)

Microsatellite markers were developed for the Australian bird-pollinated woody shrub Banksia sphaerocarpa var. caesia to study gene flow among populations in a highly fragmented landscape. Eight loci were developed, and in a sample of 40 individuals from one population, the number of alleles per locus ranged from five to 21 and observed heterozygosities ranged from 0.385 to 0.914. All eight loci showed independent inheritance. Analysis of open-pollinated progeny arrays confirmed Mendelian inheritance at seven loci, while null alleles were suspected at the remaining locus.     

Isolation and characterization of microsatellite loci from the endangered Mediterranean sponge Spongia agaricina (Demospongiae: Dictyoceratida)

The abundance of the bath sponge Spongia agaricina has decreased drastically in recent years and it is now considered an endangered species under Annex 3 of Bern and Barcelona conventions. We describe eight microsatellite markers and present data on their allelic variation and utility as high resolution genetic markers. We analyzed 36 individuals from two populations and found that the number of alleles per locus ranged between 1 and 7. Observed heterozygosity ranged from 0 to 0.72. We found deviations from Hardy–Weinberg expectations for some loci. We exclusively detected null alleles for those loci that deviated from Hardy–Weinberg expectations. Also, distributions of allele frequencies differed significantly between the two populations, making them suitable for population genetic analyses.     

Analysis of population structure in autotetraploid species.

Population structure parameters commonly used for diploid species are reexamined for the particular case of tetrasomic inheritance (autotetraploid species). Recurrence equations that describe the evolution of identity probabilities for neutral genes in an "island model" of population structure are derived assuming tetrasomic inheritance. The expected equilibrium value of FST is computed. In contrast to diploids, the correlation of genes between individuals within populations with respect to genes between populations (FST) may vary among loci due to the particular segregation patterns expected under tetrasomic inheritance and is consequently inappropriate for estimating demographic parameters in such populations. We thus define a new parameter (rho) and derive its relationship with Nm. This relationship is shown to be independent from both the selfing rate and the proportion of double reduction. Finally, the statistical procedure required to evaluate these parameters using data on gene frequencies distribution among autotetraploid populations is developed.     

Genetic diversity and population structure of an insular tree, Santalum austrocaledonicum in New Caledonian archipelago

We present a study of the genetic diversity and structure of a tropical tree in an insular system. Santalum austrocaledonicum is endemic to the archipelago of New Caledonia and is exploited for oil extraction from heartwood. A total of 431 individuals over 17 populations were analysed for eight polymorphic microsatellite loci. The number of alleles per locus ranged from 3 to 33 and the observed heterozygosity per population ranged from 0.01 in Mare to 0.74 in Ile des Pins. The genetic diversity was lowest in the most recent islands, the Loyautes, and highest in the oldest island, Grande Terre, as well as the nearby small Ile des Pins. Significant departures from panmixia were observed for some loci-population combinations (per population FIS = 0-0.03 on Grande-Terre and Ile des Pins, and 0-0.67 on Loyautes). A strong genetic differentiation among all islands was observed (FST = 0.22), and the amount of differentiation increased with geographic distance in Iles Loyaute and in Grande Terre. At both population and island levels, island age and isolation seem to be the main factors influencing the amount of genetic diversity. In particular, populations from recent islands had large average FIS that could not be entirely explained by null alleles or a Wahlund effect. This result suggests that, at least in some populations, selfing occurred extensively. Conclusively, our results indicate a strong influence of insularity on the genetic diversity and structure of Santalum austrocaledonicum.     

Global Population Genetic Structure and Male-Mediated Gene Flow in the Green Turtle (Chelonia Mydas): RFLP Analyses of Anonymous Nuclear Loci

We introduce an approach for the analysis of Mendelian polymorphisms in nuclear DNA (nDNA), using restriction fragment patterns from anonymous single-copy regions amplified by the polymerase chain reaction, and apply this method to the elucidation of population structure and gene flow in the endangered green turtle, Chelonia mydas. Seven anonymous clones isolated from a total cell DNA library were sequenced to generate primers for the amplification of nDNA fragments. Nine individuals were screened for restriction site polymorphisms at these seven loci, using 40 endonucleases. Two loci were monomorphic, while the remainder exhibited a total of nine polymorphic restriction sites and three size variants (reflecting 600-base pair (bp) and 20-bp deletions and a 20-bp insertion). A total of 256 turtle specimens from 15 nesting populations worldwide were then scored for these polymorphisms. Genotypic proportions within populations were in accord with Hardy-Weinberg expectations. Strong linkage disequilibrium observed among polymorphic sites within loci enabled multisite haplotype assignments. Estimates of the standardized variance in haplotype frequency among global collections (FST = 0.17), within the Atlantic-Mediterranean (FST = 0.13), and within the Indian-Pacific (FST = 0.13), revealed a moderate degree of population substructure. Although a previous study concluded that nesting populations appear to be highly structured with respect to female (mitochondrial DNA) lineages, estimates of Nm based on nDNA data from this study indicate moderate rates of male-mediated gene flow. A positive relationship between genetic similarity and geographic proximity suggests historical connections and/or contemporary gene flow between particular rookery populations, likely via matings on overlapping feeding grounds, migration corridors or nonnatal rookeries.     

Simultaneous Statistical Inference for Epigenetic Data

Epigenetic research leads to complex data structures. Since parametric model assumptions for the distribution of epigenetic data are hard to verify we introduce in the present work a nonparametric statistical framework for two-group comparisons. Furthermore, epigenetic analyses are often performed at various genetic loci simultaneously. Hence, in order to be able to draw valid conclusions for specific loci, an appropriate multiple testing correction is necessary. Finally, with technologies available for the simultaneous assessment of many interrelated biological parameters (such as gene arrays), statistical approaches also need to deal with a possibly unknown dependency structure in the data. Our statistical approach to the nonparametric comparison of two samples with independent multivariate observables is based on recently developed multivariate multiple permutation tests. We adapt their theory in order to cope with families of hypotheses regarding relative effects. Our results indicate that the multivariate multiple permutation test keeps the pre-assigned type I error level for the global null hypothesis. In combination with the closure principle, the family-wise error rate for the simultaneous test of the corresponding locus/parameter-specific null hypotheses can be controlled. In applications we demonstrate that group differences in epigenetic data can be detected reliably with our methodology.     

Background selection and population differentiation

A general analytical formula is derived, which predicts the effects of background selection on population differentiation at a neutral locus as a result of its linkage with selected loci of deleterious mutations. The theory is based on the assumptions of random mating, multiplicative fitness, and weak selection in hermaphrodite plants in the island model of population structure. The analytical results show that Fst at the neutral locus increases as a result of the effects of background selection, regardless of the dependence or independence among linked background selective loci. The increment in Fst is closely related to the magnitude of linkage disequilibria between the neutral locus and selected loci, and can be estimated by the ratio of Fst with background selection to Fst without background selection minus one. The steady-state linkage disequilibrium between a neutral locus and a selected locus in subpopulations, primarily attained by gene flow, decreases with the recombination rate, and can be enhanced when there are dependence among linked selected loci. Monte Carlo computer simulations with two- and three-locus models show that the analytical formulae perform well under general conditions. Application of the present theory may aid in analyzing the genome-wide mapping of the effect of background selection in terms of Fst.     

Microsatellite loci isolated from the tropical tree Hymenaea courbaril L. (Fabaceae)

We developed 11 microsatellite markers for Hymenaea courbaril for the purpose of studying spatial genetic structure and gene flow. The microsatellite loci were screened in 44 trees from two populations. All loci were polymorphic, exhibiting between two and 16 alleles, and levels of expected heterozygosity from 0.174 to 0.909. Departures from Hardy-Weinberg equilibrium were detected for all loci in one population. The estimated null allele frequency is low or moderate. No locus combinations exhibited linkage disequilibrium.