ISOLATION BY DISTANCE

A novel theory to explain species diversity in landscapes: positive frequency dependence and habitat suitability

Theories to explain the diversity of species have required that individual species occupy unique niches and/or vary in their response to environmental factors. Positive interactions within a species, although common in communities, have not been thought to maintain species diversity because in non-spatial models the more abundant species always outcompetes the rarer species. Here, we show, using a stochastic spatial model, that positive intraspecific interactions such as those caused by positive frequency dependence and/or priority effects, can maintain species diversity if interactions between individuals are primarily local and the habitat contains areas that cannot be colonized by any species, such as boulders or other physical obstructions. When intraspecific interactions are primarily neutral, species diversity will eventually erode to a single species. When the landscape is homogeneous (i.e. does not contain areas that cannot be colonized by any species), the presence of strong intraspecific interactions will not maintain diversity.

A simulated annealing approach to define the genetic structure of populations

We present a new approach for defining groups of populations that are geographically homogeneous and maximally differentiated from each other. As a by-product, it also leads to the identification of genetic barriers between these groups. The method is based on a simulated annealing procedure that aims to maximize the proportion of total genetic variance due to differences between groups of populations (spatial analysis of molecular variance; samova). Monte Carlo simulations were used to study the performance of our approach and, for comparison, the behaviour of the Monmonier algorithm, a procedure commonly used to identify zones of sharp genetic changes in a geographical area. Simulations showed that the samova algorithm indeed finds maximally differentiated groups, which do not always correspond to the simulated group structure in the presence of isolation by distance, especially when data from a single locus are available. In this case, the Monmonier algorithm seems slightly better at finding predefined genetic barriers, but can often lead to the definition of groups of populations not differentiated genetically. The samova algorithm was then applied to a set of European roe deer populations examined for their mitochondrial DNA (mtDNA) HVRI diversity. The inferred genetic structure seemed to confirm the hypothesis that some Italian populations were recently reintroduced from a Balkanic stock, as well as the differentiation of groups of populations possibly due to the postglacial recolonization of Europe or the action of a specific barrier to gene flow.

Comparing relative rates of pollen and seed gene flow in the island model using nuclear and organelle measures of population structure

We describe a method for comparing nuclear and organelle population differentiation (F(ST)) in seed plants to test the hypothesis that pollen and seed gene flow rates are equal. Wright's infinite island model is used, with arbitrary levels of self-fertilization and biparental organelle inheritance. The comparison can also be applied to gene flow in animals. Since effective population sizes are smaller for organelle genomes than for nuclear genomes and organelles are often uniparentally inherited, organelle F(ST) is expected to be higher at equilibrium than nuclear F(ST) even if pollen and seed gene flow rates are equal. To reject the null hypothesis of equal seed and pollen gene flow rates, nuclear and organelle F(ST)'s must differ significantly from their expected values under this hypothesis. Finite island model simulations indicate that infinite island model expectations are not greatly biased by finite numbers of populations (>/=100 subpopulations). The power to distinguish dissimilar rates of pollen and seed gene flow depends on confidence intervals for fixation index estimates, which shrink as more subpopulations and loci are sampled. Using data from the tropical tree Corythophora alta, we rejected the null hypothesis that seed and pollen gene flow rates are equal but cannot reject the alternative hypothesis that pollen gene flow is 200 times greater than seed gene flow.

Slow mitochondrial DNA sequence evolution in the Anthozoa (Cnidaria)

Mitochondrial genes have been used extensively in population genetic and phylogeographical analyses, in part due to a high rate of nucleotide substitution in animal mitochondrial DNA (mtDNA). Nucleotide sequences of anthozoan mitochondrial genes, however, are virtually invariant among conspecifics, even at third codon positions of protein-coding sequences. Hence, mtDNA markers are of limited use for population-level studies in these organisms. Mitochondrial gene sequence divergence among anthozoan species is also low relative to that exhibited in other animals, although higher level relationships can be resolved with these markers. Substitution rates in anthozoan nuclear genes are much higher than in mitochondrial genes, whereas nuclear genes in other metazoans usually evolve more slowly than, or similar to, mitochondrial genes. Although several mechanisms accounting for a slow rate of sequence evolution have been proposed, there is not yet a definitive explanation for this observation. Slow evolution and unique characteristics may be common in primitive metazoans, suggesting that patterns of mtDNA evolution in these organisms differ from that in other animal systems.

Allozyme diversity in natural populations of Viola palmensis Webb & Berth. (Violaceae) from La Palma (Canary Islands): implications for conservation genetics

Genetic diversity was measured by allozyme electrophoresis in eight natural populations of the threatened Canarian endemic Viola palmensis Webb & Berth. (Violaceae). Nineteen alleles corresponding to 11 gene loci were detected. High levels of genetic diversity were found, ranging from 36.3 to 45.4 % for the percentage of polymorphic loci (P), from 1.45 to 1.60 for the average number of alleles per locus (A) and from 0.128 to 0.200 for the expected heterozygosity (H(e)). Between 85.5 and 96.6 % of genetic variability was apportioned within populations. As a whole, populations were not at Hardy-Weinberg equilibrium, with a deficit of heterozygous individuals attributable to the existence of genetic structuring in the populations analysed. The levels of interpopulation genetic differentiation were low (mean F(ST) = 0.100), while genetic identity pair-wise comparisons were high (mean I = 0.973) suggesting considerable levels of gene flow among populations. No relationship was detected between genetic differentiation and geographical distances between populations. An outcrossing insect-mediated breeding system might contribute to pollen dispersion of this species. For conservation genetics we suggest in situ preservation areas are defined that are free of disturbance and that include populations with the highest genetic diversity.

Phylogeographic breaks without geographic barriers to gene flow

The spatial distribution of genetic markers can be useful both in estimating patterns of gene flow and in reconstructing biogeographic history, particularly when gene genealogies can be estimated. Genealogies based on nonrecombining genetic units such as mitochondrial and chloroplast DNA often consist of geographically separated clades that come into contact in narrow regions. Such phylogeographic breaks are usually assumed to be the result of long-term barriers to gene flow. Here I show that deep phylogeographic breaks can form within a continuously distributed species even when there are no barriers to gene flow. The likelihood of observing phylogeographic breaks increases as the average individual dispersal distance and population size decrease. Those molecular markers that are most likely to show evidence of real geographic barriers are also most likely to show phylogeographic breaks that formed without any barrier to gene flow. These results might provide an explanation as to why some species, such as the greenish warblers (Phylloscopus trochiloides), have phylogeographic breaks in mitochondrial or chloroplast DNA that do not coincide with sudden changes in other traits.

Historical vicariance and postglacial colonization effects on the evolution of genetic structure in Lophocereus, a Sonoran Desert columnar cactus

Distinguishing the historical effects of gene migration and vicariance on contemporary genetic structure is problematic without testable biogeographic hypotheses based on preexisting geological and environmental evidence. The availability of such hypotheses for North America's Sonoran Desert has contributed to our understanding of the effect of historical vicariance and dispersal events on the diversification of this region's vertebrate biota but have not yet been applied to its flora. In this paper we describe a detailed allozyme analysis of the population genetic structure and phylogeography of the Sonoran Desert columnar cactus, Lophocereus schottii (senita). Inferred phylogroup distributions reflect two historical vicariance events: (1) a middle Pliocene northward transgression of the Sea of Cortéz that is reflected in well-supported Baja California peninsular and continental phylogroups but not in current taxonomic treatments of the species; and (2) a late Pliocene transpeninsular seaway across southern Baja that is reflected in tentative support for peninsular and southern Cape Region phylogroups corresponding to taxonomic varieties L. schottii var. schottii and L. schottii var. australis, respectively. A middle Pleistocene midpeninsular seaway hypothesized to explain congruent phylogroup distributions in several vertebrate taxa is not reflected in L. s. var. schottii, nor is the distinction of a third variety, L. s. var. tenuis, from continental populations of L. s. var. schottii. Linear regression of pairwise estimates of interpopulation differentiation (M and F(ST)/[1 - F(ST)]) on interpopulation geographic distance revealed significant evidence of isolation by distance within peninsular and continental phylogroups but not between them, consistent with historical vicariance between but not within these regions. We also found significant evidence of isolation by distance between putative L. s. var. schottii and L. s. var. australis phylogroups, suggesting that reproductive isolation between peninsular and Cape Region forms is incomplete. Within peninsular, but not continental, phylogroups, northward range expansion from southern Pleistocene refugia is reflected in significant declines in genetic variation with increasing latitude and in an area phenogram in which populations are progressively nested from south (ancestral) to north (descendant) along the Baja peninsula. Although the geographic concordance of phylogenetic topologies suggests that ancient vicariance events, and not dispersal, have primarily influenced the biogeographic distributions of Baja's vertebrate biota, the phylogeographic structure of L. schottii suggests that Sonoran Desert plant species may exhibit genetic signatures of postglacial range expansion and gene flow as well as vicariance.

Population subdivision and gene flow among wild orangutans

Genetic variability among populations of orangutans from Borneo and Sumatra was assessed using seven SSR loci. Most SSR loci were highly polymorphic and their allele frequencies exhibited substantial variation across subpopulations. While significant genetic subdivision was observed among the island populations, genetic distance did not increase with geographic distance and sufficient gene flow persists to prevent marked genetic subdivision. Since it is unlikely that the Bornean Orangutans dispersed naturally among locations separated by such formidable geographic barriers, human assistance might already have altered their genetic structure. Our data suggests that there may be at least two subspecific clades of orangutans within Borneo while Central Kalimantan animals may have become more genetically related to animals in Sumatra due to human intervention.

Fine-scale genetic structure in populations ofQuercus variabilis(Fagaceae) from southern Korea

Quercus variabilis Blume (Fagaceae) is a deciduous broad-leaved tree, and an important forest element among the hillsides of southern Korea. To date, there are contrasting results with respect to fine-scale spatial genetic structure among adults in populations of several oak species; some studies have shown evidence of significant within-population spatial genetic structure, while others found weak or little evidence of fine-scale genetic structuring within populations. We used allozyme loci, Wright's F statistics, and multilocus spatial autocorrelation statistics to examine the distribution of genetic diversity and the spatial genetic structure within three undisturbed populations at a landscape level (72.5 ha, 500 × 1450 m) on Dolsan Island, South Korea. The spatial autocorrelation analyses of adults showed little evidence of fine-scale genetic structuring within populations, which could be due to random mortality among related seedlings, resulting in extensive thinning within maternal half-sib groups. Alternatively, low genetic differentiation between adjacent populations (mean FST= 0.023) and little within-population spatial genetic structure suggest probable secondary acorn movement by animals. Our results are very similar to those observed in three populations of Quercus acutissima at a landscape level (15 ha, 250 × 600 m) on Oenaro Island, South Korea. Together, these studies describe relatively subtle differences in genetic structure among adjacent populations of oaks on southern islands in Korea.Key words: allozymes, Quercus variabilis, Fagaceae, landscape level, multiple populations, spatial genetic structure.

A spatially explicit neutral model of beta-diversity in tropical forests

To represent species turnover in tropical rain forest, we use a neutral model where a tree's fate is not affected by what species it belongs to, seeds disperse a limited distance from their parents, and speciation is in equilibrium with random extinction. We calculate the similarity function, the probability F(r) that two trees separated by a distance r belong to the same species, assuming that the dispersal kernel P(r), the distribution of seeds about their parents and the prospects of mortality and reproduction, are the same for all trees regardless of their species. If P(r) is radially symmetric Gaussian with mean-square dispersal distance sigma, F(r) can be expressed in closed form. If P(r) is a radially symmetric Cauchy distribution, then, in two-dimensional space, F(r) is proportional to 1/r for large r. Analytical results are compared with individual-based simulations, and the relevance to field observations is discussed.

Effective Population Size and Population Subdivision in Demographically Structured Populations

A fast-timescale approximation is applied to the coalescent process in a single population, which is demographically structured by sex and/or age. This provides a general expression for the probability that a pair of alleles sampled from the population coalesce in the previous time interval. The effective population size is defined as the reciprocal of twice the product of generation time and the coalescence probability. Biologically explicit formulas for effective population size with discrete generations and separate sexes are derived for a variety of different modes of inheritance. The method is also applied to a nuclear gene in a population of partially self-fertilizing hermaphrodites. The effects of population subdivision on a demographically structured population are analyzed, using a matrix of net rates of movement of genes between different local populations. This involves weighting the migration probabilities of individuals of a given age/sex class by the contribution of this class to the leading left eigenvector of the matrix describing the movements of genes between age/sex classes. The effects of sex-specific migration and nonrandom distributions of offspring number on levels of genetic variability and among-population differentiation are described for different modes of inheritance in an island model. Data on DNA sequence variability in human and plant populations are discussed in the light of the results.

Population genetic structure of the cyclic snowshoe hare (Lepus americanus) in southwestern Yukon, Canada

Spatial population structure has important ecological and evolutionary consequences. Little is known about the population structure of snowshoe hares (Lepus americanus), despite their ecological importance in North American boreal forests. We used seven variable microsatellite DNA loci to determine the spatial genetic structure of snowshoe hares near Kluane Lake, Yukon during a cyclic population peak. We sampled 317 hares at 12 sites separated by distances ranging from 3 to 140 km, and used 46 additional samples from Alaska and Montana. The level of genetic variation was high (13.4 alleles/locus, 0.67 expected heterozygosity) and the distribution of alleles and genotypes was not homogeneous across the sites. The degree of differentiation was low among Yukon sites (FST = 0.015) and between Yukon and Alaska (FST = 0.012), but the Montana site was highly differentiated (FST = 0.20). A weak pattern of isolation by distance was found over the Yukon study area, with an indication that local genetic drift may be important in shaping the regional genetic structure. Landscape barriers expected to influence gene flow did not consistently affect genetic structure, although there was evidence for a partial barrier effect of Kluane Lake. The high level of inferred gene flow confirms that snowshoe hare dispersal is widespread, successful and equal between the sexes. A stepping-stone model of gene flow, potentially influenced by the synchronous density cycle, appears to best explain the observed genetic structure. Our results suggest that despite their dramatic fluctuations in density, snowshoe hares in the northern boreal forest have a large evolutionary effective population size and are not strongly subdivided by either physical or social barriers to gene flow.

Inferring the phylogeny of disjunct populations of the azure-winged magpie Cyanopica cyanus from mitochondrial control region sequences

The azure-winged magpie (AWM), Cyanopica cyanus, is found in Asia and Iberia. This remarkable disjunct distribution has been variously explained by either the sixteenth-century introduction of birds into Iberia from the Far East, or by the loss of individuals from the central part of their range as a result of Pleistocene glaciations. We have used the mitochondrial control region to undertake a molecular phylogenetic analysis of the AWM, with sequences examined from individuals collected from across the current distribution range and incorporating representatives of all currently defined subspecies. The Western birds are genetically distinct from their Asian congeners and their divergence is basal in the phylogenetic tree. This indicates that the AWM is native to Iberia and not the result of a recent introduction from Asia. In Asia, two major mitochondrial DNA lineages were identified. These correspond to an Inland Asia group and a Pacific Seaboard group, and are separated topographically by the Da Hingan Ling mountains and the Yellow Sea. Molecular clock estimates suggest that these divergences are associated with Pleistocene glaciations. Furthermore, our data do not support the current classification of the AWM into 10 subspecies, as defined based on morphology and geographical distribution.

Geographical and altitudinal population genetic structure of two dung fly species with contrasting mobility and temperature preference

Local adaptation of populations requires some degree of spatio-temporal isolation. Previous studies of the two dung fly species Scathophaga stercoraria and Sepsis cynipsea have revealed low levels of geographic and altitudinal genetic differentiation in quantitative life history and morphological traits, but instead high degrees of phenotypic plasticity. These patterns suggest that gene flow is extensive despite considerable geographic barriers and large spatio-temporal variation in selection on body size and related traits. In this study we addressed this hypothesis by investigating genetic differentiation of dung fly populations throughout Switzerland based on the same 10 electrophoretic loci in each species. Overall, we found no significant geographic differentiation of populations for either species. This is inconsistent with the higher rates of gene flow expected due to better flying capacity of the larger S. stercoraria. However, heterozygote deficiencies within populations indicated structuring on a finer scale, seen for several loci in S. cynipsea, and for the locus PGM (Phosphoglucomutase) in S. stercoraria. Additionally, S. cynipsea showed a tendency towards a greater gene diversity at higher altitudes, mediated primarily by the locus MDH (malate dehydrogenase), at which a second allele was only present in populations above 1000 m. This may be caused by increased environmental stress at higher altitudes in this warm-adapted species. MDH might thus be a candidate locus subject to thermal selection in this species, but this remains to be corroborated by direct evidence. In S. stercoraria, no altitudinal variation was found.

A modified stepping-stone model of population structure in red drum, Sciaenops ocellatus (Sciaenidae), from the northern Gulf of Mexico

Genetic studies of population or 'stock' structure in exploited marine fishes typically are designed to determine whether geographic boundaries useful for conservation and management planning are identifiable. Implicit in many such studies is the notion that subpopulations or stocks, if they exist, have fixed territories with little or no gene exchange between them. Herein, we review our long-term genetic studies of red drum (Sciaenops ocellatus), an estuarine-dependent sciaenid fish in the Gulf of Mexico and western Atlantic Ocean. Significant differences in frequencies of mitochondrial DNA haplotypes and of alleles at nuclear-encoded microsatellites occur among red drum sampled across the northern Gulf of Mexico. The spatial distribution of the genetic variation, however, follows a pattern of isolation-by-distance consistent with the hypothesis that gene flow occurs among subpopulations and is an inverse (and continuous) function of geographic distance. However, successful reproduction and recruitment of red drum depend on estuarine habitats that have geographically discrete boundaries. We hypothesize that population structure in red drum follows a modified one-dimensional, linear stepping-stone model where gene exchange occurs primarily (but not exclusively) between adjacent bays and estuaries distributed linearly along the coastline. Gene flow does occur among estuaries that are not adjacent but probabilities of gene exchange decrease as a function of geographic distance. Implications of our hypothesis are discussed in terms of inferences drawn from patterns of isolation-by-distance and relative to conservation and management of estuarine-dependent species like red drum. Based on estimates of the ratio of genetic effective population size and census size in red drum, observed patterns of gene flow in red drum may play a significant role in recruitment.

Implications of movement patterns for gene flow in black rat snakes (Elaphe obsoleta)

Gene flow is fundamental to evolutionary processes but knowledge about movements of individuals and their offspring necessary for gene flow is scant. We investigated potential ecological components of genetic connectivity within a population of black rat snakes (Elaphe obsoleta) by radio-tracking 82 individuals. Because adult black rat snakes are highly faithful to their hibernaculum, gene flow between hibernaculum populations has to occur through mating between members of different hibernacula or through juvenile dispersal. The present study was the first to assess the spatial dispersion of a complete network of hibernacula. The mean distance between the nearest-neighbour hibernacula was 811 m, which was less than the mean distances that reproductive males and females were found from their hibernacula during the mating season. Estimates of maximum distances individuals were from their hibernacula during the mating season indicated that, on average, a female was likely to mate with males that came from two hibernacula away from the female's own hibernaculum. Both males and females appeared to contribute actively to gene flow by moving more and increasing their distance from their hibernacula during the mating season. In addition, on average, females nested closer to a hibernaculum other than the one they attended, thereby potentially increasing the likelihood that their offspring would join hibernacula other than their mothers'. Thus, spatial and movement patterns of male and female black rat snakes are consistent with genetic evidence of extensive out-breeding among local hibernaculum populations.

Genetic differentiation within and between populations of a hermaphroditic freshwater planarian

Dispersal of individuals is an important factor that can influence genetic differentiation between populations. The hermaphroditic freshwater planarian Schmidtea polychroainhabits shallow regions of lakes and streams, in which they appear to be continuously distributed. In the present study we used three highly polymorphic markers for analysing small-scale and large-scale genetic structure within one, and between four natural lake populations. Genetic differentiation could already be observed between samples collected at least 13 m apart, but not between neighbouring samples, and was most pronounced between samples from different lakes. Probably due to the high variance in F(ST)values, a significant correlation between genetic differentiation and geographic distance could not be observed. These results show that individual dispersal of S. polychroa is limited, but that there is gene flow between subpopulations from the same lake. They further suggest that long-distance dispersal and gene flow between lakes, if present, is not a common process in S. polychroa.

Phenotypic and allozyme variation in Mediterranean and desert populations of wild barley, Hordeum spontaneum Koch

Populations of wild barley, Hordeum spontaneum Koch, were collected in two distinct climatic regions, desert and Mediterranean. Plants from five desert and five Mediterranean populations were compared and contrasted for extent and structure of phenotypic variation. These same 10 and one other population from each region were analyzed for allozyme variation. In a field trial of phenotypic diversity, two phenological and 14 morphological traits were examined. Study of allozyme variation was performed using eight enzyme systems encoding for 13 loci. Plants from the desert and Mediterranean regions were significantly different in seven of 16 phenotypic traits, exhibited a high (30%) interregional component of phenotypic variation, and showed a high degree of segregation on a principal component scattergram indicating ecotypic differentiation. Mediterranean populations were twice as variable as desert populations in reproductive growth parameters (stem and spike length) and grain filling (spikelet weight), but half as variable for onset of reproduction. The extent and structure of phenotypic and allozyme variation did not match. The Mediterranean and desert populations did not differ in amount of allozyme variation as estimated by mean number of alleles per locus, effective number of alleles, polymorphism, and gene diversity (n(a), n(e), P, and H(e)), did not segregate on the basis of population genetic distances, and exhibited a low proportion of interregion allozyme diversity (2%). No effect of selection on allozyme distribution was detected. Our results suggest that the adaptation of plants originating from desert and Mediterranean environments is reflected in phenotypic but not in allozyme variation.

Population structure of a lodgepole pine (Pinus contorta) and jack pine (P. banksiana) complex as revealed by random amplified polymorphic DNA

We studied the population structure of a lodgepole (Pinus contorta Dougl.) and jack pine (Pinus banksiana Lamb.) complex in west central Alberta and neighboring areas by assessing random amplified polymorphic DNA (RAPD) variability in 23 lodgepole pine, 9 jack pine, and 8 putative hybrid populations. Of 200 random primers screened, 10 that amplified 39 sharp and reproducible RAPDs were chosen for the study. None of the 39 RAPDs were unique to the parental species. RAPD diversity ranged from 0.085 to 0.190 among populations and averaged 0.143 for lodgepole pine, 0.156 for jack pine, 0.152 for hybrids, and 0.148 for all 40 populations. The estimated population differentiation based on G(ST) was 0.168 for hybrids, 0.162 for lodgepole pine, 0.155 for jack pine, and 0.247 across all 40 populations. Cluster analysis of genetic distances generally separated jack pine from lodgepole pine and hybrids, but no division could be identified that further separated lodgepole pine from hybrids. The observed weak to mild trend of "introgression by distance" in the complex and neighbouring areas was consistent with the view that introgressive hybridization between lodgepole and jack pines within and outside the hybrid zone may have been through secondary contact and primary intergradation, respectively.

The coalescent in a continuous, finite, linear population

In this article we present a model for analyzing patterns of genetic diversity in a continuous, finite, linear habitat with restricted gene flow. The distribution of coalescent times and locations is derived for a pair of sequences sampled from arbitrary locations along the habitat. The results for mean time to coalescence are compared to simulated data. As expected, mean time to common ancestry increases with the distance separating the two sequences. Additionally, this mean time is greater near the center of the habitat than near the ends. In the distant past, lineages that have not undergone coalescence are more likely to have been at opposite ends of the population range, whereas coalescent events in the distant past are biased toward the center. All of these effects are more pronounced when gene flow is more limited. The pattern of pairwise nucleotide differences predicted by the model is compared to data collected from sardine populations. The sardine data are used to illustrate how demographic parameters can be estimated using the model.

Comparison of genetic diversity estimates within and among populations of maritime pine using chloroplast simple-sequence repeat and amplified fragment length polymorphism data

We compared the genetic variation of Pinus pinaster populations using amplified fragment length polymorphism (AFLP) and chloroplast simple-sequence repeat (cpSSR) loci. Populations' levels of diversity within groups were found to be similar with AFLPs, but not with cpSSRs. The high interlocus variance associated with the AFLP loci could account for the lack of differences in the former. Although AFLPs revealed much lower genetic diversity than cpSSRs, the levels of among-population differentiation found with the two types of marker were similar, provided that loci showing fewer than four null-homozygotes, in any population, were pruned from the AFLP data. Moreover, the French and Portuguese populations were clearly differentiated from each other, with both markers. The Mantel test showed that the genetic distance matrix calculated using the AFLP data was correlated with the matrix derived from the cpSSRs. Because of the concordance found between markers we conclude that gene flow was indeed the predominant force shaping nuclear and chloroplastic genetic variation of the populations within regions, at the geographical scale studied.

Pollen production, microsporangium dehiscence and pollen flow in Himalayan cedar (Cedrus deodara Roxb. ex d. Don)

Microsporangium dehiscence, pollen production and dispersal were studied in Himalayan cedar (Cedrus deodara) during 1998 and 1999. Microsporangium dehiscence showed diurnal periodicity and was found to be related to air temperature and relative air humidity, with a strobilus taking 2 d to dehisce completely in warmer conditions and 3 d in cooler ones. The frequency of flowering in C. deodara was highly variable during the two successive years; however, cyclical production of pollen grains was observed in 50% of the trees. The maximum concentration of pollen grains in the air was found between 1200 and 1600 h, and this period was also noted to be the best time for pollination. Studying migration of pollen grains from isolated single trees in three directions showed that migration was not uniform in all directions. Long-distance transport of pollen grains was observed in the downhill direction. However, in the uphill and horizontal directions grains could travel only up to 97.5 and l95.1 m, respectively, and the frequency of pollen grains to the source frequency at these distances was only 1.9 and 2-5%, respectively. The results suggest that an isolation barrier of 190 m may be considered as a minimum for the management of deodar seed orchards.

The genetic structure of Escherichia coli populations in primary and secondary habitats

Escherichia coli were recovered from the members of two two-person households and their associated septic tanks. The E. coli were isolated using selective and non-selective isolation techniques, characterized using the method of multi-locus enzyme electrophoresis and screened for the presence of virulence factors associated with extra-intestinal disease by using PCR. The growth rate-temperature relationships of strains from the two habitats were also determined. Temporal variation explained 25% of the observed electrophoretic type (ET) diversity in the humans. Among-host variation accounted for 29% of the observed allelic diversity. In one household, ET diversity of the E. coli population in the septic tank was significantly lower than ET diversity in the humans providing the inputs to the septic tank. Molecular analysis of variance revealed that, on average, strains recovered from the septic tank of this household were genetically distinct from strains recovered from the humans providing the faecal inputs to the septic tank. Further, the growth rate-temperature response of strains differed between strains isolated from the septic tank and strains isolated from the humans. Septic tank isolates grew better at low temperatures than strains isolated from humans, but more slowly at high temperatures compared to the human isolates. By contrast, no real differences in ET diversity, allelic diversity, or the growth characteristics of strains could be detected between strains from the humans and strains from the septic tank of the other household. The results of this study suggest there are strains of E. coli that are better "adapted" to conditions found in the external environment compared to strains isolated from the gastrointestinal habitat. Further, the finding that the numerically dominant clones and clonal diversity in secondary habitats can differ substantially from those found in the source populations will confound efforts to identify the sources of faecal pollution in the environment.

Origin of Fueguian-Patagonians: an approach to population history and structure using R matrix and matrix permutation methods

A complicated history of isolation between Fueguian and Patagonian groups (originated by the appearance of the Straits of Magellan) as much as differences in population structure and life strategies constitute important factors in the clustering pattern of those groups. The aim of this work was to test several hypotheses about population structure and history of Fueguian-Patagonians to propose a model that incorporates predictions for future studies. R matrix methods and matrix permutation analyses were performed upon a data matrix of craniofacial measurements of 441 skulls divided into nine samples pertaining to six Patagonian and three Fueguian populations. Association of biological distances with three matrices representing several settlement patterns was tested using matrix permutation tests. Results of R matrix study show that the minimum genetic distance obtained confirms separation between Fueguians and Patagonians. Moreover, an analysis of residual variances from the expected regression line confirms admixture between Andean and Pampean populations and Araucanian groups, consistent with ethnohistorical observations. A model representing a long history of isolation between Fueguian and Patagonians, rather than a model emphasizing differences in life-strategies, presented the best correlation with the biological distance matrix. Because similar results were already obtained in archaeological, molecular, and morphological studies, a model for the settlement of Tierra del Fuego is proposed. It is summarized by four main hypotheses that can be tested independently by different disciplines in the future.

Inbreeding and relatedness coefficients: what do they measure?

This paper reviews and discusses what is known about the relationship between identity in state, allele frequency, inbreeding coefficients, and identity by descent in various uses of these terms. Generic definitions of inbreeding coefficients are given, as ratios of differences of probabilities of identity in state. Then some of their properties are derived from an assumption in terms of differences between distributions of coalescence times of different genes. These inbreeding coefficients give an approximate measurement of how much higher the probability of recent coalescence is for some pair of genes relative to another pair. Such a measure is in general not equivalent to identity by descent; rather, it approximates a ratio of differences of probabilities of identity by descent. These results are contrasted with some other formulas relating identity, allele frequency, and inbreeding coefficients. Additional assumptions are necessary to obtain most of them, and some of these assumptions are not always correct, for example when there is localized dispersal. Therefore, definitions based on such formulas are not always well-formulated. By contrast, the generic definitions are both well-formulated and more broadly applicable.

Two-generation analysis of pollen flow across a landscape. IV. Estimating the dispersal parameter

The distance of pollen movement is an important determinant of the neighborhood area of plant populations. In earlier studies, we designed a method for estimating the distance of pollen dispersal, on the basis of the analysis of the differentiation among the pollen clouds of a sample of females, spaced across the landscape. The method was based solely on an estimate of the global level of differentiation among the pollen clouds of the total array of sampled females. Here, we develop novel estimators, on the basis of the divergence of pollen clouds for all pairs of females, assuming that an independent estimate of adult population density is available. A simulation study shows that the estimators are all slightly biased, but that most have enough precision to be useful, at least with adequate sample sizes. We show that one of the novel pairwise methods provides estimates that are slightly better than the best global estimate, especially when the markers used have low exclusion probability. The new method can also be generalized to the case where there is no prior information on the density of reproductive adults. In that case, we can jointly estimate the density itself and the pollen dispersal distance, given sufficient sample sizes. The bias of this last estimator is larger and the precision is lower than for those estimates based on independent estimates of density, but the estimate is of some interest, because a meaningful independent estimate of the density of reproducing individuals is difficult to obtain in most cases.

Population extinction and recolonization in human demographic history

For over 30 years, a debate has raged among anthropologists about the origins of anatomically modern humans. At first the debate centered on fossil evidence, but in the past 10-15 years population geneticists have entered the fray. One model, the multiregional evolution model, posits a gradual transition from Homo erectus to anatomically modern humans throughout the Old World. In contrast, the recent African origin model hypothesizes that anatomically modern humans arose from a small, isolated population in Africa, spread out of Africa, and replaced indigenous H. erectus populations in Eurasia and Australasia. A primary objection, from a population genetics perspective, to the multiregional model is that the genetic data suggest a small effective population size for humans. This effective size is on the order of 10,000. Effective population size has a complex relationship with census size, but it has been argued that, assuming that effective size is roughly equal to the number of breeding individuals in human populations under standard demographic conditions, 10,000 breeding individuals could not have occupied much of the Old World throughout the Pleistocene and remained a cohesive species via gene flow. However, this argument is not valid if one considers population extinction and recolonization, which might have played an important role in human history during the Pleistocene. With population extinction and recolonization, the inbreeding effective population size can be small and the census size extremely large. In this paper, I will show that under conditions of population extinction and recolonization, an effective population size of 10,000 suggested by genetic data is compatible with a large census size consistent with the multiregional model.

Genetic fragmentation of populations of the fluvial sculpinCottus nozawae(Pisces: Cottidae) at the southern margin of its native range

The influence of fragmentation on local genetic population structure of the fluvial sculpin Cottus nozawae was studied by comparing populations in two regions: the Tohoku District (the southern margin of the species' native range) and Hokkaido Island (the center of its native range), Japan, based on restriction fragment length polymorphism analysis of mitochondrial DNA amplified by means of the polymerase chain reaction. Relatively high genetic divergences were found among populations from different regions. An analysis of molecular variance showed that ΦSTvalues among samples from the same river system in the Tohoku District were higher than among samples from Hokkaido Island. Heterogeneities in haplotype composition among tributary populations and among distinctive clusters were also found in the Tohoku District river system. Cottus nozawae in the Tohoku District are distributed patchily in headwater streams and may be thermally isolated. As a result, there appears to be little gene flow among tributary populations from different areas that exhibit divergence due to random genetic drift and have small effective population sizes.

Reconstructing recent divergence: evaluating nonequilibrium population structure in New Zealand chinook salmon

Newly established or perturbed populations are often the focus of conservation concerns but they pose special challenges for population genetics because drift-migration equilibrium is unlikely. To advance our understanding of the evolution of such populations, we investigated structure and gene flow among populations of chinook salmon that formed via natural straying following introduction to New Zealand in the early 1900s. We examined 11 microsatellite loci from samples collected in several sites and years to address two questions: (i) what population differentiation has arisen in the approximately 30 generations since salmon were introduced to New Zealand, relative to temporal variation within populations; and (ii) what are the approximate effective population sizes and amounts of gene flow in these populations? These questions are routinely addressed in studies of indigenous populations, but less often in the case of new populations and rarely with consideration of equilibrium assumptions. We show that despite the recent introduction, continued gene flow and high temporal variability among samples, detectable population structure has arisen among the New Zealand populations, consistent with their colonization pattern and isolation by geographical distance. Furthermore, we use simple individual-based simulations and estimates of effective population sizes to estimate the effective gene flow among drainages under likely nonequilibrium conditions. Similar methodology may be broadly applicable to other studies of population structure and phenotypic evolution under similar nonequilibrium, high gene flow conditions.

Genetic effective size of a wild primate population: influence of current and historical demography

A comprehensive assessment of the determinants of effective population size (N(e)) requires estimates of variance in lifetime reproductive success and past changes in census numbers. For natural populations, such information can be best obtained by combining longitudinal data on individual life histories and genetic marker-based inferences of demographic history. Independent estimates of the variance effective size (N(ev), obtained from life-history data) and the inbreeding effective size (N((eI), obtained from genetic data) provide a means of disentangling the effects of current and historical demography. The purpose of this study was to assess the demographic determinants of N(e) in one of the most intensively studied natural populations of a vertebrate species: the population of savannah baboons (Papio cynocephalus) in the Amboseli Basin, southern Kenya. We tested the hypotheses that N(eV) < N < N(eI) (where N = population census number) due to a recent demographic bottleneck. N(eV) was estimated using a stochastic demographic model based on detailed life-history data spanning a 28-year period. Using empirical estimates of age-specific rates of survival and fertility for both sexes, individual-based simulations were used to estimate the variance in lifetime reproductive success. The resultant values translated into an N(eV)/N estimate of 0.329 (SD = 0.116, 95% CI = 0.172-0.537). Historical N(eI), was estimated from 14-locus microsatellite genotypes using a coalescent-based simulation model. Estimates of N(eI) were 2.2 to 7.2 times higher than the contemporary census number of the Amboseli baboon population. In addition to the effects of immigration, the disparity between historical N(eI) and contemporary N is likely attributable to the time lag between the recent drop in census numbers and the rate of increase in the average probability of allelic identity-by-descent. Thus, observed levels of genetic diversity may primarily reflect the population's prebottleneck history rather than its current demography.

Evidence of low gene flow in a neotropical clustered tree species in two rainforest stands of French Guiana

The spatial genetic structure of the neotropical, clustered tree species Vouacapoua americana (Aublet) was studied in two natural forest stands (Paracou and Nouragues) in French Guiana. Using eight microsatellite loci, V. americana is characterized by a marked genetic structure at small spatial distances (under 30-60 m), in agreement with the limited seed dispersal by rodent species. Gene flow through pollen is also shown to be mainly restricted to less than 100 m. This result suggests that most pollination events (mediated through small insects) are probably limited to within-patches of individuals, which might explain the high genetic differentiation among patches (F(ST) = 0.11) separated by less than 2 km. We also assume that stronger genetic structure in Paracou is likely to be due to lower seed dispersal by rodents, large spatial distances separating patches, or a recent recolonization event.

Late Quaternary distributional stasis in the submediterranean mountain plant Anthyllis montana L. (Fabaceae) inferred from ITS sequences and amplified fragment length polymorphism markers

Anthyllis montana is a submediterranean, herbaceous plant of the southern and central European mountains. The internal transcribed spacer (ITS) regions of nuclear ribosomal DNA were sequenced from multiple accessions of the species and several closely related taxa. In addition, amplified fragment length polymorphism (AFLP) was analysed from 71 individuals of A. montana collected in 20 localities, mainly in the Pyrenees, Alps, Italian Peninsula and Balkans. Our ITS phylogeny showed a sequential branching pattern in A. montana, implying a western Mediterranean origin followed by an eastward migration. ITS clock calibrations suggest that speciation of A. montana took place at the Pliocene-Pleistocene boundary, while intraspecific divergence dates to Late Quaternary times (i.e. 0.7 million years ago). The AFLP analyses revealed a major genetic (west/east) subdivision within A. montana, probably caused by the massive glaciation of the Alps during this latter time period. The present-day absence of A. montana from vast parts of the Alps, which appear ecologically suitable for the species, together with the finding of evenly distributed AFLP variability within each of the two western and eastern lineages identified, is taken as evidence for a largely static Late Quaternary history without large-scale migration. High levels of AFLP variation observed among populations, together with weak or absent patterns of isolation by distance, seem to be in accord with long-term population insularization and distributional stasis. However, recent small-scale migration and a narrow hybrid zone between western and eastern lineages need to be postulated to explain the intermediate genetic composition of individuals from the Maritime Alps, a well-known suture-zone for other plant and animal species.

Genetic markers in the study of Anisakis typica (Diesing, 1860): larval identification and genetic relationships with other species of Anisakis Dujardin, 1845 (Nematoda: Anisakidae)

Genetic variation at 21 gene-enzyme systems was studied in a sample of an adult population of Anisakis typica (Diesing, 1860) recovered in the dolphin Sotalia fluviatilis from the Atlantic coast of Brazil. The characteristic alleles, detected in this population, made it possible to identify as A. typica, Anisakis larvae with a Type I morphology (sensu Berland, 1961) from various fishes: Thunnus thynnus and Auxis thazard from Brazil waters, Trachurus picturatus and Scomber japonicus from Madeiran waters, Scomberomorus commerson, Euthynnus affinis, Sarda orientalis and Coryphaena hippurus from the Somali coast of the Indian Ocean, and Merluccius merluccius from the Eastern Mediterranean. Characteristic allozymes are given for the identification, at any life-stage and in both sexes, of A. typica and the other Anisakis species so far studied genetically. The distribution of A. typica in warmer temperate and tropical waters is confirmed; the definitive hosts so far identified for this species belong to delphinids, phocoenids and pontoporids. The present findings represent the first established records of intermediate/paratenic hosts of A. typica and extend its range to Somali waters of the Indian Ocean and to the Eastern Mediterranean Sea. A remarkable genetic homogeneity was observed in larval and adult samples of A. typica despite their different geographical origin; interpopulation genetic distances were low, ranging from D(Nei)=0.004 (Eastern Mediterranean versus Somali) to D(Nei)=0.010 (Brazilian versus Somali). Accordingly, indirect estimates of gene flow gave a rather high average value of Nm = 6.00. Genetic divergence of A. typica was, on average, D(Nei)=1.12 from the members of the A. simplex complex (A. simplex s.s, A. pegreffii, A. simplex C) and D(Nei)=1.41 from A. ziphidarum, which all share Type I larvae; higher values were found from both A. physeteris (D(Nei)=2.77)

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.

The estimation of population differentiation with microsatellite markers

Microsatellite markers are routinely used to investigate the genetic structuring of natural populations. The knowledge of how genetic variation is partitioned among populations may have important implications not only in evolutionary biology and ecology, but also in conservation biology. Hence, reliable estimates of population differentiation are crucial to understand the connectivity among populations and represent important tools to develop conservation strategies. The estimation of differentiation is c from Wright's FST and/or Slatkin's RST, an FST -analogue assuming a stepwise mutation model. Both these statistics have their drawbacks. Furthermore, there is no clear consensus over their relative accuracy. In this review, we first discuss the consequences of different temporal and spatial sampling strategies on differentiation estimation. Then, we move to statistical problems directly associated with the estimation of population structuring itself, with particular emphasis on the effects of high mutation rates and mutation patterns of microsatellite loci. Finally, we discuss the biological interpretation of population structuring estimates.

Neutral evolution in spatially continuous populations

We introduce a general recursion for the probability of identity in state of two individuals sampled from a population subject to mutation, migration, and random drift in a two-dimensional continuum. The recursion allows for the interactions induced by density-dependent regulation of the population, which are inevitable in a continuous population. We give explicit series expansions for large neighbourhood size and for low mutation rates respectively and investigate the accuracy of the classical Malécot formula for these general models. When neighbourhood size is small, this formula does not give the identity even over large scales. However, for large neighbourhood size, it is an accurate approximation which summarises the local population structure in terms of three quantities: the effective dispersal rate, sigma(e); the effective population density, rho(e); and a local scale, kappa, at which local interactions become significant. The results are illustrated by simulations.

Great bustard population structure in central Spain: concordant results from genetic analysis and dispersal study

We found significant sex differences in the mtDNA genetic structure and dispersal patterns of great bustards in a population of 11 breeding groups, "leks", in central Spain. The analysis of genetic distances showed that the female population was divided into three groups of leks separated by ca. 50 km, whereas male haplotypes were randomly distributed among leks. Genetic distances among pairs of leks were positively correlated with geographical distances in females but not in males. While female haplotype distributions were homogeneous among leks at close distances, differences in male genetic structure were highly variable even between two close leks. These results from genetic analyses were concordant with those from a radiotracking study on natal dispersal. Natal dispersal distances were higher in males than in females. Also, the frequency of movement of a female between two leks was positively correlated with their genetic affinity and geographical proximity. In males, the frequency of movement was correlated with geographical proximity but not with genetic affinity. Males dispersed among genetically unrelated leks, contributing to keep nuclear genetic diversity in the population, whereas females tended to be philopatric. These results suggest that isolation-by-distance influences the distribution of maternal lineages at a regional level.