Inferring colonization history from analyses of spatial genetic structure within populations of Pinus strobus and Quercus rubra

Many factors interact to determine genetic structure within populations including adult density, the mating system, colonization history, natural selection, and the mechanism and spatial patterns of gene dispersal. We examined spatial genetic structure within colonizing populations of Quercus rubra seedlings and Pinus strobus juveniles and adults in an aspen-white pine forest in northern Michigan, USA. A 20-year spatially explicit demographic study of the forest enables us to interpret the results in light of recent colonization of the site for both species. We assayed 217 Q. rubra seedlings and 171 P. strobus individuals at 11 polymorphic loci using nine allozyme systems. Plant genotypes and locations were used in an analysis of spatial genetic structure. Q. rubra and P. strobus showed similar observed levels of heterozygosity, but Q. rubra seedlings have less heterozygosity than expected. Q. rubra seedlings show spatial genetic clumping of individuals on a scale to 25 m and levels of genetic relatedness expected from the clumped dispersion of half-siblings. In contrast, P. strobus has low levels of genetic relatedness at the smallest distance class and positive spatial genetic structure at scales < 10 m within the plot. The low density of adult Q. rubra outside the study plot and limited, spatially clumped rodent dispersal of acorns is likely responsible for the observed pattern of spatial genetic structure and the observed heterozygote deficit (i.e. a Wahlund effect). We attribute weaker patterns observed in P. strobus to the longer dispersal distance of seeds and the historical overlap of seed shadows from adults outside of the plot coupled with the overlap of seed shadows from younger, more recently established reproductive adults. The study demonstrates the utility of long-term demographic data in interpreting mechanisms responsible for generating contemporary patterns of genetic structure within populations.     

Comparative analysis of the within-population genetic structure in wild cherry (Prunus avium L.) at the self-incompatibility locus and nuclear microsatellites

Gametophytic self-incompatibility (SI) systems in plants exhibit high polymorphism at the SI controlling S-locus because individuals with rare alleles have a higher probability to successfully pollinate other plants than individuals with more frequent alleles. This process, referred to as frequency-dependent selection, is expected to shape number, frequency distribution, and spatial distribution of self-incompatibility alleles in natural populations. We investigated the genetic diversity and the spatial genetic structure within a Prunus avium population at two contrasting gene loci: nuclear microsatellites and the S-locus. The S-locus revealed a higher diversity (15 alleles) than the eight microsatellites (4-12 alleles). Although the frequency distribution of S-alleles differed significantly from the expected equal distribution, the S-locus showed a higher evenness than the microsatellites (Shannon's evenness index for the S-locus: E = 0.91; for the microsatellites: E = 0.48-0.83). Also, highly significant deviations from neutrality were found for the S-locus whereas only minor deviations were found for two of eight microsatellites. A comparison of the frequency distribution of S-alleles in three age-cohorts revealed no significant differences, suggesting that different levels of selection acting on the S-locus or on S-linked sites might also affect the distribution and dynamics of S-alleles. Autocorrelation analysis revealed a weak but significant spatial genetic structure for the multilocus average of the microsatellites and for the S-locus, but could not ascertain differences in the extent of spatial genetic structure between these locus types. An indirect estimate of gene dispersal, which was obtained to explain this spatial genetic pattern, indicated high levels of gene dispersal within our population (sigma(g) = 106 m). This high gene dispersal, which may be partly due to the self-incompatibility system itself, aids the effective gene flow of the microsatellites, thereby decreasing the contrast between the neutral microsatellites and the S-locus.     

Genetic structure in chloroplast and nuclear microsatellites in Rosa rugosa around sea straits in northern Japan

Sea‐strait geological history and seawater dispersal patterns can affect the genetic structure of coastal plants. We investigated microsatellite (simple sequence repeat, SSR) variations in chloroplast (cp) and nuclear (n) DNA of the coastal shrub R osa rugosa, native to northeastern Asia, and three species closely related to it. In northern Japan, we sampled R. rugosa along coasts around the Soya and Tsugaru straits with different histories of the Quaternary land formation as a result of different water depths. Both cpSSR and nSSR variations suggest rare hybridization between R . rugosa and the closely related species. Variations in one of two cpSSRs showed genetic differentiation between eastern and western Hokkaido. This genetic structure may result from introgression from R. davurica in eastern Hokkaido and/or isolation between the eastern and western coastlines of Hokkaido as a result of the Quaternary land formation in the Soya Strait. On the other hand, variations in 10 nSSRs were geographically homogeneous with weak isolation by distance along coastlines. These results suggest that extensive gene flow has been homogenizing the genetic structure of R . rugosa in northern Japan.     

Temporal and spatial genetic structure in Vitellaria paradoxa (shea tree) in an agroforestry system in southern Mali

Ten microsatellite loci were used to investigate the impact of human activity on the spatial and temporal genetic structure of Vitellaria paradoxa (Sapotaceae), a parkland tree species in agroforestry systems in southern Mali. Two stands (forest and fallow) and three cohorts (adults, juveniles and natural regeneration) in each stand were studied to: (i) compare their levels of genetic diversity (gene diversity, HE; allelic richness, Rs; and inbreeding, FIS); (ii) assess their genetic differentiation (FST); and (iii) compare their levels of spatial genetic structuring. Gene diversity parameters did not vary substantially among stands or cohorts, and tests for bottleneck events were nonsignificant. The inbreeding coefficients were not significantly different from zero in most cases (FIS = -0.025 in forest and 0.045 in fallow), suggesting that the species is probably outbreeding. There was a weak decrease in F(IS) with age, suggesting inbreeding depression. Differentiation of stands within each cohort was weak (FST = 0.026, 0.0005, 0.010 for adults, juveniles and regeneration, respectively), suggesting extensive gene flow. Cohorts within each stand were little differentiated (FST = -0.001 and 0.001 in forest and fallow, respectively). The spatial genetic structure was more pronounced in fallow than in forest where adults showed no spatial structuring. In conclusion, despite the huge influence of human activity on the life cycle of Vitellaria paradoxa growing in parkland systems, the impact on the pattern of genetic variation at microsatellite loci appears rather limited, possibly due to the buffering effect of extensive gene flow between unmanaged and managed populations.     

Fine‐scale spatial genetic structure of common and declining bumble bees across an agricultural landscape

Land‐use changes have threatened populations of many insect pollinators, including bumble bees. Patterns of dispersal and gene flow are key determinants of species' ability to respond to land‐use change, but have been little investigated at a fine scale (<10 km) in bumble bees. Using microsatellite markers, we determined the fine‐scale spatial genetic structure of populations of four common Bombus species (B. terrestris, B. lapidarius, B. pascuorum and B. hortorum) and one declining species (B. ruderatus) in an agricultural landscape in Southern England, UK. The study landscape contained sown flower patches representing agri‐environment options for pollinators. We found that, as expected, the B. ruderatus population was characterized by relatively low heterozygosity, number of alleles and colony density. Across all species, inbreeding was absent or present but weak (F_IS = 0.01–0.02). Using queen genotypes reconstructed from worker sibships and colony locations estimated from the positions of workers within these sibships, we found that significant isolation by distance was absent in B. lapidarius, B. hortorum and B. ruderatus. In B. terrestris and B. pascuorum, it was present but weak; for example, in these two species, expected relatedness of queens founding colonies 1 m apart was 0.02. These results show that bumble bee populations exhibit low levels of spatial genetic structure at fine spatial scales, most likely because of ongoing gene flow via widespread queen dispersal. In addition, the results demonstrate the potential for agri‐environment scheme conservation measures to facilitate fine‐scale gene flow by creating a more even distribution of suitable habitats across landscapes.     

Population genetic structure of Bombus terrestris in Europe: Isolation and genetic differentiation of Irish and British populations

The genetic structure of the earth bumblebee (Bombus terrestris L.) was examined across 22 wild populations and two commercially reared populations using eight microsatellite loci and two mitochondrial genes. Our study included wild bumblebee samples from six populations in Ireland, one from the Isle of Man, four from Britain and 11 from mainland Europe. A further sample was acquired from New Zealand. Observed levels of genetic variability and heterozygosity were low in Ireland and the Isle of Man, but relatively high in continental Europe and among commercial populations. Estimates of F_st revealed significant genetic differentiation among populations. Bayesian cluster analysis indicated that Irish populations were highly differentiated from British and continental populations, the latter two showing higher levels of admixture. The data suggest that the Irish Sea and prevailing south westerly winds act as a considerable geographical barrier to gene flow between populations in Ireland and Britain; however, some immigration from the Isle of Man to Ireland was detected. The results are discussed in the context of the recent commercialization of bumblebees for the European horticultural industry.     

Fine-scale geographical structure of genetic diversity in inland wild beet populations

Introgression arising from crop-to-wild gene flow provides novel sources of genetic variation in plant species complexes. Hybridization within the Beta vulgaris species complex is of immediate concern; crop lineages ( B .  vulgaris ssp. vulgaris ) hybridize easily with their wild relatives ( B .  vulgaris ssp. maritima ) thereby threatening wild beet gene diversity with genetic swamping. Hybridization 'hotspots' occur in European seed production areas because inland ruderal wild beets occur and reproduce in sympatry with cultivated beets. We studied gene flow occurring between seed-producing cultivars and ruderal wild B .  vulgaris in southwestern France to determine whether feral beets, arising from unharvested cultivated seed, represent an opportunity for crop-to-wild gene flow. We surveyed 42 inland ruderal beet populations located near seed production fields for nucleo-cytoplasmic variation and used a cytoplasmic marker diagnostic of cultivated lines. Occurrence of cultivated-type cytoplasm within ruderal populations clearly reflected events of crop seed escape. However, we found no genetic signatures of nuclear cultivated gene introgression, which suggests past introgression of cultivated cytoplasm into a wild nuclear background through seed escape rather than recent direct pollen flow. Overall, patterns of genetic structure suggested that inland ruderal wild beet populations act as a metapopulation, with founding events involving a few sib groups, followed by low rates of seed or pollen gene flow after populations are established. Altogether, our results indicate that a long-lived seed bank plays a key role in maintaining cultivated-type cytoplasm in the wild and highlight the need for careful management of seed production areas where wild and cultivated relatives co-occur.     

Inferring recruitment history from spatial genetic structure within populations of the colonizing tree Albizia julibrissin (Fabaceae)

Comparative analyses of spatial genetic structure (SGS) among species, populations, or cohorts give insight into the genetic consequences of seed dispersal in plants. We analysed SGS of a weedy tree in populations with known and unknown recruitment histories to first establish patterns in populations with single vs. multiple founders, and then to infer possible recruitment scenarios in populations with unknown histories. We analysed SGS in six populations of the colonizing tree Albizia julibrissin Durazz. (Fabaceae) in Athens, Georgia. Study sites included two large populations with multiple, known founders, two small populations with a single, known founder, and two large populations with unknown recruitment histories. Eleven allozyme loci were used to genotype 1385 individuals. Insights about the effects of colonization history from the SGS analyses were obtained from correlograms and Sp statistics. Distinct differences in patterns of SGS were identified between populations with multiple founders vs. a single founder. We observed significant, positive SGS, which decayed with increasing distance in the populations with multiple colonists, but little to no SGS in populations founded by one colonist. Because relatedness among individuals is estimated relative to a local reference population, which usually consists of those individuals sampled in the study population, SGS in populations with high background relatedness, such as those with a single founder, may be obscured. We performed additional analyses using a regional reference population and, in populations with a single founder, detected significant, positive SGS at all distances, indicating that these populations consist of highly related descendants and receive little seed immigration. Subsequent analyses of SGS in size cohorts in the four large study populations showed significant SGS in both juveniles and adults, probably because of a relative lack of intraspecific demographic thinning. SGS in populations of this colonizing tree is pronounced and persistent and is determined by the number and relatedness of founding individuals and adjacent seed sources. Patterns of SGS in populations with known histories may be used to indirectly infer possible colonization scenarios for populations where it is unknown.     

Social and spatial effects on genetic variation between foraging flocks in a wild bird population

Social interactions are rarely random. In some instances, animals exhibit homophily or heterophily, the tendency to interact with similar or dissimilar conspecifics, respectively. Genetic homophily and heterophily influence the evolutionary dynamics of populations, because they potentially affect sexual and social selection. Here, we investigate the link between social interactions and allele frequencies in foraging flocks of great tits (Parus major) over three consecutive years. We constructed co‐occurrence networks which explicitly described the splitting and merging of 85,602 flocks through time (fission–fusion dynamics), at 60 feeding sites. Of the 1,711 birds in those flocks, we genotyped 962 individuals at 4,701 autosomal single nucleotide polymorphisms (SNPs). By combining genomewide genotyping with repeated field observations of the same individuals, we were able to investigate links between social structure and allele frequencies at a much finer scale than was previously possible. We explicitly accounted for potential spatial effects underlying genetic structure at the population level. We modelled social structure and spatial configuration of great tit fission–fusion dynamics with eigenvector maps. Variance partitioning revealed that allele frequencies were strongly affected by group fidelity (explaining 27%–45% of variance) as individuals tended to maintain associations with the same conspecifics. These conspecifics were genetically more dissimilar than expected, shown by genomewide heterophily for pure social (i.e., space‐independent) grouping preferences. Genomewide homophily was linked to spatial configuration, indicating spatial segregation of genotypes. We did not find evidence for homophily or heterophily for putative socially relevant candidate genes or any other SNP markers. Together, these results demonstrate the importance of distinguishing social and spatial processes in determining population structure.     

Population genetic structure,genetic diversity,and natural history of the South American species of Nothofagus subgenus Lophozonia (Nothofagaceae) inferred from nuclear microsatellite data

The effect of glaciation on the levels and patterns of genetic variation has been well studied in the Northern Hemisphere. However, although glaciation has undoubtedly shaped the genetic structure of plants in the Southern Hemisphere, fewer studies have characterized the effect, and almost none of them using microsatellites. Particularly, complex patterns of genetic structure might be expected in areas such as the Andes, where both latitudinal and altitudinal glacial advance and retreat have molded modern plant communities. We therefore studied the population genetics of three closely related, hybridizing species of Nothofagus (N. obliqua, N. alpina, and N. glauca, all of subgenus Lophozonia; Nothofagaceae) from Chile. To estimate population genetic parameters and infer the influence of the last ice age on the spatial and genetic distribution of these species, we examined and analyzed genetic variability at seven polymorphic microsatellite DNA loci in 640 individuals from 40 populations covering most of the ranges of these species in Chile. Populations showed no significant inbreeding and exhibited relatively high levels of genetic diversity (H_E = 0.502–0.662) and slight, but significant, genetic structure (R_ST = 8.7–16.0%). However, in N. obliqua, the small amount of genetic structure was spatially organized into three well‐defined latitudinal groups. Our data may also suggest some introgression of N. alpina genes into N. obliqua in the northern populations. These results allowed us to reconstruct the influence of the last ice age on the genetic structure of these species, suggesting several centers of genetic diversity for N. obliqua and N. alpina, in agreement with the multiple refugia hypothesis.     

Costs and benefits of admixture between foreign genotypes and local populations in the field

Admixture is the hybridization between populations within one species. It can increase plant fitness and population viability by alleviating inbreeding depression and increasing genetic diversity. However, populations are often adapted to their local environments and admixture with distant populations could break down local adaptation by diluting the locally adapted genomes. Thus, admixed genotypes might be selected against and be outcompeted by locally adapted genotypes in the local environments. To investigate the costs and benefits of admixture, we compared the performance of admixed and within‐population F1 and F2 generations of the European plant Lythrum salicaria in a reciprocal transplant experiment at three European field sites over a 2‐year period. Despite strong differences between site and plant populations for most of the measured traits, including herbivory, we found limited evidence for local adaptation. The effects of admixture depended on experimental site and plant population, and were positive for some traits. Plant growth and fruit production of some populations increased in admixed offspring and this was strongest with larger parental distances. These effects were only detected in two of our three sites. Our results show that, in the absence of local adaptation, admixture may boost plant performance, and that this is particularly apparent in stressful environments. We suggest that admixture between foreign and local genotypes can potentially be considered in nature conservation to restore populations and/or increase population viability, especially in small inbred or maladapted populations.     

Association between host's genetic diversity and parasite burden in damselflies

Recent research indicates that low genetic variation in individuals can increase susceptibility to parasite infection, yet evidence from natural invertebrate populations remains scarce. Here, we studied the relationship between genetic heterozygosity, measured as AFLP‐based inbreeding coefficient f_AFLP, and gregarine parasite burden from eleven damselfly, Calopteryx splendens, populations. We found that in the studied populations, 5–92% of males were parasitized by endoparasitic gregarines (Apicomplexa: Actinocephalidae). Number of parasites ranged from none to 47 parasites per male, and parasites were highly aggregated in a few hosts. Mean individual f_AFLP did not differ between populations. Moreover, we found a positive association between individual's inbreeding coefficient and parasite burden. In other words, the more homozygous the individual, the more parasites it harbours. Thus, parasites are likely to pose strong selection pressure against inbreeding and homozygosity. Our results support the heterozygosity‐fitness correlation hypothesis, which suggests the importance of heterozygosity for an individual's pathogen resistance.     

Spatial and temporal genetic structure of Symbiodinium populations within a common reef‐building coral on the Great Barrier Reef

The dinoflagellate photosymbiont Symbiodinium plays a fundamental role in defining the physiological tolerances of coral holobionts, but little is known about the dynamics of these endosymbiotic populations on coral reefs. Sparse data indicate that Symbiodinium populations show limited spatial connectivity; however, no studies have investigated temporal dynamics for in hospite Symbiodinium populations following significant mortality and recruitment events in coral populations. We investigated the combined influences of spatial isolation and disturbance on the population dynamics of the generalist Symbiodinium type C2 (ITS1 rDNA) hosted by the scleractinian coral Acropora millepora in the central Great Barrier Reef. Using eight microsatellite markers, we genotyped Symbiodinium in a total of 401 coral colonies, which were sampled from seven sites across a 12‐year period including during flood plume–induced coral bleaching. Genetic differentiation of Symbiodinium was greatest within sites, explaining 70–86% of the total genetic variation. An additional 9–27% of variation was explained by significant differentiation of populations among sites separated by 0.4–13 km, which is consistent with low levels of dispersal via water movement and historical disturbance regimes. Sampling year accounted for 6–7% of total genetic variation and was related to significant coral mortality following severe bleaching in 1998 and a cyclone in 2006. Only 3% of the total genetic variation was related to coral bleaching status, reflecting generally small (8%) reductions in allelic diversity within bleached corals. This reduction probably reflected a loss of genotypes in hospite during bleaching, although no site‐wide changes in genetic diversity were observed. Combined, our results indicate the importance of disturbance regimes acting together with limited oceanographic transport to determine the genetic composition of Symbiodinium types within reefs.     

Plant‐mediated interactions between two herbivores differentially affect a subsequently arriving third herbivore in populations of wild cabbage

• Plants are part of biodiverse communities and frequently suffer from attack by multiple herbivorous insects. Plant responses to these herbivores are specific for insect feeding guilds: aphids and caterpillars induce different plant phenotypes. Moreover, plants respond differentially to single or dual herbivory, which may cascade into a chain of interactions in terms of resistance to other community members. Whether differential responses to single or dual herbivory have consequences for plant resistance to yet a third herbivore is unknown.
• We assessed the effects of single or dual herbivory by Brevicoryne brassicae aphids and/or Plutella xylostella caterpillars on resistance of plants from three natural populations of wild cabbage to feeding by caterpillars of Mamestra brassicae. We measured plant gene expression and phytohormone concentrations to illustrate mechanisms involved in induced responses.
• Performance of both B. brassicae and P. xylostella was reduced when feeding simultaneously with the other herbivore, compared to feeding alone. Gene expression and phytohormone concentrations in plants exposed to dual herbivory were different from those found in plants exposed to herbivory by either insect alone. Plants previously induced by both P. xylostella and B. brassicae negatively affected growth of the subsequently arriving M. brassicae. Furthermore, induced responses varied between wild cabbage populations.
• Feeding by multiple herbivores differentially activates plant defences, which has plant‐mediated negative consequences for a subsequently arriving herbivore. Plant population‐specific responses suggest that plant populations adapt to the specific communities of insect herbivores. Our study contributes to the understanding of plant defence plasticity in response to multiple insect attacks.

     

Long‐distance dispersal in a fire‐ and livestock‐protected savanna

Savannas are highly diverse and dynamic environments that can shift to forest formations due to protection policies. Long‐distance dispersal may shape the genetic structure of these new closed forest formations. We analyzed eight microsatellite loci using a single‐time approach to understand contemporary pollen and effective seed dispersal of the tropical tree, Copaifera langsdorffii Desf. (Fabaceae), occurring in a Brazilian fire‐ and livestock‐protected savanna. We sampled all adult trees found within a 10.24 ha permanent plot, young trees within a subplot of 1.44 ha and open‐pollinated seeds. We detected a very high level of genetic diversity among the three generations in the studied plot. Parentage analysis revealed high pollen immigration rate (0.64) and a mean contemporary pollen dispersal distance of 74 m. In addition, half‐sib production was 1.8 times higher than full‐sibs in significant higher distances, indicating foraging activity preference for different trees at long distances. There was a significant and negative correlation between diameter at breast height (DBH) of the pollen donor with the number of seeds (r = ?0.640, P‐value = 0.032), suggesting that pollen donor trees with a higher DBH produce less seeds. The mean distance of realized seed dispersal (recruitment kernel) was 135 m due to the large home range dispersers (birds and mammals) in the area. The small magnitude of spatial genetic structure found in young trees may be a consequence of overlapping seed shadows and increased tree density. Our results show the positive side of closed canopy expansion, where animal activities regarding pollination and seed dispersal are extremely high.     

Seed dispersal and realized gene flow of two forest orchids in a fragmented landscape

• Species with vast production of dust‐like windborne seeds, such as orchids, should not be limited by seed dispersal. This paradigm, however, does not fit recent studies showing that many sites suitable for orchids are unoccupied and most seeds land close to their maternal plant. To explore this issue, we studied seed dispersal and gene flow of two forest orchid species, Epipactis atrorubens and Cephalanthera rubra, growing in a fragmented landscape of forested limestone hills in southwest Bohemia, Czech Republic.
• We used a combination of seed trapping and plant genotyping methods (microsatellite DNA markers) to quantify short‐ and long‐distance dispersal, respectively. In addition, seed production of both species was estimated.
• We found that most seeds landed very close to maternal plants (95% of captured seeds were within 7.2 m) in both species, and dispersal distance was influenced by forest type in E. atrorubens. In addition, C. rubra showed clonal reproduction (20% of plants were of clonal origin) and very low fruiting success (only 1.6% of plants were fruiting) in comparison with E. atrorubens (25.7%). Gene flow was frequent up to 2 km in C. rubra and up to 125 km in E. atrorubens, and we detected a relatively high dispersal rate among regions in both species.
• Although both species occupy similar habitats and have similar seed dispersal abilities, C. rubra is notably rarer in the study area. Considerably low fruiting success in this species likely limits its gene flow to longer distances and designates it more sensitive to habitat loss and fragmentation.

     

Postglacial colonization history reflects in the genetic structure of natural populations of Festuca rubra in Europe

We conducted a large‐scale population genetic survey of genetic diversity of the host grass Festuca rubra s.l., which fitness can be highly dependent on its symbiotic fungus Epichloë festucae, to evaluate genetic variation and population structure across the European range. The 27 studied populations have previously been found to differ in frequencies of occurrence of the symbiotic fungus E. festucae and ploidy levels. As predicted, we found decreased genetic diversity in previously glaciated areas in comparison with nonglaciated regions and discovered three major maternal genetic groups: southern, northeastern, and northwestern Europe. Interestingly, host populations from Greenland were genetically similar to those from the Faroe Islands and Iceland, suggesting gene flow also between those areas. The level of variation among populations within regions is evidently highly dependent on the postglacial colonization history, in particular on the number of independent long‐distance seed colonization events. Yet, also anthropogenic effects may have affected the population structure in F. rubra. We did not observe higher fungal infection rates in grass populations with lower levels of genetic variability. In fact, the fungal infection rates of E. festucae in relation to genetic variability of the host populations varied widely among geographical areas, which indicate differences in population histories due to colonization events and possible costs of systemic fungi in harsh environmental conditions. We found that the plants of different ploidy levels are genetically closely related within geographic areas indicating independent formation of polyploids in different maternal lineages.