High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Water availability drives signatures of local adaptation in whitebark pine (Pinus albicaulis Engelm.) across fine spatial scales of the Lake Tahoe Basin,USA

Patterns of local adaptation at fine spatial scales are central to understanding how evolution proceeds, and are essential to the effective management of economically and ecologically important forest tree species. Here, we employ single and multilocus analyses of genetic data (n = 116 231 SNPs) to describe signatures of fine‐scale adaptation within eight whitebark pine (Pinus albicaulis Engelm.) populations across the local extent of the environmentally heterogeneous Lake Tahoe Basin, USA. We show that despite highly shared genetic variation (F_ST = 0.0069), there is strong evidence for adaptation to the rain shadow experienced across the eastern Sierra Nevada. Specifically, we build upon evidence from a common garden study and find that allele frequencies of loci associated with four phenotypes (mean = 236 SNPs), 18 environmental variables (mean = 99 SNPs), and those detected through genetic differentiation (n = 110 SNPs) exhibit significantly higher signals of selection (covariance of allele frequencies) than could be expected to arise, given the data. We also provide evidence that this covariance tracks environmental measures related to soil water availability through subtle allele frequency shifts across populations. Our results replicate empirical support for theoretical expectations of local adaptation for populations exhibiting strong gene flow and high selective pressures and suggest that ongoing adaptation of many P. albicaulis populations within the Lake Tahoe Basin will not be constrained by the lack of genetic variation. Even so, some populations exhibit low levels of heritability for the traits presumed to be related to fitness. These instances could be used to prioritize management to maintain adaptive potential. Overall, we suggest that established practices regarding whitebark pine conservation be maintained, with the additional context of fine‐scale adaptation.     

Testing hypotheses driving genetic structure in the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla

Habitat fragmentation has often been implicated in the decline of many species. For habitat specialists and/or sedentary species, loss of habitat can result in population isolation and lead to negative genetic effects. However, factors other than fragmentation can often be important and also need to be considered when assessing the genetic structure of a species. We genotyped individuals from 13 populations of the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla in Florida to test three alternative hypotheses regarding the effects that habitat fragmentation might have on genetic structure. A map of potential habitat developed from recent satellite imagery suggested that Brown‐headed Nuthatch populations in southern Florida occupied smaller and more isolated habitat patches (i.e. were more fragmented) than populations in northern Florida. We also genotyped individuals from a small, isolated Brown‐headed Nuthatch population on Grand Bahama Island. We found that populations associated with more fragmented habitat in southern Florida had lower allelic richness than populations in northern Florida (P = 0.02), although there were no differences in heterozygosity. Although pairwise estimates of F_ST were low overall, values among southern populations were generally higher than northern populations. Population assignment tests identified K = 3 clusters corresponding to a northern cluster, a southern cluster and a unique population in southeast Florida; using sampling localities as prior information revealed K = 7 clusters, with greater structure only among southern Florida populations. The Bahamas population showed moderate to high differentiation compared with Florida populations. Overall, our results suggest that fragmentation could affect gene flow in Brown‐headed Nuthatch populations and is likely to become more pronounced over time.     

Polyploid origin,genetic diversity and population structure in the tetraploid sea lavender <Emphasis Type="Italic">Limonium narbonense</Emphasis> Miller (Plumbaginaceae) from eastern Spain

Limonium narbonense Miller is a fertile tetraploid species with a sporophytic self-incompatibility system. This sea lavender is found in coastal salt marshes which have been under intense human pressure during the past decades resulting in significant habitat fragmentation. Eleven microsatellite loci specifically designed for this species were amplified in 135 individuals from five populations. These markers were used to investigate the polyploid nature, the levels of genetic diversity and population structure in this species. L. narbonense showed high levels of genetic diversity (A = 7.82, P = 100% H _T = 0.446), consistent with its likely autotetraploid origin revealed in this study and obligate outcrossing breeding system. Inbreeding (F _IS) values were low in the three southern populations (mean F _IS = 0.062), and higher in the northern populations (mean F _IS = 0.184). Bayesian analysis of population structure revealed that populations could be grouped into two genetic clusters, one including three southern populations and the other the two northernmost ones. Individuals from the two northernmost populations showed higher admixture of the two genetic clusters than individuals from the three southern ones. A thorough analysis of microsatellite electrophoretic patterns suggests an autotetraploid origin for L. narbonense. The genetic structure revealed in this study is attributed to a recent migration from the southern area. This result suggests a net gene flow from the south to the north, likely facilitated by migratory movements of birds visiting the temporary flooded ponds occupied by L. narbonense.     

Isolation and characterization of microsatellite loci in the Cape fynbos heath Erica coccinea (Ericaceae)

Enriched genomic libraries were used to isolate and characterize dinucleotide microsatellite loci in Erica coccinea, a South African Cape fynbos heath species with distinct resprouter and seeder populations. Microsatellites were required to investigate the effect of the contrasting demographic pattern driven by these two post-fire responses in the population genetic structure of seeder and resprouter forms within this species. Eight microsatellite loci were characterised and amplified a total of 106 alleles in 2 samples each of 30 individuals from 1 resprouter and 1 seeder population. Mean allele numbers were 7.88 and 11.0 for the resprouter and seeder population, respectively. Both populations showed similar average observed and expected heterozygosity levels, H _O(resprouter) = 0.683, H _O(seeder) = 0.696; H _E(resprouter) = 0.726, H _E(seeder) = 0.756, and average positive inbreeding coefficients F _IS(resprouter) = 0.058, F _IS(seeder) = 0.080. This set of microsatellite loci will be used to conduct a population genetic survey of seeder and resprouter populations throughout the range of the species. Cross-species transferability was also assayed in four other South African and four European species of the genus Erica, supporting their potential use for population genetic analyses.     

Diversity and genetic structure of the endangered cycad Dioon sonorense (Zamiaceae) from Sonora, Mexico: Evolutionary and conservation implications

The diversity and genetic structure of the cycad Dioon sonorense (De Luca, Sabato & Vázq. Torres) Chemnick, T.J. Greg. & S. Salas-Mor. were evaluated from the four known populations in tropical dry forests of Sonora, Mexico. The estimates of genetic diversity in the 19 loci analysed were: mean number of alleles per locus (A = 2.0), percentage of polymorphic loci (P = 81.6) and expected heterozygosity (H_E = 0.314). The subpopulations were found to be in Hardy–Weinberg equilibrium at the local level (F_is), even though we found a global excess of homozygotes (F_it = 0.13). Genetic variation attributed to differences among populations was 15%. Our results indicate that the geographical isolation caused by the historical effects of the Pleistocene among populations has probably generated clinal variation of the allelic frequencies at two loci, in relation to their latitudinal distribution. We compare our results with the long-lived southern hemisphere trees Fitzroya and Nothofagus in relationship to the Pleistocene glaciations. We also infer that the populations have become fragmented recently due to increasing pressure of habitat conversion, disturbance and illegal extraction of plants and seeds for the international horticultural trade. We therefore recommend three aspects: (i) propagation of this species for reintroduction, (ii) declare the localities of its northern distribution as a protected area, and (iii) an amendment of the current International Union for Conservation of Nature Red List categorisation for this species to ‘Critically Endangered’.     

Novel Y‐chromosome short tandem repeats in Sus scrofa and their variation in European wild boar and domestic pig populations

Y‐chromosome markers are important tools for studying male‐specific gene flow within and between populations, hybridization patterns and kinship. However, their use in non‐human mammals is often hampered by the lack of Y‐specific polymorphic markers. We identified new male‐specific short tandem repeats (STRs) in Sus scrofa using the available genome sequence. We selected four polymorphic loci (5–10 alleles per locus), falling in one duplicated and two single‐copy regions. A total of 32 haplotypes were found by screening 211 individuals from eight wild boar populations across Europe and five domestic pig populations. European wild boar were characterized by significantly higher levels of haplotype diversity compared to European domestic pigs (H_D = 0.904 ± 0.011 and H_D = 0.491 ± 0.077 respectively). Relationships among STR haplotypes were investigated by combining them with single nucleotide polymorphisms at two linked genes (AMELY and UTY) in a network analysis. A differentiation between wild and domestic populations was observed (F_ST = 0.229), with commercial breeds sharing no Y haplotype with the sampled wild boar. Similarly, a certain degree of geographic differentiation was observed across Europe, with a number of local private haplotypes and high diversity in northern populations. The described Y‐chromosome markers can be useful to track male inheritance and gene flow in wild and domestic populations, promising to provide insights into evolutionary and population genetics in Sus scrofa.     

Development of EST‐derived microsatellite markers for Arabidopsis lyrata subspecies petraea (L.)

A set of expressed sequence tag–simple sequence repeat (EST‐SSR) loci has been developed for Arabidopsis lyrata ssp. petraea. From 768 root cDNA clones, 126 microsatellites, including di‐, tri‐, tetra‐ and pentanucleotide repeat motifs were identified and primers were designed to 24 EST‐SSRs. Eleven loci were subsequently screened on 150 individuals sampled from five natural populations, which revealed three to nine alleles per locus (mean 5.36) and expected heterozygosity (H_E) estimates ranging from 0.046 to 0.698. Significant deviations from random mating were observed at 10 EST‐SSR loci, likely due to inbreeding (global F_IS = 0.151) and population structure (global F_ST = 0.246).     

Genetic isolation in an endemic African habitat specialist

The Chestnut‐banded Plover Charadrius pallidus is a Near‐Threatened shorebird species endemic to mainland Africa. We examined levels of genetic differentiation between its two morphologically and geographically distinct subspecies, C. p. pallidus in southern Africa (population size 11 000–16 000) and C. p. venustus in eastern Africa (population size 6500). In contrast to other plover species that maintain genetic connectivity over thousands of kilometres across continental Africa, we found profound genetic differences between remote sampling sites. Phylogenetic network analysis based on four nuclear and two mitochondrial gene regions, and population genetic structure analyses based on 11 microsatellite loci, indicated strong genetic divergence, with 2.36% mitochondrial sequence divergence between individuals sampled in Namibia (southern Africa) and those of Kenya and Tanzania (eastern Africa). This distinction between southern and eastern African populations was also supported by highly distinct genetic clusters based on microsatellite markers (global F_ST = 0.309,  = 0.510, D = 0.182). Behavioural factors that may promote genetic differentiation in this species include habitat specialization, monogamous mating behaviour and sedentariness. Reliance on an extremely small number of saline lakes for breeding and limited dispersal between populations are likely to promote reproductive and genetic isolation between eastern and southern Africa. We suggest that the two Chestnut‐banded Plover subspecies may warrant elevation to full species status. To assess this distinction fully, additional sample collection will be needed, with analysis of genetic and phenotypic traits from across the species’ entire breeding range.     

Microsatellite loci in the house fly,Musca domestica L. (Diptera: Muscidae)

Genomic libraries from house flies enriched for (CA)₁₅ and (CAG)₁₀ repeats were constructed by using biotinylated probes. Twenty‐five loci were isolated and evaluated for polymorphisms in wild flies representing two geographically diverse populations. Fourteen of 19 dinucleotide loci, and one of six trinucleotide loci were polymorphic. One hundred and twenty‐seven alleles were detected, 39 of which were private. Average number of alleles per polymorphic locus was 8.4 ± 2.5 and average heterozygosity was 72 ± 4%. F_ST by the private allele method was 0.73. Three of 15 loci showed significant heterozygote deficiencies, attributed to null alleles. Five of 15 loci were amplified in the face fly, Musca autumnalis.     

Twelve polymorphic expressed sequence tags‐derived markers for Plasmopara halstedii,the causal agent of sunflower downy mildew

Twelve expressed sequence tags‐derived markers were isolated from Plasmopara halstedii (Oomycetes), the causal agent of sunflower downy mildew. A total of 25 single nucleotide polymorphisms and five indels were detected by single‐strand conformation polymorphism analysis and developed for high‐throughput genotyping of 32 isolates. There was a high level of genetic diversity (H_E = 0.484). Observed heterozygosity ranged from 0 to 0.143 indicating that P. halstedii is probably a selfing species. These markers were also useful in detecting significant genetic variations among French populations (F_ST = 0.193) and between French and Russian populations (F_ST = 0.23). Cross‐amplification tests on three closely related species indicated that no loci amplified in other Oomycete species.     

Genetic variation amongst viviparid snails in the genus Mekongia (Prosobranchia: Viviparidae) in Thailand

Allozyme variation in viviparid snails from the genus Mekongia in Thailand were examined across the different species, subspecies and geographical locations (river drainage systems). Using horizontal starch gel electrophoresis, 11 presumed allozyme loci (eight polymorphic) from eight enzyme systems were screened. Heterozygosity was moderately low (Hexp = 0.000–0.109, mean = 0.037). One population of Mekongia pongensis (Nong Khai) was monomorphic at all 11 examined loci in contrast to the other two populations, suggesting bottleneck within this population. Populations were more differentiated in the Mekongia sphaericula complex (F_ST = 0.587) than in either the Mekongia swainsoni complex (F_ST = 0.161) or M. pongensis species (F_ST = 0.073). Mekongia sphaericula sphaericula and Mekongia sphaericula extensa exhibited fixed allele differences at two loci, a high genetic distance (D = 0.265–0.300) and a potential polyphyletic relationship, suggesting two distinct lineages (species).     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.     

Genetic structure and demographic history of Lymantria dispar (Linnaeus, 1758) (Lepidoptera: Erebidae) in its area of origin and adjacent areas

We analyzed the population genetic structure and demographic history of 20 Lymantria dispar populations from Far East Asia using microsatellite loci and mitochondrial genes. In the microsatellite analysis, the genetic distances based on pairwise F_ST values ranged from 0.0087 to 0.1171. A NeighborNet network based on pairwise F_ST genetic distances showed that the 20 regional populations were divided into five groups. Bayesian clustering analysis (K = 3) demonstrated the same groupings. The populations in the Korean Peninsula and adjacent regions, in particular, showed a mixed genetic pattern. In the mitochondrial genetic analysis based on 98 haplotypes, the median‐joining network exhibited a star shape that was focused on three high‐frequency haplotypes (Haplotype 1: central Korea and adjacent regions, Group 1; Haplotype 37: southern Korea, Group 2; and Haplotype 90: Hokkaido area, Group 3) connected by low‐frequency haplotypes. The mismatch distribution dividing the three groups was unimodal. In the neutral test, Tajima's D and Fu's FS tests were negative. We can thus infer that the Far East Asian populations of L. dispar underwent a sudden population expansion. Based on the age expansion parameter, the expansion time was inferred to be approximately 53,652 years before present (ybp) for Group 1, approximately 65,043 ybp for Group 2, and approximately 76,086 ybp for Group 3. We propose that the mixed genetic pattern of the inland populations of Far East Asia is due to these expansions and that the inland populations of the region should be treated as valid subspecies that are distinguishable from other subspecies by genetic traits.     

Unique evolutionary trajectories in repeated adaptation to hydrogen sulphide‐toxic habitats of a neotropical fish (Poecilia mexicana)

Replicated ecological gradients are prime systems to study processes of molecular evolution underlying ecological divergence. Here, we investigated the repeated adaptation of the neotropical fish Poecilia mexicana to habitats containing toxic hydrogen sulphide (H₂S) and compared two population pairs of sulphide‐adapted and ancestral fish by sequencing population pools of >200 individuals (Pool‐Seq). We inferred the evolutionary processes shaping divergence and tested the hypothesis of increase of parallelism from SNPs to molecular pathways. Coalescence analyses showed that the divergence occurred in the face of substantial bidirectional gene flow. Population divergence involved many short, widely dispersed regions across the genome. Analyses of allele frequency spectra suggest that differentiation at most loci was driven by divergent selection, followed by a selection‐mediated reduction of gene flow. Reconstructing allelic state changes suggested that selection acted mainly upon de novo mutations in the sulphide‐adapted populations. Using a corrected Jaccard index to quantify parallel evolution, we found a negligible proportion of statistically significant parallel evolution of J_corr = 0.0032 at the level of SNPs, divergent genome regions (J_corr = 0.0061) and genes therein (J_corr = 0.0091). At the level of metabolic pathways, the overlap was J_corr = 0.2545, indicating increasing parallelism with increasing level of biological integration. The majority of pathways contained positively selected genes in both sulphide populations. Hence, adaptation to sulphidic habitats necessitated adjustments throughout the genome. The largely unique evolutionary trajectories may be explained by a high proportion of de novo mutations driving the divergence. Our findings favour Gould's view that evolution is often the unrepeatable result of stochastic events with highly contingent effects.     

Noninvasive individual and species identification of jaguars (Panthera onca), pumas (Puma concolor) and ocelots (Leopardus pardalis) in Belize,Central America using cross‐species microsatellites and faecal DNA

There is a great need to develop efficient, noninvasive genetic sampling methods to study wild populations of multiple, co‐occurring, threatened felids. This is especially important for molecular scatology studies occurring in challenging tropical environments where DNA degrades quickly and the quality of faecal samples varies greatly. We optimized 14 polymorphic microsatellite loci for jaguars (Panthera onca), pumas (Puma concolor) and ocelots (Leopardus pardalis) and assessed their utility for cross‐species amplification. Additionally, we tested their reliability for species and individual identification using DNA from faeces of wild felids detected by a scat detector dog across Belize in Central America. All microsatellite loci were successfully amplified in the three target species, were polymorphic with average expected heterozygosities of H_E = 0.60 ± 0.18 (SD) for jaguars, H_E = 0.65 ± 0.21 (SD) for pumas and H_E = 0.70 ± 0.13 (SD) for ocelots and had an overall PCR amplification success of 61%. We used this nuclear DNA primer set to successfully identify species and individuals from 49% of 1053 field‐collected scat samples. This set of optimized microsatellite multiplexes represents a powerful tool for future efforts to conduct noninvasive studies on multiple, wild Neotropical felids.     

Isolation and characterization of polymorphic microsatellite loci for the study of paternity and population structure in the little penguin Eudyptula minor

We isolated and characterized eight novel microsatellite loci in the little penguin Eudyptula minor, using nonradioactive polymerase chain reaction‐based techniques to screen GA and GAAA repeats from enriched genomic DNA libraries. All eight loci were polymorphic and seven were variable in our main study population (mean H_E = 0.613, mean N_A = 7.14). Cross‐amplification using a microsatellite primer developed in Spheniscus demersus (African penguin) yielded one additional polymorphic locus. This locus combined with six of the little penguin loci is suitable for paternity assignment in little penguins (exclusion probability for seven unlinked loci = 0.993).     

Isolation and characterization of microsatellite markers in musk duck (Biziura lobata: Aves), and their application to other waterfowl species

We isolated and characterized 11 novel microsatellite loci to study paternity in the Australian musk duck (Biziura lobata), using nonradioactive PCR‐based techniques to screen GA and GAAA repeats enriched genomic DNA libraries. Nine of 11 loci showed no evidence of null alleles and were variable (mean H_E = 0.825, mean number of alleles = 9). This set of nine loci is suitable for paternity assignment (exclusion probability for nine unlinked loci = 0.9999). We also demonstrated that many of these loci cross‐amplify in various other waterfowl species.