Mitochondrial phylogeography and subspecies of the wide‐ranging sub‐Saharan leopard tortoise Stigmochelys pardalis (Testudines: Testudinidae) – a case study for the pitfalls of pseudogenes and GenBank sequences

The leopard tortoise (Stigmochelys pardalis) is the most widely distributed sub‐Saharan tortoise species, with a range extending from the Horn of Africa all over eastern Africa to the Republic of South Africa, Namibia and southernmost Angola. Using 1938 bp of mitochondrial DNA (cyt b gene, partial ND4 gene plus adjacent tRNA genes) from a nearly range‐wide sampling, we examined its phylogeographic structure and compared our findings with previously published GenBank sequences. We identified seven major clades that are largely parapatrically distributed. A few records of distinct haplotypes at the same locality or in close proximity could be the result of translocation of tortoises by man. The greatest diversity occurs in the south of the species’ range, with five out of the seven clades. Testing for isolation‐by‐distance suggests that the observed phylogeographic structure is the result of restricted geographical gene flow and not of historical vicariance. This is in sharp contrast to wide‐ranging thermophilic reptiles from the western Palaearctic, whose phylogeographic structure was significantly shaped by Pleistocene range interruptions, but also by earlier dispersal and vicariant events. Most cyt b sequences of S. pardalis from GenBank turned out to be nuclear pseudogenes, or to be of chimerical origin from such pseudogenes and authentic mitochondrial sequences, which argues for caution regarding uncritical usage of GenBank sequences. The recent revalidation of the two subspecies of S. pardalis was based on such a chimerical sequence that was erroneously identified with the subspecies S. p. babcocki. Furthermore, according to our data, the distribution of mitochondrial clades does match neither the traditional subspecies ranges nor the pronounced geographical size variation of leopard tortoises. We conclude that there is no rationale for recognizing subspecies within S. pardalis.     

Mitochondrial DNA sequence variability in red-legged partridge, <Emphasis Type="Italic">Alectoris rufa</Emphasis>, Spanish populations and the origins of genetic contamination from <Emphasis Type="Italic">A. chukar</Emphasis>

The analysis of 135 mitochondrial D-loop sequences of the Iberian autochthonous red-legged partridge (Alectoris rufa) from wild population hunting bags from various locations and fowl runs in Spain yielded 37 different haplotypes. Among these, three haplotypes correspond to chukar partridges (Alectoris chukar), indicating genetic introgression from birds illegally introduced for restocking: three individuals carrying such haplotypes where found in natural populations, one appeared among those sampled on a mass reproduction farm and the remaining 10 in another fowl-run. The geographical origin of the contaminating chukar haplotypes could be assigned to the most easterly area of the chukar partridge geographical distribution in China. Molecular diversity parameters in the A. rufa samples indicate a considerable amount of genetic variation. Φ_ST showed significant differences among populations that are not explained by geographical distance alone. Particularly, one northern population (Palencia) shows a certain degree of genetic differentiation that could reflect a previously suggested subspecies division. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Colour polymorphism in Salamandra salamandra (Amphibia: Urodela), revealed by a lack of genetic and environmental differentiation between distinct phenotypes

The existence of two or more distinctly coloured phenotypes among individuals of an interbreeding population is known as colour polymorphism. In amphibians, this phenomenon is pervasive among anurans, but rare or absent among salamanders and caecilians, respectively. Here, we examine whether various distinct phenotypes of Salamandra salamandra in North Spain, used as a basis to describe the subspecies S. s. bernardezi and S. s. alfredschmidti, indeed warrant separate taxonomic status or that these co‐occur and belong to a single taxon. Based on a sample of 1147 individuals from 27 local populations, six phenotype classes were designated. Although two phenotypes that are attributable to S. s. alfredschmidti show some degree of geographical restriction, these co‐occur with those representing typical S. s. bernardezi. A fifth phenotype class could not be unambiguously attributed to either subspecies due to an overlap in previously suggested diagnostic characteristics. Mitochondrial (cytochrome b) and nuclear (β‐fibrinogen) DNA analyses revealed S. s. alfredschmidti to be nested within several subclades of S. s. bernardezi, without displaying unique lineages. Furthermore, no significant divergence was recovered by means of niche overlap analyses. As a result, we revoke the subspecies status of S. s. alfredschmidti, which should be regarded as a junior synonym of S. s. bernardezi. The current findings confirm the existence of colour polymorphism in S. salamandra and the family Salamandridae, which provides exciting possibilities for future research.     

Biometrical analysis reveals major differences between the two subspecies of the European rabbit

The climatic oscillations that have occurred in the last few million years have strongly affected species distribution ranges. Highly divergent genetic lineages arose, some of which correspond to recognized subspecies that currently occupy small geographical areas. Understanding the implications of the genetic differences between these subspecies is crucial for proper conservation of Evolutionarily Significant Units. We use the two European rabbit subspecies, Oryctolagus cuniculus cuniculus and O. c. algirus, in the Iberian Peninsula as a model to investigate the repercussions at the biometric level of their largely recognized genetic differentiation. To accomplish this we analysed the ear and hind foot length, and the body mass of 999 adult rabbits from 27 locations across the distribution range of both subspecies in their native range, the Iberian Peninsula. Our results show biometric differences between the two subspecies, also explained by geographical location and sex, O. c. algirus being lighter and having shorter ear and hind foot lengths. We examine these findings under an evolutionary framework, and discuss their implications for current conservation efforts. Future research should focus on the ecological implications of these biometric differences, namely potential different habitat use and anti‐predatory strategies in the species' native range.     

Mitochondrial DNA variation in the caramote prawn <Emphasis Type="Italic">Penaeus (Melicertus) kerathurus</Emphasis> across a transition zone in the Mediterranean Sea

In this study we analysed mitochondrial DNA variation in Penaeus kerathurus prawns collected from seven locations along a transect across the Siculo–Tunisian region in order to verify if any population structuring exists over a limited geographical scale and to delineate the putative transition zone with sufficient accuracy. Partial DNA sequences of COI and 16S genes were analysed. In contrast to the highly conservative 16S gene, the COI sequences exhibited sufficient diversity for population analysis. The COI gene revealed low levels of haplotype and nucleotide diversities. The size of the annual landings of this commercial species suggests large population sizes. Hence, the low genetic diversity detected in this study could indicate a possible reduction in effective population sizes in the past. We detected significant genetic differentiation between eastern and western populations likely due to restricted gene flow across the Siculo–Tunisian boundary. We discuss the different evolutionary forces that may have shaped the genetic variation and suggest that the genetic divide is probably maintained by present-day dispersal limitation. R. Zitari-Chatti and N. Chatti are contributed equally to the work.     

Phylogeny and biogeography of the Japanese rhinoceros beetle,Trypoxylus dichotomus (Coleoptera: Scarabaeidae) based on SNP markers

The Japanese rhinoceros beetle Trypoxylus dichotomus is one of the largest beetle species in the world and is commonly used in traditional Chinese medicine. Ten subspecies of T. dichotomus and a related Trypoxylus species (T. kanamorii) have been described throughout Asia, but their taxonomic delimitations remain problematic. To clarify issues such as taxonomy, and the degree of genetic differentiation of Trypoxylus populations, we investigated the genetic structure, genetic variability, and phylogeography of 53 specimens of Trypoxylus species from 44 locations in five Asian countries (China, Japan, Korea, Thailand, and Myanmar). Using specific‐locus amplified fragment sequencing (SLAF‐seq) techniques, we developed 330,799 SLAFs over 114.16M reads, in turn yielding 46,939 high‐resolution single nucleotide polymorphisms (SNPs) for genotyping. Phylogenetic analysis of SNPs indicated the presence of three distinct genetic groups, suggesting that the various subspecies could be treated as three groups of populations. PCA and ADMIXTURE analysis also identified three genetic clusters (North, South, West), which corresponded to their locations, suggesting that geographic factors were important in maintaining within population homogeneity and between population divergence. Analyses of SNP data confirmed the monophyly of certain subspecies on islands, while other subspecies (e.g., T. d. septentrionalis) were found to be polyphyletic and nested in more than one lineage. AMOVA demonstrated high level of differentiation among populations/groups. Also, pairwise F_ST values revealed high differentiation, particularly between South and West, as well as between North and South. Despite the differentiation, measurable gene flow was inferred between genetic clusters but at varying rates and directions. Our study demonstrated that SLAF‐seq derived markers outperformed 16S and COII sequences and provided improved resolution of the genetic differentiation of rhinoceros beetle populations from a large part of the species’ range.     

Subterranean termite phylogeography reveals multiple postglacial colonization events in southwestern Europe

A long‐standing goal of evolutionary biology is to understand how paleoclimatic and geological events shape the geographical distribution and genetic structure within and among species. Using a diverse set of markers (cuticular hydrocarbons, mitochondrial and nuclear gene sequences, microsatellite loci), we studied Reticulitermes grassei and R. banyulensis, two closely related termite species in southwestern Europe. We sought to clarify the current genetic structure of populations that formed following postglacial dispersal from refugia in southern Spain and characterize the gene flow between the two lineages over the last several million years. Each marker type separately provided a fragmented picture of the evolutionary history at different timescales. Chemical analyses of cuticular hydrocarbons and phylogenetic analyses of mitochondrial and nuclear genes showed clear separation between the species, suggesting they diverged following vicariance events in the Late Miocene. However, the presence of intermediate chemical profiles and mtDNA introgression in some Spanish colonies suggests ongoing gene flow. The current genetic structure of Iberian populations is consistent with alternating isolation and dispersal events during Quaternary glacial periods. Analyses of population genetic structure revealed postglacial colonization routes from southern Spain to France, where populations underwent strong genetic bottlenecks after traversing the Pyrenees resulting in parapatric speciation.     

Population genetic structure and phylogeographical pattern of rice grasshopper, <Emphasis Type="Italic">Oxya hyla intricata</Emphasis>, across Southeast Asia

The rice grasshopper, Oxya hyla intricata, is a rice pest in Southeast Asia. In this study, population genetic diversity and structure of this Oxya species was examined using both DNA sequences and AFLP technology. The samples of 12 populations were collected from four Southeast Asian countries, among which 175 individuals were analysed using mitochondrial DNA cytochrome c oxidase subunit I (COI) sequences, and 232 individuals were examined using amplified fragment length polymorphisms (AFLP) to test whether the phylogeographical pattern and population genetics of this species are related to past geological events and/or climatic oscillations. No obvious trend of genetic diversity was found along a latitude/longitude gradient among different geographical groups. Phylogenetic analysis indicated three deep monophyletic clades that approximately correspond to three geographical regions separated by high mountains and a deep strait, and TCS analysis also revealed three disconnected networks, suggesting that spatial and temporal separations by vicariance, which were also supported by AMOVA as a source of the molecular variance presented among groups. Gene flow analysis showed that there had been frequent historical gene flow among local populations in different regions, but the networks exhibited no shared haplotype among populations. In conclusion, the past geological events and climatic fluctuations are the most important factor on the phylogeographical structure and genetic patterns of O. hyla intricata in Southeast Asia. Habitat, vegetation, and anthropogenic effect may also contribute to gene flow and introgression of this species. Moreover, temperature, abundant rainfall and a diversity of graminaceous species are beneficial for the migration of O. hyla intricata. High haplotype diversity, deep phylogenetic division, negative Fu’s F _s values and unimodal and multimodal distribution shapes all suggest a complicated demographic expansion pattern of these O. hyla intricata populations, which might have been caused by climatic oscillations during glacial periods in the Quaternary.     

Genetic diversity,genetic structure and phylogeography of the Iberian endemic Gypsophila struthium (Caryophyllaceae) as revealed by AFLP and plastid DNA sequences: connecting habitat fragmentation and diversification

Iberian gypsum outcrops are highly fragmented and ecologically challenging environments for plant colonization. As gypsophytes occur exclusively in such habitats, they are ideal models for the study of both the effects of habitat fragmentation and selection on population genetic diversity and structure. In this study, we used amplified fragment length polymorphism (AFLP) and plastid DNA sequences to investigate the phylogeographical history of the Iberian plant Gypsophila struthium (Caryophyllaceae), a widespread endemic restricted to Iberian gypsum outcrops. Gypsophila struthium consists of two subspecies that differ in the architecture of their inflorescence and have mostly allopatric ranges. Gypsophila struthium subsp. struthium occurs in central, eastern and south‐eastern Iberia, whereas G. struthium subsp. hispanica occurs in northern and eastern areas. AFLPs revealed low but significant genetic differentiation between the subspecies, probably as a result of a recent diversification during the Pliocene–Pleistocene. In the geographical contact zone between the taxa, the Bayesian analyses revealed populations with mixed ancestries and genetic clusters predominantly of one or the other subspecies, indicating incomplete reproductive barriers between them. Plastid DNA haplotypes revealed strong geographical structure and testified to processes of isolation by distance and continuous range expansion for some haplotype clades. The Bayesian analyses of the population structure of AFLP data and nested clade phylogeographical analysis (NCPA) of plastid haplotypes revealed that the putative ancestral range corresponded to central and eastern populations of G. struthium subsp. struthium, with those lineages contributing through more recent expansion to increased genetic diversity and structure of the south‐eastern and eastern ranges of this subspecies and to the diversification of G. struthium subsp. hispanica in northern and eastern gypsum outcrops. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2013, 173 , 654–675.     

Phylogeography of the Eurasian green woodpecker (Picus viridis)

Aim In this paper we investigate the evolutionary history of the Eurasian green woodpecker (Picus viridis) using molecular markers. We specifically focus on the respective roles of Pleistocene climatic oscillations and geographical barriers in shaping the current population genetics within this species. In addition, we discuss the validity of current species and subspecies limits. Location Western Palaearctic: Europe to western Russia, and Africa north of the Sahara. Methods We sequenced two mitochondrial genes and five nuclear introns for 17 Eurasian green woodpeckers. Multilocus phylogenetic analyses were conducted using maximum likelihood and Bayesian algorithms. In addition, we sequenced a fragment of the cytochrome b gene (cyt b, 427 bp) and of the Z‐linked BRM intron 15 for 113 and 85 individuals, respectively. The latter data set was analysed using population genetic methods. Results Our phylogenetic results support the monophyly of Picus viridis and suggest that this taxon comprises three allopatric/parapatric lineages distributed in North Africa, the Iberian Peninsula and Europe, respectively. The North African lineage split from the Iberian/European clade during the early Pleistocene (1.6–2.2 Ma). The divergence event between the Iberian and the European lineages occurred during the mid‐Pleistocene (0.7–1.2 Ma). Our results also support a post‐glacial range expansion of these two lineages from distinct refugia located in the Iberian Peninsula and possibly in eastern Europe or Anatolia, which led to the establishment of a secondary contact zone in southern France. Main conclusions Our results emphasize the crucial role of both Pleistocene climatic oscillations and geographical barriers (Strait of Gibraltar, Pyrenees chain) in shaping the current genetic structure of the Eurasian green woodpecker. Our molecular data, in combination with diagnosable plumage characters, suggest that the North African green woodpecker (Levaillant’s woodpecker) merits species rank as Picus vaillantii (Malherbe, 1847). The two European lineages could be distinguished by molecular and phenotypic characters over most of their respective geographical ranges, but they locally exchange genes in southern France. Consequently, we prefer to treat them as subspecies (P. viridis viridis, P. viridis sharpei) pending further studies.     

PHYLOGEOGRAPHY OF SPOTTED OWL (STRIX OCCIDENTALIS) POPULATIONS BASED ON MITOCHONDRIAL DNA SEQUENCES: GENE FLOW,GENETIC STRUCTURE,AND A NOVEL BIOGEOGRAPHIC PATTERN

Mitochondrial DNA control region sequences of spotted owls (Strix occidentalis) allowed us to investigate gene flow, genetic structure, and biogeographic relationships among these forest-dwelling birds of western North America Estimates of gene flow based on genetic partitioning and the phylogeography of haplotypes indicate substantial dispersal within three long-recognized subspecies. However, patterns of individual phyletic relationships indicate a historical absence of gene flow among the subspecies, which are essentially monophyletic. The pattern of haplotype coalescence enabled us to identify the approximate timing and direction of a recent episode of gene flow from the Sierra Nevada to the northern coastal ranges. The three subspecies comprise phylogenetic species, and the northern spotted owl (S. o. caurina) is sister to a clade of California (S. o. occidentalis) plus Mexican spotted owls (S o lucida); this represents a novel biogeographic pattern within birds. The California spotted owl had substantially lower nucleotide diversity than the other two subspecies; this result is inconsistent with present patterns of population density A causal explanation requires postulating a severe bottleneck or a selective sweep, either of which was confined to only one geographic region.     

Phylogeny,molecular ecology and taxonomy of southern Iberian lineages of <Emphasis Type="Italic">Triops mauritanicus</Emphasis> (Crustacea: Notostraca)

We investigated the phylogeography of the main lineages in the tadpole shrimp Triops mauritanicus Ghigi in the south-western Iberian Peninsula, using mitochondrial 12S and 16S rDNA sequences. Our results indicate that a fourth, hitherto unknown main phylogenetic lineage occurs in Iberia, so that in total, the species is divided into six distinct clades, comprising T. m. mauritanicus, T. m. simplex Ghigi, and four as yet unnamed lineages that appear to be endemic to Iberia. Percentages of sequence divergence among the main clades in T. mauritanicus reach the range reported for recognized species in other notostracan lineages. A thorough morphological investigation also revealed that the differentiation among these lineages is higher than previously thought, and that populations of three of the main clades within T. mauritanicus can be reliably separated from each other and from the remaining lineages based on the morphology of adult males. The remaining clades also show a significant level of morphological differentiation, but include a certain proportion of populations for which the additional application of molecular methods is needed for a reliable determination. The geographic distributions of 12S haplotypes are indicative of frequent dispersal events and gene flow among populations belonging to the same main lineage, but give no evidence of recent migration events among different main lineages, suggesting that there is no gene flow among the latter. Our data thus suggest that the six main lineages within T. mauritanicus represent distinct species. We therefore describe the Iberian lineages as T. baeticus Korn n. sp., T. emeritensis Korn & Pérez-Bote n. sp., T. gadensis Korn & García-de-Lomas n. sp., and T. vicentinus Korn, Machado, Cristo & Cancela da Fonseca n. sp., and reinstate T. simplex Ghigi to full species status. Our data confirm the general, previously recognized pattern of a lower dispersal probability in gonochoric Triops taxa. However, we found evidence that passive dispersal in Triops may be further complicated by a strong habitat dependence of dispersal probability, mediated by prevailing dispersal vectors.     

Present and past genetic connectivity of the Indo‐Pacific tropical eel Anguilla bicolor

Aim The objective of this study was to reveal the present population structure and infer the gene‐flow history of the Indo‐Pacific tropical eel Anguilla bicolor. Location The Indo‐Pacific region. Methods The entire mitochondrial control region sequence and the genotypes at six microsatellite loci were analysed for 234 specimens collected from eight representative localities where two subspecies have been historically designated. In order to infer the population structure, genetic differentiation estimates, analysis of molecular variance and gene‐tree reconstruction were performed. The history of migration events and population growth was assessed using neutrality tests based on allelic frequency spectrum, coalescent‐based estimation of gene flow and Bayesian demographic analysis using control region sequences. Results Population structure analysis showed genetic divergence between eels from the Indian and Pacific oceans (F_ST = 0.0174–0.0251, P < 0.05 for microsatellites; Φ_ST = 0.706, P < 0.001 for control region), while no significant variation was observed within each ocean. Two mitochondrial sublineages that do not coincide with geographical regions were found in the Indian Ocean clade of a gene tree. However, these two sublineages were not differentiated at the microsatellite markers. The estimation of mitochondrial gene‐flow history suggested allopatric isolation between the Indian and Pacific oceans, and a possible secondary contact within the Indian Ocean after an initial population splitting. Bayesian demographic history reconstruction and neutrality tests indicated population growth in each ocean after the Indo‐Pacific divergence. Main conclusions Anguilla bicolor has diverged between the Indian and Pacific oceans, which is consistent with the classical subspecies designation, but is apparently genetically homogeneous in the Indian Ocean. The analysis of gene‐flow and demographic history indicated that the two mitochondrial sublineages observed in the Indian Ocean probably represent the haplotype groups of relict ancestral populations. A comparison with a sympatric congener suggested that absolute physical barriers to gene flow may not be necessary for population divergence in eels.     

Unveiling introgression in bumblebee (Bombus terrestris) populations through mitogenome-based markers

The bumblebee, Bombus terrestris, is an important pollinator commercially used on a global scale. The exported subspecies B. t. terrestris has colonised diverse environments, in some cases displacing wild pollinators to the verge of local extinction. In this sense, the native Iberian subspecies B. t. lusitanicus may be threatened by the subspecies B. t. terrestris, naturally distributed from the Pyrenees to Central Europe but also observed in southern Spain due to escapes from commercial nests. Mitochondrial genomes have a low recombination rate and a small effective population size owing to their maternal inheritance, thus providing an accurate approach to study hybridisation events between populations. Therefore, we present the sequences of the mitogenomes of both subspecies as a molecular framework to select suitable markers to detect possible introgression events between them. We used metagenomics to obtain approximately 17 kbp of the mitogenome from both subspecies. Their mitogenomes differed in 358 bp (excluding the AT-rich region). Four mitogenomic fragments were selected to be tested as subspecific diagnostic markers. A RFLP detected in the gene nad2 (NADH dehydrogenase subunit 2) has proven to be an efficient, quick and cost-effective tool to assess the dispersion of the non-endemic subspecies into Iberian native populations. Subspecific haplotypes were observed in both morphological subspecies, suggesting introgression events in the northern natural contact area and in the new human-mediated contact area in the south of the Iberian Peninsula.     

Genetic and phenotypic variation among geographically isolated populations of the globally threatened Dupont's lark Chersophilus duponti

Identifying genetically and phenotypically distinct populations of threatened species is critical if we are to delineate appropriate plans for their conservation. We conducted an integrated analysis of population genetic structure, historical demographic events, current gene flow (all based on mtDNA sequences) and morphological variation of three geographically separated groups of populations of Dupont's lark Chersophilus duponti, located in the Iberian Peninsula (three populations), Morocco (two populations), and Tunisia (one population). Unusually, this lark species is the only one among the genus Chersophilus. Our results revealed the early historical divergence of an eastern Dupont's lark lineage (in Tunisia) and a western lineage (in Morocco and Spain), consistent with subspecies taxonomy and distribution. The western lineage subsequently split into two lineages, following the isolation of Iberian and African populations. Such pattern of historical differentiation caused great population genetic structure, with differences among geographic areas explaining more than 80% of total genetic variation. Mismatch distributions and coalescent estimates of divergence time showed that lineage divergence was associated with sudden population expansion events, which apparently took place during the last glaciation, when steppe habitats were widespread across the Mediterranean region. Extant populations from different geographic areas hardly shared any haplotype (only one out of 16 ND2 haplotypes was shared by Tunisian and Moroccan Dupont's larks), and consequently gene flow between geographic areas was found to be virtually absent. Apart from showing great genetic differentiation, Dupont's larks from different geographic areas were morphologically distinct, showing substantial variation in body size and feeding-related traits (length of feet and bill). We conclude that Dupont's lark populations isolated in the Iberian Peninsula, Morocco, and Tunisia are distinct evolutionary entities and should be considered as such in conservation plans. Such circumstance sets a daunting conservation challenge that exemplifies the need of incorporating knowledge of historical processes to our general understanding of the demography of threatened species.     

Genetic diversity of the imperilled bath sponge Spongia officinalis Linnaeus, 1759 across the Mediterranean Sea: patterns of population differentiation and implications for taxonomy and conservation

The Mediterranean bath sponge Spongia officinalis is an iconic species with high socio‐economic value and precarious future owing to unregulated harvesting, mortality incidents and lack of established knowledge regarding its ecology. This study aims to assess genetic diversity and population structure of the species at different geographical scales throughout its distribution. For this purpose, 11 locations in the Eastern Mediterranean (Aegean Sea), Western Mediterranean (Provence coast) and the Strait of Gibraltar were sampled; specimens were analysed using partial mitochondrial cytochrome oxidase subunit I (COI) sequences, along with a set of eight microsatellite loci. According to our results (i) no genetic differentiation exists among the acknowledged Mediterranean morphotypes, and hence, S. officinalis can be viewed as a single, morphologically variable species; (ii) a notable divergence was recorded in the Gibraltar region, indicating the possible existence of a cryptic species; (iii) restriction to gene flow was evidenced between the Aegean Sea and Provence giving two well‐defined regional clusters, thus suggesting the existence of a phylogeographic break between the two systems; (iv) low levels of genetic structure, not correlated to geographical distance, were observed inside geographical sectors, implying mechanisms (natural or anthropogenic) that enhance dispersal and gene flow have promoted population connectivity; (v) the genetic diversity of S. officinalis is maintained high in most studied locations despite pressure from harvesting and the influence of devastating epidemics. These findings provide a basis towards the effective conservation and management of the species.     

Patterns of nuclear and mitochondrial DNA variation in Iberian populations of <Emphasis Type="Italic">Emys orbicularis</Emphasis> (Emydidae): conservation implications

The European pond turtle (Emys orbicularis) is threatened and in decline in several regions of its natural range, due to habitat loss combined with population fragmentation. In this work, we have focused our efforts on studying the genetic diversity and structure of Iberian populations with a fine-scale sampling (254 turtles in 10 populations) and a representation from North Africa and Balearic island populations. Using both nuclear and mitochondrial markers (seven microsatellites, ∼1048 bp nDNA and ∼1500 bp mtDNA) we have carried out phylogenetic and demographic analyses. Our results show low values of genetic diversity at the mitochondrial level although our microsatellite dataset revealed relatively high levels of genetic variability with a latitudinal genetic trend decreasing from southern to northern populations. A moderate degree of genetic differentiation was estimated for Iberian populations (genetic distances, F _ST values and clusters in the Bayesian analysis). The results in this study combining mtDNA and nDNA, provide the most comprehensive population genetic data for E. orbicularis in the Iberian Peninsula. Our results suggest that Iberian populations within the Iberian–Moroccan lineage should be considered as a single subspecies with five management units, and emphasize the importance of habitat management rather than population reinforcement (i.e. captive breeding and reintroduction) in this long-lived species.     

Relative role of life-history traits and historical factors in shaping genetic population structure of sardines (<Emphasis Type="Italic">Sardina pilchardus</Emphasis>)

Background  

Marine pelagic fishes exhibit rather complex patterns of genetic differentiation, which are the result of both historical processes and present day gene flow. Comparative multi-locus analyses based on both nuclear and mitochondrial genetic markers are probably the most efficient and informative approach to discerning the relative role of historical events and life-history traits in shaping genetic heterogeneity. The European sardine (Sardina pilchardus) is a small pelagic fish with a relatively high migratory capability that is expected to show low levels of genetic differentiation among populations. Previous genetic studies based on meristic and mitochondrial control region haplotype frequency data supported the existence of two sardine subspecies (S. p. pilchardus and S. p. sardina).     

Fuzzy Boundaries: Color and Gene Flow Patterns among Parapatric Lineages of the Western Shovel-Nosed Snake and Taxonomic Implication

Accurate delineation of lineage diversity is increasingly important, as species distributions are becoming more reduced and threatened. During the last century, the subspecies category was often used to denote phenotypic variation within a species range and to provide a framework for understanding lineage differentiation, often considered incipient speciation. While this category has largely fallen into disuse, previously recognized subspecies often serve as important units for conservation policy and management when other information is lacking. In this study, we evaluated phenotypic subspecies hypotheses within shovel-nosed snakes on the basis of genetic data and considered how evolutionary processes such as gene flow influenced possible incongruence between phenotypic and genetic patterns. We used both traditional phylogenetic and Bayesian clustering analyses to infer range-wide genetic structure and spatially explicit analyses to detect possible boundary locations of lineage contact. Multilocus analyses supported three historically isolated groups with low to moderate levels of contemporary gene exchange. Genetic data did not support phenotypic subspecies as exclusive groups, and we detected patterns of discordance in areas where three subspecies are presumed to be in contact. Based on genetic and phenotypic evidence, we suggested that species-level diversity is underestimated in this group and we proposed that two species be recognized, Chionactis occipitalis and C. annulata. In addition, we recommend retention of two subspecific designations within C. annulata (C. a. annulata and C. a. klauberi) that reflect regional shifts in both genetic and phenotypic variation within the species. Our results highlight the difficultly in validating taxonomic boundaries within lineages that are evolving under a time-dependent, continuous process.