A comparison of four methods for detecting weak genetic structure from marker data

Genetic structure is ubiquitous in wild populations and is the result of the processes of natural selection, genetic drift, mutation, and gene flow. Genetic drift and divergent selection promotes the generation of genetic structure, while gene flow homogenizes the subpopulations. The ability to detect genetic structure from marker data diminishes rapidly with a decreasing level of differentiation among subpopulations. Weak genetic structure may be unimportant over evolutionary time scales but could have important implications in ecology and conservation biology. In this paper we examine methods for detecting and quantifying weak genetic structures using simulated data. We simulated populations consisting of two putative subpopulations evolving for up to 50 generations with varying degrees of gene flow (migration), and varying amounts of information (allelic diversity). There are a number of techniques available to detect and quantify genetic structure but here we concentrate on four methods: F(ST), population assignment, relatedness, and sibship assignment. Under the simple mating system simulated here, the four methods produce qualitatively similar results. However, the assignment method performed relatively poorly when genetic structure was weak and we therefore caution against using this method when the analytical aim is to detect fine-scale patterns. Further work should examine situations with different mating systems, for example where a few individuals dominate reproductive output of the population. This study will help workers to design their experiments (e.g., sample sizes of markers and individuals), and to decide which methods are likely to be most appropriate for their particular data.     

Conservation genetics of insular Podarcis lizards using partial cytochrome b sequences

Sequence data derived from a 306 bp fragment of mitochondrial cytochrome b and molecular variance estimates were used to investigate the genetic population structure of the endangered and endemic lizard Podarcis atrata of the Columbretes archipelago (Mediterranean, Spain). Our results show a very high and significant among-population genetic differentiation. F _ST values and phylogenetic analyses confirm the evolutionary distinctiveness of P. atrata populations, suggesting that the populations of these islands deserve special protection measures. The populations of the two islands Columbrete Grande and Mancolibre are less differentiated than those of Foradada and Lobo, and seem to have retained mainland haplotypes. This situation needs further attention as the origin of the mainland haplotypes is still unclear. If they are a result of recent introductions from mainland specimens, then they may represent a threat to the endemic lizards of the Columbretes islands.     

SMALL POPULATION GENETIC VARIABILITY AT LOCI UNDER STABILIZING SELECTION

Genetic variability at a locus under stabilizing selection in a finite population is investigated using analytic methods and computer simulations. Three measures are examined: the number of alleles k, heterozygosity H, and additive genetic variance Vg. A nearly-neutral theory results. The composite parameter S = NV_M/V_s (where N is the population size, V_M the variance of new mutant allelic effects and V_s the weakness of stabilizing selection) figures prominently in the results. The equilibrium heterozygosity is similar to that of strictly neutral theory, H = 4N_μc/ (1 + 4N_μc), except that μ_c = where c is about 0.5. Simulations corroborate except for very low N. Genetic variability attains similar equilibrium values at both a “lone” locus and at an “embedded” locus. This agrees with my earlier work concerning molecular clock rates. These results modify the neutralist interpretation of data concerning genetic variability and genetic distances between populations. Low H values are proportional not to N but to . This may explain the narrow observed range of H among species. Heterozygosities need not be highly correlated to genetic variances. Genetic variances are not highly dependent on population size except in very small populations which are difficult to sample without bias because the smallest populations go extinct the fastest. Nearly neutral evolution will not be easily distinguished from strictly neutral theory under the Hudson-Kreitman-Aguade inter-/intraspecific variation ratio test, since a similar effective mutation rate holds for genetic distances and D = 2μ_ct, where . As with strictly neutral theory, comparisons across loci should show D and H to be positively correlated because of the shared μ_c. But unlike neutral theory, for a given locus, comparisons across species should show D and H to be negatively correlated. There is no obvious threshold of population size below which genetic variability inevitably declines. Extinction depends on both genetic variation and natural selection. Neither theory nor observation presently indicates the measure of genetic variability (k, H, V_G or other) that best indicates vulnerability of a small population to extinction.     

Genetic differentiation and connectivity of morphological types of the broadcast‐spawning coral Galaxea fascicularis in the Nansei Islands,Japan

Population connectivity resulting from larval dispersal is essential for the maintenance or recovery of populations in marine ecosystems, including coral reefs. Studies of species diversity and genetic connectivity within species are essential for the conservation of corals and coral reef ecosystems. We analyzed mitochondrial DNA sequence types and microsatellite genotypes of the broadcast‐spawning coral, Galaxea fascicularis, from four regions in the subtropical Nansei Islands in the northwestern Pacific Ocean. Two types (soft and hard types) of nematocyst morphology are known in G. fascicularis and are significantly correlated with the length of a mitochondrial DNA noncoding sequence (soft type: mt‐L; hard type: mt‐S type). Using microsatellites, significant genetic differentiation was detected between the mitochondrial DNA sequence types in all regions. We also found a third genetic cluster (mt‐L+), and this unexpected type may be a cryptic species of Galaxea. High clonal diversity was detected in both mt‐L and mt‐S types. Significant genetic differentiation, which was found among regions within a given type (F_ST = 0.009–0.024, all Ps ≤ 0.005 in mt‐L; 0.009–0.032, all Ps ≤ 0.01 in mt‐S), may result from the shorter larval development than in other broadcast‐spawning corals, such as the genus Acropora. Nevertheless, intraspecific genetic diversity and connectivity have been maintained, and with both sexual and asexual reproduction, this species appears to have a potential for the recovery of populations after disturbance.     

Allozyme variation in the three extant populations of the narrowly endemic cycad Dioon angustifolium Miq. (Zamiaceae) from North-eastern Mexico

BACKGROUND AND AIMS: Dioon angustifolium was considered within D. edule. Recent morphometric and allozyme studies on D. edule have shown that D. angustifolium has originated from geographic isolation and is therefore considered to be a separate species. This cycad is endemic to north-eastern Mexico and is known only from three populations in the Sierra Madre Oriental mountain chain. Its populations are small when compared with its southern relative D. edule. In this study, genetic variation was determined within and between populations of D. angustifolium and the genetic consequences of habitat fragmentation and isolation of populations of this species were assessed. METHODS: Allozyme electrophoresis of 14 presumptive loci was used. The data were analysed with statistical approximations for estimating genetic diversity, structure, gene flow and recent genetic bottlenecks. KEY RESULTS: Means and standard deviations of genetic diversity estimators were: number of alleles per locus (A = 1.67 +/- 0.23), percentage of polymorphic loci (P = 52.4 +/- 23 %) and expected heterozygosity (H(E) = 0.218 +/- 0.093). The genetic variation attributable to differences among populations was 16.7 %. Mean gene flow between paired populations was Nm = 1.55 +/- 0.67, which is similar to that reported for endemic plant species of narrow geographical distribution and species with gravity-dispersed seed. A recent bottleneck is detected in the populations studied. CONCLUSIONS: Dioon angustifolium presents high levels of genetic diversity compared with other cycad species, in spite of small population sizes. The recent bottleneck effect did not effectively reduce the genetic variation to the extent of eliminating these populations. The distribution of D. angustifolium appears to be the result of historical biogeographical effects related to the Pleistocene glaciations. It is recommended that this species be catalogued in the IUCN Red List of Threatened Species and conservation efforts be made to preserve it.     

Gene flow and range expansion in a mountain‐dwelling passerine with a fragmented distribution

We studied gene flow and bottleneck events in the population history of locally isolated citril finches endemic to European mountains. For the present study, we used two genetic markers with different rates of evolution: a fast evolving mitochondrial marker (ATPase6/8) and a more slowly evolving nuclear marker (02401). Populations north of the Pyrenees showed in general fewer haplotypes and a considerable lower nucleotide and gene diversity than the Iberian populations. Unexpectedly, we found very little genetic variability in the fast evolving mitochondrial marker, arguing for a strong and relatively recent bottleneck event in the species population history. This pattern potentially reflects a sudden decrease of crucial resources during Mid‐Holocene (mountain pine, Scots pine, and black pine) and a subsequent breakdown of the population. The bottleneck could also have been caused or coincide with a selective sweep in the mitochondrion. By contrast, the slowly evolving nuclear marker showed a much higher variability. This marker probably reflects major gene flow along a potential expansion pathway from the Eastern Pyrenees, northwards to the populations of Central Europe, and southwards to the more fragmented populations of central and southern Spain. The population of the Western Pyrenees (Navarra) appears to be cut‐off from this major gene flow and our data indicate a certain degree of partial isolation, probably reflecting more ancient events (e.g. the separation in distinct refuge sites during the last glacial maximum). © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103 , 707–721.     

Effects of multiple founder populations on spatial genetic structure of reintroduced American martens

Reintroductions and translocations are increasingly used to repatriate or increase probabilities of persistence for animal and plant species. Genetic and demographic characteristics of founding individuals and suitability of habitat at release sites are commonly believed to affect the success of these conservation programs. Genetic divergence among multiple source populations of American martens (Martes americana) and well documented introduction histories permitted analyses of post‐introduction dispersion from release sites and development of genetic clusters in the Upper Peninsula (UP) of Michigan <50 years following release. Location and size of spatial genetic clusters and measures of individual‐based autocorrelation were inferred using 11 microsatellite loci. We identified three genetic clusters in geographic proximity to original release locations. Estimated distances of effective gene flow based on spatial autocorrelation varied greatly among genetic clusters (30–90 km). Spatial contiguity of genetic clusters has been largely maintained with evidence for admixture primarily in localized regions, suggesting recent contact or locally retarded rates of gene flow. Data provide guidance for future studies of the effects of permeabilities of different land‐cover and land‐use features to dispersal and of other biotic and environmental factors that may contribute to the colonization process and development of spatial genetic associations.     

Genetic management of the red squirrel,Sciurus vulgaris: a practical approach to regional conservation

The progressive decline in red squirrel (Sciurus vulgaris) numbers in Wales has led to conservation and reintroduction projects being established on the island of Anglesey. The recovery of the island’s remnant wild population was initially successful, however concern remained over potential loss of genetic diversity resulting from an observed demographic bottleneck. We used mitochondrial DNA (mtDNA) control region sequences and six microsatellite loci to assess current levels of genetic variation in the population. Samples were monomorphic for control region sequences and a historic specimen from the same area carrying a different haplotype demonstrated a loss of mtDNA diversity during the last 20 years. Inclusion of other Welsh haplotypes indicated phylogeographic structure in the region, in contrast to previous UK studies. Genotyping results showed allelic diversity and heterozygosity to be less than 50% of that recorded in other UK populations, with strong evidence for a recent genetic bottleneck. A parallel reintroduction programme on Anglesey included genetic analysis of individuals during the selection of captive breeding pairs. We present analysis of sequence and microsatellite data, and subsequent management decisions taken to maximise diversity in the founder and F1 generations. Population and Habitat Viability Analysis applied to both populations modelled future levels of heterozygosity and allelic diversity. Supplementation of the remnant and reintroduced populations with translocated squirrels was simulated as a potential management tool; results support use of this strategy to reduce loss of diversity and increase survival. The limitations of applying conservation genetic theory within small-scale management projects are discussed.