Natural populations of Arabidopsis thaliana differ in seedling responses to high-temperature stress.

Little is known about adaptive within-species variation in thermotolerance in wild plants despite its likely role in both functional adaptation at range limits and in predicting response to climate change. Heat shock protein Hsp101, rapidly heat induced in Arabidopsis thaliana, plays a central role in thermotolerance in laboratory studies, yet little is known about variation in its expression in natural populations. We explored variation in thermotolerance and Hsp101 expression in seedlings from 16 natural populations of A. thaliana sampled along an elevation and climate gradient. We tested both naive controls (maintained at 22 °C until heat stress) and thermally pre-acclimated plants (exposed to a 38 °C 3-h acclimation treatment). After acclimation, seedlings were exposed to one of two heat stresses: 42 or 45 °C. Thermotolerance was measured as post-stress seedling survival and root growth. When stressed at 45 °C, both thermotolerance and Hsp101 expression were significantly increased by pre-acclimation. However, thermotolerance did not differ between pre-acclimation and control when followed by a 42 °C stress. Immediately after heat stress, pre-acclimated seedlings contained significantly more Hsp101 than control seedlings. At 45 °C, Hsp101 expression was positively associated with survival (r(2) = 0.37) and post-stress root growth (r(2) = 0.15). Importantly, seedling survival, post-stress root growth at 45 °C and Hsp101 expression at 42 °C were significantly correlated with the home sites' first principal component of climate variation. This climate gradient mainly reflects a temperature and precipitation gradient. Thus, the extent of Hsp101 expression modulation and thermotolerance appear to be interrelated and to evolve adaptively in natural populations of A. thaliana.

A geographic cline in leaf salicylic acid with increasing elevation in Arabidopsis thaliana.

Salicylic acid (SA) occupies a key role as a hormone central to both plant resistance to bacterial pathogens and tolerance of abiotic stresses. Plants at high elevation experience colder temperatures and elevated UV levels. While it has been predicted that SA concentrations will be higher in plants from high elevation populations, few studies have addressed this question. Here, we asked how concentrations of SA vary in natural populations of Arabidopsis thaliana collected across an elevational gradient on the Iberian Peninsula. In a series of common garden experiments, we found that constitutive SA concentrations were highest in genotypes from the low elevation populations. This result was in the opposite direction from our prediction and is an exception to the general finding that phenolic compounds increase with increasing elevation. These data suggest that high constitutive SA is not associated with resistance to cold temperatures in these plants. Furthermore, we also found that leaf constitutive camalexin concentrations, an important defense against some bacterial and fungal enemies, were highest in the low elevation populations, suggesting that pathogen pressures may be important. Further examination of this elevational cline will likely provide additional insights into the interplay between phenolic compounds and biotic and abiotic stress.

Constitutive camalexin production and environmental stress response variation in Arabidopsis populations from the Iberian Peninsula.

Optimal defense theory predicts that induction of defensive secondary metabolites in plants will be inversely correlated with constitutive expression of those compounds. Here, we asked whether camalexin, an important defense against fungal and bacterial pathogens, support this prediction in structured natural populations of Arabidopsis thaliana from the Iberian Peninsula. In common garden experiments, we found that genotypes from the VIE population constitutively hyper-accumulated camalexin. Camalexin concentrations were not induced significantly when plants were exposed to a temperature of 10°C for 48h. However, they were induced when plants were exposed to 48h of infection by the virulent bacterial pathogen, Pseudomonas syringae pv. tomato DC3000. Genotypes from the VIE population with the hyper-accumulation of camalexin were significantly more resistant to bacterial growth. Induction of camalexin was negatively correlated with constitutive camalexin concentrations following log transformation and two different corrections for autocorrelation, thus supporting the tradeoff predicted by optimal defense theory. Constitutive overexpression of camalexin was not explained by the only known natural genetic polymorphism at the Accelerated Cell Death 6, ACD6, locus. Collectively, the results support an important role of camalexin in defense against P. syringae as well as significant structured variation in defense levels within wild populations.

Adaptation to spring heat and drought in northeastern Spanish Arabidopsis thaliana.

The extent to which a species' environmental range reflects adaptive differentiation remains an open question. Environmental gradients can lead to adaptive divergence when differences in stressors among sites along the gradient place conflicting demands on the balance of stress responses. The extent to which this is accomplished through stress tolerance vs stress avoidance is also an open question. We present results from a controlled environment study of 48 lineages of Arabidopsis thaliana collected along a gradient in northeastern Spain across which temperatures increase and precipitation decreases with decreasing elevation. We tested the extent to which clinal adaptive divergence in heat and drought is explained through tolerance and avoidance traits by subjecting plants to a dynamic growth chamber cycle of increasing heat and drought stress analogous to low elevation spring in northeastern Spain. Lineages collected at low elevation were the most fit and fitness scaled with elevation of origin. Higher fitness was associated with earlier bolting, greater early allocation to increased numbers of inflorescences, reduction in rosette leaf photosynthesis and earlier fruit ripening. We propose that this is a syndrome of avoidance through early flowering accompanied by restructuring of the organism that adapts A. thaliana to low-elevation Mediterranean climates.

Developmental instability is genetically correlated with phenotypic plasticity, constraining heritability, and fitness.

Although adaptive plasticity would seem always to be favored by selection, it occurs less often than expected. This lack of ubiquity suggests that there must be trade-offs, costs, or limitations associated with plasticity. Yet, few costs have been found. We explore one type of limitation, a correlation between plasticity and developmental instability, and use quantitative genetic theory to show why one should expect a genetic correlation. We test that hypothesis using the Landsberg erecta × Cape Verde Islands recombinant inbred lines (RILs) of Arabidopsis thaliana. RILs were grown at four different nitrogen (N) supply levels that span the range of N availabilities previously documented in North American field populations. We found a significant multivariate relationship between the cross-environment trait plasticity and the within-environment, within-RIL developmental instability across 13 traits. This genetic covariation between plasticity and developmental instability has two costs. First, theory predicts diminished fitness for highly plastic lines under stabilizing selection, because their developmental instability and variance around the optimum phenotype will be greater compared to nonplastic genotypes. Second, empirically the most plastic traits exhibited heritabilities reduced by 57% on average compared to nonplastic traits. This demonstration of potential costs in inclusive fitness and heritability provoke a rethinking of the evolutionary role of plasticity.

Clinal variation in seed traits influencing life cycle timing in Arabidopsis thaliana.

Early life-history transitions are crucial determinants of lifetime survival and fecundity. Adaptive evolution in early life-history traits involves a complex interplay between the developing plant and its current and future environments. We examined the plant's earliest life-history traits, dissecting an integrated suite of pregermination processes: primary dormancy, thermal induction of secondary dormancy, and seasonal germination response. We examined genetic variation in the three processes, genetic correlations among the processes, and the scaling of germination phenology with the source populations' climates. A spring annual life history was associated with genetic propensities toward both strong primary dormancy and heat-induced secondary dormancy, alone or in combination. Lineages with similar proportions of winter and spring annual life history have both weak primary dormancy and weak thermal dormancy induction. A genetic bias to adopt a spring annual strategy, mediated by rapid loss of primary dormancy and high thermal dormancy induction, is associated with a climatic gradient characterized by increasing temperature in summer and rainfall in winter. This study highlights the importance of considering combinations of multiple genetically based traits along a climatic gradient as adaptive strategies differentiating annual plant life-history strategies. Despite the genetic-climatic cline, there is polymorphism for life-history strategies within populations, classically interpreted as bet hedging in an unpredictable world.

Population genomics and the causes of local differentiation.

Exactly 50 years ago, a revolution in empirical population genetics began with the introduction of methods for detecting allelic variation using protein electrophoresis (Throckmorton 1962; Hubby 1963; Lewontin & Hubby 1966). These pioneering scientists showed that populations are chock-full of genetic variation. This variation was a surprise that required a re-thinking of evolutionary genetic heuristics. Understanding the causes for the maintenance of this variation became and remains a major area of research. In the process of addressing the causes, this same group of scientists documented geographical genetic structure (Prakash et al. 1969), spawning the continued accumulation of what is now a huge case study catalogue of geographical differentiation (e.g. Loveless & Hamrick 1984; Linhart & Grant 1996). Geographical differentiation is clearly quite common. Yet, a truly general understanding of the patterns in and causes of spatial genetic structure across the genome remains elusive. To what extent is spatial structure driven by drift and phylogeography vs. geographical differences in environmental sources of selection? What proportion of the genome participates? A general understanding requires range-wide data on spatial patterning of variation across the entire genome. In this issue of Molecular Ecology, Lasky et al. (2012) make important strides towards addressing these issues, taking advantage of three contemporary revolutions in evolutionary biology. Two are technological: high-throughput sequencing and burgeoning computational power. One is cultural: open access to data from the community of scientists and especially data sets that result from large collaborative efforts. Together, these developments may at last put answers within reach.

Arabidopsis thaliana populations show clinal variation in a climatic gradient associated with altitude.

• Understanding the adaptive basis of life history variation is a central goal in evolutionary ecology. The use of model species enables the combination of molecular mechanistic knowledge with ecological and evolutionary questions, but the study of life history variation in natural environments is required to merge these disciplines. • Here, we tested for clinal variation in life history and associated traits along an environmental and altitudinal gradient in the model species Arabidopsis thaliana. Seventeen natural populations of A. thaliana were geo-referenced in north-eastern Spain on a gradient in which precipitation increases but maximum spring temperature and minimum winter temperature decrease with altitude. • One hundred and eighty-nine genotypes from the 17 populations were grown under uniform controlled conditions. Variations in traits related to biomass allocation, fecundity, phenology and vegetative growth were tested for relationships with the altitude and climatic variables associated with the home sites. Above-ground mass, number of rosette leaves at bolting, developmental time and seed weight increased with the home site's altitude. Root allocation, vegetative growth during winter and number of seeds decreased with altitude. • We suggest that the differences among home sites provide clues to the variation in adaptive strategies associated with the climatic gradient. We compared these results with adaptations and clinal relationships reported for other species and with molecular mechanisms described in Arabidopsis.

Something old and something new: wedding recombinant inbred lines with traditional line cross analysis increases power to describe gene interactions.

In this paper we present a novel approach to quantifying genetic architecture that combines recombinant inbred lines (RIL) with line cross analysis (LCA). LCA is a method of quantifying directional genetic effects (i.e. summed effects of all loci) that differentiate two parental lines. Directional genetic effects are thought to be critical components of genetic architecture for the long term response to selection and as a cause of inbreeding depression. LCA typically begins with two inbred parental lines that are crossed to produce several generations such as F1, F2, and backcrosses to each parent. When a RIL population (founded from the same P1 and P2 as was used to found the line cross population) is added to the LCA, the sampling variance of several nonadditive genetic effect estimates is greatly reduced. Specifically, estimates of directional dominance, additive x additive, and dominance x dominance epistatic effects are reduced by 92%, 94%, and 56% respectively. The RIL population can be simultaneously used for QTL identification, thus uncovering the effects of specific loci or genomic regions as elements of genetic architecture. LCA and QTL mapping with RIL provide two qualitatively different measures of genetic architecture with the potential to overcome weaknesses of each approach alone. This approach provides cross-validation of the estimates of additive and additive x additive effects, much smaller confidence intervals on dominance, additive x additive and dominance x dominance estimates, qualitatively different measures of genetic architecture, and the potential when used together to balance the weaknesses of LCA or RIL QTL analyses when used alone.

A simultaneous test of trophic interaction models: which vegetation characteristic explains herbivore control over plant community mass?

Predicting herbivore control over plants (i.e. changes in plant mass due to herbivore damage) is a central goal of ecology. Progress has been limited, however, because the vegetation characteristics thought to influence herbivore control are naturally correlated and typically experimentally confounded. To address this problem, we defined eight conventional models that predict herbivore control over plant community mass, each model based on a different vegetation characteristic (i.e. host concentration, tissue nitrogen, growth rate, size, tolerance of herbivory or net primary productivity). We then used structural equation modelling to test each model against two field experiments. Our results clearly rejected all models except for a tolerance of herbivory mechanism; stems with greater access to limiting resources better tolerated herbivory, regardless of where herbivore activity was greatest. Consequently, herbivore reductions of plant community mass were greatest at low resource availability. This adds to evidence that herbivore activity poorly predicts herbivore control.

Demographic and genetic patterns of variation among populations of Arabidopsis thaliana from contrasting native environments.

BACKGROUND: Understanding the relationship between environment and genetics requires the integration of knowledge on the demographic behavior of natural populations. However, the demographic performance and genetic composition of Arabidopsis thaliana populations in the species' native environments remain largely uncharacterized. This information, in combination with the advances on the study of gene function, will improve our understanding on the genetic mechanisms underlying adaptive evolution in A. thaliana. METHODOLOGY/PRINCIPAL FINDINGS: We report the extent of environmental, demographic, and genetic variation among 10 A. thaliana populations from Mediterranean (coastal) and Pyrenean (montane) native environments in northeast Spain. Geographic, climatic, landscape, and soil data were compared. Demographic traits, including the dynamics of the soil seed bank and the attributes of aboveground individuals followed over a complete season, were also analyzed. Genetic data based on genome-wide SNP markers were used to describe genetic diversity, differentiation, and structure. Coastal and montane populations significantly differed in terms of environmental, demographic, and genetic characteristics. Montane populations, at higher altitude and farther from the sea, are exposed to colder winters and prolonged spring moisture compared to coastal populations. Montane populations showed stronger secondary seed dormancy, higher seedling/juvenile mortality in winter, and initiated flowering later than coastal populations. Montane and coastal regions were genetically differentiated, montane populations bearing lower genetic diversity than coastal ones. No significant isolation-by-distance pattern and no shared multilocus genotypes among populations were detected. CONCLUSIONS/SIGNIFICANCE: Between-region variation in climatic patterns can account for differences in demographic traits, such as secondary seed dormancy, plant mortality, and recruitment, between coastal and montane A. thaliana populations. In addition, differences in plant mortality can partly account for differences in the genetic composition of coastal and montane populations. This study shows how the interplay between variation in environmental, demographic, and genetic parameters may operate in natural A. thaliana populations.

Evolutionary time for dispersal limits the extent but not the occupancy of species' potential ranges in the tropical plant genus Psychotria (Rubiaceae).

Explaining the diversity in geographic range sizes among species is a central goal of ecological and evolutionary studies. We tested species age as an explanation of range size variation within a group of understory shrubs in the Neotropics (Psychotria subgenus Psychotria, Rubiaceae). We distinguish between range occupancy (filling an occupied area) and range extent (maximum distances dispersed). We used Bayesian relaxed-clock dating of molecular sequence data to estimate the relative age of species, and we used species distribution modeling to predict species' potential ranges. If the range sizes of species are limited by time for dispersal, we hypothesize that older species should have (1) larger realized range occupancies and realized range extents than younger species, (2) filled a greater proportion of their potential range occupancies, and (3) colonized a greater proportion of their potential range extents. We found (1) a significant but weak positive relationship between species age versus both realized range occupancy and realized range extent, (2) no relationship between species age and filling of potential range occupancies, but (3) that older species had colonized a significantly greater proportion of their potential range extents than younger species. Our results indicate that a time-for-dispersal effect can limit the extent of ranges of species but not necessarily their occupancies.

Inflorescences contribute more than rosettes to lifetime carbon gain in Arabidopsis thaliana (Brassicaceae).

A metamorphosis from rosette to inflorescence in many annuals shifts photosynthetic tissue from a two-dimensional array in the soil boundary layer during cool months to a three-dimensional structure in the troposphere as spring progresses. We propose that this shift allows escape from both self-shading and an increasingly stressful boundary layer microclimate, permitting continued increases in growth. As a first step in exploring this hypothesis, we compared the lifetime C gain, water loss, and instantaneous water use efficiency (WUE) of five Arabidopsis thaliana genotypes by measuring gas exchange across the life cycle. On average, the inflorescence contributed 55% (± 5% SE) of lifetime C gain, but only 25% of lifetime water loss. Mean inflorescence WUE was nearly fourfold that of the rosette. The inflorescence continued to fix C after rosette senescence. The percentage inflorescence: total C gain varied among genotypes, from 36% to 93%. Genotypes differed in WUE for both structures. We suggest that local climates may have selected for divergence in these traits. For many annuals and winter annuals, understanding C and water budgets and their evolution must include measures of both rosette and inflorescence gas exchange.

Plastic trait integration across a CO2 gradient in Arabidopsis thaliana.

Shifts across environments in patterns of trait integration may govern or alter adaptive responses. Changes in resource supply rates may be an especially important cause of plasticity of trait integration because they can lead to shifts in colimitation and coregulation of traits. Traditional evolutionary genetic characterization of trait integration relies on covariance analyses. Structural equation modeling (SEM) can complement such analyses. The SEM provides insights into causal structure not possible with a covariance analysis, thereby providing mechanistic understanding of shifts in integration and suggesting likely foci of selection in changing environments. We tested for changes in trait integration by growing 35 genotypes of Arabidopsis thaliana (Brassicaceae; mouse-eared cress) from throughout the species' range in four atmospheric CO2 concentrations: 250 (past), 355 (approximately recent CO2), 530, and 710 (future) microM M(-1). SEM revealed significant shifts in the integration of N, C, and H2O use and their effects on reproductive dry mass across the CO2 gradient. The low CO2 stress of 250 microM M(-1) had the most divergent integration structures. Standardized total effects of C assimilation, water loss, and early N mass changed in sign across the C supply gradient, and the total effect of quantum yield decreased from significant to nonsignificant values across the gradient. Transpiration exhibited significant genetic variation and is thus a candidate target for selection and adaptation under novel growth CO2 concentrations. The strength of the correlation between C assimilation and transpiration declined by 19% from 250 to 710 microM M(-1), indicating a partial decoupling of their current mutual evolutionary constraint in the atmosphere of the future. Structural equation analysis of functional integration provides unique insights into the mechanisms through which changes in limiting resources can alter the nature of trait integration.

Flower color-flower scent associations in polymorphic Hesperis matronalis (Brassicaceae).

Floral scent emission rate and composition of purple and white flower color morphs of Hesperis matronalis (Brassicaceae) were determined for two populations and, for each, at two times of day using dynamic headspace collection and GC-MS. The floral volatile compounds identified for this species fell into two main categories, terpenoids and aromatics. Principal component analysis of 30 compounds demonstrated that both color morphs emitted more scent at dusk than at dawn. Color morphs varied in chemical composition of scent, but this differed between populations. The white morphs exhibited significant differences between populations, while the purple morphs did not. In the white morphs, one population contains color-scent associations that match expectations from classical pollination syndrome theory, where the flowers have aromatic scents, which are expected to maximize night-flying moth pollinator attraction; in the second population, white morphs were strongly associated with terpenoid compounds. The potential impact that pollinators, conserved biosynthetic pathways, and the genetics of small colonizing populations may have in determining population-specific associations between floral color and floral scent are discussed.

EVOLUTIONARY PREDICTABILITY IN NATURAL POPULATIONS: DO MATING SYSTEM AND NONADDITIVE GENETIC VARIANCE INTERACT TO AFFECT HERITABILITIES IN PLANTAGO LANCEOLATA?

Quantitative genetics has been an immensely powerful tool in manipulating the phenotypes of domesticated plants and animals. Much of the predictive power of quantitative genetics depends on the breeder's control over the context in which phenotype and mating are being expressed. In the natural world, these contexts are often difficult to describe, let alone control. We are left, therefore, with a poor understanding of the limits of quantitative genetics in natural populations. One of the crucial contextual elements for assessing breeding value is the genetic background in which an individual's genes are being assessed. When interacting genes are polymorphic within a population, the degree of mating among relatives can influence the correlations among mates and the predictions of a response to selection. Population structure can strongly influence the degree to which dominance and epistasis influences additive genetic variance and heritability. The extent of inbreeding can also influence heritabilities through its effect on the environmental component of phenotypic variance. The applicability of standard quantitative genetic breeding designs to the measurement of heritabilities in natural populations therefore depends in part on: (1) the mating system of the population; and (2) the importance of gene interactions in determining phenotypic variation. We tested for an effect of mating structure on the partitioning of phenotypic variance and heritability by comparing two breeding designs in a common environment. Both breeding designs used 139 pollen parents taken from mapped locations in a population of Plantago lanceolata L., and crossed to 280 seed parents from the same population. One design was random-mating, the second was biased toward near-neighbor matings to an extent determined by field measure of pollen-mediated gene flow distances. The offspring were grown randomly mixed in a common garden. Nine traits were measured: central corm diameter, number of leaves, area of the most recently fully expanded leaf, density of hairs (cm-2 ) on the leaves, dry weight per unit leaf area, photosynthetic capacity, transpiration rates, water use efficiency, and reproductive dry weight. Heritabilities and variance components from the two designs were compared using randomization tests. None of the variance components or the heritabilities differed significantly between breeding designs at the 0.05 level. The test could distinguish differences between the heritabilities measured in the two breeding designs as small as 0.11, on average. Thus, for the degree of inbreeding normally exhibited in P. lanceolata there is insufficient gene interaction present within populations to influence the partitioning of variance between additive and nonadditive components or to influence heritability estimates to a meaningful extent. We suggest that for Plantago other sources of variation in heritability estimates, such as maternal effects and genotype × environment interactions, are more important influences than the interaction between inbreeding and gene interactions, and standard heritability estimate based on random breeding is as accurate as one taking the natural mating structure into account.

A SIMULATION OF WRIGHT'S SHIFTING-BALANCE PROCESS: MIGRATION AND THE THREE PHASES.

Wright partitioned the shifting-balance process into three phases. Phase one is the shift of a deme within a population to the domain of a higher adaptive peak from that of the historical peak. Phase two is mass selection within a deme towards that higher peak. Phase three is the conversion of additional demes to the higher peak. The migration rate between demes is critical for the existence of phases one and three. Phase one requires small effective population sizes, hence low migration rates. Phase three is optimal under high migration rates that spread the most-fit genotype from deme to deme. Thus, a population-wide peak shift requires intermediate levels of migration. By altering the rates of phases one and three, migration affects the predominant direction of mass selection within a population. This study examines the degree to which migration, through its effects on phases one and three, determines the probability of a simulated population arriving at its genotypic optimum after 12,000 generations. These simulations reveal that there is a range of migration rates for which an entire population might be expected to shift to a higher peak. Below m = 0.001 peak shifts occur frequently (phases I and II) but are not successfully exported out of subpopulations (phase III), and above 0.01 peak shifts within demes (phase I and II), required to initiate phase III, become increasingly uncommon. Because it is unlikely that real populations will have uniform migration rates from generation to generation, the probable effects of varying migration rates on broadening the range of conditions producing peak shifts are discussed.

A LIFE-HISTORY BASED STUDY OF POPULATION GENETIC STRUCTURE: SEED BANK TO ADULTS IN PLANTAGO LANCEOLATA.

We explored the extent to which the soil seed bank differed genetically and spatially in comparison to two actively growing stages in a natural population of Plantago lanceolata. All seed-bank seeds, seedlings, and adults of P. lanceolata within eight subunits in a larger population were mapped, subjected to starch gel electrophoresis, and allozyme analysis in 1988. Gel electrophoresis was also used to estimate the mating system in two years, 1986 and 1988. The spatial distributions of seeds, seedlings, and adults were highly coincident. Allele frequencies of the dormant seeds differed significantly from those of the adults for four of the five polymorphic loci. In addition, a comparison of the genotype frequencies of the three life-history stages indicated that the seed bank had an excess of homozygotes. Homozygosity, relative to Hardy-Weinberg expectations, decreased during the life cycle (for seed bank, seedlings, and adults respectively: Fit = 0.19, 0.09, 0.01; Fis = 0.14, 0.04, -0.12). Spatial genetic differentiation increased sixfold during the life cycle: (for seed bank, seedling and adults: Fs1 ∗∗∗ = 0.02, 0.05, 0.12). The apparent selfing rate was 0.01 in 1986 and 0.09 in 1988. These selfing rates are not large enough to account for the elevated homozygosity of the seed bank. Inbreeding depression, overdominance for fitness, and a "temporal Wahlund's effect" are discussed as possible mechanisms that could generate high homozygosity in the seed bank, relative to later life-history stages. In Plantago lanceolata, the influence of the mating system and the "genetic memory" of the seed bank are obscured by the time plants reach the reproductive stage.

Leptokurtic pollen-flow, non-leptokurtic gene-flow in a wind-pollinated herb, Plantago lanceolata L.

The purpose of this study was to simultaneously measure pollen dispersal distance and actual pollen-mediated gene-flow distance in a wind-pollinated herb, Plantago lanceolata. The pollen dispersal distribution, measured as pollen deposition in a wind tunnel, is leptokurtic, as expected from previous studies of wind-pollinated plants. Gene-flow, measured as seeds produced on rows of male-sterile inflorescences in the wind tunnel, is non-leptokurtic, peaking at an intermediate distance. The difference between the two distributions results from the tendency of the pollen grains to cluster. These pollen clusters are the units of gene dispersal, with clusters of intermediate and large size contributing disproportionately to gene-flow. Since many wind-pollinated species show pollen clustering (see text), the common assumption for wind-pollinated plants that gene-flow is leptokurtic requires re-examination. Gene-flow was also measured in an artifical outdoor population of male-steriles, containing a single pollen source plant in the center of the array. The gene flow distribution is significantly platykurtic, and has the same general properties outdoors, where wind speed and turbulence are uncontrolled, as it does in the wind tunnel. I estimated genetic neighborhood size based on my measure of gene-flow in the outdoor population. The estimate shows that populations of Plantago lanceolata will vary in effective number from a few tens of plants to more than five hundred plants, depending on the density of the population in question. Thus, the measured pollen-mediated gene-flow distribution and population density will interact to produce effective population sizes ranging from those in which there is no random genetic drift to those in which random genetic drift plays an important role in determining gene frequencies within and among populations. Despite the platykurtosis in the distribution, pollen-mediated gene dispersal distances are still quite limited, and considerable within and among-population genetic differentiation is to be expected in this species.