Population differentiation in Trinidadian guppies (Poecilia reticulata): patterns and problems

Populations of the guppy, Poecilia reficulata , in N. Trinidad exhibit marked population differentiation in allozyme frequencies. Here we investigate six further populations electrophoretically at 25 enzyme-coding loci to examine patterns in geographical structuring, genotypic distributions and genetic diversity. With one exception, possibly related to an experimental introduction, populations divided broadly (dendrogram of Nei's mean genetic identity, Ī) in accordance with proposed ancestral colonization. Most populations were in Hardy Weinberg equilibrium, though some significant deficits in heterozygotes were detected. Incorporating information from published data, markedly hjgher levels of genetic diversity (mean observed heterozygosity, _o) were recorded in lowland [_o= 0.0382 ± k0406 ( s.e .), n = 9] compared with upstream populations [_o= 0.0112 ± 00034 (S.E.), n = 9]. Patterns are discussed in relation to historical and present-day evolutionary forces.     

Some morphometric and biochemical features of ready-to-migrate silver and pre-migratory yellow stages of the shortfin eel of south-eastern Australian waters

The mean total length ( L _T), mass and age of ready to migrate female silver shortfin eels Anguilla australis from the Hopkins River estuary and the mouth of the Merri River in south-eastern Australia, were 83·2 ± 1·2 cm, 1051 ± 51 g, and 17·2 ± 1·79 years, respectively. The eye index ( I _E) of the silver shortfin eels was < 5·2 (mean 7·64 ± 0·29) and differed significantly from that of the yellow shortfin eels collected from two other sites. The I _E increased with L _T (mm) and was related by log I _E= 2·656 log L _T6·925. The per cent moisture, protein and ash content of the liver of silver shortfin eels was significantly lower than in yellow shortfin eels, but lipid content was significantly higher in the former (35·5 ± 2·0%). The mean mass μg mg lipid ⁻) of saturates (230·4 ± 2·6 v . 181·7 ±2·6), monoenes (367·4 ± 6·3 v . 290·8 ± 8·9) and PUFA (177·3 ± 5·3 v . 159·7 ± 4·6) in muscle was significantly higher, and the great majority of individual fatty acids was found also in higher quantities in silver shortfin eels. In the liver, the PUFA found in the highest quantity was 22:6n-3, except in shortfin eels from Hopkins River estuary, and the amount of 18:2n-6 in the liver of silver shortfin eels was significantly higher than that in yellow shortfin eels but the reverse was true of 20:4n-6. In both muscle and liver tissues the saturate 16:0 and the monoene 18:ln-9 collectively accounted for >50% of all the fatty acids in the lipid.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.     

Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

Temporal effective size estimates of a managed walleye Sander vitreus population and implications for genetic-based management

The goal of this research was to use the long-term fishery data set and DNA from archived scales of walleye Sander vitreus in Escanaba Lake, WI, U.S.A., to improve the understanding of the underlying mechanism(s) influencing genetic diversity in naturally recruiting populations. The introduced population of S. vitreus in Escanaba Lake has a low mean effective population size ( N _E) between 124·6 and 185·5 despite a mean census size ( N _C) of 4659 ( N _E/ N _C c. 0·04), suggesting an accelerated rate of genetic drift between 1952 and 2002. These values are smaller than the median N _E range of several studies suggesting typical N _E/ N _C ratios of 0·11–0·16 in a wide range of taxa. N _E increased steadily during the past two sampled decades (1992 and 2002) and was consistent with a lowering of the variance in S. vitreus reproductive success, possibly linked to a large, sustained exploitation (mean 28%) rate. Variance in reproductive success is one of the most important factors influencing N _E in species, like S. vitreus , which have a potential for large fecundities and large juvenile mortalities (type III survivorship). The N _B estimates across six sequential cohorts (age classes of S. vitreus , assayed from 1994 to 1999) was consistent with estimates of N _E reported for 1992–2002. These results, coupled with in-depth census and exploitation data, show that the genetic characteristics of Escanaba Lake S. vitreus have changed substantially and that management activities, such as supplemental stocking and harvest practices, have profoundly influenced the genetic dynamics of S. vitreus in this lake.     

Analysis of the enzyme polymorphism in a plaice, Pleuronectes platessa L., population from the north-west coast of Brittany, France

A study of enzyme polymorphism by electrophoretical techniques permitted the genetic structure and the homogeneity degree of the plaice population present in the Abers Wrac'h and Benoit along the north-west coast of Brittany to be determined. Thirteen enzymic systems encoded by 21 loci were analysed and nine of these loci were polymorphic at the 5% level. The comparison of allelic frequencies and of observed heterozygosity h _obs did not show significant differences between these two individual groups. Nei's genetic identity ( I = 0.998) and distance ( D = 0.002) as estimated from allelic frequencies showed that the plaice of the two estuaries are very similar.
Thus, the plaice of A. Wrac'h and of A. Benoit belong to the same population in which the average level of polymorphism per locus P is 47% and the observed average heterozygosity per individual H _obs is 0.163 ± 0.094.     

Allozyme Diversity in Populations of Cymbidium goeringii (Orchidaceae)

Abstract: Using 14 allozyme loci, we investigated levels of genetic diversity within populations, and degree of genetic divergence among 24 populations of Cymbidium goeringii (Orchidaceae) in Korea and Japan. Cymbidium goeringii maintains high levels of genetic diversity both at population (mean expected heterozygosity, H _e = 0.238) and species levels (0.260). Means of H _e found in 24 populations were not significantly different from each other. About 90 % of the total variation in the species is common to all populations (mean G _ST = 0.108). No unique allele was found in any population. The indirect estimate of gene flow based on the mean G _ST was high ( Nm = 2.06). Nei's genetic identities for pairs of populations had high values (mean = 0.974 [SD = 0.013]). The Mantel-Z test showed a significant correlation between genetic distance and geographic distance. However, the mean G _ST value between 17 populations in Korea and seven Japanese populations was relatively low (0.029), even though the land connection between the southern Korean peninsula and southern Japanese archipelagos has not existed since the middle Pleistocene. Large numbers of small seeds of C. goeringii might travel long distances by wind from populations to populations both in Korea and Japan, increasing genetic diversity within populations and maintaining low genetic differentiation among populations.     

Genetic differentiation between Australian and North American populations of the monarch butterfly Danaus plexippus (L.) (Lepidoptera: Nymphalidae): an exploration using allozyme electrophoresis

Allozyme analysis was used to address the question of the source of the Australian populations of the monarch butterfly Danaus plexippus (L.). The study had three major aims: (1) To compare the levels of diversity of Australian and Hawaiian populations with potential source populations. (2) To determine whether eastern and western North American populations were sufficiently divergent for the Australian populations to be aligned to a source population. (3) To compare the differentiation among regions in Australia and North America to test the prediction of greater genetic structure in Australia, as a consequence of reduced migratory behaviour. The reverse was found, with F _ST values an order of magnitude lower in Australia than in North America. Predictably, Australian and Hawaiian populations had lower allelic diversity, but unexpected higher heterozygosity values than North American populations. It was not possible to assign the Australian populations to a definitive source, although the high levels of similarity of Australian populations to each other suggest a single colonization event. The possibility that the Australian populations have not been here long enough to reach equilibrium is discussed. © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 75 , 437–452.     

Genetic structure within and among regional populations of the Eurasian badger (Meles meles) from Denmark and the Netherlands

The Eurasian badger Meles meles has a wide distribution area ranging from Japan to Ireland. In western Europe badger habitats are severely disturbed by anthropogenic factors, leading to fragmentation into subpopulations and formation of a metapopulation substructuring of once continuous panmictic populations. We have examined the genetic structure of Dutch and Danish badger populations on a relatively small scale (within countries) and a larger scale (between countries). The levels of genetic variation of populations were moderate and did not differ significantly among populations (overall H _O=0.30, overall H _E=0.34). Considerable genetic differentiation between the Dutch and Danish populations was found (overall F _ST=0.32, mean pairwise Dutch–Danish F _ST=0.42), indicating a large-scale substructuring of these western European badger populations. Further analysis showed that the Danish badger population can be substructured into three clusters [ P _{( k =3)}=0.99], but the Dutch populations cluster into one more or less panmictic population [ P _{( k =1)}=0.99] with little or no substructuring. The presence of migration barriers, such as roads, together with the peninsular geography of Denmark, may have led to this structuring of badger populations. In contrast, measures that improve migration and connection to other populations from neighboring countries may have prevented substructuring of the Dutch badger population.     

The distribution of random amplified polymorphic DNA (RAPD) diversity amongst populations of Isotoma petraea (Lobeliaceae)

RAPDs were generated from plants of six populations of Isotoma petraea F. Muell. The species occurs on rock outcrops in southern and western Australia, with populations exhibiting different breeding systems, including complete autogamy, varying levels of outbreeding and complex hybridity. Non-metric multidimensional scaling (nMDS) analysis of the random amplified polymorphic DNA (RAPD) data set clearly resolved all populations. The Pigeon Rock population, which is home to both complex hybrid and structural homozygote plants, was divided into those two groups by the nMDS analysis. There was little diversity in highly autogamous populations, but levels were higher in the outbred Yackeyackine population. All complex hybrid populations and plants possessed numerous genetic system-specific RAPDs, some of which were shown to be held in fixed heterozygosity. Estimating G _ST using RAPDs has been problematical due to their dominance, and analytical methods usually rely on knowledge of the selfing rate or assume Hardy–Weinberg equilibrium. This assumption does not hold when populations exhibit fixed heterozygosity, and an alternative method, Shannon's Index, was used to partition genetic diversity. The distribution of genetic diversity fit expectations for an inbreeding species, with most of the variation (87.5%) occurring between populations. This compares to an average RAPD-based G _ST of 59.6% for inbreeding species generally and 15.5% for outbreeding species.     

Microsatellite analysis of genetic population structure in an endangered salmonid: the coastal cutthroat trout (Oncorhynchus clarki clarki)

The genetic population structure of coastal cutthroat trout ( Oncorhynchus clarki clarki ) in Washington state was investigated by analysis of variation in allele frequencies at six highly polymorphic microsatellite loci for 13 anadromous populations, along with one outgroup population from the Yellowstone subspecies ( O. clarki bouvieri) (mean heterozygosity = 67%; average number of alleles per locus = 24). Tests for genetic differentiation revealed highly significant differences in genotypic frequencies for pairwise comparisons between all populations within geographical regions and overall population subdivision was substantial ( F _ST = 0.121, R _ST = 0.093), with 44.6% and 55.4% of the among-population diversity being attributable to differences between streams ( F _SR = 0.054) and between regions ( F _RT = 0.067), respectively. Analysis of genetic distances and geographical distances did not support a simple model of isolation by distance for these populations. With one exception, neighbour-joining dendrograms from the Cavalli-Sforza and Edwards' chord distances and maximum likelihood algorithms clustered populations by physiogeographic region, although overall bootstrap support was relatively low (53%). Our results suggest that coastal cutthroat trout populations are ultimately structured genetically at the level of individual streams. It appears that the dynamic balance between gene flow and genetic drift in the subspecies favours a high degree of genetic differentiation and population subdivision with the simultaneous maintenance of high heterozygosity levels within local populations. Results are discussed in terms of coastal cutthroat trout ecology along with implications for the designation of evolutionarily significant units pursuant to the US Endangered Species Act of 1973 and analogous conservation units.     

The patterns of genetic diversity in Carpentaria acuminata (Arecaceae), and rainforest history in northern Australia

Carpentaria acuminata occurs in monsoon rainforest and is endemic to the Northern Territory, Australia. The genetic diversity of C. acuminata populations was surveyed across the geographical range of the species using isozyme analysis. Genetic diversity within C. acuminata populations ( H _E = 0.143) was typical of rainforest species and woody angiosperms generally. Genetic diversity was not correlated with rainforest patch size. However, there was significant heterogeneity among populations ( F _ST = 0.379), with infrequent effective gene flow among populations ( Nm = 0.39). Genetic diversity was negatively correlated with increasing distance between neighbouring C. acuminata populations, but geographical distance was not a good predictor of genetic similarity. C. acuminata is a favoured food of mobile frugivores such as Torres Strait pigeons and flying foxes. The decreased diversity with decreasing density of populations indicated that seed dispersal by frugivores has been important for the maintenance of diversity in this species. Populations known to have originated on relatively young, Holocene landforms were not necessarily genetically depauperate. Gene flow by pollen is apparently limited because C. acuminata populations are significantly inbred regardless of genetic diversity ( F = 0.641). The distribution and diversity of rare alleles, i.e. those occurring in few populations, is consistent with the theory of rainforest contraction during the Pleistocene.     

Population genetic structure of wild rice Oryza glumaepatula distributed in the Amazon flood area influenced by its life-history traits

Wild plant species develop their own way of living to adapt to the specific environment of their habitats. Their life-history traits strongly affect the genetic structure of the population. The wild species Oryza glumaepatula Steud. growing in the Amazon basin seems to have characteristic life-history traits suited for the flood condition. At the vegetative growth stage, the culms frequently break at internodes. With no roots anchoring on the ground, plant bodies floating in the water move downriver by water current and wind. To examine the association between the life-history traits and genetic population structure of Amazonian O. glumaepatula , we analysed allozyme variability at 29 loci of 16 enzymes using 37 populations from five regions. Allozymes were not so variable (total gene diversity H _E = 0.044) compared with Asian wild rice, O. rufipogon Griff. The bottleneck effect and rare opportunity of interspecies gene flow may prevent the development of allozyme variability. Population genotypes tended to be differentiated among geographically isolated regions. Observed heterozygosities were much lower than expected heterozygosities, or gene diversity ( H _O = 0.003 for whole population) and F _IS over polymorphic loci was 0.931, indicating that O. glumaepatula has developed an inbreeding system. But, the intrapopulation gene diversity ( H _S) was higher than interpopulation gene diversity ( D _ST), as generally observed in outbreeding populations. The migration ability of O. glumaepatula makes long-distance seed dispersal possible. This might have led to frequent gene flow among populations.     

Population genetic structure in Myrtus communis L. in a chronically fragmented landscape in the Mediterranean: can gene flow counteract habitat perturbation?

Ancient managed landscapes provide ideal opportunities to assess the consequences of habitat fragmentation on the patterns of genetic diversity and gene flow in long-lived plant species. Using amplified fragment length polymorphism (AFLP) and allozyme markers, we quantified seed-mediated gene flow and population genetic diversity and structure in 14 populations of Myrtus communis (myrtle), a common endozoochorous shrub species of forest patches in lowland agricultural Mediterranean areas. Overall, allozyme diversity for myrtle was low (P₉₅   =   25%; A   =   1.411; H_e = 0.085) compared to other known populations, and a significant portion of populations (57%) had lower levels of allelic diversity and/or heterozygosity than expected at random, as shown by simulated resampling of the whole diversity of the landscape. We found significant correlations between allozyme variability and population size and patch isolation, but no significant inbreeding in any population. Genetic differentiation among populations for both allozyme and AFLP markers was significant (Φ_ST = 0.144 and Φ_ST = 0.142, respectively) but an isolation-by-distance pattern was not detected. Assignment tests on AFLP data indicated a high immigration rate in the populations ( ca. 20–22%), likely through effective seed dispersal across the landscape by birds and mammals. Our results suggest that genetic isolation is not the automatic outcome of habitat destruction since substantial levels of seed-mediated gene flow are currently detectable. However, even moderate rates of gene flow seem insufficient in this long-lived species to counteract the genetic erosion and differentiation imposed by chronic habitat destruction.     

Genetic variation in BoLA microsatellite loci in Portuguese cattle breeds

Major histocompatibility complex (MHC) typing based on microsatellites can be a valuable approach to understanding the selective processes occurring at linked or physically close MHC genes and can provide important information on variability and relationships of populations. Using microsatellites within or in close proximity with bovine lymphocyte antigen (BoLA) genes, we investigated the polymorphisms in the bovine MHC, known as the BoLA, in eight Portuguese cattle breeds. Additional data from non-BoLA microsatellite loci were also used to compare the variability between these regions. Diversity was higher in BoLA than in non-BoLA microsatellites, as could be observed by the number of alleles, allelic richness and observed heterozygosity. Brava de Lide, a breed selected for aggressiveness and nobility, presented the lowest values of observed heterozygosity and allelic richness in both markers. Results from neutrality tests showed few statistically significant differences between the observed Hardy–Weinberg homozygosity ( F ) and the expected homozygosity ( F _E), indicating the apparent neutrality of the BoLA microsatellites within the analysed breeds. Nevertheless, we detected a trend of lower values of observed homozygosity compared with the expected one. We also detected some differences in the levels of allelic variability among the four BoLA microsatellites. Our data showed a higher number of alleles at the BoLA-DRB3 locus than at the BoLA-DRBP1 locus. These differences could be related to their physical position in the chromosome and may reflect functional requirements for diversity.     

Small population size and extremely low levels of genetic diversity in island populations of the platypus, Ornithorhynchus anatinus

Genetic diversity generally underpins population resilience and persistence. Reductions in population size and absence of gene flow can lead to reductions in genetic diversity, reproductive fitness, and a limited ability to adapt to environmental change increasing the risk of extinction. Island populations are typically small and isolated, and as a result, inbreeding and reduced genetic diversity elevate their extinction risk. Two island populations of the platypus, Ornithorhynchus anatinus, exist; a naturally occurring population on King Island in Bass Strait and a recently introduced population on Kangaroo Island off the coast of South Australia. Here we assessed the genetic diversity within these two island populations and contrasted these patterns with genetic diversity estimates in areas from which the populations are likely to have been founded. On Kangaroo Island, we also modeled live capture data to determine estimates of population size. Levels of genetic diversity in King Island platypuses are perilously low, with eight of 13 microsatellite loci fixed, likely reflecting their small population size and prolonged isolation. Estimates of heterozygosity detected by microsatellites (H(E)= 0.032) are among the lowest level of genetic diversity recorded by this method in a naturally outbreeding vertebrate population. In contrast, estimates of genetic diversity on Kangaroo Island are somewhat higher. However, estimates of small population size and the limited founders combined with genetic isolation are likely to lead to further losses of genetic diversity through time for the Kangaroo Island platypus population. Implications for the future of these and similarly isolated or genetically depauperate populations are discussed.     

Does Species Range and Rarity Affect Population Genetics? A Case Study of Four Graptophyllum Species from Queensland,Australia

There are four species of Graptophyllum (Acanthaceae) in Australia which span a gradient in rarity and conservation status from the most endangered, Graptophyllum reticulatum (endangered), to G. ilicifolium (vulnerable), to G. excelsum (rare), to the least endangered G. spinigerum (widespread). Graptophyllum reticulatum and G. ilicifolium are each only known from three locations. All species were found to have quite high levels of genetic diversity (H_e, A, A_p, P) and there was no significant difference among species in genetic diversity ( P > 0.05), thus, genetic diversity was unrelated to conservation status or species geographic range. All species had high allelic fixation values ( F  ), indicating that populations were effectively inbred regardless of conservation status and geographic range. The genetic diversity among populations (F_ST) was not related to rarity and contrasted among the most restricted species, where diversity among populations was twice as high among G. reticulatum populations (F_ST, 0.22) compared with among G. ilicifolium populations (F_ST, 0.11).     

Genetic diversity and effective size of the Atlantic salmon Salmo salar L. inhabiting the River Eo (Spain) following a stock collapse

Historic angling records suggest the occurrence of a drastic decline in the River Eo (Spain) Atlantic salmon population size during the past two decades, as a result of overexploitation and habitat deterioration. In recent years, the population has been apparently recovering, and the present study is aimed to report information on the level of genetic diversity and the effective size of the current population as these may have immediate consequences for its conservation. Eighty-six salmon from two temporal groups (1998–1999 and 2004–2006), representing three generations, were genotyped using a panel of eight microsatellites. Inspite of the recent decline in census numbers and the detection of the signs of a population bottleneck, the population exhibits a high level of genetic diversity, similar to that from other populations, and almost unchanged during the period of study [average allelic richness ( A ) = 14·0 and 13·9, and average heterozygosity ( H _e) = 0·843 and 0·851 in 1998–1999 and 2004–2006, respectively]. The effective population size ( N _e) estimated by two different temporal methods showed a consistent value around 80 salmon, whereas the estimates from the linkage disequilibrium (LD) method provided a value around 165 individuals for either sample. The recent growing number of salmon, as indicated by fisheries records, the relatively large estimates of the ratio N _e/ N (with range 0·23–0·44 for the temporal estimates and 0·31–0·59 for the LD estimates) and the high levels of diversity found suggest that the population has not been greatly affected by the historical census declines and can be expected to recover in the future.     

Evolutionary implications of a high selfing rate in the freshwater snail Lymnaea truncatula

Abstract Self-compatible hermaphroditic organisms that mix self-fertilization and outcrossing are of great interest for investigating the evolution of mating systems. We investigate the evolution of selfing in Lymnaea truncatula , a self-compatible hermaphroditic freshwater snail. We first analyze the consequences of selfing in terms of genetic variability within and among populations and then investigate how these consequences along with the species ecology (harshness of the habitat and parasitism) might govern the evolution of selfing. Snails from 13 localities (classified as temporary or permanent depending on their water availability) were sampled in western Switzerland and genotyped for seven microsatellite loci. F_IS (estimated on adults) and progeny array analyses (on hatchlings) provided similar selfing rate estimates of 80%. Populations presented a low polymorphism and were highly differentiated (F_ST= 0.58). Although the reproductive assurance hypothesis would predict higher selfing rate in temporary populations, no difference in selfing level was observed between temporary and permanent populations. However, allelic richness and gene diversity declined in temporary habitats, presumably reflecting drift. Infection levels varied but were not simply related to either estimated population selfing rate or to differences in heterozygosity. These findings and the similar selfing rates estimated for hatchlings and adults suggest that within-population inbreeding depression is low in L. truncatula.     

Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce

Conifers are characterized by a large genome size and a rapid decay of linkage disequilibrium, most often within gene limits. Genome scans based on noncoding markers are less likely to detect molecular adaptation linked to genes in these species. In this study, we assessed the effectiveness of a genome-wide single nucleotide polymorphism (SNP) scan focused on expressed genes in detecting local adaptation in a conifer species. Samples were collected from six natural populations of white spruce ( Picea glauca ) moderately differentiated for several quantitative characters. A total of 534 SNPs representing 345 expressed genes were analysed. Genes potentially under natural selection were identified by estimating the differentiation in SNP frequencies among populations ( F _ST) and identifying outliers, and by estimating local differentiation using a Bayesian approach. Both average expected heterozygosity and population differentiation estimates ( H _E = 0.270 and F _ST = 0.006) were comparable to those obtained with other genetic markers. Of all genes, 5.5% were identified as outliers with F _ST at the 95% confidence level, while 14% were identified as candidates for local adaptation with the Bayesian method. There was some overlap between the two gene sets. More than half of the candidate genes for local adaptation were specific to the warmest population, about 20% to the most arid population, and 15% to the coldest and most humid higher altitude population. These adaptive trends were consistent with the genes' putative functions and the divergence in quantitative traits noted among the populations. The results suggest that an approach separating the locus and population effects is useful to identify genes potentially under selection. These candidates are worth exploring in more details at the physiological and ecological levels.