Population structure,diversity, and phylogeography in the near‐threatened Eurasian black vultures Aegypius monachus (Falconiformes; Accipitridae) in Europe: insights from microsatellite and mitochondrial DNA variation

The Eurasian black vulture (Aegypius monachus) has experienced a severe decline during the last two centuries and is globally classified as near‐threatened. This has led to the extinction of many traditional breeding areas in Europe and resulted in the present patchy distribution (Iberian and Balkan peninsulas) in the Western Palearctic. In the present study, we describe the current genetic status of the European populations using both mitochondrial cytochrome b sequences and nuclear microsatellite markers, comparing with those found in Asia (Mongolia and Caucasus region). Although, mitochondrial (mt)DNA revealed a relatively low genetic variability (haplotype diversity), no evidence of genome‐wide genetic erosion exists because nuclear diversity exhibits normal levels and strong differentiation. A highly philopatric dispersal behaviour must be invoked to explain the existence of a clear pattern that revealed by the phylogeographic analysis, which indicates a sharp East–West clinal distribution and an allopatric differentiation. The distribution of mtDNA haplotypes one in the Iberian population and two in Balkan population and the significance divergence at nuclear loci fulfill the definitions of those populations as evolutionary significant units. We discuss how management strategies should aim at the maintenance (or increase) of current genetic variability levels, suggesting that independent conservation plans are urgently required to protect these two breeding European populations from extinction. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 859–872.     

Genetic diversity and ecological niche modelling of the restricted Recordia reitzii (Verbenaceae) from southern Brazilian Atlantic forest

Genetic diversity analyses, coupled with ecological niche modelling (ENM) of species with a restricted distribution, may provide valuable information for understanding diversification patterns in endangered areas. We analyzed the genetic diversity of Recordia reitzii, a tree restricted to the threatened and highly fragmented Brazilian Atlantic forest, using three intergenic cpDNA spacers and ten microsatellite (SSR) loci. To assess the historical processes that may have influenced the distribution of extant R. reitzii populations, the current potential distributions of R. reitzii and Recordia boliviana, a closely related species, were modelled and projected onto the Last Glacial Maximum (LGM) and Last Interglacial (LIG) periods. Niche divergence was quantified between these two. The cpDNA and SSR data showed a north–south pattern of the diversity distribution and structured populations, suggesting that gene flow is probably limited. According to our data, R. reitzii exhibits low genetic diversity, which may be a result of a founder or distribution‐reduction effect, narrow distribution or small population size. The ecological niche models showed a wider palaeodistribution during the LIG and a retraction during the LGM for both species. Tests of niche divergence and conservatism indicated that bioclimatic factors might have influenced the diversification of these Recordia species. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 176 , 332–348.     

Long‐term isolation of the coastal plant Calystegia soldanella (Convolvulaceae) in ancient freshwater Lake Biwa,Japan

Lake Biwa is an ancient freshwater lake that was formed approximately 4 Mya and harbours many coastal plants that commonly inhabit the seashore. We used chloroplast DNA haplotype analysis using two spacer sequences and simple sequence repeat (SSR) analysis using eight nuclear microsatellite markers to detect genomic signatures indicating long‐term isolation of inland populations of Calystegia soldanella in Lake Biwa from coastal populations. We used 348 samples from 63 populations for haplotype analysis and 478 samples from 27 populations for SSR analysis covering the inland and coastal distribution of the species. We detected seven haplotypes, and the distribution pattern of these haplotypes was geographically highly structured between Lake Biwa and the coast. Nuclear SSR analysis also supported genetic differentiation between Lake Biwa and coastal populations (analyses of molecular variance, 43%), and the grouping of Lake Biwa and coastal populations by a Neighbour‐joining tree. In addition, genetic diversity of the inland populations (mean H_E = 0.153) was significantly lower than that of coastal populations (mean H_E = 0.328). These results suggested that inland populations at Lake Biwa have been isolated from coastal populations for a very long time. The inland populations most likely experienced a bottleneck effect, resulting in sufficient in situ genetic divergence to clearly distinguish them from coastal populations. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 51–66.     

Phylogeography of the humbug damselfish,Dascyllus aruanus (Linnaeus, 1758): evidence of Indo‐Pacific vicariance and genetic differentiation of peripheral populations

The phylogeographical structure of coral‐associated reef fishes may have been severely affected, more than species from deeper habitats, by habitat loss during periods of low sea level. The humbug damselfish, Dascyllus aruanus, is widely distributed across the Indo‐West Pacific, and exclusively inhabits branching corals. We used mitochondrial cytochrome b sequence and seven microsatellite loci on D. aruanus samples (260 individuals) from 13 locations across the Indo‐West Pacific to investigate its phylogeographical structure distribution‐wide. A major genetic partition was found between the Indian and Pacific Ocean populations, which we interpret as the result of geographical isolation on either side of the Indo‐Pacific barrier during glacial periods. The peripheral populations of the Red Sea and the Society Islands exhibited lower genetic diversity, and substantial genetic differences with the other populations, suggesting relative isolation. Thus, vicariance on either side of the Indo‐Pacific barrier and peripheral differentiation are considered to be the main drivers that have shaped the phylogeographical patterns presently observed in D. aruanus. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 931–942.     

Barcoding,molecular taxonomy,and exploration of the diversity of shrews (Soricomorpha: Soricidae) on Mount Nimba (Guinea)

We tested the efficiency of cytochrome oxidase I (COI)‐barcoding as a taxonomic tool to discriminate and identify sympatric shrew species on Mount Nimba (Guinea). We identified 148 specimens at the species level using morphological characters and comparison with type specimens, including several taxa from Mount Nimba. We identified ten morphospecies and tested aspects of genetic diversity and monophyly using genetic data from three mitochondrial (16S, cytochrome b, and COI) and one nuclear marker (the breast cancer gene, BRCA). Nine morphospecies were validated under the phylogenetic and genetic species concepts, including the recently diverged species Crocidura buettikoferi, Crocidura theresae, and Crocidura grandiceps. Under the same concepts, our analyses revealed the presence of two cryptic species amongst animals identified as Crocidura muricauda. We then tested the efficiency of barcoding thanks to commonly used phenetic methods, with the 148 specimens representing 11 potentially valid species based on morphological and molecular data. We show that COI‐barcoding is a powerful tool for shrew identification and can be used for taxonomic surveys. The comparison of genetic divergence values shows the presence of a barcoding gap (i.e. difference between the highest intraspecific and the lowest interspecific genetic divergence values). Given that only a few COI sequences are available for Afrotropical shrews, our work is an important step forward toward their enrichment. We also tested the efficiency of the three other sequenced markers and found that cytochrome b is as efficient as COI for barcoding shrews. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 166 , 672–687.     

Lack of association between winter coat colour and genetic population structure in the Japanese hare,Lepus brachyurus (Lagomorpha: Leporidae)

Seasonal changes in fur colour in some mammalian species have long attracted the attention of biologists, especially in species showing population variation in these seasonal changes. Genetic differences among populations that show differences in seasonal changes in coat colour have been poorly studied. Because the Japanese hare (Lepus brachyurus) has two allopatric morphotypes that show remarkably different coat colours in winter, we examined the population genetic structure of the species using partial sequences of the SRY gene and six autosomal genes: three coat colour‐related genes (ASIP, TYR, and MC1R) and three putatively neutral genes (TSHB, APOB, and SPTBN1). The phylogenetic tree of SRY sequences exhibited two distinct lineages that diverged approsimately 1 Mya. Although the two lineages exhibited a clear allopatric distribution, it was not consistent with the distribution of morphotypes. In addition, six nuclear gene sequences failed to reveal genetic differences between morphotypes. Population network trees for 11 expedient populations divided the populations into four groups. Genetic structure analysis revealed an admixture of four genetic clusters in L. brachyurus, two of which showed large genetic differences. Our results suggest ancient vicariance in L. brachyurus, and we detected no genetic differences between the two morphotypes. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 761–776.     

Molecular phylogeography of the microturbellarian Monocelis lineata (Platyhelminthes: Proseriata) in the North‐East Atlantic

Monocelis lineata is a complex of cryptic species (three in the Mediterranean and one in the Atlantic) widespread in midlittoral habitats. Throughout the range, populations with or without an ocular pigmented shield are found. We investigated the genetic structure of the North‐East Atlantic populations with the aim of shedding light on their phylogeography and reconstructing possible patterns of recolonization after the Würmian glaciation. Fourteen samples were investigated using cytochrome c oxidase subunit I (COI) and 13 by inter‐simple sequence repeats (ISSRs). COI did not exhibit a clear pattern of decreased genetic diversity along a latitudinal gradient. Populations from Ferrol (Spain), Doolin (Ireland), and Helsingør (Denmark) showed a higher genetic variability, whereas a reduction in the number of haplotypes was found at the northernmost edge of the distribution and in northern Ireland and Scotland. Two genetically differentiated areas (southern Europe and south‐western Ireland versus northern Atlantic) were revealed by ISSR data. The results obtained provided evidence of three refugia (Iberian Peninsula, south‐western Ireland, and North Sea), and the occurrence of secondary contacts that shaped the genetic variability of some of the populations examined. Two different recolonization pathways of north‐western Europe during the post‐Würmian glaciations have been detected. Furthermore, ISSR analysis provided evidence of genetic divergence among populations with and without pigmented eyespot, suggesting the action of ecological differentiation. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103 , 117–135.     

Mitochondrial DNA and nuclear microsatellites reveal high diversity and genetic structure in an avian top predator,the white‐tailed sea eagle,in central Europe

We analysed 123 white‐tailed sea eagles (Haliaeetus albicilla) from (primarily central) Europe with respect to variability and differentiation based on 499 bp of the mitochondrial control region and genotypes at seven unlinked nuclear microsatellites. Variability was high (overall expected heterozygosity, haplotype and nucleotide diversity being 0.70, 0.764 and 0.00698, respectively) and both marker systems showed a subdivision into two main genetic clusters (microsatellites) or haplogroups (mtDNA). In line with earlier analyses focusing on populations from northern and eastern Europe, as well as from Asia, we found a high level of admixture in Europe and no signs of a bottleneck – despite a severe decline of white‐tailed sea eagle populations during the 20^th century. Europe is thus a global stronghold for this species not only with respect to the number of breeding pairs but also regarding the proportion of species‐wide genetic diversity. Our dense sampling revealed a possibly clinal variation within central Europe from north‐west to south‐east that was reflected by the distribution of mtDNA haplotypes as well as the two microsatellite‐based clusters. This population differentiation in central Europe probably originated from a geographically structured postglacial colonization and was later enhanced by recent demographic fluctuations. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 727–737.     

Gut microbial diversity across a contact zone for California voles: Implications for lineage divergence of hosts and mitonuclear mismatch in the assembly of the mammalian gut microbiome

Gut microbial diversity is thought to reflect the co‐evolution of microbes and their hosts as well as current host‐specific attributes such as genetic background and environmental setting. To explore interactions among these parameters, we characterized variation in gut microbiome composition of California voles (Microtus californicus) across a contact zone between two recently diverged lineages of this species. Because this contact zone contains individuals with mismatched mitochondrial‐nuclear genomes (cybrids), it provides an important opportunity to explore how different components of the genotype contribute to gut microbial diversity. Analyses of bacterial 16S rRNA sequences and joint species distribution modelling revealed that host genotypes and genetic differentiation among host populations together explained more than 50% of microbial community variation across our sampling transect. The ranked importance (most to least) of factors contributing to gut microbial diversity in our study populations were: genome‐wide population differentiation, local environmental conditions, and host genotypes. However, differences in microbial communities among vole populations (β‐diversity) did not follow patterns of lineage divergence (i.e., phylosymbiosis). Instead, among‐population variation was best explained by the spatial distribution of hosts, as expected if the environment is a primary source of gut microbial diversity (i.e., dispersal limitation hypothesis). Across the contact zone, several bacterial taxa differed in relative abundance between the two parental lineages as well as among individuals with mismatched mitochondrial and nuclear genomes. Thus, genetic divergence among host lineages and mitonuclear genomic mismatches may also contribute to microbial diversity by altering interactions between host genomes and gut microbiota (i.e., hologenome speciation hypothesis).     

Phylogeography and genetic structure of the threatened Canarian Juniperus cedrus (Cupressaceae)

We used plastid sequences (trnL, trnL‐trnF, petN‐psbM and trnT‐trnL) to infer the phylogenetic relationships and inter‐island connections of the Canarian Juniperus cedrus, and AFLP fingerprints to assess its genetic diversity patterns. Maximum Likelihood, Maximum Parsimony and Bayesian methods suggest independent colonization events for the three Macaronesian junipers and support the monophyly of J. cedrus. Plastid sequences reveal a low genetic diversity (three haplotypes) and do not provide sufficient information to resolve its temporal and geographical origin. AFLPs indicate a greater isolation in J. cedrus than in other Macaronesian trees with similar distributions and dispersal syndromes. Gran Canaria harbours the least genetically diverse population, which justifies immediate conservation actions. This island and Tenerife also show independent genetic structure, meaning that genetic exchange from other islands should be avoided in eventual reinforcements. Populations from La Palma and La Gomera show the highest genetic diversity levels and number of polymorphic AFLPs, probably because a lower incidence of felling has allowed a less dramatic influence of genetic bottlenecks. We suggest that management efforts should prioritize populations from these islands to preserve the evolutionary potential of the species, but we also stress the importance of knowledge of the evolutionary history, genetic structure and ecological interactions in conservation strategies. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 175 , 376–394.     

Colonization history of Mallorca Island by the European rabbit,Oryctolagus cuniculus,and the Iberian hare,Lepus granatensis (Lagomorpha: Leporidae)

The Mediterranean islands have a long history of human‐mediated introductions resulting in frequent replacements of their fauna and flora. Although these histories are sometimes well documented or may be inferred from paleontological studies, the use of phylogenetic and population genetic reconstruction methods provides a complementary perspective for answering questions related to the history of insular species. In the present study, we infer the colonization history of Mallorca (Balearic Islands) by the European rabbit (Oryctolagus cuniculus) and the Iberian hare (Lepus granatensis) using sequence variation of the mitochondrial DNA control region from continental and insular specimens (total of 489 sequences). Additionally, the taxonomic identity of Mallorcan L. granatensis was confirmed using a diagnostic nuclear marker. For both Mallorcan rabbits and hares, genetic diversity was comparable to the continental populations, suggesting the introduction of multiple lineages. Two Mallorcan haplogroups were found in hares, which likely correspond to two introduction events. Rabbits from Mallorca were identified as belonging to the subspecies Oryctolagus cuniculus cuniculus, and may have been originated both from Iberian and French populations. The molecular estimates of the timing of the colonization events of the Mallorcan lagomorphs are consistent with human‐mediated introductions by early settlers on the islands. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 748–760.     

Long‐term ‘islands’ in the landscape: low gene flow,effective population size and genetic divergence in the shrub Hakea oldfieldii (Proteaceae)

The genetic structure of disjunct populations is determined by founding genetic properties, demographic processes, gene flow, drift and local selection. We aim to identify the genetic consequences of natural population disjunction at regional and local scales in Hakea oldfieldii using nuclear and plastid markers to investigate long‐term effective population sizes and gene flow, and patterns of diversity and divergence, among populations. Regional divergence was significant as shown by a consistent pattern in principal coordinates, neighbor‐joining and Bayesian analyses, but divergence at the local level was also significant with localized distribution of plastid haplotypes and populations clustering separately in Bayesian analyses. Historical, recent and first‐generation gene flow was low, suggesting that recent habitat fragmentation has not reduced gene migration significantly. Genetic bottlenecks were detected in three populations. Long‐term effective population size was significantly correlated with the number of alleles/locus and observed heterozygosity, but not with census population size, suggesting that the loss of diversity is associated with long‐term changes rather than recent fragmentation. Inbreeding coefficients were significant in only three populations, suggesting that the loss of diversity is linked to drift and bottlenecks associated with demographic processes (local extinction by fires) rather than inbreeding. Historical disjunction as a result of specific ecological requirements, contraction of habitats following drying during the Pleistocene, low gene flow and changes in population size are likely to have been important forces driving divergence through isolation by distance and drift. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 319–334.     

Spatial genetic structure in Beta vulgaris subsp. maritima and Beta macrocarpa reveals the effect of contrasting mating system,influence of marine currents,and footprints of postglacial recolonization routes

Understanding the factors that contribute to population genetic divergence across a species' range is a long‐standing goal in evolutionary biology and ecological genetics. We examined the relative importance of historical and ecological features in shaping the present‐day spatial patterns of genetic structure in two related plant species, Beta vulgaris subsp. maritima and Beta macrocarpa. Using nuclear and mitochondrial markers, we surveyed 93 populations from Brittany (France) to Morocco – the southern limit of their species' range distribution. Whereas B. macrocarpa showed a genotypic structure and a high level of genetic differentiation indicative of selfing, the population genetic structure of B. vulgaris subsp. maritima was consistent with an outcrossing mating system. We further showed (1) a strong geographic clustering in coastal B. vulgaris subsp. maritima populations that highlighted the influence of marine currents in shaping different lineages and (2) a peculiar genetic structure of inland B. vulgaris subsp. maritima populations that could indicate the admixture of distinct evolutionary lineages and recent expansions associated with anthropogenic disturbances. Spatial patterns of nuclear diversity and differentiation also supported a stepwise recolonization of Europe from Atlantic‐Mediterranean refugia after the last glacial period, with leading‐edge expansions. However, cytoplasmic diversity was not impacted by postglacial recolonization: stochastic long‐distance seed dispersal mediated by major oceanic currents may mitigate the common patterns of reduced cytoplasmic diversity observed for edge populations. Overall, the patterns we documented here challenge the general view of reduced genetic diversity at the edge of a species' range distribution and provide clues for understanding how life‐history and major geographic features interact to shape the distribution of genetic diversity.     

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.     

Land masses and oceanic currents drive population structure of Heritiera littoralis,a widespread mangrove in the Indo‐West Pacific

Phylogeographic forces driving evolution of sea‐dispersed plants are often influenced by regional and species characteristics, although not yet deciphered at a large spatial scale for many taxa like the mangrove species Heritiera littoralis. This study aimed to assess geographic distribution of genetic variation of this widespread mangrove in the Indo‐West Pacific region and identify the phylogeographic factors influencing its present‐day distribution. Analysis of five chloroplast DNA fragments’ sequences from 37 populations revealed low genetic diversity at the population level and strong genetic structure of H. littoralis in this region. The estimated divergence times between the major genetic lineages indicated that glacial level changes during the Pleistocene epoch induced strong genetic differentiation across the Indian and Pacific Oceans. In comparison to the strong genetic break imposed by the Sunda Shelf toward splitting the lineages of the Indian and Pacific Oceans, the genetic differentiation between Indo‐Malesia and Australasia was not so prominent. Long‐distance dispersal ability of H. littoralis propagules helped the species to attain transoceanic distribution not only across South East Asia and Australia, but also across the Indian Ocean to East Africa. However, oceanic circulation pattern in the South China Sea was found to act as a barrier creating further intraoceanic genetic differentiation. Overall, phylogeographic analysis in this study revealed that glacial vicariance had profound influence on population differentiation in H. littoralis and caused low genetic diversity except for the refugia populations near the equator which might have persisted through glacial maxima. With increasing loss of suitable habitats due to anthropogenic activities, these findings therefore emphasize the urgent need for conservation actions for all populations throughout the distribution range of H. littoralis.     

Patterns and processes in the evolutionary history of parrotfishes (Family Labridae)

Phylogenetic reconstruction of the evolutionary relationships among 61 of the 70 species of the parrotfish genera Chlorurus and Scarus (Family Labridae) based on mitochondrial and nuclear gene sequences retrieved 15 well‐supported clades with mid Pliocene/Pleistocene diversification. Twenty‐two reciprocally monophyletic sister‐species pairs were identified: 64% were allopatric, and the remainder were sympatric. Age of divergence was similar for allopatric and sympatric species pairs. Sympatric sister pairs displayed greater divergence in morphology, ecology, and sexually dimorphic colour patterns than did allopatric pairs, suggesting that both genetic drift in allopatric species pairs and ecologically adaptive divergence between members of sympatric pairs have played a role in diversification. Basal species typically have small geographical ranges and are restricted to geographically and ecologically peripheral reef habitats. We found little evidence that a single dominant process has driven diversification, nor did we detect a pattern of discrete, sequential stages of diversification in relation to habitat, ecology, and reproductive biology. The evolution of Chlorurus and Scarus has been complex, involving a number of speciation processes. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 107 , 529–557.     

Phylogeography of the Poecilimon luschani species group (Orthoptera,Tettigoniidae): a radiation strictly correlated with climatic transitions in the Pleistocene

The Quaternary biogeography of Anatolia has received considerable interest recently. Here, the genealogical history of the Anatolio?Balkan lineage of the Poecilimon luschani species group was evaluated. Using concatenated data from 16S rDNA and cytochrome c oxidase subunit I (COI) sequences, the timings of inter‐ and intraspecies radiations were estimated. The demographic history of the populations was estimated using a data set established from COI sequences. Genetic diversity was very high in almost all of the populations studied. Fixation indices suggested extreme divergence of P. luschani. A molecular chronogram estimated a radiation history for the species/subspecies over a period ranging from 1.323 to 0.440 Myr. Demographic analyses applied to 11 populations suggested departure in population size for most of the local populations. The following conclusions were reached: (1) P. luschani originated from an Anatolio‐Aegean ancestral stock and extended its range to the Balkans through Dardanelles during the Early Pleistocene; (2) the Mid‐Pleistocene Transition, the lengthening of glacial period from 41 to 100 Kyr and the initiation of intense glaciation periods are the three main events corresponding to the main nodes of the chronogram; (3) altitudinal heterogeneity played a buffer role during the glacial cycles, allowing populations to cope with severe environmental changes; (4) the effects of Pleistocene climate cycles on populations differ according to altitudinal and latitudinal location in Anatolia, and (5) habitat preferences, such as altitudinal range, may easily shift because of changes in environmental conditions. © 2014 The Linnean Society of London     

Cryptic speciation along a bathymetric gradient

The deep ocean supports a highly diverse and mostly endemic fauna, yet little is known about how or where new species form in this remote ecosystem. How speciation occurs is especially intriguing in the deep sea because few obvious barriers exist that would disrupt gene flow. Geographic and bathymetric patterns of genetic variation can provide key insights into how and where new species form. We quantified the population genetic structure of a protobranch bivalve, Neilonella salicensis, along a depth gradient (2200–3800 m) in the western North Atlantic using both nuclear (28S and calmodulin intron) and mitochondrial (cytochrome c oxidase subunit I) loci. A sharp genetic break occurred for each locus between populations above 2800 m and below 3200 m, defining two distinct clades with no nuclear or mitochondrial haplotypes shared between depth regimes. Bayesian phylogenetic analyses provided strong support for two clades, separated by depth, within N. salicensis. Although no morphological divergence was apparent, we suggest that the depth‐related population genetic and phylogenetic divergence is indicative of a cryptic species. The frequent occurrence of various stages of divergence associated with species formation along bathymetric gradients suggests that depth, and the environmental gradients that attend changes in depth, probably play a fundamental role in the diversification of marine organisms, especially in deep water. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 897–913.