Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species

High‐throughput DNA sequencing facilitates the analysis of large portions of the genome in nonmodel organisms, ensuring high accuracy of population genetic parameters. However, empirical studies evaluating the appropriate sample size for these kinds of studies are still scarce. In this study, we use double‐digest restriction‐associated DNA sequencing (ddRADseq) to recover thousands of single nucleotide polymorphisms (SNPs) for two physically isolated populations of Amphirrhox longifolia (Violaceae), a nonmodel plant species for which no reference genome is available. We used resampling techniques to construct simulated populations with a random subset of individuals and SNPs to determine how many individuals and biallelic markers should be sampled for accurate estimates of intra‐ and interpopulation genetic diversity. We identified 3646 and 4900 polymorphic SNPs for the two populations of A. longifolia, respectively. Our simulations show that, overall, a sample size greater than eight individuals has little impact on estimates of genetic diversity within A. longifolia populations, when 1000 SNPs or higher are used. Our results also show that even at a very small sample size (i.e. two individuals), accurate estimates of F_ST can be obtained with a large number of SNPs (≥1500). These results highlight the potential of high‐throughput genomic sequencing approaches to address questions related to evolutionary biology in nonmodel organisms. Furthermore, our findings also provide insights into the optimization of sampling strategies in the era of population genomics.     

Distance-based tests for homogeneity of multivariate dispersions

Summary The traditional likelihood‐based test for differences in multivariate dispersions is known to be sensitive to nonnormality. It is also impossible to use when the number of variables exceeds the number of observations. Many biological and ecological data sets have many variables, are highly skewed, and are zero‐inflated. The traditional test and even some more robust alternatives are also unreasonable in many contexts where measures of dispersion based on a non‐Euclidean dissimilarity would be more appropriate. Distance‐based tests of homogeneity of multivariate dispersions, which can be based on any dissimilarity measure of choice, are proposed here. They rely on the rotational invariance of either the multivariate centroid or the spatial median to obtain measures of spread using principal coordinate axes. The tests are straightforward multivariate extensions of Levene's test, with P‐values obtained either using the traditional F‐distribution or using permutation of either least‐squares or LAD residuals. Examples illustrate the utility of the approach, including the analysis of stabilizing selection in sparrows, biodiversity of New Zealand fish assemblages, and the response of Indonesian reef corals to an El Niño. Monte Carlo simulations from the real data sets show that the distance‐based tests are robust and powerful for relevant alternative hypotheses of real differences in spread.     

Sample Size Effects on Estimates of Population Genetic Structure: Implications for Ecological Restoration

The field of ecological restoration is growing rapidly, and the sourcing of suitable seed is a major issue. Information on the population genetic structure of a species can provide valuable information to aid in defining seed collection zones. For a practical contribution from genetics, a rapid approach to delineating seed collection zones using genetic markers (amplified fragment length polymorphisms [AFLPs]) has been developed. Here, we test the effects of sampling regime on the efficacy of this method. Genetic data were collected for an outcrossing seeder, Daviesia divaricata ssp. divaricata, an important species in urban bushland restoration in Perth, Western Australia. The effect of sample size and number of AFLP markers on estimates of genetic variation and population structure was examined in relation to implications for sourcing material for restoration. Three different sample sizes were used (n= 8, 15, and 30) from six urban bushland remnants. High levels of genetic diversity were observed in D. divaricata (87.4% polymorphic markers), with significant population differentiation detected among sampled populations (Θ^B= 0.1386, p < 0.001). Although sample size does not appear to affect the spatial pattern in principle co‐ordinates analysis (PCA) plots, the number of polymorphic loci increased with sample size and estimates of population subdivision (F_ST and Θ^B) and associated confidence intervals decreased with increasing sample size. We recommend using a minimum of 30 plants for sourcing seed for restoration projects.     

Using maximum entropy to predict the potential distribution of an invasive freshwater snail

相似文献   

Global patterns of population genetic differentiation in seed plants

Evaluating the factors that drive patterns of population differentiation in plants is critical for understanding several biological processes such as local adaptation and incipient speciation. Previous studies have given conflicting results regarding the significance of pollination mode, seed dispersal mode, mating system, growth form and latitudinal region in shaping patterns of genetic structure, as estimated by F_ST values, and no study to date has tested their relative importance together across a broad scale. Here, we assembled a 337‐species data set for seed plants from publications with data on F_ST from nuclear markers and species traits, including variables pertaining to the sampling scheme of each study. We used species traits, while accounting for sampling variables, to perform phylogenetic multiple regressions. Results demonstrated that F_ST values were higher for tropical, mixed‐mating, non‐woody species pollinated by small insects, indicating greater population differentiation, and lower for temperate, outcrossing trees pollinated by wind. Among the factors we tested, latitudinal region explained the largest portion of variance, followed by pollination mode, mating system and growth form, while seed dispersal mode did not significantly relate to F_ST. Our analyses provide the most robust and comprehensive evaluation to date of the main ecological factors predicted to drive population differentiation in seed plants, with important implications for understanding the basis of their genetic divergence. Our study supports previous findings showing greater population differentiation in tropical regions and is the first that we are aware of to robustly demonstrate greater population differentiation in species pollinated by small insects.     

Standardized sampling plan for Aphis gossypii based on the cotton cultivar,plant phenology and crop size

A standardized sampling plan is the starting point for developing a decision‐making system for pest control. Aphis gossypii (Hemiptera: Aphididae) is a destructive sap‐feeding pest on cotton worldwide. However, research addressing cotton cultivar, plant phenology and field size with the aim of developing a sampling plan for A. gossypii has not been done. Therefore, in this study, we developed a standardized sampling for A. gossypii as a function of these factors. To accomplish this, A. gossypii densities in four experimental cotton cultivars were sampled weekly during year one to determine the ideal aphid characteristic to sample (by individual or colony). During year one and two, A. gossypii densities were sampled weekly in the same cultivars to determine sampling unit, sampling technique and the number of samples for an A. gossypii sampling plan. Using the sample number determined, the sampling time was recorded for cotton field size of 1, 5, 10, 50, 100 and 150 ha in order to estimate the sampling cost. In cotton, the count of individuals was the best characteristic for the assessment of A. gossypii. Leaves of the most apical branches for the vegetative and reproductive cotton plant stage were the best sampling units. The best sampling technique was direct counting. The cotton cultivar did not affect the development of the sampling plan. The A. gossypii sampling plan involved the evaluation of 58 samples per zone and required 20 min (<0.35 min/sample) for the evaluation of these samples. However, the walking time between samples was the main factor responsible for the total sampling time and cost in cotton fields, and this factor strongly depends on the size of the cotton field.     

Sampling Hubbell's neutral theory of biodiversity

In the context of neutral theories of community ecology, a novel genealogy‐based framework has recently furnished an analytic extension of Ewens’ sampling multivariate abundance distribution, which also applies to a random sample from a local community. Here, instead of taking a multivariate approach, we further develop the sampling theory of Hubbell's neutral spatially implicit theory and derive simple abundance distributions for a random sample both from a local community and a metacommunity. Our result is given in terms of the average number of species with a given abundance in any randomly extracted sample. Contrary to what has been widely assumed, a random sample from a metacommunity is not fully described by the Fisher log‐series, but by a new distribution. This new sample distribution matches the log‐series expectation at high biodiversity values (θ > 1) but clearly departs from it for species‐poor metacommunities (θ < 1). Our theoretical framework should be helpful in the better assessment of diversity and testing of the neutral theory by using abundance data.     

EDENetworks: A user‐friendly software to build and analyse networks in biogeography,ecology and population genetics

The recent application of graph‐based network theory analysis to biogeography, community ecology and population genetics has created a need for user‐friendly software, which would allow a wider accessibility to and adaptation of these methods. EDENetworks aims to fill this void by providing an easy‐to‐use interface for the whole analysis pipeline of ecological and evolutionary networks starting from matrices of species distributions, genotypes, bacterial OTUs or populations characterized genetically. The user can choose between several different ecological distance metrics, such as Bray‐Curtis or Sorensen distance, or population genetic metrics such as F_ST or Goldstein distances, to turn the raw data into a distance/dissimilarity matrix. This matrix is then transformed into a network by manual or automatic thresholding based on percolation theory or by building the minimum spanning tree. The networks can be visualized along with auxiliary data and analysed with various metrics such as degree, clustering coefficient, assortativity and betweenness centrality. The statistical significance of the results can be estimated either by resampling the original biological data or by null models based on permutations of the data.     

Hylobius abietis L. feeding on the novel host Pinus brutia Ten. increases emission of volatile organic compounds

Plants respond to feeding by herbivorous insects by producing volatile organic chemicals, which mediate interactions between herbivores and plants. Yet, few studies investigated whether such plant responses to herbivory differ between historical host and novel plants. Here, we investigated whether herbivory by the pine weevil Hylobius abietis causes a release of volatile organic chemicals from a novel tree Pinus brutia and compared the relative amounts of volatiles released from herbivore's historical hosts and P. brutia. We collected volatiles emitted from P. brutia seedlings that were either subjected to feeding by H. abietis or no feeding. Our results indicated that feeding increased emission of volatile compounds, composed of monoterpenes and sesquiterpenes, and that the emission was several fold higher in the damaged seedlings than in undamaged seedlings. In particular, emission of monoterpenes and sesquiterpenes increased by 4.4‐and 10‐fold in the damaged plants, respectively. Strikingly, individual monoterpenes and sesquiterpenes showed much greater dissimilarity between damaged and undamaged seedlings. Furthermore, several minor monoterpenes showed negative relationships with the weevil gnawed area. We discussed these results with the results of previous studies focused on historical host plants of H. abietis and hypothesized the ecological relevance and importance of our results pertaining relevance to the plant–herbivory interactions.     

Sampling related individuals within ponds biases estimates of population structure in a pond‐breeding amphibian

Effective conservation and management of pond‐breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding‐site connectivity. Population‐level studies of pond‐breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within‐pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome‐wide SNPs generated using a double‐digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within‐pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.     

QST–FST comparisons with unbalanced half‐sib designs

Q_ST, a measure of quantitative genetic differentiation among populations, is an index that can suggest local adaptation if Q_ST for a trait is sufficiently larger than the mean F_ST of neutral genetic markers. A previous method by Whitlock and Guillaume derived a simulation resampling approach to statistically test for a difference between Q_ST and F_ST, but that method is limited to balanced data sets with offspring related as half‐sibs through shared fathers. We extend this approach (i) to allow for a model more suitable for some plant populations or breeding designs in which offspring are related through mothers (assuming independent fathers for each offspring; half‐sibs by dam); and (ii) by explicitly allowing for unbalanced data sets. The resulting approach is made available through the R package QstFstComp.     

Validation of a cheap and simple nondestructive method for obtaining AFLPs and DNA sequences (mitochondrial and nuclear) in amphibians

The use of nondestructive methods for obtaining DNA from amphibians (e.g. buccal swabs) allows genetic studies to be performed without affecting the survival of the studied individuals. In this study, we compared two methods of nondestructive DNA sampling, buccal swabs and interdigital membrane or toe‐clipping, in several amphibian species of different size: Rhinella spinulosa, R. atacamensis, six species of the genus Telmatobius and Pleurodema thaul. We evaluated the integrity of the DNA extracted by sequencing fragments of mitochondrial and nuclear genes and by generating amplified fragment length polymorphisms markers (AFLPs). In all cases, we obtained an adequate amount of DNA (mean range 55–298 ng/μL). We obtained identical DNA sequences from buccal swab and interdigital membrane/toe‐clip for all individuals. The differences in the coding of AFLP markers between the tissues were similar to those reported for replicas of the same type of sample in similar analyses in other species of amphibians. In conclusion, the use of buccal swabs is a trustworthy and inexpensive method to obtain DNA for mitochondrial and nuclear sequencing and AFLP analyses. Given the types of markers evaluated, buccal swabs may be used for phylogenetic, phylogeographic and population genetic studies, even in small amphibians (<33 mm).     

wisepair: a computer program for individual matching in genetic tracking studies

Individual‐based data sets tracking organisms over space and time are fundamental to answering broad questions in ecology and evolution. A ‘permanent’ genetic tag circumvents a need to invasively mark or tag animals, especially if there are little phenotypic differences among individuals. However, genetic tracking of individuals does not come without its limits; correctly matching genotypes and error rates associated with laboratory work can make it difficult to parse out matched individuals. In addition, defining a sampling design that effectively matches individuals in the wild can be a challenge for researchers. Here, we combine the two objectives of defining sampling design and reducing genotyping error through an efficient Python‐based computer‐modelling program, wisepair . We describe the methods used to develop the computer program and assess its effectiveness through three empirical data sets, with and without reference genotypes. Our results show that wisepair outperformed similar genotype matching programs using previously published from reference genotype data of diurnal poison frogs (Allobates femoralis) and without‐reference (faecal) genotype sample data sets of harbour seals (Phoca vitulina) and Eurasian otters (Lutra lutra). In addition, due to limited sampling effort in the harbour seal data, we present optimal sampling designs for future projects. wisepair allows for minimal sacrifice in the available methods as it incorporates sample rerun error data, allelic pairwise comparisons and probabilistic simulations to determine matching thresholds. Our program is the lone tool available to researchers to define parameters a priori for genetic tracking studies.     

Monitoring an Arizona Ponderosa Pine Restoration: Sampling Efficiency and Multivariate Analysis of Understory Vegetation

As monitoring plans for the restoration of Pinus ponderosa forests in the southwestern United States evolve toward examining multifactor ecosystem responses to ecological restoration, designing efficient sampling procedures for understory vegetation will become increasingly important. The objective of this study was to compare understory composition and diversity among thin/burn and control treatments in a P. ponderosa restoration, while simultaneously examining the effects of sampling design and multivariate analyses on which conclusions were based. Using multi‐response permutation procedures (MRPP), we tested the null hypothesis of no difference in understory species composition among treatments using different data matrices (e.g., frequency and cover) for two different sampling methods. Treatment differences were subtle and were detected by an intensive 50, 1‐m² subplot sampling method for all data matrices but were not detected by a less intensive point‐intercept sampling method for any matrix. Sampling methods examined in this study controlled results of multivariate analyses more than the data matrices used to summarize data generated by a sampling method. We partitioned data into plant life form and native/exotic species categories for MRPP, and this partitioning isolated plant groups most responsible for treatment differences. We also examined the effects of number of 1‐m² subplots sampled on mean‐species‐richness/m² estimates and found that estimates based on 10 subplots and based on 50 subplots were highly correlated (r = 0.99). Species–area curves indicated that the 50, 1‐m² subplot sampling method detected the common species of sites but failed to detect the majority of rare species. Additional sampling‐design studies are needed to develop single sampling designs that produce multifactor data on plant composition, diversity, and spatial patterns amenable to multivariate analyses as part of monitoring plans of vegetation responses to ecological restoration.     

Incorporating functional dissimilarities into sample‐based rarefaction curves: from taxon resampling to functional resampling

Question: Indices of functional diversity have been seen as the key for integrating information on species richness with measures that focus on those components of community composition related to ecosystem functioning. For comparing species richness among habitats on an equal‐effort basis, so‐called sample‐based rarefaction curves may be used. Given a study area that is sampled for species presence and absence in N plots, sample‐based rarefaction generates the expected number of accumulated species as the number of sampled plots increases from 1 to N. Accordingly, the question for this study is: can we construct a ‘functional rarefaction curve’ that summarizes the expected functional dissimilarity between species when n plots are drawn at random from a larger pool of N plots? Methods: In this paper, we propose a parametric measure of functional diversity that is obtained by combining sample‐based rarefaction techniques that are usually applied to species richness with Rao's quadratic diversity. For a given set of N presence/absence plots, the resulting measure summarizes the expected functional dissimilarity at an increasingly larger cumulative number of plots n (n≤N). Results and Conclusions: Due to its parametric nature, the proposed measure is progressively more sensitive to rare species with increasing plot number, thus rendering this measure adequate for comparing the functional diversity of species assemblages that have been sampled with variable effort.     

Epigenetic variation predicts regional and local intraspecific functional diversity in a perennial herb

The ecological significance of epigenetic variation has been generally inferred from studies on model plants under artificial conditions, but the importance of epigenetic differences between individuals as a source of intraspecific diversity in natural plant populations remains essentially unknown. This study investigates the relationship between epigenetic variation and functional plant diversity by conducting epigenetic (methylation‐sensitive amplified fragment length polymorphisms, MSAP) and genetic (amplified fragment length polymorphisms, AFLP) marker–trait association analyses for 20 whole‐plant, leaf and regenerative functional traits in a large sample of wild‐growing plants of the perennial herb Helleborus foetidus from ten sampling sites in south‐eastern Spain. Plants differed widely in functional characteristics, and exhibited greater epigenetic than genetic diversity, as shown by per cent polymorphism of MSAP fragments (92%) or markers (69%) greatly exceeding that for AFLP ones (41%). After controlling for genetic structuring and possible cryptic relatedness, every functional trait considered exhibited a significant association with at least one AFLP or MSAP marker. A total of 27 MSAP (13.0% of total) and 12 AFLP (4.4%) markers were involved in significant associations, which explained on average 8.2% and 8.0% of trait variance, respectively. Individual MSAP markers were more likely to be associated with functional traits than AFLP markers. Between‐site differences in multivariate functional diversity were directly related to variation in multilocus epigenetic diversity after multilocus genetic diversity was statistically accounted for. Results suggest that epigenetic variation can be an important source of intraspecific functional diversity in H. foetidus, possibly endowing this species with the capacity to exploit a broad range of ecological conditions despite its modest genetic diversity.     

Trophic niche ontogeny and palaeoecology of early Toarcian Stenopterygius (Reptilia: Ichthyosauria)

Reconstructing ecological niche shifts during ontogeny in extinct animals with no living analogues is difficult without exceptional fossil collections. Here we demonstrate how a previously identified ontogenetic shift in the size and shape of the dentition in the early Toarcian ichthyosaur Stenopterygius quadriscissus accurately predicts a particular dietary shift. The smallest S. quadriscissus fed on small, burst‐swimming fishes, with a steady shift towards faster moving fish and cephalopods with increasing body size. Larger adult specimens appear to have been completely reliant on cephalopods, with fish completely absent from gut contents shortly after onset of sexual maturity. This is consistent with a previously proposed ontogenetic niche shift based on tooth shape and body size, corroborating the idea that dental ontogeny may be a useful predictor of dietary shifts in marine reptiles. Applying the theoretical framework used here to other extinct species will improve the resolution of palaeoecological reconstructions, where appropriate sample sizes exist.     

An ecological genetic delineation of local seed‐source provenance for ecological restoration

An increasingly important practical application of the analysis of spatial genetic structure within plant species is to help define the extent of local provenance seed collection zones that minimize negative impacts in ecological restoration programs. Here, we derive seed sourcing guidelines from a novel range‐wide assessment of spatial genetic structure of 24 populations of Banksia menziesii (Proteaceae), a widely distributed Western Australian tree of significance in local ecological restoration programs. An analysis of molecular variance (AMOVA) of 100 amplified fragment length polymorphism (AFLP) markers revealed significant genetic differentiation among populations (Φ_PT = 0.18). Pairwise population genetic dissimilarity was correlated with geographic distance, but not environmental distance derived from 15 climate variables, suggesting overall neutrality of these markers with regard to these climate variables. Nevertheless, Bayesian outlier analysis identified four markers potentially under selection, although these were not correlated with the climate variables. We calculated a global R‐statistic using analysis of similarities (ANOSIM) to test the statistical significance of population differentiation and to infer a threshold seed collection zone distance of ~60 km (all markers) and 100 km (outlier markers) when genetic distance was regressed against geographic distance. Population pairs separated by >60 km were, on average, twice as likely to be significantly genetically differentiated than population pairs separated by <60 km, suggesting that habitat‐matched sites within a 30‐km radius around a restoration site genetically defines a local provenance seed collection zone for B. menziesii. Our approach is a novel probability‐based practical solution for the delineation of a local seed collection zone to minimize negative genetic impacts in ecological restoration.     

Low levels of hybridization between sympatric Arctic char (Salvelinus alpinus) and Dolly Varden char (Salvelinus malma) highlights their genetic distinctiveness and ecological segregation

Understanding the extent of interspecific hybridization and how ecological segregation may influence hybridization requires comprehensively sampling different habitats over a range of life history stages. Arctic char (Salvelinus alpinus) and Dolly Varden (S. malma) are recently diverged salmonid fishes that come into contact in several areas of the North Pacific where they occasionally hybridize. To better quantify the degree of hybridization and ecological segregation between these taxa, we sampled over 700 fish from multiple lake (littoral and profundal) and stream sites in two large, interconnected southwestern Alaskan lakes. Individuals were genotyped at 12 microsatellite markers, and genetic admixture (Q) values generated through Bayesian‐based clustering revealed hybridization levels generally lower than reported in a previous study (<0.6% to 5% of samples classified as late‐generation hybrids). Dolly Varden and Arctic char tended to make different use of stream habitats with the latter apparently abandoning streams for lake habitats after 2–3 years of age. Our results support the distinct biological species status of Dolly Varden and Arctic char and suggest that ecological segregation may be an important factor limiting opportunities for hybridization and/or the ecological performance of hybrid char.