中华猕猴桃和美味猕猴桃自然居群遗传结构及其种间杂交渐渗

 利用9对SSR引物对中华猕猴桃（Actinidia chinensis）和美味猕猴桃（A. deliciosa）两近缘种的5个同域分布复合体和各自1个非同域分布居群进行了遗传多样性、居群遗传结构的分析以及种间杂交渐渗的探讨。结果表明：1）两物种共有等位基因比例高达81.13%,物种特有等位基因较少（中华猕猴桃：13.27%,美味猕猴桃：5.61%）,但共享等位基因表型频率在两近缘种间存在差异,而且与各同域复合体中两物种样本的交错程度或间距存在关联;2）两种猕猴桃均具有极高遗传多样性,美味猕猴桃的遗传多样性（H_o＝0 .749, PIC＝0.818）都略高于中华猕猴桃（H_o＝0.686,PIC＝0.799）;3）两猕 猴桃物种均具有较低的Nei’s居群遗传分化度,但AMOVA分析结果揭示种内异域居群间（F_ST=0.091 5）和同域复合体种间（F_ST=0.111 5）均存在一定程度的遗传分化;中华猕猴桃居群遗传分化（G_ST=0.086; F_ST=0.212 1）高于美味猕猴桃（G_ST= 0.080;F_ST=0.142 0）;4）同域分布复合体两物种间的遗传分化（G_ST=0.020）低于物种内异域居群间的遗传分化（中华猕猴桃：G_ST=0.086; 美味猕猴桃：G_ST=0.080）,同域复合体物种间的基因流（N_m＝7.89 -29.75）远远高于 同种异域居群间（中华猕猴桃：N_m =2.663; 美味猕猴桃：N_m=2.880）; 5）居群UPGMA聚类揭示在同一地域的居群优先聚类,个体聚类结果显示多数个体聚在各自居群组内,但各地理居群并不按地理距离的远近聚类,这与Mantel相关性检测所揭示的居群间遗传距离与地理距离没有显著性相关的结果一致。进一步分析表明两种猕猴桃的遗传多样性和居群遗传结构不仅受其广域分布、远交、晚期分化等生活史特性的影响,同时还与猕猴桃的染色体基数高 (x=29)、倍性复杂和种间杂交等因素密切相关,其中两种猕猴桃的共享祖先多态性和同域分布种间杂交基因渗透对两猕猴桃的居群遗传结构产生了重要影响。     

Genomewide introgressive hybridization patterns in wild Atlantic salmon influenced by inadvertent gene flow from hatchery releases

Many salmonid fish populations are threatened by genetic homogenization, primarily due to introgressive hybridization with hatchery‐reared conspecifics. By applying genomewide analysis using two molecular marker types (1986 SNPs and 17 microsatellites), we assessed the genetic impacts of inadvertent gene flow via straying from hatchery releases on wild populations of Atlantic salmon in the Gulf of Finland, Baltic Sea, over 16 years (1996–2012). Both microsatellites and SNPs revealed congruent population genetic structuring, indicating that introgression changed the genetic make‐up of wild populations by increasing genetic diversity and reducing genetic divergence. However, the degree of genetic introgression varied among studied populations, being higher in the eastern part and lower in the western part of Estonia, which most likely reflects the history of past stocking activities. Using kernel smoothing and permutation testing, we detected considerable heterogeneity in introgression patterns across the genome, with a large number of regions exhibiting nonrandom introgression widely dispersed across the genome. We also observed substantial variation in nonrandom introgression patterns within populations, as the majority of genomic regions showing elevated or reduced introgression were not consistently detected among temporal samples. This suggests that recombination, selection and stochastic processes may contribute to complex nonrandom introgression patterns. Our results suggest that (i) some genomic regions in Atlantic salmon are more vulnerable to introgressive hybridization, while others show greater resistance to unidirectional gene flow; and (ii) the hybridization of previously separated populations leads to complex and dynamic nonrandom introgression patterns that most likely have functional consequences for indigenous populations.     

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.     

Long live the alien: is high genetic diversity a pivotal aspect of crested porcupine (Hystrix cristata) long‐lasting and successful invasion?

Studying the evolutionary dynamics of an alien species surviving and continuing to expand after several generations can provide fundamental information on the relevant features of clearly successful invasions. Here, we tackle this task by investigating the dynamics of the genetic diversity in invasive crested porcupine (Hystrix cristata) populations, introduced to Italy about 1500 years ago, which are still growing in size, distribution range and ecological niche. Using genome‐wide RAD markers, we describe the structure of the genetic diversity and the demographic dynamics of the H. cristata invasive populations and compare their genetic diversity with that of native African populations of both H. cristata and its sister species, H. africaeaustralis. First, we demonstrate that genetic diversity is lower in both the invasive Italian and the North Africa source range relative to other native populations from sub‐Saharan and South Africa. Second, we find evidence of multiple introduction events in the invasive range followed by very limited gene flow. Through coalescence‐based demographic reconstructions, we also show that the bottleneck at introduction was mild and did not affect the introduced genetic diversity. Finally, we reveal that the current spatial expansion at the northern boundary of the range is following a leading‐edge model characterized by a general reduction of genetic diversity towards the edge of the expanding range. We conclude that the level of genome‐wide diversity of H. cristata invasive populations is less important in explaining its successful invasion than species‐specific life‐history traits or the phylogeographic history in the native source range.     

Evaluation of introgressive hybridization among Cervidae in Japan's Kinki District via two novel genetic markers developed from public NGS data

Hybridization and backcrossing of native populations with introduced species can lead to introgression and genetic alteration. In this study, we evaluated introgression in 43 deer from a potential hybrid zone around Okinoshima Island, Kinki District, Japan. This region witnessed the migration of a hybrid population (cross between the Formosan sika deer [Cervus nippon taiouanus] and other deer species) that could potentially breed with the native Japanese sika deer (C. n. centralis). We used an existing genetic marker for the mitochondrial cytochrome b gene and two novel markers for nuclear DNA, developed using publicly available next‐generation sequencing data. We identified one mainland deer with a mitochondrial haplotype identical to that of the Formosan sika deer as well as nuclear heterozygous sequences identical to those of Formosan and Japanese sika deer. This suggests that the mainland deer is a hybrid offspring of the Okinoshima population and native deer. However, only Japanese sika deer sequences were found in the other 42 samples, indicating limited introgression. Nevertheless, hybridization pre‐ and postintroduction in the Okinoshima population could cause multispecies introgression among Japanese sika deer, negatively affecting genetic integrity. We developed a simple test based on polymerase chain reaction–restriction fragment length polymorphism to detect introgression in natural populations. Our method can accelerate genetic monitoring of Japanese sika deer in Kinki District. In conclusion, to prevent further introgression and maintain genetic integrity of Japanese sika deer, we recommend establishing fences around Okinoshima Island to limit migration, besides a continued genetic monitoring of the native deer.     

Conservation genomics in the fight to help the recovery of the critically endangered Siamese crocodile Crocodylus siamensis

Endangered species are often characterized by low genetic diversity and it is imperative for conservation efforts to incorporate the knowledge obtained from genetic studies for effective management. However, despite the promise of technological advances in sequencing, application of genome‐wide data to endangered populations remains uncommon. In the present study we pursued a holistic conservation‐genomic approach to inform a field‐based management programme of a Critically Endangered species, the Siamese crocodile Crocodylus siamensis. Using thousands of single nucleotide polymorphisms from throughout the genome, we revealed signals of introgression from two other crocodile species within our sample of both wild and captive‐bred Siamese crocodiles from Cambodia. Our genetic screening of the Siamese crocodiles resulted in the subsequent re‐introduction of 12 individuals into the wild as well as the selection of four individuals for captive breeding programmes. Comparison of intraspecific genetic diversity revealed an alarmingly low contemporary effective population size in the wild (<50) with evidence of a recent bottleneck around Tonle Sap Lake. We also projected a probable future extinction in the wild (within fewer than five generations) in this population in the absence of re‐introduction efforts. However, an increase in the number of potential breeders through re‐introductions, including the one resulting from this project, could counter this trend. Our results have been implemented in ongoing re‐introduction and captive breeding programmes, with major implications for the conservation management of Siamese crocodiles, and provide a blueprint for the rescue effort of other “terminally ill” populations of critically endangered species.     

Investigating temporal changes in hybridization and introgression in a predominantly bimodal hybridizing population of invasive sika (Cervus nippon) and native red deer (C. elaphus) on the Kintyre Peninsula,Scotland

We investigated temporal changes in hybridization and introgression between native red deer (Cervus elaphus) and invasive Japanese sika (Cervus nippon) on the Kintyre Peninsula, Scotland, over 15 years, through analysis of 1513 samples of deer at 20 microsatellite loci and a mtDNA marker. We found no evidence that either the proportion of recent hybrids, or the levels of introgression had changed over the study period. Nevertheless, in one population where the two species have been in contact since ～1970, 44% of individuals sampled during the study were hybrids. This suggests that hybridization between these species can proceed fairly rapidly. By analysing the number of alleles that have introgressed from polymorphic red deer into the genetically homogenous sika population, we reconstructed the haplotypes of red deer alleles introduced by backcrossing. Five separate hybridization events could account for all the recently hybridized sika‐like individuals found across a large section of the Peninsula. Although we demonstrate that low rates of F1 hybridization can lead to substantial introgression, the progress of hybridization and introgression appears to be unpredictable over the short timescales.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Cultivation shapes genetic novelty in a globally important invader

Acacia saligna is a species complex that has become invasive in a number of countries worldwide where it has caused substantial environmental and economic impacts. Understanding genetic and other factors contributing to its success may allow managers to limit future invasions of closely related species. We used three molecular markers to compare the introduced range (South Africa) to the native range (Western Australia). Nuclear markers showed that invasive populations are divergent from native populations and most closely related to a cultivated population in Western Australia. We also found incongruence between nuclear and chloroplast data that, together with the long history of cultivation of the species, suggest that introgressive hybridization (coupled with chloroplast capture) may have occurred within A. saligna. While we could not definitively prove introgression, the genetic distance between cultivated and native A. saligna populations was comparable to known interspecific divergences among other Acacia species. Therefore, cultivation, multiple large‐scale introductions and possibly introgressive hybridization have rapidly given rise to the divergent genetic entity present in South Africa. This may explain the known global variation in invasiveness and inaccuracy of native bioclimatic models in predicting potential distributions.     

Effect of microsatellite selection on individual and population genetic inferences: an empirical study using cross‐specific and species‐specific amplifications

Although whole‐genome sequencing is becoming more accessible and feasible for nonmodel organisms, microsatellites have remained the markers of choice for various population and conservation genetic studies. However, the criteria for choosing microsatellites are still controversial due to ascertainment bias that may be introduced into the genetic inference. An empirical study of red deer (Cervus elaphus) populations, in which cross‐specific and species‐specific microsatellites developed through pyrosequencing of enriched libraries, was performed for this study. Two different strategies were used to select the species‐specific panels: randomly vs. highly polymorphic markers. The results suggest that reliable and accurate estimations of genetic diversity can be obtained using random microsatellites distributed throughout the genome. In addition, the results reinforce previous evidence that selecting the most polymorphic markers leads to an ascertainment bias in estimates of genetic diversity, when compared with randomly selected microsatellites. Analyses of population differentiation and clustering seem less influenced by the approach of microsatellite selection, whereas assigning individuals to populations might be affected by a random selection of a small number of microsatellites. Individual multilocus heterozygosity measures produced various discordant results, which in turn had impacts on the heterozygosity‐fitness correlation test. Finally, we argue that picking the appropriate microsatellite set should primarily take into account the ecological and evolutionary questions studied. Selecting the most polymorphic markers will generally overestimate genetic diversity parameters, leading to misinterpretations of the real genetic diversity, which is particularly important in managed and threatened populations.     

Founder events predict changes in genetic diversity during human‐mediated range expansions

Intentional or accidental introduction of species to new locations is predicted to result in loss of genetic variation and increase the likelihood of inbreeding, thus reducing population viability and evolutionary potential. However, multiple introductions and large founder numbers can prevent loss of genetic diversity and may therefore facilitate establishment success and range expansion. Based on a meta‐analysis of 119 introductions of 85 species of plants and animals, we here show a quantitative effect of founding history on genetic diversity in introduced populations. Both introduction of large number of individuals and multiple introduction events significantly contribute to maintaining or even increasing genetic diversity in introduced populations. The most consistent loss of genetic diversity is seen in insects and mammals, whereas introduced plant populations tend to have higher genetic variation than native populations. However, loss or gain of genetic diversity does not explain variation in the extent to which plant or animal populations become invasive outside of their native range. These results provide strong support for predictions from population genetics theory with respect to patterns of genetic diversity in introduced populations, but suggest that invasiveness is not limited by genetic bottlenecks.     

Limited gene exchange between two sister species of leaf beetles within a hybrid zone in the Alps

Analysing genomic variation within and between sister species is a first step towards understanding species boundaries. We focused on two sister species of cold‐resistant leaf beetles, Gonioctena quinquepunctata and G. intermedia, whose ranges overlap in the Alps. A previous study of DNA sequence variation had revealed multiple instances of mitochondrial genome introgression in this region, suggesting recent hybridization between the two species. To evaluate the extent of gene exchange resulting from these hybridization events, we sampled individuals of both species inside and outside the hybrid zone and analysed genomic variation among them using RAD‐seq markers. Individual levels of introgression in the nuclear genome were estimated first by defining species‐specific SNPs (displaying a fixed difference between species) a priori and second by using model‐based methods. Both types of analyses indicated little gene exchange, if any, between species at the level of the nuclear genome. Whereas the first method suggested slightly more gene flow, we argue that it has likely overestimated introgression in the phylogeographic context of this study. We conclude that strong intrinsic barriers prevent genetic exchange at the level of the nuclear genome between the two species. The apparent discrepancy observed between introgression occurring in the nuclear and mitochondrial genomes could be explained by selection acting in favour of the latter. Also, these results have consequences for the phylogeographic study of each species, since we can assume that genetic diversity in the overlapping portion of their ranges is not the product of introgression.     

Genetic diversity and population structure of Quercus fabri Hance in China revealed by genotyping‐by‐sequencing

Analysis of genetic diversity and population structure among Quercus fabri populations is essential for the conservation and utilization of Q. fabri resources. Here, the genetic diversity and structure of 158 individuals from 13 natural populations of Quercus fabri in China were analyzed using genotyping‐by‐sequencing (GBS). A total of 459,564 high‐quality single nucleotide polymorphisms (SNPs) were obtained after filtration for subsequent analysis. Genetic structure analysis revealed that these individuals can be clustered into two groups and the structure can be explained mainly by the geographic barrier, showed gene introgression from coastal to inland areas and high mountains could significantly hinder the mutual introgression of genes. Genetic diversity analysis indicated that the individual differences within groups are greater than the differences between the two groups. These results will help us better understand the genetic backgrounds of Q. fabri.     

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada,California, USA

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross‐breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short‐lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.     

Negligible nuclear introgression despite complete mitochondrial capture between two species of chipmunks

The idea that species boundaries can be semipermeable to gene flow is now widely accepted but the evolutionary importance of introgressive hybridization remains unclear. Here we examine the genomic contribution of gene flow between two hybridizing chipmunk species, Tamias ruficaudus and T. amoenus. Previous studies have shown that ancient hybridization has resulted in complete fixation of introgressed T. ruficaudus mitochondrial DNA (mtDNA) in some populations of T. amoenus, but the extent of nuclear introgression is not known. We used targeted capture to sequence over 10,500 gene regions from multiple individuals of both species. We found that most of the nuclear genome is sorted between these species and that overall genealogical patterns do not show evidence for introgression. Our analysis rules out all but very minor levels of interspecific gene flow, indicating that introgressive hybridization has had little impact on the overall genetic composition of these species outside of the mitochondrial genome. Given that much of the evidence for introgression in animals has come from mtDNA, our results underscore that unraveling the importance introgressive hybridization during animal speciation will require a genome‐wide perspective that is still absent for many species.     

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito‐genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito‐genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex‐biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt‐b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt‐b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial–nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt‐b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito‐nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex‐biased dispersal. A weak pattern of female‐biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex‐biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris‐type mito‐genome in thermal adaptation may also contribute to the introgression between the two species.     

Widespread mito‐nuclear discordance with evidence for introgressive hybridization and selective sweeps in Lycaeides

We investigated the extent and potential cause(s) of mitochondrial introgression within the polytypic North American Lycaeides species complex (Lepidoptera). By comparing population genetic structure based on mitochondrial DNA (COI and COII) and nuclear DNA (251 polymorphic amplified fragment length polymorphism markers), we detected substantial mito‐nuclear discordance, primarily involving a single mitochondrial haplotype (h01), which is likely due to mitochondrial introgression between differentiated Lycaeides populations and/or species. We detected reduced mitochondrial genetic diversity relative to nuclear genetic diversity in populations where mitochondrial haplotype h01 occurs, suggesting that the spread of this haplotype was facilitated by selection. We found no evidence that haplotype h01 is associated with increased fitness (in terms of survival to eclosion, fresh adult weight, and adult longevity) in a polymorphic Lycaeides melissa population. However, we did find a positive association between mitochondrial haplotype h01 and infection by the endoparasitic bacterium Wolbachia in one out of three lineages tested. Linkage disequilibrium between mitochondrial haplotype h01 and Wolbachia infection status may have resulted in indirect selection favouring the spread of haplotype h01 in at least one lineage of North American Lycaeides. These results illustrate the potential for introgressive hybridization to produce substantial mito‐nuclear discordance and demonstrate that an individual's mitochondrial and nuclear genome may have strikingly different evolutionary histories resulting from non‐neutral processes and intrinsic differences in the inheritance and biology of these genomes.     

Searching for the genetic footprint of ancient and recent hybridization

Determining the long‐term consequences of hybridization remains a central quest for evolutionary biologists. A particular challenge is to establish whether and to what extent widespread hybridization results in gene flow (introgression) between parental taxa. In this issue of Molecular Ecology, Jordan et al. ( 2018 ) search for evidence of gene flow between two closely related species of Geum (Rosaceae), which hybridize readily in contemporary populations and where hybrid swarms have been recorded for at least 200 years (Ruhsam, Hollingsworth, & Ennos, 2013 ). The authors find mixed evidence of ancient introgression when analysing allopatric populations. Intriguingly, when analysing populations of a region where the two species occur either mixed in the same population or in close proximity, and where hybrids are presently common, Jordan and colleagues find that the majority of randomly sampled individuals analysed (92/96) show no evidence of introgression (defined as individuals with admixture coefficients of <1%). The few individuals identified as hybrids are shown to likely be F1 or early‐generation backcrosses, indicating that even in sympatric regions, hybridization does not penetrate beyond a few generations. Based on their findings, Geum seems to be an example of little to no introgression despite contemporary hybridization.     

Minimum sample sizes for invasion genomics: Empirical investigation in an invasive whitefly

Analysis of population genetics provides insights into the evolutionary processes, among which the sample size choice is per se a crucial issue in the analysis. Genome‐wide high‐throughput techniques based on RADseq have been increasingly used in studies on the population genomics of invasive species. However, there is little information available regarding optimal sample sizes for analyzing population genomics of invasive species. In this study, we first use type IIB endonucleases restriction site‐associated DNA (2b‐RAD) to mine thousands of single nucleotide polymorphisms (SNPs) for native and introduced populations in Q1 clade (SPB and 17JN) and Q2 clade (ISQ and UAS0601) of the whitefly, Bemisia tabaci (Gennadius) MED (also known as B. tabaci biotype Q). Then, we used resampling techniques to create simulated populations with a random subset of individuals and 3,000 SNPs to determine how many individuals should be sampled for accurate estimates of intra‐ and interpopulation genetic diversity. We calculated the intrapopulation genetic diversity parameters (unbiased expected heterozygosity, observed heterozygosity, and the number of effect alleles) and pairwise genetic differentiation F_ST; finally, an ad hoc statistic, ΔK, was used to determine the optimal value. Our results showed that a sample size greater than four individuals (n ≥ 4) has little impact on estimates of genetic diversity within whitefly populations; moreover, precise estimate of F_ST can be easily achieved at a very small simple size (n = 3 or 4). Our results will provide in‐depth understanding of the optimization of sampling schemes in population genomics of invasive species.