Genetic drift outweighs natural selection at toll‐like receptor (TLR) immunity loci in a re‐introduced population of a threatened species

During population establishment, genetic drift can be the key driver of changes in genetic diversity, particularly while the population is small. However, natural selection can also play a role in shaping diversity at functionally important loci. We used a well‐studied, re‐introduced population of the threatened Stewart Island robin (N = 722 pedigreed individuals) to determine whether selection shaped genetic diversity at innate immunity toll‐like receptor (TLR) genes, over a 9‐year period of population growth following establishment with 12 genetic founders. We found no evidence for selection operating with respect to TLR diversity on first‐year overwinter survival for the majority of loci, genotypes and alleles studied. However, survival of individuals with TLR4_BE genotype was significantly improved: these birds were less than half as likely to die prior to maturity compared with all other TLR4 genotypes. Furthermore, the population frequency of this genotype, at a two‐fold excess over Hardy–Weinberg expectation, was increased by nonrandom mating. Near‐complete sampling and full pedigree and reproductive data enabled us to eliminate other potential causes of these patterns including inbreeding, year effects, density dependence, selection on animals at earlier life history stages or genome‐level association of the TLR4_E allele with ‘good genes’. However, comparison of observed levels of gene diversity to predictions under simulated genetic drift revealed results consistent with neutral expectations for all loci, including TLR4. Although selection favoured TLR4_BE heterozygotes in this population, these effects were insufficient to outweigh genetic drift. This is the first empirical study to show that genetic drift can overwhelm natural selection in a wild population immediately following establishment.     

Isolation and characterization of microsatellite loci from the apple maggot fly Rhagoletis pomonella (Diptera: Tephritidae)

We report the isolation and development of 81 novel primers for amplifying microsatellite loci in the Rhagoletis pomonella sibling species complex, and the sequencing, characterization and analysis of basic population genetic parameters for nine of these genes. We also report the successful cross‐species amplification of several of these loci. The R. pomonella sibling species complex is a textbook example of genetic differentiation in sympatry via host‐plant shifting. Microsatellite markers can be useful for mapping host‐plant‐associated adaptations in Rhagoletis that generate reproductive isolation and facilitate speciation, as well as for resolving the genetic structure and evolutionary history of fly populations.     

Development and characterization of novel microsatellite markers for the common earthworm (Lumbricus terrestris L.)

We developed and characterized 10 highly polymorphic microsatellite loci from an SSR‐enriched genomic DNA library of the common earthworm (Lumbricus terrestris L). Characterization of these loci using 32 individuals revealed high levels of genetic diversity, five to 18 alleles per locus and a high observed and expected heterozygosity. These loci will be used for paternity analysis and population genetic studies of the co‐evolution between L. terrestris and its parasites.     

Variable nuclear markers for Melanoplus oregonensis identified from the screening of a genomic library

We report the isolation of seven single‐copy nuclear loci from the montane grasshopper Melanoplus oregonensis. With an average length of approximately 930 bp and between five and 67 variable sites, these loci are useful for population genetic analyses within M. oregonensis. Amplification of nuclear loci in additional species of western Melanoplus grasshoppers suggests that they will also be useful for a variety of population genetic, phylogeographic and phylogenetic studies.     

Microsatellite multiplexes for high‐throughput genotyping of French grunts (Haemulon flavolineatum,Pisces: Haemulidae) and their utility in other grunt species

Grunts of the genus Haemulon are abundant in the western Atlantic and are an important component of regional fisheries. Management efforts for these and other Caribbean reef fish require an understanding of population structure and the extent of dispersal among populations. We characterized 11 polymorphic microsatellite loci for French grunts (Haemulon flavolineatum) and grouped them into three PCR multiplexes to facilitate high‐throughput genotyping. These loci will be useful for population genetic studies of French grunts. Cross‐species amplifications suggest these loci are also suitable for population genetic studies of H. aurolineatum, H. plumieri, and H. sciurus.     

Isolation and characterization of Ligustrum micranthum (Oleaceae) microsatellite loci using paired‐end Illumina reads

Twenty‐six microsatellite loci were developed and characterized for Ligustrum micranthum, a species endemic to the Ogasawara Islands, Japan. The genetic structure of this species must be clarified in order to restore the island's ecosystem. A total of 8511 primer pairs were designed from de novo sequencing. Of the 48 primer pairs selected, amplification and polymorphisms were tested using one population each from the Chichijima and Hahajima Islands of the Ogasawara Islands. Twenty‐six microsatellite loci were successfully amplified and the number of alleles for these loci ranged from five to 31 per locus, and the mean expected heterozygosities were 0.858 and 0.849, respectively. No significant deviation from the Hardy–Weinberg equilibrium was observed in either population, and no significant linkage disequilibrium was detected between any locus pair. The microsatellite loci reported in this study can be used in future studies to evaluate the genetic structure and mating system of L. micranthum.     

Characterization of comparative genome‐derived simple sequence repeats for acanthopterygian fishes

Simple sequence repeats (SSRs) have become one of the most popular molecular markers for population genetic studies. The application of SSR markers has often been limited to source species because SSR loci are too labile to be maintained in even closely related species. However, a few extremely conserved SSR loci have been reported. Here, we tested for the presence of conserved SSR loci in acanthopterygian fishes, which include over 14 000 species, by comparing the genome sequences of four acanthopterygian fishes. We also examined the comparative genome‐derived SSRs (CG‐SSRs) for their transferability across acanthopterygian fishes and their applicability to population genetic analysis. Forty‐six SSR loci with conserved flanking regions were detected and examined for their transferability among seven nonacanthopterygian and 27 acanthopterygian fishes. The PCR amplification success rate in nonacanthopterygian fishes was low, ranging from 2.2% to 21.7%, except for Lophius litulon (Lophiiformes; 80.4%). Conversely, the rate in most acanthopterygian fishes exceeded 70.0%. Sequencing of these 46 loci revealed the presence of SSRs suitable for scoring while fragment analysis of 20 loci revealed polymorphisms in most of the acanthopterygian fishes. Population genetic analysis of Cottus pollux (Scorpaeniformes) and Sphaeramia orbicularis (Perciformes) using CG‐SSRs showed that these populations did not deviate from linkage equilibrium or Hardy–Weinberg equilibrium. Furthermore, almost no loci showed evidence of null alleles, suggesting that CG‐SSRs have strong resolving power for population genetic analysis. Our findings will facilitate the use of these markers in species in which markers remain to be identified.     

Genetic variation and structure in the Mediterranean shrubs Myrtus communis and Pistacia lentiscus in different landscape contexts

Studies concerning different habitat configurations can provide insights into the complex interactions between species’ life‐history traits and the environment and can help to predict patterns in population genetics. In this study, we compared patterns of genetic variation in two Mediterranean shrub species (Myrtus communis and Pistacia lentiscus) that co‐occur in populations within three contrasting landscape contexts: continuous, fragmented‐connected and fragmented‐isolated populations. Analysing variation at microsatellites loci, our results revealed weak responses to the landscape contexts. We rather found a population‐specific response in both study species. However, despite both study species sharing similar levels of genetic diversity, Myrtus displayed higher levels of homozygosity and genetic differentiation among populations, stronger patterns of within‐population spatial genetic structure, lower values of mutation‐scaled effective population size and stronger evidence for recent genetic bottlenecks than Pistacia. This result highlights the influence of past events (e.g. historical connectivity, fluctuations in population size) and local factors (e.g. microhabitat availability for recruitment, habitat quality, plant density, native fauna) and that the landscape configuration per se (i.e. fragment size and/or isolation) might not completely determine the species’ genetic patterns.     

Does Hyphal‐Tip Ensure the Same Allelic Composition at SSR Loci as Monosporic Isolates of Sclerotinia sclerotiorum?

The proper characterization of individual is a basic stage in population genetic studies. In Sclerotinia sclerotiorum, genetic uniformity of an individual can be obtained by isolation of single ascospore; however, hyphal‐tip isolates are commonly used in genetic studies. The aim of this study was to assess whether hyphal‐tip isolates of S. sclerotiorum can be used as surrogate of monoascosporic (monosporic) isolates. Twenty‐eight isolates of S. sclerotiorum were collected from common bean plants with white mold symptoms and were purified by hyphal‐tip or single ascospore. The correspondence between hyphal‐tip and monosporic isolates was assessed through the allelic composition at 10 microsatellite (SSR) loci of the isolates obtained by both methods. For the SSR loci comprised of dinucleotide repeats in 92% of the cases, the difference (di) between the amplicon size values for hyphal‐tip and monosporic isolates was no more than one base pair. For the loci comprised of tetra or pentanucleotide repeats in 89% of the cases, di was no more than one base pair. The same allelic profile was found for hyphal‐tip or single ascospore isolates of S. sclerotiorum. When monosporic isolates cannot be easily obtained, hyphal‐tip can safeguard the genotypic identity of S. sclerotiorum isolates.     

Isolation and characterization of microsatellite markers in Dendrobium officinale,an endangered herb endemic to China

As an invaluable herb, Dendrobium officinale has been in severe danger since 1950 because of human exploitation and habitat deterioration. In order to investigate the genetic diversity and structure of this species, we isolated and characterized 12 microsatellite loci derived from two microsatellite‐enriched libraries. Twenty‐two individuals from Leye population were analysed. These loci were polymorphic and displayed three to 12 alleles per locus. The observed and expected heterozygosities ranged from 0.150 to 0.624 and from 0.162 to 0.605, respectively. These are the first microsatellite loci characterized from D. officinale that will contribute to research on the conservation, genetic diversity, population structure and individual authentication.     

Population genomics of rapidly invading lionfish in the Caribbean reveals signals of range expansion in the absence of spatial population structure

Range expansions driven by global change and species invasions may have significant genomic, evolutionary, and ecological implications. During range expansions, strong genetic drift characterized by repeated founder events can result in decreased genetic diversity with increased distance from the center of the historic range, or the point of invasion. The invasion of the Indo‐Pacific lionfish, Pterois volitans, into waters off the US East Coast, Gulf of Mexico, and Caribbean Sea provides a natural system to study rapid range expansion in an invasive marine fish with high dispersal capabilities. We report results from 12,759 single nucleotide polymorphism loci sequenced by restriction enzyme‐associated DNA sequencing for nine P. volitans sampling areas in the invaded range, including Florida and other sites throughout the Caribbean, as well as mitochondrial control region D‐loop data. Analyses revealed low to no spatially explicit metapopulation genetic structure, which is partly consistent with previous finding of little structure within ocean basins, but partly divergent from initial reports of between‐basin structure. Genetic diversity, however, was not homogeneous across all sampled sites. Patterns of genetic diversity correlate with invasion pathway. Observed heterozygosity, averaged across all loci within a population, decreases with distance from Florida while expected heterozygosity is mostly constant in sampled populations, indicating population genetic disequilibrium correlated with distance from the point of invasion. Using an F_ST outlier analysis and a Bayesian environmental correlation analysis, we identified 256 and 616 loci, respectively, that could be experiencing selection or genetic drift. Of these, 24 loci were shared between the two methods.     

Polymorphic microsatellite markers for studies of Aedes aegypti (Diptera: Culicidae), the vector of dengue and yellow fever

A significant challenge to population genetic studies of the dengue vector, Aedes aegypti, has been the lack of polymorphic microsatellite loci. In an effort to develop useful markers, we evaluated the genetic variation at 17 microsatellite loci identified in the A. aegypti genome. Nine loci with at least five alleles were identified in field‐collected specimens from Thailand. An additional two loci carried five alleles if samples from an A. aegypti laboratory colony were included. Our results greatly increase the number of highly variable markers available for the study of the genetics and the population structure of this medically important species.     

Polymorphic microsatellite loci in the bush‐cricket Metrioptera roeseli

Metrioptera roeseli is a common bush‐cricket in Europe which prefers grassland as habitat. We developed microsatellite loci in order to study the population genetic structure and the effects of different habitat types on the genetic structure of the bush‐cricket in a rural landscape. In this paper, we report on six polymorphic microsatellite loci with seven to 16 alleles per locus in 20 individuals from the Lahn‐Dill‐Bergland in Hesse, Germany.     

More precisely biased: increasing the number of markers is not a silver bullet in genetic bottleneck testing

In response to our review of the use of genetic bottleneck tests in the conservation literature (Peery et al. 2012, Molecular Ecology, 21 , 3403–3418), Hoban et al. (2013, Molecular Ecology, in press) conducted population genetic simulations to show that the statistical power of genetic bottleneck tests can be increased substantially by sampling large numbers of microsatellite loci, as they suggest is now possible in the age of genomics. While we agree with Hoban and co‐workers in principle, sampling large numbers of microsatellite loci can dramatically increase the probability of committing type 1 errors (i.e. detecting a bottleneck in a stable population) when the mutation model is incorrectly assumed. Using conservative values for mutation model parameters can reduce the probability of committing type 1 errors, but doing so can result in significant losses in statistical power. Moreover, we believe that practical limitations associated with developing large numbers of high‐quality microsatellite loci continue to constrain sample sizes, a belief supported by a literature review of recent studies using next generation sequencing methods to develop microsatellite libraries. conclusion, we maintain that researchers employing genetic bottleneck tests should proceed with caution and carefully assess both statistical power and type 1 error rates associated with their study design.     

Accounting for missing data in the estimation of contemporary genetic effective population size (Ne)

Theoretical models are often applied to population genetic data sets without fully considering the effect of missing data. Researchers can deal with missing data by removing individuals that have failed to yield genotypes and/or by removing loci that have failed to yield allelic determinations, but despite their best efforts, most data sets still contain some missing data. As a consequence, realized sample size differs among loci, and this poses a problem for unbiased methods that must explicitly account for random sampling error. One commonly used solution for the calculation of contemporary effective population size (N_e) is to calculate the effective sample size as an unweighted mean or harmonic mean across loci. This is not ideal because it fails to account for the fact that loci with different numbers of alleles have different information content. Here we consider this problem for genetic estimators of contemporary effective population size (N_e). To evaluate bias and precision of several statistical approaches for dealing with missing data, we simulated populations with known N_e and various degrees of missing data. Across all scenarios, one method of correcting for missing data (fixed‐inverse variance‐weighted harmonic mean) consistently performed the best for both single‐sample and two‐sample (temporal) methods of estimating N_e and outperformed some methods currently in widespread use. The approach adopted here may be a starting point to adjust other population genetics methods that include per‐locus sample size components.     

Drifting fronds and drifting alleles: range dynamics,local dispersal and habitat isolation shape the population structure of the estuarine seaweed Fucus ceranoides

Aim The seaweed Fucus ceranoides is restricted to spatially discrete estuarine habitats and lacks planktonic dispersal phases; it is therefore expected to exhibit strong population differentiation. Its cold‐temperate affinities and mtDNA variation imply that the northern part of the species’ range, where F. ceranoides is now ubiquitous, was recently colonized after the onset of the last deglaciation, potentially resulting in areas with greater genetic homogeneity. Here we examine the population structure of F. ceranoides to test these predictions, emphasizing the contrasting genetic signatures of limited dispersal in refugial versus recently colonized regions. Location North‐eastern Atlantic estuaries from Portugal to Norway. Methods A total of 504 individuals from 21 estuarine sites spanning the entire range of F. ceranoides were sampled and genotyped for nine microsatellite loci. Population structure was inferred from several genotypic and allele‐frequency analyses. Geographical patterns of genetic diversity were used to reconstruct the historical biogeography of the species. Results Genetic diversity and regional population differentiation showed a consistent decline with increasing latitude. Southernmost populations harboured most of the endemic variation, whereas the northern populations (> 55° N) were almost fixed for the same alleles across loci. In southern and central regions of its distribution, F. ceranoides showed striking population subdivision, with many of the sampled estuaries corresponding to coherent genetic units that were easily discriminated from one another with standard clustering methods. Main conclusions The geographical pattern of genetic diversity supports the long‐term refugial status of Iberia and a post‐glacial range expansion of F. ceranoides into previously glaciated latitudes. Despite the species’ capacity to colonize newly available habitats, the genetic structure of F. ceranoides outside the recently (re)colonized range reveals that gene flow between populations is extremely low. This study provides a remarkable example of how infrequent and spatially limited dispersal can have contrasting effects at the scales of meta‐population (connectivity) versus range dynamics (habitat tracking), and of how dispersal restrictions can result in either genetic divergence or homogeneity depending on the maturity and demographic conditions of the populations.     

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site‐associated DNA sequencing (RAD‐Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F_ST 0.008; range 0.00–0.018) were concordant with previous microsatellite‐based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age‐1 and age‐2 families of full and half‐siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non‐native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.     

Identification of microsatellite loci in Collinsia verna (Veronicaceae)

We developed eight polymorphic microsatellite loci for Collinsia verna (Veronicaceae). In a sample of 18–35 individuals from a single population, we found two to 15 alleles per locus (mean 8.3). We also tested these loci for cross‐amplification in all 22 species in the tribe Collinseae. Overall, more than half the species in the tribe amplified one microsatellite while three species most closely related to C. verna (Collinsia violacea, Collinsia parviflora and Collinsia grandiflora) amplified multiple microsatellite loci. These microsatellite loci will be used in future studies of mating system in this tribe and other quantitative genetic and population genetic studies.