POPULATION GENETIC STRUCTURE AND GENE FLOW ACROSS ARID VERSUS MESIC ENVIRONMENTS: A COMPARATIVE STUDY OF TWO PARAPATRIC SENECIO SPECIES FROM THE NEAR EAST

To shed light on the potential effects of xeric/arid versus mesic environments on plant population genetic structure and patterns of gene flow, we have compared allozyme and cpDNA haplotype variation in populations of two closely related, highly outcrossed, and largely wind-dispersed winter annuals of Senecio (Asteraceae). The species form a distinctive zone of parapatric distribution in the Near East by differing in their ecogeographical regimes. Senecio vernalis mainly thrives in the mesic Mediterranean life zone of Israel, whereas S. glaucus inhabits either xeric maritime or arid (semi-) desert sites. Significant differences in allozymic population subdivision among S. vernalis (θ_n = 0.04; Nm_n = 5.85) and S. glaucus (θ_n = 0.12; Nm_n = 1.85) largely resulted from topogeographical substructuring present within the latter species. Because of the similarity of within-region estimates of population structure for S. glaucus with those measured among populations of S. vernalis, it appears unlikely that ecological “aridity” factors per se are important in influencing levels of population differentiation in these species. Based on hierarchical F-statistics and tests of isolation by distance, we further conclude that geographical topologies influence the level and mode of nuclear gene flow (via pollen and/or seed) among and within subsets of S. glaucus populations, although without providing a complete barrier to interregional dispersal (dNm_reg = 2.16) and without promoting allopatric differentiation via drift. The allozymic data further suggested that S. vernalis and S. glaucus form a zone of secondary contact in the Near East, accompanied by an almost complete interspecific barrier to nuclear gene flow (dnNm_sp = 0.253). However, to account for the considerable sharing of cpDNA haplotypes, both at the intra- and interspecific level, it is necessary to invoke either (1) selection acting against alien nuclear but not cytoplasmic DNA; or (2) the sporadic immigration of cpDNA via seed with large homogenizing effects on cytoplasmic population structure over time.     

POPULATION STRUCTURE AND LEVELS OF GENE FLOW IN THE MEDITERRANEAN LAND SNAIL ALBINARIA CORRUGATA (PULMONATA: CLAUSILIIDAE)

The amount of gene flow among local populations partly determines the relative importance of genetic drift and natural selection in the differentiation of such populations. Land snails, because of their limited powers for dispersal, may be particularly likely to show such differentiation. In this study, we directly estimate gene flow in Albinaria corrugata, a sedentary, rock-dwelling gastropod from Crete, by mark-recapture studies. In the same area, 23 samples were taken and studied electrophoretically for six polymorphic enzyme loci. The field studies indicate that the population structure corresponds closely to the stepping-stone model: demes are present on limestone boulders that are a few meters apart, and dispersal takes place mainly between adjacent demes. Average deme size (N) is estimated at 29 breeding individuals and the proportion of migrants per generation at 0.195 (Nm = 5.7). We find no reason to assume long-distance dispersal, apart from dispersal along occasional stretches of suitable habitat. Genetic subdivision of the population, as derived from F_ST values, corresponds to the direct estimate only at the lowest spatial level (distance between sample sites < 10 m), where values for Nm of 5.4 and 17.6 were obtained. In contrast, at the larger spatial scales, F_ST values give gene-flow estimates that are incompatible with the expected amount of gene flow at these scales. We explain these discrepancies by arguing that gene flow is in fact extremely limited, making correct estimates of Nm from F_ST impossible at the larger spatial scales. In view of these low levels of gene flow, it is concluded that both genetic drift and natural selection may play important roles in the genetic differentiation of this species, even at the lowest spatial scales.     

DYNAMICS OF POLYPLOID FORMATION IN TRAGOPOGON (ASTERACEAE): RECURRENT FORMATION,GENE FLOW,AND POPULATION STRUCTURE

Polyploidy is a major feature of angiosperm evolution and diversification. Most polyploid species have formed multiple times, yet we know little about the genetic consequences of recurrent formations. Among the clearest examples of recurrent polyploidy are Tragopogon mirus and T. miscellus (Asteraceae), each of which has formed repeatedly in the last ～80 years from known diploid progenitors in western North America. Here, we apply progenitor‐specific microsatellite markers to examine the genetic contributions to each tetraploid species and to assess gene flow among populations of independent formation. These data provide fine‐scale resolution of independent origins for both polyploid species. Importantly, multiple origins have resulted in considerable genetic variation within both polyploid species; however, the patterns of variation detected in the polyploids contrast with those observed in extant populations of the diploid progenitors. The genotypes detected in the two polyploid species appear to represent a snapshot of historical population structure in the diploid progenitors, rather than modern diploid genotypes. Our data also indicate a lack of gene flow among polyploid plants of independent origin, even when they co‐occur, suggesting potential reproductive barriers among separate lineages in both polyploid species.     

MOLECULAR PHYLOGEOGRAPHY, RETICULATION, AND LINEAGE SORTING IN MEDITERRANEAN SENECIO SECT. SENECIO (ASTERACEAE)

Abstract The Mediterranean species complex of Senecio serves to illustrate evolutionary processes that are likely to confound phylogenetic inference, including rapid diversification, gene tree‐species tree discordance, reticulation, interlocus concerted evolution, and lack of complete lineage sorting. Phylogeographic patterns of chloroplast DNA (cpDNA) haplotype variation were studied by sampling 156 populations (502 individuals) across 18 species of the complex, and a species phylogeny was reconstructed based on sequences from the internal transcribed spacer (ITS) regions of nuclear ribosomal DNA. For a subset of species, randomly amplified polymorphic DNAs (RAPDs) provided reference points for comparison with the cpDNA and ITS datasets. Two classes of cpDNA haplotypes were identified, with each predominating in certain parts of the Mediterranean region. However, with the exception of S. gallicus, intraspecific phylogeographic structure is limited, and only a few haplotypes detected were species‐specific. Nuclear sequence divergence is low, and several unresolved phylogenetic groupings are suggestive of near simultaneous diversification. Two well‐supported ITS clades contain the majority of species, amongst which there is a pronounced sharing of cpDNA haplotypes. Our data are not capable of diagnosing the relative impact of reticulation versus insufficient lineage sorting for the entire complex. However, there is firm evidence that S. flavus subsp. breviflorus and S. rupestris have acquired cpDNA haplotypes and ITS sequences from co‐occurring species by reticulation. In contrast, insufficient lineage sorting is a viable hypothesis for cpDNA haplotypes shared between S. gallicus and its close relatives. We estimated the minimum coalescent times for these haplotypes by utilizing the inferred species phylogeny and associated divergence times. Our data suggest that ancestral cpDNA polymorphisms may have survived for ca. 0.4–1.0 million years, depending on molecular clock calibrations.     

GENETIC POPULATION STRUCTURE AND LEVELS OF GENE FLOW IN THE STREAM DWELLING WATERSTRIDER,AQUARIUS (=GERRIS) REMIGIS (HEMIPTERA: GERRIDAE)

Gene flow, in combination with selection and drift, determines levels of differentiation among local populations. In this study we estimate gene flow in a stream dwelling, flightless waterstrider, Aquarius remigis. Twenty-eight Aquarius remigis populations from Quebec, Ontario, New Brunswick, Iowa, North Carolina, and California were genetically characterized at 15 loci using starch gel electrophoresis. Sampling over two years was designed for a hierarchical analysis of population structure incorporating variation among sites within streams, streams within watersheds, watersheds within regions, and regions within North America. Hierarchical F statistics indicated that only sites within streams maintained enough gene flow to prevent differentiation through drift (Nm = 27.5). Above the level of sites within streams gene flow is highly restricted (Nm ≤ 0.5) and no correlation is found between genetic and geographic distances. This agrees well with direct estimates of gene flow based on mark and recapture data, yielding an N_e of approximately 170 individuals. Previous assignment of subspecific status to Californian A. remigis is not supported by genetic distances between those populations and other populations in North America. Previous suggestion of specific status for south-eastern A. remigis is supported by genetic distances between North Carolina populations and other populations in North America, and a high proportion of region specific alleles in the North Carolina populations. However, because of the high degree of morphological and genetic variability throughout the range of this species, the assignment of specific or subspecific status to parts of the range may be premature.     

SPECIES-WIDE POPULATION STRUCTURE IN A SOUTHEASTERN U.S. FRESHWATER FISH,HETERANDRIA FORMOSA: GENE FLOW AND BIOGEOGRAPHY

The phylogeography of the freshwater fish fauna of the southeastern United States has almost achieved paradigm status in evolutionary biology (Avise 1992), and the major geographic features responsible for shaping species distributions are well-characterized. Nevertheless, variation among species in distributions of allele or haplotype frequencies suggests that species-specific processes (e.g., migration) may also play a role in establishing those distributions. There has also been relatively little investigation into how population structure may differ among subregions in the Southeast, for example, on the Florida peninsula versus the U.S. mainland to the northwest and/or northeast. The geology of the peninsula is such that both physical and biotic fluctuations may have been (and still be) particularly important in establishing the population structure of freshwater taxa. This possibility leads to two interesting questions in population genetics. (1) Does gene flow in freshwater species of the region better approximate a one- or two-dimensional pattern? (2) Are populations on the peninsula farther from migration-genetic drift equilibrium than their counterparts on the mainland? These questions were addressed by examining the population strucuture of a livebearing fish, Heterandria formosa; several features of the biology of the species make it particularly likely that recent gene flow has been important in its evolution. I surveyed electrophoretic variation in 34 populations distributed throughout the species range. The phylogeographic patterns observed are in general concordance with those found in other species, although with some differences. A two-dimensional hypothesis of gene flow on the Florida peninsula better explains the data than does a one-dimensional one. There is no evidence that populations on the peninsula are farther from migration-drift equilibrium than those to the northwest. Populations in the northeast have lower genetic diversity than those to the south and west and show no isolation by distance; those results are consistent with a recent range expansion into the northeast, although smaller historical effective population sizes could also explain the pattern.     

GENE FLOW,REFUGIA, AND EVOLUTION OF GEOGRAPHIC VARIATION IN THE SONG SPARROW (MELOSPIZA MELODIA)

We surveyed mtDNA restriction-site variation in song sparrows taken from across their continental range. Despite marked geographic variation in size and plumage color, mtDNA variation was not geographically structured. Subspecies were not identifiable by mtDNA analysis. We suggest that postglaciation dispersal scattered mtDNA haplotypes across the continent, explaining the lack of mtDNA geographic patterns. Evolution of size and plumage coloration has probably proceeded faster than mtDNA evolution, leading to the well-structured continental pattern of morphological variation. We suggest that the nonordered geographic distribution of haplotypes reflects the recency of population establishment following completion of range expansion. Dispersal distance was estimated from the mtDNA data at 6.1 km per generation, an order of magnitude greater than that (0.3 km) estimated from demographic data. Island samples were not especially different from continental ones. Rooting the haplotype cladogram with a putative primitive haplotype identified Newfoundland and the Queen Charlotte Islands as potential sites of recent refugia. We question whether study of geographic variation in song sparrows leads to insights concerning speciation.     

GENETIC STRUCTURE AND MALE-MEDIATED GENE FLOW IN THE GHOST BAT (MACRODERMA GIGAS)

The Australian ghost bat is a large, opportunistic carnivorous species that has undergone a marked range contraction toward more mesic, tropical sites over the past century. Comparison of mitochondrial DNA (mtDNA) control region sequences and six nuclear microsatellite loci in 217 ghost bats from nine populations across subtropical and tropical Australia revealed strong population subdivision (mtDNA φ_ST = 0.80; microsatellites UR_ST = 0.337). Low-latitude (tropical) populations had higher heterozygosity and less marked phylogeographic structure and lower subdivision among sites within regions (within Northern Territory [NT] and within North Queensland [NQ]) than did populations at higher latitudes (subtropical sites; central Queensland [CQ]), although sampling of geographically proximal breeding sites is unavoidably restricted for the latter. Gene flow among populations within each of the northern regions appears to be male biased in that the difference in population subdivision for mtDNA and microsatellites (NT φ_ST = 0.39, UR_ST = 0.02; NQ φ_ST = 0.60, UR_ST = ?0.03) is greater than expected from differences in the effective population size of haploid versus diploid loci. The high level of population subdivision across the range of the ghost bat contrasts with evidence for high gene flow in other chiropteran species and may be due to narrow physiological tolerances and consequent limited availability of roosts for ghost bats, particularly across the subtropical and relatively arid regions. This observation is consistent with the hypothesis that the contraction of the species' range is associated with late Holocene climate change. The extreme isolation among higher-latitude populations may predispose them to additional local extinctions if the processes responsible for the range contraction continue to operate.     

ECOLOGICAL AND HISTORICAL ASSOCIATIONS OF GENE FLOW IN DARTERS (TELEOSTEI: PERCIDAE)

Life history should relate to gene flow (Nm) through its influence on dispersal and effective population size. Comparative studies designed to elucidate this relationship must contend with historical events that can yield misleading estimates of gene flow and statistical problems associated with inclusion of life-history traits correlated with phylogeny. We studied the relationships of life-history characters and gene flow in 15 species of darters, a monophyletic group of stream fishes. Populations of coexisting species were sampled in three geographic regions with different Pleistocene glaciation histories. Gene flow was estimated indirectly from allozymes using two methods, 8 and private alleles. Isolation-by-distance was also tested using regression of pairwise estimates of gene flow (M?) on distance. Theta and private-alleles methods produced congruent estimates of Nm, except in a study region hypothesized to have been historically fragmented and then united following Pleistocene glaciation. A relatively weak association between life-history traits and Nm (based on θ) was observed when species from the historically fragmented region were included in stepwise regression analysis, because Nm was low despite life-history differences among taxa in this region. Excluding observations from this region produced stronger associations between clutch size and Nm (r² = 0.57), and between female size, egg size, and Nm (r² = 0.95). Additional analyses that corrected for female body size and phylogenetic nonindependence agreed that darters with high fecundity and small eggs exhibited high gene flow, whereas darters with small clutches and large eggs had low gene flow. The latter combination of life-history traits primarily is exhibited in species from headwater habitats where parental investment presumably confers survivorship on offspring. Reduced gene flow and genetic divergence among demes appear to be evolutionary consequences of this strategy.     

THE EFFECTS OF THE RELATIVE GEOGRAPHIC SCALES OF GENE FLOW AND SELECTION ON MORPH FREQUENCIES IN THE WALKING-STICK TIMEMA CRISTINAE

Gene frequencies in large populations are determined by a balance between selection and gene flow between neighborhoods of different selection regimes. This balance is affected by the area of the patches of a given selection regime relative to the gene-flow distance. If patches are small relative to gene-flow distance, similarity in the total area occupied by different patch types is a crucial condition for the stability of polymorphisms. However, if patches are larger than the gene-flow distance, then the relative area of different patch types is less important because of reduced gene flow resulting from isolation by distance. Two morphs (striped and unstriped) of the walking-stick Timema cristinae were each strongly associated with patches of distinct species of food plants on which they are most cryptic. The frequency of a morph was high on the plant on which it is most cryptic when either: (1) the area occupied by the food plant (patch) was very large; (2) the patch was completely isolated from other patches; or (3) the patch was larger than adjacent patches. Results (1) and (2) are consistent with isolation-by-distance models, and result (3) is consistent with Levene's multiple-niche polymorphism model.     

LONG-DISTANCE GENE FLOW IN THE SEDENTARY BUTTERFLY,EUPHILOTES ENOPTES (LEPIDOPTERA: LYCAENIDAE)

The relationship between gene flow and geographic proximity has been assessed for many insect species, but dispersal distances are poorly known for most of these. Thus, we are able to assess the concordance between vagility and gene flow for only a few species. In this study, I documented variation at six allozyme loci among Washington and Oregon populations of the sedentary, patchily distributed, lycaenid butterfly, Euphilotes enoptes (Boisduval) to assess whether the relationship between gene flow and geographic distance is consistent with the dispersal biology of this species. Both a phenogram based on genetic distances between populations and a regression analysis of gene flow estimates on geographic distances showed a pattern consistent with genetic isolation by distance. Many estimates of gene flow among pairs of populations separated by more than 100 km exceeded the equivalent of 10 individuals exchanged per generation, a value much greater than would be predicted from the limited dispersal ability of this species. However, based on the allozyme data, genetic neighborhood size was estimated to be approximately 39 individuals, a value that is consistent with poor vagility. The results of this study speak to the power of stepping-stone gene flow among populations and are compared to the results of other studies that have examined the relationship between dispersal and gene flow in sedentary insects.     

PHENOLOGICAL ISOLATION,GENE FLOW AND DEVELOPMENTAL DIFFERENCES AMONG LOW- AND HIGH-ELEVATION POPULATIONS OF EUPHILOTES ENOPTES (LEPIDOPTERA: LYCAENIDAE)

Populations of the specialist herbivore, Euphilotes enoptes (Lepidoptera: Lycaenidae), along three elevational transects in the mountains of central Washington state, differed markedly in the phenology of adult flight. In spite of this apparent limitation to gene flow, six allozyme loci revealed substantial gene exchange among populations along these gradients. The elevational difference, and thus the phenological difference, between populations has not influenced the extent of gene flow between them. Because the direct exchange of genes between low- and high-elevation populations is very unlikely, gene flow between them has probably occurred in a stepwise fashion via intermediate populations. It is hypothesized that such gene flow has been biased in an uphill direction due to the combined effects of source size and oviposition behavior. Adult emergence times of populations in the same region are positively correlated with elevation in a nonlinear fashion, consistent with the hypothesis that gene flow from low-elevation populations has been swamping selection at higher altitudes.     

CORRELATION OF PAIRWISE GENETIC AND GEOGRAPHIC DISTANCE MEASURES: INFERRING THE RELATIVE INFLUENCES OF GENE FLOW AND DRIFT ON THE DISTRIBUTION OF GENETIC VARIABILITY

Attempts to relate estimates of regional F_ST to gene flow and drift via Wright's (1931) equation F_ST ≈ 1/ (4Nm + 1) are often inappropriate because most natural sets of populations probably are not at equilibrium (McCauley 1993), as assumed by the island model upon which the equation is based, or ineffective because the influences of gene flow and drift are confounded in the product Nm. Evaluations of the association between genetic (F_ST) and geographic distances separating all pairwise populations combinations in a region allows one to test for regional equilibrium, to evaluate the relative influences of gene flow and drift on population structure both within and between regions, and to visualize the behavior of the association across all degrees of geographic separation. Tests of the model using microsatellite data from 51 populations of eastern collared lizards (Crotaphytus collaris collaris) collected from four distinct geographical regions gave results highly consistent with predicted patterns of association based on regional differences in various historical and ecological factors that affect the amount of drift and gene flow. The model provides a prerequisite for and an alternative to regional F_ST analyses, which often simply assume regional equilibrium, thus potentially leading to erroneous and misleading inferences regarding regional population structure.     

GENETIC STRUCTURE OF THE BLUE RIDGE DUSKY SALAMANDER (DESMOGNATHUS ORESTES): INFERENCES FROM ALLOZYMES, MITOCHONDRIAL DNA, AND BEHAVIOR

Abstract The plethodontid salamander Desmognathus orestes, a member of the D. ochrophaeus species complex, is distributed in southwestern Virginia, eastern Tennessee, and western North Carolina. Previous allozyme analyses indicate that D. orestes consists of two distinct groups of populations (D. orestes‘B’ and D. orestes‘C’) with extensive intergradation and probable gene flow between these two groups. Spatially varying allele frequencies can reflect historical associations, current gene flow, or a combination of population‐level processes. To differentiate among these processes, we use multiple markers to further characterize divergence among populations of D. orestes and assess the degree of intergradation between D. orestes‘B’ and D. orestes‘C’, specifically investigating variation in allozymes, mitochondrial DNA (mtDNA), and reproductive behavior among populations. On a broad scale, the mtDNA genealogies reconstruct haplotype clades that correspond to the species identified from previous allozyme analyses. However, at a finer geographic scale, the distributions of the allozyme and mtDNA markers for D. orestes‘B’ and D. orestes‘C’ are discordant. MtDNA haplotypes corresponding to D. orestes‘B’ are more broadly distributed across western North Carolina than predicted by allozyme data, and the region of intergradation with D. orestes‘C’ indicates asymmetric gene flow of these markers. Asymmetric mating may contribute to observed discordance in nuclear versus cytoplasmic markers. Results support describing D. orestes as a single species and emphasize the importance of using multiple markers to examine fine‐scale patterns and elucidate evolutionary processes affecting gene flow when making species‐level taxonomic decisions.     

HISTORICAL SEPARATION AND PRESENT GENE FLOW THROUGH A ZONE OF SECONDARY CONTACT IN PONDEROSA PINE

I examined the effects of historical division and secondary contact between eastern and western varieties of ponderosa pine (Pinus ponderosa Laws Pinaceae) on extant patterns of genetic variation. Fossil and biogeographic evidence both indicate that the current point of contact between these two varieties represents secondary contact following historical separation during the Wisconsin glaciation. Current gene flow was assessed by observing the degree of introgression of paternally inherited cpDNA and maternally inherited mtDNA polymorphisms. Both seeds and pollen are wind dispersed in ponderosa pine. Introgression was primarily from west to east, the direction of the prevailing wind, for both organelles, but introgression of cpDNA far exceeded that of mtDNA. Thus pollen is the main agent of contemporary gene flow between the two varieties. Neither seeds nor pollen showed enough introgression since secondary contact to have homogenized the two gene pools. However, allozyme differentiation was minimal. This calls into question assumptions of selective neutrality for at least some of the markers. Theory predicts that nuclear markers will show a high locus-to-locus variance of F_ST following historical separation. This prediction is confirmed by the allozyme data for ponderosa pine, and may provide a useful means of identifying historical separations from allele frequency data.     

SPATIAL SCALE OF GENETIC STRUCTURE AND AN INDIRECT ESTIMATE OF GENE FLOW IN EELGRASS,ZOSTERA MARINA

In this study, the first investigation of population structure in an aquatic angiosperm, I show that populations of a marine angiosperm (eelgrass, Zostera marina) are genetically differentiated at a number of spatial scales. I find also that there is no correspondence between geographic and genetic distances separating subpopulations, an increasingly common result in spatially stratified studies of genetic structure in marine invertebrates. F-statistics, calculated for two years from electrophoretic variation at five polymorphic allozyme loci, indicate significant genetic differentiation among sampling quadrats within each of two bays (θ = 0.064-0.208), between tide zones within a bay (θ = 0.025-0.157) and between bays (θ = 0.079). Spatial autocorrelation analysis was used to explore genetic differentiation at smaller spatial scales; estimated patch sizes (within which genetic individuals are randomly associated) indicated no appeciable genetic structure at scales less than 20 m × 20 m. Calculated values of F-statistics were a function of the spatial scale from which samples were drawn: increasing the size of the “subpopulation” included in calculation of fixation indices for the same “total” sample resulted in an increase in the magnitude of f (e.g., from 0.092 to 0.181) and a decrease in θ (e.g., from 0.186 to 0.025). On the basis of the best estimate of the spatial scale of subpopulations, the effective number of migrants per generation (N_em) ranges from 1.1 to 2.8. Genetic consequences of the disturbance regime in the eelgrass habitat sampled were extreme variation between years in the allele richness and proportion of heterozygotes in a sample and a positive relationship between the extinction probability of patches and the genetic variance among them. The changes in F-statistics as a function of sampling scale and the observation that θ among sampled quadrats was positively associated with the probability of extinction among quadrats indicated that indirect estimates of gene flow (N_em) calculated from θ should be cautiously interpreted in populations that may not yet be in drift-migration equilibrium.     

POPULATION GENETICS AND COLONY STRUCTURE OF THE ARGENTINE ANT (LINEPITHEMA HUMILE) IN ITS NATIVE AND INTRODUCED RANGES

Abstract.— Introduced species often possess low levels of genetic diversity relative to source populations as a consequence of the small population sizes associated with founder events. Additionally, native and introduced populations of the same species can possess divergent genetic structuring at both large and small geographic scales. Thus, genetic systems that have evolved in the context of high diversity may function quite differently in genetically homogeneous introduced populations. Here we conduct a genetic analysis of native and introduced populations of the Argentine ant (Linepithema humile) in which we show that the population‐level changes that have occurred during introduction have produced marked changes in the social structure of this species. Native populations of the Argentine ant are characterized by a pattern of genetic isolation by distance, whereas this pattern is absent in introduced populations. These differences appear to arise both from the effects of recent range expansion in the introduced range as well as from differences in gene flow within each range. Relatedness within nests and colonies is lower in the introduced range than in the native range as a consequence of the widespread genetic similarity that typifies introduced populations. In contrast, nestmates and colony‐mates in the native range are more closely related, and local genetic differentiation is evident. Our results shed light on the problem posed for kin selection theory by the low levels of relatedness that are characteristic of many unicolonial species and suggest that the loss of genetic variation may be a common mechanism for the transition to a unicolonial colony structure.     

UNDERSTANDING THE POPULATION GENETIC STRUCTURE OF GLEDITSIA TRIACANTHOS L.: THE SCALE AND PATTERN OF POLLEN GENE FLOW

Spatial and temporal patterns of gene flow determine the extent to which populations can differentiate from one another as a result of natural selection or genetic drift. In this study, we investigated pollen-mediated gene flow in two eastern Kansas populations of the subdioecious tree species, Gleditsia triacanthos L. (Leguminosae), or honeylocust. In 2 yr at each site, we used paternity-exclusion analysis to estimate the proportion of seeds sired by immigrant pollen. We also used a single-parent and parent-pair exclusion analysis on naturally established seedlings and saplings to estimate gene flow into one site over a 12-yr period and into the second site over a 22-yr period. Results of both analyses showed high minimum estimates of pollen gene flow into each site (17%–30%). In each population, we found significantly less gene flow in years of high fruit production than in years of low fruit production, but in one population, we observed little variation in gene-flow rates among age classes of seedlings and saplings. The level of pollen gene flow showed weak negative dependence on the relative isolation distances of the maternal trees sampled (140–240 m at one site vs. 85–120 m at the second site), and gene-flow estimates from naturally established juveniles were very similar at the two sites. Within populations, a multiple regression model showed that maximum-likelihood estimates of male fertility were negatively associated with distances between mates and positively associated with male size as measured by stem diameter. In neither population, however, did the regression explain more than 16% of the total variation in male fertilities.     

THE INFLUENCE OF POPULATION SIZE AND ISOLATION ON GENE FLOW BY POLLEN IN SILENE ALBA

In a series of experiments conducted over two seasons, we used arrays of experimental populations to examine the effects of flower number and distance between patches on gene flow by pollen. For this study we used the dioecious, short-lived perennial plant Silene alba (Caryophyllaceae). This species lives in disturbed roadside and agricultural habitats and displays a weedy population dynamic with high colonization and extinction rates. The motivation for the study was to understand what factors may be influencing genetic connectedness among newly colonized populations within a regional metapopulation. By using experimental populations composed of genotypes homozygous at a diagnostic locus, it was possible to identify explicitly pollen movement into a focal patch as a function of flower number and distance to the nearest neighboring patch. Overall, the mean immigration rate (measured as the fraction of seeds sired by males outside the focal patch) at 20 m was just over 47%, whereas at 80 m immigration rates were less than 6%. In addition, by knowing the context in which each of these gene-flow events occurred, it was possible to understand some of the factors that influenced the exchange of genes. Both the number of flowers in the focal population (target) and in the neighboring populations (source) had a significant effect on the frequency of gene flow. Our experimental data also demonstrate that factors that influence gene flow at one spatial scale may not act in the same way at another. Specifically, the influence of target size and the relative size of the target and source patches on rates of gene flow depended on whether the patches were separated by 20 m or 80 m. These data suggest that the patterns of gene flow within a metapopulation system can be complex and may vary within a growing season.     

LANDSCAPE GENOMICS OF POPULUS TRICHOCARPA: THE ROLE OF HYBRIDIZATION,LIMITED GENE FLOW,AND NATURAL SELECTION IN SHAPING PATTERNS OF POPULATION STRUCTURE

Populus trichocarpa is an ecologically important tree across western North America. We used a large population sample of 498 accessions over a wide geographical area genotyped with a 34K Populus SNP array to quantify geographical patterns of genetic variation in this species (landscape genomics). We present evidence that three processes contribute to the observed patterns: (1) introgression from the sister species P. balsamifera, (2) isolation by distance (IBD), and (3) natural selection. Introgression was detected only at the margins of the species’ distribution. IBD was significant across the sampled area as a whole, but no evidence of restricted gene flow was detected in a core of drainages from southern British Columbia (BC). We identified a large number of F_ST outliers. Gene Ontology analyses revealed that F_ST outliers are overrepresented in genes involved in circadian rhythm and response to red/far‐red light when the entire dataset is considered, whereas in southern BC heat response genes are overrepresented. We also identified strong correlations between geoclimate variables and allele frequencies at F_ST outlier loci that provide clues regarding the selective pressures acting at these loci.