How multilocus genotypic pattern helps to understand the history of selfing populations: a case study in Medicago truncatula

The occurrence of populations exhibiting high genetic diversity in predominantly selfing species remains a puzzling question, since under regular selfing genetic diversity is expected to be depleted at a faster rate than under outcrossing. Fine-scale population genetics approaches may help to answer this question. Here we study a natural population of the legume Medicago truncatula in which both the fine-scale spatial structure and the selfing rate are characterized using three different methods. Selfing rate estimates were very high ( approximately 99%) irrespective of the method used. A clear pattern of isolation by distance reflecting small seed dispersal distances was detected. Combining genotypic data over loci, we could define 34 multilocus genotypes. Among those, six highly inbred genotypes (lines) represented more than 75% of the individuals studied and harboured all the allelic variation present in the population. We also detected a large set of multilocus genotypes resembling recombinant inbred lines between the most frequent lines occurring in the population. This finding illustrates the importance of rare recombination in redistributing available allelic diversity into new genotypic combinations. This study shows how multilocus and fine-scale spatial analyses may help to understand the population history of self-fertilizing species, especially to make inferences about the relative role of foundation/migration and recombination events in such populations.     

MAINTENANCE OF POLYGENIC VARIATION IN SPATIALLY STRUCTURED POPULATIONS: ROLES FOR LOCAL MATING AND GENETIC REDUNDANCY

Although heritable genetic variation is critical to the evolutionary process, we know little about how it is maintained. Obviously, mutation-selection balance must play a role, but there is considerable doubt over whether it can account for heritabilities as high as 0.5, which are commonly found in natural populations. Most models of mutation-selection balance assume panmictic populations. In this paper we use Monte Carlo simulations to examine the effect of isolation by distance on the variation maintained by mutation in a polygenic trait subject to optimizing selection. We show that isolation by distance can substantially increase the total variation maintained in continuous populations over a wide range of dispersal patterns, but only if more than one genotype produces the optimal phenotype (genetic redundancy). Isolation by distance alone has only a slight effect on the variation maintained in the total population for neutral alleles. The combined effect of isolation by distance and genetic redundancy, however, allows the maintenance of substantial variation despite strong stabilizing selection. The mechanism is straightforward. Isolation by distance allows mutation and drift to operate independently in different parts of the population. Because of their independent evolutionary histories, different parts of the population independently draw from the available set of redundant genotypes. Because the genotypes are redundant, selection does not discriminate among them, and they will persist until eliminated by drift. The population as a whole maintains many distinct genotypes. We show that this process allows mutation to maintain high levels of variation, even under strong stabilizing selection, and that over a moderate range of dispersal patterns the amount of variation maintained in the entire population is independent of both the strength of selection and the variance of the dispersal distance. Furthermore, we show that individual heterozygosity is increased in locally mating populations when selection is strong. Finally, our simulations provide a rough picture of how selection and the dispersal pattern influence the spatial distribution of genetic and phenotypic variation.     

Gene flow in a facultative apomictic poacea, the savanna grass Hyparrhenia diplandra

The genetics of the poacea Hyparrhenia diplandra was studied in four natural populations from an ecological station in West Africa, where it makes up 80% of grasses from wet savanna and constitutes a dense continuum of randomly distributed individuals. DNA content and cytogenetical observations suggest it is an allotetraploid. Using two highly variable microsatellites (heterozygosity H = 0.615-0.616), we show that this species is an apomict with rare sexual reproduction events that account for approximately 0.5% of seeds pollinated in the wild. Hexaploid individuals were also produced, corroborating the observation of aberrant genotypes in the wild. The spatial extent of asexual clones in the field was low in comparison with the predominance of apomixis, thus indicating a low dispersal of seeds from their parent. Heterozygosity and departure from Hardy-Weinberg predictions were similar in the four populations, revealing a high apparent selfing rate s = 0.599 among sexually produced seeds. This is an overestimate since we could not distinguish true selfing from reciprocal outcrosses between neighboring individuals from the same apomictic clone. Gene flow by pollen could be substantial, possibly explaining the absence of isolation by distance in the studied area.     

Genetic Diversity and Population Structure in Polygonum cespitosum: Insights to an Ongoing Plant Invasion

Molecular markers can help elucidate how neutral evolutionary forces and introduction history contribute to genetic variation in invaders. We examined genetic diversity, population structure and colonization patterns in the invasive Polygonum cespitosum, a highly selfing, tetraploid Asian annual introduced to North America. We used nine diploidized polymorphic microsatellite markers to study 16 populations in the introduced range (northeastern North America), via the analyses of 516 individuals, and asked the following questions: 1) Do populations have differing levels of within-population genetic diversity? 2) Do populations form distinct genetic clusters? 3) Does population structure reflect either geographic distances or habitat similarities? We found low heterozygosity in all populations, consistent with the selfing mating system of P. cespitosum. Despite the high selfing levels, we found substantial genetic variation within and among P. cespitosum populations, based on the percentage of polymorphic loci, allelic richness, and expected heterozygosity. Inferences from individual assignment tests (Bayesian clustering) and pairwise F _ST values indicated high among-population differentiation, which indicates that the effects of gene flow are limited relative to those of genetic drift, probably due to the high selfing rates and the limited seed dispersal ability of P. cespitosum. Population structure did not reflect a pattern of isolation by distance nor was it related to habitat similarities. Rather, population structure appears to be the result of the random movement of propagules across the introduced range, possibly associated with human dispersal. Furthermore, the high population differentiation, genetic diversity, and fine-scale genetic structure (populations founded by individuals from different genetic sources) in the introduced range suggest that multiple introductions to this region may have occurred. High genetic diversity may further contribute to the invasive success of P. cespitosum in its introduced range.     

Spatial genetic structure within two populations of a self-incompatible perennial,Chionographis japonica var.japonica (Liliaceae)

Intrapopulational spatial genetic structure was examined in two populations ofChionographis japonica var.japonica, a self-incompatible perennial, by spatial autocorrelation analysis of enzyme polymorphism. Although most spatial autocorrelation indices (Moran'sI) in the shortes distance class were significantly positive, most in the other distance classes did not significantly deviate from the values expected from random distributions of genotypes in both populations. This contrasts with a spatial genetic pattern previously reported for a population of the predominantly selfing congener,C. japonica var.kurohimensis, indicating that pollen-mediated gene flow highly impedes genetic substructuring within populations of outcrossingC. japonica var.japonica. Genetic similarity in very proximate distance found in outcrossingC. japonica var.japonica is probably due to restricted dispersal of seeds.     

Effects of seed bank disturbance on the fine-scale genetic structure of populations of the rare shrub Grevillea macleayana

Dispersal in most plants is mediated by the movement of seeds and pollen, which move genes across the landscape differently. Grevillea macleayana is a rare, fire-dependent Australian shrub with large seeds lacking adaptations for dispersal; yet it produces inflorescences adapted to pollination by highly mobile vertebrates (eg birds). Interpreting fine-scale genetic structure in the light of these two processes is confounded by the recent imposition of anthropogenic disturbances with potentially contrasting genetic consequences: (1) the unusual foraging behaviour of exotic honeybees and 2. widespread disturbance of the soil-stored seedbank by road building and quarrying. To test for evidence of fine-scale genetic structure within G. macleayana populations and to test the prediction that such structure might be masked by disturbance of the seed bank, we sampled two sites in undisturbed habitat and compared their genetic structure with two sites that had been strongly affected by road building using a test for spatial autocorrelation of genotypes. High selfing levels inferred from genotypes at all four sites implies that pollen dispersal is limited. Consistent with this, we observed substantial spatial clustering of genes at 10 m or less in the two undisturbed populations and argue that this reflects the predicted effects of both high selfing levels and limited seed dispersal. In contrast, at the two sites disturbed by road building, spatial autocorrelation was weak. This suggests there has been mixing of the seed bank, counteracting the naturally low dispersal and elevated selfing due to honeybees. Pollination between near neighbours with reduced relatedness potentially has fitness consequences for G. macleayana in disturbed sites.     

Molecular evidence for limited dispersal of vegetative propagules in the epiphytic lichen Lobaria pulmonaria

Propagation, whether sexual or asexual, is a fundamental step in the life cycle of every organism. In lichenized fungi, a great variety of vegetative propagules have evolved in order for the symbiotic partners to disperse simultaneously. For lichens with the ability of sexual and asexual reproduction, the relative contribution of vegetative dispersal is unknown but could, nonetheless, be inferred by studying genotype distribution. The genetic structure of three Lobaria pulmonaria (Lobariaceae) populations from Switzerland was investigated based on the observed variation at six microsatellite loci. All three populations had a clustered distribution of identical genotypes at small spatial scales. The maximum distance between identical genotypes was 230 m. At a distance of 350 m from a source tree, seemingly suitable habitat patches were too far apart to be colonized. Some multilocus genotypes were frequent within local populations but no genotypes were shared among populations. The restricted occurrences of common genotypes as well as the clustered distributions are evidence for a limited dispersal of vegetative propagules in L. pulmonaria. Gene flow among isolated populations will ultimately depend on the capacity of long-distance dispersal and thus probably depend on sexual reproduction.     

Influence of landscape on the population genetic structure of the alpine butterfly parnassius smintheus (Papilionidae)

Four microsatellite DNA markers were developed which were used to examine the relationship between landscape and population genetic structure among a set of populations of the butterfly Parnassius smintheus located in the foothills of the Canadian Rockies. Detailed information on the dispersal of adult butterflies among this same set of populations was available. Simple and partial Mantel tests were used to examine the relationships between genetic distances, predicted rates of dispersal, and a number of landscape variables, all measured pairwise for 17 sample sites. Nei's standard genetic distance was negatively correlated with predicted dispersal. We observed a significant pattern of isolation by distance at a very small spatial scale. The distance between sites that was through forest was a stronger predictor of genetic distance than the distance through open meadow, indicating a significant effect of landscape on population genetic structure beyond that of simple isolation by distance. Our results suggest that rises in the tree-line in alpine areas, caused by global warming, will lead to reduced gene flow among populations of P. smintheus.     

Fine-scale spatial genetic structure with nonuniform distribution of individuals.

This paper presents the first theoretical study of spatial genetic structure within nonuniformly distributed continuous plant populations. A novel individual-based model of isolation by distance was constructed to simulate genetic evolution within such populations. We found larger values of spatial genetic autocorrelations in highly clumped populations than in uniformly distributed populations. Most of this difference was caused by differences in mean dispersal distances, but aggregation probably also produced a slight increase in spatial genetic structure. Using an appropriate level of approximation of the continuous distribution of individuals in space, we assessed the potential effects of density, seed and pollen dispersal, generation overlapping, and overdominance selection at an independent locus, on fine-scale genetic structure, by varying them separately in a few particular cases with extreme clumping. When selfing was allowed, all these input variables influenced both aggregation and spatial genetic structure. Most variations in spatial genetic structure were closely linked to variations in clumping and/or local density. When selfing was not allowed, spatial genetic structure was lower in most cases.     

Gene flow, dispersal, and nested clade analysis among populations of the stonefly Peltoperla tarteri in the southern Appalachians

We examined the genetic structure and phylogeography of populations of the stonefly Peltoperla tarteri in the Southern Appalachians to determine the extent and likely mechanism for dispersal of this stream insect. A 454-base-pair (bp) portion of the mitochondrial control region was sequenced from a minimum of 20 individuals from eight populations. Pairwise FST and exact tests showed high levels of differentiation among almost all populations except those on the same stream. amova analysis detected significant genetic differentiation between streams within drainages (phi(SD) = 0.14, P < 0.001), and there was a slight positive correlation between aquatic distance and genetic distance (r = 0.295, P = 0.03). According to nested clade analysis, the present day pattern of genetic variation in P. tarteri is the result of a historical range expansion coupled with restricted gene flow with isolation by distance. Together, these analyses suggest that adult dispersal is limited and that movement by larvae is the primary dispersal mechanism for P. tarteri.     

Consequences of multiple mating‐system shifts for population and range‐wide genetic structure in a coastal dune plant

Evolutionary transitions from outcrossing to selfing can strongly affect the genetic diversity and structure of species at multiple spatial scales. We investigated the genetic consequences of mating‐system shifts in the North American, Pacific coast dune endemic plant Camissoniopsis cheiranthifolia (Onagraceae) by assaying variation at 13 nuclear (n) and six chloroplast (cp) microsatellite (SSR) loci for 38 populations across the species range. As predicted from the expected reduction in effective population size (N_e) caused by selfing, small‐flowered, predominantly selfing (SF) populations had much lower nSSR diversity (but not cpSSR) than large‐flowered, predominantly outcrossing (LF) populations. The reduction in nSSR diversity was greater than expected from the effects of selfing on N_e alone, but could not be accounted for by indirect effects of selfing on population density. Although selfing should reduce gene flow, SF populations were not more genetically differentiated than LF populations. We detected five clusters of nSSR genotypes and three groups of cpSSR haplotypes across the species range consisting of parapatric groups of populations that usually (but not always) differed in mating system, suggesting that selfing may often initiate ecogeographic isolation. However, lineage‐wide genetic variation was not lower for selfing clusters, failing to support the hypothesis that selection for reproductive assurance spurred the evolution of selfing in this species. Within three populations where LF and SF plants coexist, we detected genetic differentiation among diverged floral phenotypes suggesting that reproductive isolation (probably postzygotic) may help maintain the striking mating‐system differentiation observed across the range of this species.     

Effects of kin-structured seed dispersal on the genetic structure of the clonal dioecious shrub ilex leucoclada

Nonrandom patterns of gene dispersal have been identified as possible causes of genetic structuring within populations. Attempts to model these patterns have generally focused solely on the effects of isolation by distance, but the processes involved are more complex than such modeling suggests. Here, we extend considerations of gene dispersal processes beyond simple isolation by distance effects by directly evaluating the effects of kin-structured gene dispersal mediated by the group dispersal of related seeds within fruits (i.e., kin-structured seed dispersal) by birds on genetic structure in Ilex leucoclada, a clonal dioecious shrub. To examine the genetic structure patterns, we established two 30x30 m plots (one with immature soils in old-growth forest and one in secondary forest, designated IM and SC, respectively) with different I. leucoclada stem densities. In these two plots 145 and 510 stems were found, representing 78 and 85 genets, respectively, identified by analyzing their genotypes at eight microsatellite loci. The clonal structure was stronger in the SC plot than in the IM plot. Correlograms of coancestry for genets in both plots exhibited significant, positive, high values in the shortest distance class, indicating the presence of strong genetic structure. However, Sp statistics revealed that the pattern of the genetic structure differed between the plots. In addition, to estimate the family structure within fruits, we sampled forty fruits, in total, from 15 randomly selected plants in the area around the IM and SC plots, and found that 80% of the fruits were multiseeded and 42-100% of the multiseeded fruits contained at least one pair of full sibs. Simulations based on these estimates demonstrated that the group dispersal of related seeds produced through correlated mating both within and across fruits, but not unstructured half-sib dispersal, could generate the observed magnitude and trends of genetic structure found in the IM plot. Furthermore, in addition to kin-structured seed dispersal, isolation by distance processes is also likely to promote genetic substructuring in the SC plot. After discussing possible ecological factors that may have contributed to the observed genetic structure, we contrast our results with those predicted by general isolation by distance models, and propose that kin-structured seed dispersal should promote some evolutionary phenomena, and thus should be incorporated, where appropriate, in models of gene dispersal in natural plant populations.     

Two invasive populations of the dry rot fungus Serpula lacrymans show divergent population genetic structures

The dry rot fungus Serpula lacrymans is a devastating basidiomycete occurring in wooden constructions in temperate regions worldwide. In this study, we compare the genetic structures of two invasive populations from Europe and Japan. Microsatellite data from 14 loci and DNA sequences from four loci demonstrated that the two populations were highly differentiated. Significant isolation by distance effect was observed in Europe and Japan. Higher genetic variation was observed within the Japanese population than within the European population, corresponding with the observed higher richness of vegetative compatibility types in Japan, indicating that there has been a higher level of gene flow from the Asian source populations to Japan than to Europe. The European population is genetically more homogenous with only six detected vegetative compatibility types. Various tests indicate that both the European and the Japanese populations have gone through population bottlenecks prior to population expansion. No identical multi‐locus genotypes were observed within Japan and very few within Europe, indicating limited clonal dispersal. Deviations from Hardy Weinberg expectations were observed both in Europe and Japan and heterozygote excess were observed at several loci, especially in Europe. Possible explanations for this pattern are discussed.     

An unusual case of seed dispersal in an invasive aquatic; yellow flag iris (<Emphasis Type="Italic">Iris pseudacorus</Emphasis>)

Invasive aquatic plants typically reproduce vegetatively, but there have been conflicting hypotheses of invasive aquatic yellow flag iris dispersing primarily by seed versus rhizome fragmentation. We performed genetic analysis of 20 aquatic yellow flag iris populations across the Pacific Northwest, USA, with leaf tissue taken from plants between 2 and 5 m apart. We found 167 unique genotypes in 171 plants, and we never found genetically identical plants from different populations. We found that 99.1 % of seed is viable. Our results support that this obligately outcrossing invasive disperses almost entirely by seed, not rhizome fragmentation. We found no significant relationship between genetic and geographic distance across the Pacific Northwest, suggesting recent long distance dispersal and/or multiple founding events. Bayesian analysis shows the presence of two major genotypic clusters within our collections which also suggests more than one distinct genetic source for the invasion. These processes have led to genetically distinct populations that can be geographically close. Our findings are unusual for an aquatic invasive, and inform yellow flag iris managers of two things: (1) to limit dispersal, development of mature seed in the field should be prevented; and (2) if classical biological control is proposed, an agent guild that limits seed production would be effective for managing most dispersal.     

Evaluating the influence of different aspects of habitat fragmentation on mating patterns and pollen dispersal in the bird-pollinated Banksia sphaerocarpa var. caesia

Habitat fragmentation can significantly affect mating and pollen dispersal patterns in plant populations, although the differential effects of the various aspects of fragmentation are poorly understood. In this study, we used eight microsatellite loci to investigate the effect of fragmentation on the mating system and pollen dispersal within one large and eight small population remnants of Banksia sphaerocarpa var. caesia, a bird-pollinated shrub in the southern agricultural region of Western Australia. The large population had a much larger neighbourhood size and lower selfing rate, maternal pollen pool differentiation and within-plot mean pollen dispersal distance than the small populations. Outcrossing was consistently high and ranged from 85.7% ± 2.6 to 98.5% ± 0.9, and mating patterns suggested nearest-neighbour pollination. Pollen immigration into small populations ranged from 2.8% ± 1.8 to 16.5% ± 3.2. Using the small populations, we tested for correlations between various fragmentation variables and mating system and pollen dispersal parameters. We found significant negative linear relationships between population isolation and outcrossing rate; population shape and neighbourhood size; and conspecific density and mean pollen dispersal distance. There were significant positive linear relationships between population shape and pollen pool differentiation and between population size and number of different fathers per seed crop. Our results suggest that birds may use a series of fragmented populations as a vegetation corridor while foraging across the landscape and that population connectivity is a critical determinant of pollinator visitation. Our results also suggest that the effect of a linear population shape on the mating system and pollen dispersal is routinely underestimated.     

Neither connectivity nor genetic diversity matter in the conservation of a rare fern and a moss on insular erratic boulders

Erratic boulders provide habitat for rock-dwelling species and contribute to the biodiversity of landscapes. In the calcareous Swiss lowlands, siliceous erratic boulders are exclusive habitat islands for the regionally critically endangered fern Asplenium septentrionale, about 20 bryophyte species and numerous lichens. Focusing on island biogeographical processes, we analysed the conservation genomics of A. septentrionale and the moss Hedwigia ciliata on insular erratic boulders in the Swiss lowlands and the adjacent “mainland” in siliceous mountains. We genotyped both species using double digest restriction associated DNA sequencing (ddRAD). For the tetraploid A. septentrionale, abundant identical multilocus genotypes within populations suggested prevalent intragametophytic selfing, and six out of eight boulder populations consisting of a single multilocus genotype each indicated single spore founder events. The genetic structure of A. septentrionale mainland populations coincided with Pleistocene glacial refugia. Four genetic lineages of H. ciliata were identified, and populations consisting of a single multilocus genotype were less common than in A. septentrionale. For both taxa, multilocus genotype diversity on boulders was lower than in mainland populations. The absence of common genetic groups among boulder populations, and the absence of isolation by distance patterns, suggested colonisation of boulders through independent long-distance dispersal events. Successful boulder colonisation of A. septentrionale seems to be rare, while colonisation by H. ciliata appears to be more frequent. We conclude that pivotal principles of conservation biology, such as connectivity and genetic diversity, are of less importance for the studied cryptogams on insular erratic boulders because of long-distance dispersal, intragametophytic selfing and polyploidy.

     

Apparent widespread gene flow in the predominantly flightless planthopper Tumidagena minuta

1. To determine whether dispersal biology can predict the pattern of population‐genetic variation among insect populations accurately, allozyme variation was assayed for populations of a saltmarsh planthopper, Tumidagena minuta, in which > 99% of the adults are flightless. 2. The pattern of genetic isolation by distance in T. minuta was compared with that in other insects, to determine whether it was similar to isolation by distance in other sedentary insects. 3. In contrast to predictions, the pattern of isolation by distance in T. minuta was most similar to that seen in the most mobile insects in a recent review of population‐genetic variation in insects. Furthermore, population‐genetic subdivision over a spatial scale of > 400 km was weak. 4. Possible causes of the apparent contradiction between dispersal biology and population‐genetic structure in this species are discussed. The results for T. minuta highlight the fact that although mobility is generally correlated with gene flow in insects, studies of population‐genetic variation must be combined with direct studies of dispersal to understand fully the degree to which populations exchange individuals.     

The Limits of Mean-Field Heterozygosity Estimates under Spatial Extension in Simulated Plant Populations

Computational models of evolutionary processes are increasingly required to incorporate multiple and diverse sources of data. A popular feature to include in population genetics models is spatial extension, which reflects more accurately natural populations than does a mean field approach. However, such models necessarily violate the mean field assumptions of classical population genetics, as do natural populations in the real world. Recently, it has been questioned whether classical approaches are truly applicable to the real world. Individual based models (IBM) are a powerful and versatile approach to achieve integration in models. In this study an IBM was used to examine how populations of plants deviate from classical expectations under spatial extension. Populations of plants that used three different mating strategies were placed in a range of arena sizes giving crowded to sparse occupation densities. Using a measure of population density, the pollen communication distance (P(cd)), the deviation exhibited by outbreeding populations differed from classical mean field expectations by less than 5% when P(cd) was less than 1, and over this threshold value the deviation significantly increased. Populations with an intermediate mating strategy did not have such a threshold and deviated directly with increasing isolation between individuals. Populations with a selfing strategy were influenced more by the mating strategy than by increased isolation. In all cases pollen dispersal was more influential than seed dispersal. The IBM model showed that mean field calculations can be reasonably applied to natural outbreeding plant populations that occur at a density in which individuals are less than the average pollen dispersal distance from their neighbors.     

Isolation by distance and gene flow in the Eurasian badger (Meles meles) at both a local and broad scale

Eurasian badgers, Meles meles, have been shown to possess limited genetic population structure within Europe; however, field studies have detected high levels of philopatry, which are expected to increase population structure. Population structure will be a consequence of both contemporary dispersal and historical processes, each of which is expected to be evident at a different scale. Therefore, to gain a greater understanding of gene flow in the badger, we examined microsatellite diversity both among and within badger populations, focusing on populations from the British Isles and western Europe. We found that while populations differed in their allelic diversity, the British Isles displayed a similar degree of diversity to the rest of western Europe. The lower genetic diversity occurring in Ireland, Norway and Scotland was more likely to have resulted from founder effects rather than contemporary population density. While there was significant population structure (F ST = 0.19), divergence among populations was generally well explained by geographic distance (P < 0.0001) across the entire range studied of more than 3000 km. Transient effects from the Pleistocene appear to have been replaced by a strong pattern of genetic isolation by distance across western Europe, suggestive of colonization from a single refugium. Analysis of individuals within British populations through Mantel tests and spatial autocorrelation demonstrated that there was significant local population structure across 3-30 km, confirming that dispersal is indeed restricted. The isolation by distance observed among badger populations across western Europe is likely to be a consequence of this restricted local dispersal.