Mutation at high rates reduces spatial structure within populations

Genetic and epigenetic divergence between disturbed and undisturbed subpopulations of a Mediterranean shrub: a 20-year field experiment

Little is known on the potential of ecological disturbance to cause genetic and epigenetic changes in plant populations. We take advantage of a long‐term field experiment initiated in 1986 to study the demography of the shrub Lavandula latifolia, and compare genetic and epigenetic characteristics of plants in two adjacent subplots, one experimentally disturbed and one left undisturbed, 20 years after disturbance. Experimental setup was comparable to an unreplicated ‘Before‐After‐Control‐Impact’ (BACI) design where a single pair of perturbed and control areas were compared. When sampled in 2005, plants in the two subplots had roughly similar ages, but they had established in contrasting environments: dense conspecific population (‘Undisturbed’ subpopulation) versus open area with all conspecifics removed (‘Disturbed’ subpopulation). Plants were characterized genetically and epigenetically using amplified fragment length polymorphism (AFLP) and two classes of methylation‐sensitive AFLP (MSAP) markers. Subpopulations were similar in genetic diversity but differed in epigenetic diversity and multilocus genetic and epigenetic characteristics. Epigenetic divergence between subpopulations was statistically unrelated to genetic divergence. Bayesian clustering revealed an abrupt linear boundary between subpopulations closely coincident with the arbitrary demarcation line between subplots drawn 20 years back, which supports that genetic and epigenetic divergence between subpopulations was caused by artificial disturbance. There was significant fine‐scale spatial structuring of MSAP markers in both subpopulations, which in the Undisturbed one was indistinguishable from that of AFLP markers. Genetic differences between subpopulations could be explained by divergent selection alone, while the concerted action of divergent selection and disturbance‐driven appearance of new methylation variants in the Disturbed subpopulation is proposed to explain epigenetic differences. This study provides the first empirical evidence to date suggesting that relatively mild disturbances could leave genetic and epigenetic signatures on the next adult generation of long‐lived plants.

Microsatellites for microbiologists

Microsatellites are repeating sequences of 2-6base pairs of DNA. Currently, they are used as molecular markers in many organisms, specifically in genetic studies analyzing kinship and population structure. In addition, they can be used to study gene duplication and/or deletion. Although they are used in studies on microorganisms including fungi, bacteria, protists, and archaea, it appears that these genetic markers are not being utilized to their full microbiological potential. Microsatellites have many advantages over other genetic markers currently in use as they are in general species specific, and therefore, cross-contamination by nontarget organisms is rare. Furthermore, microsatellites are suitable for use with fast and cheap DNA extraction methods, with ancient DNA or DNA from hair and fecal samples used in noninvasive sampling, making them widely available as a genetic marker. Microsatellites have already proven to be a useful tool for evolutionary studies of pathogenic microorganisms such as Candida albicans and Helicobacter pylori, and the onset of new sequencing techniques (such as 454, PACBIO, and mini-ion sequencing) means the ability to detect such markers will become less time consuming and cheaper, thus further expanding their potential to answer important microbial ecology questions.

Comparative spatial genetics and epigenetics of plant populations: heuristic value and a proof of concept

Despite the recent upsurge of interest on natural epigenetic variation of nonmodel organisms, factors conditioning the spatial structure of epigenetic diversity in wild plant populations remain virtually unexplored. We propose that information on processes shaping natural epigenetic variation can be gained using the spatial structure of genetic diversity as null model. Departures of epigenetic isolation-by-distance (IBD) patterns from genetic IBD patterns for the same sample, particularly differences in slope of similarity-distance regressions, will reflect the action of factors that operate specifically on epigenetic variation, including imperfect transgenerational inheritance and responsiveness to environmental factors of epigenetic marks. As a proof of concept, we provide a comparative analysis of spatial genetic and epigenetic structure of 200 mapped individuals of the perennial herb Helleborus foetidus. Plants were fingerprinted using nuclear microsatellites, amplified fragment length polymorphisms (AFLP) and methylation-sensitive AFLP markers. Expectations from individual-level IBD patterns were tested by means of kinship-distance regressions. Both genetic and epigenetic similarity between H. foetidus individuals conformed to theoretical expectations under individual-level IBD models. Irrespective of marker type, there were significant negative linear relationships between the kinship coefficient for plant pairs and their spatial separation. Regression slopes were significantly steeper for epigenetic markers. Epigenetic similarity between individuals was much greater than genetic similarity at shortest distances, such epigenetic 'kinship excess' tending to decrease as plant separation increased. Results suggest that moderate-to-high heritability and responsiveness to local environments are major drivers of epigenetic spatial structure in H. foetidus, and illustrate the heuristic value of comparing genetic and epigenetic spatial structure for formulating and testing hypotheses on forces shaping epigenetic diversity in wild plant populations.

Are microsatellite fragment lengths useful for population-level studies? The case of Polygala lewtonii (Polygalaceae)

PREMISE OF THE STUDY

Microsatellites, despite being commonly used population-level markers, contain biases because scoring relies solely on fragment length. Their complexity can lead to homoplasy, the effects of which are poorly understood. We investigated the impact of using fragment lengths, repeats, or flanking region sequences on common population-level analyses.

METHODS

Five polymorphic microsatellite markers amplified across the central Florida scrub endemic Polygala lewtonii (Polygalaceae) and its close, widespread congener P. polygama. We genotyped 147 individuals of P. lewtonii and 156 of P. polygama, and sequenced the amplicons of four markers across all individuals. We ran basic statistics, spatial clustering analysis, historical demographics, and migration tests.

RESULTS

One population of intermediate morphology was genetically clearly identified as P. polygama, making it the southernmost population of this species. Statistics were comparable between the fragment length and repeat numbers, with some notable differences. Flanking regions exhibited surprisingly high polymorphism between species, and between geographically distant conspecific populations.

DISCUSSION

The increasing use of markers developed in one species and amplified in another is only a good practice if precautions are taken, notably the sequencing of the fragments between species and populations. Flanking region sequences are a useful marker at the interspecific level.

Assessing Genetic Structure in Common but Ecologically Distinct Carnivores: The Stone Marten and Red Fox

The identification of populations and spatial genetic patterns is important for ecological and conservation research, and spatially explicit individual-based methods have been recognised as powerful tools in this context. Mammalian carnivores are intrinsically vulnerable to habitat fragmentation but not much is known about the genetic consequences of fragmentation in common species. Stone martens (Martes foina) and red foxes (Vulpes vulpes) share a widespread Palearctic distribution and are considered habitat generalists, but in the Iberian Peninsula stone martens tend to occur in higher quality habitats. We compared their genetic structure in Portugal to see if they are consistent with their differences in ecological plasticity, and also to illustrate an approach to explicitly delineate the spatial boundaries of consistently identified genetic units. We analysed microsatellite data using spatial Bayesian clustering methods (implemented in the software BAPS, GENELAND and TESS), a progressive partitioning approach and a multivariate technique (Spatial Principal Components Analysis-sPCA). Three consensus Bayesian clusters were identified for the stone marten. No consensus was achieved for the red fox, but one cluster was the most probable clustering solution. Progressive partitioning and sPCA suggested additional clusters in the stone marten but they were not consistent among methods and were geographically incoherent. The contrasting results between the two species are consistent with the literature reporting stricter ecological requirements of the stone marten in the Iberian Peninsula. The observed genetic structure in the stone marten may have been influenced by landscape features, particularly rivers, and fragmentation. We suggest that an approach based on a consensus clustering solution of multiple different algorithms may provide an objective and effective means to delineate potential boundaries of inferred subpopulations. sPCA and progressive partitioning offer further verification of possible population structure and may be useful for revealing cryptic spatial genetic patterns worth further investigation.

Permeability of the landscape matrix between amphibian breeding sites

For organisms that reproduce in discrete habitat patches, land cover between patches (known as the matrix) is important for dispersal among breeding sites. Models of patchy populations often incorporate information on the permeability of the matrix to dispersal, sometimes based on expert opinion. I estimated the relative resistance to gene flow of land cover types and barriers using F(ST) calculated from microsatellite markers in two amphibians, within an 800-km(2) area in northern Switzerland. The species included a frog (Rana temporaria: 996 individuals, 48 populations, seven markers) and a newt (Triturus alpestris: 816 individuals, 41 populations, seven markers). Open fields and urban areas were more resistant to gene flow than forested land; roads and highways also reduced permeability. Results were similar for the two species. However, differences in resistance among matrix elements were relatively low: gene flow through urban areas was reduced by only 24-42% relative to forest; a divided highway reduced gene flow by 11-40% and was 7-8 times more resistant than a secondary road. These data offer an empirically based alternative to expert opinion for setting relative resistance values in landscape models.

Coastal pollution limits pelagic larval dispersal

The ecological impact of large coastal human populations on marine ecosystems remains relatively unknown. Here, we examine the population structure of Patiria miniata, the bat star, and correlate genetic distances with a model based on flow rates and proximity to P. miniata populations for the four major stormwater runoff and wastewater effluent sources of the Southern California Bight. We show that overall genetic connectivity is high (F(ST)~0.005); however, multivariate analyses show that genetic structure is highly correlated with anthropogenic inputs. The best models included both stormwater and wastewater variables and explained between 26.55 and 93.69% of the observed structure. Additionally, regressions between allelic richness and distance to sources show that populations near anthropogenic pollution have reduced genetic diversity. Our results indicate that anthropogenic runoff and effluent are acting as barriers to larval dispersal, effectively isolating a high gene flow species that is virtually free of direct human impact.

Microsatellites reveal substantial among-population genetic differentiation and strong inbreeding in the relict fern Dryopteris aemula

BACKGROUND AND AIMS

A previous study detected no allozyme diversity in Iberian populations of the buckler-fern Dryopteris aemula. The use of a more sensitive marker, such as microsatellites, was thus needed to reveal the genetic diversity, breeding system and spatial genetic structure of this species in natural populations.

METHODS

Eight microsatellite loci for D. aemula were developed and their cross-amplification with other ferns was tested. Five polymorphic loci were used to characterize the amount and distribution of genetic diversity of D. aemula in three populations from the Iberian Peninsula and one population from the Azores.

KEY RESULTS

Most microsatellite markers developed were transferable to taxa close to D. aemula. Overall genetic variation was low (H(T) = 0.447), but was higher in the Azorean population than in the Iberian populations of this species. Among-population genetic differentiation was high (F(ST) = 0.520). All loci strongly departed from Hardy-Weinberg equilibrium. In the population where genetic structure was studied, no spatial autocorrelation was found in any distance class.

CONCLUSIONS

The higher genetic diversity observed in the Azorean population studied suggested a possible refugium in this region from which mainland Europe has been recolonized after the Pleistocene glaciations. High among-population genetic differentiation indicated restricted gene flow (i.e. lack of spore exchange) across the highly fragmented area occupied by D. aemula. The deviations from Hardy-Weinberg equilibrium reflected strong inbreeding in D. aemula, a trait rarely observed in homosporous ferns. The absence of spatial genetic structure indicated effective spore dispersal over short distances. Additionally, the cross-amplification of some D. aemula microsatellites makes them suitable for use in other Dryopteris taxa.

Distinctively variable sequence-based nuclear DNA markers for multilocus phylogeography of the soybean- and rice-infecting fungal pathogen Rhizoctonia solani AG-1 IA

A series of multilocus sequence-based nuclear DNA markers was developed to infer the phylogeographical history of the Basidiomycetous fungal pathogen Rhizoctonia solani AG-1 IA infecting rice and soybean worldwide. The strategy was based on sequencing of cloned genomic DNA fragments (previously used as RFLP probes) and subsequent screening of fungal isolates to detect single nucleotide polymorphisms (SNPs). Ten primer pairs were designed based on these sequences, which resulted in PCR amplification of 200-320 bp size products and polymorphic sequences in all markers analyzed. By direct sequencing we identified both homokaryon and heterokaryon (i.e. dikaryon) isolates at each marker. Cloning the PCR products effectively estimated the allelic phase from heterokaryotic isolates. Information content varied among markers from 0.5 to 5.9 mutations per 100 bp. Thus, the former RFLP codominant probes were successfully converted into six distinctively variable sequence-based nuclear DNA markers. Rather than discarding low polymorphism loci, the combination of these distinctively variable anonymous nuclear markers would constitute an asset for the unbiased estimate of the phylogeographical parameters such as population sizes and divergent times, providing a more reliable species history that shaped the current population structure of R. solani AG-1 IA.

As the raven flies: using genetic data to infer the history of invasive common raven (Corvus corax) populations in the Mojave Desert

Common raven (Corvus corax) populations in Mojave Desert regions of southern California and Nevada have increased dramatically over the past five decades. This growth has been attributed to increased human development in the region, as ravens have a commensal relationship with humans and feed extensively at landfills and on road-killed wildlife. Ravens, as a partially subsidized predator, also represent a problem for native desert wildlife, in particular threatened desert tortoises (Gopherus agassizii). However, it is unclear whether the more than 15-fold population increase is due to in situ population growth or to immigration from adjacent regions where ravens have been historically common. Ravens were sampled for genetic analysis at several local sites within five major areas: the West Mojave Desert (California), East Mojave Desert (southern Nevada), southern coastal California, northern coastal California (Bay Area), and northern Nevada (Great Basin). Analyses of mtDNA control region sequences reveal an increased frequency of raven 'Holarctic clade' haplotypes from south to north inland, with 'California clade' haplotypes nearly fixed in the California populations. There was significant structuring among regions for mtDNA, with high F(ST) values among sampling regions, especially between the Nevada and California samples. Analyses of eight microsatellite loci reveal a mostly similar pattern of regional population structure, with considerably smaller, but mostly significant, values. The greater mtDNA divergences may be due to lower female dispersal relative to males, lower N(e), or effects of high mutation rates on maximal values of F(ST). Analyses indicate recent population growth in the West Mojave Desert and a bottleneck in the northern California populations. While we cannot rule out in situ population growth as a factor, patterns of movement inferred from our data suggest that the increase in raven populations in the West Mojave Desert resulted from movements from southern California and the Central Valley. Ravens in the East Mojave Desert are more similar to ones from northern Nevada, indicating movement between those regions. If this interpretation of high gene flow into the Mojave Desert is correct, then efforts to manage raven numbers by local control may not be optimally effective.

Distribution and fine-scale spatial-genetic structure in British wild cherry (Prunus avium L.)

Insights into the within-population spatial-genetic structure (SGS) of forest tree species, where little is known regarding seed and pollen dispersal patterns, enhance understanding of their ecology and provide information of value in conservation and breeding. This study utilised 13 polymorphic simple sequence repeat loci to investigate the impact of asexual recruitment, management regime and tree size on the development of SGS in wild cherry (Prunus avium L). Only 246 genotypes were identified in the 551 trees sampled, reflecting significant levels of clonal reproduction in both managed and unmanaged populations. Naturally regenerated wild cherry was spatially aggregated under both management regimes. However, in the managed population, sexually derived trees accounted for a greater proportion of the smaller size classes, whereas vegetatively produced trees dominated the smaller size classes in the unmanaged population. High overall SGS values (Sp 0.030-Sp 0.045) were observed when considering only sexually derived genets and kinship coefficients were significant up to the 120 m distance class for both populations. The inclusion of clonal ramets in the analysis significantly increased the overall SGS (Sp 0.089-Sp 0.119) as well as kinship coefficients in the 40-80 m distance classes, illustrating the dramatic impact of vegetative propagation on SGS in this species. Increased spatial aggregation and regeneration appeared to be concomitant with increased SGS in the 40 m distance class in the unmanaged population. Neighbourhood size estimates were relatively small for both populations and kinship coefficients were found to decline with distance under both management regimes, suggesting that common mechanisms may restrict gene dispersal in wild cherry.

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.

Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers

Recent improvements in genetic analysis and genotyping methods have resulted in a rapid expansion of the power of molecular markers to address ecological questions. Microsatellites have emerged as the most popular and versatile marker type for ecological applications. The rise of commercial services that can isolate microsatellites for new study species and genotype samples at reasonable prices presents ecologists with the unprecedented ability to employ genetic approaches without heavy investment in specialized equipment. Nevertheless, the lack of accessible, synthesized information on the practicalities and pitfalls of using genetic tools impedes ecologists' ability to make informed decisions on using molecular approaches and creates the risk that some will use microsatellites without understanding the steps needed to evaluate the quality of a genetic data set. The first goal of this synthesis is to provide an overview of the strengths and limitations of microsatellite markers and the risks, cost and time requirements of isolating and using microsatellites with the aid of commercial services. The second goal is to encourage the use and consistent reporting of thorough marker screening to ensure high quality data. To that end, we present a multistep screening process to evaluate candidate loci for inclusion in a genetic study that is broadly targeted to both novice and experienced geneticists alike.

Geographic differentiation in the house fly estimated by microsatellite and mitochondrial variation

Gene flow over very large geographic scales has been investigated in few species. Examples include Drosophila melanogaster, Drosophila subobscura, Drosophila simulans, and the Mediterranean fruit fly (Ceratitis capitata). The cosmopolitan house fly, a highly vagile, fecund, colonizing species offers an additional exemplar. Genotypes at seven microsatellite loci were scored in 14 widely separated natural house fly populations from the Nearctic, neotropics, Afrotropics, Palearctic, and Asia. Allelic diversities and heterozygosities differed significantly among populations. Averaged over all populations, Weir and Cockerham's theta = 0.13 and RST = 0.20. Pairwise genetic distance measures were uncorrelated with geographic distance. Microsatellite frequencies were compared with mitochondrial data from 13 of the same populations in which theta = 0.35 and Nei's GST = 0.72. Mitochondrial variation indicated up to threefold greater indices of genetic differentiation than the microsatellites. We were unable to draw any biogeographical inferences from these results or from tree or network topologies constructed from the genetic data. It is likely that high microsatellite diversities, mutation rates, and homoplasy greatly compromised their usefulness in estimating gene flow. House fly colonization dynamics include a large number of primary and secondary colonizations coupled with substantial genetic drift, but no detectable bottlenecks.

Estimating dispersal from short distance spatial autocorrelation

A series of theoretical studies has formed a strong connection between spatial statistics observed in populations and summary measures of the amount of dispersal. Synthesized, these developments allow dispersal to be indirectly estimated from standing spatial patterns of genetic variation under a range of conditions broad enough to be likely met in most populations of either plants or animals. The spatial correlations at the shortest distances are particularly robust to range of conditions and have disproportionately high statistical power. This review integrates theoretical results in a way that maximizes robustness and flexibility in the use of short distance autocorrelation to estimate Wright's neighborhood size, or the total variance in dispersal distances. Empirical guidelines are developed that are meant to be as practical and broad as possible. The guidelines focus on Moran's I-statistics for diploid genotypes converted to allele frequencies, but are also extended to or compared with several other approaches.