Vertical and horizontal photobiont transmission within populations of a lichen symbiosis

Lichens are widespread symbioses and play important roles in many terrestrial ecosystems. The genetic structure of lichens is the result of the association between fungal and algal populations constituting the lichen thallus. Using eight fungus- and seven alga-specific highly variable microsatellite markers on within-population spatial genetic data from 62 replicate populations across Europe, North America, Asia and Africa, we investigated the contributions of vertical and horizontal transmission of the photobiont to the genetic structure of the epiphytic lichen Lobaria pulmonaria. Based on pairwise comparisons of multilocus genotypes defined separately for the mycobiont and for the photobiont, we inferred the transmission mode of the photobiont and the relative contribution of somatic mutation and recombination. After constraining the analysis of one symbiont to pairs of individuals with genetically identical symbiotic partners, we found that 77% of fungal and 70% of algal pairs were represented by clones. Thus, the predominant dispersal mode was by means of symbiotic vegetative propagules (vertical transmission), which dispersed fungal and algal clones co-dependently over a short distance, thus shaping the spatial genetic structure up to distances of 20m. Evidence for somatic mutation generating genetic diversity was found in both symbionts, accounting for 30% of pairwise comparisons in the alga and 15% in the fungus. While the alga did not show statistically significant evidence of recombination, recombination accounted for 7.7% of fungal pairs with identical algae. This implies that, even in a mostly vegetatively reproducing species, horizontal transmission plays a role in shaping the symbiotic association, as shown in many coral and other symbioses in nature.     

Effects of multiple founder populations on spatial genetic structure of reintroduced American martens

Reintroductions and translocations are increasingly used to repatriate or increase probabilities of persistence for animal and plant species. Genetic and demographic characteristics of founding individuals and suitability of habitat at release sites are commonly believed to affect the success of these conservation programs. Genetic divergence among multiple source populations of American martens (Martes americana) and well documented introduction histories permitted analyses of post‐introduction dispersion from release sites and development of genetic clusters in the Upper Peninsula (UP) of Michigan <50 years following release. Location and size of spatial genetic clusters and measures of individual‐based autocorrelation were inferred using 11 microsatellite loci. We identified three genetic clusters in geographic proximity to original release locations. Estimated distances of effective gene flow based on spatial autocorrelation varied greatly among genetic clusters (30–90 km). Spatial contiguity of genetic clusters has been largely maintained with evidence for admixture primarily in localized regions, suggesting recent contact or locally retarded rates of gene flow. Data provide guidance for future studies of the effects of permeabilities of different land‐cover and land‐use features to dispersal and of other biotic and environmental factors that may contribute to the colonization process and development of spatial genetic associations.     

Analysis of spatial genetic structure in an expanding Pinus halepensis population reveals development of fine-scale genetic clustering over time

We analysed the change of spatial genetic structure (SGS) of reproductive individuals over time in an expanding Pinus halepensis population. To our knowledge, this is the first empirical study to analyse the temporal component of SGS by following the dynamics of successive cohorts of the same population over time, rather than analysing different age cohorts at a single time. SGS is influenced by various factors including restricted gene dispersal, microenvironmental selection, mating patterns and the spatial pattern of reproductive individuals. Several factors that affect SGS are expected to vary over time and as adult density increases. Using air photo analysis, tree-ring dating and molecular marker analysis we reconstructed the spread of reproductive individuals over 30 years beginning from five initial individuals. In the early stages, genotypes were distributed randomly in space. Over time and with increasing density, fine-scale (< 20 m) SGS developed and the magnitude of genetic clustering increased. The SGS was strongly affected by the initial spatial distribution and genetic variation of the founding individuals. The development of SGS may be explained by fine-scale environmental heterogeneity and possibly microenvironmental selection. Inbreeding and variation in reproductive success may have enhanced SGS magnitude over time.     

Molecular genetic diversity in populations of the stingless bee Plebeia remota: A case study

Genetic diversity is a major component of the biological diversity of an ecosystem. The survival of a population may be seriously threatened if its genetic diversity values are low. In this work, we measured the genetic diversity of the stingless bee Plebeia remota based on molecular data obtained by analyzing 15 microsatellite loci and sequencing two mitochondrial genes. Population structure and genetic diversity differed depending on the molecular marker analyzed: microsatellites showed low population structure and moderate to high genetic diversity, while mitochondrial DNA (mtDNA) showed high population structure and low diversity in three populations. Queen philopatry and male dispersal behavior are discussed as the main reasons for these findings.     

Non-random distribution of individual genetic diversity along an environmental gradient

Improving our knowledge of the links between ecology and evolution is especially critical in the actual context of global rapid environmental changes. A critical step in that direction is to quantify how variation in ecological factors linked to habitat modifications might shape observed levels of genetic variability in wild populations. Still, little is known on the factors affecting levels and distribution of genetic diversity at the individual level, despite its vital underlying role in evolutionary processes. In this study, we assessed the effects of habitat quality on population structure and individual genetic diversity of tree swallows (Tachycineta bicolor) breeding along a gradient of agricultural intensification in southern Québec, Canada. Using a landscape genetics approach, we found that individual genetic diversity was greater in poorer quality habitats. This counter-intuitive result was partly explained by the settlement patterns of tree swallows across the landscape. Individuals of higher genetic diversity arrived earlier on their breeding grounds and settled in the first available habitats, which correspond to intensive cultures. Our results highlight the importance of investigating the effects of environmental variability on individual genetic diversity, and of integrating information on landscape structure when conducting such studies.     

Population structure, genetic diversity, and clone formation in Quercus chrysolepis (Fagaceae)

Stands of canyon live oak (Quercus chrysolepis, Fagaceae) are maintained for fuelwood, fire management, recreation, and as habitat for wildlife. Information about the link between the oak's reproductive ecology and its extent of genetic diversity is important in developing land management policies that will maintain the long-term viability of populations. Basal sprouting is the primary means of reproduction following fire or cutting, and stands frequently include groups of visibly connected trees in a clustered distribution that suggests cloning. We determined the extent to which clusters of trees were clonal and defined the spatial pattern and diversity of genotypes for six populations across nearly the entire east-west extent of the San Bernardino Mountains in southern California. We mapped over 100 trees at each of five sites and genotyped each tree for allozymes at seven polymorphic loci. We identified clones using these multilocus genotypes and detected an average of 34.4 ± 7.3 (SD) clones per site, most of which had unique genotypes. In general, clustered trees belong to single clones and most clones consist of few trees (mean = 3.4 ± 0.6 trees per clone). However, clone size increased significantly with increased individual heterozygosity, suggesting that selection may favor highly heterozygous clones. Clonal diversity and evenness were high relative to reports for most other clonal species; an average of 97% of clones had distinct genotypes, and Simpson's index of diversity averaged 0.95 ± 0.02. Population genetic analyses of 319 clones from six sites revealed high genetic diversity within sites (mean HS = 0.443). Only a small proportion of the total genetic diversity was explained by variation among sites (mean GST = 0.018), which is consistent with high gene flow among sites (Nm = 9.5). We found no significant substructure among plots within sites, and fixation indices within sites were generally small, suggesting that either little inbreeding occurs, and/or few inbred progeny survive. However, spatial autocorrelation analysis of clones indicated fine-scale genetic structure at distances under 4 m, possibly due to limited seed dispersal. Our data suggest that guidelines for seed collection of canyon live oak for use in restoration can be specified in a manner similar to that recommended for conifer species within the region studied.     

Extremely low effective population sizes, genetic structuring and reduced genetic diversity in a threatened bumblebee species, Bombus sylvarum (Hymenoptera: Apidae)

Habitat fragmentation may severely affect survival of social insect populations as the number of nests per population, not the number of individuals, represents population size, hence they may be particularly prone to loss of genetic diversity. Erosion of genetic diversity may be particularly significant among social Hymenoptera such as bumblebees (Bombus spp.), as this group may be susceptible to diploid male production, a suggested direct cost of inbreeding. Here, for the first time, we assess genetic diversity and population structuring of a threatened bumblebee species (Bombus sylvarum) which exists in highly fragmented habitat (rather than oceanic) islands. Effective population sizes, estimated from identified sisterhoods, were very low (range 21-72) suggesting that isolated populations will be vulnerable to loss of genetic variation through drift. Evidence of significant genetic structuring between populations (theta = 0.084) was found, but evidence of a bottleneck was detected in only one population. Comparison across highly fragmented UK populations and a continental population (where this species is more widespread) revealed significant differences in allelic richness attributable to a high degree of genetic diversity in the continental population. While not directly related to population size, this is perhaps explained by the high degree of isolation between UK populations relative to continental populations. We suggest that populations now existing on isolated habitat islands were probably linked by stepping-stone populations prior to recent habitat loss.     

Mollusc grazing may constrain the ecological niche of the old forest lichen Pseudocyphellaria crocata

This study reports on mollusc grazing of two epiphytic cyanobacterial lichens (Pseudocyphellaria crocata and Lobaria pulmonaria) transplanted within three Picea abies-dominated boreal rain forest stands (clear-cut, young and old forests) in west central Norway. Grazing was particularly high in transplants located in the old forest and was almost absent in clear-cut transplants. Grazing marks were absent on natural thalli on nearby spruce twigs (required creeping distance for mollusc from the ground >4 m). Transplantation of lichens from twigs to artificial transplantation frames reduced the creeping distance to 1.2 m, and caused a significant increase in grazing damage in P. crocata. Given a paired choice under transplantation, molluscs consistently preferred P. crocata and avoided L. pulmonaria, implying species-specific differences in herbivore defence. Pseudocyphellaria crocata has a much lower content of the medullary depsidones stictic and constictic acid than L. pulmonaria. Heavy grazing occurred in the P. crocata thalli lowest in these two depsidones. The upper part of the medulla hosting the photobiont was the preferred fodder for grazing molluscs. Molluscs avoided the yellow soralia in P. crocata (localised pulvinic acid), suggesting a role for pulvinic acid in preventing grazing of detached soredia and early establishment stages. The preference of P. crocata for thin spruce twigs is probably a result of a lower grazing pressure on twigs compared to e.g. deciduous stems that frequently support the better defended L. pulmonaria. Ongoing climate changes with increased annual rainfall and milder winters have presumably increased mollusc grazing, particularly in SW parts of Norway which have more species of lichen-feeding molluscs than the boreal sites studied. These temperate areas lacking natural spruce populations have recently experienced reported extinctions of the poorly defended P. crocata from rocks and deciduous stems prone to mollusc grazing. Lichen-feeding molluscs have likely played a role in these extinctions, causing spruce twigs in Atlantic boreal forests to be a last strong foothold for P. crocata in Scandinavia.     

Variable extent of sex‐biased dispersal in a strongly polygynous mammal

For mammals with a polygynous mating system, dispersal is expected to be male‐biased. However, with the increase in empirical studies, discrepancies are arising between the expected and observed direction/extent of the bias in dispersal. In this study, we assessed sex‐biased dispersal in red deer (Cervus elaphus) on 13 estates from the Scottish Highlands. A total of 568 adult individuals were genotyped at 21 microsatellite markers and sequenced for 821 bp of the mitochondrial control region. Estimates of population structure with mitochondrial sequences were eight times larger than that obtained with microsatellite data (F_{st′‐mt DNA} = 0.831; F_{st′‐micros} = 0.096) indicating overall male‐biased dispersal in the study area. Comparisons of microsatellite data between the sexes indicated a predominance of male‐biased dispersal in the study area but values of F_ST and relatedness were only slighter larger for females. Individual‐based spatial autocorrelation analysis generated a similar pattern of relatedness across geographical distances for both sexes, with differences only significant at two distance intervals (25–30 and 70–112 km). Patterns of relatedness differed between estates, male biased‐dispersal was detected in eight estates but no sex‐biased dispersal was found in the remaining five. Neither population density nor landscape cover was found to be associated with the patterns of relatedness found across the estates. Differences in management strategies that could influence age structure, sex ratio and dispersal behaviour are proposed as potential factors influencing the relatedness patterns observed. This study provides new insights on dispersal of a strongly polygynous mammal at geographical scales relevant for management and conservation.     

Contrasting genetic structures of two parasitic nematodes,determined on the basis of neutral microsatellite markers and selected anthelmintic resistance markers

For the first time, the neutral genetic relatedness of natural populations of Trichostrongylid nematodes was investigated in relation to polymorphism of the β‐tubulin gene, which is selected for anthelminthic treatments. The aim of the study was to assess the contribution of several evolutionary processes: migration and genetic drift by neutral genetic markers and selection by anthelminthic treatments on the presence of resistance alleles at β‐tubulin. We studied two nematode species (Teladorsagia circumcincta and Haemonchus contortus) common in temperate climatic zones; these species are characterized by contrasting life history traits. We studied 10 isolated populations of goat nematode parasites: no infected adult goat had been exchanged after the herds were established. Beta‐tubulin polymorphism was similar in these two species. One and two β‐tubulin alleles from T. circumcincta and H. contortus respectively were shared by several populations. Most of the β‐tubulin alleles were ‘private’ alleles. No recombination between alleles was detected in BZ‐resistant alleles from T. circumcincta and H. contortus. The T. circumcincta populations have not diverged much since their isolation (F_ST <0.08), whereas H. contortus displayed marked local genetic differentiation (F_ST ranging from 0.08 to 0.18). These findings suggest that there are severe bottlenecks in the H. contortus populations, possibly because of their reduced abundance during unfavourable periods and their high reproductive rate, which allows the species to persist even after severe population reduction. Overall, the data reported contradict the hypothesis of the origin of β‐tubulin resistance alleles in these populations from a single mutational event, but two other hypotheses (recurrent mutation generating new alleles in isolated populations and the introduction of existing alleles) emerge as equally likely.     

Stability of population structure and genetic diversity across generations assessed by microsatellites among sympatric populations of landlocked Atlantic salmon (Salmo salar L.)

It may often be necessary to perform genetic analyses of temporal replicates to estimate the significance of spatial variation independently from that of temporal variation in order to ensure the reliability of estimates of a defined population structure. Nevertheless, temporal studies of genetic diversity remain scarce in the literature relative to the plethora of empirical studies of population structure. In vertebrates, a limited number of studies have specifically assessed the temporal stability of population structure for more than one generation. In this study, we performed a microsatellite analysis of DNA obtained from archived scales to compare the population structure among four sympatric landlocked populations of Atlantic salmon ( Salmo salar ) over a time frame of three to five generations. The same patterns of allele frequency distribution, θ, R _ST and genetic distance estimates were observed among populations for two time periods, confirming the temporal stability of the population structure. Despite population declines and stocking during this period, no statistically significant changes in intrapopulation genetic diversity were apparent. This study illustrates the feasibility and usefulness of microsatellite analysis of temporal samples, not only to infer changes of intrapopulation genetic diversity, but also to assess the stability of population structure over a time frame of several generations.     

Colonization history and genetic diversity: adaptive potential in early stage invasions

The introduction of Anolis cristatellus from the multiple species anole community of Puerto Rico in the Greater Antilles to the island of Dominica in the Lesser Antilles, with its solitary endemic anole, provides an example of a very recent, timed, single colonization. We investigate the geographic origin and adaptive potential of the Dominican population using a range of methods including mtDNA phylogeography, nuclear microsatellite variation and multiple paternity studies, as well as heritability estimates, common garden experiments and comparative geographic studies of quantitative scalation traits. Phylogeographic analysis of NADH2 and microsatellite studies suggests that the Dominican population arose from a set of individuals from the central west area of Puerto Rico within their endemic range. The multiple‐individual inoculation, together with sperm storage and evidence of multiple paternity indicate genetic variability and suggest the potential for adaptation by natural selection. Estimates of heritability, common garden experiments and broad sense Q_ST/F_ST ratios, linked to replicated comparisons along elevational transects go some way to suggesting that the invasive populations may be adapting by natural selection, in parallel with the endemic anole, in the brief period since their introduction.     

Spatial genetic structure of a small rodent in a heterogeneous landscape

Gene flow in natural populations may be strongly influenced by landscape features. The integration of landscape characteristics in population genetic studies may thus improve our understanding of population functioning. In this study, we investigated the population genetic structure and gene flow pattern for the common vole, Microtus arvalis, in a heterogeneous landscape characterised by strong spatial and temporal variation. The studied area is an intensive agricultural zone of approximately 500 km² crossed by a motorway. We used individual-based Bayesian methods to define the number of population units and their spatial borders without prior delimitation of such units. Unexpectedly, we determined a single genetic unit that covered the entire area studied. In particular, the motorway considered as a likely barrier to dispersal was not associated with any spatial genetic discontinuity. Using computer simulations, we demonstrated that recent anthropogenic barriers to effective dispersal are difficult to detect through analysis of genetic variation for species with large effective population sizes. We observed a slight, but significant, pattern of isolation by distance over the whole study site. Spatial autocorrelation analyses detected genetic structuring on a local scale, most probably due to the social organisation of the study species. Overall, our analysis suggests intense small-scale dispersal associated with a large effective population size. High dispersal rates may be imposed by the strong spatio-temporal heterogeneity of habitat quality, which characterises intensive agroecosystems.     

Reduced genetic diversity, increased isolation and multiple introductions of invasive giant hogweed in the western Swiss Alps

The giant hogweed ( Heracleum mantegazzianum ) has successfully invaded 19 European countries as well as parts of North America. It has become a problematic species due to its ability to displace native flora and to cause public health hazards. Applying population genetics to species invasion can help reconstruct invasion history and may promote more efficient management practice. We thus analysed levels of genetic variation and population genetic structure of H. mantegazzianum in an invaded area of the western Swiss Alps as well as in its native range (the Caucasus), using eight nuclear microsatellite loci together with plastid DNA markers and sequences. On both nuclear and plastid genomes, native populations exhibited significantly higher levels of genetic diversity compared to invasive populations, confirming an important founder event during the invasion process. Invasive populations were also significantly more differentiated than native populations. Bayesian clustering analysis identified five clusters in the native range that corresponded to geographically and ecologically separated groups. In the invaded range, 10 clusters occurred. Unlike native populations, invasive clusters were characterized by a mosaic pattern in the landscape, possibly caused by anthropogenic dispersal of the species via roads and direct collection for ornamental purposes. Lastly, our analyses revealed four main divergent groups in the western Swiss Alps, likely as a consequence of multiple independent establishments of H. mantegazzianum .     

A review of techniques for spatial modeling in geographical, conservation and landscape genetics

Most evolutionary processes occur in a spatial context and several spatial analysis techniques have been employed in an exploratory context. However, the existence of autocorrelation can also perturb significance tests when data is analyzed using standard correlation and regression techniques on modeling genetic data as a function of explanatory variables. In this case, more complex models incorporating the effects of autocorrelation must be used. Here we review those models and compared their relative performances in a simple simulation, in which spatial patterns in allele frequencies were generated by a balance between random variation within populations and spatially-structured gene flow. Notwithstanding the somewhat idiosyncratic behavior of the techniques evaluated, it is clear that spatial autocorrelation affects Type I errors and that standard linear regression does not provide minimum variance estimators. Due to its flexibility, we stress that principal coordinate of neighbor matrices (PCNM) and related eigenvector mapping techniques seem to be the best approaches to spatial regression. In general, we hope that our review of commonly used spatial regression techniques in biology and ecology may aid population geneticists towards providing better explanations for population structures dealing with more complex regression problems throughout geographic space.     

Genetic consequences of an invasion through a patchy environment — the cynipid gallwasp Andricus quercuscalicis (Hymenoptera: Cynipidae)

Over the last 300–400 years, the cynipid gallwasp Andricus quercuscalicis has invaded northern and western Europe following human introduction of an obligate host plant, the Turkey oak ( Quercus cerris ) from south-eastern Europe. In the introduced range, distances between Turkey oak patches are greater and the numbers of oaks in each patch are far lower than in its native range. These changes in spatial distribution of Turkey oak are predicted to result in high genetic subdivision of A. quercuscalicis in invaded areas relative to its native range. Allozyme electrophoresis was used to examine genetic variation in 823 gall wasps from 39 populations of A. quercuscalicis. No new electrophoretic variants were found in the invaded range and both allelic diversity and mean heterozygosity decreased significantly with distance from the native range. Spatial autocorrelation analysis and values of Wright's F_st indicate that differences in allele frequences between populations were substantially greater in the invaded range than in the native range. Spatial autocorrelation analysis also suggests that changes in allele frequencies across Europe are unlikely to be the results of selection, but rather of strong directional migration followed by limited gene flow between populations. Patterns of genetic differentiation across Europe suggest that populations of A. quercuscalicis have been founded sequentially from the east through a process of random genetic subsampling and not by source-sink colonization directly from the native range.     

Approximate Bayesian computation reveals the factors that influence genetic diversity and population structure of foxsnakes

Contemporary geographical range and patterns of genetic diversity within species reflect complex interactions between multiple factors acting across spatial and temporal scales, and it is notoriously difficult to disentangle causation. Here, we quantify patterns of genetic diversity and genetic population structure using mitochondrial DNA sequences (101 individuals, cytochrome b) and microsatellites (816 individuals, 12 loci) and use Approximate Bayesian computation methods to test competing models of the demographic history of eastern and western foxsnakes. Our analyses indicate that post-glacial colonization and past population declines, probably caused by the infilling of deciduous forest and cooler temperatures since the mid-Holocene, largely underpin large-scale genetic patterns for foxsnakes. At finer geographical scales, our results point to more recent anthropogenic habitat loss as having accentuated genetic population structure by causing further declines and fragmentation.     

The patterns of genetic diversity in Carpentaria acuminata (Arecaceae), and rainforest history in northern Australia

Carpentaria acuminata occurs in monsoon rainforest and is endemic to the Northern Territory, Australia. The genetic diversity of C. acuminata populations was surveyed across the geographical range of the species using isozyme analysis. Genetic diversity within C. acuminata populations ( H _E = 0.143) was typical of rainforest species and woody angiosperms generally. Genetic diversity was not correlated with rainforest patch size. However, there was significant heterogeneity among populations ( F _ST = 0.379), with infrequent effective gene flow among populations ( Nm = 0.39). Genetic diversity was negatively correlated with increasing distance between neighbouring C. acuminata populations, but geographical distance was not a good predictor of genetic similarity. C. acuminata is a favoured food of mobile frugivores such as Torres Strait pigeons and flying foxes. The decreased diversity with decreasing density of populations indicated that seed dispersal by frugivores has been important for the maintenance of diversity in this species. Populations known to have originated on relatively young, Holocene landforms were not necessarily genetically depauperate. Gene flow by pollen is apparently limited because C. acuminata populations are significantly inbred regardless of genetic diversity ( F = 0.641). The distribution and diversity of rare alleles, i.e. those occurring in few populations, is consistent with the theory of rainforest contraction during the Pleistocene.