Species delimitation and geography

Despite the importance of the geographical arrangement of populations for the inference of species boundaries, only a few approaches that integrate spatial information into species delimitation have thus far been developed. Persistent differentiation of sympatric groups of individuals is the best criterion for species status. Species delimitation becomes more prone to error if allopatric metapopulations are considered because it is often difficult to assess whether observed differences between allopatric metapopulations would be sufficient to prevent the fusion of these metapopulations upon contact. We propose a novel approach for testing the hypothesis that the multilocus genetic distances between individuals or populations belonging to two different candidate species are not larger than expected based on their geographical distances and the relationship of genetic and geographical distances within the candidate species. A rejection of this null hypothesis is an argument for classifying the two studied candidate species as distinct species. Case studies show that the proposed tests are suitable to distinguish between intra‐ and interspecific differentiation. The regression approach proposed here is more appropriate for testing species hypotheses with regard to isolation by distance than (partial) Mantel tests. Our tests assume a linear relationship between genetic and (transformed) geographical distances. This assumption can be compromised by a high genetic variability within populations as found in a case study with microsatellite markers.     

Genetic differentiation between individuals

I describe a method of analysis of genetic differentiation which is suitable for the comparison of genetic and demographic estimates of the ‘neighborhood size’– more precisely the product of density and second moment of dispersal distance σ²– in continuous populations sampled at the smallest scale. This method is based on results of models of isolation by distance common to a wide variety of dispersal distances. The performance of this method is tested by simulation for some highly leptokurtic dispersal distributions, and it is applied to a previous study of a kangaroo‐rat (Dipodomys spectabilis) population. In this case, genetic and demographic estimates are within a factor of two from each other. Thus, in line with some previous examples, this study shows that a better agreement may be attained than is usually recognized between genetic and demographic estimates.     

Living on the edge: the role of geography and environment in structuring genetic variation in the southernmost populations of a tropical oak

Understanding the factors determining genetic diversity and structure in peripheral populations is a long‐standing goal of evolutionary biogeography, yet little empirical information is available for tropical species. In this study, we combine information from nuclear microsatellite markers and niche modelling to analyse the factors structuring genetic variation across the southernmost populations of the tropical oak Quercus segoviensis. First, we tested the hypothesis that genetic variability decreases with population isolation and increases with local habitat suitability and stability since the Last Glacial Maximum (LGM). Second, we employed a recently developed multiple matrix regression with randomisation (MMRR) approach to study the factors associated with genetic divergence among the studied populations and test the relative contribution of environmental and geographic isolation to contemporary patterns of genetic differentiation. We found that genetic diversity was negatively correlated with average genetic differentiation with other populations, indicating that isolation and limited gene flow have contributed to erode genetic variability in some populations. Considering the relatively small size of the study area (<120 km), analyses of genetic structure indicate a remarkable inter‐population genetic differentiation. Environmental dissimilarity and differences in current and past climate niche suitability and their additive effects were not associated with genetic differentiation after controlling for geographic distance, indicating that local climate does not contribute to explain spatial patterns of genetic structure. Overall, our data indicate that geographic isolation, but not current or past climate, is the main factor determining contemporary patterns of genetic diversity and structure within the southernmost peripheral populations of this tropical oak.     

Interactions between environmental factors can hide isolation by distance patterns: a case study of Ctenomys rionegrensis in Uruguay

Identifying the factors responsible for the structuring of genetic diversity is of fundamental importance for biodiversity conservation. However, arriving at such understanding is difficult owing to the many factors involved and the potential interactions between them. Here, we present an example of how such interactions can preclude us from arriving at a complete characterization of the demographic history and genetic structure of a species. Ctenomys rionegrensis is a species with restricted dispersal abilities and, as such, should exhibit an isolation by distance (IBD) pattern, which previous studies were unable to uncover. It was therefore concluded that this species underwent a recent population expansion. Using a novel hierarchical Bayesian method, we show that the inability to detect the IBD pattern is due to the interaction between elevation and geographical distance. We posit that populations in low areas suffer periodic floods that may reduce local population sizes, increasing genetic drift, a process that masks the effect of distance on genetic differentiation. Our results do not refute the possibility that the populations of C. rionegrensis underwent a recent population expansion but they indicate that an alternative scenario described by a metapopulation model at or near migration-drift equilibrium cannot be excluded either.     

Limitations to the inference of gene flow at regional geographic scales–an example from the Pieris napi group (Lepidoptera: Pieridae) in Europe

We used hierarchical and pairwise F-statistics to describe genetic differentiation and infer gene flow (M) on local and regional scales within and among parapatric European butterfly taxa in the Pieris napi (L.) group. Within-population allozyme variability is consistently high, and local effective population sizes are inferred to be in the thousands of individuals. The pairwise analysis yields an average neighbourhood area of radius 3.5 km. Among populations within most regions, differentiation is low and M > 2 effective individuals population^-1generation^-1. Pairwise comparisons within the brilannica group show a disjunction indicating that it is out of equilibrium, perhaps as a result of secondary contact between highland and lowland groups. Comparison between meridionalis groups on mainland Italy and Corsica yields M > 12; this is surely loo high and lack of equilibrium resulting from initial colonization is suspected. The hierarchical analysis indicates that 23 ≤0020M≤ 88 among the taxa napi, bryoniae and meridionalis that meet in hybrid zones; no effective gene flow barrier exists among them. This high estimate could also result from recent primary contact, but such a genetic barrier should produce the ‘edge effects' seen in population genetic simulations, and no evidence of this was found among geographically close samples of napi and bryoniae populations from Switzerland. Studies of gene flow among geographic regions are greatly limited by the equilibrium assumption, though studies of local differentiation are much less so. Population studies of gene flow on local scales at regional boundaries provide limited means of testing the equilibrium assumption, and both regional and local analyses provide testable predictions about local population structure. When the equilibrium assumption is not upheld, local patterns at regional boundaries can provide historical information about primary vs. secondary contact.     

Isolation by environment

The interactions between organisms and their environments can shape distributions of spatial genetic variation, resulting in patterns of isolation by environment (IBE) in which genetic and environmental distances are positively correlated, independent of geographic distance. IBE represents one of the most important patterns that results from the ways in which landscape heterogeneity influences gene flow and population connectivity, but it has only recently been examined in studies of ecological and landscape genetics. Nevertheless, the study of IBE presents valuable opportunities to investigate how spatial heterogeneity in ecological processes, agents of selection and environmental variables contributes to genetic divergence in nature. New and increasingly sophisticated studies of IBE in natural systems are poised to make significant contributions to our understanding of the role of ecology in genetic divergence and of modes of differentiation both within and between species. Here, we describe the underlying ecological processes that can generate patterns of IBE, examine its implications for a wide variety of disciplines and outline several areas of future research that can answer pressing questions about the ecological basis of genetic diversity.     

Quantifying the roles of ecology and geography in spatial genetic divergence

Investigating the properties of ecological landscapes that influence gene flow among populations can provide key insights into the earliest stages of biological divergence. Both ecological and geographical factors can reduce gene flow, which can lead to population divergence, but we know little of the relative strengths of these phenomena in nature. Here, we use a novel application of structural equation modelling to quantify the contributions of ecological and geographical isolation to spatial genetic divergence in 17 species of Anolis lizards. Our comparative analysis shows that although both processes contributed significantly, geographical isolation explained substantially more genetic divergence than ecological isolation (36.3 vs. 17.9% of variance respectively), suggesting that despite the proposed ubiquity of ecological divergence, non‐ecological factors play the dominant role in the evolution of spatial genetic divergence.     

Small-Scale Fragmentation Effects on Local Genetic Diversity in Two Phyllostomid Bats with Different Dispersal Abilities in Panama

Habitat fragmentation is one of the greatest threats to biodiversity. Despite their importance for conservation, the genetic consequences of small-scale habitat fragmentation for bat populations are largely unknown. In this study, we linked genetic with ecological and demographic data to assess the effects of habitat fragmentation on two species of phyllostomid bats ( Uroderma bilobatum and Carollia perspicillata ) that differ in their dispersal abilities and demographic response to fragmentation. We hypothesized that population differentiation and the effect of habitat fragmentation on levels of genetic diversity will be a function of the species' mobility. We sequenced mtDNA from 232 bats caught on 11 islands in Gatún Lake, Panamá, isolated from the mainland for ca 90 yr, and in adjacent, continuous forest on the mainland. Populations of both species showed significant genetic differentiation ( F _ST). Consistent with our prediction, population subdivision was lower in the highly mobile U. bilobatum ( F _ST= 0.01) compared to the less vagile C. perspicillata ( F _ST= 0.06), and only the latter species showed a pattern indicative of isolation by distance and, in addition, an effect of fragmentation. Genetic erosion as a result of fragmentation was also only detectable in the less mobile species, C. perspicillata , where haplotype diversity was lower in island compared to mainland populations. Our results suggest that some Neotropical bat species are prone to loss of genetic variation in response to anthropogenic small-scale habitat fragmentation. In this context, our findings point toward mobility as a good predictor of a species' vulnerability to fragmentation and altered population genetic structure.     

Dispersal and gene flow in the rare, parasitic Large Blue butterfly Maculinea arion

Dispersal is crucial for gene flow and often determines the long‐term stability of meta‐populations, particularly in rare species with specialized life cycles. Such species are often foci of conservation efforts because they suffer disproportionally from degradation and fragmentation of their habitat. However, detailed knowledge of effective gene flow through dispersal is often missing, so that conservation strategies have to be based on mark–recapture observations that are suspected to be poor predictors of long‐distance dispersal. These constraints have been especially severe in the study of butterfly populations, where microsatellite markers have been difficult to develop. We used eight microsatellite markers to analyse genetic population structure of the Large Blue butterfly Maculinea arion in Sweden. During recent decades, this species has become an icon of insect conservation after massive decline throughout Europe and extinction in Britain followed by reintroduction of a seed population from the Swedish island of Öland. We find that populations are highly structured genetically, but that gene flow occurs over distances 15 times longer than the maximum distance recorded from mark–recapture studies, which can only be explained by maximum dispersal distances at least twice as large as previously accepted. However, we also find evidence that gaps between sites with suitable habitat exceeding ～20 km induce genetic erosion that can be detected from bottleneck analyses. Although further work is needed, our results suggest that M. arion can maintain fully functional metapopulations when they consist of optimal habitat patches that are no further apart than ～10 km.     

A note on the nutrition of Rhodella

Rhodella maculata Evans has been grown in axenic culture. The alga is euryhaline but will not grow in acid conditions. It requires an exogenous source of vitamin B₁₂, is unable to grow in the dark on acetate or glucose and has a limited capability for utilising organic nitrogen sources.     

Genetic diversity and genetic structure of populations of Ranunculus japonicus Thunb. (Ranunculaceae)

Abstract In order to clarify the genetic diversity and population structure of Ranunculus japonicus , allozymic analysis was conducted on 60 populations in southwestern Japan. Considerable genetic variati ons were detected among the populations of R. japonicus . The genetic diversities within species ( H _es = 0.215) and within populations ( H _ep = 0.172) were slightly higher than those of other perennial herbs with widespread distribution and outcrossing plants. Significantly higher values of fixation index were detected in some populations, which might have arisen from restricted mating partners. The majority of genetic variation (approx. 80%) resided within a population and a moderate level of genetic differentiation ( G _ST = 0.203) was observed among populations. The F _ST value (0.203) suggests the existence of a substantial population structure in this species. The highly significant correlation between geographic distance and F _ST values indicates that isolation by distance has played an important role in the construction of the genetic structure of this species.     

Temporal genetic stability in natural populations of the waterflea Daphnia magna in response to strong selection pressure

Studies monitoring changes in genetic diversity and composition through time allow a unique understanding of evolutionary dynamics and persistence of natural populations. However, such studies are often limited to species with short generation times that can be propagated in the laboratory or few exceptional cases in the wild. Species that produce dormant stages provide powerful models for the reconstruction of evolutionary dynamics in the natural environment. A remaining open question is to what extent dormant egg banks are an unbiased representation of populations and hence of the species’ evolutionary potential, especially in the presence of strong environmental selection. We address this key question using the water flea Daphnia magna, which produces dormant stages that accumulate in biological archives over time. We assess temporal genetic stability in three biological archives, previously used in resurrection ecology studies showing adaptive evolutionary responses to rapid environmental change. We show that neutral genetic diversity does not decline with the age of the population and it is maintained in the presence of strong selection. In addition, by comparing temporal genetic stability in hatched and unhatched populations from the same biological archive, we show that dormant egg banks can be consulted to obtain a reliable measure of genetic diversity over time, at least in the multidecadal time frame studied here. The stability of neutral genetic diversity through time is likely mediated by the buffering effect of the resting egg bank.     

A comparison of individual‐based genetic distance metrics for landscape genetics

A major aim of landscape genetics is to understand how landscapes resist gene flow and thereby influence population genetic structure. An empirical understanding of this process provides a wealth of information that can be used to guide conservation and management of species in fragmented landscapes and also to predict how landscape change may affect population viability. Statistical approaches to infer the true model among competing alternatives are based on the strength of the relationship between pairwise genetic distances and landscape distances among sampled individuals in a population. A variety of methods have been devised to quantify individual genetic distances, but no study has yet compared their relative performance when used for model selection in landscape genetics. In this study, we used population genetic simulations to assess the accuracy of 16 individual‐based genetic distance metrics under varying sample sizes and degree of population genetic structure. We found most metrics performed well when sample size and genetic structure was high. However, it was much more challenging to infer the true model when sample size and genetic structure was low. Under these conditions, we found genetic distance metrics based on principal components analysis were the most accurate (although several other metrics performed similarly), but only when they were derived from multiple principal components axes (the optimal number varied depending on the degree of population genetic structure). Our results provide guidance for which genetic distance metrics maximize model selection accuracy and thereby better inform conservation and management decisions based upon landscape genetic analysis.     

EXAMINING THE FULL EFFECTS OF LANDSCAPE HETEROGENEITY ON SPATIAL GENETIC VARIATION: A MULTIPLE MATRIX REGRESSION APPROACH FOR QUANTIFYING GEOGRAPHIC AND ECOLOGICAL ISOLATION

Understanding the effects of landscape heterogeneity on spatial genetic variation is a primary goal of landscape genetics. Ecological and geographic variables can contribute to genetic structure through geographic isolation, in which geographic barriers and distances restrict gene flow, and ecological isolation, in which gene flow among populations inhabiting different environments is limited by selection against dispersers moving between them. Although methods have been developed to study geographic isolation in detail, ecological isolation has received much less attention, partly because disentangling the effects of these mechanisms is inherently difficult. Here, I describe a novel approach for quantifying the effects of geographic and ecological isolation using multiple matrix regression with randomization. I explored the parameter space over which this method is effective using a series of individual‐based simulations and found that it accurately describes the effects of geographic and ecological isolation over a wide range of conditions. I also applied this method to a set of real‐world datasets to show that ecological isolation is an often overlooked but important contributor to patterns of spatial genetic variation and to demonstrate how this analysis can provide new insights into how landscapes contribute to the evolution of genetic variation in nature.     

Gene flow between nascent species: geographic,genotypic and phenotypic differentiation within and between Aquilegia formosa and A. pubescens

Speciation can be described as a reduction, and the eventual cessation, in the ability to interbreed. Thus, determining how gene flow differs within and between nascent species can illuminate the relative stage the taxa have attained in the speciation process. Aquilegia formosa and A. pubescens are fully intercompatible, yet occur in different habitats and have flowers specialized for pollination by hummingbirds and hawkmoths, respectively. Using 79 SNP loci, we genotyped nearly 1000 individuals from populations of both species in close proximity to each other and from putative hybrid zones. The species shared all but one SNP polymorphism, and on average, allele frequencies differed by only 0.14. However, the species were clearly differentiated using Structure, and admixed individuals were primarily identified at putative hybrid zones. PopGraph identified a highly integrated network among all populations, but populations of each species and hybrid zones occupied distinct regions in the network. Using either conditional graph distance (cGD) or Fst/(1‐Fst), we found significant isolation by distance (IBD) among populations. Within species, IBD was strong, indicating high historic gene flow. IBD extended approximately 100 km in A. pubescens and 30 km in A. formosa. However, IBD between the species was very weak and extended only a few km beyond hybrid zones, suggesting little recent gene flow. The extensive sharing of SNP polymorphisms between these species suggests that they are very early in the speciation process while the low signal of IBD suggests that they have largely ceased gene exchange.     

Microsatellite analysis of genetic variation among and within Alpine marmot populations in the French Alps

The genetic structure of the Alpine marmot, Marmota marmota, was studied by an analysis of five polymorphic microsatellite loci. Eight locations were sampled in the French Alps, one from Les Ecrins valley (n = 160), another from La Sassière valley (n = 289) and the six others from the Maurienne valley (n = 139). Information on social group structure was available for both Les Ecrins and La Sassière but not for the other samples. The high levels of genetic diversity observed are at odds with the results obtained using microsatellites, minisatellites and allozymes on Alpine marmots from Germany, Austria and Switzerland. Strong deficits in heterozygotes were found in Les Ecrins and La Sassière. They are caused by a Wahlund effect due to the family structure (i.e. differentiation between the family groups). The family groups exhibit excess of heterozygotes rather than deficits. This may be caused by outbreeding and this is compatible with recent results from the genetics of related social species when information on the social structure is taken into account. The observed outbreeding could be the result of females mating with transient males or males coming from neighbouring colonies. Both indicate that the species may not be as monogamous as is usually believed. The results are also compatible with a male-biased dispersal but do not allow us to exclude some female migration. We also found a significant correlation between geographical and genetic distance indicating that isolation by distance could be an issue in marmots. This study is the first that analysed populations of marmots taking into account the social structure within populations and assessing inbreeding at different levels (region, valley, population, and family groups). Our study clearly demonstrated that the sampling strategy and behavioural information can have dramatic effects on both the results and interpretation of the genetic data.     

The names of Spain: a study of the isonymy structure of Spain

In order to estimate the isonymy structure of Spain, we studied surname distribution in 283 Spanish towns based on 3.625 million telephone users selected from 6.328 million users, downloaded from a commercial CD-ROM which contains all 13 million users in the country. Since in Spain the surname is made by the paternal and the maternal surname, it was possible to classify surnames according to parental origin. Two matrices of isonymy distances, one for paternal and one for maternal surnames, were constructed and tested for correlation with geographic distance. For the whole of Spain, Euclidean distance was significantly but weakly correlated with geographic distance both for paternal and maternal surnames, with r = 0.205 +/- 0.013 and r = 0.263 +/- 0.012, respectively. Two dendrograms of the 283 sampled towns were built from the two matrices of Euclidean distance. They are largely colinear. Four main clusters identified by the dendrograms are correlated with geography. Given the surname structure of Spain, we were able to calculate from isonymy and for each town 1). total or expressed inbreeding, 2). random or expected inbreeding, and 3). local inbreeding. Total inbreeding, F(IT), was highest in the North Atlantic regions and lowest along the Mediterranean Coast. The lowest levels were found in Andalusia, Catalunyia, Valencia, and Navarra. Random inbreeding, F(ST), had a similar geographical pattern. Local inbreeding, F(IS), was relatively uniform in the whole of Spain. In towns, random inbreeding dominates over local inbreeding. From the analysis, it emerges that the northwestern area of Spain is the most inbred.     

History,geography and host use shape genomewide patterns of genetic variation in the redheaded pine sawfly (Neodiprion lecontei)

Divergent host use has long been suspected to drive population differentiation and speciation in plant‐feeding insects. Evaluating the contribution of divergent host use to genetic differentiation can be difficult, however, as dispersal limitation and population structure may also influence patterns of genetic variation. In this study, we use double‐digest restriction‐associated DNA (ddRAD) sequencing to test the hypothesis that divergent host use contributes to genetic differentiation among populations of the redheaded pine sawfly (Neodiprion lecontei), a widespread pest that uses multiple Pinus hosts throughout its range in eastern North America. Because this species has a broad range and specializes on host plants known to have migrated extensively during the Pleistocene, we first assess overall genetic structure using model‐based and model‐free clustering methods and identify three geographically distinct genetic clusters. Next, using a composite‐likelihood approach based on the site frequency spectrum and a novel strategy for maximizing the utility of linked RAD markers, we infer the population topology and date divergence to the Pleistocene. Based on existing knowledge of Pinus refugia, estimated demographic parameters and patterns of diversity among sawfly populations, we propose a Pleistocene divergence scenario for N. lecontei. Finally, using Mantel and partial Mantel tests, we identify a significant relationship between genetic distance and geography in all clusters, and between genetic distance and host use in two of three clusters. Overall, our results indicate that Pleistocene isolation, dispersal limitation and ecological divergence all contribute to genomewide differentiation in this species and support the hypothesis that host use is a common driver of population divergence in host‐specialized insects.