Cryptic speciation in a benthic isopod from Patagonian and Falkland Island waters and the impact of glaciations on its population structure
© Leese et al; licensee BioMed Central Ltd. 2008
Received: 26 May 2008
Accepted: 19 December 2008
Published: 19 December 2008
The Falkland Islands and Patagonia are traditionally assigned to the Magellan Biogeographic Province. Most marine species in Falkland waters are also reported from southern Patagonia. It remains unclear if relatively immobile, marine benthic, shallow-water species maintain gene flow, and by what mechanism. Recurrent fluctuations in sea level during glacial cycles are regarded as a possible mechanism that might have allowed genetic exchange between the regions. However, the realized genetic exchange between the Falkland Islands and Patagonia has never been estimated.
This study analyses the genetic structure of three populations of the marine shallow-water isopod Serolis paradoxa (Fabricius, 1775) from the Falkland Islands and southern Patagonia (central Strait of Magellan and the Atlantic opening) applying seven nuclear microsatellites and a fragment of the mitochondrial 16S rRNA gene. Both marker systems report highest genetic diversity for the population from the central Strait of Magellan and lowest for the Falkland Islands. The estimated effective population sizes were large for all populations studied. Significant differentiation was observed among all three populations. The magnitude of differentiation between Patagonia and the Falkland Islands (16S: uncorrected p-distance 2.1%; microsatellites: standardized F'ST > 0.86) was an order of magnitude higher than between populations from within Patagonia. This indicates that there is currently no effective gene flow for nominal S. paradoxa between these two regions and it has been absent for time exceeding the last glacial maximum. We argue that specimens from the Strait of Magellan and the Falkland Islands very likely represent two distinct species that separated in the mid-Pleistocene (about 1 MY BP).
The results of this study indicate limited gene flow between distant populations of the brooding isopod Serolis paradoxa. The patterns of genetic diversity suggest that the only recently inundated Strait of Magellan was colonized by different source populations, most likely from Atlantic and Pacific coastal waters. Our results demonstrate that more systematic testing of shared faunal inventory and realized genetic exchange between Patagonia and the Falkland Islands is needed before a consensus concerning the position of the Falkland Islands relative to the Magellan zoogeographic province can be reached.
In this context, the relatively few reports of species endemic to the Falkland Islands are not unexpected [8, 20–22]. On the whole, evidence from marine species supports that the Falkland Islands form a part of the Atlantic Magellan Biogeographic Province and that migration of species between the continental South America and the Falkland Island is occurring repeatedly. However, recent molecular studies have shown that unrecognized cryptic species may be more common than previously assumed [15, 17, 23–33]. They indicate that morphological and ecological similarity may be an unreliable piece of evidence on which to base taxonomic identifications and, by extension, the definition of biogeographic provinces derived from them.
In this study we investigate spatial partitioning of intraspecific molecular polymorphisms in nominal Serolis paradoxa (Fabricius, 1775), a marine benthic shallow-water isopod, using two independent genetic marker systems. S. paradoxa was originally described from the Falkland Islands but is also frequently reported from the Strait of Magellan, the Patagonian coastal shelf (Atlantic and Pacific side), and also from regions further to the equator [5, 20]. For the current taxonomic status and synonyms of S. paradoxa see . The vertical distribution of S. paradoxa ranges from shallow waters (about 5 m, Held pers. observ., Lopaz-Gappa pers. observ.) down to 113 m . Although in the Magellan region S. paradoxa can be locally very abundant  it is often encountered at medium densities (about 1 ind/m2, Mutschke and Rios pers. comm.). Like almost all isopods, S. paradoxa lacks free-swimming distribution stages and broods its offspring in a ventral brood pouch, the marsupium, and is thus expected to be limited in its dispersal capabilities. No information on the life cycle and the duration of S. paradoxa are known. Based on significantly extended embryonic stages of Antarctic serolid isopods in comparison to non-Antarctic isopods  it can be assumed that embryonic development, maturation and brooding of S. paradoxa from the cool-temperate regions, each stage lasts several months. Altogether, S. paradoxa is expected to have very limited dispersal capacity due to these factors.
Direct measurements of dispersal and migration over large geographical distances provide a poor means of assessing effective gene flow. The small number of immigrants needed per generation to appreciably change the genetic structure of a population will not be picked up by realistic sampling schemes. Indirect genetic estimates use tools that interpret the genetic structure of a population as a result of past genetic influx and thus represent an easier and more reliable method . However, in this context historic extrinsic factors that may have exerted a structuring force must be considered when estimating present-day population structure. One extrinsic factor known to have had a major impact on genetic structure and distribution of species are glaciation events . Their influence on the marine fauna is two-fold: large-scale glaciations may directly render entire coastal habitats unavailable during glacial maxima [40–42] and also lead to a decrease in sea level of up to 130 m . The latter may disrupt inshore habitats on either side of an emerging barrier (e.g. appearance of the Panama land bridge , or connect shallow water habitats that are disjunct during periods of high sea level (additional file 1).
The focus of the present study is the genetic structure of S. paradoxa from the Falkland Islands and the Patagonian shelf. The present-day situation suggests that the deeper waters on the South American shelf may present an insurmountable barrier to S. paradoxa inhabiting shallow waters around Patagonia and the Falkland Islands. Historically, the lower sea level during glacial maxima may have connected both regions and facilitated migration between habitats thus counteracting independent lineage sorting in the two regions.
By investigating the coherence between gene pools and construction of an approximate timeline, we test whether the disruptive or unifying forces predominated and if the influence of the last glaciations exerted a major influence on the evolutionary history of S. paradoxa. We also test whether the major age difference between marine habitats in the central Strait of Magellan and around the Falkland Islands (< 14 KY BP [43–50] and millions of years , respectively) exerted a measurable influence on the genetics of populations living there today. In particular, we test for differences in genetic diversity and patterns of recent population expansions or secondary contact of colonizers from the Atlantic and the Pacific side in central Magellan Strait. We hypothesize that populations in the central Strait of Magellan are genetically less diverse than populations from the coast or the Falkland Islands due to recent range expansion into the Magellan Strait subsequent to the retreat of the glaciers after the last glacial maximum (LGM).
Materials and methods
Sampling sites (PA = Strait of Magellan near Punta Arenas, AO = Atlantic opening of the Strait of Magellan, FI = Falkland Islands).
Strait of Magellan near Punta Arenas
SCUBA diving (CH)
Opening of the Strait of Magellan to the Atlantic Ocean
2nd Cruce Estrecho
West Falkland Islands
DNA extraction, PCR, sequencing/genotyping
Total DNA was extracted from muscle tissue using the Qiagen DNeasy Mini Kit according to the standard tissue protocol. Only 100 μl of elution buffer were used to increase DNA concentration.
Microsatellite markers Spa04, Spa12, Spa34, Spa35, Spa39, Spa42 and Spa43  were applied to assess intraspecific genetic polymorphisms for all specimens from the three sampling sites. Standard 15 μl reactions consisted of 1× PCR HotMaster Buffer, 0.2 mM dNTPs, 0.5–0.75 μM of each primer (one labelled, one unlabelled), 0.03 U/μl HotMaster Taq (Eppendorf, 5-Prime), 0.5 M Betaine (Sigma Aldrich) and 5–20 ng of genomic DNA. Cycling conditions on an epgradient thermocycler (Eppendorf) were 2 min at 94°C followed by 30 to 34 cycles with 20 s at 94°C, 15 s at annealing temperature, 30 s at 65°C. A final extension step of 45 minutes at 65°C was performed to reduce in vitro artefacts due to incomplete adenylation of products [see  for details]. PCR products were visualized on 2% TBE agarose gels, diluted 1–15 fold with molgrade water (CARL ROTH) and 1 μl of the diluted product was denatured in a mixture of 14.7 μl HI-DI formamide with 0.3 μl GeneScan ROX 500 size standard (both Applied Biosystems). Allele length scoring was performed using the software GENEMAPPER 4.0 (Applied Biosystems). To minimize genotyping errors [52, 53], up to four independent reactions were performed on a subset of samples to estimate allele calling errors.
The universal primers 16Sar and 16Sbr  were used for amplification. Reactions were carried out in 25 μl volumes with 1× HotMaster reaction buffer, 0.2 mM dNTPs, 0.5 μM of each primer, 0.025 U/μl HotMaster Taq (Eppendorf, 5-Prime). Reaction conditions were: Initial denaturation for 2 min at 94°C followed by 36 cycles of 20 s at 94°, 15 s at 46°C and 80 s at 65°C plus a final elongation step of 5 min at 65°C. PCR products were purified using Qiagen QIAquick or Eppendorf Perfectprep Gel cleanup kits. Cycle-sequencing was performed in 10 μl reaction volumes using 1 μM of either 16Sar or 16Sbr primer, 1 μl of the purified template DNA and the BigDye Terminator Kit 3.1 chemistry (Applied Biosystems) according to the recommendations of the manufacturer. Reactions were purified according to the 'modified protocol' of the Qiagen DyeEx Kit. Sequencing was conducted on an ABI 3130xl sequencer.
Raw data were checked and corrected for genotyping errors using the software MICRO-CHECKER version 2.2.3  and DROPOUT version 1.3 . In addition, MICRO-CHECKER was used to test for the presence of null alleles in populations, i.e. alleles that fail to amplify due to substitutions in the primer binding regions. Corrected genotype tables were converted to specific software formats using the software MSTOOLKIT version 3.1  and CONVERT version 1.3.1 . The program ANIMALFARM version 1.0  was used to test for loci with significantly disproportionate variances that may bias allele-size based distance estimates such as Slatkin's or Rousset's RST estimates [60, 61]. Tests for Hardy Weinberg equilibrium (HWE) and linkage disequilibrium (LD) were performed using GENEPOP version 4.0.6 . Parameter settings: 10,000 dememorization steps, 50 batches, 20,000 MCMC sampling steps. HWE tests aim at testing whether there is a statistically significant deviation of genotype frequencies from those expected according to Mendelian inheritance. Linkage disequilibrium occurs when two genomic loci are not inherited independently, e.g. due to physical linkage or other processes at population level hindering independent recombination of loci.
To assess partitioning of genetic variability within individuals, subpopulations and regions, we performed hierarchical analyses of molecular variance (AMOVA) using ARLEQUIN version 3.11 . Therefore, populations PA and AO were assigned to one group, while FI constituted the other group. In addition, single and multilocus inbreeding coefficients (FIS) and pairwise population coancestry coefficients (FST, similar to Weir and Cockerham's Theta) were estimated as in  using GENEPOP. We also calculated pairwise allele-size based differentiation estimates, RST, according to  using GENEPOP. Significance was assessed by exact G tests as implemented in GENEPOP. The interpretation of the FST values from multiallelic data is problematic because their maximum values depend on the amount of within-population variation and even in the absence of any shared allele often fail to reach the theoretical maximum of 1 [65–67]. We therefore applied a standardization approach suggested by Hedrick  for calculations of GST and derived by Meirmans  for Analysis of Variance frameworks (ANOVA).
The main principle of this standardization approach is to correct the maximum possible value for FST as follows: FST(max) is calculated using GENEPOP applying the sampling bias correction suggested by Meirmans  using the Software RECODEDATA. F'ST was subsequently calculated by dividing FST by this inferred maximum value.
The standardized F'ST measure calculated range from 0 (populations equifrequent for all alleles) to 1 (populations fixed for different alleles) and therefore makes interpretation of the degree of subdivision much easier and facilitates comparing results among studies.
In addition to these ANOVA based coancestry estimates we performed individual assignment tests using the program STRUCTURE, version 2.2.2  to investigate population subdivision. The advantage of the Bayesian clustering algorithm of STRUCTURE is that no classification of populations has to be done a priori. Assuming HWE and no or only weak LD within subpopulations, STRUCTURE assigns individual genotypes probabilistically to populations and calculates the likelihood of the genotype dataset for a given number of populations (K), i.e. ln Pr (D|K) for K = 1 to K = n, using a Markov Chain Monte Carlo algorithm [69, 70]. For the S. paradoxa data set, the most likely number of populations was inferred without making assumptions on geographic origin of individuals. The number of MCMC steps needed to reach convergence was first estimated by comparing run lengths between 10,000 and 2,000,000 steps. Convergence was generally reached with <5,000 steps. Therefore, for the parameter sets 10 independent runs with a burn-in of 5,000 and subsequent 100,000 MCMC steps were performed with and without assuming recent admixture in the prior model, and considering alleles as correlated and uncorrelated. The number of clusters (K) to infer was defined from K = 1 to K = 4 to allow detection of potential cryptic subpopulations. Alpha was inferred from the data for each population separately. Results from 10 independent runs were analysed in CLUMPP, version 1.1.1  to compute a consensus membership coefficient Q-matrix from all 10 independent Q-matrices. Both the individual Q-matrix and averaged population Q-matrix were visualized using DISTRUCT, version 1.1 .
To assess estimates of the present effective population size (N e ), we applied the linkage disequilibrium method proposed by Hill (1981) , modified by Waples  to account for a bias correction when sample size is much smaller than effective population size. This method is implemented in the program LDNE, version 1.3 . Calculations of N e and the confidence intervals (CI) were estimated considering alleles with a frequency of c ≥ 0.05 and c ≥ 0.02 and ≥ 0.01, respectively.
Tests for historical population bottlenecks were performed using the program BOTTLENECK. Tests implemented in this program are based on the hypothesis that populations that have experienced recent reductions in their effective population size (N e ) show a reduction in both allelic richness and heterozygosity. In populations decreasing in size, the number of alleles (N A ) drops faster than heterozygosity  and therefore the observed heterozygosity is larger than the expected heterozygosity (H O > H E ). Conversely, in expanding populations often the number of alleles increases faster than heterozygosity until equilibrium is reached. From the comparison of both parameters, allelic diversity and heterozygosity, it is possible to make inferences regarding historical demography of a population. For each locus and population BOTTLENECK computes distribution of H E expected from the observed N A , given the sample size (n) under the assumption of mutation-drift equilibrium. This distribution is obtained through simulating the coalescent process of n genes under the three possible mutation models, i.e. a) the Infinite Allele Model (IAM), b) the Two-Phase Model (TPM), c) the Stepwise-Mutation Model (SMM). As recommended by Cornuet and Luikart  we tested several proportions of the SMM for the TPM (70–90%). Statistical significance of the parameters were inferred applying a Sign-test and a Wilcoxon-rank-test [76, 78, 79].
Assembly of forward and reverse strands and editing was performed using the software SEQMAN (Dnastar, Lasergene) and GENEIOUS version 4.0.2 (Biomatters Ltd.). Sequence alignment was performed using MUSCLE version 3.6 . The alignment required no manual correction based on secondary structure information . Sequence variation was analyzed using MEGA 4.0 . Gene diversity and nucleotide diversity according to Nei  and Theta based on the number of segregating sites, Theta (S), were calculated with ARLEQUIN version 3.11. Genetic differentiation between populations and between regions [(PA + AO) vs. FI] were assessed using an FST and AMOVA framework as implemented in ARLEQUIN. Assuming neutrality, evidence of a population expansion was tested applying Tajima's D  and Fu's F S statistic  as implemented in ARLEQUIN applying a coalescent simulation approach generating 10,000 selectively neutral samples for assessment of significance of results. A test for sudden population expansion based on the pairwise mismatch distribution was calculated using ARLEQUIN and significance was assessed by 50,000 pseudo replicates.
A statistical parsimony network with a 95% connection-probability limit was created for the 490-bp alignment using TCS version 1.21 . In addition, two outgroup sequences of the serolid isopods Cuspidoserolis luethjei and Cuspidoserolis johnstoni (GenBank accession numbers AJ269802, AJ269803; see ) were aligned to the S. paradoxa sequences using MUSCLE, resulting in a 492-bp alignment. This alignment was used to calculate a neighbor joining tree  with bootstrap support (1000 replicates) based on uncorrected p-distances using PAUP* version 4b10 .
The coalescent-based MCMC approach implemented in the software BEAST was used to date the splitting event between the different Serolis paradoxa lineages applying both a strict and a relaxed molecular clock model. The sequence model HKY85 was used for modelling sequence evolution  together with a predefined mutation rate of 0.37% per million years , based on a molecular clock for the serolid isopod Ceratoserolis trilobitoides. Dating times and confidence intervals (CI) were filtered using TRACER version 1.4 .
Total number of specimens scored for each locus (N S ), number of different alleles (N A ), inbreeding coefficient (F IS ), observed heterozygosity (H O ) and expected heterozygosity (H E ) for the seven microsatellites and three populations of Serolis paradoxa.
Mean N A /locus
Mean N A per location
Mean H O
Mean H E
Results of ANIMALFARM confirmed that none of the loci contributed disproportionally to distance-based differentiation estimates after Bonferroni or Sidak adjustment of the significance level.
Hierarchical analyses of molecular variance (AMOVA) among Serolis paradoxa populations within and between two regions using 7 microsatellite markers.
Component of differentiation
ΦCT = 0.310
Among populations within regions
ΦSC = -0.024
Among individuals within populations
ΦIS = 0.119
ΦIT = 0.407
Genetic differentiation among populations of Serolis paradoxa from three stations as assessed by F-statistics (FST, lower diagonal) and R-statistics (RST, upper diagonal), based on seven polymorphic microsatellite loci.
Standardized pairwise FST calculates  in this study showed very strong pairwise population differentiation between Patagonia and the Falklands (PA vs FI: 0.86; AO vs FI: 0.91), and much smaller values among the Magellan Strait populations (PA vs AO: 0.063). These values demonstrate that both regions are almost fixed for different alleles at the seven loci investigated. When removing locus Spa39, which is biased due to the presence of null-alleles in population FI, the standardized values did not change, however, the non-standardized FST values were almost twice as high (data not shown).
Estimating the present effective population sizes using the LD approach , we consistently received negative Ne estimates with confidence intervals ranging from negative values to infinity, thus indicating very large population sizes . When testing for recent demographic contractions or expansions by looking for deviations from mutation-drift equilibrium under different mutation models using BOTTLENECK we found a significant heterozygosity deficiency under particular mutation models. For population AO there was a significant heterozygosity deficit under both SMM and TPM models (additional file 3), which provides strong evidence for recent population expansion. For FI the evidence for recent population expansion was somewhat weaker: a significant heterozygosity deficit was detected only using the SMM and the Wilcoxon test (P = 0.0195, additional file 3). Thus results of BOTTLENECK do not provide evidence for a similarly drastic decline and subsequent recovery in population size for PA. For population AO under a TPM and a strict SMM, the significant excess of heterozygosity may indicate that this population is expanding presently. Although the evolutionary dynamics of microsatellites are not fully understood [94, 95] it is commonly accepted that the IAM model is not an appropriate descriptor of the mutational dynamics of microsatellite markers and hence that its application often leads to unrealistic conclusions.
In summary, the results from the microsatellite analyses provide evidence for moderate differentiation between the two Patagonian populations and very strong subdivision between populations from Patagonia and the Falkland Islands. Genetic diversity was highest in the center of the Strait of Magellan, lower near its opening towards the Atlantic Ocean and lowest around the Falkland Islands. All populations showed a significant heterozygosity deficit corroborated by high FIS values (Table 2) which may be indicative for inbreeding of local populations.
Distribution of the n = 71 16S rDNA sequences on the two sampling locations and GenBank accession numbers.
Hierarchical analyses of molecular variance (AMOVA) among Serolis paradoxa populations within and between two regions based on the 16S rDNA data.
Component of differentiation
ΦCT = 0.905
Among populations within regions
ΦSC = -0.494
ΦST = 0.952
Genetic diversity and neutrality indices for the 16S rDNA data sets.
Magellan Strait (PA)
Atlantic Opening (AO)
Falkland Islands (FI)
0.0014 ± 0.0012
0.0006 ± 0.0008
0.5755 ± 0.0952
0.2600 ± 0.1202
1.038 ± 0.5900
0.823 ± 0.5237
-0.945 (P = 0.200)
-1.240 (P = 0.085)
-1.729 (P = 0.009)
-1.819 (P = 0.059)
-1.827 (P = 0.034)
-2.889 (P = 0.001)
Dating of the time to the most recent common ancestor (tMRCA) between both S. paradoxa lineages using a strict molecular clock for the mutation rate differed for the two monophyletic lineages. For the Patagonian taxa as ingroup the tMRCA inferred was 0.948 MY (5% CI 0.344 MY, 95% CI 1.658 MY) and for the Falkland Islands lineage 0.643 MY (5% CI 0.136 MY, 95% CI 1.207 MY). Thus, from both inferences, evidence for a splitting event in the mid-Pleistocene is supported.
The genetic variability within the nominal species Serolis paradoxa turned out to have extensive spatial structure. The differences in mutation rates and coalescent dynamics of the two marker systems help describe present-day population structure and reconstruct historical demographic processes.
Two genetically distinct lineages
There is strong evidence of divergence between populations from Patagonia and the Falkland Islands, supported by microsatellite and mitochondrial data. The dominant feature of the intraspecific variability of mitochondrial DNA data for the Patagonian populations (PA and AO) and the Falkland Island population is that populations form two shallow subnetworks, corresponding to the two geographic regions. The nuclear microsatellite markers support the geographic partitioning of variation with high and significant FST, RST differentiation values and strong support from Bayesian cluster analyses (Figure 2).
The geographic positions of our sampling locations along an East-West axis might suggest testing for isolation by distance effects (IBD). However, in the context of this study, the IBD is an inappropriate method and it is unlikely that this would become more meaningful even if more intermediate sampling locations were available. This is because the central Strait of Magellan became available for (re)colonization only very recently, approximately 9–14 KY BP [46, 48]. This rapid range expansion is typically accompanied by loss of alleles and an excess of homozygosity  which violates a mutation-drift equilibrium assumed by the IBD model. Investigating distance effects on the distribution of intraspecific variance inside the Magellan Strait offers a means to trace the recolonization of this young habitat and would be appropriate for IBD but this requires more fine-scaled sampling and is outside the scope of this paper.
Absence of effective gene flow between the Falkland Islands and Patagonia is strongly suggested by nearly fixed population specific differences in fast evolving microsatellites and the perfect congruence of haplotype identity and geography for the 16S rDNA data. The long branch connecting the two groups of haplotypes (Figure 3) and their reciprocal monophyly (Figure 4) indicates complete lineage sorting in both groups. The magnitude of genetic differentiation between 16S genotypes is on the order of magnitude known for reproductively isolated species [26–28, 97]. Speciation ultimately involves the irreversible disruption of a once contiguous gene pool into two . The recognition of species thus centers around direct or indirect evidence for gene flow between them. Our data from two independent molecular markers are in line with the expectations of two independently evolving lineages. The patterns and magnitude of the remaining differences do not suggest the presence of additional cryptic species inside (PA and AO) vs. (FI), they indicate, however, that gene flow is restricted even within the Strait of Magellan. The congruence between both marker systems supports that the 16S rRNA gene tree reflects the species tree rather than being a result of shared ancestral polymorphisms  or other processes affecting mitochondrial genes (see  for review).
Evolutionary history of nominal Serolis paradoxa
Southern Hemisphere glaciations differently affected both regions: The Falkland Islands were little affected by glacial advances , thus S. paradoxa was able to survive by following the rising and falling sea levels. In Patagonia, however, major parts of today's distribution of nominal S. paradoxa became unavailable due to ice coverage and/or low sea levels. Western Patagonia was covered by a contiguous ice shield similar to the Antarctic Peninsula today and the central Strait of Magellan was inundated only after the LGM, approximately 14-9 KY BP . Contrary to the situation around the Falkland Islands where S. paradoxa was presumably continuously present over evolutionary times, this species was forced to immigrate into the Strait of Magellan only recently after the retreat of the glaciers. Surprisingly, genetic diversity estimates for population PA from central Magellan Strait (Table 7) indicate that the population has the highest genetic diversity and shows almost no signs for recent population expansion (Table 7, additional file 3, 4), although colonization of a new habitat is often accompanied by a loss of genetic diversity (founder event). In comparison, population FI is less diverse for the 16S rDNA with one dominant haplotype only (HT7) and reveals strong evidence for recent population expansion. Diversity estimates decline from population PA in the west to population FI in the east, which seems counterintuitive as the effects of past glaciations are likely to have been much more severe for PA than for AO and FI. This apparent contradiction may, however, be explained by the fact that the Magellan Strait was recolonized after the LGM not only from the Atlantic but also the Pacific side, thus receiving allelic diversity from different source populations. In the contact zone in the central Magellan Strait, this scenario explains the inflated genetic diversity estimates for PA.
In summary, our data are in agreement with the following scenario: Populations of an ancestral species were separated geographically and evolved in allopatry (Falkland Islands vs. Pacific and Atlantic side of Patagonia). Applying a rate for the accumulation of substitutions in 16S rDNA estimated by Held  for the serolid isopod Ceratoserolis trilobitoides (Eights, 1833) with a rate of 0.37% per MY for transitions and transversions, the time of divergence was estimated to have occurred several hundreds of thousands of years before present. Thus the initial separation of lineages predates the last glaciations and took place in mid-Pleistocene (average estimates 0.643 – 0.948 BP). In view of the strong genomic signatures of differentiation between Patagonia and the Falkland Islands we must therefore reject the hypothesis that low sea levels during glacial periods led to significantly elevated levels of gene flow between populations of S. paradoxa due to greater proximity of shallow-water habitats. A similar argument applies to potential migration between Patagonia and the Falkland Islands via passive rafting on drifting substrates. Although there are major directional ocean currents that frequently transport substrates suitable for transportation of even rather immobile species [101–103] this apparently played no role in the recent evolutionary history of S. paradoxa. This species exclusively inhabits soft-bottom shallow waters and is frequently half-buried in the sediment (Held pers. observ.). Its capability to colonize new island habitats and maintain genetic continuity across barriers to dispersal and over evolutionary times is therefore small. Further sampling effort should focus on sampling specimens from the West Falkland Islands. It cannot be excluded that members of both lineages live in sympatry today.
Reliability and systematic bias in differentiation estimates
The equilibrium FST estimate for totally isolated populations based on microsatellites can reach the maximum value of FST = 1 only theoretically. Due to the high mutation rate of microsatellites [94, 104, 105] and often a restricted allelic spectrum (, but see ), the intrapopulation variability is generally very high in particular after a long time of independent evolution of large populations. Applying Meirman's standardization approach for pairwise FST, differentiation between PA and FI is 0.86, between AO and FI 0.91 and among the Magellan Strait populations PA and AO 0.063 and thus about three times larger than without this correction. These values underpin that populations from Patagonia and the Falkland Islands are almost fixed for different alleles at the seven loci investigated. The results point out the importance of the recently introduced standardization approach [67, 68] in order to allow for easier comparison and interpretation of the data. Differentiation estimates of the 16S rDNA yield comparable results. Differentiation was significant between all three populations. Although populations AO and PA shared the most common haplotypes, FST estimates between PA and AO revealed much higher differentiation than inferred using microsatellite data. The most plausible explanation is that the fourfold smaller effective population size of mitochondrial DNA  lead to much stronger effects of genetic drift, resulting in higher differentiation estimates.
In principle, differentiation estimates can also be biased due to comparing samples obtained in different years (PA: 1997, AO: 2003, FI: 2004). However, as only few years, corresponding to even fewer generations of S. paradoxa, separate the samples and no major disturbances in the regions were reported for the time in-between the samplings. Thus, we regard this a negligible issue.
Concerning the dating of the split between the two lineages, it must be stated that genetic distances between two lineages increase much faster than predicted by molecular clocks if populations experience population bottlenecks . Thus, the realistic tMRCA between the two lineages might be shorter than the estimated mean using the molecular clock. In addition, it is not entirely certain if the molecular clock can be applied to S. paradoxa. The time estimates are based on 16S rRNA substitution rates commonly used for other Crustacea .
Taxonomic and conservation status of the newly delimited species
The genetic data strongly suggest that nominal Serolis paradoxa (Fabricius, 1775) consists of two reproductively isolated species one of which occurs in Patagonia while the other is presumably confined to shallow waters around the Falkland Islands. As the type was originally described by Fabricius as Oniscus paradoxum Fabricius, 1775 from the Falkland Island the species from Patagonia is in need of formal description and a scientific name.
The occurrence of cryptic species has important implications for the conservation of biodiversity in general . If a cryptic species is not recognized, unique and endangered local faunas cannot be efficiently protected. However, the estimates of effective population size for both species contained inside nominal Serolis paradoxa imply that both are highly abundant and neither needs to be considered endangered.
In summary, our data prove low but significant differentiation among populations within the Strait of Magellan and the absence of effective gene flow among populations from the Strait of Magellan and the Falkland Islands. In fact, specimens from both regions belong to two cryptic lineages that probably diverged in the mid-Pleistocene and may already represent reproductively isolated species. The 16S rDNA data supports a genetically rich central Strait of Magellan population, an intermediate population near the Atlantic opening of the Strait of Magellan and a genetically depauperate Falkland Island population. The results are in line with the expectations of colonization of the central Strait of Magellan from both sides of Patagonia after the last glacial maximum approximately 9-14 KY after deglaciation of the habitat and rise of sea levels.
While the fauna of the Falkland Islands has often been accepted to share most of their faunal inventory with Patagonia our results indicate that shallow water species with low mobility may in fact be strongly differentiated populations of one species or even reproductively isolated species.
We thank Erika Mutschke and Carlos Rios (Universidad de Magallanes, Punta Arenas) for providing us with material from the 2nd Cruce estrecho and helpful information, 2003. We also thank Andrea Eschbach and Shobhit Agrawal (AWI Bremerhaven) for technical assistance in this study. This work was supported by a DFG grant HE-3391/3 to CH, NSF grant OPP-0132032 to H.W. Detrich III, and a DAAD scholarship to FL and AK. This is publication number 22 from the ICEFISH Cruise of 2004.
- Avise JC: Phylogeography. The history and formation of species. 2000, Cambridge: Harvard University Press
- Brandt A: Origin of Antarctic Isopoda (Crustacea, Malacostraca). Marine Biology. 1992, 113: 415-423.
- Wägele J: Notes on Antarctic and South American Serolidae (Crustacea, Isopoda) with remarks on the phylogenetic biogeography and a description of new genera. Zoologische Jahrbücher Abteilung für Systematik Ökologie und Geographie der Tiere. 1994, 121: 3-69.
- Knox GA, Lowry JK: A comparison between the benthos of the Southern Ocean and the North Polar Ocean with special reference of the Amphipoda and the Polychaeta. Polar Oceans. Edited by: Dunbar MJ. 1977, Calgary: Arctic Institute of North America, 423-462.
- Lancellotti DA, Vasquez JA: Zoogeografía de macroinvertebrados bentónicos de la costa de Chile: contribución para la conservación marina. Revista Chilena de Historia Natural. 2002, 73: 99-129.
- McDowall RM: Falkland Islands biogeography: converging trajectories in the Southern Ocean. Journal of Biogeography. 2005, 32: 49-62.
- Longhurst AR: Ecological geography of the sea. 2007, Burlington, MA, USA: Academic Press
- Montiel San Martin A: Biodiversity, zoogeography and ecology of polychaetes from the Magellan region and adjacent areas. Berichte zur Polarforschung. 2005, 505: 1-112.
- Adie RJ: The position of the Falkland Islands in a reconstruction of Gondwanaland. Geological Magazine. 1952, 89: 401-410.
- Mitchell C, Taylor GK, Cox KG, Shaw J: Are the Falkland Islands a rotated microplate?. Nature. 1986, 319: 131-134.
- Marshall JEA: The Falkland Islands: A key element in Gondwana palaeogeography. Tectonics. 1994, 13: 499-514.
- Knowlton N: Sibling species in the sea. Annual Review of Ecology and Systematics. 1993, 24: 189-216.
- Palumbi SR: Genetic divergence, reproductive isolation and marine speciation. Annual Review of Ecology and Systematics. 1994, 25: 547-572.
- Bohonak AJ: Dispersal, gene flow, and population structure. Quarterly Review of Biology. 1999, 74: 21-43.PubMed
- Knowlton N: Molecular genetic analyses of species boundaries in the sea. Hydrobiologia. 2000, 420: 73-90.
- Gallardo MH, Penaloza L, Clasing E: Gene flow and allozymic population structure in the clam Venus antiqua (King of Broderip), (Bivalvia, Veneriidae) from Southern Chile. Journal of Experimental Marine Biology and Ecology. 1998, 230: 193-205.
- Hunter RL, Halanych KM: Evaluating connectivity in the brooding brittle star Astrotoma agassizii across the drake passage in the Southern Ocean. Journal of Heredity. 2008, 99: 137-48.PubMed
- Pérez-Barros P, D'Amato M, Guzmán NV, Lovrich GA: Taxonomic status of two South American sympatric squat lobsters, Munida gregaria and Munida subrugosa (Crustacea: Decapoda: Galatheidae), challenged by DNA sequence information. Biological Journal of the Linnean Society. 2008, 94: 421-434.
- Thornhill DJ, Mahon AR, Norenburg JL, Halanych KM: Open-ocean barriers to dispersal: a test case with the Antarctic Polar Front and the ribbon worm Parborlasia corrugatus (Nemertea: Lineidae). Molecular Ecology. 2008, 17: 5104-5117.PubMed
- Brandt A: Colonization of the Antarctic shelf by the Isopoda (Crustacea, Malacostraca). Berichte zur Polarforschung. 1991, 98: 1-240.
- Brandt A, Linse K, Mühlenhardt-Siegel U: Biogeography of Crustacea and Mollusca of the Subantarctic and Antarctic regions. Scientia Marina. 1999, 63: 383-389.
- Linse K: New records of shelled marine molluscs at Bouvet Island and preliminary assessment of their biogeographic affinities. Polar Biology. 2006, 29: 120-127.
- Allcock AL, Brierley AS, Thorpe JP, Rodhouse PG: Restricted gene flow and evolutionary divergence between geographically separated populations of the Antarctic octopus Paraledone turqueti. Marine Biology. 1997, 129: 97-102.
- Held C: Phylogeny and biogeography of serolid isopods (Crustacea, Isopoda, Serolidae) and the use of ribosomal expansion segments in molecular systematics. Mol Phylogenet Evol. 2000, 15 (2): 165-178.PubMed
- Page TJ, Linse K: More evidence of speciation and dispersal across the Antarctic Polar Front through molecular systematics of Southern Ocean Limatula (Bivalvia: Limidae). Polar Biology. 2002, 25: 818-826.
- Held C: Molecular evidence for cryptic speciation within the widespread Antarctic crustacean Ceratoserolis trilobitoides (Crustacea, Isopoda). Antarctic biology in a global context. Edited by: Huiskes AH, Gieskes WW, Rozema J, Schorno RM, van der Vies SM, Wolff WJ. 2003, Leiden: Backhuys Publishers, 135-139.
- Held C, Wägele J-W: Cryptic speciation in the giant Antarctic isopod Glyptonotus antarcticus (Isopoda: Valvifera: Chaetiliidae). Scientia Marina. 2005, 69: 175-181.
- Raupach MJ, Wägele JW: Distinguishing cryptic species in Antarctic Asellota (Crustacea: Isopoda) – a preliminary study of mitochondrial DNA in Acanthaspidia drygalskii. Antarctic Science. 2006, 18: 191-198.
- Linse K, Cope T, Lörz A-N, Sands C: Is the Scotia Sea a centre of Antarctic marine diversification? Some evidence of cryptic speciation in the circum-Antarctic bivalve Lissarca notorcadensis (Arcoidea: Philobryidae). Polar Biology. 2007, 30: 1059-1068.
- Wares JP, Daley S, Wetzer R, Toonen RJ: An evaluation of cryptic lineages of Idotea balthica (Isopoda: Idoteidae): Morphology and microsatellites. Journal of Crustacean Biology. 2007, 27: 643-648.
- Wilson NG, Hunter RL, Lockhart SJ, Halanych KM: Multiple lineages and absence of panmixia in the "circumpolar" crinoid Promachocrinus kerguelensis from the Atlantic sector of Antarctica. Marine Biology. 2007, 152: 895-904.
- Leese F, Kop A, Agrawal S, Held C: Isolation and characterization of microsatellite markers from the marine isopods Serolis paradoxa and Septemserolis septemcarinata (Crustacea: Peracarida). Molecular Ecology Resources. 2008, 8: 818-821.PubMed
- Mahon AR, Arango CP, Halanych KM: Genetic diversity of Nymphon (Arthropoda: Pycnogonida: Nymphonidae) along the Antarctic Peninsula with a focuson Nymphon australe Hodgson 1902. Marine Biology. 2008, 155: 315-323.
- Bastida R, Torti MR: Los isópodos Serolidae de la Argentina. Clave para su reconocimiento. Physis Sección A. 1973, 32: 19-46.
- Gappa L, Sueiro M: The subtidal macrobenthic assemblages of Bahia San Sebastian (Tierra del Fuego, Argentina). Polar Biology. 2007, 30: 679-687.
- Rios C, Mutschke E, Morrison E: Biodiversidad bentónica sublitoral en el estrecho de Magallanes, Chile. Revista de Biologia Marina y Oceanografia. 2003, 38: 1-12.
- Wägele J: On the reproductive biology of Ceratoserolis trilobitoides (Crustacea: Isopoda): Latitudinal variation of fecundity and embryonic development. Polar Biology. 1987, 7: 11-24.
- Neigel JE: Population genetics and demography of marine species. Marine biodiversity: Patterns and processes. Edited by: Ormond RF, Gage JD, Angel MV. 1997, Cambridge: Cambridge University Press, 274-292.
- Hewitt GM: The genetic legacy of the quaternary ice ages. Nature. 2000, 405: 907-913.PubMed
- Anderson JB, Shipp SS, Lowe AL, Weller JS, Mosola AB: The Antarctic ice sheet during the last glacial maximum and its subsequent retreat history: a review. Quaternary Science Reviews. 2002, 21: 49-70.
- Huybrechts P: Sea-level changes at the LGM from ice-dynamic reconstructions of the Greenland and Antarctic ice sheets during the glacial cycles. Quaternary Science Reviews. 2002, 21: 203-231.
- Ingolfsson O: Quaternary glacial and climate history of Antarctica. Quaternary Glaciations – Extent and Chronology, Part III. Edited by: Ehlers J, Gibbard PL. 2004, Elsevier, 3-43.
- Polanski J: The maximum glaciation in the Argentine Cordillera. Geological Society of America, Inc. Special Paper. 1965, 84: 453-472.
- Flint RF, Fidalgo F: Glacial drift in the eastern Argentine Andes between latitude 41° 10'S and latitude 43° 10'S. Geological Society of America Bulletin. 1969, 80: 1043-1052.
- Hulton N, Sugden D, Payne A, Clapperton C: Glacier modeling and the climate of Patagonia during the last glacial maximum. Quaternary Research. 1994, 42: 1-19.
- Clapperton CM, Sugden DE, Kaufman DS, McCulloch RD: The last glaciation in central Magellan Strait, southernost Chile. Quaternary Research. 1995, 44: 133-148.
- Benn DI, Clapperton CM: Glacial Sediment-Landform Associations and Paleoclimate during the Last Glaciation, Strait of Magellan, Chile. Quaternary Research. 2000, 54: 13-23.
- McCulloch RD, Bentley MJ, Purves RS, Hulton NRJ, Sugden DE, Clapperton CM: Climatic inferences from glacial and palaeoecological evidence at the last glacial termination, southern South America. Journal of Quaternary Sciences. 2000, 15: 409-417.
- Rostami K, Peltier WR, Mangini A: Quaternary marine terraces, sea-level changes and uplift history of Patagonia, Argentina: comparisons with predictions of the ICE-4G (VM2) model of the global process of glacial isostatic adjustment. Quaternary Science Reviews. 2000, 19: 1495-1525.
- Hulton NRJ, Purves RS, McCulloch RD, Sugden DE, Bentley MJ: The Last Glacial Maximum and deglaciation in southern South America. Quaternary Science Reviews. 2002, 21: 233-241.
- Clapperton CM: Quaternary glaciations in the Southern Ocean and the Antarctic Peninsula area. Quaterary Science Reviews. 1990, 9: 229-252.
- Pompanon F, Bonin A, Bellemain E, Taberlet P: Genotyping errors: causes, consequences and solutions. Nature Reviews Genetics. 2005, 6: 847-859.PubMed
- Hoffman JI, Amos W: Microsatellite genotyping errors: detection approaches, common sources and consequences for paternal exclusion. Molecular Ecology. 2005, 14: 599-612.PubMed
- Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, Weighting, and Phylogenetic Utility of Mitochondrial Gene-Sequences and A Compilation of Conserved Polymerase Chain-Reaction Primers. Annals of the Entomological Society of America. 1994, 87: 651-701.
- van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Notes. 2004, 4: 535-538.
- McKelvey KS, Schwartz MS: DROPOUT: A program to identify problem loci and samples for noninvasive genetic samples in a capture-mark-recapture framework. Molecular Ecology Notes. 2005, 5: 716-718.
- Park SDE: Trypanotolerance in west African cattle and the population genetic effects of selection. 2001, University of Dublin
- Glaubitz JC: CONVERT: A user-friendly program to reformat diploid genotypc data for commonly used population genetic software packages. Molecular Ecology Notes. 2004, 4: 309-310.
- Landry PA, Koskinen MT, Primmer CR: Deriving evolutionary relationships among populations using microsatellites and (deltamu)(2): all loci are equal, but some are more equal than others. Genetics. 2002, 161: 1339-47.PubMed CentralPubMed
- Slatkin M: A measure of population subdivision based on microsatellite allele frequencies. Genetics. 1995, 139: 457-462.PubMed CentralPubMed
- Rousset F: Equilibrium values of measures of population subdivision for stepwise mutation processes. Genetics. 1996, 142: 1357-1362.PubMed CentralPubMed
- Rousset F: GENEPOP '007: a complete reimplementation of the GENEPOP software for Windows and Linux. Molecular Ecology Notes. 2008, 8: 103-106.
- Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online. 2005, 1: 47-50.PubMed Central
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370.
- Hedrick P: Highly variable loci and their interpretation in evolution and conservation. Evolution. 1999, 53: 313-318.
- Whitlock MC, McCauley DE: Indirect measures of gene flow and migration: FST not equal to 1/(4Nm + 1). Heredity. 1999, 82 (Pt 2): 117-25.PubMed
- Hedrick P: A standardized genetic differentiation measure. Evolution. 2005, 59: 1633-1638.PubMed
- Meirmans PG: Using the AMOVA framework to estimate a standardized genetic differentiation measure. Evolution. 2006, 60: 2399-402.PubMed
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMed
- Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-1587.PubMed CentralPubMed
- Jakobsson M, Rosenberg NA: CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007, 23: 1801-1806.PubMed
- Rosenberg NA: Distruct: a program for the graphical display of population structure. Molecular Ecology Notes. 2004, 4: 137-138.
- Hill WG: Estimation of effective population size from data on linkage disequilibrium. Genetical Research. 1981, 38: 209-216.
- Waples RS: A bias correction for estimates of effective population size based on linkage disequilibrium at unlinked gene loci. Conservation Genetics. 2006, 7: 167-184.
- Waples RS, Do C: LDNe: a program for estimating effective population size from data on linkage disequilibrium. Molecular Ecology Resources. 2008, 8: 753-756.PubMed
- Piry S, Luikart G, Cornuet JM: BOTTLENECK: A computer program for detecting recent reductions in the effective population size using allele frequency data. Journal of Heredity. 1999, 90: 502-503.
- Nei M, Maruyama T, Chakraborty R: The bottleneck effect and genetic variability in populations. Evolution. 1975, 29: 1-10.
- Cornuet JM, Luikart G: Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics. 1996, 144: 2001-14.PubMed CentralPubMed
- Luikart G, Allendorf FW, Cornuet JM, Sherwin WB: Distortion of allele frequency distributions provides a test for recent population bottlenecks. Journal of Heredity. 1998, 89: 238-247.PubMed
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004, 32: 1792-1797.PubMed CentralPubMed
- Gutell RR: Collection of small subunit (16S- and 16S-like) ribosomal RNA structures. Nucleic Acids Research. 1993, 21: 3051-3054.PubMed CentralPubMed
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599.PubMed
- Nei M: Molecular evolutionary genetics. 1987, New York: Columbia University Press
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMed CentralPubMed
- Fu Y-X: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-923.PubMed CentralPubMed
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Molecular Ecology. 2000, 9: 1657-1659.PubMed
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMed
- Swofford DL: PAUP*. Phylogenetic analysis using parsimony (* and other methods). 1998
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology. 2007, 7: 214-PubMed CentralPubMed
- Hasegawa M, Kishino H, Yano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. Journal of Molecular Evolution. 1985, 22: 160-174.PubMed
- Held C: No evidence for slow-down of molecular substitution rates at subzero temperatures in Antarctic serolid isopods (Crustacea, Isopoda, Serolidae). Polar Biology. 2001, 24: 497-501.
- Rambaut A, Drummond AJ: Tracer v1.4. 2007, [http://beast.bio.ed.ac.uk/Tracer]
- Rice WR: Analyzing tables of statistical tests. Evolution. 1989, 43: 223-225.
- Ellegren H: Microsatellites: Simple sequences with complex evolution. Nature Reviews Genetics. 2004, 5: 435-445.PubMed
- Buschiazzo E, Gemmell NJ: The rise, fall and renaissance of microsatellites in eukaryotic genomes. Bioessays. 2006, 28: 1040-50.PubMed
- Ibrahim KM, Nichols RA, Hewitt GM: Spatial patterns of genetic variation generated by different forms of dispersal during range expansion. Heredity. 1996, 77: 282-291.
- Schubart CD, Diesel R, Hedges SB: Rapid evolution to terrestrial life in Jamaican crabs. Nature. 1998, 393: 363-365.
- Coyne JA, Orr HA: Speciation. 2004, Sinauer Associates, Inc
- Moore WS: Inferring phylogenies from mtDNA variation: Mitochondrial-gene trees versus nuclear-gene trees. Evolution. 1995, 49: 718-726.
- Ballard JW, Whitlock MC: The incomplete natural history of mitochondria. Molecular Ecology. 2004, 13: 729-744.PubMed
- Helmuth B, Veit RR, Holberton R: Long-distance dispersal of a subantarctic brooding bivalve (Gaimardia trapesina) by kelp-rafting. Marine Biology. 1994, 120: 421-426.
- Thiel M, Gutow L: The ecology of rafting in the marine environment – I – The floating substrata. Oceanography And Marine Biology: An Annual Review. 2005, 42: 181-263.
- Leese FL, Agrawal SA, Held CH: An exception to the rule? Long-distance dispersal of the brooding benthic isopods Septemserolis septemcarinata from remote Southen Ocean islands. Journal of Biogeography. 2008, under review
- Weber JL, Wong C: Mutation of human short tandem repeats. Human Molecular Genetics. 1993, 2: 1123-1128.PubMed
- Schug MD, Hutter CM, Wetterstrand KA, Gaudette MS, Mackay TFC, Aquadro CF: The mutation rates of di-, tri- and tetranucleotide repeats in Drosophila melanogaster. Molecular Biology and Evolution. 1998, 15: 1751-1760.PubMed
- Nauta MJ, Weissing FJ: Constraints on allele size at microsatellite loci: implications for genetic differentiation. Genetics. 1996, 143: 1021-1032.PubMed CentralPubMed
- Primmer CR, Ellegren H, Saino N, Moller AP: Directional evolution in germline microsatellite mutations. Nature Genetics. 1996, 13: 391-393.PubMed
- Birky CW, Maruyama T, Fuerst P: An approach to population and evolutionary genetic theory for genes in mitochondria and chloroplasts and some results. Genetics. 1983, 103: 513-527.PubMed CentralPubMed
- Bickford D, Lohman DJ, Sohdi NS, Ng PK, Meier R, Winker K, Ingram KK, Das I: Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 2007, 22 (3): 148-155.PubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.