The limits of mtDNA phylogeography: complex patterns of population history in a highly structured Iberian lizard are only revealed by the use of nuclear markers

Phylogeographic analyses based on the sole use of the mitochondrial DNA (mtDNA) molecule reveal only a small part of the evolutionary history of a species or a set of related species. In this study, we have combined the application of slow- and fast-evolving nuclear markers (proteins and microsatellites, respectively) together with the analysis of two-gene genealogies to further understand the history of the Iberian endemic Schreiber's green lizard, Lacerta schreiberi, a species for which a well established phylogeographical scenario is available. In sharp contrast with the observation of four divergent and almost allopatric mtDNA clades, our nuclear data revealed how two groups of populations diverged, persisted and began to admix along the mountains of the Iberian Central System. In addition, the combination of mtDNA and nuclear data showed how the core area of the species distribution responded to ice ages, first by relatively old processes of population expansion to the south followed by episodes of contraction that are at the origin of present-day isolates, and more recently by a postglacial expansion to the Iberian Northwest where new habitats were made available after climatic amelioration. Taken together with recently published results for a variety of other organisms, our results suggest that complex processes of fragmentation, expansion and admixture can only be properly addressed through the use of several and complementary types of molecular markers. Finally, we also suggest that southern European refugia are both hotspots and melting pots of genetic diversity.     

Miocene to Pliocene vegetation reconstruction and climate estimates in the Iberian Peninsula from pollen data

Pollen analysis of Miocene and Pliocene sediments from the Iberian Peninsula shows a progressive reduction in plant diversity through time caused by the disappearance of thermophilous and high-water requirement plants. In addition, an increase in warm-temperate (mesothermic), seasonal-adapted “Mediterranean” taxa, high-elevation conifers and herbs (mainly Artemisia) occurred during the Middle and Late Miocene and Pliocene. This has mainly been interpreted as a response of the vegetation to global and regional processes, including climate cooling related to the development of the East Antarctic Ice Sheet and then the onset of the Arctic Ice Sheet, uplift of regional mountains related to the Alpine uplift and the progressive movement of Eurasia towards northern latitudes as a result of the northwards subduction of Africa. The development of steppe-like vegetation in southern Iberia is ancient and probably started during the Oligocene. The onset of a contrasted seasonality in temperature during the Mid-Pliocene superimposed on the pre-existing seasonality in precipitation, the annual length of which increased southward. The Mediterranean climatic rhythm (summer drought) began about 3.4 Ma and caused the individualization of modern Mediterranean ecosystems. Quaternary-type Mediterranean climatic fluctuations started at 2.6 Ma (Gelasian) resulting in repeated steppe vs. forest alternations. A latitudinal climatic gradient between the southern and the northern parts of the Iberian Peninsula existed since the Middle Miocene.     

Population genetic structure in the Iberian cyprinid fish Iberochondrostoma lemmingii (Steindachner, 1866): disentangling species fragmentation and colonization processes

Interplay between the complex geography, hydrogeomorphological history, past climatic changes, and anthropogenic pressures is likely responsible for the current diversity and species' distribution of freshwater fishes in the Iberian Peninsula. To further disentangle the evolutionary processes promoting the diversification of endemic Iberian Cyprinids through time and space, we explored the patterns of genetic diversity of the Iberian arched‐mouth nase, Iberochondrostoma lemmingii (Steindachner, 1866), using molecular markers rendering at different timescales: the mitochondrial gene cytochrome b and seven microsatellite loci. Both markers showed significant differentiation of populations though the relative genetic distances among populations were different between markers. Mitochondrial DNA results indicate the isolation of hydrographic basins as the main driver of population differentiation, with Tejo as the centre of diversification. The results also support connections between Tejo, Guadiana, and Guadalquivir, with levels of divergence suggesting an earlier severance of Guadalquivir, whereas Guadiana and Tejo maintained connections until a more recent past. Establishment of more peripherial populations in small southern basins (Quarteira and Almargem) could have been ruled by founder events. However, the analysis of present‐day genetic configuration suggested by microsatellite data implies, for the first time, the involvement of other factors in the evolution of arched‐mouth Iberian nase populations. Relative low genetic distances between inter‐basin populations (Tejo and Guadiana) and the lack of concordance between differentiation and geography suggest a possible influence of human‐mediated translocations in the population genetic patterns of I. lemmingii. High intra‐basin differentiation levels were found within Tejo and Guadiana and may be associated with factors intrinsic to the species (e.g. low dispersal capability) or natural and/or artificial barriers to gene flow. The low vagility of the species appears to be an important factor influencing the evolutionary processes shaping the phylogeographical patterns of I. lemmingii, which could be relevant for the conservation of this threatened species. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 559–572.     

Steppes, savannahs, forests and phytodiversity reservoirs during the Pleistocene in the Iberian Peninsula

A palaeobotanical analysis of the Pleistocene floras and vegetation in the Iberian Peninsula shows the existence of patched landscapes with Pinus woodlands, deciduous and mixed forests, parklands (savannah-like), shrublands, steppes and grasslands. Extinctions of Arctotertiary woody taxa are recorded during the Early and Middle Pleistocene, but glacial refugia facilitated the survival of a number of temperate, Mediterranean and Ibero-North African woody angiosperms. The responses of Iberian vegetation to climatic changes during the Pleistocene have been spatially and temporarily complex, including rapid changes of vegetation in parallel to orbital and suborbital variability, and situations of multi-centennial resilience or accommodation to climatic changes. Regional characteristics emerged as soon as for the Middle Pleistocene, if not earlier: Ericaceae in the Atlantic coast indicating wetter climate, thermo-mediterranean elements in the south as currently, and broad-leaf trees in the northeastern. Overall, steppe landscapes and open Pinus woodlands prevailed over many continental regions during the cold spells of the Late Pleistocene. The maintenance of a high phytodiversity during the glacials was linked to several refuge zones in the coastal shelves of the Mediterranean and intramountainous valleys. Northern Iberia, especially on coastal areas, was also patched with populations of tree species, and this is not only documented by palaeobotanical data (pollen and charcoal) but also postulated by phylogeographical models.     

Assessing the Conservation Status of an Iberian Moth Using Pseudo-Absences

ABSTRACT Knowing the distribution of endangered species is of substantial importance for conservation. We considered a useful approach for modeling species distribution when managing information from atlases and museums but when absence data is not available. By modeling the distribution for Graellsia isabelae, a threatened moth species, we assessed its current conservation status and identified its most relevant distribution explanatory variables using Geographic Information System and Generalized Linear Models. The distribution model was built from 136 occurrence records and 25 digitized explanatory variables at a 10 × 10 km resolution. Model predictions from logistic-regressed pseudo-absences, obtained from a presence-only method (Ecological-Niche Factor Analysis), explained 96.23% of the total deviance. We found that the best predictor variables were summer precipitation, aridity, and mean elevation. With respect to host plants, the presence of G. isabelae associated mainly with Scots pine (Pinus sylvestris) and Austrian pine (P. nigra). The finding of 8 areas, exclusively in the eastern Iberian territory, and a larger unoccupied habitat in the western Iberian Peninsula indicates that this species is probably not in equilibrium with its environment by historical factors. Sites of Community Importance under protection do not seem sufficient to maintain current populations, necessitating the protection of suitable neighboring habitats. Our methodology is useful to manage the conservation status of species for which reliable absence data is not available. It is possible to determine those variables that most affect the distribution of species as well as the potential suitable areas with the purpose of evaluating protected areas, connectivity among populations, and possible reintroductions.     

Genetic signature of recent glaciation on populations of a near-shore marine fish species (Syngnathus leptorhynchus)

Continental glaciation has played a major role in shaping the present-day phylogeography of freshwater and terrestrial species in the Northern Hemisphere. Recent work suggests that coastal glaciation during ice ages may have also had a significant impact on marine species. The bay pipefish, Syngnathus leptorhynchus , is a near-shore Pacific coast fish species with an exceptionally wide latitudinal distribution, ranging from Bahia Santa Maria, Baja California to Prince William Sound, Alaska. Survey data indicate that S. leptorhynchus is experiencing a range expansion at the northern limit of its range, consistent with colonization from southern populations. The present study uses six novel microsatellite markers and mitochondrial DNA (mtDNA) sequence data to study the present-day population genetic structure of four coastal populations of S. leptorhynchus . Deficits in mtDNA and nuclear DNA diversity in northern populations from regions glaciated during the last glacial maximum (LGM) [ c . 18 000 years before present ( bp )] suggest that these populations were effected by glacial events. Direct estimates of population divergence times derived from both isolation and isolation-with-migration models of evolution are also consistent with a postglacial phylogenetic history of populations north of the LGM. Sequence data further indicate that a population at the southern end of the species range has been separated from the three northern populations since long before the last interglacial event ( c . 130 000 years bp ), suggesting that topographical features along the Pacific coast may maintain population separation in regions unimpacted by coastal glaciation.     

Postglacial dispersal of the European rabbit (Oryctolagus cuniculus) on the Iberian peninsula reconstructed from nested clade and mismatch analyses of mitochondrial DNA genetic variation

Nested clade analysis was applied to cytochrome b restriction site data previously obtained on 20 natural populations of the European rabbit across the Iberian Peninsula to test the hypothesis of postglacial dispersal from two main refugia, one in the northeast and the other in the southwest. Apart from historical fragmentation that resulted in geographic discontinuity of two distinct mitochondrial DNA (mtDNA) clades A and B, patterns of haplotype genetic variability have been shaped mostly by restricted gene flow via isolation by distance. The distribution of tip versus interior haplotypes suggests that dispersal occurred from both the southwestern and northeastern groups. Dispersal from the southwest had a north and northwest direction, whereas from the northeast it had mostly a western and southern orientation, with subsequent overlap in a southeastern-northwestern axis across the Iberian Peninsula. The analysis of the pairwise mismatch distribution of a 179-181-bp fragment of the mtDNA control region, for seven of those populations, further supports the idea that major patterns of dispersal were in the direction of central Iberia. Additionally, rabbit populations do not show signs of any significant loss of genetic diversity in the recent past, implying that they maintained large population sizes and structure throughout the ice ages. This is congruent with the fact that the Iberian Peninsula was itself a glacial refugium during Quaternary ice ages. Nonetheless, climatic oscillations of this period, although certainly milder than in northern Europe, were sufficient to affect the range distributions of Iberian organisms.     

Cranial variation in the Iberian Peninsula and the Balearic Islands: Inferences about the history of the population

A multivariate analysis of four prehistoric and nine historic populations from the Iberian Peninsula and Balearic Islands with large sample sizes (n > 30 individuals for the neurocranium and n > 15 for the facial skeleton) is presented, considering 874 male and 557 female skulls and using 20 craniometric measurements. Cluster analyses have been undertaken using the squared Euclidean distance as a measure of proximity and the average linkage between groups (UPGMA), and neighbor-joining algorithms as a branching method, and a bootstrap analysis was used to assess the robustness of the clustering topology. The study was complemented with a principal coordinate analysis and with the application of the Mantel test to measure the degree of correspondence between the information furnished by the female and the male samples. The analyses show that the main source of morphometric variability in the Iberian Peninsula is the Basque population. The second source of variation is provided by two populations (Muslims and Jews), different from the rest from an archaeological and cultural point of view, and can probably be attributed to influences from sub-Saharan Africa. The massive deportations of the Jews in 1492 and of the Moors between the 15th and 17th centuries may have erased this source of variability from the present population of the Iberian Peninsula. The remaining studied populations, including samples from Castile, Cantabria, Andalusia, Catalonia and Balearic Islands, are grouped together, showing a notable morphological homogeneity, despite their temporal and geographic heterogeneity. These results are in general agreement with those obtained in synthetic maps, by analyzing multiple genetic markers. In such studies, the Basque population is described as the main source of genetic variability, not only in the Iberian Peninsula, but also in Western Europe. © 1996 Wiley-Liss, Inc.     

Evolution of rattlesnakes (Viperidae; Crotalus) in the warm deserts of western North America shaped by Neogene vicariance and Quaternary climate change

During Pleistocene, the Laurentide ice sheet rearranged and diversified biotic distributions in eastern North America, yet had minimal physical impact in western North America where lineage diversification is instead hypothesized to result from climatic changes. If Pleistocene climatic fluctuations impacted desert species, the latter would reflect patterns of restricted gene flow concomitant with indications of demographic bottlenecks. Accordingly, molecular evidence for refugia should be present within these distributions and for subsequent range expansions as conditions improved. We sought answers to these questions by evaluating mitochondrial DNA (mtDNA) sequences from four species of rattlesnakes [Crotalus mitchellii (speckled rattlesnake), Crotalus cerastes (sidewinder), Crotalus tigris (tiger rattlesnake), Crotalus ruber (red diamond rattlesnake)] with distributions restricted to desert regions of southwestern North America. We inferred relationships using parsimony and maximum likelihood, tested intraspecific clades for population expansions, applied an isolation-with-migration model to determine bi-directional migration rates (m) among regions, and inferred divergence times for species and clades by applying a semiparametric penalized likelihood approach to our molecular data. Evidence for significant range expansion was present in two of eight regions in two species (Crotalus mitchellii pyrrhus, C. tigris region north). Two species (C. cerastes, C. mitchellii) showed a distribution concomitant with northward displacement of Baja California from mainland México, followed by vicariant separation into subclades. Effects of Pleistocene climate fluctuations were found in the distributions of all four species. Three regional diversification patterns were identified: (i) shallow genetic diversity that resulted from Pleistocene climatic events (C. tigris, C. ruber); (ii) deep Pleistocene divisions indicating allopatric segregation of subclades within refugia (C. mitchellii, C. cerastes); and (iii) lineage diversifications that extended to Pliocene or Late Miocene (C. mitchellii, C. cerastes). Clade-diversifying and clade-constraining effects impacted the four species of rattlesnakes unequally. We found relatively high levels of molecular diversification in the two most broadly distributed species (C. mitchellii, C. cerastes), and lower levels of genetic diversification in the two species (C. tigris, C. ruber) whose ranges are relatively more restricted. Furthermore, in several cases, the distributions of subspecies were not congruent with our molecular information. We suggest regional conservation perspectives for southwestern deserts cannot rely upon subspecies as biodiversity surrogates, but must instead employ a molecular and deep historical perspective as a primary mechanism to frame biodiversity reserves within this region.     

The geography of speciation in narrow‐range endemics of the ‘Haenydra’ lineage (Coleoptera,Hydraenidae, Hydraena)

Aim We test whether species of western Mediterranean aquatic Coleoptera of the ‘Haenydra’ lineage (Hydraenidae, Hydraena) originated through: (1) successive periods of dispersal and speciation, (2) range fragmentation by random vicariance, or (3) range fragmentation by geographic isolation owing to a general reduction of population density. Location Europe. Methods To discriminate between scenarios we use contrasting predictions of the relationship between phylogenetic and geographic distance. The phylogeny was based on 3 kb of four mitochondrial and two nuclear gene fragments of about half of the known species of ‘Haenydra’, including most western Mediterranean taxa. Divergences were estimated using a molecular clock. The relationship between phylogenetic and geographic distance was tested using bivariate plots, Mantel tests and comparison of the observed phylogeny with the one minimizing geographic distances between species, as measured using Euclidean minimum spanning trees (EMSTs). Results The monophyly of ‘Haenydra’ was strongly supported, although its phylogenetic placement was not resolved. ‘Haenydra’ was estimated to be of late Miocene age, with most species originating during the Pleistocene. In two clades (Hydraena tatii and Hydraena emarginata clades) there was a significant association between geographic and phylogenetic distance, and the reconstructed phylogeny was identical to that obtained through the EMST, demonstrating a strong non‐randomness of the geographic distribution of the species. In two other clades (Hydraena iberica and Hydraena bitruncata clades) there was no association between geographic and phylogenetic distance, and the observed phylogeny was not the one minimizing geographic distances. In one of the clades this seems to be due to a secondary, recent range expansion of one species (H. iberica), which erased the geographic signal of their distributions. Main conclusions We show that it is possible to obtain strong evidence of stasis of the geographic ranges of narrow‐range endemic species through the study of their phylogenetic relationships and current distributions. In at least two of the studied clades, current species seem to have originated through the fragmentation of a more widely distributed species, without further range movements. A process of range expansion and fragmentation may have occurred repeatedly within the ‘Haenydra’ lineage, contributing to the accumulation of narrow‐range endemics in Mediterranean Pleistocene refugia.     

Phylogenetic analysis of the Lacerta agilis subspecies complex

The Sand lizard Lacerta agilis inhabits a huge area across Eurasia with several subspecies. Nine subspecies are generally approved, L. a. agilis, L. a. argus, L. a. bosnica; L. a. chersonensis, L. a. exigua, L. a. grusinica, L. a. iorinensis, L. a. brevicaudata and L. a. boemica, but several more have been described. The emergence of this large number of subspecies is connected with the phylogeographic history of this species, defined by geographic and climatic processes. A study revealing phylogenetic relationships of this species was already conducted with a broad taxon sampling and coverage. However, the latter study was based solely on the cytochrome b gene and, furthermore, the Balkan Peninsula inhabited by the subspecies L. a. bosnica was underrepresented. This peninsula is a centre of European herpetofaunal endemism with high levels of phenotypic and genotypic variation. Therefore, the inclusion of the subspecies L. a. bosnica is important to clarify the overall view of the phylogenetic relations within the Lacerta agilis subspecies complex and to investigate the level of population differentiation within this highly diverse area. Thus, the aim of the present study was a more thorough analysis of the Balkan Peninsula with a broader taxon sampling. Furthermore, we extended the available datasets by adding the mitochondrial d-loop region and by further samples of different areas of the distribution range apart from the Balkan Peninsula. Our study reveals that the Balkan Peninsula is apparently inhabited by several differentiated lineages, whereby the Central Greek lineage might deserve subspecies status. Furthermore, the distribution area of the two subspecies L. a. agilis and L. a. argus should be revised, as the results of our study raise doubt about the commonly assumed distribution of both subspecies. As the most important outcome our results support that L. a. boemica deserves species status.     

The automation and evaluation of nested clade phylogeographic analysis

Nested clade phylogeographic analysis (NCPA) is a popular method for reconstructing the demographic history of spatially distributed populations from genetic data. Although some parts of the analysis are automated, there is no unique and widely followed algorithm for doing this in its entirety, beginning with the data, and ending with the inferences drawn from the data. This article describes a method that automates NCPA, thereby providing a framework for replicating analyses in an objective way. To do so, a number of decisions need to be made so that the automated implementation is representative of previous analyses. We review how the NCPA procedure has evolved since its inception and conclude that there is scope for some variability in the manual application of NCPA. We apply the automated software to three published datasets previously analyzed manually and replicate many details of the manual analyses, suggesting that the current algorithm is representative of how a typical user will perform NCPA. We simulate a large number of replicate datasets for geographically distributed, but entirely random-mating, populations. These are then analyzed using the automated NCPA algorithm. Results indicate that NCPA tends to give a high frequency of false positives. In our simulations we observe that 14% of the clades give a conclusive inference that a demographic event has occurred, and that 75% of the datasets have at least one clade that gives such an inference. This is mainly due to the generation of multiple statistics per clade, of which only one is required to be significant to apply the inference key. We survey the inferences that have been made in recent publications and show that the most commonly inferred processes (restricted gene flow with isolation by distance and contiguous range expansion) are those that are commonly inferred in our simulations. However, published datasets typically yield a richer set of inferences with NCPA than obtained in our random-mating simulations, and further testing of NCPA with models of structured populations is necessary to examine its accuracy.     

Climatic niche and neutral genetic diversity of the six Iberian pine species: a retrospective and prospective view

Quaternary climatic fluctuations have left contrasting historical footprints on the neutral genetic diversity patterns of existing populations of different tree species. We should expect the demography, and consequently the neutral genetic structure, of taxa less tolerant to particular climatic extremes to be more sensitive to long‐term climate fluctuations. We explore this hypothesis here by sampling all six pine species found in the Iberian Peninsula (2464 individuals, 105 populations), using a common set of chloroplast microsatellite markers, and by looking at the association between neutral genetic diversity and species‐specific climatic requirements. We found large variation in neutral genetic diversity and structure among Iberian pines, with cold‐enduring mountain species (Pinus uncinata, P. sylvestris and P. nigra) showing substantially greater diversity than thermophilous taxa (P. pinea and P. halepensis). Within species, we observed a significant positive correlation between population genetic diversity and summer precipitation for some of the mountain pines. The observed pattern is consistent with the hypotheses that: (i) more thermophilous species have been subjected to stronger demographic fluctuations in the past, as a consequence of their maladaptation to recurrent glacial cold stages; and (ii) altitudinal migrations have allowed the maintenance of large effective population sizes and genetic variation in cold‐tolerant species, especially in more humid regions. In the light of these results and hypotheses, we discuss some potential genetic consequences of impending climate change.     

Persistent genetic signatures of colonization in Brachionus manjavacas rotifers in the Iberian Peninsula

Recent phylogeographical assessments have consistently shown that continental zooplankton display high levels of population subdivision, despite the high dispersal capacity of their diapausing propagules. As such, there is an apparent paradox between observed cosmopolitanism in the zooplankton that is associated with long-distance dispersal, and strong phylogeographical structures at a regional scale. Such population dynamics, far from migration-drift equilibrium, have been shown in the rotifer species complex Brachionus plicatilis, a group of over a dozen species inhabiting salt lakes and coastal lagoons worldwide. Here we present the mitochondrial DNA phylogeography of one of these species, Brachionus manjavacas, in the Iberian Peninsula, where it often co-occurs with the morphologically similar species B. plicatilis sensu stricto. We obtained sequences from 233 individuals from diapausing eggs and clonal cultures from 16 lakes in the Iberian Peninsula, and a Tunisian lake. Two strongly supported deep mitochondrial DNA clades were found (A and B). Phylogenetic and nested clade analysis showed that clade A has a strong phylogeographical structure, with a strong similarity of phylogeographical patterns between B. manjavacas clade A and B. plicatilis s.s. These include (i) signatures of allopatric fragmentation between central and southern populations, and (ii) range expansions in the Iberian Peninsula, both likely to have occurred during the Pleistocene. We find evidence for a glacial refugium in the Guadiana basin. Clades A and B co-occurred in several of these lakes because of range expansion and secondary contact between both clades. The co-occurrence between B. plicatilis s.s. and B. manjavacas is not recent, and both species might have experienced similar environmental challenges during the Pleistocene. The strong correlation of genetic and geographical distance found suggests that historical events can lead to such correlation, mirroring the effects of 'isolation by distance' in equilibrium populations.     

The rise and fall of the mountain hare (Lepus timidus) during Pleistocene glaciations: expansion and retreat with hybridization in the Iberian Peninsula

The climatic fluctuations during glaciations have affected differently arctic and temperate species. In the northern hemisphere, cooling periods induced the expansion of many arctic species to the south, while temperate species were forced to retract in southern refugia. Consequently, in some areas the alternation of these species set the conditions for competition and eventually hybridization. Hares in the Iberian Peninsula appear to illustrate this phenomenon. Populations of Iberian hare (Lepus granatensis), brown hare (Lepus europaeus) and broom hare (Lepus castroviejoi) in Northern Iberia harbour mitochondrial haplotypes from the mountain hare (Lepus timidus), a mainly boreal and arctic species presently absent from the peninsula. To understand the history of this past introgression we analysed sequence variation and geographical distribution of mitochondrial control region and cytochrome b haplotypes of L. timidus origin found in 378 specimens of these four species. Among 124 L. timidus from the Northern Palaearctic and the Alps we found substantial nucleotide diversity (2.3%) but little differentiation between populations. Based on the mismatch distribution of the L. timidus sequences, this could result from an expansion at a time of temperature decrease favourable to this arctic species. The nucleotide diversity of L. timidus mtDNA found in Iberian L. granatensis, L. europaeus and L. castroviejoi (183, 70 and 1 specimens, respectively) was of the same order as that in L. timidus over its range (1.9%), suggesting repeated introgression of multiple lineages. The structure of the coalescent of L. granatensis sequences indicates that hybridization with L. timidus was followed by expansion of the introgressed haplotypes, as expected during a replacement with competition, and occurred when temperatures started to rise, favouring the temperate species. Whether a similar scenario explains the introgression into Iberian L. europaeus remains unclear but it is possible that it hybridized with already introgressed L. granatensis.     

Gene flow, dispersal, and nested clade analysis among populations of the stonefly Peltoperla tarteri in the southern Appalachians

We examined the genetic structure and phylogeography of populations of the stonefly Peltoperla tarteri in the Southern Appalachians to determine the extent and likely mechanism for dispersal of this stream insect. A 454-base-pair (bp) portion of the mitochondrial control region was sequenced from a minimum of 20 individuals from eight populations. Pairwise FST and exact tests showed high levels of differentiation among almost all populations except those on the same stream. amova analysis detected significant genetic differentiation between streams within drainages (phi(SD) = 0.14, P < 0.001), and there was a slight positive correlation between aquatic distance and genetic distance (r = 0.295, P = 0.03). According to nested clade analysis, the present day pattern of genetic variation in P. tarteri is the result of a historical range expansion coupled with restricted gene flow with isolation by distance. Together, these analyses suggest that adult dispersal is limited and that movement by larvae is the primary dispersal mechanism for P. tarteri.     

Post-glacial colonization of Western Europe brown bears from a cryptic Atlantic refugium out of the Iberian Peninsula

The European brown bear (Ursus arctos) shows a particular phylogeography that has been used to illustrate the model for contraction-expansion dynamics related to glacial refugia in Southern European peninsulas. Recent studies, however, have nuanced the once generally accepted paradigm, indicating the existence of cryptic refugia for some species further north. In this paper we collected available data on chronology and mitochondrial haplotypes from Western European brown bears, adding new sequences from present day individuals from the Cantabrian (North Iberia) area, in order to reconstruct the dynamics of the species in the region. Both genetics and chronology show that the Iberian Pleistocene lineages were not the direct ancestors of the Holocene ones, the latter entering the Peninsula belatedly (around 10,000 years BP) with respect to other areas such as the British Isles. We therefore propose the existence of a cryptic refugium in continental Atlantic Europe, from where the bears would expand as the ice receded. The delay in the recolonization of the Iberian Peninsula could be due to the orographic characteristics of the Pyrenean-Cantabrian region and to the abundant presence of humans in the natural entrance to the Peninsula.     

Mitochondrial DNA phylogeography of Lissotriton boscai (Caudata, Salamandridae): evidence for old, multiple refugia in an Iberian endemic

In Europe, southern peninsulas served as refugia during cold periods in the Pleistocene, acting both as centres of origin of endemisms and as sources from which formerly glaciated areas were recolonized during interglacial periods. Previous studies have revealed that within the main refugial areas, intraspecific lineages often survived in allopatric refugia. We analysed two mitochondrial markers (nad4, control region, approximately 1.4 kb) in 103 individuals representing the entire distribution of Lissotriton boscai, a newt endemic to the western Iberian Peninsula. We inferred the evolutionary history of the species through phylogenetic, phylogeographic and historical demographic analyses. The results revealed unexpected, deep levels of geographically structured genetic variability. We identified two main evolutionary lineages, each containing three well-supported clades. The first historical split involved populations from central-southwestern coastal Portugal and the ancestor of all the remaining populations around 5.8 million years ago. Both lineages were subsequently fragmented into different population groups between 2.5 and 1.2 million years ago. According to nested clade analysis, at lower hierarchical levels the patterns suggest restricted gene flow with isolation by distance, whereas at higher levels the clades exhibit signatures of contiguous range expansion. Bayesian Skyline Plots show recent bottlenecks, followed by demographic expansions in all lineages. The significant genetic structure found is consistent with long-term survival of populations in allopatric refugia, supporting the 'refugia-within-refugia' scenario for southern European peninsulas. The comparison of our results with other co-distributed species highlights the generality of this hypothesis for the Iberian herpetofauna and suggests that Mediterranean refuges had more relevance for the composition and distribution of present biodiversity patterns than currently acknowledged. We briefly discuss the taxonomic and conservation implications of our results.     

Monitoring of marine benthic communities and taxonomic resolution: an approach through diverse habitats and substrates along the Southern Iberian coastline

Studies conducted along the southern Iberian coastline validate macrobenthic community analyses at taxonomic levels higher than that of species. Twelve studies on littoral benthic communities, carried out by the same research team, were selected spanning both a variety of sampling strategies (spatial, temporal, spatio-temporal) and substrate/habitat types (sediment, rock, algae). In order to establish differences between the results obtained at the taxonomic levels of species, family and order, similarities among stations were calculated using Spearman’s coefficient for ranges. A subset of three studies was selected to investigate possible differences in ‘best-explaining’ environmental variables with taxonomic level. The environmental variables selected at species level were the same as those found at levels of family and order. It is concluded that studies at the different levels of taxonomic resolution (species, family, order) lead to similar results both with regard to relative community distributions and the environmental variables associated with these. The importance of this result for monitoring similar benthic communities is discussed.     

Vegetation history and human activity during the last 6000 years on the central Catalan coast (northeastern Iberian Peninsula)

In this paper, we present the results of four pollen diagrams obtained from the northeastern coast of the Iberian Peninsula. These data, together with a set of 11 radiocarbon dates, allow us to make some suggestions about human activity in this area mainly during the last 6000 years. We have established four main stages of this activity. Phase I (7000–3000 B.P.) shows some sporadic human clearance without qualitative and/or quantitative changes inside the natural woodland, except for slight oscillations in AP values. These clearances are synchronous with Neolithic and Bronze Age settlements. Phase II (3000–1500/1300 B.P.) demonstrates a different human action on the landscape along the Catalan coast. Sampling sites located in rich agricultural plains and close to urban centres show continuous woodland clearance during the Iberian period and especially during Roman times. In contrast, at boring sites far from these towns, the irregular and sporadic woodland disturbances continue, although these become more common. Phase III (1500/1300–850 B.P.) shows further woodland clearance fires along the coast. Pollen, archaeological and historical evidence allow us to relate this to the introduction of grazing in the coastal area, connected with seasonal transhumance between the littoral and mountain regions. Phase IV (850–300/150 B.P.) is characterised by the final clearance of woodland and the development of olive farming. During this phase the agrarian Mediterranean landscape was definitively formed.