Why sampling scheme matters: the effect of sampling scheme on landscape genetic results

There has been a recent trend in genetic studies of wild populations where researchers have changed their sampling schemes from sampling pre-defined populations to sampling individuals uniformly across landscapes. This reflects the fact that many species under study are continuously distributed rather than clumped into obvious “populations”. Once individual samples are collected, many landscape genetic studies use clustering algorithms and multilocus genetic data to group samples into subpopulations. After clusters are derived, landscape features that may be acting as barriers are examined and described. In theory, if populations were evenly sampled, this course of action should reliably identify population structure. However, genetic gradients and irregularly collected samples may impact the composition and location of clusters. We built genetic models where individual genotypes were either randomly distributed across a landscape or contained gradients created by neighbor mating for multiple generations. We investigated the influence of six different sampling protocols on population clustering using program STRUCTURE, the most commonly used model-based clustering method for multilocus genotype data. For models where individuals (and their alleles) were randomly distributed across a landscape, STRUCTURE correctly predicted that only one population was being sampled. However, when gradients created by neighbor mating existed, STRUCTURE detected multiple, but different numbers of clusters, depending on sampling protocols. We recommend testing for fine scale autocorrelation patterns prior to sample clustering, as the scale of the autocorrelation appears to influence the results. Further, we recommend that researchers pay attention to the impacts that sampling may have on subsequent population and landscape genetic results. The U.S. Government's right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

Source verification of mis-identified Arabidopsis thaliana accessions

A major strength of Arabidopsis thaliana as a model lies in the availability of a large number of naturally occurring inbred lines. Recent studies of A. thaliana population structure, using thousands of accessions from stock center and natural collections, have revealed a robust pattern of isolation by distance at several spatial scales, such that genetically identical individuals are generally found close to each other. However, some individual accessions deviate from this pattern. While some of these may be the products of rare long‐distance dispersal events, many deviations may be the result of mis‐identification, in the sense that the data regarding location of origin data are incorrect. Here, we aim to identify such discrepancies. Of the 5965 accessions examined, we conclude that 286 deserve special attention as being potentially mis‐identified. We describe these suspicious accessions and their possible origins, and advise caution with regard to their use in experiments in which accurate information on geographic origin is important. Finally, we discuss possibilities for maintaining the integrity of stock lines.     

Natural genetic variation of Arabidopsis thaliana is geographically structured in the Iberian peninsula

To understand the demographic history of Arabidopsis thaliana within its native geographical range, we have studied its genetic structure in the Iberian Peninsula region. We have analyzed the amount and spatial distribution of A. thaliana genetic variation by genotyping 268 individuals sampled in 100 natural populations from the Iberian Peninsula. Analyses of 175 individuals from 7 of these populations, with 20 chloroplast and nuclear microsatellite loci and 109 common single nucleotide polymorphisms, show significant population differentiation and isolation by distance. In addition, analyses of one genotype from 100 populations detected significant isolation by distance over the entire Iberian Peninsula, as well as among six Iberian subregions. Analyses of these 100 genotypes with different model-based clustering algorithms inferred four genetic clusters, which show a clear-cut geographical differentiation pattern. On the other hand, clustering analysis of a worldwide sample showed a west–east Eurasian longitudinal spatial gradient of the commonest Iberian genetic cluster. These results indicate that A. thaliana genetic variation displays significant regional structure and consistently support the hypothesis that Iberia has been a glacial refugium for A. thaliana. Furthermore, the Iberian geographical structure indicates a complex regional population dynamics, suggesting that this region contained multiple Pleistocene refugia with a different contribution to the postglacial colonization of Europe.     

Development of microsatellite markers for the critically endangered Limonium dufourii (Girard) Kuntze (Plumbaginaceae)

Limonium dufourii is an endemic plant from the eastern Mediterranean coast of Spain with a triploid chromosome number and apomictic reproduction. We have isolated and characterized 13 polymorphic microsatellite loci from an enriched library in order to investigate its population genetic structure. Simple sequence repeat (SSR) loci were screened in 120 individuals from the six extant populations of this species. They show an average of 5.76 alleles per locus, ranging from 2 to 18, with seven loci exhibiting heterozygosities larger than 0.60. Three loci present one single allele in each individual, whereas one locus presents three alleles in every individual analysed.     

Effects of sample size,number of markers,and allelic richness on the detection of spatial genetic pattern

The influence of study design on the ability to detect the effects of landscape pattern on gene flow is one of the most pressing methodological gaps in landscape genetic research. To investigate the effect of study design on landscape genetics inference, we used a spatially‐explicit, individual‐based program to simulate gene flow in a spatially continuous population inhabiting a landscape with gradual spatial changes in resistance to movement. We simulated a wide range of combinations of number of loci, number of alleles per locus and number of individuals sampled from the population. We assessed how these three aspects of study design influenced the statistical power to successfully identify the generating process among competing hypotheses of isolation‐by‐distance, isolation‐by‐barrier, and isolation‐by‐landscape resistance using a causal modelling approach with partial Mantel tests. We modelled the statistical power to identify the generating process as a response surface for equilibrium and non‐equilibrium conditions after introduction of isolation‐by‐landscape resistance. All three variables (loci, alleles and sampled individuals) affect the power of causal modelling, but to different degrees. Stronger partial Mantel r correlations between landscape distances and genetic distances were found when more loci were used and when loci were more variable, which makes comparisons of effect size between studies difficult. Number of individuals did not affect the accuracy through mean equilibrium partial Mantel r, but larger samples decreased the uncertainty (increasing the precision) of equilibrium partial Mantel r estimates. We conclude that amplifying more (and more variable) loci is likely to increase the power of landscape genetic inferences more than increasing number of individuals.     

Planning and Assessment of Activity Budget Studies Employing Instantaneous Sampling

To obtain accurate estimates of activity budget parameters, samples must be unbiased and precise. Many researchers have considered how biased data may affect their ability to draw conclusions and examined ways to decrease bias in sampling efforts, but few have addressed the implications of not considering estimate precision. We propose a method to assess whether the number of instantaneous samples collected is sufficient to obtain precise activity budget parameter estimates. We draw on sampling theory to determine the number of observations per animal required to reach a desired bound on the error of estimation based on a stratified random sample, with individual animals acting as strata. We also discuss the optimal balance between the number of individuals sampled and the number of observations sampled per individual for a variety of sampling conditions. We present an empirical dataset on pronghorn (Antilocapra americana) as an example of the utility of the method. The required numbers of observation to reach precise estimates for pronghorn varied between common and rare behaviors, but precise estimates were achieved with <255 observations per individual for common behaviors. The two most apparent factors affecting the required number of observations for precise estimates were the number of individuals sampled and the complexity of the activity budget. This technique takes into account variation associated with individual activity budgets and population variation in activity budget parameter estimates, and helps to ensure that estimates are precise. The method can also be used for planning future sampling efforts.     

Genetic basis of adaptation: flowering time in Arabidopsis thaliana

 We have mapped QTLs (quantitative trait loci) for an adaptive trait, flowering time, in a selfing annual, Arabidopsis thaliana. To obtain a mapping population we made a cross between an early-summer, annual strain, Li-5, and an individual from a late over-wintering natural population, Naantali. From the backcross to Li-5 298 progeny were grown, of which 93 of the most extreme individuals were genotyped. The data were analysed with both interval mapping and composite interval mapping methods to reveal one major and six minor QTLs, with at least one QTL on each of the five chromosomes. The QTL on chromosome 4 was a major one with an effect of 17.3 days on flowering time and explaining 53.4% of the total variance. The others had effects of at most 6.5 days, and they accounted for only small portions of the variance. Epistasis was indicated between one pair of the QTLs. The result of finding one major QTL and little epistasis agrees with previous studies on flowering time in Arabidopsis thaliana and other species. That several QTLs were found was expected considering the large number of possible candidate loci. In the light of the suggested genetic models of gene action at the candidate loci, epistasis was to be expected. The data showed that major QTLs for adaptive traits can be detected in non-domesticated species. Received: 15 January 1997/Accepted: 21 February 1997     

Assessing polar bear (Ursus maritimus) population structure in the Hudson Bay region using SNPs

Defining subpopulations using genetics has traditionally used data from microsatellite markers to investigate population structure; however, single‐nucleotide polymorphisms (SNPs) have emerged as a tool for detection of fine‐scale structure. In Hudson Bay, Canada, three polar bear (Ursus maritimus) subpopulations (Foxe Basin (FB), Southern Hudson Bay (SH), and Western Hudson Bay (WH)) have been delineated based on mark–recapture studies, radiotelemetry and satellite telemetry, return of marked animals in the subsistence harvest, and population genetics using microsatellites. We used SNPs to detect fine‐scale population structure in polar bears from the Hudson Bay region and compared our results to the current designations using 414 individuals genotyped at 2,603 SNPs. Analyses based on discriminant analysis of principal components (DAPC) and STRUCTURE support the presence of four genetic clusters: (i) Western—including individuals sampled in WH, SH (excluding Akimiski Island in James Bay), and southern FB (south of Southampton Island); (ii) Northern—individuals sampled in northern FB (Baffin Island) and Davis Strait (DS) (Labrador coast); (iii) Southeast—individuals from SH (Akimiski Island in James Bay); and (iv) Northeast—individuals from DS (Baffin Island). Population structure differed from microsatellite studies and current management designations demonstrating the value of using SNPs for fine‐scale population delineation in polar bears.     

Isolation and characterization of 12 tetra‐repeated microsatellite loci from the Formosan grass lizard (Takydromus formosanus)

Twelve tetra‐repeated microsatellite loci were developed for the island‐endemic grass lizard Takydromus formosanus (Squamata: Lacertidae). We characterized these loci by genotyping 51 individuals sampled from six localities. The number of alleles per locus ranged between 12 and 23, whereas the observed heterozygosity ranged between 0.686 and 0.941. None of these loci showed gametic disequilibrium or sex linkage. With such high polymorphism, we believe that these loci should be suitable for fine‐scaled analysis of population structure, hybridizations or individual identification of grass lizards.     

Evidence for multiple introductions of Centaurea stoebe micranthos (spotted knapweed,Asteraceae) to North America

Invasive species’ success may depend strongly on the genetic resources they maintain through the invasion process. We ask how many introductions have occurred in the North American weed Centaurea stoebe micranthos (Asteraceae), and explore whether genetic diversity and population structure have changed as a result of introduction. We surveyed individuals from 15 European native range sites and 11 North American introduced range sites at six polymorphic microsatellite loci. No significant difference existed in the total number of alleles or in the number of private alleles found in each range. Shannon–Weaver diversity of phenotype frequencies was also not significantly different between the ranges, while expected heterozygosity was significantly higher in the invasive range. Population structure was similar between the native range and the invasive range, and isolation by distance was not significant in either range. Traditional assignment methods did not allocate any North American individuals to the sampled European populations, while Bayesian assignment methods grouped individuals into nine genetic clusters, with three of them shared between North America and Europe. Invasive individuals tended to have genetically admixed profiles, while natives tended to assign more strongly to a single cluster. Many North American individuals share assignment with Romania and Bulgaria, suggesting two separate invasions that have undergone gene flow in North America. Samples from three other invasive range sites were genetically distinct, possibly representing three other unique introductions. Multiple introductions and the maintenance of high genetic diversity through the introduction process may be partially responsible for the invasive success of C. stoebe micranthos.     

THE EVOLVING GENETIC HISTORY OF A POPULATION OF LATHYRUS SYLVESTRIS: EVIDENCE FROM TEMPORAL AND SPATIAL GENETIC STRUCTURE

We analyze patterns of genetic microdifferentiation within a natural population of Lathyrus sylvestris, a perennial herb with both sexual reproduction and clonal growth. In a population from the northern foothills of the Pyrénées in southwestern France, a combined demographic and genetic investigation enabled the study not only of spatial genetic structure of the population, but also of the history of the population's spatial genetic structure over time. Excavation of all individuals allowed identification of clonemates. Age of each individual was determined by counting annual growth rings in the taproot, a method tested with individuals of known age planted in experimental gardens. Each individual was mapped, and genotypes of all individuals were determined using allozyme markers for a number of polymorphic loci. Distribution patterns and spatial genetic structure, both for all individuals and for different age classes, were analyzed using spatial autocorrelation statistics (Geary's Index, Moran's Index). Patterns of gene flow within the population were also studied using F-statistics and tests for random associations of alleles. Because age, allele frequencies, and location were known for each individual, it was possible to study how spatial genetic structure changed over time. Results from all these diverse approaches are consistent with one another, and clearly show the following: (1) founder effects, with the study transect being first colonized by individuals at either end of the transect that were homozygous for different alleles at one marker locus; (2) a difference in spatial distribution of individuals originated from sexual reproduction (seedlings) and from clonal growth (connected individuals); (3) restricted gene flow, due to inbreeding among related, clumped individuals; and (4) increase in heterozygote deficit within the youngest cohort of individuals. The results indicate that genetic differentiation in time was much less marked than differentiation in space. Nevertheless, the results revealed that the studied population is experiencing demographic and genetic variation in time, suggesting that it is not at equilibrium. On the one hand, spatial structuring is becoming less marked due to the recombination of founder genotypes; on the other hand, as establishment of new individuals increases, a new spatial structure emerges due to mating between relatives.     

Identifying the source of species invasions: sampling intensity vs. genetic diversity

Population geneticists and community ecologists have long recognized the importance of sampling design for uncovering patterns of diversity within and among populations and in communities. Invasion ecologists increasingly have utilized phylogeographical patterns of mitochondrial or chloroplast DNA sequence variation to link introduced populations with putative source populations. However, many studies have ignored lessons from population genetics and community ecology and are vulnerable to sampling errors owing to insufficient field collections. A review of published invasion studies that utilized mitochondrial or chloroplast DNA markers reveals that insufficient sampling could strongly influence results and interpretations. Sixty per cent of studies sampled an average of less than six individuals per source population, vs. only 45% for introduced populations. Typically, far fewer introduced than source populations were surveyed, although they were sampled more intensively. Simulations based on published data forming a comprehensive mtDNA haplotype data set highlight and quantify the impact of the number of individuals surveyed per source population and number of putative source populations surveyed for accurate assignment of introduced individuals. Errors associated with sampling a low number of individuals are most acute when rare source haplotypes are dominant or fixed in the introduced population. Accuracy of assignment of introduced individuals is also directly related to the number of source populations surveyed and to the degree of genetic differentiation among them ( F _ST). Incorrect interpretations resulting from sampling errors can be avoided if sampling design is considered before field collections are made.     

Ecological and historical determinants of population genetic structure and diversity in the Mediterranean shrub Rosmarinus officinalis (Lamiaceae)

Population genetic studies of widespread Mediterranean shrubs are scarce compared with those of trees and narrow endemics or studies from phylogeographical perspectives, despite the key role these species may play in Mediterranean ecosystems. Knowledge on the effect of ecological factors in shaping their genetic patterns is also limited. In this study we investigate genetic diversity and population structure across 18 populations of Rosmarinus officinalis, a Mediterranean shrubland plant. Populations were sampled along two elevational gradients, one each on calcareous and siliceous soils in a mountain system in the eastern Iberian Peninsula, to decipher the effect of ecological factors on the genetic diversity and structure based on 11 microsatellite loci. We found overall high levels of genetic diversity and weak population structure. Genetic diversity increased with elevation, whereas population differentiation was stronger among populations growing on siliceous soils. The nested analysis of elevational gradients within soil types revealed that these general patterns were mostly driven by siliceous populations, whereas calcareous populations were more homogeneous along elevational belts. Bayesian analysis of population structure revealed genetic membership of lowland and high‐elevation populations to different genetic clusters and a higher admixture of intermediate‐elevation populations to both clusters. High‐elevation populations were less differentiated from a hypothetical ancestral cluster, suggesting the persistence of their gene pool during the Pleistocene glaciations. In contrast, lowland populations resulted from more recent divergence. We propose that life‐history and reproductive traits mostly contribute to explain the high levels of genetic diversity and weak population structure, whereas ecological and historical factors mostly contribute to the stronger differentiation of siliceous populations and a rapid expansion of R. officinalis on calcareous soils possibly mediated by human landscape transformations, © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2016, 180 , 50–63.     

Identifying the minimum number of microsatellite loci needed to assess population genetic structure: A case study in fly culturing

Small, isolated populations are constantly threatened by loss of genetic diversity due to drift. Such situations are found, for instance, in laboratory culturing. In guarding against diversity loss, monitoring of potential changes in population structure is paramount; this monitoring is most often achieved using microsatellite markers, which can be costly in terms of time and money when many loci are scored in large numbers of individuals. Here, we present a case study reducing the number of microsatellites to the minimum necessary to correctly detect the population structure of two Drosophila nigrosparsa populations. The number of loci was gradually reduced from 11 to 1, using the Allelic Richness (AR) and Private Allelic Richness (PAR) as criteria for locus removal. The effect of each reduction step was evaluated by the number of genetic clusters detectable from the data and by the allocation of individuals to the clusters; in the latter, excluding ambiguous individuals was tested to reduce the rate of incorrect assignments. We demonstrate that more than 95% of the individuals can still be correctly assigned when using eight loci and that the major population structure is still visible when using two highly polymorphic loci. The differences between sorting the loci by AR and PAR were negligible. The method presented here will most efficiently reduce genotyping costs when small sets of loci (“core sets”) for long-time use in large-scale population screenings are compiled.     

Spatial structuring of the population genetics of a European subterranean termite species

In population genetics studies, detecting and quantifying the distribution of genetic variation can help elucidate ecological and evolutionary processes. In social insects, the distribution of population‐level genetic variability is generally linked to colony‐level genetic structure. It is thus especially crucial to conduct complementary analyses on such organisms to examine how spatial and social constraints interact to shape patterns of intraspecific diversity. In this study, we sequenced the mitochondrial COII gene for 52 colonies of the subterranean termite Reticulitermes grassei (Isoptera: Rhinotermitidae), sampled from a population in southwestern France. Three haplotypes were detected, one of which was found exclusively in the southern part of the study area (near the Pyrenees). After genotyping 6 microsatellite loci for 512 individual termites, we detected a significant degree of isolation by distance among individuals over the entire range; however, the cline of genetic differentiation was not continuous, suggesting the existence of differentiated populations. A spatial principal component analysis based on allele frequency data revealed significant spatial autocorrelation among genotypes: the northern and southern groups were strongly differentiated. This finding was corroborated by clustering analyses; depending on the randomized data set, two or three clusters, exhibiting significant degrees of differentiation, were identified. An examination of colony breeding systems showed that colonies containing related neotenic reproductives were prevalent, suggesting that inbreeding may contribute to the high level of homozygosity observed and thus enhance genetic contrasts among colonies. We discuss the effect of evolutionary and environmental factors as well as reproductive and dispersal modes on population genetic structure.     

The influence of trap type on evaluating population structure of the semifossorial and social rodent <Emphasis Type="Italic">Octodon degus</Emphasis>

Trap type may influence captures of individuals in different age-sex categories in small mammal studies, resulting in biased population and demographic information. We deployed 4 live trap types at burrow systems of the rodent, Octodon degus Molina, 1782, in central Chile to determine trap efficacy in capturing individuals of 6 demographic categories. We captured 2672 individuals in 17 709 trap days (15.1% trapping success). Tomahawks were the most efficient trap capturing half of individuals during both years, followed by mesh Sherman traps, large Sherman traps, and medium Sherman traps in 2005. All trap types equally sampled sexes. Large and medium Sherman traps provided similar demographic structure, where half of the individuals captured were pups; Tomahawk traps sampled more adults than pups. Relative captures of pups were similar across different trap types, suggesting that pups are equally sampled by each of the deployed trap types. Relative captures of adults were lower in Sherman traps, suggesting that this age class avoided solid-walled traps. For Octodon degus, the sole use of Tomahawk traps may produce sufficient, unbiased demographic data. Only 4 trap mortalities occurred (0.15%). Researchers may minimize trap mortality without compromising sufficient demographic sampling by trapping during peak animal activity.     

Paternity analysis in Ophiopogon xylorrhizus Wang et Tai (Liliaceae s.l.): selfing assures reproductive success

We assessed the utility of codominant allozyme markers for paternity analysis in natural populations of Ophiopogon xylorrhizus (Liliaceae s. l.) by means of likelihood‐based approach. In the three independent natural stands of 46 individual plants with 354 seeds of O. xylorrhizus (Liliaceae s. l), we assigned 319 seeds to a single individual with the highest LOD score. Among these, 179 seeds had a significant D‐value. An individual (acting either as female or male) mateds with an average of one other individual on average to reproduce, except when it self‐pollinated. If an individual that neither donated pollen to other individuals nor selfed, being strictly a pollen acceptor, the number of offspring produced was half that produced by as many offspring as individuals which both provided and accepted pollen (4.1 vs. 9.5). Significant positive relationships were found in all the three stands between the number of self‐pollen an individual produced and the number of seeds that individual set. Male outcrossing reproductive success was unevenly varied, both in amount and spatial location. The individuals located in the centre of the stand and the individuals closest to each other had higher male outcrossing reproductive success than those at the edge of the stand. Male selfing reproductive success showed a similar pattern, with an even distribution of female reproductive output. Most mating events were took place within 20 m. The distance between mates mirrored the active ranges of pollinators, and shaped the fine‐scale population spatial genetic spatial structure. The results indicate that selfing assures reproductive success in O. xylorrhizus.     

Measuring mast seeding behavior: relationships among population variation,individual variation and synchrony

Mast seeding, or masting, is the variable production of flowers, seeds, or fruit across years more or less synchronously by individuals within a population. A critical issue is the extent to which temporal variation in seed production over a collection of individuals can be viewed as arising from a combination of individual variation and synchrony among individuals. Studies of masting typically quantify such variation in terms of the coefficient of variation (CV). In this paper we examine mathematically how the population CV relates to the mean individual CV and synchrony, concluding that the relationship is a complex one which cannot isolate an overall measure of synchrony, and involves additional factors, principally the number of plants sampled and the mean productivity per plant. Our development suggests some simple approximate relationships of population CV to individual variability, synchrony and the number of individuals. These were found to fit quite well when applied to data from 59 studies which included seed production at the individual level.     

The biological control of Pomacea canaliculata population by rice-duck mutualism in paddy fields

Duck has been used as a non-chemical control method against Pomacea canaliculata Lamarck, but little is known about its principles that underlie the control of snail populations. An indoor experiment was initially used to observe the predation potential of ducks, followed by replicated field trials. In the indoor studies, ducks effectively preyed on juvenile snails, but had a weak predatory effect on large snails and egg clusters. In the field, application of a rice-duck mutualism system significantly reduced the numbers of snails (especially number of immature individuals), number of snail egg clusters and snail damage to rice plants. The controlling effect was longer and more stable than the chemical application, resulting in a better yield than with the pentachlorophenol sodium and tea seed powder treatment. Our experimental results also suggested that the snail age structure in the rice-duck mutualism plots was shifted towards older snails by ducks preying, indicating a trend towards population decline, and ducks caused snails to oviposit on sites not ideal for hatchling establishment. Throughout the studies, it is suggested that a rice-duck mutualism system could be used for controlling P. canaliculata in organic rice production.