Microsatellites reveal a high population structure in Triatoma infestans from Chuquisaca, Bolivia

Background

For Chagas disease, the most serious infectious disease in the Americas, effective disease control depends on elimination of vectors through spraying with insecticides. Molecular genetic research can help vector control programs by identifying and characterizing vector populations and then developing effective intervention strategies.

Methods and Findings

The population genetic structure of Triatoma infestans (Hemiptera: Reduviidae), the main vector of Chagas disease in Bolivia, was investigated using a hierarchical sampling strategy. A total of 230 adults and nymphs from 23 localities throughout the department of Chuquisaca in Southern Bolivia were analyzed at ten microsatellite loci. Population structure, estimated using analysis of molecular variance (AMOVA) to estimate F_ST (infinite alleles model) and R_ST (stepwise mutation model), was significant between western and eastern regions within Chuquisaca and between insects collected in domestic and peri-domestic habitats. Genetic differentiation at three different hierarchical geographic levels was significant, even in the case of adjacent households within a single locality (R _ST = 0.14, F _ST = 0.07). On the largest geographic scale, among five communities up to 100 km apart, R _ST = 0.12 and F _ST = 0.06. Cluster analysis combined with assignment tests identified five clusters within the five communities.

Conclusions

Some houses are colonized by insects from several genetic clusters after spraying, whereas other households are colonized predominately by insects from a single cluster. Significant population structure, measured by both R _ST and F _ST, supports the hypothesis of poor dispersal ability and/or reduced migration of T. infestans. The high degree of genetic structure at small geographic scales, inferences from cluster analysis and assignment tests, and demographic data suggest reinfesting vectors are coming from nearby and from recrudescence (hatching of eggs that were laid before insecticide spraying). Suggestions for using these results in vector control strategies are made.     

Population structure of the malaria vector Anopheles sinensis (Diptera: Culicidae) in China: two gene pools inferred by microsatellites

Background

Anopheles sinensis is a competent malaria vector in China. An understanding of vector population structure is important to the vector-based malaria control programs. However, there is no adequate data of A. sinensis population genetics available yet.

Methodology/Principal Findings

This study used 5 microsatellite loci to estimate population genetic diversity, genetic differentiation and demographic history of A. sinensis from 14 representative localities in China. All 5 microsatellite loci were highly polymorphic across populations, with high allelic richness and heterozygosity. Hardy–Weinberg disequilibrium was found in 12 populations associated with heterozygote deficits, which was likely caused by the presence of null allele and the Wahlund effect. Bayesian clustering analysis revealed two gene pools, grouping samples into two population clusters; one includes six and the other includes eight populations. Out of 14 samples, six samples were mixed with individuals from both gene pools, indicating the coexistence of two genetic units in the areas sampled. The overall differentiation between two genetic pools was moderate (F _ST = 0.156). Pairwise differentiation between populations were lower within clusters (F _ST = 0.008–0.028 in cluster I and F _ST = 0.004–0.048 in cluster II) than between clusters (F _ST = 0.120–0.201). A reduced gene flow (Nm = 1–1.7) was detected between clusters. No evidence of isolation by distance was detected among populations neither within nor between the two clusters. There are differences in effective population size (Ne = 14.3-infinite) across sampled populations.

Conclusions/Significance

Two genetic pools with moderate genetic differentiation were identified in the A. sinensis populations in China. The population divergence was not correlated with geographic distance or barrier in the range. Variable effective population size and other demographic effects of historical population perturbations could be the factors affecting the population differentiation. The structured populations may limit the migration of genes under pressures/selections, such as insecticides and immune genes against malaria.     

Population genetics and fine-scale genetic structure of Rhizopogon roseolus in the Tottori sand dune

We investigated the population genetics and fine-scale genetic structure of Rhizopogon roseolus. A total of 173 R. roseolus sporocarps were collected from two stands in the Tottori sand dune. We developed and applied five novel polymorphic microsatellite (SSR; simple sequence repeat) markers for sporocarp genotyping. In total, we identified 110 genets, most of which were small in size. Spatial autocorrelation analyses revealed a significantly positive genetic structure in short-distance classes. The inbreeding coefficient value was significant in both stands (F_IS = 0.18), while the F_ST value (F_ST = 0.020) indicated little genetic differentiation between the two populations. The majority of alleles were distributed in both stands with similar frequencies. These results suggest that short-distance spore dispersal plays a dominant role in generating new genets, and eventually increases the frequency of inbreeding in the Tottori sand dune, whereas rare gene flow between the two stands, possibly associated with spore dispersal by mycophagous animals, could reduce genetic differentiation.     

Three Decades of Farmed Escapees in the Wild: A Spatio-Temporal Analysis of Atlantic Salmon Population Genetic Structure throughout Norway

Each year, hundreds of thousands of domesticated farmed Atlantic salmon escape into the wild. In Norway, which is the world’s largest commercial producer, many native Atlantic salmon populations have experienced large numbers of escapees on the spawning grounds for the past 15–30 years. In order to study the potential genetic impact, we conducted a spatio-temporal analysis of 3049 fish from 21 populations throughout Norway, sampled in the period 1970–2010. Based upon the analysis of 22 microsatellites, individual admixture, F_ST and increased allelic richness revealed temporal genetic changes in six of the populations. These changes were highly significant in four of them. For example, 76% and 100% of the fish comprising the contemporary samples for the rivers Vosso and Opo were excluded from their respective historical samples at P = 0.001. Based upon several genetic parameters, including simulations, genetic drift was excluded as the primary cause of the observed genetic changes. In the remaining 15 populations, some of which had also been exposed to high numbers of escapees, clear genetic changes were not detected. Significant population genetic structuring was observed among the 21 populations in the historical (global F_ST = 0.038) and contemporary data sets (global F_ST = 0.030), although significantly reduced with time (P = 0.008). This reduction was especially distinct when looking at the six populations displaying temporal changes (global F_ST dropped from 0.058 to 0.039, P = 0.006). We draw two main conclusions: 1. The majority of the historical population genetic structure throughout Norway still appears to be retained, suggesting a low to modest overall success of farmed escapees in the wild; 2. Genetic introgression of farmed escapees in native salmon populations has been strongly population-dependent, and it appears to be linked with the density of the native population.     

Uniform selection as a primary force reducing population genetic differentiation of cavitation resistance across a species range

Background

Cavitation resistance to water stress-induced embolism determines plant survival during drought. This adaptive trait has been described as highly variable in a wide range of tree species, but little is known about the extent of genetic and phenotypic variability within species. This information is essential to our understanding of the evolutionary forces that have shaped this trait, and for evaluation of its inclusion in breeding programs.

Methodology

We assessed cavitation resistance (P ₅₀), growth and carbon isotope composition in six Pinus pinaster populations in a provenance and progeny trial. We estimated the heritability of cavitation resistance and compared the distribution of neutral markers (F _ST) and quantitative genetic differentiation (Q _ST), for retrospective identification of the evolutionary forces acting on these traits.

Results/Discussion

In contrast to growth and carbon isotope composition, no population differentiation was found for cavitation resistance. Heritability was higher than for the other traits, with a low additive genetic variance (h² _ns = 0.43±0.18, CV_A = 4.4%). Q _ST was significantly lower than F _ST, indicating uniform selection for P ₅₀, rather than genetic drift. Putative mechanisms underlying Q_STST are discussed.     

Population genetics of an ecosystem-defining reef coral Pocillopora damicornis in the Tropical Eastern Pacific

Background

Coral reefs in the Tropical Eastern Pacific (TEP) are amongst the most peripheral and geographically isolated in the world. This isolation has shaped the biology of TEP organisms and lead to the formation of numerous endemic species. For example, the coral Pocillopora damicornis is a minor reef-builder elsewhere in the Indo-West Pacific, but is the dominant reef-building coral in the TEP, where it forms large, mono-specific stands, covering many hectares of reef. Moreover, TEP P. damicornis reproduces by broadcast spawning, while it broods mostly parthenogenetic larvae throughout the rest of the Indo-West Pacific. Population genetic surveys for P. damicornis from across its Indo-Pacific range indicate that gene flow (i.e. larval dispersal) is generally limited over hundreds of kilometers or less. Little is known about the population genetic structure and the dispersal potential of P. damicornis in the TEP.

Methodology

Using multilocus microsatellite data, we analyzed the population structure of TEP P. damicornis among and within nine reefs and test for significant genetic structure across three geographically and ecologically distinct regions in Panama.

Principal Findings/Conclusions

We detected significant levels of population genetic structure (global R_ST = 0.162), indicating restricted gene flow (i.e. larvae dispersal), both among the three regions (R_RT = 0.081) as well as within regions (R_SR = 0.089). Limited gene flow across a distinct environmental cline, like the regional upwelling gradient in Panama, indicates a significant potential for differential adaptation and population differentiation. Individual reefs were characterized by unexpectedly high genet diversity (avg. 94%), relatively high inbreeding coefficients (global F_IS = 0.183), and localized spatial genetic structure among individuals (i.e. unique genets) over 10 m intervals. These findings suggest that gene flow is limited in TEP P. damicornis populations, particularly among regions, but even over meter scales within populations.     

Genetic Diversity in Introduced Golden Mussel Populations Corresponds to Vector Activity

We explored possible links between vector activity and genetic diversity in introduced populations of Limnoperna fortunei by characterizing the genetic structure in native and introduced ranges in Asia and South America. We surveyed 24 populations: ten in Asia and 14 in South America using the mitochondrial cytochrome c oxidase subunit I (COI) gene, as well as eight polymorphic microsatellite markers. We performed population genetics and phylogenetic analyses to investigate population genetic structure across native and introduced regions. Introduced populations in Asia exhibit higher genetic diversity (H _E = 0.667–0.746) than those in South America (H _E = 0.519–0.575), suggesting higher introduction effort for the former populations. We observed pronounced geographical structuring in introduced regions, as indicated by both mitochondrial and nuclear markers based on multiple genetic analyses including pairwise Ф_ST, F _ST, Bayesian clustering method, and three-dimensional factorial correspondence analyses. Pairwise F _ST values within both Asia (F _ST = 0.017–0.126, P = 0.000–0.009) and South America (F _ST = 0.004–0.107, P = 0.000–0.721) were lower than those between continents (F _ST = 0.180–0.319, P = 0.000). Fine-scale genetic structuring was also apparent among introduced populations in both Asia and South America, suggesting either multiple introductions of distinct propagules or strong post-introduction selection and demographic stochasticity. Higher genetic diversity in Asia as compared to South America is likely due to more frequent propagule transfers associated with higher shipping activities between source and donor regions within Asia. This study suggests that the intensity of human-mediated introduction vectors influences patterns of genetic diversity in non-indigenous species.     

Global phylogeography with mixed-marker analysis reveals male-mediated dispersal in the endangered scalloped hammerhead shark (Sphyrna lewini)

Background

The scalloped hammerhead shark, Sphyrna lewini, is a large endangered predator with a circumglobal distribution, observed in the open ocean but linked ontogenetically to coastal embayments for parturition and juvenile development. A previous survey of maternal (mtDNA) markers demonstrated strong genetic partitioning overall (global Φ_ST = 0.749) and significant population separations across oceans and between discontinuous continental coastlines.

Methodology/Principal Findings

We surveyed the same global range with increased sample coverage (N = 403) and 13 microsatellite loci to assess the male contribution to dispersal and population structure. Biparentally inherited microsatellites reveal low or absent genetic structure across ocean basins and global genetic differentiation (F _ST = 0.035) over an order of magnitude lower than the corresponding measures for maternal mtDNA lineages (Φ_ST = 0.749). Nuclear allelic richness and heterozygosity are high throughout the Indo-Pacific, while genetic structure is low. In contrast, allelic diversity is low while population structure is higher for populations at the ends of the range in the West Atlantic and East Pacific.

Conclusions/Significance

These data are consistent with the proposed Indo-Pacific center of origin for S. lewini, and indicate that females are philopatric or adhere to coastal habitats while males facilitate gene flow across oceanic expanses. This study includes the largest sampling effort and the most molecular loci ever used to survey the complete range of a large oceanic predator, and findings emphasize the importance of incorporating mixed-marker analysis into stock assessments of threatened and endangered shark species.     

Multi-Gene Analysis Reveals a Lack of Genetic Divergence between Calanus agulhensis and C. sinicus (Copepoda; Calanoida)

The discrimination and taxonomic identification of marine species continues to pose a challenge despite the growing number of diagnostic metrics and approaches. This study examined the genetic relationship between two sibling species of the genus Calanus (Crustacea; Copepoda; Calanidae), C. agulhensis and C. sinicus, using a multi-gene analysis. DNA sequences were determined for portions of the mitochondrial cytochrome c oxidase I (mtCOI); nuclear citrate synthase (CS), and large subunit (28S) rRNA genes for specimens collected from the Sea of Japan and North East (NE) Pacific Ocean for C. sinicus and from the Benguela Current and Agulhas Bank, off South Africa, for C. agulhensis. For mtCOI, C. sinicus and C. agulhensis showed similar levels of haplotype diversity (H_d = 0.695 and 0.660, respectively) and nucleotide diversity (π = 0.003 and 0.002, respectively). Pairwise F_ST distances for mtCOI were significant only between C. agulhensis collected from the Agulhas and two C. sinicus populations: the Sea of Japan (F_ST = 0.152, p<0.01) and NE Pacific (F_ST = 0.228, p<0.005). Between the species, F_ST distances were low for both mtCOI (F_ST = 0.083, p = 0.003) and CS (F_ST = 0.050, p = 0.021). Large subunit (28S) rRNA showed no variation between the species. Our results provide evidence of the lack of genetic distinction of C. sinicus and C. agulhensis, raise questions of whether C. agulhensis warrants status as a distinct species, and indicate the clear need for more intensive and extensive ecological and genetic analysis.     

Tertiary Origin and Pleistocene Diversification of Dragon Blood Tree (Dracaena cambodiana-Asparagaceae) Populations in the Asian Tropical Forests

Background

The origin of extraordinarily rich biodiversity in tropical forests is often attributed to evolution under stable climatic conditions over a long period or to climatic fluctuations during the recent Quaternary period. Here, we test these two hypotheses using Dracaena cambodiana, a plant species distributed in paleotropical forests.

Methods

We analyzed nucleotide sequence data of two chloroplast DNA (cpDNA: atpB-rbcL and trnD-trnT) regions and genotype data of six nuclear microsatellites from 15 populations (140 and 363 individuals, respectively) distributed in Indochina Peninsular and Hainan Island to infer the patterns of genetic diversity and phylogeographic structure. The population bottleneck and genetic drift were estimated based upon nuclear microsatellites data using the software programs BOTTLENECK and 2MOD. The lineage divergence times and past population dynamics based on cpDNA data were estimated using coalescent-based isolation-with-migration (IMa) and BEAST software programs.

Results

A significant phylogeographic structure (N _ST = 0.876, G _ST = 0.796, F _ST-SSR =0.329, R _ST =0.449; N _ST>G _ST, R _ST>F _ST-SSR, P<0.05) and genetic differentiation among populations were detected. Bottleneck analyses and Bayesian skyline plot suggested recent population reduction. The cpDNA haplotype network revealed the ancestral populations from the southern Indochina region expanded to northward. The most recent ancestor divergence time of D. cambodiana dated back to the Tertiary era and rapid diversification of terminal lineages corresponded to the Quaternary period.

Conclusions

The results indicated that the present distribution of genetic diversity in D. cambodiana was an outcome of Tertiary dispersal and rapid divergence during the Quaternary period under limited gene flow influenced by the uplift of Himalayan-Tibetan Plateau and Quaternary climatic fluctuations respectively. Evolutionary processes, such as extinction-recolonization during the Pleistocene may have contributed to the fast diversification in D. cambodiana.     

Genetic population structure in the Antarctic benthos: insights from the widespread amphipod, Orchomenella franklini

Currently there is very limited understanding of genetic population structure in the Antarctic benthos. We conducted one of the first studies of microsatellite variation in an Antarctic benthic invertebrate, using the ubiquitous amphipod Orchomenella franklini (Walker, 1903). Seven microsatellite loci were used to assess genetic structure on three spatial scales: sites (100 s of metres), locations (1–10 kilometres) and regions (1000 s of kilometres) sampled in East Antarctica at Casey and Davis stations. Considerable genetic diversity was revealed, which varied between the two regions and also between polluted and unpolluted sites. Genetic differentiation among all populations was highly significant (F _ST = 0.086, R _ST = 0.139, p<0.001) consistent with the brooding mode of development in O. franklini. Hierarchical AMOVA revealed that the majority of the genetic subdivision occurred across the largest geographical scale, with N_em≈1 suggesting insufficient gene flow to prevent independent evolution of the two regions, i.e., Casey and Davis are effectively isolated. Isolation by distance was detected at smaller scales and indicates that gene flow in O. franklini occurs primarily through stepping-stone dispersal. Three of the microsatellite loci showed signs of selection, providing evidence that localised adaptation may occur within the Antarctic benthos. These results provide insights into processes of speciation in Antarctic brooders, and will help inform the design of spatial management initiatives recently endorsed for the Antarctic benthos.     

Tracking the emergence of a new breed using 49,034 SNP in sheep

Domestic animals are unique in that they have been organised into managed populations called breeds. The strength of genetic divergence between breeds may vary dependent on the age of the breed, the scenario under which it emerged and the strength of reproductive isolation it has from other breeds. In this study, we investigated the Gulf Coast Native breed of sheep to determine if it contains lines of animals that are sufficiently divergent to be considered separate breeds. Allele sharing and principal component analysis (PCA) using nearly 50,000 SNP loci revealed a clear genetic division that corresponded with membership of either the Florida or Louisiana Native lines. Subsequent analysis aimed to determine if the strength of the divergence exceeded that found between recognised breed pairs. Genotypes from 14 breeds sampled from Europe and Asia were used to obtain estimates of pair-wise population divergence measured as F _ST. The divergence separating the Florida and Louisiana Native (F _ST = 6.2%) was approximately 50% higher than the average divergence separating breeds developed within the same region of Europe (F _ST = 4.2%). This strongly indicated that the two Gulf Coast Native lines are sufficiently different to be considered separate breeds. PCA using small SNP sets successfully distinguished between the Florida and Louisiana Native animals, suggesting that allele frequency differences have accumulated across the genome. This is consistent with a population history involving geographic separation and genetic drift. Suggestive evidence was detected for divergence at the poll locus on sheep chromosome 10; however drift at neutral markers has been the largest contributor to the genetic separation observed. These results document the emergence of populations that can be considered separate breeds, with practical consequences for bio-conservation priorities, animal registration and the establishment of separate breed societies.     

Genetic structure of the tree peony (Paeonia rockii) and the Qinling Mountains as a geographic barrier driving the fragmentation of a large population

Background

Tree peonies are great ornamental plants associated with a rich ethnobotanical history in Chinese culture and have recently been used as an evolutionary model. The Qinling Mountains represent a significant geographic barrier in Asia, dividing mainland China into northern (temperate) and southern (semi–tropical) regions; however, their flora has not been well analyzed. In this study, the genetic differentiation and genetic structure of Paeonia rockii and the role of the Qinling Mountains as a barrier that has driven intraspecific fragmentation were evaluated using 14 microsatellite markers.

Methodology/Principal Findings

Twenty wild populations were sampled from the distributional range of P. rockii. Significant population differentiation was suggested (F_ST value of 0.302). Moderate genetic diversity at the population level (H_S of 0.516) and high population diversity at the species level (H_T of 0.749) were detected. Significant excess homozygosity (F_IS of 0.076) and recent population bottlenecks were detected in three populations. Bayesian clusters, population genetic trees and principal coordinate analysis all classified the P. rockii populations into three genetic groups and one admixed Wenxian population. An isolation-by-distance model for P. rockii was suggested by Mantel tests (r = 0.6074, P<0.001) and supported by AMOVA (P<0.001), revealing a significant molecular variance among the groups (11.32%) and their populations (21.22%). These data support the five geographic boundaries surrounding the Qinling Mountains and adjacent areas that were detected with Monmonier''s maximum-difference algorithm.

Conclusions/Significance

Our data suggest that the current genetic structure of P. rockii has resulted from the fragmentation of a formerly continuously distributed large population following the restriction of gene flow between populations of this species by the Qinling Mountains. This study provides a fundamental genetic profile for the conservation and responsible exploitation of the extant germplasm of this species and for improving the genetic basis for breeding its cultivars.     

Linking seed dispersal and genetic structure of trees: a biogeographical approach

Aim Natural and human‐induced differences in frugivore assemblages can influence the seed dispersal distances of trees. An important issue in seed dispersal systems is to understand whether differences in seed dispersal distances also affect the genetic structure of mature trees. One possible approach to test for a relationship between seed dispersal and the genetic structure of mature trees is to compare the genetic structure of two closely related tree species between two biogeographical regions that differ in frugivore assemblages and seed dispersal distances. Previous studies on two Commiphora species revealed that Commiphora guillauminii in Madagascar has a much lower seed dispersal distance than Commiphora harveyi in South Africa. We tested whether the lower seed dispersal distance might have caused decreased gene flow, resulting in a stronger genetic structure in Madagascar than in South Africa. Location Madagascar and South Africa. Methods Using amplified fragment length polymorphism markers we investigated the genetic structure of 134 trees in Madagascar and 158 trees in South Africa at a local and a regional spatial scale. Results In concordance with our hypothesis, kinship analysis suggests that gene flow was restricted mostly to 3 km in Madagascar and to 30 km in South Africa. At the local spatial scale, the genetic differentiation among groups of trees within sample sites was marginally significantly higher in Madagascar (F_ST = 0.069) than in South Africa (F_ST = 0.021). However, at a regional spatial scale genetic differentiation was lower in Madagascar (F_ST = 0.053) than in South Africa (F_ST = 0.163). Main conclusions Our results show that lower seed dispersal distances of trees were linked to higher genetic differentiation of trees only at a local spatial scale. This suggests that seed dispersal affects the genetic population structure of trees at a local, but not at a regional, spatial scale.     

Quantitative genetics of body size and timing of maturation in two nine-spined stickleback (Pungitius pungitius) populations

Due to its influence on body size, timing of maturation is an important life-history trait in ectotherms with indeterminate growth. Comparison of patterns of growth and maturation within and between two populations (giant vs. normal sized) of nine-spined sticklebacks (Pungitius pungitius) in a breeding experiment revealed that the difference in mean adult body size between the populations is caused by differences in timing of maturation, and not by differential growth rates. The fish in small-sized population matured earlier than those from large-sized population, and maturation was accompanied by a reduction in growth rate in the small-sized population. Males matured earlier and at smaller size than females, and the fish that were immature at the end of the experiment were larger than those that had already matured. Throughout the experimental period, body size in both populations was heritable (h² = 0.10–0.64), as was the timing of maturation in the small-sized population (h² = 0.13–0.16). There was a significant positive genetic correlation between body size and timing of maturation at 140 DAH, but not earlier (at 80 or 110 DAH). Comparison of observed body size divergence between the populations revealed that Q _ST exceeded F _ST at older ages, indicating adaptive basis for the observed divergence. Hence, the results suggest that the body size differences within and between populations reflect heritable genetic differences in the timing of maturation, and that the observed body size divergence is adaptive.     

Estimates of Genetic Differentiation Measured by F(ST) Do Not Necessarily Require Large Sample Sizes When Using Many SNP Markers

Population genetic studies provide insights into the evolutionary processes that influence the distribution of sequence variants within and among wild populations. F_ST is among the most widely used measures for genetic differentiation and plays a central role in ecological and evolutionary genetic studies. It is commonly thought that large sample sizes are required in order to precisely infer F_ST and that small sample sizes lead to overestimation of genetic differentiation. Until recently, studies in ecological model organisms incorporated a limited number of genetic markers, but since the emergence of next generation sequencing, the panel size of genetic markers available even in non-reference organisms has rapidly increased. In this study we examine whether a large number of genetic markers can substitute for small sample sizes when estimating F_ST. We tested the behavior of three different estimators that infer F_ST and that are commonly used in population genetic studies. By simulating populations, we assessed the effects of sample size and the number of markers on the various estimates of genetic differentiation. Furthermore, we tested the effect of ascertainment bias on these estimates. We show that the population sample size can be significantly reduced (as small as n = 4–6) when using an appropriate estimator and a large number of bi-allelic genetic markers (k>1,000). Therefore, conservation genetic studies can now obtain almost the same statistical power as studies performed on model organisms using markers developed with next-generation sequencing.     

Genetic and morphological variability in the Raddia brasiliensis complex (Poaceae: Bambusoideae)

Raddia brasiliensis forms a species complex with the recently segregated R. megaphylla, R. lancifolia, R. soderstromii and R. stolonifera, occurring in the Atlantic rainforest, Brazil. Allozymic markers were used in 272 individuals of 14 populations of this group to investigate its genetic variability, correlating this with morphological variability, and testing the proposed taxonomy based on multivariate morphometrics. Genetic variability was low in almost all populations (P = 22.2–66.7, A = 1.3–2.0, H _e=0.04–0.17). R. brasiliensis showed a very high endogamy (F _IS = 0.329). Values for genetic and morphological structuring were very high to high (F _ST = 0.43, A _MRPP = 0.22; F _ST = 0.19, A _MRPP = 0.10 and F _ST = 0.18, A _MRPP = 0.39), respectively, for R. brasiliensis, R. soderstromii and R. megaphylla. The lowest genetic identity between populations was also found in R. brasiliensis, and the highest morphological differentiation was found between populations of R. megaphylla. Allozymic and morphological data were congruent and complementary, and confirm that we are dealing with five distinct species as previously circumscribed.     

Genome-wide analysis of single nucleotide polymorphisms uncovers population structure in Northern Europe

Background

Genome-wide data provide a powerful tool for inferring patterns of genetic variation and structure of human populations.

Principal Findings

In this study, we analysed almost 250,000 SNPs from a total of 945 samples from Eastern and Western Finland, Sweden, Northern Germany and Great Britain complemented with HapMap data. Small but statistically significant differences were observed between the European populations (F_ST = 0.0040, p<10⁻⁴), also between Eastern and Western Finland (F_ST = 0.0032, p<10⁻³). The latter indicated the existence of a relatively strong autosomal substructure within the country, similar to that observed earlier with smaller numbers of markers. The Germans and British were less differentiated than the Swedes, Western Finns and especially the Eastern Finns who also showed other signs of genetic drift. This is likely caused by the later founding of the northern populations, together with subsequent founder and bottleneck effects, and a smaller population size. Furthermore, our data suggest a small eastern contribution among the Finns, consistent with the historical and linguistic background of the population.

Significance

Our results warn against a priori assumptions of homogeneity among Finns and other seemingly isolated populations. Thus, in association studies in such populations, additional caution for population structure may be necessary. Our results illustrate that population history is often important for patterns of genetic variation, and that the analysis of hundreds of thousands of SNPs provides high resolution also for population genetics.     

Population genetics of two asexually and sexually reproducing psocids species inferred by the analysis of mitochondrial and nuclear DNA sequences

Background

The psocids Liposcelis bostrychophila and L. entomophila (Psocoptera: Liposcelididae) are found throughout the world and are often associated with humans, food stores and habitations. These insects have developed high levels of resistance to various insecticides in grain storage systems. However, the population genetic structure and gene flow of psocids has not been well categorized, which is helpful to plan appropriate strategies for the control of these pests.

Methodology/Principal Findings

The two species were sampled from 15 localities in China and analyzed for polymorphisms at the mitochondrial DNA (Cytb) and ITS (ITS1-5.8S-ITS2) regions. In total, 177 individual L. bostrychophila and 272 individual L. entomophila were analysed. Both Cytb and ITS sequences showed high genetic diversity for the two species with haplotype diversities ranged from 0.154±0.126 to 1.000±0.045, and significant population differentiation (mean F _ST = 0.358 for L. bostrychophila; mean F _ST = 0.336 for L. entomophila) was also detected among populations investigated. A Mantel test indicated that for both species there was no evidence for isolation-by-distance (IBD). The neutrality test and mismatch distribution statistics revealed that the two species might have undergone population expansions in the past.

Conclusion

Both L. bostrychophila and L. entomophila displayed high genetic diversity and widespread population genetic differentiation within and between populations. The significant population differentiation detected for both psocids may be mainly due to other factors, such as genetic drift, inbreeding or control practices, and less by geographic distance since an IBD effect was not found.