Spatial and Space-Time Correlations in Systems of Subpopulations with Genetic Drift and Migration

The geographic distribution of genetic variation is an important theoretical and experimental component of population genetics. Previous characterizations of genetic structure of populations have used measures of spatial variance and spatial correlations. Yet a full understanding of the causes and consequences of spatial structure requires complete characterization of the underlying space-time system. This paper examines important interactions between processes and spatial structure in systems of subpopulations with migration and drift, by analyzing correlations of gene frequencies over space and time. We develop methods for studying important features of the complete set of space-time correlations of gene frequencies for the first time in population genetics. These methods also provide a new alternative for studying the purely spatial correlations and the variance, for models with general spatial dimensionalities and migration patterns. These results are obtained by employing theorems, previously unused in population genetics, for space-time autoregressive (STAR) stochastic spatial time series. We include results on systems with subpopulation interactions that have time delay lags (temporal orders) greater than one. We use the space-time correlation structure to develop novel estimators for migration rates that are based on space-time data (samples collected over space and time) rather than on purely spatial data, for real systems. We examine the space-time and spatial correlations for some specific stepping stone migration models. One focus is on the effects of anisotropic migration rates. Partial space-time correlation coefficients can be used for identifying migration patterns. Using STAR models, the spatial, space-time, and partial space-time correlations together provide a framework with an unprecedented level of detail for characterizing, predicting and contrasting space-time theoretical distributions of gene frequencies, and for identifying features such as the pattern of migration and estimating migration rates in experimental studies of genetic variation over space and time.     

The influence of contemporary and historic landscape features on the genetic structure of the sand dune endemic,Cirsium pitcheri (Asteraceae)

Narrow endemics are at risk from climate change because of their restricted habitat preferences, lower colonization ability and dispersal distances. Landscape genetics combines new tools and analyses that allow us to test how both past and present landscape features have facilitated or hindered previous range expansion and local migration patterns, and thereby identifying potential limitations to future range shifts. We have compared current and historic habitat corridors in Cirsium pitcheri, an endemic of the linear dune ecosystem of the Great Lakes, to determine the relative contributions of contemporary migration and post-glacial range expansion on genetic structure. We used seven microsatellite loci to characterize the genetic structure for 24 populations of Cirsium pitcheri, spanning the center to periphery of the range. We tested genetic distance against different measures of geographic distance and landscape permeability, based on contemporary and historic landscape features. We found moderate genetic structure (Fst=0.14), and a north–south pattern to the distribution of genetic diversity and inbreeding, with northern populations having the highest diversity and lowest levels of inbreeding. High allelic diversity, small average pairwise distances and mixed genetic clusters identified in Structure suggest that populations in the center of the range represent the point of entry to the Lake Michigan and a refugium of diversity for this species. A strong association between genetic distances and lake-level changes suggests that historic lake fluctuations best explain the broad geographic patterns, and sandy habitat best explains local patterns of movement.     

Livestock trade history, geography, and parasite strains: the mitochondrial genetic structure of Echinococcus granulosus in Argentina

A sample of 114 isolates of Echinococcus granulosus (Cestoda: Taeniidae) collected from different host species and sites in Argentina has been sequenced for 391 bp from the mitochondrial cytochrome c oxidase subunit I gene to analyze genetic variability and population structure. Nine different haplotypes were identified, 5 of which correspond to already characterized strains. Analysis of molecular variance and nested clade analysis of the distribution of haplotypes among localities within 3 main geographic regions indicate that geographic differentiation accounts for the overall pattern of genetic variability in E. granulosus populations. Significant geographic differentiation is also present when the sheep strain alone is considered. Our results suggest that geographic patterns are not due to actual restricted gene flow between regions but are rather a consequence of past history, probably related to the time and origin of livestock introduction in Argentina.     

TESTING MODELS OF MIGRATION AND ISOLATION AMONG POPULATIONS OF CHINOOK SALMON (ONCORHYNCHUS TSCHAWYTSCHA)

The chinook salmon (Oncorhynchus tschawytscha) is a behaviorally, morphologically, and ecologically variable species distributed over a large geographic range. Although previous genetic surveys have revealed considerable genetic differences among populations with different life history types and from different major river drainages, it is not clear to what degree these genetically distinct populations are connected by low levels of gene flow. The work described in this paper addresses this question by surveying DNA restriction site variation at six nuclear genes from nine populations encompassing most of the species's North American range, and then attempting to fit the patterns of variation observed at these genes to five different evolutionary models using computer simulations of the coalescent process. Two commonly used constant population size models, one hypothesizing no migration among populations and one hypothesizing equal rates of migration among populations, were found to be statistically inconsistent with the observed patterns of variation. The other three models, which involved either recent divergence among populations coupled with large changes in populations size, unequal migration rates among populations, or selection, were all found to be consistent with the observed patterns of variation. Assuming selective neutrality, these results suggest that either the populations have all descended from a common ancestral population within the last ~50,000 years and have all suffered large declines in effective population size since that time, or that they have a more ancient divergence time but are connected by low levels of gene flow. These conclusions rest on the assumption of selective neutrality. With the methods employed, it was not possible to simultaneously test hypotheses of both selective neutrality and population structure.     

Disruption of historical patterns of isolation by distance in coastal steelhead

When natural populations exchange migrants at a rate proportional to their geographic distance, population genetics theory leads to the expectation of a pattern of isolation-by-distance (IBD), whereby geographic and genetic distance are correlated. However, the presence or absence of such patterns in modern populations may not fully reflect the historical relationships among those populations. Thus, historical samples, collected prior to modern human impacts, can often provide a critical baseline for comparison with modern populations. Steelhead, the anadromous form of rainbow trout, Oncorhynchus mykiss, are native to western North America and are endangered or threatened throughout most of California, near the southern extent of their native range. Population samples of steelhead collected in 1897 and 1909 in Central California rivers provided the opportunity to evaluate the historical genetic composition and population structure of these threatened fish. Here we show that these steelhead populations had a historically strong correlation between genetic and geographic distance that has been virtually erased in modern populations, suggesting that current relationships among modern steelhead populations are no longer reflective of natural migratory pathways. This demonstrates the critical role of migration in maintaining population relationships of threatened species and highlights the importance of natural history museums in providing historical baseline information.     

Hitchhiking effect of a beneficial mutation spreading in a subdivided population

A central problem in population genetics is to detect and analyze positive natural selection by which beneficial mutations are driven to fixation. The hitchhiking effect of a rapidly spreading beneficial mutation, which results in local removal of standing genetic variation, allows such an analysis using DNA sequence polymorphism. However, the current mathematical theory that predicts the pattern of genetic hitchhiking relies on the assumption that a beneficial mutation increases to a high frequency in a single random-mating population, which is certainly violated in reality. Individuals in natural populations are distributed over a geographic space. The spread of a beneficial allele can be delayed by limited migration of individuals over the space and its hitchhiking effect can also be affected. To study this effect of geographic structure on genetic hitchhiking, we analyze a simple model of directional selection in a subdivided population. In contrast to previous studies on hitchhiking in subdivided populations, we mainly investigate the range of sufficiently high migration rates that would homogenize genetic variation at neutral loci. We provide a heuristic mathematical analysis that describes how the genealogical structure at a neutral locus linked to the locus under selection is expected to change in a population divided into two demes. Our results indicate that the overall strength of genetic hitchhiking--the degree to which expected heterozygosity decreases--is diminished by population subdivision, mainly because opportunity for the breakdown of hitchhiking by recombination increases as the spread of the beneficial mutation across demes is delayed when migration rate is much smaller than the strength of selection. Furthermore, the amount of genetic variation after a selective sweep is expected to be unequal over demes: a greater reduction in expected heterozygosity occurs in the subpopulation from which the beneficial mutation originates than in its neighboring subpopulations. This raises a possibility of detecting a "hidden" geographic structure of population by carefully analyzing the pattern of a selective sweep.     

INTEGRATIVE TESTING OF HOW ENVIRONMENTS FROM THE PAST TO THE PRESENT SHAPE GENETIC STRUCTURE ACROSS LANDSCAPES

Tests of the genetic structure of empirical populations typically focus on the correlative relationships between population connectivity and geographic and/or environmental factors in landscape genetics. However, such tests may overlook or misidentify the impact of candidate factors on genetic structure, especially when connectivity patterns differ between past and present populations because of shifting environmental conditions over time. Here we account for the underlying demographic component of population connectivity associated with a temporarily dynamic landscape in tests of the factors structuring population genetic variation in an Australian lizard, Lerista lineopunctulata, from 24 nuclear loci. Correlative tests did not support significant effect from factors associated with a static contemporary landscape. However, spatially explicit demographic modeling of genetic differentiation shows that changes in environmental conditions (as estimated from paleoclimatic data) and corresponding distributional shifts from the past to present landscape significantly structures genetic variation. Results from model‐based inference (i.e., from an integrative modeling approach that generates spatially explicit expectations that are tested with approximate Bayesian computation) contrasts with those from correlative analyses, highlighting the importance of expanding the landscape genetic perspective to tests the links between pattern and process, revealing how factors shape patterns of genetic variation within species.     

Among population differentiation at nuclear genes in native Scots pine (Pinus sylvestris L.) in Scotland

In the Scottish Highlands, Scots pine is at the north-western extreme of its wide natural distribution. Here, the remaining native populations are patchily distributed in highly variable environments, from the more continental, drier eastern Highlands to the milder, wetter Atlantic Ocean coast. As these pinewoods are the remnants of a naturally established forest, they form a valuable system for analysis of genetic and adaptive variation in heterogeneous environments. Using samples from across the Scottish population, we analysed data from nuclear and mitochondrial genes to assess patterns of within and between population genetic variation. Within population diversity levels were high, and significant genetic differentiation among pairs of Scottish populations at relatively small spatial scales was present at several nuclear loci. At these loci, no differentiation had been found among continental populations, even those separated by large geographic distances. Overall, no clear clustering of Scottish samples was found in population structure analysis suggesting that geographically distant populations with high intra-population nucleotide diversity are not strongly isolated or diverged from each other. Scottish populations lacked a mitotype that is widespread in eastern and north-eastern Europe, indicating that pines from that area may not have participated in the most recent colonisation of the British Isles.     

Inter‐island divergence within Drosophila mauritiana,a species of the D. simulans complex: Past history and/or speciation in progress?

Speciation with gene flow may be more common than generally thought, which makes detailed understanding of the extent and pattern of genetic divergence between geographically isolated populations useful. Species of the Drosophila simulans complex provide a good model for speciation and evolutionary studies, and hence understanding their population genetic structure will increase our understanding of the context in which speciation has occurred. Here, we describe genetic diversity and genetic differentiation of two distant populations of D. mauritiana (Mauritius and Rodrigues Islands) at mitochondrial and nuclear loci. We surveyed the two populations for their mitochondrial haplotypes, eight nuclear genes and 18 microsatellite loci. A new mitochondrial type is fixed in the Rodrigues population of D. mauritiana. The two populations are highly differentiated, their divergence appears relatively ancient (100,000 years) compared to the origin of the species, around 0.25MYA, and they exhibit very limited gene flow. However, they have similar levels of divergence from their sibling, D. simulans. Both nuclear genes and microsatellites revealed contrasting demographic histories between the two populations, expansion for the Mauritius population and stable population size for the Rodrigues Island population. The discovery of pronounced geographic structure within D. mauritiana combined to genetic structuring and low gene flow between the two island populations illuminates the evolutionary history of the species and clearly merits further attention in the broad context of speciation.     

Evidence for Inbreeding and Genetic Differentiation among Geographic Populations of the Saprophytic Mushroom Trogia venenata from Southwestern China

During the past 40 years, more than 400 Sudden Unexplained Deaths (SUDs) have occurred in Yunnan, southwestern China. Epidemiological and toxicological analyses suggested that a newly discovered mushroom called Trogia venenata was the leading culprit for SUDs. At present, relatively little is known about the genetics and natural history of this mushroom. In this study, we analyzed the sequence variation at four DNA fragments among 232 fruiting bodies of T. venenata collected from seven locations. Our ITS sequence analyses confirmed that all the isolates belonged to the same species. The widespread presence of sequence heterozygosity within many strains at each of three protein-coding genes suggested that the fruiting bodies were diploid, dikaryotic or heterokaryotic. Within individual geographic populations, we found significant deviations of genotype frequencies from Hardy-Weinberg expectations, with the overall observed heterozygosity lower than that expected under random mating, consistent with prevalent inbreeding within local populations. The geographic populations were overall genetically differentiated. Interestingly, while a positive correlation was found between population genetic distance and geographic distance, there was little correlation between genetic distance and barium concentration difference for the geographic populations. Our results suggest frequent inbreeding, geographic structuring, and limited gene flow among geographic populations of T. venenata from southwestern China.     

Likelihood analysis of ongoing gene flow and historical association

Abstract.— We develop a Monte Carlo-based likelihood method for estimating migration rates and population divergence times from data at unlinked loci at which mutation rates are sufficiently low that, in the recent past, the effects of mutation can be ignored. The method is applicable to restriction fragment length polymorphisms (RFLPs) and single nucleotide polymorphisms (SNPs) sampled from a subdivided population. The method produces joint maximum-likelihood estimates of the migration rate and the time of population divergence, both scaled by population size, and provides a framework in which to test either for no ongoing gene flow or for population divergence in the distant past. We show the method performs well and provides reasonably accurate estimates of parameters even when the assumptions under which those estimates are obtained are not completely satisfied. Furthermore, we show that, provided that the number of polymorphic loci is sufficiently large, there is some power to distinguish between ongoing gene flow and historical association as causes of genetic similarity between pairs of populations.     

Population Genetics and Phylogeography of Bufo gargarizans in China

The Chinese toad Bufo gargarizans has an extensive range, covering nearly all of China. This study addresses the population genetic structure of this toad and the relation to geographic location. Partial sequences of the mtDNA control region were obtained from 14 populations, 29 haplotypes were defined, and 73 variable sites were found, including 66 parsimony informative sites. All population genetic analyses indicated no clear geographic pattern in the distribution of haplotypes. The genetic divergence between the populations was significant. Phylogeographic analyses suggested that past fragmentation and/or long-distance colonization seemed to have shaped the present-day distribution of the haplotypes of B. gargarizans.     

Genogenography of the aboriginal population of Marii El (from data on immunobiochemical polymorphism)

The geographic distribution of the frequencies of genes related to the immunological and biochemical polymorphism was studied in the Maris, who are the indigenous population of the Marii El Republic. Data on the frequencies of 33 alleles of 10 loci (ABO, TF, GC, PI, HP, AHS, F13B, ACP1, PGM1, and GLO1) in five raions (districts) of Marii El were obtained. Computer interpolation maps were constructed for all alleles. The maps allows to predict the distribution of the alleles throughout Marii El. A map of the reliability of the cartographic prediction was drawn. For the first time, the reliability of predicted gene frequencies were taken into account in constructing and interpreting the maps of gene frequencies. For the entire set of the studied genes, parameters of heterozygosity (HS) and gene diversity (GST) were estimated. Cartographic correlation analysis was performed to reveal the relationship between gene frequencies and geographic coordinates. It was found that 42% of the studied genes predominantly correlated with latitude and 9% with longitude. It was assumed that the genetic structure of Mari populations had been mainly determined by latitude-related factors. A map of Nei's genetic distances between the overall Mari gene pool and the local populations revealed a central core, which was close to the "average Mari" gene pool, and a periphery, which was genetically distant from it. Suggestions on the microevolution of the Mari gene pool were advanced. Maps of the genes with the most characteristic genetic relief (ABO*B, ACP*A, TF*D, GC*1F, PI*M2, HP*1F, and F13B*3) are shown. These maps exhibit a high correlation with the maps of principal components.     

GENETIC VARIATION IN TIME AND SPACE: MICROSATELLITE ANALYSIS OF EXTINCT AND EXTANT POPULATIONS OF ATLANTIC SALMON

Information on genetic composition of past and present populations may be obtained by analyzing DNA from archival samples. A study is presented on the genetic population structure of extant and extinct local populations of Atlantic salmon from 1913 to 1989 using dried scales as a source of DNA. Variation at six microsatellite loci was studied. Tests for differentiation among populations and among time series within populations showed that population structure was stable over time. This was also confirmed by a neighbor-joining dendrogram, which showed a clear clustering of samples from individual rivers that covered a time span of up to 76 years. These results suggest that salmon populations evolve as semi-independent units connected by modest amounts of gene flow. Additionally, a clear association between geographic and genetic distance was found. This relationship has otherwise been difficult to establish in several recent studies. The discrepancy may be due to impact of human activities on the genetic structure of present populations, whereas old samples represent populations in a more unaffected state. However, other explanations related to differences in the sampling of past and present populations may be equally valid.     

Distinguishing migration from isolation: a Markov chain Monte Carlo approach

A Markov chain Monte Carlo method for estimating the relative effects of migration and isolation on genetic diversity in a pair of populations from DNA sequence data is developed and tested using simulations. The two populations are assumed to be descended from a panmictic ancestral population at some time in the past and may (or may not) after that be connected by migration. The use of a Markov chain Monte Carlo method allows the joint estimation of multiple demographic parameters in either a Bayesian or a likelihood framework. The parameters estimated include the migration rate for each population, the time since the two populations diverged from a common ancestral population, and the relative size of each of the two current populations and of the common ancestral population. The results show that even a single nonrecombining genetic locus can provide substantial power to test the hypothesis of no ongoing migration and/or to test models of symmetric migration between the two populations. The use of the method is illustrated in an application to mitochondrial DNA sequence data from a fish species: the threespine stickleback (Gasterosteus aculeatus).     

Phylogeography of the crown-of-thorns starfish: genetic structure within the Pacific species

Previous studies of the population genetic structure of the corallivorous crown-of-thorns starfish (COTS) Acanthaster planci in the Pacific Ocean showed high levels of gene flow that were assumed to reflect a high dispersal potential. However, the phylogeographic analyses of the Pacific crown-of-thorns starfish species of this study, using the highly variable mitochondrial control region and the most complete geographic coverage to date, contradict this view. Results show high levels of overall genetic structure (Φ_ST = 0.198), suggesting a complex history of range restrictions and expansions, a pattern that we hypothesize results from changes in topography and oceanography associated with sea-level changes. However, results also show signatures of ongoing gene flow between populations isolated in the past and high levels of genetic connectivity even among distant populations. Combined, these results indicate that while there are significant limits to genetic exchange among populations among Pacific Ocean populations of the crown-of-thorns starfish, the high larval dispersal potential of this species is often achieved as well.     

线粒体基因和核基因揭示自然屏障和第四纪更新世气候振荡影响双斑乙蠊谱系地理格局

[目的]明确双斑乙蠊Sigmella biguttata种群间遗传分化程度,并揭示该种地理分布格局成因.[方法]PCR扩增双斑乙蠊19个地理种群284头个体的线粒体基因COⅠ,COⅡ和ND1以及核基因ITS的序列;使用MEGA v.7.0,DnaSP v.5.0和Arlequin v.3.5软件分析双斑乙蠊地理种群的遗...     

Patterns of genetic diversity in southern and southeastern Araucaria angustifolia (Bert.) O. Kuntze relict populations

Habitat fragmentation and a decrease in population size may lead to a loss in population genetic diversity. For the first time, the reduction in genetic diversity in the northernmost limit of natural occurence (southeastern Brazil) of Araucaria angustifolia in comparison with populations in the main area of the species continuous natural distribution (southern Brazil), was tested. The 673 AFLPs markers revealed a high level of genetic diversity for the species (Ht = 0.27), despite anthropogenic influence throughout the last century, and a decrease of H in isolated populations of southeastern Brazil (H = 0.16), thereby indicating the tendency for higher genetic diversity in remnant populations of continuous forests in southern Brazil, when compared to natural isolated populations in the southeastern region. A strong differentiation among southern and southeastern populations was detected (AMOVA variance ranged from 10%-15%). From Bayesian analysis, it is suggested that the nine populations tested form five "genetic clusters" (K = 5). Five of these populations, located in the northernmost limit of distribution of the species, represent three "genetic clusters". These results are in agreement with the pattern of geographic distribution of the studied populations.     

Isolation by distance and vicariance drive genetic structure of a coral reef fish in the Pacific Ocean

We studied the genetic diversity of a coral reef fish species to investigate the origin of the differentiation. A total of 727 Acanthurus triostegus collected from 15 locations throughout the Pacific were analyzed for 20 polymorphic loci. The genetic structure showed limited internal disequilibrium within each population; 3.7% of the loci showed significant Hardy-Weinberg disequilibrium, mostly associated with Adh*, and we subsequently removed this locus from further analysis of geographic pattern. The genetic structure of A. triostegus throughout the tropical Pacific Ocean revealed a strong geographic pattern. Overall, there was significant population differentiation (multilocus F(ST) = 0.199), which was geographically structured according to bootstraps of neighbor-joining analysis on Nei's unbiased genetic distances and AMOVA analysis. The genetic structure revealed five geographic groups in the Pacific Ocean: western Pacific (Guam, Philippines, Palau, and Great Barrier Reef); central Pacific (Solomons, New Caledonia, and Fiji); and three groups made up of the eastern populations, namely Hawaiian Archipelago (north), Marquesas (equatorial), and southern French Polynesia (south) that incorporates Clipperton Island located in the northeastern Pacific. In addition, heterozygosity values were found to be geographically structured with higher values grouped within Polynesian and Clipperton populations, which exhibited lower population size. Finally, the genetic differentiation (F(ST)) was significantly correlated with geographic distance when populations from the Hawaiian and Marquesas archipelagos were separated from all the other locations. These results show that patterns of differentiation vary within the same species according to the spatial scale, with one group probably issued from vicariance, whereas the other followed a pattern of isolation by distance. The geographic pattern for A. triostegus emphasizes the diversity of the evolutionary processes that lead to the present genetic structure with some being more influential in certain areas or according to a particular spatial scale.