Contrasting patterns of genetic diversity at three different genetic markers in a marine mammal metapopulation

Many studies use genetic markers to explore population structure and variability within species. However, only a minority use more than one type of marker and, despite increasing evidence of a link between heterozygosity and individual fitness, few ask whether diversity correlates with population trajectory. To address these issues, we analysed data from the Steller's sea lion, Eumetiopias jubatus , where three stocks are distributed over a vast geographical range and where both genetic samples and detailed demographic data have been collected from many diverse breeding colonies. To previously published mitochondrial DNA (mtDNA) and microsatellite data sets, we have added new data for amplified fragment length polymorphism (AFLP) markers, comprising 238 loci scored in 285 sea lions sampled from 23 natal rookeries. Genotypic diversity was low relative to most vertebrates, with only 37 loci (15.5%) being polymorphic. Moreover, contrasting geographical patterns of genetic diversity were found at the three markers, with Nei's gene diversity tending to be higher for AFLPs and microsatellites in rookeries of the western and Asian stocks, while the highest mtDNA values were found in the eastern stock. Overall, and despite strongly contrasting demographic histories, after applying phylogenetic correction we found little correlation between genetic diversity and either colony size or demography. In contrast, we were able to show a highly significant positive relationship between AFLP diversity and current population size across a range of pinniped species, even though equivalent analyses did not reveal significant trends for either microsatellites or mtDNA.     

Glacial refugia and the phylogeography of Steller's sea lion (Eumatopias jubatus) in the North Pacific

Mitochondrial DNA sequence data were used to examine the phylogeographic history of Steller's sea lions (Eumetopias jubatus) in relation to the presence of Plio-Pleistocene insular refugia. Cytochrome b and control region sequences from 336 Steller's sea lions reveal phylogenetic lineages associated with continental refugia south of the ice sheets in North America and Eurasia. Phylogenetic analysis suggests the genetic structure of E. jubatus is the result of Pleistocene glacial geology, which caused the elimination and subsequent reappearance of suitable rookery habitat during glacial and interglacial periods. The cyclic nature of geological change produced a series of independent population expansions, contractions and isolations that had analogous results on Steller's sea lions and other marine and terrestrial species. Our data show evidence of four glacial refugia in which populations of Steller's sea lions diverged. These events occurred from approximately 60,000 to 180,000 years BP and thus preceded the last glacial maximum.     

Microsatellite diversity, population subdivision and gene flow in the Lipizzan horse

Blood samples of 561 Lipizzan horses from subpopulations (studs) of seven European countries representing a large fraction of the breed's population were used to examine the genetic diversity, population subdivision and gene flow in the breed. DNA analysis based on 18 microsatellite loci revealed that genetic diversity (observed heterozygosity = 0.663, gene diversity = 0.675 and the mean number of alleles = 7.056) in the Lipizzan horse is similar to other horse breeds as well as to other domestic animal species. The genetic differentiation between Lipizzan horses from different studs, although moderate, was apparent (pairwise F(ST) coefficients ranged from 0.021 to 0.080). Complementary findings explaining the genetic relationship among studs were revealed by genetic distance and principal component analysis. One genetic cluster consisted of the subpopulations of Austria, Italy and Slovenia, which represent the classical pool of Lipizzan horse breeding. A second cluster was formed by the Croatian, Hungarian and Slovakian subpopulations. The Romanian subpopulation formed a separate unit. The largest genetic differentiation was found between the Romanian and Italian subpopulation. Genetic results are consistent with the known breeding history of the Lipizzan horse. Correct stud assignment was obtained for 80.9% and 92.1% of Lipizzan horses depending on the inclusion or exclusion of migrant horses, respectively. The results of the present study will be useful for the development of breeding strategies, which consider classical horse breeding as well as recent achievements of population and conservation genetics.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.     

Genetic diversity and genetic structure of populations of Ranunculus japonicus Thunb. (Ranunculaceae)

Abstract In order to clarify the genetic diversity and population structure of Ranunculus japonicus , allozymic analysis was conducted on 60 populations in southwestern Japan. Considerable genetic variati ons were detected among the populations of R. japonicus . The genetic diversities within species ( H _es = 0.215) and within populations ( H _ep = 0.172) were slightly higher than those of other perennial herbs with widespread distribution and outcrossing plants. Significantly higher values of fixation index were detected in some populations, which might have arisen from restricted mating partners. The majority of genetic variation (approx. 80%) resided within a population and a moderate level of genetic differentiation ( G _ST = 0.203) was observed among populations. The F _ST value (0.203) suggests the existence of a substantial population structure in this species. The highly significant correlation between geographic distance and F _ST values indicates that isolation by distance has played an important role in the construction of the genetic structure of this species.     

Fine-scale population structure and riverscape genetics of brook trout (Salvelinus fontinalis) distributed continuously along headwater channel networks

Linear and heterogeneous habitat makes headwater stream networks an ideal ecosystem in which to test the influence of environmental factors on spatial genetic patterns of obligatory aquatic species. We investigated fine-scale population structure and influence of stream habitat on individual-level genetic differentiation in brook trout (Salvelinus fontinalis) by genotyping eight microsatellite loci in 740 individuals in two headwater channel networks (7.7 and 4.4 km) in Connecticut, USA. A weak but statistically significant isolation-by-distance pattern was common in both sites. In the field, many tagged individuals were recaptured in the same 50-m reaches within a single field season (summer to fall). One study site was characterized with a hierarchical population structure, where seasonal barriers (natural falls of 1.5-2.5 m in height during summer base-flow condition) greatly reduced gene flow and perceptible spatial patterns emerged because of the presence of tributaries, each with a group of genetically distinguishable individuals. Genetic differentiation increased when pairs of individuals were separated by high stream gradient (steep channel slope) or warm stream temperature in this site, although the evidence of their influence was equivocal. In a second site, evidence for genetic clusters was weak at best, but genetic differentiation between individuals was positively correlated with number of tributary confluences. We concluded that the population-level movement of brook trout was limited in the study headwater stream networks, resulting in the fine-scale population structure (genetic clusters and clines) even at distances of a few kilometres, and gene flow was mitigated by 'riverscape' variables, particularly by physical barriers, waterway distance (i.e. isolation-by-distance) and the presence of tributaries.     

Natural and anthropogenic determinants of genetic structure in the largest remaining population of the endangered golden-brown mouse lemur, Microcebus ravelobensis

Genetic differentiation between natural populations is best understood as a result of both natural and anthropogenic factors. Genetic studies on large populations still living under relatively undisturbed conditions are extremely valuable to disentangle these influences. The effect of three natural (geographic distance, landscape, dispersal) factors and two anthropogenic factors (road, savannah) on gene flow was analyzed in the largest remaining forest region in the range of the endangered golden-brown mouse lemur in Madagascar. A total of 187 individuals from 12 sites were sampled and genotyped at eight polymorphic microsatellite loci. All sites exhibited similar levels of genetic variation. The level of genetic differentiation was low to moderate with pairwise F(ST) values ranging from -0.002 to 0.12, but most were significant and all sites exhibited high self-assignment rates. A spatial autocorrelation analysis was performed at two geographic scales revealing a pattern of isolation-by-distance and suggesting that no clear differences exist between male and female local dispersal. Two Bayesian approaches revealed that a stretch of savannah represented a significant barrier to movement, whereas the influence of the road on gene flow was less clear. Finally, we found that landscape characteristics, in particular altitude, play a role in the functional connectivity of the sites. The study underlines the importance of studies in relatively undisturbed conditions for the interpretation of population genetics data in fragmented environments. The results are discussed in terms of their conservation relevance for forest-dwelling animals such as most primate species.     

Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco)

A recent study of mitochondrial phylogeography of tawny owls (Strix aluco) in western Europe suggested that this species survived the Pleistocene glaciations in three allopatric refugia located in Iberia, Italy, and the Balkans, and the latter was likely the predominant source of postglacial colonization of northern Europe. New data from seven microsatellite loci from 184 individual owls distributed among 14 populations were used to assess the genetic congruence between nuclear and mitochondrial DNA (mtDNA) markers. Microsatellites corroborated the major phylogeographical conclusions reached on the basis of the mtDNA sequences, but also showed important differences leading to novel inferences. Microsatellites corroborated the three major refugia and supported the Balkan origin of northern populations. When corrected for differences in effective population size, microsatellites and mtDNA yielded generally congruent overall estimates of population structure (N*ST=0.12 vs. RST=0.16); however, there was substantial heterogeneity in the RST among the seven nuclear loci that was not correlated with heterozygosity. Populations representing the Balkans postglacial expansion interact with populations from the other two refugia forming two clines near the Alps and the Pyrenees. In both cases, the apparent position of the contact zones differed substantially between markers due to the genetic composition of populations sampled in northern Italy and Madrid. Microsatellite data did not corroborate the lower genetic diversity of northern, recently populated regions as was found with mtDNA; this discrepancy was taken as evidence for a recent bottleneck recovery. Finally, this study suggests that congruence among genetic markers should be more likely in cases of range expansion into new areas than when populations interact across contact zones.     

Dispersal and genetic structure in the American marten, Martes americana

Natal dispersal in a vagile carnivore, the American marten (Martes americana), was studied by comparing radio-tracking data and microsatellite genetic structure in two populations occupying contrasting habitats. The genetic differentiation determined among groups of individuals using F(ST) indices appeared to be weak in both landscapes, and showed no increase with geographical distance. Genetic structure investigated using pairwise genetic distances between individuals conversely showed a pattern of isolation by distance (IBD), but only in the population occurring in a homogeneous high-quality habitat, therefore showing the advantage of individual-based analyses in detecting within-population processes and local landscape effects. The telemetry study of juveniles revealed a leptokurtic distribution of dispersal distances in both populations, and estimates of the mean squared parent-offspring axial distance (sigma2) inferred both from the genetic pattern of IBD and from the radio-tracking survey showed that most juveniles make little contribution to gene flow.     

Microsatellite variation within and among recently fragmented populations of the golden lion tamarin (Leontopithecus rosalia)

Four variable microsatellite loci were used toexamine the genetic diversity and differentiation of golden lion tamarins (Leontopithecus rosalia) in four populations recently isolated by habitat fragmentation. Using Rst estimates of genetic differentiation, a considerable genetic divergence was detected among these populations, with an averagedifferentiation of 31%. Significant differences in allele number among these populations were found. However, the heterozygosity among these populations was not statistically different. These results suggestthat loss in allele diversity was faster than loss in heterozygosity. Conservation implications, particularly for golden lion tamarins, are then discussed. Loss of allelic diversity might be as serious a concern to endangered species as heterozygosity or inbreeding.     

Extremely low genetic variability and highly structured local populations of Arabidopsis thaliana at higher latitudes

The genetic diversity and population structure of Arabidopsis thaliana populations from Norway were studied and compared to a worldwide sample of A. thaliana to investigate the demographic history and elucidate possible colonization routes of populations at the northernmost species limit. We genotyped 282 individuals from 31 local populations using 149 single nucleotide polymorphism markers. A high level of population subdivision (F(ST) = 0.85 ± 0.007) was found indicating that A. thaliana is highly structured at the regional level. Significant relationships between genetic and geographical distances were found, suggesting an isolation by distance mode of evolution. Genetic diversity was much lower, and the level of linkage disequilibrium was higher in populations from the north (65-68°N) compared to populations from the south (59-62°N); this is consistent with a northward expansion pattern. A neighbour-joining tree showed that populations from northern Norway form a separate cluster, while the remaining populations are distributed over a few minor clusters. Minimal gene flow seems to have occurred between populations in different regions, especially between the geographically distant northern and southern populations. Our data suggest that northern populations represent a homogenous group that may have been established from a few founders during northward expansions, while populations in the central part of Norway constitute an admixed group established by founders of different origins, most probably as a result of human-mediated gene flow. Moreover, Norwegian populations appeared to be homogenous and isolated compared to a worldwide sample of A. thaliana, but they are still grouped with Swedish populations, which may indicate common colonization histories.     

Population genetic structure of Chilo suppressalis (Walker) (Lepidoptera: Crambidae): strong subdivision in China inferred from microsatellite markers and mtDNA gene sequences

Chilo suppressalis (Walker) displays significant geographical differences in ecological preference that may be congruent with patterns of molecular variation. To test this, we collected and analysed 381 individuals of this species from cultivated rice at 18 localities in China during the rice-growing season of 2005–2006. We used four microsatellite DNA markers and four mitochondrial DNA gene fragments. We found that this species is highly differentiated, coupled with an estimated population expansion date of at least 60 000  bp . Phylogenetic analyses, Bayesian clustering, and phylogeographical analyses of statistical parsimony haplotype network consistently divided the populations into three clades: a central China (CC) clade, a northern plus northeastern China (NN) clade and a southwestern China (SW) clade. Analysis of molecular variance indicated a high level of geographical differentiation at different hierarchical levels [ F _ST for microsatellite markers, COI, COII, 16S and ND1 is 0.06004 ( P  < 0.0001), 0.27607 ( P  < 0.0001), 0.22949 ( P  < 0.0001), 0.19485 ( P  < 0.0001) and 0.29285 ( P  < 0.0001), respectively]. Isolation by distance appeared among the samples from within China ( r  = 0.404, P  = 0.0002); N_em values estimated using a coalescent-based method were small (< 2 migrants per generation), suggesting that the observed levels of differentiation are a result of migration–drift equilibrium. Our results imply that the genetic differentiation of this borer, which is approximately in accordance with its observed number of generations per year in different Chinese geographical regions, is probably attributed to climatic and/or geological events (e.g. the last glacial maximum) and subsequently strengthened by the domestication of rice.     

Genetic survey of two endangered grasshopper subspecies, Prionotropis hystrix rhodanica and Prionotropis hystrix azami (Orthoptera, Pamphagidae): within- and between-population dynamics at the regional scale

The fragmentation and destruction of natural habitats is one of the major causes of extinction of wild flora and fauna due to human activities. The area around the Mediterranean Basin is particularly affected: this region is one of the world’s richest in terms of biodiversity and is undergoing substantial modification due to extensive changes in land use. We investigated the demographic state (through a genetic survey) of an endangered grasshopper, the geographical distribution of which is restricted to arid open land or glades in south-eastern France. We compared two subspecies: one is severely threatened by the fragmentation of its specific and restricted habitat (Prionotropis hystrix rhodanica in the ‘Crau’ steppe) whereas the other, P.␣h.␣azami, occurs more widely in glades or clearings in the evergreen Mediterranean forest. Genetic diversity at eight microsatellite loci was highly structured, indicating substantial isolation of populations. Within-population diversity was high at all but one site. Genetic drift was the major force involved in this genetic structure, with very little gene flow at the regional scale of observation, consistent with both the limited dispersal of this flightless species and the patchy configuration of its habitat. No significant differences in the extent of genetic diversity or population dynamics were observed between the two subspecies. The observed genetic pattern calls into question the validity of the subspecies status, and shows no strong impact of recent anthropogenic modification of the Crau steppe. The evolution of the Mediterranean xeric and open habitats under both natural and anthropogenic processes is discussed.     

Genetic population structure and contemporary dispersal patterns of a recent European invader, the Chinese mitten crab, Eriocheir sinensis

Genetic studies of recently established populations are challenging because the assumption of equilibrium underlying many analyses is likely to be violated. Using microsatellites, we investigated determinants of genetic structure and migration among invasive European-Chinese mitten crab populations, applying a combination of traditional population genetic analyses and nonequilibrium Bayesian methods. Consistent with their recent history, invasive populations showed much lower levels of genetic diversity than a native Chinese population, indicative of recent bottlenecks. Population differentiation was generally low but significant and especially pronounced among recently established populations. Significant differentiation among cohorts from the same geographical location (River Thames) suggests the low effective population size and associated strong genetic drift that would be anticipated from a very recent colonization. An isolation-by-distance pattern appears to be driven by an underlying correlation between geographical distance and population age, suggesting that cumulative homogenizing gene flow reduces founder bottleneck-associated genetic differentiation between longer-established populations. This hypothesis was supported by a coalescent analysis, which supported a drift + gene flow model as more likely than a model excluding gene flow. Furthermore, admixture analysis identified several recent migrants between the UK and Continental European population clusters. Admixture proportions were significantly predicted by the volume of shipping between sites, indicating that human-mediated transport remains a significant factor for dispersal of mitten crabs after the initial establishment of populations. Our study highlights the value of nonequilibrium methods for the study of invasive species, and also the importance of evaluating nonequilibrium explanations for isolation by distance patterns.     

Compatible genetic and ecological estimates of dispersal rates in insect (Coenagrion mercuriale: Odonata: Zygoptera) populations: analysis of 'neighbourhood size' using a more precise estimator

Genetic and demographic estimates of dispersal are often thought to be inconsistent. In this study, we use the damselfly Coenagrion mercuriale (Odonata: Zygoptera) as a model to evaluate directly the relationship between estimates of dispersal rate measured during capture-mark-recapture fieldwork with those made from the spatial pattern of genetic markers in linear and two-dimensional habitats. We estimate the 'neighbourhood size' (Nb) - the product of the mean axial dispersal rate between parent and offspring and the population density - by a previously described technique, here called the regression method. Because C. mercuriale is less philopatric than species investigated previously by the regression method we evaluate a refined estimator that may be more applicable for relatively mobile species. Results from simulations and empirical data sets reveal that the new estimator performs better under most situations, except when dispersal is very localized relative to population density. Analysis of the C. mercuriale data extends previous results which demonstrated that demographic and genetic estimates of Nb by the regression method are equivalent to within a factor of two at local scales where genetic estimates are less affected by habitat heterogeneity, stochastic processes and/or differential selective regimes. The corollary is that with a little insight into a species' ecology the pattern of spatial genetic structure provides quantitative information on dispersal rates and/or population densities that has real value for conservation management.     

Extensive gene flow over Europe and possible speciation over Eurasia in the ectomycorrhizal basidiomycete Laccaria amethystina complex

Biogeographical patterns and large-scale genetic structure have been little studied in ectomycorrhizal (EM) fungi, despite the ecological and economic importance of EM symbioses. We coupled population genetics and phylogenetic approaches to understand spatial structure in fungal populations on a continental scale. Using nine microsatellite markers, we characterized gene flow among 16 populations of the widespread EM basidiomycete Laccaria amethystina over Europe (i.e. over 2900 km). We also widened our scope to two additional populations from Japan (10(4) km away) and compared them with European populations through microsatellite markers and multilocus phylogenies, using three nuclear genes (NAR, G6PD and ribosomal DNA) and two mitochondrial ribosomal genes. European L. amethystina populations displayed limited differentiation (average F(ST) = 0.041) and very weak isolation by distance (IBD). This panmictic European pattern may result from effective aerial dispersal of spores, high genetic diversity in populations and mutualistic interactions with multiple hosts that all facilitate migration. The multilocus phylogeny based on nuclear genes confirmed that Japanese and European specimens were closely related but clustered on a geographical basis. By using microsatellite markers, we found that Japanese populations were strongly differentiated from the European populations (F(ST) = 0.416), more than expected by extrapolating the European pattern of IBD. Population structure analyses clearly separated the populations into two clusters, i.e. European and Japanese clusters. We discuss the possibility of IBD in a continuous population (considering some evidence for a ring species over the Northern Hemisphere) vs. an allopatric speciation over Eurasia, making L. amethystina a promising model of intercontinental species for future studies.     

Microgeographic genetic structure and gene flow in Hibiscus moscheutos (Malvaceae) populations

Microgeographic genetic variation in populations of a wetland macrophyte, Hibiscus moscheutos L. (Malvaceae), was investigated using allozyme polymorphism. The species is a self-compatible insect-pollinated perennial, and seeds are water dispersed (hydrochory). Six hundred plants were analyzed from eight brackish and two freshwater populations within the Rhode River watershed/estuarine system. The genetic structure of the populations was assessed by fixation indices and spatial autocorrelation analyses. The degree of genetic differentiation among sites and gene flow between all paired combinations of sites (M ) was analyzed using three hypothetical gene flow models. Fixation indices indicated almost complete panmixia within populations, and spatial autocorrelations showed that genotypes were randomly distributed within sites, most likely the result of water dispersal of seeds. Allele frequencies were significantly different among sites, and estimated FST indicated moderate genetic differentiation (_ = 0.062). Genetic differences between populations were mostly explained by a gene flow model that accounted for the location of populations relative to the tidal stream. The importance of hydrochory in affecting spatial genetic structure was thus suggested both within and among H. moscheutos populations.