Application of high-throughput sequencing technologies to microsatellite genotyping (SSRseq) has been shown to remove many of the limitations of electrophoresis-based methods and to refine inference of population genetic diversity and structure. We present here a streamlined SSRseq development workflow that includes microsatellite development, multiplexed marker amplification and sequencing, and automated bioinformatics data analysis. We illustrate its application to five groups of species across phyla (fungi, plant, insect and fish) with different levels of genomic resource availability. We found that relying on previously developed microsatellite assay is not optimal and leads to a resulting low number of reliable locus being genotyped. In contrast, de novo ad hoc primer designs gives highly multiplexed microsatellite assays that can be sequenced to produce high quality genotypes for 20–40 loci. We highlight critical upfront development factors to consider for effective SSRseq setup in a wide range of situations. Sequence analysis accounting for all linked polymorphisms along the sequence quickly generates a powerful multi-allelic haplotype-based genotypic dataset, calling to new theoretical and analytical frameworks to extract more information from multi-nucleotide polymorphism marker systems.
This study assesses whether the effective number of breeders (N-b) can be estimated using a time and cost-effective protocol using genetic sibship reconstruction from a single sample of young-of-the-year (YOY) for the purposes of Atlantic salmon Salmo salar population monitoring. N-b was estimated for 10 consecutive reproductive seasons for S. salar in the River Nivelle, a small population located at the rear-edge of the species distribution area in France, chronically under its conservation limit and subjected to anthropogenic and environmental changes. Subsampling of real and simulated data showed that accurate estimates of N-b can be obtained from YOY genotypes, collected at moderate random sampling intensity, achievable using routine juvenile electrofishing protocols. Spatial bias and time elapsed since spawning were found to affect estimates, which must be accounted for in sampling designs. N-b estimated in autumn for S. salar in the River Nivelle was low and variable across years from 23 (95% C.I. 14-41) to 75 (53-101) and was not statistically correlated with the estimated number of returning adults, but it was positively correlated with the estimated number of YOY at age 9months. N-b was found to be lower for intermediate levels of redd aggregation, suggesting that the strength of the competition between males to access females affects reproductive success variance depending on redd spatial configuration. Thus, environmental factors such as habitat availability and quality for spawning and YOY development predominate over demographic ones (number of returning adults) in driving long-term population viability for S. salar in the River Nivelle. This study showcases N-b as an integrated parameter, encompassing demographic and ecological information about a reproductive event, relevant to the assessment of both short-term effects of management practices and long-term population conservation status.
BACKGROUND AND AIMS:Polyploidy in plants has been studied extensively. In many groups, two or more cytotypes represent separate biological entities with distinct distributions, histories and ecology. This study examines the distribution and origins of cytotypes of Alnus glutinosa in Europe, North Africa and western Asia.METHODS:A combined approach was used involving flow cytometry and microsatellite analysis of 12 loci in 2200 plants from 209 populations combined with species distribution modelling using MIROC and CCSM climatic models, in order to analyse (1) ploidy and genetic variation, (2) the origin of tetraploid A. glutinosa, considering A. incana as a putative parent, and (3) past distributions of the species.KEY RESULTS:The occurrence of tetraploid populations of A. glutinosa in Europe is determined for the first time. The distribution of tetraploids is far from random, forming two geographically well-delimited clusters located in the Iberian Peninsula and the Dinaric Alps. Based on microsatellite analysis, both tetraploid clusters are probably of autopolyploid origin, with no indication that A. incana was involved in their evolutionary history. A projection of the MIROC distribution model into the Last Glacial Maximum (LGM) showed that (1) populations occurring in the Iberian Peninsula and North Africa were probably interconnected during the LGM and (2) populations occurring in the Dinaric Alps did not exist throughout the last glacial periods, having retreated southwards into lowland areas of the Balkan Peninsula.CONCLUSIONS:Newly discovered tetraploid populations are situated in the putative main glacial refugia, and neither of them was likely to have been involved in the colonization of central and northern Europe after glacial withdrawal. This could mean that neither the Iberian Peninsula nor the western part of the Balkan Peninsula served as effective refugial areas for northward post-glacial expansion of A. glutinosa.
Understanding the consequences of habitat disturbance on mating patterns although pollen and seed dispersal in forest trees has been a long‐standing theme of forest and conservation genetics. Forest ecosystems face global environmental pressures from timber exploitation to genetic pollution and climate change, and it is therefore essential to comprehend how disturbances may alter the dispersal of genes and their establishment in tree populations in order to formulate relevant recommendations for sustainable resource management practices and realistic predictions of potential adaptation to climate change by means of range shift or expansion (Kremer et al. ). However, obtaining reliable evidence of disturbance‐induced effects on gene dispersal processes from empirical evaluation of forest tree populations is difficult. Indeed, tree species share characteristics such as high longevity, long generation time and large reproductive population size, which may impede the experimenter's ability to assess parameters at the spatial and time scales at which any change may occur (Petit and Hampe ). It has been suggested that appropriate study designs should encompass comparison of populations before and after disturbance as well as account for demonstrated variation in conspecific density, that is, the spatial distribution of mates, and forest density, including all species and relating to alteration in landscape openness (Bacles & Jump ). However, more often than not, empirical studies aiming to assess the consequences of habitat disturbance on genetic processes in tree populations assume rather than quantify a change in tree densities in forests under disturbance and generally fail to account for population history, which may lead to inappropriate interpretation of a causal relationship between population genetic structure and habitat disturbance due to effects of unmonitored confounding variables (Gauzere et al. 2013). In this issue, Shohami and Nathan () take advantage of the distinctive features of the fire‐adapted wind‐pollinated Aleppo pine Pinus halepensis (Fig. 1) to provide an elegant example of best practice. Thanks to long‐term monitoring of the study site, a natural stand in Israel, Shohami and Nathan witnessed the direct impact of habitat disturbance, here taking the shape of fire, on conspecific and forest densities and compared pre‐ and postdisturbance mating patterns estimated from cones of different ages sampled on the same surviving maternal individuals (Fig. 2). This excellent study design is all the more strong that Shohami and Nathan took further analytical steps to account for confounding variables, such as historical population genetic structure and possible interannual variation in wind conditions, thus giving high credibility to their findings of unequivocal fire‐induced alteration of mating patterns in P. halepensis. Most notably, the authors found, at the pollen pool level, a disruption of local genetic structure which, furthermore, they were able to attribute explicitly to enhanced pollen‐mediated gene immigration into the low‐density fire‐disturbed stand. This cleverly designed research provides a model approach to be followed if we are to advance our understanding of disturbance‐induced dispersal and genetic change in forest trees.
Predicting likely species responses to an alteration of their local environment is key to decision‐making in resource management, ecosystem restoration and biodiversity conservation practice in the face of global human‐induced habitat disturbance. This is especially true for forest trees which are a dominant life form on Earth and play a central role in supporting diverse communities and structuring a wide range of ecosystems. In Europe, it is expected that most forest tree species will not be able to migrate North fast enough to follow the estimated temperature isocline shift given current predictions for rapid climate warming. In this context, a topical question for forest genetics research is to quantify the ability for tree species to adapt locally to strongly altered environmental conditions (Kremer et al. ). Identifying environmental factors driving local adaptation is, however, a major challenge for evolutionary biology and ecology in general but is particularly difficult in trees given their large individual and population size and long generation time. Empirical evaluation of local adaptation in trees has traditionally relied on fastidious long‐term common garden experiments (provenance trials) now supplemented by reference genome sequence analysis for a handful of economically valuable species. However, such resources have been lacking for most tree species despite their ecological importance in supporting whole ecosystems. In this issue of Molecular Ecology, De Kort et al. () provide original and convincing empirical evidence of local adaptation to temperature in black alder, Alnus glutinosa L. Gaertn, a surprisingly understudied keystone species supporting riparian ecosystems. Here, De Kort et al. () use an innovative empirical approach complementing state‐of‐the‐art landscape genomics analysis of A. glutinosa populations sampled in natura across a regional climate gradient with phenotypic trait assessment in a common garden experiment (Fig. ). By combining the two methods, De Kort et al. () were able to detect unequivocal association between temperature and phenotypic traits such as leaf size as well as with genetic loci putatively under divergent selection for temperature. The research by De Kort et al. () provides valuable insight into adaptive response to temperature variation for an ecologically important species and demonstrates the usefulness of an integrated approach for empirical evaluation of local adaptation in nonmodel species (Sork et al. ).
Populations located at the rear-edge of a species' distribution may have disproportionate ecological and evolutionary importance for biodiversity conservation in a changing global environment. Yet genetic studies of such populations remain rare. This study investigates the evolutionary history of North-African low latitude marginal populations of Alnus glutinosa Gaertn., a European tree species that plays a significant ecological role as a keystone of riparian ecosystems. We genotyped 551 adults from 19 populations located across North Africa at 12 microsatellite loci and applied a coalescent-based simulation approach to reconstruct the demographic and evolutionary history of these populations. Surprisingly, Moroccan trees were tetraploids demonstrating a strong distinctiveness of these populations within a species otherwise known as diploid. Best-fitting models of demographic reconstruction revealed the relict nature of Moroccan populations that were found to have withstood past climate change events and to be much older than Algerian and Tunisian populations. This study highlights the complex demographic history that can be encountered in rear-edge distribution margins that here consist of both old stable climate relict and more recent populations, distinctively diverse genetically both quantitatively and qualitatively. We emphasize the high evolutionary and conservation value of marginal rear-edge populations of a keystone riparian species in the context of on-going climate change in the Mediterranean region.
PREMISE OF THE STUDY:Microsatellite markers were developed using second-generation sequencing in Solanum rostratum as a tool to study the reproductive biology and genetic structure of this invasive species.METHODS AND RESULTS:Thirteen microsatellites were successfully discovered and amplified in a single multiplexed PCR. All loci showed genetic variation in S. rostratum. Cross-amplification in five closely related taxa was successful for a subset of loci.CONCLUSIONS:The set of 13 microsatellite markers developed here provides a time-effective and cost-effective genetic tool to study the reproductive biology of S. rostratum. The demonstrated transferability of the PCR multiplex to related taxa also highlights its usefulness for evolutionary studies across Solanum sect. Androceras.
Despite longstanding research, how anthropogenic disturbance affects the genetics of tree populations remains poorly understood. Although empirical evidence often conflicts with theoretical expectations, little progress has been made in refining experimental design or in reformulating theoretical hypotheses. Such progress is, however, essential to understand how forest tree species can tolerate anthropogenic disturbance. Further advances in forest fragmentation genetics research will require that processes driving reproduction and recruitment in fragmented populations are assessed from a tree's perspective instead of experimental convenience, using a multidisciplinary approach to explain the spatiotemporal dynamics of gene dispersal. In this opinion article we aim to inspire a new perspective in forest fragmentation genetics research.
Streamlining the development and genotyping of microsatellites in species for which no genetic information is available represents an important technical challenge to overcome in order to enable mainstream application of state-of-the-art population genetic analysis techniques in nonmodel organisms. Using the example of Acacia harpophylla, an acacia tree endemic of north-eastern Australia, we show that high-throughput shotgun pyrosequencing technology, so-called second-generation sequencing, reduces time and cost of microsatellite marker discovery in nonmodel organisms and of their large-scale typing in natural populations. We found that 0.5% of short sequence reads generated on 454 Genome Sequencer FLX Titanium from random genome sampling and 2.2% of reads generated with prior microsatellite enrichment yielded microsatellite markers with designed polymerase chain reaction (PCR) primers, suggesting that enrichment increases efficiency of pyrosequencing when microsatellite discovery is the primary goal. Using stringent selection criteria to facilitate downstream PCR multiplex design, we identified 1435 microsatellite loci with designed primers from a total of 200,908 short sequence reads. From a subset of 96 loci tested for amplification, 38 were validated for population genetics applications, leading to the optimization of a cost-effective multiplex PCR protocol for the simultaneous typing of nine microsatellites in natural populations of A. harpophylla.
Recent developments in sequencing technologies and bioinformatics analyses provide an unprecedented opportunity for cost and time effective high quality microsatellite marker discovery in nonmodel organisms for which no genomic information is available. Here, we use shotgun pyrosequencing of a microsatellite-enriched library to develop, for the first time, microsatellite markers for Alnus glutinosa, a keystone tree species of European riparian woodland communities. From a total of 17 855 short sequences, we identified 590 perfect microsatellites from which 392 had designed primers. A subset of 48 loci were tested for amplification, 12 of which were polymorphic in A. glutinosa. These 12 loci were successfully coamplified in a single multiplex polymerase chain reaction experiment and validated for population genetics applications. In addition, 10 and 8 of these microsatellites were found to be transferable to the related A. incana and A. cordata species. The developed multiplex of 12 microsatellite markers therefore provides new opportunities for experimental evolutionary and forest genetics research in Alnus.
The spotted gum species complex represents a group of four eucalypt hardwood taxa that have a native range that spans the east coast of Australia, with a morphological cline from Victoria to northern Queensland. Of this group, Corymbia citriodora subsp. variegata (CCV) is widespread in south-eastern Queensland and northern New South Wales. It is currently the most commonly harvested native hardwood in Queensland. However, little basic knowledge of the reproductive biology of the species is available to inform genetic improvement and resource management programmes. Here we take an integrative approach, using both field and molecular data, to identify ecological factors important to mating patterns in native populations of CCV. Field observation of pollinator visitation and flowering phenology of 20 trees showed that foraging behaviour of pollinator guilds varies depending on flowering phenology and canopy structure. A positive effect of tree mean flowering effort was found on insect visitation, while bat visitation was predicted by tree height and by the number of trees simultaneously bearing flowers. Moreover, introduced honeybees were observed frequently, performing 73% of detected flower visits. Conversely, nectar-feeding birds and mammals were observed sporadically with lorikeets and honeyeaters each contributing to 11% of visits. Fruit bats, represented solely by the grey-headed flying fox, performed less than 2% of visits. Genotyping at six microsatellite markers in 301 seeds from 17 families sampled from four of Queensland's native forests showed that CCV displays a mixed-mating system that is mostly outcrossing (t m = 0.899 ± 0.021). Preferential effective pollination from near-neighbours was detected by means of maximum-likelihood paternity analysis with up to 16% of reproduction events resulting from selfing. Forty to 48% of fertilising pollen was also carried from longer distance (>60 m). Marked differences in foraging behaviour and visitation frequency between observed pollinator guilds suggests that the observed dichotomy of effective pollen movement in spotted gums may be due to frequent visit from introduced honeybees favouring geitonogamy and sporadic visits from honeyeaters and fruit bats resulting in potential long-distance pollinations.
Paternity analysis based on microsatellite marker genotyping was used to infer contemporary genetic connectivity by pollen of three population remnants of the wind-pollinated, wind-dispersed tree Fraxinus excelsior, in a deforested Scottish landscape. By deterministically accounting for genotyping error and comparing a range of assignment methods, individual-based paternity assignments were used to derive population-level estimates of gene flow. Pollen immigration into a 300 ha landscape represents between 43 and 68% of effective pollination, mostly depending on assignment method. Individual male reproductive success is unequal, with 31 of 48 trees fertilizing one seed or more, but only three trees fertilizing more than ten seeds. Spatial analysis suggests a fat-tailed pollen dispersal curve with 85% of detected pollination occurring within 100 m, and 15% spreading between 300 and 1900 m from the source. Identification of immigrating pollen sourced from two neighbouring remnants indicates further effective dispersal at 2900 m. Pollen exchange among remnants is driven by population size rather than geographic distance, with larger remnants acting predominantly as pollen donors, and smaller remnants as pollen recipients. Enhanced wind dispersal of pollen in a barren landscape ensures that the seed produced within the catchment includes genetic material from a wide geographic area. However, gene flow estimates based on analysis of non-dispersed seeds were shown to underestimate realized gene immigration into the remnants by a factor of two suggesting that predictive landscape conservation requires integrated estimates of post-recruitment gene flow occurring via both pollen and seed.
The theoretical impacts of anthropogenic habitat degradation on genetic resources have been well articulated. Here we use a simulation approach to assess the magnitude of expected genetic change, and review 31 studies of 23 neotropical tree species to assess whether empirical case studies conform to theory. Major differences in the sensitivity of measures to detect the genetic health of degraded populations were obvious. Most studies employing genetic diversity (nine out of 13) found no significant consequences, yet most that assessed progeny inbreeding (six out of eight), reproductive output (seven out of 10) and fitness (all six) highlighted significant impacts. These observations are in line with theory, where inbreeding is observed immediately following impact, but genetic diversity is lost slowly over subsequent generations, which for trees may take decades. Studies also highlight the ecological, not just genetic, consequences of habitat degradation that can cause reduced seed set and progeny fitness. Unexpectedly, two studies examining pollen flow using paternity analysis highlight an extensive network of gene flow at smaller spatial scales (less than 10 km). Gene flow can thus mitigate against loss of genetic diversity and assist in long-term population viability, even in degraded landscapes. Unfortunately, the surveyed studies were too few and heterogeneous to examine concepts of population size thresholds and genetic resilience in relation to life history. Future suggested research priorities include undertaking integrated studies on a range of species in the same landscapes; better documentation of the extent and duration of impact; and most importantly, combining neutral marker, pollination dynamics, ecological consequences, and progeny fitness assessment within single studies.
In endangered species, it is critical to analyse the level at which populations interact (i.e. dispersal) as well as the levels of inbreeding and local adaptation to set up conservation policies. These parameters were investigated in the endangered species Parnassia palustris living in contrasted habitats. We analysed population structure in 14 populations of northern France for isozymes, cpDNA markers and phenotypic traits related to fitness. Within population genetic diversity and inbreeding coefficients were not correlated to population size. Populations seem not to have undergone severe recent bottleneck. Conversely to pollen migration, seed migration seems limited at a regional scale, which could prevent colonization of new sites even if suitable habitats appear. Finally, the habitat type affects neither within-population genetic diversity nor genetic and phenotypic differentiation among populations. Thus, even if unnoticed local adaptation to habitats exists, it does not influence gene flow between populations.