BACKGROUND
phylogenetic relationships among cetacean taxa are contended at many different levels.
however, robust phylogenies are necessary for gaining insight into the evolutionary histories of these taxa and can help in understanding speciation of highly mobile taxa in an environment with seemingly few barriers to movement.
studies using molecular markers, often mitochondrial dna  sequences, to elucidate phylogenetic relationships among cetacean taxa have just as often created new controversy as resolved standing controversy  <cit> .
at deeper evolutionary levels, nuclear molecular data supported early conclusions based on morphological data that the hippopotamids are the sister lineage to cetaceans within the artiodactyla  <cit> , although, more recent morphological analyses do no support this relationship and instead identify raoellids as the sister group to cetaceans  <cit> .
within the cetacea, a previously controversial issue was the placement of sperm whales ; while initial mtdna sequence analyses suggested the sperm whale lineage is more closely related to the baleen whales rather than the rest of the toothed whales  <cit> , analysis of mtdna cytochrome b data  <cit> , multiple mtdna sequences and morphology  <cit>  and nuclear data  <cit>  all support the traditional view of monophyletic suborders.
for some cetacean lineages, such as the beaked whales, mitochondrial sequences have proven excellent markers for resolving phylogenies and even pinpointing new species  <cit> , although nuclear markers offer new evolutionary insights  <cit> .
for the endangered right whale species, mitochondrial markers render phylogenies congruent with those inferred from nuclear data  <cit> .
one group where relationships among taxa still remain unresolved is the family delphinidae, and in particular the subfamily delphininae.
this subfamily encompasses eleven nominal species : tursiops truncatus , t. aduncus , stenella frontalis , s. attenuata , s. longirostris , s. clymene , s. coeruleoalba , delphinus delphis linnaeus,  <dig>  d. capensis gray,  <dig>  lagenodelphis hosei fraser,  <dig>  and sousa chinensis .
morphological data have most often excluded the genus sotalia gray,  <dig> from this subfamily  in concordance with the mtdna cytochrome b phylogeny of leduc et al.
<cit> , but a recent molecular analysis suggested the genus should be included  <cit> .
however, interpreted in the most conservative manner, the combined analysis of mtdna and nuclear loci in caballero et al.
<cit>  supports a sister relationship between the lineages leading to the subfamily delphininae and the genus sotalia but is not necessarily strong evidence for inclusion of sotalia in the subfamily; a subtle but important distinction.
a number of morphological and genetic studies have been conducted in an attempt to resolve evolutionary relationships among taxa within the family delphinidae and/or subfamily delphininae  <cit> .
leduc et al.
<cit>  and may-collado et al.
<cit>  used mtdna cytochrome b sequences to reconstruct relationships among the taxa within the family delphinidae.
although the leduc et al.
<cit>  phylogeny included representatives of all species within the subfamily delphininae, each species was represented by only a few sequences ; the mitochondrial cytochrome b locus was unable to completely resolve branching order within the delphininae  <cit> .
the lack of resolution in the cytochrome b phylogeny suggests the delphininae are the product of a recent radiation; divergence among the numerous taxa is small, rendering resolution of branching order difficult.
the leduc et al.
<cit>  cytochrome b phylogeny also suggests polyphyly of the delphinine genera stenella and tursiops and this result for stenella is also supported by may-collado et al.
<cit> .
so far, nuclear and combined analyses concerning the delphininae raise more questions and provide support for only a few internal nodes within the subfamily  <cit> .
the polyphyletic genera and lack of branching order resolution among many of the delphinine taxa point to a need for a new approach to discerning the evolutionary relationships among the species in this subfamily and perhaps a revision of the subfamily .
we attempted to resolve phylogenetic relationships among members of the subfamily delphininae using two approaches.
both use complete taxon sampling for the subfamily  and incorporate multiple individuals per species to capture the intraspecific variation inherent within the species.
first, phylogenetic analysis was performed using mitochondrial dna control region sequences.
the control region was chosen for comparison because it is commonly relied upon for studies of cetacean systematics  <cit>  and species identification  <cit> .
the higher mutation rate may allow the control region locus to resolve relationships that the cytochrome b gene sequences could not  <cit> .
however, mtdna sequences represent a single locus gene tree and phylogenetic reconstruction of species trees can be greatly improved through multi-locus analyses.
to address this problem, we employed an alternate approach to phylogenetic reconstruction targeting multiple polymorphic markers from anonymous sites across the genome.
amplified fragment length polymorphism  analysis is a powerful molecular technique combining a restriction fragment length polymorphism  assay and dna amplification via the polymerase chain reaction   <cit> .
hundreds of genomic markers can be generated from the assay's restriction enzyme digest and two rounds of fragment amplification via pcr  <cit> .
the aflp method provides a multi-locus approach and may overcome problems in phylogenetic reconstruction resulting from incomplete lineage sorting  <cit> .
aflps have proven successful for resolving phylogenetic histories at both shallow  <cit>  and deeper species level hierarchies  <cit> .
within delphinids, these markers are powerful enough to reveal differentiation between two sets of closely related taxa, delphinus delphis and d. capensis, and offshore and coastal morphotypes of t. truncatus  <cit> .
in addition, koopman  <cit>  demonstrated substantial phylogenetic signal in aflp data and congruence between the aflp markers and other nuclear data .
here we compare and contrast the efficacy of mtdna sequences and aflp markers for reconstructing evolutionary histories in a group that has undergone a rapid radiation.
methods
sampling - mtdna
a total of  <dig> control region mtdna haplotypes  representing eleven delphinine species and seven outgroup taxa were utilized for this analysis.
fourteen of the fifteen outgroup taxa haplotypes were downloaded from genbank as well as  <dig> tursiops aduncus haplotypes,  <dig> south australian t. truncatus haplotypes, and one indopacific t. truncatus haplotype .
the remaining haplotypes resulted from sequences obtained in our lab from  <dig> individual samples.
six hundred and six of these samples were sequenced for previous studies:  <dig> delphinus spp.
samples,  <dig> s. frontalis samples, and  <dig> tursiops truncatus samples  <cit> .
unpublished sequences from genbank:
charlton k, mckechnie sw: diversity of mitochondrial dna control region of bottlenose dolphins  from southern australian waters.
biological sciences, monash university;  <dig>
grohmann l, bokermann i, unseld m, hiesel r, malek o, giese a, brennicke a: whale meat from protected species is still being sold on japanese markets.
berlin, germany;  <dig>
ji gq, yang g, liu s, zhou ky: a study on the variability of the mitochondrial dna control region of bottlenose dolphins  in chinese waters.
institute of genetic resources, college of life sciences, nanjing university;  <dig>
yamagiwa d: grampus griseus  mitochondrial d-loop region dna.
graduate school, agricultural life sciences, the university of tokyo, veterinary anatomy;  <dig>
tissue  was obtained from free ranging dolphins following the methods of gorgone et al.
<cit>  or from sampling of dead, stranded individuals; a minimum of five individuals per species was included in an attempt to encompass some geographic variation that may be present in species with large oceanic distributions .
samples from tursiops truncatus included collections from both of the distinct morphotypes found in the northwest atlantic ocean, described as the coastal form and the offshore form  <cit> .
dna from the four t. aduncus samples not obtained from genbank and ten sousa chinensis samples was provided by the national marine fisheries service, southwest fisheries science center marine mammal and sea turtle dna archive .
total number of samples  and haplotypes  in the mtdna dataset, and total number of samples in aflp dataset .
numbers in italics represent sequences from genbank .
region abbreviations: etp = eastern tropical pacific; nep = northeast pacific; wsa = western south atlantic; gom = gulf of mexico; wna = western north atlantic; ena = eastern north atlantic; bs = black sea; sa = south australia; ipo = indo-pacific ocean;  = outgroup species
sampling - aflp
a total of  <dig> samples representing eleven delphinine species plus two outgroup species- steno bredanensis  and lagenorhynchus acutus  - were included in the aflp analysis.
a minimum of four individuals per species was used and samples from different geographic areas were used when available .
tissue  was obtained as above via remote biopsy or sampling of stranded individuals.
dna from the four t. aduncus samples and ten sousa chinensis samples was provided as above.
one hundred fourteen samples from the aflp dataset were also included in the mtdna analysis; overlapping samples were incorporated for all ingroup taxa.
dna extraction
dna from  <dig> delphinus spp.
samples,  <dig> stenella frontalis samples, and  <dig> tursiops truncatus samples was extracted for previous studies  <cit> .
dna from the remaining samples  was extracted via the proteinase k method as described in  <cit>  with the exception of the buffer volume .
dna concentrations were assessed on a hoefer dyna quant  <dig> fluorometer prior to pcr amplification and dna quality was assessed via agarose gel electrophoresis.
mtdna control region sequencing
a dna fragment approximately  <dig> bp in size containing the flanking proline trna gene and 5' end of the control region was amplified using the primers l <dig> and h <dig>  <cit> .
amplification reactions contained  <dig> ng of dna,  <dig>  um of each primer,  <dig> um dntps,  <dig>  mm mgcl <dig>  1× pcr buffer  and  <dig>  u of taq dna polymerase .
the thermal cycling profile started with a  <dig> second 94°c denaturing step followed by  <dig> cycles of 94°c for  <dig> s, 55°c for  <dig> s, and 72°c for  <dig> s. a final 10-minute extension step was added after the last cycle to extend incomplete fragments.
the product was gel purified via excision from a  <dig> % low melting point agarose gel followed by agarase digestion.
cycle sequencing was performed using abi big dye terminator® fluorescent dye chemistry versions  <dig> ,  <dig> , or  <dig>   and pcr products were sequenced in both directions.
cycle sequencing products were cleaned via ethanol precipitation according to abi protocols or with centrisep spin column strips .
products were sequenced via capillary electrophoresis on abi  <dig> prism® and abi  <dig> genetic analyzers.
sequence electropherograms for forward and reverse reads of each fragment were edited using sequence navigator  <dig> .
<dig>  and a consensus of the two directions was constructed.
the consensus sequences were truncated to control region only  and haplotypes aligned by eye in seqpup  <dig> f  <cit> .
aflp assay
the aflp assay was run according to the protocol of vos et al.
<cit>  and applied biosystems  <cit>  with modifications made in kingston and rosel  <cit> .
taqi was used as the frequent cutter enzyme rather than msei for increased resolution in c/g rich vertebrate genomes  <cit> .
twenty selective primer combinations were used to generate the aflp fragments  <cit> .
ecori selective primers were fluorescently labeled for detection on an abi  <dig> prism® genetic analyzer so only fragments containing an ecori site were detected.
aflp scoring
resulting electropherograms were scored for polymorphic peaks using genotyper®  <dig>  software .
peaks were scored as dominant markers, present or absent .
one base-pair  sized bins were created for each dominant marker category.
markers ranged in size from  <dig> bp to  <dig> bp.
we used the conservative scoring protocol described in kingston and rosel  <cit>  to protect against potential problems associated with uneven amplification among samples and poor amplification of larger fragments for degraded dna samples.
only fragments that showed even amplification across all samples were scored  and scoring was halted at the size  where poorer quality samples began to lose monomorphic peaks  <cit> .
analyses - mtdna control region
sequences were aligned by eye  and unique haplotypes identified.
a model and parameters for the phylogenetic reconstruction were determined empirically using likelihood via modeltest  <dig>   <cit> .
the akaike information criterion  indicated that the tamura-nei model of dna evolution with a gamma correction , proportion of invariable sites =  <dig> , and empirical base frequencies  was most appropriate given the data.
the alignment of control region haplotypes was analyzed in a likelihood framework using garli  <dig>   <cit> .
the model and parameters determined above in modeltest  <dig>   <cit>  were applied in garli .
two replicates were performed in order to assess convergence on a topology.
stop generation and stop time were set at  <dig> ,000; genthreshfortopoterm was set at 20000; scorethreshforterm was set at  <dig> .
the remaining options were set as default.
the analysis was bootstrapped for  <dig> iterations.
additionally, the aligned control region haplotypes were analyzed in paup*  <dig> b <dig>  <cit>  using distance methods.
using the above model and the neighbor-joining algorithm, a phylogenetic reconstruction was rendered and bootstrapped  using  <dig> iterations.
within species, between species and corrected between species genetic distances were estimated using mega  <dig>  and the tamura-nei model with parameters as described above.
the trees were rooted with a single outgroup haplotype  although six more outgroup taxa  were included in the analysis .
analyses - aflp
an initial binary data matrix was compiled for all individuals and dominant aflp markers.
this presence-absence matrix was used as the basis for all aflp analyses.
species-specific markers were defined directly from the raw binary data; a species-specific marker was shared by all individuals sampled from a particular species to the exclusion of all other taxa .
relationships among taxa and individuals were defined using a neighbor-joining phylogram  <cit>  which was built using nei-li distance derived from the binary data matrix and bootstrapped   <dig> times using paup*  <dig> b <dig>  <cit> .
distance-based methods were included because the parsimony criterion, in particular, may be inappropriate for use with dominant, anonymous markers due to the inherent faulty assumption of homology among shared absent markers and the possible parsimonious, but incorrect, reconstruction in which no markers are assigned to an ancestor at a given internal node  <cit> .
we also utilized the bayesian phylogenetic inference method through mrbayes  <dig>   <cit> .
due to the binary nature of the aflp data and the difference in probabilities between  <dig> to  <dig> and  <dig> to  <dig> state changes, we chose the restriction "noinvariantsites" option and the "noabsentsites" option.
the analysis was run over  <dig> replicates to assure convergence on a topology; each run was performed over  <dig> , <dig> generations  and burn-in was set at  <dig>  generations.
the remaining options were set as default.
non-metric multidimensional scaling analysis  was used to further clarify and visualize relationships among taxa outside the context of a bifurcating tree.
nmds is an ordination technique designed to portray relationships as defined by a jaccard similarity matrix in three-dimensional space.
a jaccard similarity matrix was created from the initial binary data using ntsys-pc  <cit> .
jaccard similarity values range from  <dig>  to  <dig>  and are based on shared presence of markers:
where a is the number of polymorphic markers shared by individuals x and y, b is the number of markers present in x but absent in y, and c is the number of markers present in y but absent in x  <cit> .
nmds plots were created using ntsys-pc  <cit> ; three sets of principal coordinates analysis values were used as the initial configuration for better fit.
the goodness of fit of the nmds model was measured using a stress value ranging from 0- <dig> .
interspecies hybrids
during the study, four spotted dolphin individuals  were genetically identified as possible interspecies hybrids due to incongruence between field identification, mtdna haplotype and placement within the phylogenetic reconstruction based on the aflp markers.
these four individuals were added to the core aflp data set in separate, identical phylogenetic analyses repeated as above .
we then utilized the program structure  <dig>  on the multilocus aflp data to assess hybrid origin and perform assignment tests for the putative hybrids  <cit> .
we first ran structure using the aflp genotypes of what we considered pure s. attenuata  and s. frontalis  samples to confirm the number of 'populations' and their composition without the putative hybrids.
these data were analyzed via the admixture model, independent alleles,  <dig>  burn-in and  <dig>  mcmc, three replicate runs at each k , lambda =  <dig>  recessive alleles option, and no prior population information.
this test was repeated with inferred lambda as well as under the no admixture model.
next, the hybrids were added and prior population info was used for the  <dig> "known" s. attenuata and s. frontalis samples; this dataset was used to identify affinities of each of the four putative hybrids to either species.
this analysis was run under both admixture and no admixture models,  <dig>  burn-in and  <dig>  mcmc, migration  <dig> .
finally, that same analysis was run without prior population info for all samples included in order to examine probabilities of alternative parentage and grandparentage.
RESULTS
mtdna control region
the  <dig> delphinid sequences -  <dig> sequenced for previous studies  and  <dig> sequenced in this study - resulted in  <dig> unique  <dig> bp mtdna control region haplotypes.
with the addition of outside sequences from genbank, a total of  <dig> haplotypes were present in the alignment .
two t. aduncus sequences  from our study matched the downloaded genbank haplotypes af <dig> and af <dig>  haplotypes from our study submitted to genbank bear accession numbers dq845437-dq <dig> and gq504058-gq <dig>  the maximum within-species diversity was equal to or exceeded net between-species distances in  <dig> of  <dig>  pairwise comparisons .
net mean between-species distance ranged from  <dig> - <dig>  while mean within-species distances ranged from  <dig>  to  <dig> .
net mean between species distance below diagonal, mean within species distance along diagonal , uncorrected mean between species distance above diagonal.
in brackets, standard error based on  <dig> bootstrap replicates.
the phylogenetic trees derived from the mtdna control region data reveal that although some species form monophyletic groups , these nodes are not unilaterally supported after bootstrap resampling .
stenella attenuata and sousa chinensis are the only species that show robust bootstrap support .
the control region trees offer little resolution with regard to branching order among species even without bootstrap support.
several d. delphis and s. coeruleoalba haplotypes fall outside any coherent species clade .
aflp
the aflp assay rendered  <dig> total polymorphic markers among  <dig> individuals representing  <dig> species.
each primer combination exhibited on average  <dig>  ±  <dig>  polymorphic markers .
each individual exhibited on average  <dig>  ±  <dig>   polymorphic markers when all primer combinations were considered.
five of the  <dig> polymorphic characters were species-specific markers as they demonstrated species-level synapomorphies.
both aflp phylogenies are more resolved than the phylogenies derived from the mtdna data, although branching order among some of the deeper nodes is not determined .
the bayesian analysis of the aflp data resolves the deeper nodes best.
both d. capensis and s. longirostris are nested within taxa, d. delphis and s. clymene, respectively.
the two spotted dolphin species, s. frontalis and s. attenuata, are joined with very high support, although on a species level, only s. frontalis is monophyletic upon bootstrapping.
these sister species do not cluster exclusively with other stenella species, rendering the genus polyphyletic .
tursiops truncatus coastal and offshore morphotypes from the atlantic ocean form two distinct groups; the node joining these two sister taxa exhibits excellent support.
an unsupported node clustering the rest of the subfamily suggests sousa chinensis is the sister taxon to the rest of the delphinines in the distance analysis.
however, the placement of sousa chinensis is the most notable point where the two analyses disagree; the bayesian tree suggests this species clusters with the spotted dolphins and t. truncatus.
the entire subfamily clusters to the exclusion of steno bredanensis with excellent bootstrap support in both reconstruction methods.
non-metric multidimensional scaling analysis reveals large-scale differences between the sister taxa s. frontalis and s. attenuata and the rest of the delphininae .
the genus tursiops is slightly less distant from the remaining delphinines.
stenella coeruleoalba, l. hosei, and sousa chinensis form distinct species groups, while the overlapping s. longirostris and s. clymene clusters associate closely with delphinus spp.
although s. longirostris and s. clymene cannot be distinguished on the delphinine nmds plot due to the variance among s. clymene, when the nmds analysis was run for only these taxa, the two species formed distinct groups .
the same has been demonstrated for sister taxa pairs d. delphis, d. capensis, as well as coastal and offshore morphotypes of t. truncatus  <cit> .
tursiops aduncus, although represented by only four individuals, does not form a tight species cluster.
two of the four t. aduncus representatives associate more closely with the t. truncatus group than any other group.
interspecies hybrids
initial structure analyses confirmed two 'populations'  for k = 1-5: - <dig> , - <dig> , - <dig> , - <dig> , - <dig> , respectively) from the aflp data  corresponding to the two discrete species .
probabilities of identity of the four putative hybrids based on the aflp data were opposite the species identifications inferred from the mtdna data, but congruent with field identifications based on external morphological characters .
an additional hybrid, for which the mtdna, field data and aflp data are all congruent, may also have been identified.
* posterior probability that each individual has recent ancestry  from the species opposite that identified by the aflp analysis.
two structure assignment and ancestry analyses were performed, one with prior population information assigned and the second with putative hybrid population identity set as unknown to examine alternate parentage contributions.
last row identifies fifth possible hybrid discovered by the analysis.
discussion
mitochondrial versus aflp phylogenies
traditional methods of phylogenetic inference using a single gene, or even several genes, often yield a limited picture of evolutionary history for closely related species.
by focusing on a scale so fine and specific, discrepancies between loci are the norm  <cit> .
the difficulty of finding a gene that reflects complete lineage sorting, yet is variable enough to reliably untangle relationships between closely related taxa that are the product of a recent radiation event is apparent in many phylogenetic studies  <cit> .
with the advent of more powerful analyses that consider hundreds of nuclear markers, the resolution and power of phylogenetic analyses has drastically increased  <cit> .
here we compared the ability of two different molecular marker types, mtdna sequences and aflp markers, to reconstruct the evolutionary history among a group of recent and likely rapidly radiated taxa.
the comparison of mitochondrial and aflp phylogenies reveals congruence in some cases but discrepancies in other inferences concerning the relationships among members of the subfamily delphininae.
overall, the mtdna control region phylogeny, similar to previously published cytochrome b phylogenies  <cit> , offers little power for resolving relationships among these delphinine taxa .
in the control region, the high levels of within species diversity, often equal to or exceeding levels of between-species divergence , may significantly interfere with the ability to construct a robust phylogeny.
in some cases, pairwise between-species distances are less than the distance between two individuals of the same species.
the aflp reconstructions provide stronger support for many relationships within the subfamily and untangle other relationships that mtdna leaves unresolved .
discrepancies between mitochondrial and nuclear phylogenies are not unique to this study  <cit> .
although both sets of data exhibit a pattern of high levels of intraspecific variation and low levels of interspecies divergence, phylogenetic inference from the mtdna data in particular is hindered by this pattern.
one explanation for this pattern is a recent, rapid radiation, with incomplete lineage sorting of mtdna control region haplotypes among the delphinine species  <cit> .
the occurrence of the marked haplotypes that do not cluster together with conspecifics  could be explained by such a hypothesis.
the fact that all of the "misplaced" haplotypes in the mtdna phylogeny group with the expected species in the aflp tree indicates they were not simply misidentified samples.
although this pattern might also arise if there were a recent nuclear duplication of this mitochondrial gene  <cit> , there is no evidence to date of such a control region duplication in cetaceans.
whatever the cause, the mtdna control region tree alone is insufficient for robust phylogenetic inference in this subfamily.
although the resolution for individual species groups is better in the aflp trees, the deeper internodes are short and not consistently well-supported.
this short internode congruence with the mtdna pattern suggests rapid radiation over other alternatives.
in the mtdna analysis, few species form monophyletic groups in the control region phylogeny and only two species, sousa chinensis and stenella attenuata, render monophyletic species clades consistently supported.
while it is possible to use the mtdna control region sequence for species identification if an unknown falls into one of these groups, or even a monophyletic species clade lacking bootstrap support , interpretation is much less clear if an unknown haplotype falls outside a coherent monophyletic group.
use of the aflp markers appears to be a more consistent choice for species identification when using phylogenetic methods.
it has proven successful in identification of maple species and individuals  <cit> .
even in the dolphin species, which have considerably lower levels of nuclear variability  <cit> , we found no two delphinine individuals to have identical aflp profiles; this specificity suggests the method may prove useful for identifying individuals.
however, the methodology is considerably more labor intensive than dna sequencing of an mtdna fragment, and may prove less useful in cases of highly degraded tissue samples.
given the low divergence among taxa  and high levels of diversity within species, it is not entirely surprising the single locus mtdna phylogenetic reconstructions offer little resolution.
due to the prevalence of mitochondrial control region data collected, as well as the assertion that mitochondrial markers alone are excellent in phylogeographic context  <cit> , we felt it necessary to demonstrate this evolutionarily interesting example in which the data fail.
because there is so much morphological variation within delphinine species, one goal of this study was thorough taxon sampling.
since the great majority of the variation in the control region is limited to the region we amplified  <cit> , one might be concerned about resolving power for so many individuals over the short stretch of dna.
if the number of taxa is pruned along the branches of the tree, we still observe the same poor resolution .
interspecies hybrids
interspecies hybrids among delphinid taxa have been recorded, but most documented cases come from animals in captivity rather than in the wild  <cit> .
in our study, the combined mtdna and aflp data, along with corresponding morphological field identifications, indicate at least four inter-species hybrids between the two spotted dolphin species.
these two species, although morphologically similar  <cit> , exhibit ten fixed differences among mtdna control region haplotypes  <cit> .
the s. attenuata group is one of only two species exhibiting strong bootstrap support in the mtdna trees.
although the resolution is generally poor in the mtdna phylogenies, we can at least reliably discern s. attenuata from all other delphinines based on mitochondrial data alone.
in addition, none of the putative hybrid haplotypes are unique, i.e., they are shared by other individuals of that species in the data set.
in an aflp phylogenetic reconstruction, the hybrid sample exhibiting a s. frontalis mtdna control region haplotype sits in a monophyletic clade with the s. attenuata, while the three hybrid samples bearing s. attenuata mtdna control region haplotypes fall inside the bootstrap-supported s. frontalis clade .
the mostly nuclear aflp data are concordant with field identification of the samples but conflict with mtdna haplotypes.
this pattern suggests we have detected hybridization in the wild in both directions, with possible backcrossing into the paternal species.
it is important to note we would likely not be able to detect backcrossing of hybrids into the maternal species, at least with mtdna sequence incongruence.
the structure assignment and ancestry tests on the aflp data confirm the aflp phylogenetic results .
the results of the probability of assignment of each putative hybrid to species strongly conflicts with the mtdna haplotype identity, suggesting hybrids are backcrossed many generations into the paternal species.
we also see an interesting result with the s. attenuata individual sa <dig>  this individual shows a low probability of assignment to the s. attenuata species but rather than a high probability of assignment in s. frontalis, this individual exhibits high probability of mixed grandparentage .
this individual may represent the category of hybrids mentioned in the previous paragraph, those backcrossed into the maternal species.
even some of the unambiguous individuals show non-zero probabilities for extra-species grandparentage, suggesting low levels of allelic introgression may be widespread.
hybridization and its evolutionary role have been recently revisited in the literature now that larger nuclear datasets are increasingly available for comparison to mtdna phylogenies  <cit> .
shaw  <cit>  demonstrated the phenomenon of mtdna gene flow and in some cases complete introgression  across species boundaries in the hawaiian cricket genus laupala.
the apparent importance of interspecific gene flow in this system led shaw  <cit>  to issue a caveat about potentially misleading patterns of mtdna variation among closely related species complexes.
among species of darwin's finches in the galápagos archipelago, sympatric introgressive hybridization has also played an important role in the adaptive radiation of species  <cit> .
it is now recognized that interspecific hybridization in the wild is not uncommon in rapidly radiating groups  <cit> .
a reader might form an intuitive hypothesis regarding the detection of hybridization in wild populations: hybridization levels must be relatively high if even a handful of hybrids are detected by chance.
therefore, species identities should be lost over time due to gene flow.
while this may be a possible fate for some lineages through the course of evolutionary history, there is a myriad of recent literature documenting cases of just the opposite: divergence with gene flow and long-term maintenance of species boundaries in the face of secondary contact  <cit> .
the highly labile spotted dolphin species may demonstrate a stable hybrid zone across yet to be quantified gradients or patches in marine variables such as salinity, temperature, depth, or prey distribution.
differential introgression of loci across the genome between species due to hybridization is possible and common; gene flow and species divergence are not always mutually exclusive trajectories  <cit> .
although we cannot make definitive conclusions about the role of hybridization in the evolution of the delphinine species, it is important to note the evolutionary similarities among the delphininae, the laupala, and the darwin's finches.
the delphininae, like these other taxa, are likely the product of a recent, rapid radiation event.
in addition, many of the species in the subfamily are distributed across ocean basins.
in contrast with island species, it is more difficult to discern obvious barriers to gene flow in marine species.
however, isolated allopatric or parapatric populations exist within delphinine species that are in turn sympatric with other delphinine species  <cit> .
these features outline a system in which rare, interspecies, introgressive hybridization could possibly play an evolutionary significant role among the delphinines.
at the very least, the occurrence of interspecies hybridization in the wild draws into question the practice of species identification for some delphinines based on mtdna sequences alone.
systematic relationships inferred from aflp data
the taxonomy and systematics of the family delphinidae have been unresolved for centuries.
rice  <cit>  recognized  <dig> species, but new species of sotalia and orcaella have recently been described  <cit>  and additional species will likely be described in coming years  <cit> .
it is clear that the biodiversity of this family remains underestimated at this time.
at the subfamily level, as few as two and as many as five subfamilies have been proposed  <cit> , but the author suggests that subfamily designations are "best held in abeyance" pending further studies  <cit> .
the subfamily delphininae, however, has maintained itself across the variant classifications with major conflicts of membership confined to the genera grampus, steno, sousa and sotalia  <cit> .
it seems generally agreed that grampus belongs in the subfamily globicephalinae  <cit> , that sousa belongs within the delphininae  <cit> , and sotalia may be at least a sister lineage  <cit> .
however, to date, no comprehensive morphological or molecular study has had sufficient power to satisfactorily resolve the relationships among the species in the delphininae.
most analyses do not have dense within taxon sampling or complete taxon sampling for the entire group  <cit> .
some recent single-marker analyses aimed at larger phylogenetic groups may appear to provide some resolution of internal nodes within the subfamily delphininae  <cit> .
however, when a single haplotype represents each species, we are simply looking at a gene tree among the single lineages chosen by chance from those species' wide geographic distributions.
as figure  <dig> demonstrates, the stochastic sampling of a single haplotype per species could render a myriad of different results, some more resolved than others, but none alone particularly accurate at estimating the true course of the species' evolutionary history.
to avoid the possibility of this particular confounding possibility, the comprehensive within species sampling scheme for all the delphininae was implemented in our study.
as mentioned above, even complete taxon sampling encompassing geographic variability within species does not improve the mtdna phylogeny based on control region sequences.
however, the bayesian analysis of the aflp data, which provide many loci spread across the genome, provides a more resolved phylogeny for the subfamily in which many relationships are congruent with previously described relationships based on morphological evidence  <cit> .
we believe the evidence points to a real incongruence between the gene history of the mitochondrion and the evolutionary history of these recently diverged species.
morphological data support a sister taxon relationship between the two spotted dolphin species.
stenella frontalis and s. attenuata overlap in every morphological character analyzed by perrin et al.
<cit>  with the exceptions of total vertebral count and color pattern .
although there are no fixed differences in skull morphometrics, the two species can be differentiated using simultaneous discriminant analysis of multiple characters  <cit> .
neither the mtdna control region nor the most comprehensive cytochrome b  <cit>  phylogenies indicate these two species share a most recent common ancestor relative to other stenella species.
in the mtdna control region phylogeny, s. attenuata haplotypes are so distinct that they are the only widely sampled species within the delphininae supported by a bootstrap value approaching 100%, but they never cluster with s. frontalis.
in contrast, the aflp data recover a relationship congruent with the morphological evidence with bootstrap support of 100%, suggesting that the morphological overlap may be a result of homology rather than convergence .
the detection of hybridization between the two species  also argues for a recent shared evolutionary history.
the aflp phylogeny is the first molecular study to suggest that s. clymene and s. longirostris are more closely related to each other than any other taxa in the genus stenella.
this relationship is supported by morphological data as well  <cit> .
stenella clymene was only officially recognized as a valid species in  <dig> when perrin et al.
<cit>  examined a series of skulls and photos of s. clymene, s. longirostris and s. coeruleoalba.
prior to this, the uncertainty in the validity of the species likely resulted from the fact that external color patterns in s. clymene resemble s. longirostris while the shape of the skull more closely resembles s. coeruleoalba  <cit> .
by examining a series of skulls of all three taxa, perrin et al.
<cit>  determined that s. clymene is a valid species and concluded it is most closely related to s. longirostris.
interestingly, only s. clymene and s. longirostris are known to exhibit aerial spinning behavior and perrin  <cit>  has suggested this may represent a synapomorphic character.
the aflp data support this close relationship, in fact, with the s. longirostris samples being nested within the s. clymene clade, suggesting a recent common ancestry for these two taxa.
this result contrasts with the leduc et al.
cytochrome b phylogeny  <cit>  that grouped s. clymene with s. coeruleoalba.
the authors suggested s. clymene may be of hybrid origin from parental species s. coeruleoalba and s. longirostris.
however, the aflp phylogeny does not place s. clymene intermediate to the other two species and in fact, groups s. clymene with s. longirostris and delphinus to the exclusion of s. coeruleoalba, thereby arguing against such a hybrid origin.
in addition, perrin and mead  <cit>  suggested a close relationship between s. clymene, s. longirostris and delphinus, s. coeruleoalba and lagenodelphis hosei.
the multidimensional scaling analysis of the aflp data mirrors the morphological inferences ; the phylogenetic analysis supports a relationship among s. clymene and s. longirostris and delphinus, with s. coeruleoalba as the next most closely related species to this group.
however, l. hosei is sister to all the delphinines except t. aduncus in both aflp trees, a result in conflict with the morphological data as well as the mtdna and nuclear data in caballero et al.
<cit> .
the caballero et al.
analysis of combined genes suggests, like the morphology, a close relationship between l. hosei and d. delphis and s. longirostris.
however, the support for this node is moderate  and five of the eleven species in the subfamily are not included in their molecular analysis  <cit> .
further work is needed to understand the conflicting position of l. hosei.
contrary to morphological evidence, most molecular studies  <cit>  and the aflp phylogenies all suggest that the genus tursiops is polyphyletic.
in the aflp tree, the monophyletic coastal and offshore t. truncatus clades show substantial divergence, but share a common ancestor to the exclusion of all other delphinine taxa .
the divergence between the coastal and offshore morphotypes of t. truncatus has been previously documented  <cit> .
this divergence is greater than the divergence seen between two recognized species d. delphis and d. capensis in both mtdna and aflp data and may represent species level differentiation  <cit> .
there is no support for a close relationship between t. truncatus and t. aduncus, the other member of the genus; surprising given the morphological similarities between the two species.
recent genetic studies have suggested the existence of additional species within the genus tursiops  <cit> .
given the polyphyletic nature revealed by molecular analyses, this genus deserves further investigation.
in addition, the aflp data do not provide strong evidence for monophyly of the genus stenella.
the validity of this genus as a coherent evolutionary lineage has been questioned previously on both morphological and molecular grounds  <cit> .
while nuclear, mtdna sequence data and aflp data do not support monophyly for this genus, the relationships among stenella species and other members of the subfamily do differ significantly between the data types.
the bayesian analysis of the aflp data is most in line with morphological inferences.
the bayesian analysis groups s. longirostris and s. clymene with the two delphinus species while the two spotted dolphins species are found together in a well supported group with t. truncatus and sousa chinensis .
these relationships are congruent with morphological evidence, which supports a close relationship between members of the genus tursiops, s. attenuata, and s. frontalis  <cit> .
perrin et al.
<cit>  describe a suite of cranial characters that differentiate t. truncatus, s. frontalis and s. attenuata from other delphinines.
stenella frontalis and t. truncatus also share a similar ground coloration pattern  <cit> .
all molecular data sets do support monophyly of the genus delphinus  <cit> .
in the aflp phylogenies, the d. capensis samples are nested within the d. delphis clade.
the recent divergence between the two delphinus species has been demonstrated in previous studies using both mtdna and morphological data  <cit> .
the fact that d. capensis is nested within d. delphis in the genome-wide survey as well as in the mtdna phylogenies further supports the hypothesis of incipient speciation in d. capensis; reciprocal monophyly has not yet been attained, while at the same time there is no evidence of hybridization between these species even where they are sympatric  <cit> .
one major inconsistency in the aflp analyses is the placement of sousa chinensis.
the distance-based approach positions this species on a fairly long branch but an unsupported node placed outside of the remaining delphininae taxa .
the bayesian analysis of the aflp data positions sousa firmly within the delphininae cluster, grouping it with t. truncatus and the two spotted dolphin species with posterior probability of 91% .
while this placement is the main difference between the two topologies, the incongruence is not drastic.
the main backbone of the distance-based tree is composed mostly short internodes and little bootstrap support.
while the bayesian inference offers greater resolution, the linking of sousa as a sister lineage to the pair stenella frontalis and s. attenuata alone bears a low posterior probability.
the genus sousa has historically been considered a more primitive taxon grouped with genera generally considered outside the subfamily delphininae .
as noted by leduc et al.
<cit> , arnold and heinsohn  <cit>  suggested similarities in morphological characters among sousa, tursiops and stenella, perceiving more derived characters in sousa than previous investigators.
the bayesian analysis of the aflp data would support this hypothesis.
neither analysis, however, suggests that this genus belongs outside the subfamily, supporting other molecular studies using mtdna and nuclear genes  <cit> .
investigating further the relationship between the sotalia lineage and the delphininae with complete species coverage and dense within species sampling may help clarify the placement of the genus sousa  <cit> .
CONCLUSIONS
rapid radiation and implications for phylogeny reconstruction
although there are still ambiguities to be resolved concerning relationships among some of the taxa in the subfamily delphininae, the comparison of extensive mtdna and aflp datasets and their resultant phylogenies offers considerable insight into this enigmatic group.
all molecular data to date support the hypothesis of a recent, rapid radiation in the evolutionary history of these taxa.
the support for a rapid radiation in this subfamily is of considerable interest.
how did it come about?
changes in sea level and concomitant changes in water temperatures are often suggested as having played a significant role in speciation of cetaceans  <cit> .
such changes could explain the atlantic endemics s. clymene and s. frontalis, as described earlier.
why these species did not expand outside the atlantic is not clear.
however, currents around the cape of good hope would favor movement into the atlantic rather than out.
it is more difficult to explain the rest of the delphininae, where many species are both sympatric and distributed circumglobally.
it is clear, however, that the oceans support a diverse array of similar dolphin species and hence must have produced conditions necessary for the diversification seen in this subfamily.
careful, comprehensive analyses of habitat and diet, as has been recently done with beaked whales  <cit> , may provide insight into the different processes and pressures that produced these closely related dolphin taxa.
the pattern of high intraspecific variation and low interspecies divergence exhibited by the delphinines can be problematic when dealing with single-locus mtdna phylogenetic reconstruction.
bootstrapping of mtdna datasets results in nearly completely unresolved phylogenies that are highly unstable when changes in outgroup, or even ingroup, taxa are made .
hence, inferences about the evolutionary relationships among the delphinine species cannot be made from these reconstructions.
the multi-locus approach using aflp markers offers far greater resolving power in the face of this pattern.
given the reduced effective population size of a mitochondrial marker, we often expect better phylogenetic resolution from mtdna data.
in this case, we find the opposite.
since the mitochondrion bears the evolutionary history of a single molecule nested within a species, the signal embedded, while often powerful, must be interpreted with caution.
the aflp markers are a cross section of markers across the genome, providing the phylogenetic reconstruction with a signal that is integrated across sites likely to be neutral, selected, linked to selected sites, and everything in between; the signal incorporates heterogeneity in gene histories.
it is from this attribute that the assay may draw resolving power.
finally, the aflp data when coupled with mtdna sequence data provide evidence of interspecies hybridization.
if interspecies hybridization plays a part in delphinine evolution, even if rare in frequency, extreme caution must be used when inferring phylogeny from mtdna loci in the absence of corresponding multi-locus nuclear data.
considering the extensive process of developing sequence-based, single-copy nuclear markers useful for phylogenetic reconstruction on this scale  <cit> , the multi-locus aflp approach offers us a powerful tool with which to begin addressing the problems associated with phylogenetic inference in closely related, recently diverged taxa.
authors' contributions
sek, lda, and per all carried out aspects of the molecular labwork, data analyses, and manuscript composition.
sek ran the bulk of the aflp data collection and analyses, while lda and per ran much of the sequence data collection and analyses.
lda completed the structure analyses.
sek executed the garli analyses.
per was the primary force behind the study concept.
all authors read and approved the final manuscript.
supplementary material
additional file 1
phylogenetic analysis of truncated dataset.
phylogenetic tree  inferred from a pruned dataset of mtdna control region haplotypes  using the maximum likelihood tool garli.
tree is rooted with the outgroup lagenorhynchus acutus.
ml bootstrap values  above 50% are denoted above supported nodes.
click here for file
acknowledgements
we would like to thank the many people involved with the sampling effort for this project: a. hohn, l. hansen, a. gorgone, p. martone, r. lopiccolo, g. lovewell, k. touhey, r. baird - noaa fisheries sefsc, beaufort, nc; l. youngner - georgia dnr; c. potter - smithsonian institution; t. jefferson and ocean park, aberdeen, hong kong, s. chivers, p. morin, k. roberston - noaa fisheries swfsc, la jolla, ca; j. nicolas and biopsy crews - noaa fisheries nefsc woods hole, ma; s. barco - virginia marine science museum; e. zolman, w. mcfee nos - charleston, sc; w. hoggard, k. mullin, w. jones - noaa fisheries sefsc, pascagoula, ms; r. griffin - mote marine lab; t. martinez - noaa fisheries sefsc, miami, fl; rob nawojchik - mystic aquarium; va, nc, sc, ga stranding network members and volunteers.
thanks to a. sellas not only for sampling effort, but also for contributions  in the lab.
thanks to a. westgate for his contributions and effort  in the lab.
we appreciate the sequencing efforts of s. leonard and l. flinn.
thanks to s. sonsthagen for assistance in the execution of the garli bootstrap analysis.

 phylogenetic relationships among cetacean taxa are contended at many different levels. however, robust phylogenies are necessary for gaining insight into the evolutionary histories of these taxa and can help in understanding speciation of highly mobile taxa in an environment with seemingly few barriers to movement. studies using molecular markers, often mitochondrial dna sequences, to elucidate phylogenetic relationships among cetacean taxa have just as often created new controversy as resolved standing controversy. at deeper evolutionary levels, nuclear molecular data supported early conclusions based on morphological data that the hippopotamids are the sister lineage to cetaceans within the artiodactyla, although, more recent morphological analyses do no support this relationship and instead identify raoellids as the sister group to cetaceans. within the cetacea, a previously controversial issue was the placement of sperm whales; while initial mtdna sequence analyses suggested the sperm whale lineage is more closely related to the baleen whales rather than the rest of the toothed whales, analysis of mtdna cytochrome b data, multiple mtdna sequences and morphology and nuclear data all support the traditional view of monophyletic suborders. for some cetacean lineages, such as the beaked whales, mitochondrial sequences have proven excellent markers for resolving phylogenies and even pinpointing new species, although nuclear markers offer new evolutionary insights. for the endangered right whale species, mitochondrial markers render phylogenies congruent with those inferred from nuclear data. one group where relationships among taxa still remain unresolved is the family delphinidae, and in particular the subfamily delphininae. this subfamily encompasses eleven nominal species: tursiops truncatus, t. aduncus, stenella frontalis, s. attenuata, s. longirostris, s. clymene, s. coeruleoalba, delphinus delphis linnaeus, d. capensis gray, lagenodelphis hosei fraser, and sousa chinensis. morphological data have most often excluded the genus sotalia gray, from this subfamily in concordance with the mtdna cytochrome b phylogeny of leduc et al., but a recent molecular analysis suggested the genus should be included. however, interpreted in the most conservative manner, the combined analysis of mtdna and nuclear loci in caballero et al. supports a sister relationship between the lineages leading to the subfamily delphininae and the genus sotalia but is not necessarily strong evidence for inclusion of sotalia in the subfamily; a subtle but important distinction. a number of morphological and genetic studies have been conducted in an attempt to resolve evolutionary relationships among taxa within the family delphinidae and/or subfamily delphininae. leduc et al. and may-collado et al. used mtdna cytochrome b sequences to reconstruct relationships among the taxa within the family delphinidae. although the leduc et al. phylogeny included representatives of all species within the subfamily delphininae, each species was represented by only a few sequences; the mitochondrial cytochrome b locus was unable to completely resolve branching order within the delphininae. the lack of resolution in the cytochrome b phylogeny suggests the delphininae are the product of a recent radiation; divergence among the numerous taxa is small, rendering resolution of branching order difficult. the leduc et al. cytochrome b phylogeny also suggests polyphyly of the delphinine genera stenella and tursiops and this result for stenella is also supported by may-collado et al.. so far, nuclear and combined analyses concerning the delphininae raise more questions and provide support for only a few internal nodes within the subfamily. the polyphyletic genera and lack of branching order resolution among many of the delphinine taxa point to a need for a new approach to discerning the evolutionary relationships among the species in this subfamily and perhaps a revision of the subfamily. we attempted to resolve phylogenetic relationships among members of the subfamily delphininae using two approaches. both use complete taxon sampling for the subfamily and incorporate multiple individuals per species to capture the intraspecific variation inherent within the species. first, phylogenetic analysis was performed using mitochondrial dna control region sequences. the control region was chosen for comparison because it is commonly relied upon for studies of cetacean systematics and species identification. the higher mutation rate may allow the control region locus to resolve relationships that the cytochrome b gene sequences could not. however, mtdna sequences represent a single locus gene tree and phylogenetic reconstruction of species trees can be greatly improved through multi-locus analyses. to address this problem, we employed an alternate approach to phylogenetic reconstruction targeting multiple polymorphic markers from anonymous sites across the genome. amplified fragment length polymorphism analysis is a powerful molecular technique combining a restriction fragment length polymorphism assay and dna amplification via the polymerase chain reaction. hundreds of genomic markers can be generated from the assay's restriction enzyme digest and two rounds of fragment amplification via pcr. the aflp method provides a multi-locus approach and may overcome problems in phylogenetic reconstruction resulting from incomplete lineage sorting. aflps have proven successful for resolving phylogenetic histories at both shallow and deeper species level hierarchies. within delphinids, these markers are powerful enough to reveal differentiation between two sets of closely related taxa, delphinus delphis and d. capensis, and offshore and coastal morphotypes of t. truncatus. in addition, koopman demonstrated substantial phylogenetic signal in aflp data and congruence between the aflp markers and other nuclear data. here we compare and contrast the efficacy of mtdna sequences and aflp markers for reconstructing evolutionary histories in a group that has undergone a rapid radiation. methods sampling - mtdna a total of control region mtdna haplotypes representing eleven delphinine species and seven outgroup taxa were utilized for this analysis. fourteen of the fifteen outgroup taxa haplotypes were downloaded from genbank as well as tursiops aduncus haplotypes, south australian t. truncatus haplotypes, and one indopacific t. truncatus haplotype. the remaining haplotypes resulted from sequences obtained in our lab from individual samples. six hundred and six of these samples were sequenced for previous studies: delphinus spp. samples, s. frontalis samples, and tursiops truncatus samples. unpublished sequences from genbank: charlton k, mckechnie sw: diversity of mitochondrial dna control region of bottlenose dolphins from southern australian waters. biological sciences, monash university; grohmann l, bokermann i, unseld m, hiesel r, malek o, giese a, brennicke a: whale meat from protected species is still being sold on japanese markets. berlin, germany; ji gq, yang g, liu s, zhou ky: a study on the variability of the mitochondrial dna control region of bottlenose dolphins in chinese waters. institute of genetic resources, college of life sciences, nanjing university; yamagiwa d: grampus griseus mitochondrial d-loop region dna. graduate school, agricultural life sciences, the university of tokyo, veterinary anatomy; tissue was obtained from free ranging dolphins following the methods of gorgone et al. or from sampling of dead, stranded individuals; a minimum of five individuals per species was included in an attempt to encompass some geographic variation that may be present in species with large oceanic distributions. samples from tursiops truncatus included collections from both of the distinct morphotypes found in the northwest atlantic ocean, described as the coastal form and the offshore form. dna from the four t. aduncus samples not obtained from genbank and ten sousa chinensis samples was provided by the national marine fisheries service, southwest fisheries science center marine mammal and sea turtle dna archive. total number of samples and haplotypes in the mtdna dataset, and total number of samples in aflp dataset. numbers in italics represent sequences from genbank. region abbreviations: etp = eastern tropical pacific; nep = northeast pacific; wsa = western south atlantic; gom = gulf of mexico; wna = western north atlantic; ena = eastern north atlantic; bs = black sea; sa = south australia; ipo = indo-pacific ocean; = outgroup species sampling - aflp a total of samples representing eleven delphinine species plus two outgroup species- steno bredanensis and lagenorhynchus acutus - were included in the aflp analysis. a minimum of four individuals per species was used and samples from different geographic areas were used when available. tissue was obtained as above via remote biopsy or sampling of stranded individuals. dna from the four t. aduncus samples and ten sousa chinensis samples was provided as above. one hundred fourteen samples from the aflp dataset were also included in the mtdna analysis; overlapping samples were incorporated for all ingroup taxa. dna extraction dna from delphinus spp. samples, stenella frontalis samples, and tursiops truncatus samples was extracted for previous studies. dna from the remaining samples was extracted via the proteinase k method as described in with the exception of the buffer volume. dna concentrations were assessed on a hoefer dyna quant fluorometer prior to pcr amplification and dna quality was assessed via agarose gel electrophoresis. mtdna control region sequencing a dna fragment approximately bp in size containing the flanking proline trna gene and 5' end of the control region was amplified using the primers l and h. amplification reactions contained ng of dna, um of each primer, um dntps, mm mgcl 1× pcr buffer and u of taq dna polymerase. the thermal cycling profile started with a second 94°c denaturing step followed by cycles of 94°c for s, 55°c for s, and 72°c for s. a final 10-minute extension step was added after the last cycle to extend incomplete fragments. the product was gel purified via excision from a % low melting point agarose gel followed by agarase digestion. cycle sequencing was performed using abi big dye terminator® fluorescent dye chemistry versions,, or and pcr products were sequenced in both directions. cycle sequencing products were cleaned via ethanol precipitation according to abi protocols or with centrisep spin column strips. products were sequenced via capillary electrophoresis on abi prism® and abi genetic analyzers. sequence electropherograms for forward and reverse reads of each fragment were edited using sequence navigator. and a consensus of the two directions was constructed. the consensus sequences were truncated to control region only and haplotypes aligned by eye in seqpup f. aflp assay the aflp assay was run according to the protocol of vos et al. and applied biosystems with modifications made in kingston and rosel. taqi was used as the frequent cutter enzyme rather than msei for increased resolution in c/g rich vertebrate genomes. twenty selective primer combinations were used to generate the aflp fragments. ecori selective primers were fluorescently labeled for detection on an abi prism® genetic analyzer so only fragments containing an ecori site were detected. aflp scoring resulting electropherograms were scored for polymorphic peaks using genotyper® software. peaks were scored as dominant markers, present or absent. one base-pair sized bins were created for each dominant marker category. markers ranged in size from bp to bp. we used the conservative scoring protocol described in kingston and rosel to protect against potential problems associated with uneven amplification among samples and poor amplification of larger fragments for degraded dna samples. only fragments that showed even amplification across all samples were scored and scoring was halted at the size where poorer quality samples began to lose monomorphic peaks. analyses - mtdna control region sequences were aligned by eye and unique haplotypes identified. a model and parameters for the phylogenetic reconstruction were determined empirically using likelihood via modeltest. the akaike information criterion indicated that the tamura-nei model of dna evolution with a gamma correction, proportion of invariable sites =, and empirical base frequencies was most appropriate given the data. the alignment of control region haplotypes was analyzed in a likelihood framework using garli. the model and parameters determined above in modeltest were applied in garli. two replicates were performed in order to assess convergence on a topology. stop generation and stop time were set at,000; genthreshfortopoterm was set at 20000; scorethreshforterm was set at. the remaining options were set as default. the analysis was bootstrapped for iterations. additionally, the aligned control region haplotypes were analyzed in paup* b using distance methods. using the above model and the neighbor-joining algorithm, a phylogenetic reconstruction was rendered and bootstrapped using iterations. within species, between species and corrected between species genetic distances were estimated using mega and the tamura-nei model with parameters as described above. the trees were rooted with a single outgroup haplotype although six more outgroup taxa were included in the analysis. analyses - aflp an initial binary data matrix was compiled for all individuals and dominant aflp markers. this presence-absence matrix was used as the basis for all aflp analyses. species-specific markers were defined directly from the raw binary data; a species-specific marker was shared by all individuals sampled from a particular species to the exclusion of all other taxa. relationships among taxa and individuals were defined using a neighbor-joining phylogram which was built using nei-li distance derived from the binary data matrix and bootstrapped times using paup* b. distance-based methods were included because the parsimony criterion, in particular, may be inappropriate for use with dominant, anonymous markers due to the inherent faulty assumption of homology among shared absent markers and the possible parsimonious, but incorrect, reconstruction in which no markers are assigned to an ancestor at a given internal node. we also utilized the bayesian phylogenetic inference method through mrbayes. due to the binary nature of the aflp data and the difference in probabilities between to and to state changes, we chose the restriction "noinvariantsites" option and the "noabsentsites" option. the analysis was run over replicates to assure convergence on a topology; each run was performed over, generations and burn-in was set at generations. the remaining options were set as default. non-metric multidimensional scaling analysis was used to further clarify and visualize relationships among taxa outside the context of a bifurcating tree. nmds is an ordination technique designed to portray relationships as defined by a jaccard similarity matrix in three-dimensional space. a jaccard similarity matrix was created from the initial binary data using ntsys-pc. jaccard similarity values range from to and are based on shared presence of markers: where a is the number of polymorphic markers shared by individuals x and y, b is the number of markers present in x but absent in y, and c is the number of markers present in y but absent in x. nmds plots were created using ntsys-pc; three sets of principal coordinates analysis values were used as the initial configuration for better fit. the goodness of fit of the nmds model was measured using a stress value ranging from 0-. interspecies hybrids during the study, four spotted dolphin individuals were genetically identified as possible interspecies hybrids due to incongruence between field identification, mtdna haplotype and placement within the phylogenetic reconstruction based on the aflp markers. these four individuals were added to the core aflp data set in separate, identical phylogenetic analyses repeated as above. we then utilized the program structure on the multilocus aflp data to assess hybrid origin and perform assignment tests for the putative hybrids. we first ran structure using the aflp genotypes of what we considered pure s. attenuata and s. frontalis samples to confirm the number of'populations' and their composition without the putative hybrids. these data were analyzed via the admixture model, independent alleles, burn-in and mcmc, three replicate runs at each k, lambda = recessive alleles option, and no prior population information. this test was repeated with inferred lambda as well as under the no admixture model. next, the hybrids were added and prior population info was used for the "known" s. attenuata and s. frontalis samples; this dataset was used to identify affinities of each of the four putative hybrids to either species. this analysis was run under both admixture and no admixture models, burn-in and mcmc, migration. finally, that same analysis was run without prior population info for all samples included in order to examine probabilities of alternative parentage and grandparentage.  mtdna control region the delphinid sequences - sequenced for previous studies and sequenced in this study - resulted in unique bp mtdna control region haplotypes. with the addition of outside sequences from genbank, a total of haplotypes were present in the alignment. two t. aduncus sequences from our study matched the downloaded genbank haplotypes af and af haplotypes from our study submitted to genbank bear accession numbers dq845437-dq and gq504058-gq the maximum within-species diversity was equal to or exceeded net between-species distances in of pairwise comparisons. net mean between-species distance ranged from - while mean within-species distances ranged from to. net mean between species distance below diagonal, mean within species distance along diagonal, uncorrected mean between species distance above diagonal. in brackets, standard error based on bootstrap replicates. the phylogenetic trees derived from the mtdna control region data reveal that although some species form monophyletic groups, these nodes are not unilaterally supported after bootstrap resampling. stenella attenuata and sousa chinensis are the only species that show robust bootstrap support. the control region trees offer little resolution with regard to branching order among species even without bootstrap support. several d. delphis and s. coeruleoalba haplotypes fall outside any coherent species clade. aflp the aflp assay rendered total polymorphic markers among individuals representing species. each primer combination exhibited on average ± polymorphic markers. each individual exhibited on average ± polymorphic markers when all primer combinations were considered. five of the polymorphic characters were species-specific markers as they demonstrated species-level synapomorphies. both aflp phylogenies are more resolved than the phylogenies derived from the mtdna data, although branching order among some of the deeper nodes is not determined. the bayesian analysis of the aflp data resolves the deeper nodes best. both d. capensis and s. longirostris are nested within taxa, d. delphis and s. clymene, respectively. the two spotted dolphin species, s. frontalis and s. attenuata, are joined with very high support, although on a species level, only s. frontalis is monophyletic upon bootstrapping. these sister species do not cluster exclusively with other stenella species, rendering the genus polyphyletic. tursiops truncatus coastal and offshore morphotypes from the atlantic ocean form two distinct groups; the node joining these two sister taxa exhibits excellent support. an unsupported node clustering the rest of the subfamily suggests sousa chinensis is the sister taxon to the rest of the delphinines in the distance analysis. however, the placement of sousa chinensis is the most notable point where the two analyses disagree; the bayesian tree suggests this species clusters with the spotted dolphins and t. truncatus. the entire subfamily clusters to the exclusion of steno bredanensis with excellent bootstrap support in both reconstruction methods. non-metric multidimensional scaling analysis reveals large-scale differences between the sister taxa s. frontalis and s. attenuata and the rest of the delphininae. the genus tursiops is slightly less distant from the remaining delphinines. stenella coeruleoalba, l. hosei, and sousa chinensis form distinct species groups, while the overlapping s. longirostris and s. clymene clusters associate closely with delphinus spp. although s. longirostris and s. clymene cannot be distinguished on the delphinine nmds plot due to the variance among s. clymene, when the nmds analysis was run for only these taxa, the two species formed distinct groups. the same has been demonstrated for sister taxa pairs d. delphis, d. capensis, as well as coastal and offshore morphotypes of t. truncatus. tursiops aduncus, although represented by only four individuals, does not form a tight species cluster. two of the four t. aduncus representatives associate more closely with the t. truncatus group than any other group. interspecies hybrids initial structure analyses confirmed two'populations' for k = 1-5: -, -, -, -, -, respectively) from the aflp data corresponding to the two discrete species. probabilities of identity of the four putative hybrids based on the aflp data were opposite the species identifications inferred from the mtdna data, but congruent with field identifications based on external morphological characters. an additional hybrid, for which the mtdna, field data and aflp data are all congruent, may also have been identified. * posterior probability that each individual has recent ancestry from the species opposite that identified by the aflp analysis. two structure assignment and ancestry analyses were performed, one with prior population information assigned and the second with putative hybrid population identity set as unknown to examine alternate parentage contributions. last row identifies fifth possible hybrid discovered by the analysis. discussion mitochondrial versus aflp phylogenies traditional methods of phylogenetic inference using a single gene, or even several genes, often yield a limited picture of evolutionary history for closely related species. by focusing on a scale so fine and specific, discrepancies between loci are the norm. the difficulty of finding a gene that reflects complete lineage sorting, yet is variable enough to reliably untangle relationships between closely related taxa that are the product of a recent radiation event is apparent in many phylogenetic studies. with the advent of more powerful analyses that consider hundreds of nuclear markers, the resolution and power of phylogenetic analyses has drastically increased. here we compared the ability of two different molecular marker types, mtdna sequences and aflp markers, to reconstruct the evolutionary history among a group of recent and likely rapidly radiated taxa. the comparison of mitochondrial and aflp phylogenies reveals congruence in some cases but discrepancies in other inferences concerning the relationships among members of the subfamily delphininae. overall, the mtdna control region phylogeny, similar to previously published cytochrome b phylogenies, offers little power for resolving relationships among these delphinine taxa. in the control region, the high levels of within species diversity, often equal to or exceeding levels of between-species divergence, may significantly interfere with the ability to construct a robust phylogeny. in some cases, pairwise between-species distances are less than the distance between two individuals of the same species. the aflp reconstructions provide stronger support for many relationships within the subfamily and untangle other relationships that mtdna leaves unresolved. discrepancies between mitochondrial and nuclear phylogenies are not unique to this study. although both sets of data exhibit a pattern of high levels of intraspecific variation and low levels of interspecies divergence, phylogenetic inference from the mtdna data in particular is hindered by this pattern. one explanation for this pattern is a recent, rapid radiation, with incomplete lineage sorting of mtdna control region haplotypes among the delphinine species. the occurrence of the marked haplotypes that do not cluster together with conspecifics could be explained by such a hypothesis. the fact that all of the "misplaced" haplotypes in the mtdna phylogeny group with the expected species in the aflp tree indicates they were not simply misidentified samples. although this pattern might also arise if there were a recent nuclear duplication of this mitochondrial gene, there is no evidence to date of such a control region duplication in cetaceans. whatever the cause, the mtdna control region tree alone is insufficient for robust phylogenetic inference in this subfamily. although the resolution for individual species groups is better in the aflp trees, the deeper internodes are short and not consistently well-supported. this short internode congruence with the mtdna pattern suggests rapid radiation over other alternatives. in the mtdna analysis, few species form monophyletic groups in the control region phylogeny and only two species, sousa chinensis and stenella attenuata, render monophyletic species clades consistently supported. while it is possible to use the mtdna control region sequence for species identification if an unknown falls into one of these groups, or even a monophyletic species clade lacking bootstrap support, interpretation is much less clear if an unknown haplotype falls outside a coherent monophyletic group. use of the aflp markers appears to be a more consistent choice for species identification when using phylogenetic methods. it has proven successful in identification of maple species and individuals. even in the dolphin species, which have considerably lower levels of nuclear variability, we found no two delphinine individuals to have identical aflp profiles; this specificity suggests the method may prove useful for identifying individuals. however, the methodology is considerably more labor intensive than dna sequencing of an mtdna fragment, and may prove less useful in cases of highly degraded tissue samples. given the low divergence among taxa and high levels of diversity within species, it is not entirely surprising the single locus mtdna phylogenetic reconstructions offer little resolution. due to the prevalence of mitochondrial control region data collected, as well as the assertion that mitochondrial markers alone are excellent in phylogeographic context, we felt it necessary to demonstrate this evolutionarily interesting example in which the data fail. because there is so much morphological variation within delphinine species, one goal of this study was thorough taxon sampling. since the great majority of the variation in the control region is limited to the region we amplified, one might be concerned about resolving power for so many individuals over the short stretch of dna. if the number of taxa is pruned along the branches of the tree, we still observe the same poor resolution. interspecies hybrids interspecies hybrids among delphinid taxa have been recorded, but most documented cases come from animals in captivity rather than in the wild. in our study, the combined mtdna and aflp data, along with corresponding morphological field identifications, indicate at least four inter-species hybrids between the two spotted dolphin species. these two species, although morphologically similar, exhibit ten fixed differences among mtdna control region haplotypes. the s. attenuata group is one of only two species exhibiting strong bootstrap support in the mtdna trees. although the resolution is generally poor in the mtdna phylogenies, we can at least reliably discern s. attenuata from all other delphinines based on mitochondrial data alone. in addition, none of the putative hybrid haplotypes are unique, i.e., they are shared by other individuals of that species in the data set. in an aflp phylogenetic reconstruction, the hybrid sample exhibiting a s. frontalis mtdna control region haplotype sits in a monophyletic clade with the s. attenuata, while the three hybrid samples bearing s. attenuata mtdna control region haplotypes fall inside the bootstrap-supported s. frontalis clade. the mostly nuclear aflp data are concordant with field identification of the samples but conflict with mtdna haplotypes. this pattern suggests we have detected hybridization in the wild in both directions, with possible backcrossing into the paternal species. it is important to note we would likely not be able to detect backcrossing of hybrids into the maternal species, at least with mtdna sequence incongruence. the structure assignment and ancestry tests on the aflp data confirm the aflp phylogenetic results. the results of the probability of assignment of each putative hybrid to species strongly conflicts with the mtdna haplotype identity, suggesting hybrids are backcrossed many generations into the paternal species. we also see an interesting result with the s. attenuata individual sa this individual shows a low probability of assignment to the s. attenuata species but rather than a high probability of assignment in s. frontalis, this individual exhibits high probability of mixed grandparentage. this individual may represent the category of hybrids mentioned in the previous paragraph, those backcrossed into the maternal species. even some of the unambiguous individuals show non-zero probabilities for extra-species grandparentage, suggesting low levels of allelic introgression may be widespread. hybridization and its evolutionary role have been recently revisited in the literature now that larger nuclear datasets are increasingly available for comparison to mtdna phylogenies. shaw demonstrated the phenomenon of mtdna gene flow and in some cases complete introgression across species boundaries in the hawaiian cricket genus laupala. the apparent importance of interspecific gene flow in this system led shaw to issue a caveat about potentially misleading patterns of mtdna variation among closely related species complexes. among species of darwin's finches in the galápagos archipelago, sympatric introgressive hybridization has also played an important role in the adaptive radiation of species. it is now recognized that interspecific hybridization in the wild is not uncommon in rapidly radiating groups. a reader might form an intuitive hypothesis regarding the detection of hybridization in wild populations: hybridization levels must be relatively high if even a handful of hybrids are detected by chance. therefore, species identities should be lost over time due to gene flow. while this may be a possible fate for some lineages through the course of evolutionary history, there is a myriad of recent literature documenting cases of just the opposite: divergence with gene flow and long-term maintenance of species boundaries in the face of secondary contact. the highly labile spotted dolphin species may demonstrate a stable hybrid zone across yet to be quantified gradients or patches in marine variables such as salinity, temperature, depth, or prey distribution. differential introgression of loci across the genome between species due to hybridization is possible and common; gene flow and species divergence are not always mutually exclusive trajectories. although we cannot make definitive conclusions about the role of hybridization in the evolution of the delphinine species, it is important to note the evolutionary similarities among the delphininae, the laupala, and the darwin's finches. the delphininae, like these other taxa, are likely the product of a recent, rapid radiation event. in addition, many of the species in the subfamily are distributed across ocean basins. in contrast with island species, it is more difficult to discern obvious barriers to gene flow in marine species. however, isolated allopatric or parapatric populations exist within delphinine species that are in turn sympatric with other delphinine species. these features outline a system in which rare, interspecies, introgressive hybridization could possibly play an evolutionary significant role among the delphinines. at the very least, the occurrence of interspecies hybridization in the wild draws into question the practice of species identification for some delphinines based on mtdna sequences alone. systematic relationships inferred from aflp data the taxonomy and systematics of the family delphinidae have been unresolved for centuries. rice recognized species, but new species of sotalia and orcaella have recently been described and additional species will likely be described in coming years. it is clear that the biodiversity of this family remains underestimated at this time. at the subfamily level, as few as two and as many as five subfamilies have been proposed, but the author suggests that subfamily designations are "best held in abeyance" pending further studies. the subfamily delphininae, however, has maintained itself across the variant classifications with major conflicts of membership confined to the genera grampus, steno, sousa and sotalia. it seems generally agreed that grampus belongs in the subfamily globicephalinae, that sousa belongs within the delphininae, and sotalia may be at least a sister lineage. however, to date, no comprehensive morphological or molecular study has had sufficient power to satisfactorily resolve the relationships among the species in the delphininae. most analyses do not have dense within taxon sampling or complete taxon sampling for the entire group. some recent single-marker analyses aimed at larger phylogenetic groups may appear to provide some resolution of internal nodes within the subfamily delphininae. however, when a single haplotype represents each species, we are simply looking at a gene tree among the single lineages chosen by chance from those species' wide geographic distributions. as figure demonstrates, the stochastic sampling of a single haplotype per species could render a myriad of different results, some more resolved than others, but none alone particularly accurate at estimating the true course of the species' evolutionary history. to avoid the possibility of this particular confounding possibility, the comprehensive within species sampling scheme for all the delphininae was implemented in our study. as mentioned above, even complete taxon sampling encompassing geographic variability within species does not improve the mtdna phylogeny based on control region sequences. however, the bayesian analysis of the aflp data, which provide many loci spread across the genome, provides a more resolved phylogeny for the subfamily in which many relationships are congruent with previously described relationships based on morphological evidence. we believe the evidence points to a real incongruence between the gene history of the mitochondrion and the evolutionary history of these recently diverged species. morphological data support a sister taxon relationship between the two spotted dolphin species. stenella frontalis and s. attenuata overlap in every morphological character analyzed by perrin et al. with the exceptions of total vertebral count and color pattern. although there are no fixed differences in skull morphometrics, the two species can be differentiated using simultaneous discriminant analysis of multiple characters. neither the mtdna control region nor the most comprehensive cytochrome b phylogenies indicate these two species share a most recent common ancestor relative to other stenella species. in the mtdna control region phylogeny, s. attenuata haplotypes are so distinct that they are the only widely sampled species within the delphininae supported by a bootstrap value approaching 100%, but they never cluster with s. frontalis. in contrast, the aflp data recover a relationship congruent with the morphological evidence with bootstrap support of 100%, suggesting that the morphological overlap may be a result of homology rather than convergence. the detection of hybridization between the two species also argues for a recent shared evolutionary history. the aflp phylogeny is the first molecular study to suggest that s. clymene and s. longirostris are more closely related to each other than any other taxa in the genus stenella. this relationship is supported by morphological data as well. stenella clymene was only officially recognized as a valid species in when perrin et al. examined a series of skulls and photos of s. clymene, s. longirostris and s. coeruleoalba. prior to this, the uncertainty in the validity of the species likely resulted from the fact that external color patterns in s. clymene resemble s. longirostris while the shape of the skull more closely resembles s. coeruleoalba. by examining a series of skulls of all three taxa, perrin et al. determined that s. clymene is a valid species and concluded it is most closely related to s. longirostris. interestingly, only s. clymene and s. longirostris are known to exhibit aerial spinning behavior and perrin has suggested this may represent a synapomorphic character. the aflp data support this close relationship, in fact, with the s. longirostris samples being nested within the s. clymene clade, suggesting a recent common ancestry for these two taxa. this result contrasts with the leduc et al. cytochrome b phylogeny that grouped s. clymene with s. coeruleoalba. the authors suggested s. clymene may be of hybrid origin from parental species s. coeruleoalba and s. longirostris. however, the aflp phylogeny does not place s. clymene intermediate to the other two species and in fact, groups s. clymene with s. longirostris and delphinus to the exclusion of s. coeruleoalba, thereby arguing against such a hybrid origin. in addition, perrin and mead suggested a close relationship between s. clymene, s. longirostris and delphinus, s. coeruleoalba and lagenodelphis hosei. the multidimensional scaling analysis of the aflp data mirrors the morphological inferences; the phylogenetic analysis supports a relationship among s. clymene and s. longirostris and delphinus, with s. coeruleoalba as the next most closely related species to this group. however, l. hosei is sister to all the delphinines except t. aduncus in both aflp trees, a result in conflict with the morphological data as well as the mtdna and nuclear data in caballero et al.. the caballero et al. analysis of combined genes suggests, like the morphology, a close relationship between l. hosei and d. delphis and s. longirostris. however, the support for this node is moderate and five of the eleven species in the subfamily are not included in their molecular analysis. further work is needed to understand the conflicting position of l. hosei. contrary to morphological evidence, most molecular studies and the aflp phylogenies all suggest that the genus tursiops is polyphyletic. in the aflp tree, the monophyletic coastal and offshore t. truncatus clades show substantial divergence, but share a common ancestor to the exclusion of all other delphinine taxa. the divergence between the coastal and offshore morphotypes of t. truncatus has been previously documented. this divergence is greater than the divergence seen between two recognized species d. delphis and d. capensis in both mtdna and aflp data and may represent species level differentiation. there is no support for a close relationship between t. truncatus and t. aduncus, the other member of the genus; surprising given the morphological similarities between the two species. recent genetic studies have suggested the existence of additional species within the genus tursiops. given the polyphyletic nature revealed by molecular analyses, this genus deserves further investigation. in addition, the aflp data do not provide strong evidence for monophyly of the genus stenella. the validity of this genus as a coherent evolutionary lineage has been questioned previously on both morphological and molecular grounds. while nuclear, mtdna sequence data and aflp data do not support monophyly for this genus, the relationships among stenella species and other members of the subfamily do differ significantly between the data types. the bayesian analysis of the aflp data is most in line with morphological inferences. the bayesian analysis groups s. longirostris and s. clymene with the two delphinus species while the two spotted dolphins species are found together in a well supported group with t. truncatus and sousa chinensis. these relationships are congruent with morphological evidence, which supports a close relationship between members of the genus tursiops, s. attenuata, and s. frontalis. perrin et al. describe a suite of cranial characters that differentiate t. truncatus, s. frontalis and s. attenuata from other delphinines. stenella frontalis and t. truncatus also share a similar ground coloration pattern. all molecular data sets do support monophyly of the genus delphinus. in the aflp phylogenies, the d. capensis samples are nested within the d. delphis clade. the recent divergence between the two delphinus species has been demonstrated in previous studies using both mtdna and morphological data. the fact that d. capensis is nested within d. delphis in the genome-wide survey as well as in the mtdna phylogenies further supports the hypothesis of incipient speciation in d. capensis; reciprocal monophyly has not yet been attained, while at the same time there is no evidence of hybridization between these species even where they are sympatric. one major inconsistency in the aflp analyses is the placement of sousa chinensis. the distance-based approach positions this species on a fairly long branch but an unsupported node placed outside of the remaining delphininae taxa. the bayesian analysis of the aflp data positions sousa firmly within the delphininae cluster, grouping it with t. truncatus and the two spotted dolphin species with posterior probability of 91%. while this placement is the main difference between the two topologies, the incongruence is not drastic. the main backbone of the distance-based tree is composed mostly short internodes and little bootstrap support. while the bayesian inference offers greater resolution, the linking of sousa as a sister lineage to the pair stenella frontalis and s. attenuata alone bears a low posterior probability. the genus sousa has historically been considered a more primitive taxon grouped with genera generally considered outside the subfamily delphininae. as noted by leduc et al., arnold and heinsohn suggested similarities in morphological characters among sousa, tursiops and stenella, perceiving more derived characters in sousa than previous investigators. the bayesian analysis of the aflp data would support this hypothesis. neither analysis, however, suggests that this genus belongs outside the subfamily, supporting other molecular studies using mtdna and nuclear genes. investigating further the relationship between the sotalia lineage and the delphininae with complete species coverage and dense within species sampling may help clarify the placement of the genus sousa.  rapid radiation and implications for phylogeny reconstruction although there are still ambiguities to be resolved concerning relationships among some of the taxa in the subfamily delphininae, the comparison of extensive mtdna and aflp datasets and their resultant phylogenies offers considerable insight into this enigmatic group. all molecular data to date support the hypothesis of a recent, rapid radiation in the evolutionary history of these taxa. the support for a rapid radiation in this subfamily is of considerable interest. how did it come about? changes in sea level and concomitant changes in water temperatures are often suggested as having played a significant role in speciation of cetaceans. such changes could explain the atlantic endemics s. clymene and s. frontalis, as described earlier. why these species did not expand outside the atlantic is not clear. however, currents around the cape of good hope would favor movement into the atlantic rather than out. it is more difficult to explain the rest of the delphininae, where many species are both sympatric and distributed circumglobally. it is clear, however, that the oceans support a diverse array of similar dolphin species and hence must have produced conditions necessary for the diversification seen in this subfamily. careful, comprehensive analyses of habitat and diet, as has been recently done with beaked whales, may provide insight into the different processes and pressures that produced these closely related dolphin taxa. the pattern of high intraspecific variation and low interspecies divergence exhibited by the delphinines can be problematic when dealing with single-locus mtdna phylogenetic reconstruction. bootstrapping of mtdna datasets results in nearly completely unresolved phylogenies that are highly unstable when changes in outgroup, or even ingroup, taxa are made. hence, inferences about the evolutionary relationships among the delphinine species cannot be made from these reconstructions. the multi-locus approach using aflp markers offers far greater resolving power in the face of this pattern. given the reduced effective population size of a mitochondrial marker, we often expect better phylogenetic resolution from mtdna data. in this case, we find the opposite. since the mitochondrion bears the evolutionary history of a single molecule nested within a species, the signal embedded, while often powerful, must be interpreted with caution. the aflp markers are a cross section of markers across the genome, providing the phylogenetic reconstruction with a signal that is integrated across sites likely to be neutral, selected, linked to selected sites, and everything in between; the signal incorporates heterogeneity in gene histories. it is from this attribute that the assay may draw resolving power. finally, the aflp data when coupled with mtdna sequence data provide evidence of interspecies hybridization. if interspecies hybridization plays a part in delphinine evolution, even if rare in frequency, extreme caution must be used when inferring phylogeny from mtdna loci in the absence of corresponding multi-locus nuclear data. considering the extensive process of developing sequence-based, single-copy nuclear markers useful for phylogenetic reconstruction on this scale, the multi-locus aflp approach offers us a powerful tool with which to begin addressing the problems associated with phylogenetic inference in closely related, recently diverged taxa. authors' contributions sek, lda, and per all carried out aspects of the molecular labwork, data analyses, and manuscript composition. sek ran the bulk of the aflp data collection and analyses, while lda and per ran much of the sequence data collection and analyses. lda completed the structure analyses. sek executed the garli analyses. per was the primary force behind the study concept. all authors read and approved the final manuscript. supplementary material additional file 1 phylogenetic analysis of truncated dataset. phylogenetic tree inferred from a pruned dataset of mtdna control region haplotypes using the maximum likelihood tool garli. tree is rooted with the outgroup lagenorhynchus acutus. ml bootstrap values above 50% are denoted above supported nodes. click here for file acknowledgements we would like to thank the many people involved with the sampling effort for this project: a. hohn, l. hansen, a. gorgone, p. martone, r. lopiccolo, g. lovewell, k. touhey, r. baird - noaa fisheries sefsc, beaufort, nc; l. youngner - georgia dnr; c. potter - smithsonian institution; t. jefferson and ocean park, aberdeen, hong kong, s. chivers, p. morin, k. roberston - noaa fisheries swfsc, la jolla, ca; j. nicolas and biopsy crews - noaa fisheries nefsc woods hole, ma; s. barco - virginia marine science museum; e. zolman, w. mcfee nos - charleston, sc; w. hoggard, k. mullin, w. jones - noaa fisheries sefsc, pascagoula, ms; r. griffin - mote marine lab; t. martinez - noaa fisheries sefsc, miami, fl; rob nawojchik - mystic aquarium; va, nc, sc, ga stranding network members and volunteers. thanks to a. sellas not only for sampling effort, but also for contributions in the lab. thanks to a. westgate for his contributions and effort in the lab. we appreciate the sequencing efforts of s. leonard and l. flinn. thanks to s. sonsthagen for assistance in the execution of the garli bootstrap analysis.
however, phylogenies for taxa that have undergone recent, rapid radiation events often remain unresolved. the dolphin subfamily delphininae is a group that has likely resulted from a recent and rapid radiation. the aflp phylogeny is considerably better resolved and more congruent with relationships inferred from morphological data. the aflp data indicate a close relationship between the two spotted dolphin species and recent ancestry between stenella clymene and s. longirostris. for closely related, recently diverged taxa, a multi-locus genome-wide survey is likely the most comprehensive approach currently available for phylogenetic inference. the placement of the lagenodelphis hosei lineage is ambiguous: phenetic analysis of the aflp data is consistent with morphological expectations but the phylogenetic analysis is not. alternative methodologies for discerning evolutionary relationships under these conditions are desirable. despite several efforts, the evolutionary relationships among the species in the subfamily remain unclear.
 many molecular phylogenetic analyses rely on dna sequence data obtained from single or multiple loci, particularly mitochondrial dna loci. however, phylogenies for taxa that have undergone recent, rapid radiation events often remain unresolved. alternative methodologies for discerning evolutionary relationships under these conditions are desirable. the dolphin subfamily delphininae is a group that has likely resulted from a recent and rapid radiation. despite several efforts, the evolutionary relationships among the species in the subfamily remain unclear.  here, we compare a phylogeny estimated using mitochondrial dna control region sequences to a multi-locus phylogeny inferred from polymorphic genomic markers obtained from amplified fragment length polymorphism analysis. the two sets of phylogenies are largely incongruent, primarily because the mtdna tree provides very poor resolving power; very few species' nodes in the tree are supported by bootstrap resampling. the aflp phylogeny is considerably better resolved and more congruent with relationships inferred from morphological data. both phylogenies support paraphyly for the genera stenella and tursiops. the aflp data indicate a close relationship between the two spotted dolphin species and recent ancestry between stenella clymene and s. longirostris. the placement of the lagenodelphis hosei lineage is ambiguous: phenetic analysis of the aflp data is consistent with morphological expectations but the phylogenetic analysis is not.  for closely related, recently diverged taxa, a multi-locus genome-wide survey is likely the most comprehensive approach currently available for phylogenetic inference. 