A power analysis of microsatellite-based statistics for inferring past population growth

Insights into the Y chromosome human diversity in Uruguay

BACKGROUND

With regard to the origin of its population and microevolutionary processes, Uruguay exhibits distinctive features that distinguish it from other countries in Latin America, while at the same time sharing several similarities. In this article, we will focus on the variability of paternal genetic lineages in two geographical regions with different histories that can be considered as examples of distinct populations for the continent. In general terms, the genetic diversity is a result of different demographic processes related to the American conquest and colonisation. These resulted in distinct ancestral components which vary geographical and depend on the distribution by sex within these components. In Uruguay, native maternal haplogroups are significantly more frequent in the North. Although there are several studies about the geneticvariability of Uruguay, little is known about male genetic lineages.

AIMS

The aim of this work is to present an updated study of the male genetic variability of the Uruguayan population.

METHODS

We analyzed 13 biallelic markers and 27 STRs located in the male-specific region of the Y chromosome for 157 males: 98 from the capital, Montevideo, and 59 from Tacuarembó.

RESULTS

Almost all haplogroups found in both locations are European (99% and 93.2% respectively). One Sub-Saharan African haplogroup was found in Montevideo (1%) and 2 in Tacuarembó (3%), while Native haplogroups were found only in Tacuarembó, evidencing a strong sex-biased admixture. By crossing genetic and genealogical information we could relate European haplogroups with different waves and times of migrations.

DISCUSSION

Network analysis indicated a very diverse male population, suggesting that European migrants came from heterogeneous geographic locations and in different waves. Tacuarembó has closer population affinities with Iberian populations while Montevideo is more diverse. Male population expansion expansion, can be explained by the large number of migrants that arrived during the XIX century and the first half of the XX century.

CONCLUSIONS

The Uruguayan male gene pool is the result of several migration waves with diverse origins, with strong sex-biased admixture that can be explained by the European migration, the violence against the indigenous males, and the segregation of the Africansadmixture that can be explained due to European migration, violence against Natives, and segregation against African males.admixture that can be explained due to European migration, violence against Natives, and segregation against African males.admixture that can be explained due to European migration, violence against Natives, and segregation against African males.admixture that can be explained due to European migration, violence against Natives, and segregation of hte Africans.

The estimates of effective population size based on linkage disequilibrium are virtually unaffected by natural selection

The effective population size (N_e) is a key parameter to quantify the magnitude of genetic drift and inbreeding, with important implications in human evolution. The increasing availability of high-density genetic markers allows the estimation of historical changes in N_e across time using measures of genome diversity or linkage disequilibrium between markers. Directional selection is expected to reduce diversity and N_e, and this reduction is modulated by the heterogeneity of the genome in terms of recombination rate. Here we investigate by computer simulations the consequences of selection (both positive and negative) and recombination rate heterogeneity in the estimation of historical N_e. We also investigate the relationship between diversity parameters and N_e across the different regions of the genome using human marker data. We show that the estimates of historical N_e obtained from linkage disequilibrium between markers (N_eLD) are virtually unaffected by selection. In contrast, those estimates obtained by coalescence mutation-recombination-based methods can be strongly affected by it, which could have important consequences for the estimation of human demography. The simulation results are supported by the analysis of human data. The estimates of N_eLD obtained for particular genomic regions do not correlate, or they do it very weakly, with recombination rate, nucleotide diversity, proportion of polymorphic sites, background selection statistic, minor allele frequency of SNPs, loss of function and missense variants and gene density. This suggests that N_eLD measures mainly reflect demographic changes in population size across generations.

The inference of the demographic history of populations is of great relevance in evolutionary biology. This inference can be made from genomic data using coalescence methods or linkage disequilibrium methods. However, the assessment of these methods is usually made assuming neutrality (absence of selection). Here we show by computer simulations and analyses of human data that the estimates of historical effective population size obtained from linkage disequilibrium between markers are virtually unaffected by natural selection, either positive or negative. In contrast, estimates obtained by coalescence mutation-recombination-based methods can be strongly affected by it, which could have important consequences for recent estimations of human demography.

Genetic structure as a response to anthropogenic and extreme weather disturbances of a coastal dune dwelling spider, Arctosa sanctaerosae

The continued increase in the number of tourists visiting the Northern Gulf Coast (NGC), USA, in the last century, and the resulting sprawl of large cities along the coast, has degraded and fragmented the available habitat of Arctosa sanctaerosae, a wolf spider endemic to the secondary dunes of the white sandy beaches of the NGC. In addition to anthropogenic disturbance to this coastal region, hurricanes are an additional and natural perturbation to the ecosystem. The data presented here explore the status of populations of this species spanning the entire known range and the factors influencing population demography including anthropogenic disturbance and severe tropical storms. Using microsatellite markers, we were able to document the genetic structure of A. sanctaerosae, including current and historic patterns of migration. These results combined with ecological and census data reveal the characteristics that have influenced population persistence: ecological variables affecting the recovery of the population clusters after severe tropical storms, genetic fragmentation due to anthropogenic disturbance, and their interaction. These findings demonstrate the significance that the high traffic beach communities of the NGC and their impact on the once intact contiguous dune ecosystem have on recovery after severe tropical storms. Contemporary modeling methods that compare current and historic levels of gene flow suggest A. sanctaerosae has experienced a single, contiguous population subdivision, and the isolates reduced in size since the onset of commercial development of the NGC. These results point to the need for monitoring of the species and increased protection for this endangered habitat.

Extreme genetic diversity in asexual grass thrips populations

The continuous generation of genetic variation has been proposed as one of the main factors explaining the maintenance of sexual reproduction in nature. However, populations of asexual individuals may attain high levels of genetic diversity through within-lineage diversification, replicate transitions to asexuality from sexual ancestors and migration. How these mechanisms affect genetic variation in populations of closely related sexual and asexual taxa can therefore provide insights into the role of genetic diversity for the maintenance of sexual reproduction. Here, we evaluate patterns of intra- and interpopulation genetic diversity in sexual and asexual populations of Aptinothrips rufus grass thrips. Asexual A. rufus populations are found throughout the world, whereas sexual populations appear to be confined to few locations in the Mediterranean region. We found that asexual A. rufus populations are characterized by extremely high levels of genetic diversity, both in comparison with their sexual relatives and in comparison with other asexual species. Migration is extensive among asexual populations over large geographic distances, whereas close sexual populations are strongly isolated from each other. The combination of extensive migration with replicate evolution of asexual lineages, and a past demographic expansion in at least one of them, generated high local clone diversities in A. rufus. These high clone diversities in asexual populations may mimic certain benefits conferred by sex via genetic diversity and could help explain the extreme success of asexual A. rufus populations.

Population history and pathways of spread of the plant pathogen Phytophthora plurivora

Human activity has been shown to considerably affect the spread of dangerous pests and pathogens worldwide. Therefore, strict regulations of international trade exist for particularly harmful pathogenic organisms. Phytophthora plurivora, which is not subject to regulations, is a plant pathogen frequently found on a broad range of host species, both in natural and artificial environments. It is supposed to be native to Europe while resident populations are also present in the US. We characterized a hierarchical sample of isolates from Europe and the US and conducted coalescent-, migration, and population genetic analysis of sequence and microsatellite data, to determine the pathways of spread and the demographic history of this pathogen. We found P. plurivora populations to be moderately diverse but not geographically structured. High levels of gene flow were observed within Europe and unidirectional from Europe to the US. Coalescent analyses revealed a signal of a recent expansion of the global P. plurivora population. Our study shows that P. plurivora has most likely been spread around the world by nursery trade of diseased plant material. In particular, P. plurivora was introduced into the US from Europe. International trade has allowed the pathogen to colonize new environments and/or hosts, resulting in population growth.

COMPUTER PROGRAMS: onesamp: a program to estimate effective population size using approximate Bayesian computation

The estimation of effective population size from one sample of genotypes has been problematic because most estimators have been proven imprecise or biased. We developed a web-based program, onesamp that uses approximate Bayesian computation to estimate effective population size from a sample of microsatellite genotypes. onesamp requires an input file of sampled individuals' microsatellite genotypes along with information about several sampling and biological parameters. onesamp provides an estimate of effective population size, along with 95% credible limits. We illustrate the use of onesamp with an example data set from a re-introduced population of ibex Capra ibex.

Population genetic structure of sexual and parthenogenetic damselflies inferred from mitochondrial and nuclear markers

It has been postulated that obligate asexual lineages may persist in the long term if they escape from negative interactions with either sexual lineages or biological enemies; and thus, parthenogenetic populations will be more likely to occur in places that are difficult for sexuals to colonize, or those in which biological interactions are rare, such as islands or island-like habitats. Ischnura hastata is the only known example of natural parthenogenesis within the insect order Odonata, and it represents also a typical example of geographic parthenogenesis, as sexual populations are widely distributed in North America, whereas parthenogenetic populations of this species have only been found at the Azores archipelago. In order to gain insight in the origin and distribution of parthenogenetic I. hastata lineages, we have used microsatellites, mitochondrial and nuclear DNA sequence data, to examine the population genetic structure of this species over a wide geographic area. Our results suggest that sexual populations of I. hastata in North America conform to a large subdivided population that has gone through a recent spatial expansion. A recent single long distance dispersal event, followed by a demographic expansion, is the most parsimonious hypothesis explaining the origin of the parthenogenetic population of this species in the Azores islands.

Alternatives to the Wright-Fisher model: the robustness of mitochondrial Eve dating

Methods of calculating the distributions of the time to coalescence depend on the underlying model of population demography. In particular, the models assuming deterministic evolution of population size may not be applicable to populations evolving stochastically. Therefore the study of coalescence models involving stochastic demography is important for applications. One interesting approach which includes stochasticity is the O'Connell limit theory of genealogy in branching processes. Our paper explores how many generations are needed for the limiting distributions of O'Connell to become adequate approximations of exact distributions. We perform extensive simulations of slightly supercritical branching processes and compare the results to the O'Connell limits. Coalescent computations under the Wright-Fisher model are compared with limiting O'Connell results and with full genealogy-based predictions. These results are used to estimate the age of the so-called mitochondrial Eve, i.e., the root of the mitochondrial polymorphisms of the modern humans based on the DNA from humans and Neanderthal fossils.

An examination of genetic diversity and effective population size in Atlantic salmon populations

Effective population size (Ne) is an important parameter in the conservation of genetic diversity. Comparative studies of empirical data that gauge the relative accuracy of Ne methods are limited, and a better understanding of the limitations and potential of Ne estimators is needed. This paper investigates genetic diversity and Ne in four populations of wild anadromous Atlantic salmon (Salmo salar L.) in Europe, from the Rivers Oir and Scorff (France) and Spey and Shin (Scotland). We aimed to understand present diversity and historical processes influencing current population structure. Our results showed high genetic diversity for all populations studied, despite their wide range of current effective sizes. To improve understanding of high genetic diversity observed in the populations with low effective size, we developed a model predicting present diversity as a function of past demographic history. This suggested that high genetic diversity could be explained by a bottleneck occurring within recent centuries rather than by gene flow. Previous studies have demonstrated the efficiency of coalescence models to estimate Ne. Using nine subsets from 37 microsatellite DNA markers from the four salmon populations, we compared three coalescence estimators based on single and dual samples. Comparing Ne estimates confirmed the efficiency of increasing the number and variability of microsatellite markers. This efficiency was more accentuated for the smaller populations. Analysis with low numbers of neutral markers revealed uneven distributions of allelic frequencies and overestimated short-term Ne. In addition, we found evidence of artificial stock enhancement using native and non-native origin. We propose estimates of Ne for the four populations, and their applications for salmon conservation and management are discussed.

The distribution of coalescence times and distances between microsatellite alleles with changing effective population size

We investigate the effects of past changes of the effective population size on the present allelic diversity at a microsatellite marker locus. We first derive the analytical expression of the generating function of the joint probabilities of the time to the Most Recent Common Ancestor for a pair of alleles and of their distance (the difference in allele size). We give analytical solutions in the case of constant population size and the geometrical mutation model. Otherwise, numerical inversion allows the distributions to be calculated in general cases. The effects of population expansion or decrease and the possibility to detect an ancient bottleneck are discussed. The method is extended to samples of three and four alleles, which allows investigating the covariance structure of the frequencies f(k) of pairs of alleles with a size difference of k motifs, and suggesting some approaches to the estimation of past demography.

Correlation analyses reveal a substantial influence of allelic gaps on the investigation of genetic diversity of modern human populations with microsatellites

High intra-population genetic diversity and multiple measures of genetic variability at STR loci are useful in inferring past evolutionary history. However, STRs, categorized by their repeat motif size, differ in a number of aspects, requiring separate analyses. We analyzed 783 STRs in 36 worldwide populations to examine marker suitability as well as correlations between various measurements, to evaluate the extent of genomic diversity present in modern human populations. The loci were grouped by type and analyzed separately for each population group. Genetic variation defined by gene diversity and allele size variance, shows different trends of variation across four types of STRs. Additionally, there is little variation of genetic diversity, but there is decreased allelic size variance with increasing repeat motifs. A poor correlation between genetic diversity and allelic size variance across loci in all groups for Di-STRs is probably caused by the presence of allelic size gaps. In contrast, allelic size variance, genetic diversity, and number of alleles are strongly correlated with both tri- and tetra-STRs. The positive correlation of allelic size variance and presence of gaps within the range of allelic sizes in Di-STRs alone explains these observations. An unexpected high imbalance index (beta) at Di-STRs due to high allelic size variance also supports this assertion.

Microsatellite data suggest significant population structure and differentiation within the malaria vector Anopheles darlingi in Central and South America

Background

Anopheles darlingi is the most important malaria vector in the Neotropics. An understanding of A. darlingi's population structure and contemporary gene flow patterns is necessary if vector populations are to be successfully controlled. We assessed population genetic structure and levels of differentiation based on 1,376 samples from 31 localities throughout the Peruvian and Brazilian Amazon and Central America using 5–8 microsatellite loci.

Results

We found high levels of polymorphism for all of the Amazonian populations (mean R_S = 7.62, mean H_O = 0.742), and low levels for the Belize and Guatemalan populations (mean R_S = 4.3, mean H_O = 0.457). The Bayesian clustering analysis revealed five population clusters: northeastern Amazonian Brazil, southeastern and central Amazonian Brazil, western and central Amazonian Brazil, Peruvian Amazon, and the Central American populations. Within Central America there was low non-significant differentiation, except for between the populations separated by the Maya Mountains. Within Amazonia there was a moderate level of significant differentiation attributed to isolation by distance. Within Peru there was no significant population structure and low differentiation, and some evidence of a population expansion. The pairwise estimates of genetic differentiation between Central America and Amazonian populations were all very high and highly significant (F_ST = 0.1859 – 0.3901, P < 0.05). Both the D_A and F_ST distance-based trees illustrated the main division to be between Central America and Amazonia.

Conclusion

We detected a large amount of population structure in Amazonia, with three population clusters within Brazil and one including the Peru populations. The considerable differences in N_e among the populations may have contributed to the observed genetic differentiation. All of the data suggest that the primary division within A. darlingi corresponds to two white gene genotypes between Amazonia (genotype 1) and Central America, parts of Colombia and Venezuela (genotype 2), and are in agreement with previously published mitochondrial COI gene sequences interpreted as incipient species. Overall, it appears that two main factors have contributed to the genetic differentiation between the population clusters: physical distance between the populations and the differences in effective population sizes among the subpopulations.

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.

On the origin and spread of the Scab disease of apple: out of central Asia

BACKGROUND

Venturia inaequalis is an ascomycete fungus responsible for apple scab, a disease that has invaded almost all apple growing regions worldwide, with the corresponding adverse effects on apple production. Monitoring and predicting the effectiveness of intervention strategies require knowledge of the origin, introduction pathways, and population biology of pathogen populations. Analysis of the variation of genetic markers using the inferential framework of population genetics offers the potential to retrieve this information.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we present a population genetic analysis of microsatellite variation in 1,273 strains of V. inaequalis representing 28 orchard samples from seven regions in five continents. Analysis of molecular variance revealed that most of the variation (88%) was distributed within localities, which is consistent with extensive historical migrations of the fungus among and within regions. Despite this shallow population structure, clustering analyses partitioned the data set into separate groups corresponding roughly to geography, indicating that each region hosts a distinct population of the fungus. Comparison of the levels of variability among populations, along with coalescent analyses of migration models and estimates of genetic distances, was consistent with a scenario in which the fungus emerged in Central Asia, where apple was domesticated, before its introduction into Europe and, more recently, into other continents with the expansion of apple growing. Across the novel range, levels of variability pointed to multiple introductions and all populations displayed signatures of significant post-introduction increases in population size. Most populations exhibited high genotypic diversity and random association of alleles across loci, indicating recombination both in native and introduced areas.

CONCLUSIONS/SIGNIFICANCE

Venturia inaequalis is a model of invasive phytopathogenic fungus that has now reached the ultimate stage of the invasion process with a broad geographic distribution and well-established populations displaying high genetic variability, regular sexual reproduction, and demographic expansion.

Genetic diversity in diploid vs. tetraploidRorippa amphibia(Brassicaceae)

The frequency of polyploidy increases with latitude in the Northern Hemisphere, especially in deglaciated, recently colonized areas. The cause or causes of this pattern are largely unknown, but a greater genetic diversity of individual polyploid plants due to a doubled genome and/or a hybrid origin is seen as a likely factor underlying selective advantages related to life in extreme climates and/or colonization ability. A history of colonization in itself, as well as a recent origin, and possibly a limited number of polyploidization events would all predict less genetic diversity in polyploids than in diploids. The null hypothesis of higher gene diversity in polyploids has to date hardly been quantified and is here tested in self-incompatible Rorippa amphibia (Brassicaceae). The species occurs in diploid and tetraploid forms and displays clear geographical polyploidy in Europe. On the basis of eight microsatellite loci it can be concluded that the level of gene diversity is higher in tetraploids than in diploids, to an extent that is expected under neutral evolution when taking into account the larger effective population size in the doubled cytotype. There is thus no evidence for reduced genetic diversity in the tetraploids. The evidence presented here may mean that the tetraploids' origin is not recent, has not been affected by bottlenecks and/or that tetraploids were formed multiple times while an effect of introgression may also play a role.

Variation in estimated recombination rates across human populations

Recently it has been reported that recombination hotspots appear to be highly variable between humans and chimpanzees, and there is evidence for between-person variability in hotspots, and evolutionary transience. To understand the nature of variation in human recombination rates, it is important to describe patterns of variability across populations. Direct measurement of recombination rates remains infeasible on a large scale, and population-genetic approaches can be imprecise, and are affected by demographic history. Reports to date have suggested broad similarity in recombination rates at large genomic scales and across human populations. Here, we examine recombination rate estimates at a finer population and genomic scale: 28 worldwide populations and 107 SNPs in a 1 Mb stretch of chromosome 22q. We employ analysis of variance of recombination rate estimates, corrected for differences in effective population size using genome-wide microsatellite mutation rate estimates. We find substantial variation in fine-scale rates between populations, but reduced variation within continental groups. All effects examined (SNP-pair, region, population and interactions) were highly significant. Adjustment for effective population size made little difference to the conclusions. Observed hotspots tended to be conserved across populations, albeit at varying intensities. This holds particularly for populations from the same region, and also to a considerable degree across geographical regions. However, some hotspots appear to be population-specific. Several results from studies on the population history of humans are in accordance with our analysis. Our results suggest that between-population variation in DNA sequences may underly recombination rate variation.

Chloroplast microsatellites reveal colonization and metapopulation dynamics in the Canary Island pine

Chloroplast microsatellites are becoming increasingly popular markers for population genetic studies in plants, but there has been little focus on their potential for demographic inference. In this work the utility of chloroplast microsatellites for the study of population expansions was explored. First, we investigated the power of mismatch distribution analysis and the F(S) test with coalescent simulations of different demographic scenarios. We then applied these methods to empirical data obtained for the Canary Island pine (Pinus canariensis). The results of the simulations showed that chloroplast microsatellites are sensitive to sudden population growth. The power of the F(S) test and accuracy of demographic parameter estimates, such as the time of expansion, were reduced proportionally to the level of homoplasy within the data. The analysis of Canary Island pine chloroplast microsatellite data indicated population expansions for almost all sample localities. Demographic expansions at the island level can be explained by the colonization of the archipelago by the pine, while population expansions of different ages in different localities within an island could be the result of local extinctions and recolonization dynamics. Comparable mitochondrial DNA sequence data from a parasite of P. canariensis, the weevil Brachyderes rugatus, supports this scenario, suggesting a key role for volcanism in the evolution of pine forest communities in the Canary Islands.

Rangewide population genetic structure of the African malaria vector Anopheles funestus

Anopheles funestus is a primary vector of malaria in Africa south of the Sahara. We assessed its rangewide population genetic structure based on samples from 11 countries, using 10 physically mapped microsatellite loci, two per autosome arm and the X (N = 548), and 834 bp of the mitochondrial ND5 gene (N = 470). On the basis of microsatellite allele frequencies, we found three subdivisions: eastern (coastal Tanzania, Malawi, Mozambique and Madagascar), western (Burkina Faso, Mali, Nigeria and western Kenya), and central (Gabon, coastal Angola). A. funestus from the southwest of Uganda had affinities to all three subdivisions. Mitochondrial DNA (mtDNA) corroborated this structure, although mtDNA gene trees showed less resolution. The eastern subdivision had significantly lower diversity, similar to the pattern found in the codistributed malaria vector Anopheles gambiae. This suggests that both species have responded to common geographic and/or climatic constraints. The western division showed signatures of population expansion encompassing Kenya west of the Rift Valley through Burkina Faso and Mali. This pattern also bears similarity to A. gambiae, and may reflect a common response to expanding human populations following the development of agriculture. Due to the presumed recent population expansion, the correlation between genetic and geographic distance was weak. Mitochondrial DNA revealed further cryptic subdivision in A. funestus, not detected in the nuclear genome. Mozambique and Madagascar samples contained two mtDNA lineages, designated clade I and clade II, that were separated by two fixed differences and an average of 2% divergence, which implies that they have evolved independently for approximately 1 million years. Clade I was found in all 11 locations, whereas clade II was sampled only on Madagascar and Mozambique. We suggest that the latter clade may represent mtDNA capture by A. funestus, resulting from historical gene flow either among previously isolated and divergent populations or with a related species.

Divergence with gene flow in Anopheles funestus from the Sudan Savanna of Burkina Faso, West Africa

Anopheles funestus is a major vector of malaria across Africa. Understanding its complex and nonequilibrium population genetic structure is an important challenge that must be overcome before vector populations can be successfully perturbed for malaria control. Here we examine the role of chromosomal inversions in structuring genetic variation and facilitating divergence in Burkina Faso, West Africa, where two incipient species (chromosomal forms) of A. funestus, defined principally by rearrangements of chromosome 3R, have been hypothesized. Sampling across an approximately 300-km east-west transect largely contained within the Sudan-Savanna ecoclimatic zone, we analyzed chromosomal inversions, 16 microsatellite loci distributed genomewide, and 834 bp of the mtDNA ND5 gene. Both molecular markers revealed high genetic diversity, nearly all of which was accounted for by within-population differences among individuals, owing to recent population expansion. Across the study area there was no correlation between genetic and geographic distance. Significant genetic differentiation found between chromosomal forms on the basis of microsatellites was not genomewide but could be explained by chromosome 3R alone on the basis of loci inside and near inversions. These data are not compatible with complete reproductive isolation but are consistent with differential introgression and sympatric divergence between the chromosomal forms, facilitated by chromosome 3R inversions.

Molecular differentiation between chromosomally defined incipient species of Anopheles funestus

Anopheles funestus Giles is one of the most important vectors of malaria in sub-Saharan Africa. The population structure of this mosquito in Burkina Faso, West Africa based on chromosomal inversion data led to the description of two chromosomal forms, Kiribina and Folonzo. Because both forms co-occur in the same locales yet differ significantly, both in the frequency of inverted arrangements on chromosome arms 3R and 2R and in vectorial capacity, they were hypothesized to be emerging species with at least partial barriers to gene flow. This hypothesis would be strengthened by molecular evidence of differentiation between Kiribina and Folonzo at loci outside chromosomal inversions. We surveyed molecular variation in sympatric populations of the two forms using sequences from the mitochondrial ND5 gene and genotypes at sixteen microsatellite loci distributed across the genome. Both classes of marker revealed slight but significant differentiation between the two forms (mtDNA F(ST) = 0.023, P < 0.001; microsatellite F(ST) = 0.004, P < 0.001; R(st) = 0.009, P = 0.002). Locus-by-locus analysis of the microsatellite data showed that significant differentiation was not genome-wide, but could be attributed to five loci on chromosome 3R (F(ST) = 0.010, P < 0.001; R(st) = 0.016, P = 0.002). Importantly, three of these loci are outside of, and in linkage equilibrium with, chromosomal inversions, suggesting that differentiation between chromosomal forms extends beyond the inversions themselves. The slight overall degree of differentiation indicated by both marker classes is likely an underestimate because of recent population expansion inferred for both Folonzo and Kiribina. The molecular evidence from this study is consistent with the hypothesis of incipient speciation between Kiribina and Folonzo.

Directional mating and a rapid male population expansion in a hybrid Uruguayan population

The Uruguayan population has been considered as mainly European descent, with a negligible Native American or African contributions. Based on serological and molecular markers, recent studies demonstrate that these two populations had an important influence in the conformation of the present one. To the Northeastern region of Uruguay, a 20% Native American contribution was estimated using autosomal markers and a 62% Native American female origin based on mitochondrial markers. In this paper, we analyze four Y chromosome markers, two biallelic loci (M3 and YAP) and two microsatellites (DYS389I and DYS391), to characterize the male genetic contribution of a sample from the Northeastern city of Tacuarembó. We take different approaches to estimate the origin of male contributions to the population of Tacuarembó; Native American contribution ranges between 1.60% and 8.31%, confirming strong directional mating, which was also detected before with mitochondrial markers. Furthermore, the male population of Tacuarembó presents the characteristic of a population that suffered a bottleneck and a posterior expansion, confirmed using two microsatellite-based statistics to analyze the past population growth; patrilocality and migration could be responsible of those characteristics.

Asymptotic behavior of joint distributions of characteristics of a pair of randomly chosen individuals in discrete-time Fisher–Wright models with mutations and drift

This is a continuation of the series of articles (C.R. Rao, D.N. Shanbhag (Eds.), Handbook of Statistics 19: Stochastic Processes: Theory and Methods, Elsevier Science, Amsterdam, 2001 (Chapter 8); Math. Biosci. 175 (2002) 83; Math. Meth. Appl. Sci. 26 (2003) 1587; Adv. Appl. Probab. 36 (2004) 57) devoted to a study of the interplay between two of the main forces of population genetics, mutations and drift, in the Fisher-Wright model. We provide discrete-time versions of theorems describing asymptotic behavior of joint distributions of characteristics of a pair of individuals in this model; their continuous-time counterparts were presented in the previous papers. Furthermore, we show that imbalance index, introduced in Kimmel et al. (Genetics 148 (1998) 1921) and King et al. (Mol. Biol. Evol. 17(12) (2000) 1895) in the context of continuous-time models, may also be used in discrete-time models to detect past population growth.

An island within an island: genetic differentiation of Anopheles gambiae in São Tomé, West Africa, and its relevance to malaria vector control

Islands are choice settings for experimental studies of vector control strategies based on transgenic insects. Before considering this approach, knowledge of the population structure of the vector is essential. Genetic variation at 12 microsatellite loci was therefore studied in samples of the malaria vector Anopheles gambiae s.s., collected from six localities of São Tomé island (West Africa). The objectives were (i) to assess the demographic stability and effective population size of A. gambiae from these sites, (ii) to determine population differentiation and (iii) to relate the observed patterns of population structure with geographic, ecological and historical aspects of the vector on the island. Significant population differentiation, revealed by FST and RST statistics, was found between the southernmost site, Porto Alegre, and northern localities. The observed patterns of population substructure are probably a result of restrictions to gene flow in the less inhabited, more densely forested and mountainous south. In all localities surveyed, A. gambiae appeared to be experiencing a demographic expansion, consistent with a relatively recent (ca. 500 years) founder effect. The results are discussed with respect to current and future prospects of malaria vector control.

Clinal patterns of human Y chromosomal diversity in continental Italy and Greece are dominated by drift and founder effects

We explored the spatial distribution of human Y chromosomal diversity on a microgeographic scale, by typing 30 population samples from closely spaced locations in Italy and Greece for 9 haplogroups and their internal microsatellite variation. We confirm a significant difference in the composition of the Y chromosomal gene pools of the two countries. However, within each country, heterogeneity is not organized along the lines of clinal variation deduced from studies on larger spatial scales. Microsatellite data indicate that local increases of haplogroup frequencies can be often explained by a limited number of founders. We conclude that local founder or drift effects are the main determinants in shaping the microgeographic Y chromosomal diversity.

Features of evolution and expansion of modern humans, inferred from genomewide microsatellite markers

We study data on variation in 52 worldwide populations at 377 autosomal short tandem repeat loci, to infer a demographic history of human populations. Variation at di-, tri-, and tetranucleotide repeat loci is distributed differently, although each class of markers exhibits a decrease of within-population genetic variation in the following order: sub-Saharan Africa, Eurasia, East Asia, Oceania, and America. There is a similar decrease in the frequency of private alleles. With multidimensional scaling, populations belonging to the same major geographic region cluster together, and some regions permit a finer resolution of populations. When a stepwise mutation model is used, a population tree based on TD estimates of divergence time suggests that the branches leading to the present sub-Saharan African populations of hunter-gatherers were the first to diverge from a common ancestral population (approximately 71-142 thousand years ago). The branches corresponding to sub-Saharan African farming populations and those that left Africa diverge next, with subsequent splits of branches for Eurasia, Oceania, East Asia, and America. African hunter-gatherer populations and populations of Oceania and America exhibit no statistically significant signature of growth. The features of population subdivision and growth are discussed in the context of the ancient expansion of modern humans.

Population genetics of Caenorhabditis elegans: the paradox of low polymorphism in a widespread species

Caenorhabditis elegans has become one of the most widely used model research organisms, yet we have little information on evolutionary processes and recent evolutionary history of this widespread species. We examined patterns of variation at 20 microsatellite loci in a sample of 23 natural isolates of C. elegans from various parts of the world. One-half of the loci were monomorphic among all strains, and overall genetic variation at microsatellite loci was low, relative to most other species. Some population structure was detected, but there was no association between the genetic and geographic distances among different natural isolates. Thus, despite the nearly worldwide occurrence of C. elegans, little evidence was found for local adaptation in strains derived from different parts of the world. The low levels of genetic variation within and among populations suggest that recent colonization and population expansion might have occurred. However, the patterns of variation are not consistent with population expansion. A possible explanation for the observed patterns is the action of background selection to reduce polymorphism, coupled with ongoing gene flow among populations worldwide.

Approximate Bayesian computation in population genetics

We propose a new method for approximate Bayesian statistical inference on the basis of summary statistics. The method is suited to complex problems that arise in population genetics, extending ideas developed in this setting by earlier authors. Properties of the posterior distribution of a parameter, such as its mean or density curve, are approximated without explicit likelihood calculations. This is achieved by fitting a local-linear regression of simulated parameter values on simulated summary statistics, and then substituting the observed summary statistics into the regression equation. The method combines many of the advantages of Bayesian statistical inference with the computational efficiency of methods based on summary statistics. A key advantage of the method is that the nuisance parameters are automatically integrated out in the simulation step, so that the large numbers of nuisance parameters that arise in population genetics problems can be handled without difficulty. Simulation results indicate computational and statistical efficiency that compares favorably with those of alternative methods previously proposed in the literature. We also compare the relative efficiency of inferences obtained using methods based on summary statistics with those obtained directly from the data using MCMC.

Statistical properties of new neutrality tests against population growth

A number of statistical tests for detecting population growth are described. We compared the statistical power of these tests with that of others available in the literature. The tests evaluated fall into three categories: those tests based on the distribution of the mutation frequencies, on the haplotype distribution, and on the mismatch distribution. We found that, for an extensive variety of cases, the most powerful tests for detecting population growth are Fu's F(S) test and the newly developed R(2) test. The behavior of the R(2) test is superior for small sample sizes, whereas F(S) is better for large sample sizes. We also show that some popular statistics based on the mismatch distribution are very conservative.

Statistical properties of population differentiation estimators under stepwise mutation in a finite island model

Microsatellite loci mutate at an extremely high rate and are generally thought to evolve through a stepwise mutation model. Several differentiation statistics taking into account the particular mutation scheme of the microsatellite have been proposed. The most commonly used is R(ST) which is independent of the mutation rate under a generalized stepwise mutation model. F(ST) and R(ST) are commonly reported in the literature, but often differ widely. Here we compare their statistical performances using individual-based simulations of a finite island model. The simulations were run under different levels of gene flow, mutation rates, population number and sizes. In addition to the per locus statistical properties, we compare two ways of combining R(ST) over loci. Our simulations show that even under a strict stepwise mutation model, no statistic is best overall. All estimators suffer to different extents from large bias and variance. While R(ST) better reflects population differentiation in populations characterized by very low gene-exchange, F(ST) gives better estimates in cases of high levels of gene flow. The number of loci sampled (12, 24, or 96) has only a minor effect on the relative performance of the estimators under study. For all estimators there is a striking effect of the number of samples, with the differentiation estimates showing very odd distributions for two samples.

Non-homogeneous infinite sites model under demographic change: mathematical description and asymptotic behavior of pairwise distributions

We developed a mathematical model, which makes possible to predict joint distributions of numbers of mismatches in two or more linked regions of the genome, based on the Infinite Sites Models, under mutation-drift equilibrium as well as under various patterns of population growth. With mutation rates varying in the region, one of the predictions is different correlation between numbers of mismatches in the two regions, depending on the pattern of the past population growth (constant, slowly growing, or rapidly growing). Also, for slower growth patterns of population sizes, the coalescence tree is not necessarily 'starlike'. Thus, the joint distribution of mismatches, predicted by the model, provides additional insights into the demographic history of the populations. We also developed expectations and variances of sample statistics under different growth scenarios. As an application we used a sample of mitochondrial sequences from hypervariable regions 1 and 2 (HV1 and HV2), representing major world populations (Europeans, Asians and Africans). The patterns of joint distributions of numbers of mismatches differ markedly from one population to another. In addition, there is a considerable variability in the proportion of numbers of mismatches between HV1 and HV2 sequences. The patterns of bivariate distributions from the HV1 and HV2 data in these data are consistent with those generated by the model involving a stepwise change in population size.

Testing for genetic evidence of population expansion and contraction: an empirical analysis of microsatellite DNA variation using a hierarchical Bayesian model

The role of past climatic change in shaping the distributions of tropical rain forest vertebrates is central to long-standing hypotheses about the legacy of the Quaternary ice ages. One approach to testing such hypotheses is to use genetic data to infer the demographic history of codistributed species. Population genetic theory that relates the structure of allelic genealogies to historical changes in effective population size can be used to detect a past history of demographic expansion or contraction. The fruit bats Cynopterus sphinx and C. brachyotis (Chiroptera: Pteropodidae) exhibit markedly different distribution patterns across the Indomalayan region and therefore represent an exemplary species pair to use for such tests. The purpose of this study was to test alternative hypotheses about historical patterns of demographic expansion and contraction in C. sphinx and C. brachyotis using a coalescent-based analysis of microsatellite variation. Specifically, we used a hierarchical Bayesian model based on Markov chain Monte Carlo simulations to estimate the posterior distribution of genealogical and demographic parameters. The results revealed strong evidence for population contraction in both species. Evidence for a population contraction in C. brachyotis was expected on the basis of biogeographic considerations. However, similar evidence for population contraction in C. sphinx does not support the hypothesis that this species underwent a pronounced range expansion during the late Quaternary. Genetic evidence for population decline may reflect the consequences of habitat destruction on a more recent time scale.

Identifying main evolutionary mechanisms shaping genetic variation of Leuciscus cephalus L. 1758 (Cyprinidae) in western Greece: discordance between methods

Numerous methods can be used in intraspecific phylogeographic studies to infer the evolutionary mechanisms that shaped observed genetic variation in populations. However, these methods are scarcely used jointly, and the evolutionary outcomes they could propose are not fully compared. In this study, using a chub (Leuciscus cephalus; Cyprinidae) mitochondrial DNA data set (13 populations in Western Greece, 14 haplotypes), we compare three distinct 'historical' methods that could possibly infer relative importance of basic evolutionary mechanisms (isolation vs migration) shaping genetic variation: the nested clade analysis, the psi-test and the 'mismatch distributions'. Taking together, interpretations of these analyses allow to draw a picture of the evolutionary history of chub in Western Greece based on isolation and genetic drift for higher clades. However, results issued each method can differ for low differentiated clades. We discuss such differences and suggest that methods should be used jointly in phylogeographic studies for a better evaluation of the evolutionary mechanisms that shaped genetic variation.

DNA dinucleotide evolution in humans: fitting theory to facts

We examine length distributions of approximately 6000 human dinucleotide microsatellite loci, representing chromosomes 1-22, from the GDB database. Under the stepwise mutation model, results from theory and simulation are compared with the empirical data. In both constant and expanding population scenarios, a simple single-step model with parameters chosen to account for the observed variance of microsatellite lengths produces results inconsistent with the observed heterozygosity and the dispersion of length skewness. Complicating the model by allowing a variable mutation rate accounts for the homozygosity, and introducing a small probability of a large mutation step accounts for the dispersion in skewnesses. We discuss these results in light of the long-term evolution of microsatellites.

Evidence for recent population expansion in the evolutionary history of the malaria vectors Anopheles arabiensis and Anopheles gambiae

Gene flow in malaria vectors is usually estimated based on differentiation indices (e.g., F(ST)) in order to predict the contemporary spread of genes such as those conferring resistance to insecticides. This approach is reliant on a number of assumptions, the most crucial, and the one most likely to be violated in these species, being mutation-migration-drift equilibrium. Tests of this assumption for the African malaria vectors Anopheles gambiae and Anopheles arabiensis are the focus of this study. We analyzed variation at 18 microsatellite loci and the ND5 region of the mitochondrial genome in two populations of each species. Equilibrium was rejected by six of eight tests for the A. gambiae population from western Kenya and by three tests in eastern Kenya. In western Kenya, all departures from equilibrium were consistent with a recent population expansion, but in eastern Kenya, there were traces of a recent expansion and a bottleneck. Equilibrium was also rejected by two of the eight tests for both A. arabiensis populations; the departure from equilibrium was consistent with an expansion. These multiple-locus tests detected a genomewide effect and therefore a demographic event rather than a locus-specific effect, as would be caused by selection. Disequilibrium due to a recent expansion in these species implies that rates of gene flow, as inferred from differentiation indices, are overestimates as they include a historical component. We argue that the same effect applies to the majority of pest species due to the correlation of their demography with that of humans.