Population genetic structuring in pacu (Piaractus mesopotamicus) across the Paraná-Paraguay basin: evidence from microsatellites

The Paraná-Paraguay basin encompasses central western Brazil, northeastern Paraguay, eastern Bolivia and northern Argentina. The Pantanal is a flooded plain with marked dry and rainy seasons that, due to its soil characteristics and low declivity, has a great water holding capacity supporting abundant fish fauna. Piaractus mesopotamicus, or pacu, endemic of the Paraná-Paraguay basin, is a migratory species economically important in fisheries and ecologically as a potential seed disperser. In this paper we employ eight microsatellite loci to assess the population structure of 120 pacu sampled inside and outside the Pantanal of Mato Grosso. Our main objective was to test the null hypothesis of panmixia and to verify if there was a different structuring pattern between the Pantanal were there were no physical barriers to fish movement and the heavily impounded Paraná and Paranapanema rivers. All loci had moderate to high levels of polymorphism, the number of alleles varied from three to 18. The average observed heterozygosity varied from 0.068 to 0.911. After the Bonferroni correction three loci remained significant for deviations from Hardy-Weinberg, and for those the frequency of null alleles was estimated. F ST and R ST pairwise comparisons detected low divergence among sampling sites, and differentiation was significant only between Paranapanema and Cuiabá and Paranapanema and Taquari. No correlation between genetic distance and the natural logarithm of the geographic distance was detected. Results indicate that for conservation purposes and for restoration programs small genetic differences detected in the Cuiabá and Paranapanema rivers should be taken in consideration.

Isolation and characterization of nine microsatellite loci from the Hawaiian grouper Epinephelus quernus (Serranidae) for population genetic analyses

The availability of variable genetic markers for groupers (Serranidae) has generally been limited to mitochondrial DNA. For studies of population genetic structure, more loci are usually required; particularly useful are those that are nuclear in origin such as microsatellites. Here, we isolated and characterized 9 microsatellite loci from the endemic Hawaiian grouper Epinephelus quernus using a biotin-labeled oligonucleotide-streptavidin-coated magnetic bead approach. Of the 20 repeat-containing fragments isolated, 15 had sufficient flanking region in which to design primers. Among these, 9 produced consistent polymerase chain reaction product, and 6 were highly variable. These 6 loci were all composed of dinucleotide repeats, with the number of alleles ranging from 6 to 18, and heterozygosities from 33.3% to 91.7%. The high levels of variability observed should make these markers useful for population genetic studies of E. quernus, and potentially other epinephelines.

Population genetic structure of brown crab (Cancer pagurus) in Irish waters

The brown crab (Cancer pagurus) fishery in Ireland is one of the most important financially and socio-economically, with the species worth approximately €15m per year in the first half of the decade. Only mackerel (Scomber scombrus) and Dublin Bay prawn (Nephrops norvegicus) are of greater value. Despite this, very little research has been conducted to describe the stock structure of brown crab on a national scale. In this study a country-wide assessment of genetic population structure was carried out. Sampling was conducted from commercial fishing boats from 11/06 to 04/08 at seven sample sites representing the central Irish brown crab fisheries, with one sample site from the UK also included in the study. Six microsatellite markers, specifically developed for brown crab, were used to assess genetic diversity and estimate population differentiation parameters. Significant genetic structuring was found using F-statistics (Fst = 0.007) and exact tests, but not with Bayesian methods. Samples from the UK and Wexford were found to be genetically distinct from all other populations. Three northern populations from Malm Head and Stanton Bank were genetically similar with Fst estimates suggesting connectivity between them. Also, Stanton Bank, again on the basis of Fst estimates, appeared to be connected to populations down the west coast of Ireland, as far south as Kerry. Two Galway samples, one inside and one outside of Galway Bay, were genetically differentiated despite their close geographic proximity. It is hypothesised that a persistent northerly summer current could transport pelagic larvae from populations along the southwest and west coasts of Ireland towards Stanton Bank in the North, resulting in the apparent connectivity observed in this study.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

The effect of host social system on parasite population genetic structure: comparative population genetics of two ectoparasitic mites and their bat hosts.

BACKGROUND: The population genetic structure of a parasite, and consequently its ability to adapt to a given host, is strongly linked to its own life history as well as the life history of its host. While the effects of parasite life history on their population genetic structure have received some attention, the effect of host social system has remained largely unstudied. In this study, we investigated the population genetic structure of two closely related parasitic mite species (Spinturnix myoti and Spinturnix bechsteini) with very similar life histories. Their respective hosts, the greater mouse-eared bat (Myotis myotis) and the Bechstein's bat (Myotis bechsteinii) have social systems that differ in several substantial features, such as group size, mating system and dispersal patterns. RESULTS: We found that the two mite species have strongly differing population genetic structures. In S. myoti we found high levels of genetic diversity and very little pairwise differentiation, whereas in S. bechsteini we observed much less diversity, strongly differentiated populations and strong temporal turnover. These differences are likely to be the result of the differences in genetic drift and dispersal opportunities afforded to the two parasites by the different social systems of their hosts. CONCLUSIONS: Our results suggest that host social system can strongly influence parasite population structure. As a result, the evolutionary potential of these two parasites with very similar life histories also differs, thereby affecting the risk and evolutionary pressure exerted by each parasite on its host.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Population genetic structure and demography analysis in palinurus elephas populations from atlantic and mediterranean waters

Report for the scientific sojourn at the University of Reading, United Kingdom, from January until May 2008. The main objectives have been firstly to infer population structure and parameters in demographic models using a total of 13 microsatellite loci for genotyping approximately 30 individuals per population in 10 Palinurus elephas populations both from Mediterranean and Atlantic waters. Secondly, developing statistical methods to identify discrepant loci, possibly under selection and implement those methods using the R software environment. It is important to consider that the calculation of the probability distribution of the demographic and mutational parameters for a full genetic data set is numerically difficult for complex demographic history (Stephens 2003). The Approximate Bayesian Computation (ABC), based on summary statistics to infer posterior distributions of variable parameters without explicit likelihood calculations, can surmount this difficulty. This would allow to gather information on different demographic prior values (i.e. effective population sizes, migration rate, microsatellite mutation rate, mutational processes) and assay the sensitivity of inferences to demographic priors by assuming different priors.

Population genetic analysis of Colombian Trypanosoma cruzi isolates revealed by enzyme electrophoretic profiles

Although Colombia presents an enormous biological diversity, few studies have been conducted on the population genetics of Trypanosoma cruzi. This study was carried out with 23 Colombian stocks of this protozoa analyzed for 13 isoenzymatic loci. The Hardy-Weinberg equilibrium, the genetic diversity and heterogeneity, the genetic relationships and the possible spatial structure of these 23 Colombian stocks of T. cruzi were estimated. The majority of results obtained are in agreement with a clonal population structure. Nevertheless, two aspects expected in a clonal structure were not discovered in the Colombian T. cruzi stocks. There was an absence of given zymodemes over-represented from a geographical point of view and the presumed temporal stabilizing selective phenomena was not observed either in the Colombian stocks sampled several times through the years of the study. Some hypotheses are discussed in order to explain the results found.

Usefulness of microsatellite typing in population genetic studies of Trypanosoma cruzi

Through microsatellite analysis of 53 monoclonal populations of Trypanosoma cruzi, we found a remarkable degree of genetic polymorphism with no single multilocus genotype being observed more than once. The microsatellite profile proved to be stable during 70 generations of the CL Brener clone in culture. The microsatellite profiling presented also high diagnostic sensitivity since DNA amplifications could be achieved with less than 100 fg DNA, corresponding to half parasite total DNA content. Based on these technical attributes the microsatellite assay turns out to be an important tool for direct typing T. cruzi in biological samples. By using this approach we were able to type T. cruzi in feces of artificially infected bugs and in single cells sorted by FACS. The microsatellites have shown to be excellent markers for T. cruzi phylogenetic reconstruction. We used maximum parsimony based on the minimum number of mutational steps to build an unrooted Wagner network, which confirms previous conclusions based on the analysis of the D7 domain of the LSU rDNA gene that T. cruzi is composed by two major groups. We also obtained evidence that strains belonging to rRNA group 2 are subdivided into two genetically distant clusters, and that one of these clusters is more related to rRNA group 1/2. These results suggest different origins for these strains.

Population genetic structure of the dengue mosquito Aedes aegypti in Venezuela

The mosquito Aedes aegypti is the main vector of dengue in Venezuela. The genetic structure of this vector was investigated in 24 samples collected from eight geographic regions separated by up to 1160 km. We examined the distribution of a 359-basepair region of the NADH dehydrogenase subunit 4 mitochondrial gene among 1144 Ae. aegypti from eight collections. This gene was amplified by the polymerase chain reaction and tested for variation using single strand conformation polymorphism analysis. Seven haplotypes were detected throughout Venezuela and these were sorted into two clades. Significant differentiation was detected among collections and these were genetically isolated by distance.

Population genetic structure of the major malaria vector Anopheles darlingi (Diptera: Culicidae) from the Brazilian Amazon, using microsatellite markers

The population genetic structure of Anopheles darlingi, the major human malaria vector in the Neotropics, was examined using seven microsatellite loci from nine localities in central and western Amazonian Brazil. High levels of genetic variability were detected (5-25 alleles per locus; H E = 0.519-0.949). There was deviation from Hardy-Weinberg Equilibrium for 59.79% of the tests due to heterozygote deficits, while the analysis of linkage disequilibrium was significant for only two of 189 (1.05%) tests, most likely caused by null alleles. Genetic differentiation (F ST = 0.001-0.095; Nm = 4.7-363.8) indicates that gene flow is extensive among locations < 152 km apart (with two exceptions) and reduced, but not absent, at a larger geographic scale. Genetic and geographic distances were significantly correlated (R² = 0.893, P < 0.0002), supporting the isolation by distance (IBD) model. The overall estimate of Ne was 202.4 individuals under the linkage disequilibrium model, and 8 under the heterozygote excess model. Analysis of molecular variance showed that nearly all variation (~ 94%) was within sample locations. The UPGMA phenogram clustered the samples geographically, with one branch including 5/6 of the state of Amazonas localities and the other branch the Acre, Rondônia, and remaining Amazonas localities. Taken together, these data suggest little genetic structure for An. darlingi from central and western Amazonian Brazil. These findings also imply that the IBD model explains nearly all of the differentiation detected. In practical terms, populations of An. darlingi at distances < 152 km should respond similarly to vector control measures, because of high gene flow.

Microsatellite markers for population genetic studies of the blowfly Chrysomya putoria (Diptera: Calliphoridae)

The investigation of the genetic variation and population structure of Chrysomya species is of great interest for both basic and applied research. However, very limited genetic information is available for this genus across its geographical distribution. Here, we describe 12 polymorphic microsatellite loci isolated from Chrysomya putoria with expected heterozygosities ranging from 0.1402-0.8312. These markers are of potential applied interest for forensic entomologists and for the characterisation of the genetic structure of C. putoria from recently colonised regions, with great promise for understanding the colonisation dynamics and spread of the genus Chrysomya in the New World.

Strong reciprocity or strong ferocity? A population genetic view of the evolution of altruistic punishment.

Strong reciprocity, defined as a predisposition to help others and to punish those that are not helping, has been proposed as a potent force leading to the evolution of cooperation and altruism. However, the conditions under which strong reciprocity might be favored are not clear. Here we investigate the selective pressure on strong reciprocity by letting both limited dispersal (i.e., spatial structure) and recombination between helping and punishment jointly determine the evolutionary dynamics of strong reciprocity. Our analytical model suggests that when helping and punishment are perfectly linked traits (no recombination occurring between them), strong reciprocity can spread even when the initial frequency of strong reciprocators is close to 0 in the population (i.e., a rare mutant can invade). By contrast, our results indicate that when recombination can occur between helping and punishment (i.e., both traits coevolve) and is stronger than selection, punishment is likely to invade a population of defectors only when it gives a direct fitness benefit to the actor. Overall, our results delineate the conditions under which strong reciprocity is selected for in a spatially structured population and highlight that the forces behind its evolution involves kinship (be it genetic or cultural).

Forecasting changes in population genetic structure of alpine plants in response to global warming.

Species range shifts in response to climate and land use change are commonly forecasted with species distribution models based on species occurrence or abundance data. Although appealing, these models ignore the genetic structure of species, and the fact that different populations might respond in different ways because of adaptation to their environment. Here, we introduced ancestry distribution models, that is, statistical models of the spatial distribution of ancestry proportions, for forecasting intra-specific changes based on genetic admixture instead of species occurrence data. Using multi-locus genotypes and extensive geographic coverage of distribution data across the European Alps, we applied this approach to 20 alpine plant species considering a global increase in temperature from 0.25 to 4 °C. We forecasted the magnitudes of displacement of contact zones between plant populations potentially adapted to warmer environments and other populations. While a global trend of movement in a north-east direction was predicted, the magnitude of displacement was species-specific. For a temperature increase of 2 °C, contact zones were predicted to move by 92 km on average (minimum of 5 km, maximum of 212 km) and by 188 km for an increase of 4 °C (minimum of 11 km, maximum of 393 km). Intra-specific turnover-measuring the extent of change in global population genetic structure-was generally found to be moderate for 2 °C of temperature warming. For 4 °C of warming, however, the models indicated substantial intra-specific turnover for ten species. These results illustrate that, in spite of unavoidable simplifications, ancestry distribution models open new perspectives to forecast population genetic changes within species and complement more traditional distribution-based approaches.

Low dispersing species : population genetic processes in space and time : the genus "Trochulus" (Gastropoda: Hygromiidae) as a case study

RÉSUMÉ Une espèce est rarement composée d'une population unique. Parce que les individus ont des capacités de dispersion limitées et que les paysages sont des mosaïques d'habitats, la plupart des espèces sont plutôt composées de sous-populations connectées par la migration. Cette variation spatiale influence directement la distribution de la variabilité génétique dans et entre les populations. Durant ce travail, nous avons abordé certains des processus populationnels qui ont joué un rôle supposé dans l'apparition de nouvelles espèces au sein du genre Trochulus. Plus précisément, nous avons tenté d'évaluer les impacts respectifs de l'isolement passé (facteurs historiques) et présent (facteurs locaux). Nous avons d'abord pu montrer que les faibles capacités de dispersion des escargots terrestres ont directement influencé leur histoire évolutive à toutes les échelles spatiales et temporelles. En réduisant l'effet homogénéisant de la migration, une faible dispersion maintient dans les populations les traces génétiques d'évènements passés. A l'échelle de la distribution globale de Trochulus villosus, ces traces ont permis de reconstruire une histoire faite d'isolements et d'expansions de populations. En combinant des données génétiques avec une modélisation de la niche climatique passée, il a été possible de proposer un scénario significativement meilleur que toutes les hypothèses alternatives que nous avons testées. A l'échelle locale par contre, l'héritage historique est difficile à distinguer de la dynamique actuelle. Ce fut le cas des lignées mitochondriales du complexe sericeus-hispidus : les deux principales lignées étaient phylogénétiquement éloignées, avaient eu des démographies passées différentes et corrélaient avec des différences morphologiques. D'un autre côté, le flux de gène nucléaire était fort, contredisant l'idée de deux espèces cryptiques isolées reproductivement. Pour pouvoir conclure à la présence ou non de deux espèces, il nous a manqué des informations locales sur la dynamique des populations et les conditions écologiques que l'on trouve dans la région d'étude. Enfin, nous avons pu souligner que la connectivité entre populations d'escargots est soumise à la qualité des habitats et à leur organisation spatiale. Les escargots sont dépendants d'un habitat et s'y adaptent, comme l'indiquent la présence de «poils » uniquement sur la coquille d'espèces vivant dans des habitats humides ou la corrélation entre morphologie et habitat au sein du complexe sericeus-hispidus. Logiquement donc, les escargots migrent préférentiellement au travers d'habitats favorables comme l'a montré la réduction de flux de gènes au travers des prairies chez T. villosus (une espèce forestière). De ces données, nous pouvons supposer que les populations d'escargots en particulier, et des espèces à faible dispersion en général, ont de fortes chances d'être affectées par les changements climatiques, avec de probables implications pour leurs histoires évolutives. SUMMARY : Species rarely consists in a single population. Because individuals have limited dispersal abilities, because landscapes are habitat patchworks, most species are made of several subpopulations connected by migration. This spatial variation has consequences on the distribution of genetic diversity within and between populations, creating a structure among the populations. During the present work, we investigated some of the population processes assumed to have played an important role on the speciation within the genus Trochulus. More specifically, we questioned the respective impacts of past (historical factors) or present (local factors) population isolations. We first could show that the poor dispersal abilities of land snails have had profound impacts on their evolutionary histories at all spatial and temporal scales. Low dispersal maintains a strong signature of past events in the populations by minimising the homogenising effects of geneflow. At the scale of Trochulus villosus global distribution, they allowed to retrieve the detailed history of this species population isolations and expansions. Combining a large genetic dataset with paleo-climatic niche modelling ended up with a historical scenario significantly better than all traditional alternatives we tested. At local scale on the contrary, past events become difficult to tease apart from ongoing processes. This was the case for the divergent mitochondria) lineages within the sericeus-hispidus complex: the two principal lineages appeared to be phylogenetically distant, to have experienced different demographic histories and to correlate with morphological differences. On the other hand, nuclear (present day) geneflow was high, contradicting the idea of two reproductively isolated cryptic species. Information on the local population dynamics and environmental conditions are lacking to be able to decide whether past isolation has indeed resulted here in new species. Finally, we emphasised the importance of the habitat types present in a landscape as well as their spatial organisation for the population connectivity of land snails. These species are tightly dependent on a habitat and adapt to it as shown by thé occurrence of hair-like structures only in species living in humid environments or by the correlation between shell morphology and habitat in the sericeus-hispidus complex. As a result, land snails preferentially migrate through favourable habitats: Trochulus villosus, a forest species, had its geneflow significantly reduced across meadows. From these data, we can hypothesise that the populations of land snails in particular and of low dispersing species in general are likely to be strongly affected by the ongoing climate changes, with potential major consequences on their evolutionary histories.