BACKGROUND
since the lower middle ordovician some  <dig> million years ago, pterioidean bivalves have inhabited a remarkable diversity of marine epifaunal and semi-infaunal environments around the globe, typically confined to cryptic habitats and forming byssal attachments to various substrata.
for a much shorter part of their history-from the dawn of human culture-they have become the primary source of pearls and nacre  <cit> .
pearl fishing based on natural populations of pterioidean bivalves had expanded into a growing global industry with a current value of approximately half a billion us dollars per annum  <cit> .
in pursuit of pearls, commercial introductions and accidental transport of pterioideans beyond their native distribution ranges have greatly affected population dynamics of some species endangering local indigenous biotas .
a recent interest in physical properties of mother-of-pearl made pterioidean bivalves a model system for elucidating molecular mechanisms of biomineralization with medical applications for bone regeneration .
given the long and many-sided history of pterioideans and humans, and the current economic and ecological significance of these bivalves, surprisingly little is known about the standing alpha-diversity, distribution, and evolutionary history of the group.
the present-day diversity of the superfamily pterioidea gray,  <dig>  <cit>  encompasses four families traditionally defined by shell shape and ligament structure: pteriidae, isognomonidae, malleidae, and pulvinitidae  <cit> .
pteriids are distinguished by an obliquely ovate shell shape, an enlarged posteriorly projecting auricle, and a deep byssal notch.
many pteriid species, particularly of genera pteria and electroma, evolved specialized associations with hydroid, scleractinian, and alcyonarian substrata.
most species of pinctada  inhabit sandy and hard bottoms, some living in association with commensal fishes and crustaceans.
several species of pteria and pinctada are cultured for commercial pearl production.
most species of a malleid genus malleus are irregularly shaped to conform to narrow crevices in hard coral and rocky substrata, and develop elongated extensions of the hinge for stabilization in soft sediment.
species of the second malleid genus, vulsella, are adapted for living within sponges.
isognomonids are characterized by the distinct morphology of the ligament that contains multiple grooves  for ligamental attachments arranged sequentially along the hinge line .
species of the genus isognomon often co-occur with species of malleus in crevices and on soft muddy bottoms, but also are gregarious on mangrove roots, whereas species of crenatula, the second isognomonid genus, live inside sponges .
similarly to isognomonids, pulvinitids possess multivincular ligament, but are distinguished from isognomonids by the presence of the foramen, an opening in the right valve through which the byssus is protruded.
living pulvinitids are represented by the genus pulvinites, the sole extant species of which, pulvinites exempla, lives byssally attached to vertical hard substrata at 200- <dig> m depths.
previous studies aimed at resolving higher-level phylogeny within bivalvia agreed on the monophyly of the pterioidea and its placement within the subclass autolamellibranchiata.
however, most aspects of relationships within the superfamily and the identity of its immediate sister group remain uncertain.
in a pioneering study on the evolution and ontogeny of the bivalve shell, jackson  proposed the derivation of major pterioidean groups from a pteria¬-like paleozoic ancestor  <cit> .
in his scheme, several autolamellibranchiate lineages, including common oysters of the family ostreidae and several extinct families, evolved from isognomon, the latter ultimately arising from the pteria stem lineage.
jackson also provisionally recognized close relatedness of malleus and vulsella.
despite the implication of jackson's work for pterioidean systematics, the notion of the monophyly and taxonomic composition of pterioidean families persisted for over a century.
molecular studies aimed at establishing relationships within the subclass autolamellibranchiata or the entire bivalvia suggested that several widely accepted pterioidean taxa possibly did not constitute natural groups.
the minimum evolution tree based on the analysis of the mitochondrial cytochrome oxydase i  sequences did not resolve the relationships among pteria, pinctada, and isognomon   <cit> .
thus, this analysis failed to provide support for the monophyly of pteriidae but supported the monophyly of pinctada and isognomon.
a study with more extensive sampling and based on parsimony analysis of sequence data for the nuclear small ribosomal subunit  likewise failed to recover the monophyly of pteriidae, but also that of malleidae   <cit> .
in that study, isognomon was monophyletic, pteria was paraphyletic, and pinctada was polyphyletic.
the pteriid electroma was a sister taxon to the malleid vulsella.
the maximum likelihood analysis of the same data produced a similar, but more resolved, topology recovering monophyletic pinctada   <cit> .
in a more extensive in character-but not taxonomic-sampling analysis of combined 18s and the large ribosomal subunit sequence  data, the families pteriidae and malleidae, and the genus pteria were non-monophyletic   <cit> .
when these data were combined with the nuclear histone h <dig> gene and the mitochondrial coi gene sequences, the analysis produced a polytomy for all representatives of the pterioidea except for the isognomon and electroma/vulsella clades   <cit> .
in summary, molecular studies, varying greatly in sampling and analytical techniques, produced inconsistent hypotheses of pterioidean relationships with low resolution and support.
these findings, however, agreed on the monophyly of the superfamily and suggested that the families pteriidae and malleidae were not monophyletic.
the results were conflicting with regard to the monophyly of the genera pteria and pinctada, whereas isognomon was invariably monophyletic.
the monophyly of the family isognomonidae and the genera electroma and malleus were not tested, and representatives of the family pulvinitidae and the genus crenatula were not included in the analyses.
the only study that focused exclusively on the pterioidea was based on morphological characters and included representatives of all four families and nine genera   <cit> .
it supported the monophyly of the pterioidea and non-monophyly of three families .
all the genera that were represented by multiple species  were recovered monophyletic.
species-level sampling in that study was insufficient to test the assumption of monophyly of the remaining genera and did not include type species for most nominal supraspecific taxa.
despite these recent challenges on both molecular and morphological grounds, no consensus on pterioidean relationships was reached and no revision of the existing classification was proposed.
the identity of the pterioidean sister group remains uncertain, favoring the oysters , pen shells , or a clade composed of these two groups, currently classified together with the pterioidea in the order pterioida.
most analyses based on morphological data for extinct, extant, or both kinds of taxa, have suggested a sister relationship of the pterioidea and pinnoidea  <cit> , except for one, which instead suggested a closer relationship of the pterioidea and ostreoidea  <cit> .
likewise, molecular analyses were divided in favoring either pinnoidea  <cit>  or ostreoidea  <cit>  as a sister group of the pterioidea.
in other studies, either pinnoidea and ostreoidea formed a clade sister to the pterioidea  <cit> , or the relationships among the three superfamilies were unresolved  <cit> .
the only explicitly combined analysis of morphological and molecular data resolved the pinnoidea basal to the pterioidea/ostreoidea clade  <cit> .
the unsettled question of the pterioidean sister group is pertinent to the recent debates on the origin of oysters in paleontological literature  <cit> .
the present analysis is the first study to produce a nearly comprehensive species-level phylogeny for the pterioidea based on dna sequence data.
the principal specific goals of this study include  testing the monophyly of the pterioidea, and all its families and genera;  reviewing the consequences of these results for the homology of important morphological characters; and  resolving species-level phylogenetic relationships within the commercially important pearl oyster genus pinctada.
methods
taxonomic sampling
sources of specimens
most specimens were obtained from museum collections worldwide, with additional material gathered during fieldwork in the western atlantic , the eastern pacific , the mediterranean , and the indo-pacific  regions.
living animals were collected in mangroves and shallow intertidal waters, fixed in 85-95% ethanol, and subsequently transferred to 75-95% ethanol.
collected specimens were deposited at the american museum of natural history  and the field museum of natural history .
dna sequence data were obtained for  <dig> of the  <dig> terminal taxa.
the remaining sequences were acquired from genbank http://www.ncbi.nlm.nih.gov/genbank/index.html.
collection localities and museum catalog numbers for the specimens used in this study are listed in table  <dig>
museum collection acronyms: amnh, american museum of natural history, new york; cas, california academy of sciences, san francisco; fmnh, field museum of natural history, chicago; lacm, los angeles county museum; mnhn, muséum national d'histoire naturelle, paris; nmp, natal museum, pietermaritzburg; usnm, national museum of natural history, smithsonian institution, washington, d. c. ; umml, rosenstiel school of marine and atmospheric science, university of miami, florida .
na, not available.
outgroup taxa
given that the pterioidea is currently accepted as monophyletic and the pinnoidea or ostreoidea were established as the taxa most closely related to the pterioidea on morphological and molecular grounds, representatives of these superfamilies were chosen for the outgroup.
to ensure maximum diversity of the ostreoidean taxa, the selection of the gryphaeid hyotissa, and the ostreid ostrea and crassostrea exemplars were guided by their placement in a recently proposed phylogenetic hypothesis  <cit> .
to polarize the immediate outgroups, representatives of two other autolamellibranchiate orders, mytiloida and pectinoida were selected based upon the availability of sequence data.
mytilioda, represented by the well-known common mussel species, mytilus edulis, had been resolved as the most basal order of the subclass in comprehensive studies based on morphological  <cit> , molecular  <cit> , and combined  <cit>  data and, accordingly, was designated for cladogram rooting.
ingroup taxa
inasmuch as the pterioidea is the primary focus of this study, the taxonomic sampling of this taxon was the most dense.
all nominal families and genera of the pterioidea were included in the analyses.
the sample contained representatives of type species of all pterioidean genera except for pulvinites, which was described based on an extinct species, p. adansonii, from the late cretaceous.
success in securing multiple representatives for each genus varied because the material could not be obtained for few taxa characterized by cryptic life habit and presumed low species diversity.
initially, specimens were sorted to morphospecies using anatomical and shell morhological features.
species-level identifications were verified with reference to the extant type material and assigned to most specimens.
where possible, to avoid assumptions of the extent of intraspecific phenotypic variation, multiple potentially conspecific individuals were included in the analysis.
given the lack of recent systematic revisions for the majority of pterioidean taxa, species-level identifications were tentatively chosen with correspondence to their prevailing use in current literature.
the sample of the genus pinctada was complete, with  <dig> exemplars representing  <dig> nominal species, including the type species pinctada margaritifera.
the overrepresentation of pinctada species was necessitated by an enduring controversy regarding the taxonomic status of the pinctada imbricata/fucata/radiata species complex  <cit>  that includes important perliculture species.
sampling of pteria was nearly comprehensive; it included  <dig> exemplars of  <dig> nominal species  of estimated total of  <dig> species.
the genus isognomon was represented by  <dig> exemplars corresponding to potentially  <dig> of approximately  <dig> species total, and includes the type species, i. isognomum.
the genus malleus was represented by  <dig> exemplars, corresponding to  <dig>  species, including malleus malleus, the type species of the genus, and malleus regula, a type species of the subgenus malvufundus.
it must be pointed out that the true global species diversity of isognomon and malleus is not known due to potential high endemism and cryptic diversity associated with adaptations to hidden habitats, such as interstitial crevices and submarine caves.
the extent of species coverage for the relatively rare and inadequately studied genera vulsella, crenatula, and electroma is uncertain due to cryptic lifestyle in coral crevices and inside sponges.
a broad survey based on dry shells from over twenty major museum collections worldwide suggested that these genera contain very few species and that crenatula is likely to be monotypic.
the latter was represented in the analysis by a single specimen corresponding to the type species, c. avicularis .
the genus electroma was represented by seven exemplars corresponding to potentially four species, including e. alacorvi and electroma  physoides , the sole living member of the subgenus.
two exemplars of vulsella, corresponding to the type species v. vulsella, one of possibly only two living species, were included in the analysis.
the genus pulvinites was represented by one exemplar of the only extant pulvinitid species, p. exempla.
molecular data
locus selection
genome-wide sampling of multiple independently evolving genes is instrumental in overcoming the incongruence among trees derived from individual gene analyses, and producing a resolved and strongly supported phylogeny.
because this study involves the analysis of taxa of disparate levels of divergence, a data set was assembled that consisted of nuclear loci  and a mitochondrial ribosomal rna locus, that in combination span a broad spectrum of variation.
the ribosomal dna sequences display a substantial difference in evolutionary rates among and within the genes , and are easily amplifiable due to the presence of multiple copies per genome and the pattern of concerted evolution  <cit> .
the following dna regions were used:  the complete sequence of the small nuclear ribosomal subunit ;  the d1-d <dig> fragments of the large nuclear ribosomal subunit ;  a fragment of the large mitochondrial ribosomal subunit , and  the single-copy nuclear histone gene h <dig> , which is intron-free and highly conserved at the level of amino acids.
the total length of the alignment of the complete data set was  <dig> bp.
dna isolation and amplification
original dna sequence data have been obtained for most sampled terminals; remaining sequences were obtained from the genbank database .
dna was extracted from fresh, frozen, dried, ethanol-  or rnalater®-preserved tissues.
total genomic dna was isolated from adductor muscle, foot, or the entire body using dneasy® tissue kit .
dna amplifications were carried out using geneamp® pcr system  <dig>  using standard procedures.
the pcr primers are listed in table  <dig>  amplifications were carried out in 35- <dig> cycles using the following general temperature profile: initial denaturation for  <dig> min at 94°c, denaturation for 30- <dig> s at 94°c, annealing for 15- <dig> s at 40-65°c, and extension for 30- <dig> s at 72°c, with a final extension at 72°c for  <dig> min.
specific conditions varied across loci and taxa.
the pcr products were subsequently desalted and concentrated using an arrayit® pcr product purification kit  on a biomek®  <dig> laboratory automation workstation , and sequenced in both directions using the abi prism® bigdye™ terminator cycle sequencing reaction kit  on abi prism® 3700/3730xl dna analyzers  at amnh.
sequence fragments were assembled into contigs, and checked for errors and ambiguities against chromatograms using sequencher™  <dig>  .
the sequence coverage varied from 2× to 6× depending on the quality of the material.
the resultant sequences were checked for potential contamination by blast searches  <cit>  as implemented by the ncbi website http://www.ncbi.nlm.nih.gov.
the new sequences have been deposited in genbank under accession numbers hq329247-hq <dig> .
dna sequence homology
the ribosomal dna sequences, the principal source of character data in the present study, vary considerably in length, creating uncertainty in the inference of nucleotide-level homology.
to circumvent alignment ambiguity, all potential nucleotide homologies were evaluated in the framework of dynamic homology  <cit> , as implemented in poy  <dig>   <cit> .
in this approach nucleotide correspondences  and their transformations  are simultaneously optimized for each tree topology provided an optimality criterion and a cost function .
the integration of alignment and tree searching into a single procedure to produce globally optimal trees allows for evaluating multiple unique homology schemes on a consistent theoretical basis, as both are evaluated using the same criteria  <cit> .
in simultaneous analysis using direct optimization, morphological and molecular characters are co-optimized at the level of sequence alignment in search of a globally optimal solution.
the homology scheme corresponding to the topology of the optimal cladogram was used to produce an implied alignment  <cit> .
the implied alignment was used for statistical tests of sequence variation and other analyses requiring fixed character matrices with equal length character strings .
dna sequences were partitioned into fragments to reduce computation time and to constrain the assignment of insertions and deletions  to unambiguously homologous regions.
none of the regions were excluded from the analyses except for the flanking primer sequences.
initially, multiple alignments were constructed using clustalx  <dig> .
<dig>  <cit>  under default parameter settings.
for rdna sequences, the alignments were subsequently partitioned into fragments flanked by highly conservative  regions corresponding to elements of secondary structure.
the secondary structure features were identified with reference to those inferred for molluscan  <cit>  and other metazoan  taxa.
lastly, the gaps were removed resulting in a total of  <dig> matrices of unaligned fragments.
all histone h <dig> sequences had the same length.
a comparison of the translated h <dig> open reading frame with the published database of histone sequences  <cit>  confirmed high conservation at the level of amino acids.
in this case, the nucleotide homology was unambiguous and the h <dig> data set was treated as a single aligned fragment.
patterns of nucleotide substitutions
individual data partitions were tested for substitution saturation using a non-parametric statistical test based on an information entropy index  <cit>  implemented in dambe  <dig> .
<dig>  <cit> .
the actual number of transitions and transversions, and their ratio were obtained by optimizing substitutions on the optimal topology using macclade  <dig>   <cit> .
departure from homogeneity of base composition across all taxa was assessed using the χ <dig> test as implemented in paup*  <dig> d <dig>  <cit> .
the extent of rate variation across sites for individual data partitions and for the entire data set was estimated by the shape parameter α of the gamma distribution.
the values of α were estimated by the maximum likelihood method assuming the discrete gamma model  <cit>  with four rate categories as implemented in paup*  <dig> d <dig>  <cit> .
for close inspection of species-level divergence of the problematic pinctada imbricata/fucata/radiata species complex, the genetic distances were obtained for within and among presumed species, and then compared to the inter- and intraspecific levels of sequence divergence inferred for their congeners.
in this analysis, the distance matrices for rdna data sets were generated using the minimum evolution objective function  under the general time-reversible model of nucleotide substitution  using paup*  <dig> d <dig>  <cit> .
the choice of the substitution model was guided by the best-fit model preferred by the akaike information criterion  for the combined molecular data set under maximum likelihood  optimality criterion .
the h <dig> data set was not used in this analysis due to missing data for several critical taxa in question.
morphological data
a previously published data set based on morphological data  <cit>  was used for simultaneous analysis of morphological and molecular data.
the morphological characters and methods have been described and illustrated in detail in  <cit> .
the data set included  <dig> morphological characters scored for  <dig> exemplar species, representing all valid pterioidean genera.
as in the original study, all characters were unweighted and one multistate character  was treated as additive based on unequivocal ontogenetic evidence.
molecular data were available for all the taxa for which morphological characters were sampled with the exception of malleus anatinus.
therefore, for the simultaneous analysis of morphological and molecular data sets, the sequences of a closely related species, m. albus, were used as a surrogate for the missing molecular data for m. anatinus.
given the small number of taxa for which morphological data was available, the benefits of minimizing the amount of missing data at the expense of using a composite taxon was justified in the context of the present analysis aimed primarily at resolving supraspecific relationships.
for the morphological character matrix, see  <cit> .
character matrices and incongruence-length difference tests
several different molecular data sets varying in taxon and character sampling were examined in the present study.
the most comprehensive character matrix  included sequences for four loci scored for  <dig> exemplars .
the complete data set contained representatives of all supraspecific taxa and approximately 92% of all extant species of the pterioidea.
in this data set 43% terminals were represented by sequences for all four loci, 54%-by three loci, 2%-by two loci, and only 1%-by a short fragment of a single locus.
thus, character sampling was not satisfactory for 3% of terminal taxa , and their placement must therefore be treated with caution.
the second character matrix  was restricted to the subset of those taxa for which morphological data was available.
the reduced data set was tailored for a combined analysis to have a maximum taxonomic overlap of morphological and molecular character matrices with minimum missing data.
this data set contained  <dig> exemplars , representing all pterioidean genera and families.
morphological data was scored for all  <dig> terminals;  <dig> of these terminals had a complete complement of molecular data ,  <dig> terminals had sequence data for three loci, and one terminal  was represented by a single fragment of the 18s locus.
the simultaneous analyses of molecular and morphological data were performed under maximum parsimony optimality criterion with the same cost for transformations of morphological characters as for nucleotide substitutions.
in addition to the complete and reduced matrices of combined molecular data, five data sets corresponding to the combined nuclear data  and the four individual loci  were analyzed in the parsimony framework to explore topological differences, and the extent of resolution and support among the partitions.
to evaluate the degree of incongruence between the data sets quantitatively, the incongruence-length difference tests  were performed in paup*  <dig> d <dig>  <cit>  .
the ild tests were done simultaneously on the four individual locus matrices as well as pairwise on the complete data set.
the tests were performed under parsimony using  <dig> replicates, where each replicate consisted of  <dig> random-addition sequences , with one tree retained per iteration, followed by tree bisection and reconnection  branch swapping to completion.
phylogenetic analyses
tree searching and optimization
phylogenetic analyses of dna sequence data were performed using direct optimization  <cit> , a heuristic approximation to the optimal tree alignment method  <cit> , under parsimony optimality criterion.
a heuristic extension of generalized optimization  <cit> , direct optimization includes insertion and deletion events in addition to substitutions as character transformations, allowing for the analysis of sequences of different lengths.
additional analyses based on a static alignment and ml optimality criterion were performed to characterize the patterns of nucleotide substitutions and to explore potential sources of systematic error for inferences based on the assumptions of parsimony.
these analyses were based on the implied alignment corresponding to the single optimal cladogram obtained from the parsimony analysis of the combined data performed by direct optimization under the equal-cost regime for the combined molecular character data set.
parsimony analyses under direct optimization, as implemented in poy  <dig>   <cit> , were carried out in parallel on the american museum of natural history parallel computing cluster and a dual-processor personal computer.
a stepwise tree searching strategy that combines independent multiple starting points and efficient algorithms designed for escaping local optima was used to insure sufficient sampling of tree space.
initially, 250- <dig> wagner trees  <cit>  were generated by ras and submitted to branch swapping to completion by alternating rounds of subtree pruning and regrafting  and tbr.
the resultant optimal cladograms were subjected to  <dig> iterations of tbr parsimony ratchet  <cit> , upweighting 20% of nucleotide characters by a factor of  <dig> and retaining one optimal cladogram per iteration.
in the next step, the optimal cladograms were submitted to tree fusing  <cit> , allowing for up to  <dig> pairwise clade exchanges, which was followed by another round of tbr.
finally, the resultant optimal cladograms were submitted to a final round of tbr under exact three-dimensional direct optimization , a more efficient, but substantially more computationally intensive method.
the final equally most-parsimonious trees  were rediagnosed to ensure the correctness of tree length calculation and were used to infer implied alignment.
all tree searches were performed with unconstrained ingroup and outgroup designation  <cit> .
the ml analyses were performed using raxml  <dig> .
<dig>  <cit>  for  <dig> independent replicates.
the results were evaluated for topological congruence and convergence on the optimal likelihood value to ensure that the heuristic searches of tree space were sufficiently thorough.
there appears to be no starting-point dependence to successive approximation that could potentially compromise the accuracy of approximation compared to full-optimization of ml searches  <cit> .
the best-fit model of substitution, given the optimal cladogram and the corresponding implied alignment, was selected by modeltest  <dig>   <cit>  for the combined molecular data set using the akaike information criterion .
the model favored by the akaike information criterion  corresponded to the general time-reversible model that included estimation of the proportion of invariant sites  and assumed a gamma distributed rate variation among sites .
however, due to non-independence of p-invar and the shape parameter α of the gamma distribution  <cit> , the ml analyses were performed under the gtr + Γ model.
topological differences between the trees based on the partitioned and combined data do not present a problem for choosing among alternative models of sequence evolution, because tree topology does not strongly affect model estimation unless long internal branches are involved  <cit> .
the ensemble consistency index  and the ensemble retention index  for the optimal tree were calculated in poy  <dig>   <cit>  with the dynamic-homology characters transformed into static nucleotide-level characters by static approximation  <cit> .
the robustness of phylogenetic relationships was evaluated using the bremer  and jackknife  support methods at nucleotide level.
relative contribution of individual genes was assessed by partitioned bremer support .
the batch files for constrained searches were generated by treerot  <dig>  <cit> .
the constrained searches  were performed by paup*  <dig> d <dig>  <cit> .
for calculating jackknife support values,  <dig> resampling iterations with 36% of sites removed during each pseudoreplicate were performed using poy  <dig>   <cit> .
for the partitioned analyses, only jackknife support values were estimated.
the resulting trees were visualized using poy  <dig>   <cit>  and figtree  <dig> .
<dig>  <cit> .
the optimization and ancestral state reconstruction of morphological characters were evaluated using macclade  <dig>   <cit> .
the significance of the length difference between alternative topologies was tested by the kishino-hasegawa test using resampling of estimated log-likelihoods  bootstrap  with  <dig> replicates under ml criterion, and the wilcoxon signed-ranks test  and the kh test under parsimony criterion as implemented in paup*  <dig> d <dig>  <cit> .
the tests conducted under ml were performed using the gtr + Γ, the same model used for the ml tree searches.
the likelihood-based tests are generally considered inapplicable for post-hoc tests of the optimal versus an alternative topology  <cit> .
in the present study, these paired-sites tests were used to assess the significance of alternative outgroup taxon configurations using constrained topologies defined a priori  for the reduced molecular data set.
indel treatment and sensitivity analysis
because direct optimization integrates the implied alignment and the explicit topology, the choice of analytical parameters  is critical for phylogenetic inference, as analyses performed under alternative cost matrices can produce widely different hypotheses of relationships.
to evaluate the effect of varying alignment parameter values on cladogram topology, independent analyses of the complete data set were performed for a range of parameter combinations  under parsimony.
sensitivity analysis allows discerning between stable relationships  and unstable relationships .
twelve parameter combinations  were used in the analysis: indel/substitution ratios of  <dig>   <dig>  and 4; and tv/ts ratios of  <dig>   <dig>   <dig>  and ∞.
when the tv/ts cost ratio was other than one, the insertion-deletion cost was set relative to the cost of transversions.
transversions were never given a higher cost than indels.
to account for a likely possibility that insertions and deletions spanning multiple contiguous nucleotide positions could have been accommodated by single events, additional sets of parameters taking into account affine indel  costs  <cit>  were carried out.
these analyses explored the indel-substitution ratios of  <dig>   <dig>  and  <dig>  and tv/ts ratios of  <dig>   <dig>   <dig>  where the extension gap cost was set at a constant minimum .
this corresponded to five ratios of gap opening/gap extension costs:  <dig>   <dig>   <dig>   <dig>  and  <dig>  in addition, a cost regime that maximizes homology of both sequence fragments and individual nucleotide positions in the analysis of sequences of variable length  <cit>  was also examined.
this weighting scheme, based on the view of parsimony as a two-taxon analysis  <cit> , corresponds to the following combination of costs:  <dig> for substitutions,  <dig> for gap opening, and  <dig> for gap extension  <cit> .
in total,  <dig> parameter combinations were analyzed.
the results of the sensitivity analysis are shown as binary plots  to designate presence or absence of every clade recovered in the analysis of the complete data set under equal-costs regime.
in addition, these results are summarized as a stability tree-a consensus  of clades recovered under the entire parameter space explored.
despite the fact that character and topological congruence, as measured by the ild  <cit>  and the topological ild  respectively, has been proposed as a criterion for choosing among alternative weighting schemes, these measures are partition-dependent and their values are not comparable across different data matrices  <cit> .
because of the lack of such an objective criterion to decide among alternative weighting schemes, results from equal costs were preferred.
results under additional parameter value combinations were used to assess clade stability.
RESULTS
dna sequence divergence
a detailed description of dna sequence divergence patterns is provided in additional file 1: dna sequence divergence.
as expected, the average levels and the range of divergence, based on the estimates of p distance , increased with relative taxonomic rank and, when compared across loci, were the smallest for the 18s and largest for the 16s sequences.
with regard to the taxonomically problematic pinctada imbricata/fucata/radiata species complex, even though none of the markers could provide unambiguous cut-off levels for intra- vs. interspecific sequence divergence, in all cases the levels of divergence in pairwise comparisons between p. imbricata, fucata, and radiata exemplars were invariably lower compared to other interspecific comparisons, and were the same or marginally exceeding the limits of the intraspecific sequence divergence range inferred for other species within pinctada.
nucleotide base composition
overall base composition was homogeneous across all taxa: a χ <dig> test rejected the null hypothesis of base-composition stationarity neither for any individual locus, nor for the complete data set .
the 16s sequences had a low  proportion of cytosine and the prevalence of a + t , that has been previously reported for other molluscan taxa  <cit> .
across all loci, the sequences were to some extent gc-rich .
<dig>
the average and range  of base composition  for each data partition across taxa.
the results of the χ <dig> value  and the p value  of the χ <dig> test  are listed in the rightmost column.
patterns of nucleotide substitution
the entropy-based statistical test  <cit>  of the combined molecular data set showed no significant net effect of saturation .
it revealed that the 18s and h <dig> sequences were unaffected by multiple hits , but suggested the presence of substantial substitution saturation in the 16s sequence .
saturation of the 28s sequences can potentially become problematic when extremely asymmetrical topologies are considered , which are probably unrealistic given the results of the phylogenetic analysis .
the overall value of the shape parameter α of the gamma distribution was  <dig> , suggesting extensive heterogeneity in substitution rates across sites.
there was a considerable disparity in the levels of α across loci: the highest level of across-site variation was found in 18s , followed by h <dig> , 28s , and 16s .
it is noteworthy, that all the nuclear loci displayed substantial  variation, whereas the mitochondrial locus showed relatively moderate  variation in substitution rates among sites.
the estimates of α appear to be accurate for such range of α   <cit> .
in h <dig> sequences it was likely due to codon bias: virtually no substitutions were inferred for the second position, less than 2% of sites  experienced substitutions in the first position, whereas up to roughly 20%  of sites in the third position underwent substitution.
incongruence-length difference tests
the ild tests performed simultaneously on all four molecular data partitions showed that the complete data set was homogeneous , suggesting that combining molecular partitions in a phylogenetic analysis was not likely to reduce phylogenetic accuracy.
the pairwise ild tests, revealed, however, that only the rdna partitions were entirely congruent , whereas the h <dig> data set was congruent with 16s , but not with the nuclear rdna loci .
despite the fact that the appropriate threshold of incongruence is not known, it has been argued that a significance threshold of  <dig>  might be too conservative for the ild test  <cit> .
given the overall homogeneity of the combined data and marginal incongruence of h <dig> data, the data-set incongruence was not likely to render an analysis of the combined data misleading.
consequently, the phylogenetic inference based on the combined data was the preferred description of pterioidean relationships, but one must be cautious when interpreting the results with regard to the contribution of the h <dig> data.
the ild tests together with phylogenentic analyses of individual data partitions and the assessment of their relative contribution to support on a combined-data cladogram  were used to assess the distribution, nature, and extent of agreement among data sets.
combined analyses of molecular data
a parsimony analysis of the complete molecular data set performed under equal-costs regime produced a single, well-resolved optimal cladogram .
the superfamily pterioidea was recovered as a monophyletic, well-supported clade , which was stable under most  combinations of alignment parameters.
the few cases of non-monophyly of the superfamily were due entirely to the placement of pulvinites exempla among the outgroup taxa as a sister group to the ostreoidea .
this arrangement was recovered only under parameter combinations characterized by the highest indel/substitution ratio, suggesting that the spurious placement was due to missing data, an expected artifact given the fact that p. exempla was the only taxon for which only a small fraction of sequence data was available.
excluding p. exempla from the analysis did not affect the relationships of the remaining taxa , suggesting that the effect of the missing data for this taxon on the relationships among remaining ingroup taxa was negligible.
the suprageneric relationships were completely resolved and all polytomies were restricted to a single derived subclade within the isognomon clade.
all traditionally recognized families-with the obvious exception of the monotypic pulvinitidae-were recovered polyphyletic, whereas all genera  were definitively monophyletic.
the monophyly of the genera pulvinites and crenatula could not be tested in the present study because these genera were represented by single exemplars.
thus, the polyphyly of the families was a result of unconventional relative placement of traditionally recognized genera.
the genus-level clades were distributed in a fully pectinate pattern of relationships with the basal-most pulvinites and the others branching off in the following sequence: malleus, pteria, isognomon, pinctada, vulsella, and electroma .
the nodes along the backbone had relatively high support  and were stable under the majority of alignment parameters  with the exception of the vulsella/electroma/crenatula group.
despite the high support values  obtained under the equal-costs regime, the monophyly of this clade was recovered in only 41% of alignment parameter combinations.
the sister-group relationship of crenatula and electroma  was robust  and insensitive to parameter choice.
the monophyly of all the genera was very strongly supported  and entirely stable across alignment parameter space .
the monophyly of higher taxa among the outgroups was strongly supported  in all analyses and was insensitive to the choice of alignment parameters.
the ostreoidea was identified as an immediate outgroup to the pterioidea with high support  under equal-cost regime, but the pterioidea/ostreoidea sister-group relationship was highly sensitive to alignment ambiguity and was obtained for only four  combinations of alignment parameters.
effectively, all possible arrangements of the considered superfamilies were recovered .
the cladogram generated under the parameter set that maximizes homology of both sequence fragments and base correspondences , had a similar overall structure to the cladogram generated under uniform weighting with several important differences:  the placement of isognomon and pteria clades were exchanged;  the vulsella/electroma/crenatula clade was not recovered; and  the pinnoidea was an immediate pterioidean outgroup.
relative contributions of different data partitions, calculated as partitioned bremer support values for each node, are shown in figure  <dig>  each locus provides support for some deep and apical nodes, but not all of them.
there also appears to be a slight taxonomic bias in the distribution of the pbrs values: for example, relatively high pbrs values for most nodes within pteria were provided by 28s, whereas comparable in magnitude pbrs values for the majority of the nodes within pinctada were provided by 16s.
taken together, 18s and 28s bring close to 75% of the signal and contribute the highest on average support values; h <dig> provides the least support  and the lowest on average support values.
when the support of each gene was normalized by the number of informative characters, 18s appeared as most informative, providing more than twice the support as 28s, whereas 16s and h <dig> were the least informative markers in this data set.
a ml analysis of the complete molecular data set produced in a single optimal tree .
the values of suboptimal trees from  <dig> % of search replicates were within  <dig> % of the optimal likelihood value;  <dig> of these had identical to the optimal likelihood values to the third decimal point and differed only in the placement of several conspecific taxa.
the estimated shape parameter α of gamma distribution was  <dig> ; the substitution relative rate parameters for the model were  <dig>  ,  <dig>  ,  <dig>  ,  <dig>  ,  <dig>  ,  <dig>  .
comparison of the optimal parsimony and ml topologies showed non-substantial differences with regard to placement of several species-level taxa within pteria and malleus, and alternate arrangement of conspecific individuals of pinctada imbricata.
a parsimony analysis of the reduced molecular data set performed under equal-costs regime resulted in a single well-supported  cladogram .
in spite of a great difference in taxon sampling between the complete and reduced data sets, the higher-level relationships were nearly identical with only two notable differences: in the analysis of the reduced data set pulvinites exempla emerged as a sister taxon to the outgroup crassostrea virginica  and a pinnid, atrina rigida, was recovered as an immediate outgroup to the pterioidea.
as suggested by the analyses of the complete data set, missing data and sensitivity to alignment ambiguity coupled with the limited taxon sampling of the reduced data set probably account respectively for the instability in the position of p. exempla and the placement of the immediate outgroup.
the phylogenetic analyses of individual molecular data partitions resulted in cladograms that differ in topology, degree of resolution, and relative support measures.
however, all the analyses invariably recovered the pterioidea and major genus-level clades monophyletic.
combined analyses of morphological and molecular data
a parsimony analysis of the combined morphological and the reduced molecular data sets  performed under equal-costs regime produced a single, highly supported tree .
when compared to the tree generated for the reduced  data set analyzed under the same conditions, the addition of morphological data reconstituted monophyletic the pterioidea by placing pulvinites exempla at the base of the pterioidean subtree; the placement of the remaining taxa was identical.
similarly, the arrangements of pterioidean genera in the analyses of the combined reduced and complete molecular data sets were identical with the only significant topological difference being the identity of the immediate pterioidean sister-group: the pinnoidea in the former and the ostreoidea in the latter.
the paired-sites tests, applied to three topologies constrained to the identical arrangement of pterioidean taxa , but differing solely in the placement of the outgroup taxa, showed that the compared trees were not significantly different .
only the monophyly of the ostreoidea-pinnoidea clade could be rejected in one set of tests under parsimony .
the significance of the length difference between alternative hypotheses of relationships  were evaluated using the kishino-hasegawa  test under parsimony  and maximum likelihood , and the wilcoxon signed-ranks test .
* = significant difference, p <  <dig> .
other abbreviations: ost, ostreoidea; pin, pinnoidea; pte, pterioidea.
a parsimony analysis of the combined morphological and the complete molecular data sets  produced a well-resolved and strongly supported single optimal cladogram .
the topology was essentially the same as of the tree obtained under the same conditions for the complete molecular-only data set: the higher-level relationships along the backbone were identical.
the few differences included an alternative placement of several clades within pteria  and the rearrangement of several, probably conspecific, representatives of isognomon.
the most significant discrepancy between the analyses was the placement of electroma papilionacea.
in the combined analysis it was recovered at the base of the ) clade, whereas in the molecular-only study it was nested within it: )).
discussion
reliability of phylogenetic inference
the present study provides the first explicit hypothesis for higher-level phylogeny and species-level relationships for the bivalve superfamily pterioidea based on dna sequence data.
nearly complete species-level taxon and character coverage suggest that the random error associated with sampling was not likely to compromise the fidelity of the results inasmuch as supraspecific relationships are concerned.
uncertainty in species-level relationships does, however, remain, due to incomplete taxon sampling, differences in character sampling among congeneric individuals, and potentially insufficient phylogenetic signal of chosen markers for resolving recent divergences.
sequence alignment was straightforward and unambiguous for the protein-coding gene.
for the ribosomal sequences, the effect of the uncertainty in nucleotide positional homology on stability of inferred relationships was evaluated using the sensitivity analysis, which allowed for recognizing specific parts of the phylogeny affected by the uncertainty in nucleotide homology and evaluating the extent of its effect.
nucleotide base composition was nearly perfectly stationary across taxa overall as well as for individual loci, without evidence of taxon attraction due to similar nucleotide-composition bias .
substitution saturation in dna sequences can potentially underestimate the actual amount of evolutionary change and reduce the phylogenetic signal due to excessive multiple hits in the same nucleotide positions resulting in extensive homoplasy .
despite the fact that extensive saturation was detected for transitions in the 16s sequence, its effect on the phylogenetic inference was negligible, as shown by  comparison of the results from individual data partitions,  incorporating the ts/tv into phylogenetic analysis,  character congruence among data partitions, and  consistency in the distribution of brs values.
furthermore, excluding 16s data from the analysis resulted in a decreased resolution and support throughout the tree.
these results are consistent with theoretical and empirical findings suggesting that homoplasy does not necessarily confound phylogenetic inference  <cit>  and that a given site, which is homoplastic in a sub-part of the most-parsimonious tree, can bring phylogenetic information for another sub-part of the same tree  <cit> .
heterogeneity in substitution rates among sites within a nucleotide sequence can have considerable impact on sequence divergence  <cit> .
substantial heterogeneity in substitution rates among lineages resulting in long-branch attraction is another source of misleading inference in molecular data  <cit> .
the agreement between the results of the maximum likelihood analysis under the model allowing for variable base frequencies and parsimony suggested that the predominance of transitions and among-site heterogeneity for most loci, as well as among-lineage heterogeneity in nucleotide substitution rates  were not likely to be a source of misleading inference.
the absence of long-branch attraction despite considerable differences in branch lengths was consistent with prior observation that extensively sampled phylogenies are not sensitive to this artifact  <cit> .
the doubly uniparental inheritance   <cit> , which has been reported for a number of bivalve species  <cit> , can potentially affect the accuracy of phylogenetic reconstruction due to varying degrees of divergence and homogenization between conspecific genomes.
within the subclass autolamellibranchiata, the instances of dui were reported only from the most basal extant group, mytiloida.
even though the lack of reported instances of dui in a particular group does not necessarily imply that it is not present, the absence of strict dioeciousness in the pterioidea suggests that dui was not a likely to be a problem.
all pterioidean species for which reproductive strategies have been elucidated , are rhythmical protandrous hermaphrodites, with occasional instances of protogyny and simultaneous hermaphroditism, and mtdna appears to be maternally inherited  <cit> .
simultaneous analysis of a single data matrix combining independent character sets potentially increases the descriptive efficiency and explanatory power  <cit> , provides the basis for evaluating the relative support and taking into account the hidden support provided by individual partitions  <cit> , and is less sensitive to sampling bias  <cit> .
however, the incongruence among data partitions may indicate that they have evolved under different models of evolution and render a combined-data analysis misleading  <cit> .
despite topological differences obtained in the analyses of individual data partitions, the results of ild tests showed that the data sets were largely congruent, implying that a single phylogeny underlies the relationships inferred from individual data partitions.
differences in substitution rates across sites, the average pairwise sequence divergence values, the distribution and levels of support , and the extent of resolution obtained in the analyses of individual partitions indicated that the marker loci collectively encompassed sufficiently broad spectrum of dna sequence divergence to produce a nearly completely resolved and well-supported tree.
these results also suggested that topological differences among data partitions were probably due to characters that evolved rapidly to provide a reliable signal for sub-terminal  nodes of a tree, but produced noise for the resolution of deeper  nodes.
the overall increased resolution and support of the combined-data trees when compared to the results of partitioned analyses strongly speak for the simultaneous analysis approach.
in summary, the present molecular sequence data and the choice of methodology adequately account for sampling biases and potentially confounding  effects of various molecular evolutionary processes, providing a strong foundation for accurate, reliable, and robust reconstruction of phylogenetic relationships of the recent pterioidea.
phylogeny of the pterioidea
the present study corroborates the monophyly of the superfamily pterioidea and definitively rejects that of all nominal pterioidean families .
whereas the non-monophyly of some traditional pterioidean groups was anticipated by prior studies  <cit> , the arrangement of the genera was novel and unexpected, as well as resolved and strongly supported.
in terms of taxonomic congruence with morphological studies, the analysis of molecular data offered a considerably different phylogenetic hypothesis.
generally assumed close relatedness of commercially important genera pteria and pinctada was rejected by the molecular evidence, which was not entirely surprising given the absence of unambiguously optimized non-homoplastic characters supporting the sister-group relationship of these taxa  <cit> .
the presumed close relatedness of the genera malleus and isognomon, implied by superficial resemblance in shell morphology and similar life habits  <cit> , was also contradicted by the molecular data.
the only morphological character supporting a close relatedness of these two genera-the absence of food grooves on inner demibranchs of the gills-is symplesiomorphic in the context of the dna based topology, as this condition has been also reported from other autolamellibranchiate taxa .
contrary to previously proposed hypothesis  <cit> , this study rejected a close affinity of isognomon and pulvinites, placing the latter at the base of the pterioidean tree.
this finding was corroborated by the fact that in the morphology-based analysis  the sister-group relationship of the two genera was weakly supported and relied on a unique combination of shared homoplastic characters .
the present analysis provided the first molecular evidence for a peculiar clade, recently recognised on morphological grounds   <cit> , that consists of relatively understudied, low-diversity genera electroma, vulsella, and crenatula, that were traditionally placed in pteriidae, malleidae, and isognomonidae respectively.
the arrangement among the genera, however, varied between the two studies.
in the morphology-only analysis, where each genus was represented by a single exemplar, electroma was basal to vulsella and crenatula, the sister-group relationship between which was implied by the reduction of the pedo-byssal apparatus.
instead, in the molecular analysis, vulsella appeared basal to the electroma clade that nested the sole representative of crenatula, rendering electroma paraphyletic.
an unanticipated discovery of polyphyletic electroma requires detailed morphological and taxonomic reappraisal of this poorly known genus.
another species of electroma, typically placed in its own subgenus pterelectroma, e.  physoides, invariably grouped with the pteria clade, putting into question the taxonomic status of this subgenus.
due to high sensitivity to alignment parameters and taxon sampling, uncertainty remains regarding the identity of the sister group of the pterioidea.
the paired-sites tests, applied to the alternative topologies differing only in the placement of the outgroup taxa, showed that the trees were not significantly different.
consequently, the sister-group of the pterioidea could not be unequivocally inferred in light of the present results.
based on the results of the combined-data analyses under uniform weighting and predominant support , the ostreoidea is tentatively proposed here as the most likely pterioidean sister group.
implications for homology of morphological characters
it has long been acknowledged that extensive homoplasy is a recurrent theme in bivalve phenotypic evolution  <cit>  that can potentially mislead phylogenetic inference based on morphological data alone.
optimizing morphological character state transformations onto topology derived from the analysis of molecular data provides means to discern putative homologies and instances of homoplasy.
the molecular analyses presented here have confirmed a number of suspected convergences identified on the basis of morphology, but also uncovered a number of hitherto unanticipated ones.
the hypotheses of homology of several key characters  were evaluated in light of a robust phylogenetic hypothesis based on molecular evidence.
a duplicated outer fold of the mantle edge  has been considered a derived trait that evolved twice independently in the pterioidea based on the analysis of morphological data  <cit>  and in agreement with presumably "primitive" three-fold arrangement  <cit> .
in contrast, molecular data suggested that the four-fold condition was primitive in the pterioidea, whereas the three-fold condition was a result of a single reversal to the three-fold condition exhibited by the outgroup taxa.
this finding was not entirely unexpected, given the fact that analogous transitions in mantle evolution had occurred in other bivalve lineages  <cit> .
in the analysis of morphological data, homorhabdic ctenidia  evolved once from heterorhabdic ctenidia .
molecular results, on the other hand, implied that the gill grades had evolved iteratively and the ancestral condition could not be unequivocally resolved.
the presence of food grooves in the outer demibranchs showed the same pattern.
the loss of postlarval byssus and posterior pedo-byssal retractor muscles has been previously considered as evidence for sister-group relationship of vulsella and crenatula  <cit> .
both features are likely to be adaptations for living inside sponges, the habit shared by members of the two genera, where physical stabilization is achieved by the attachment to the host sponge.
the relationships inferred from the molecular data invariably grouped crenatula with electroma, a relationship supported by a single unambiguously optimized nonhomoplastic apomorphy: a t-shaped  typhlosole, a longitudinal fold of the intestine.
molecular results were ambiguous with regard to the immediate relatedness of the vulsella and the crenatula/electroma clades due to alignment uncertainty, and suggested a possibility that the evolution of adaptation for living inside sponges could have been two separate events.
this supposition is corroborated by prior ecological observations that the host specificity, and the nature of mutualistic association of vulsella and crenatula species with their demosponge hosts are different: the former is a facultative mutualist with keratosid sponges, whereas the latter is an obligatory mutualist with monaxonid sponges  <cit> .
species of vulsella possess a relatively thick and elongated shell with an extremely abbreviated hinge line and a narrow prismatic margin that provide a rigid "endo-skeleton" for the sponge, hence facilitating the host's growth.
species of crenatula, on the other hand, have a very thin shell with a very wide prismatic margin, which appears to be influenced by the sponge growth, and produce highly irregular shell shapes without consistent orientation.
the particular form of the multivincular ligament of crenatula allows for curvature of the extensive hinge line, which facilitates attachment along the hinge line to irregular surfaces.
the loss of postlarval byssus in distantly related malleus albus  also supports the likelihood of convergence.
the symmetry and relative placement of abdominal sense organs , paired mechanoreceptors typically found on the surface of the adductor muscle, have been previously shown to be phylogenetically informative within autolamellibranchiata  <cit>  and, specifically, within the pterioidea  <cit> .
in the morphology-based analysis, the clustering of the asos around the anus was inferred as a primitive condition in the pterioidea , whereas the dorsal displacement of the right aso had occurred twice.
according to the results of the molecular character analysis, both aspects of the aso position had evolved twice.
a set of characters that has historically been used for higher-level systematics of pterioidean and other autolamellibranchiate taxa is associated with the types of ligament and their support structures.
according to the molecular-based topology, the multivincular ligament characterizes three distantly related genera: pulvinites, isognomon, and crenatula.
morphological data alone also identified an independent origin of multivincular ligament in crenatula, but suggested a sister-group relationship of pulvinites and isognomon on the basis of shared ligament morphology and two other homoplastic characters.
it is noteworthy that despite apparent similarity between multivincular systems in pulvinites and isognomon, the pattern of ligostracum deposition and the ontogenetic sequence of origination of ligamental layers along the hinge line differ between the genera.
even though not entirely congruent, the results of the molecular and morphological studies reject the long-standing assumption of close relatedness based on a single origin of multivincular ligament in the pterioidea.
the present molecular evidence confirmed the taxonomic significance of the shape of the anal funnel, an ear-shaped membranous structure situated at the posterior extremity of the rectum that presumably facilitates expulsion of fecal pellets.
the absence of the anal funnel in pulvinites is consistent with its basal placement in the pterioidea as the same condition is displayed by the outgroup taxa.
in contrast, in the analysis based on morphological data, the lack of the funnel was inferred as a loss.
the molecular and morphological analyses agree that the anal funnel is subtriangular in outline in the more basal pterioidean clades, whereas it is characteristically lanceolate and rounded in two most derived clades, pinctada and vulsella/crenatula/electroma, respectively.
given the uncertainty in the identity of pterioidean sister group, it is instructive to explore the effect of alternative combinations of outgroup taxa on the distribution of diagnostic morphological characters using the best available hypothesis of the ingroup relationships.
to this end, morphological characters were optimized on three trees constrained to the identical arrangement of the ingroup taxa  but differing in the placement of the outgroups.
the pinnoidea relationship  was defined by five unambiguously optimized, non-homoplastic apomorphic characters that remained unchanged above this node:  the descending intestine produced towards the posteroventral side of the posterior adductor muscle;  the loss of the anterior adductor muscle ;  the absence of the pseudonymphae;  the presence of a resilifer throughout ontogeny; and  the placement of the posterior adductor muscle entirely within the boundary of the inner shell layer.
previously, steiner and hammer  <cit>  had mapped a set of morphological characters on to their maximum likelihood tree substantiating their finding of ostreoidean/pterioidean sister-group relationship by two characters: the monomyary and the presence of the anal funnel.
the present study confirmed the monomyaryan condition as a synapomorphy of the clade, but not the occurrence of the anal funnel in different superfamilies, interpreted here as a result of independent origins, as was previously suggested on the basis of morphological data  <cit> .
the ostreoidea relationship  was defined by three apomorphic characters:  the absence of the marginal mantle lobe fusion;  the absence of the eulaterofrontal cirri of gill filaments; and  the ascending intestine not produced anteriorly passed the stomach.
the pterioidea relationship  was defined by two apomorphic characters:  the absence of the ciliated discs on gill filaments and  the clustering of ducts of the digestive diverticula.
given comparable levels of homoplasy and the lack of a priori criteria to differentially weigh the characters supporting each of the possibilities, the arrangement supported by the largest number of apomorphies and yielding the shortest tree length, that is ), is the preferred hypothesis of relationship.
species-level taxonomy of pinctada
despite the fact that the focus of the present study was to establish higher-level relationships within the pterioidea, it provided significant insights into species-level taxonomy of the genus pinctada, a commercially important and the most exhaustively sampled group.
the present study suggests that pinctada imbricata, p. fucata, and p. radiata are distinct genealogical units.
the taxonomic status of the p. imbricata/fucata/radiata species complex previously remained unsettled due to extensive variation in shell characters within and among populations, extremely wide geographical distribution, and transport and hybridization by humans  <cit> .
traditionally, three distinct species were recognized corresponding to three biogeographic realms: p. imbricata in the western atlantic region, p. radiata in the eastern indian ocean and the red sea regions, and p. fucata in the indo-pacific region.
the japanese populations have frequently been regarded as a distinct species, p. martensii, or a subspecies, p. fucata martensii.
because of its major role in perliculture, the species complex had become a major research focus in recent years.
a large number of comparative genetic analyses using allozyme and dna sequence data indicated significant levels of gene flow over large geographic distances among populations of indian ocean and indo-pacific regions, resulting in low levels of genetic differentiation  <cit> .
the pattern of genetic differentiation among the chinese, japanese, and australian populations using amplified fragment length polymorphism  markers was consistent with isolation by distance  <cit> .
shell morphometric studies  <cit> , mating experiments  <cit> , and comparative karyotype analyses  <cit>  failed to establish unequivocal species delineations.
to date, no diagnostic discrete morphological characters have been identified for populations from the three geographical regions.
the present results showed shallow levels of divergence among the members of the p. imbricata/fucata/radiata species complex, comparable to the levels of divergence of conspecific exemplars of other pinctada species.
however, the three groups were reciprocally monophyletic showing an apparent clinal geographical structure.
therefore, it might be meaningful to treat these individual populations as evolutionary significant units .
taking into account the comparatively low divergence values and the lack of diagnostic morphological characters distinguishing these groups, from the taxonomic standpoint it is advocated here to provisionally recognize them at a subspecies level under the senior synonym, p. imbricata .
the final word on the taxonomic status of these taxa will depend on the future population-level analysis with dense sampling throughout the entire distribution range and using genetic markers more sensitive for resolving recent divergencies.
this study also suggested that pinctada albina and p. nigra, generally distinguished by shell color , represent two colormorphs of the same species.
this possibility has been previously suggested based on extremely low levels of genetic divergence between two exemplars of p. albina and one of p. nigra  <cit> .
this finding was consistent with other molecular studies that invariably resolved these species  as sister taxa  <cit> .
the present analysis contained three exemplars of each colormorph that were always recovered as a monophyletic group and never formed separate clades diagnosed by color.
in addition, based on the fact that p. mazatlanica invariably nested among representatives of p. margaritifera or members of both species formed an unresolved clade, the present study provides preliminary support for conspecificity of these two species.
these findings need to be investigated further and cannot be definitively tested given the present sample size.
CONCLUSIONS
the most significant conclusion from the phylogenetic analyses presented here is that the state of pterioidean systematics is no longer in chaos.
the present results establish a solid-however unconventional- framework for the higher-level taxonomy of the group, focusing future revisionary systematic effort and providing a better grasp of standing alpha-diversity.
the phylogeny is crucial for understanding processes and patterns of diversification and morphological evolution of the pterioidea through time.
to this end, the major future challenge lies in integrating the knowledge on the recent pterioideans with information on the extensive, well-documented fossil record of the group.
as throughout their history, pterioideans remain a globally distributed, relatively common group of bivalves found in most warm, shallow marine environments around the globe, and are likely to display patterns of diversity typical of other marine benthic macrofauna.
therefore, the insights gained from the phylogenetically-informed studies of pterioidean bivalves might be broadly applicable across marine invertebrates.
supplementary material
additional file 1
dna sequence divergence.
word doc file containing dna sequence divergence.
click here for file
additional file 2
partitioned analyses.
word doc file displaying partitioned analyses.
click here for file
acknowledgements
i would like to thank friends and colleagues around the world who donated or loaned specimens: r. bieler , a.j.
baldinger , p. bouchet , j. gerber , g. giribet , l. groves , m.g.
harasewych , d. herbert , k. holmes , japan pearl exporters, m. knapp , k.v.
ramana murty , e. kools , p. maestrati , b. marshall , p.g.
oliver , c.j.m.
lodeiros seijo , r. shearer , l.r.
simone , a. tommasi , k. umadevi , Á. valdés , n. voss , and r.c.
willan .
g. steiner  and t. masaoka  provided information on specimens, sequences of which were deposited in genbank.
i am indebted to e. arizmendi castillo  for providing assistance during my visit of guaymas, mexico in  <dig>  r. rorandelli, f. bartolini, m. fusi, m. vannini , and m. gomei  have provided invaluable help during fieldwork in the tyrranean sea in summer  <dig>  i would like to thank all participants of the second international marine bivalve workshop  and c. printrakoon  in particular for partnership during the workshop and continuing collaboration.
r. desalle  generously made available molecular laboratory resources at the sackler institute for comparative genomics.
w.c. wheeler  granted access and logistical support for using the parallel computing cluster.
m. leslie  provided invaluable help for the laboratory component of the project.
m. coppolino and s. watson assisted with the american museum of natural history mollusk collection.
y. iwasaki  had assisted with translations from japanese.
f. almeida, a. varón, m. weksler , and x. xia  made technical suggestions.
p.m. mikkelsen , e. strong , j. hammock , and three anonymous reviewers provided critical suggestions for improving the manuscript.
partial funding toward completion of this work was provided by nsf-peet marine bivalve project deb- <dig> , nsf dissertation improvement grant deb- <dig>  amnh annette kade graduate fellowship, the conchologists of america grant, and nsf-atol bivalve tree of life project deb-0732854/0732903/ <dig> .

 since the lower middle ordovician some million years ago, pterioidean bivalves have inhabited a remarkable diversity of marine epifaunal and semi-infaunal environments around the globe, typically confined to cryptic habitats and forming byssal attachments to various substrata. for a much shorter part of their history-from the dawn of human culture-they have become the primary source of pearls and nacre. pearl fishing based on natural populations of pterioidean bivalves had expanded into a growing global industry with a current value of approximately half a billion us dollars per annum. in pursuit of pearls, commercial introductions and accidental transport of pterioideans beyond their native distribution ranges have greatly affected population dynamics of some species endangering local indigenous biotas. a recent interest in physical properties of mother-of-pearl made pterioidean bivalves a model system for elucidating molecular mechanisms of biomineralization with medical applications for bone regeneration. given the long and many-sided history of pterioideans and humans, and the current economic and ecological significance of these bivalves, surprisingly little is known about the standing alpha-diversity, distribution, and evolutionary history of the group. the present-day diversity of the superfamily pterioidea gray, encompasses four families traditionally defined by shell shape and ligament structure: pteriidae, isognomonidae, malleidae, and pulvinitidae. pteriids are distinguished by an obliquely ovate shell shape, an enlarged posteriorly projecting auricle, and a deep byssal notch. many pteriid species, particularly of genera pteria and electroma, evolved specialized associations with hydroid, scleractinian, and alcyonarian substrata. most species of pinctada inhabit sandy and hard bottoms, some living in association with commensal fishes and crustaceans. several species of pteria and pinctada are cultured for commercial pearl production. most species of a malleid genus malleus are irregularly shaped to conform to narrow crevices in hard coral and rocky substrata, and develop elongated extensions of the hinge for stabilization in soft sediment. species of the second malleid genus, vulsella, are adapted for living within sponges. isognomonids are characterized by the distinct morphology of the ligament that contains multiple grooves for ligamental attachments arranged sequentially along the hinge line. species of the genus isognomon often co-occur with species of malleus in crevices and on soft muddy bottoms, but also are gregarious on mangrove roots, whereas species of crenatula, the second isognomonid genus, live inside sponges. similarly to isognomonids, pulvinitids possess multivincular ligament, but are distinguished from isognomonids by the presence of the foramen, an opening in the right valve through which the byssus is protruded. living pulvinitids are represented by the genus pulvinites, the sole extant species of which, pulvinites exempla, lives byssally attached to vertical hard substrata at 200- m depths. previous studies aimed at resolving higher-level phylogeny within bivalvia agreed on the monophyly of the pterioidea and its placement within the subclass autolamellibranchiata. however, most aspects of relationships within the superfamily and the identity of its immediate sister group remain uncertain. in a pioneering study on the evolution and ontogeny of the bivalve shell, jackson proposed the derivation of major pterioidean groups from a pteria¬-like paleozoic ancestor. in his scheme, several autolamellibranchiate lineages, including common oysters of the family ostreidae and several extinct families, evolved from isognomon, the latter ultimately arising from the pteria stem lineage. jackson also provisionally recognized close relatedness of malleus and vulsella. despite the implication of jackson's work for pterioidean systematics, the notion of the monophyly and taxonomic composition of pterioidean families persisted for over a century. molecular studies aimed at establishing relationships within the subclass autolamellibranchiata or the entire bivalvia suggested that several widely accepted pterioidean taxa possibly did not constitute natural groups. the minimum evolution tree based on the analysis of the mitochondrial cytochrome oxydase i sequences did not resolve the relationships among pteria, pinctada, and isognomon. thus, this analysis failed to provide support for the monophyly of pteriidae but supported the monophyly of pinctada and isognomon. a study with more extensive sampling and based on parsimony analysis of sequence data for the nuclear small ribosomal subunit likewise failed to recover the monophyly of pteriidae, but also that of malleidae. in that study, isognomon was monophyletic, pteria was paraphyletic, and pinctada was polyphyletic. the pteriid electroma was a sister taxon to the malleid vulsella. the maximum likelihood analysis of the same data produced a similar, but more resolved, topology recovering monophyletic pinctada. in a more extensive in character-but not taxonomic-sampling analysis of combined 18s and the large ribosomal subunit sequence data, the families pteriidae and malleidae, and the genus pteria were non-monophyletic. when these data were combined with the nuclear histone h gene and the mitochondrial coi gene sequences, the analysis produced a polytomy for all representatives of the pterioidea except for the isognomon and electroma/vulsella clades. in summary, molecular studies, varying greatly in sampling and analytical techniques, produced inconsistent hypotheses of pterioidean relationships with low resolution and support. these findings, however, agreed on the monophyly of the superfamily and suggested that the families pteriidae and malleidae were not monophyletic. the results were conflicting with regard to the monophyly of the genera pteria and pinctada, whereas isognomon was invariably monophyletic. the monophyly of the family isognomonidae and the genera electroma and malleus were not tested, and representatives of the family pulvinitidae and the genus crenatula were not included in the analyses. the only study that focused exclusively on the pterioidea was based on morphological characters and included representatives of all four families and nine genera. it supported the monophyly of the pterioidea and non-monophyly of three families. all the genera that were represented by multiple species were recovered monophyletic. species-level sampling in that study was insufficient to test the assumption of monophyly of the remaining genera and did not include type species for most nominal supraspecific taxa. despite these recent challenges on both molecular and morphological grounds, no consensus on pterioidean relationships was reached and no revision of the existing classification was proposed. the identity of the pterioidean sister group remains uncertain, favoring the oysters, pen shells, or a clade composed of these two groups, currently classified together with the pterioidea in the order pterioida. most analyses based on morphological data for extinct, extant, or both kinds of taxa, have suggested a sister relationship of the pterioidea and pinnoidea, except for one, which instead suggested a closer relationship of the pterioidea and ostreoidea. likewise, molecular analyses were divided in favoring either pinnoidea or ostreoidea as a sister group of the pterioidea. in other studies, either pinnoidea and ostreoidea formed a clade sister to the pterioidea, or the relationships among the three superfamilies were unresolved. the only explicitly combined analysis of morphological and molecular data resolved the pinnoidea basal to the pterioidea/ostreoidea clade. the unsettled question of the pterioidean sister group is pertinent to the recent debates on the origin of oysters in paleontological literature. the present analysis is the first study to produce a nearly comprehensive species-level phylogeny for the pterioidea based on dna sequence data. the principal specific goals of this study include testing the monophyly of the pterioidea, and all its families and genera; reviewing the consequences of these results for the homology of important morphological characters; and resolving species-level phylogenetic relationships within the commercially important pearl oyster genus pinctada. methods taxonomic sampling sources of specimens most specimens were obtained from museum collections worldwide, with additional material gathered during fieldwork in the western atlantic, the eastern pacific, the mediterranean, and the indo-pacific regions. living animals were collected in mangroves and shallow intertidal waters, fixed in 85-95% ethanol, and subsequently transferred to 75-95% ethanol. collected specimens were deposited at the american museum of natural history and the field museum of natural history. dna sequence data were obtained for of the terminal taxa. the remaining sequences were acquired from genbank http://www.ncbi.nlm.nih.gov/genbank/index.html. collection localities and museum catalog numbers for the specimens used in this study are listed in table museum collection acronyms: amnh, american museum of natural history, new york; cas, california academy of sciences, san francisco; fmnh, field museum of natural history, chicago; lacm, los angeles county museum; mnhn, muséum national d'histoire naturelle, paris; nmp, natal museum, pietermaritzburg; usnm, national museum of natural history, smithsonian institution, washington, d. c.; umml, rosenstiel school of marine and atmospheric science, university of miami, florida. na, not available. outgroup taxa given that the pterioidea is currently accepted as monophyletic and the pinnoidea or ostreoidea were established as the taxa most closely related to the pterioidea on morphological and molecular grounds, representatives of these superfamilies were chosen for the outgroup. to ensure maximum diversity of the ostreoidean taxa, the selection of the gryphaeid hyotissa, and the ostreid ostrea and crassostrea exemplars were guided by their placement in a recently proposed phylogenetic hypothesis. to polarize the immediate outgroups, representatives of two other autolamellibranchiate orders, mytiloida and pectinoida were selected based upon the availability of sequence data. mytilioda, represented by the well-known common mussel species, mytilus edulis, had been resolved as the most basal order of the subclass in comprehensive studies based on morphological, molecular, and combined data and, accordingly, was designated for cladogram rooting. ingroup taxa inasmuch as the pterioidea is the primary focus of this study, the taxonomic sampling of this taxon was the most dense. all nominal families and genera of the pterioidea were included in the analyses. the sample contained representatives of type species of all pterioidean genera except for pulvinites, which was described based on an extinct species, p. adansonii, from the late cretaceous. success in securing multiple representatives for each genus varied because the material could not be obtained for few taxa characterized by cryptic life habit and presumed low species diversity. initially, specimens were sorted to morphospecies using anatomical and shell morhological features. species-level identifications were verified with reference to the extant type material and assigned to most specimens. where possible, to avoid assumptions of the extent of intraspecific phenotypic variation, multiple potentially conspecific individuals were included in the analysis. given the lack of recent systematic revisions for the majority of pterioidean taxa, species-level identifications were tentatively chosen with correspondence to their prevailing use in current literature. the sample of the genus pinctada was complete, with exemplars representing nominal species, including the type species pinctada margaritifera. the overrepresentation of pinctada species was necessitated by an enduring controversy regarding the taxonomic status of the pinctada imbricata/fucata/radiata species complex that includes important perliculture species. sampling of pteria was nearly comprehensive; it included exemplars of nominal species of estimated total of species. the genus isognomon was represented by exemplars corresponding to potentially of approximately species total, and includes the type species, i. isognomum. the genus malleus was represented by exemplars, corresponding to species, including malleus malleus, the type species of the genus, and malleus regula, a type species of the subgenus malvufundus. it must be pointed out that the true global species diversity of isognomon and malleus is not known due to potential high endemism and cryptic diversity associated with adaptations to hidden habitats, such as interstitial crevices and submarine caves. the extent of species coverage for the relatively rare and inadequately studied genera vulsella, crenatula, and electroma is uncertain due to cryptic lifestyle in coral crevices and inside sponges. a broad survey based on dry shells from over twenty major museum collections worldwide suggested that these genera contain very few species and that crenatula is likely to be monotypic. the latter was represented in the analysis by a single specimen corresponding to the type species, c. avicularis. the genus electroma was represented by seven exemplars corresponding to potentially four species, including e. alacorvi and electroma physoides, the sole living member of the subgenus. two exemplars of vulsella, corresponding to the type species v. vulsella, one of possibly only two living species, were included in the analysis. the genus pulvinites was represented by one exemplar of the only extant pulvinitid species, p. exempla. molecular data locus selection genome-wide sampling of multiple independently evolving genes is instrumental in overcoming the incongruence among trees derived from individual gene analyses, and producing a resolved and strongly supported phylogeny. because this study involves the analysis of taxa of disparate levels of divergence, a data set was assembled that consisted of nuclear loci and a mitochondrial ribosomal rna locus, that in combination span a broad spectrum of variation. the ribosomal dna sequences display a substantial difference in evolutionary rates among and within the genes, and are easily amplifiable due to the presence of multiple copies per genome and the pattern of concerted evolution. the following dna regions were used: the complete sequence of the small nuclear ribosomal subunit; the d1-d fragments of the large nuclear ribosomal subunit; a fragment of the large mitochondrial ribosomal subunit, and the single-copy nuclear histone gene h, which is intron-free and highly conserved at the level of amino acids. the total length of the alignment of the complete data set was bp. dna isolation and amplification original dna sequence data have been obtained for most sampled terminals; remaining sequences were obtained from the genbank database. dna was extracted from fresh, frozen, dried, ethanol- or rnalater®-preserved tissues. total genomic dna was isolated from adductor muscle, foot, or the entire body using dneasy® tissue kit. dna amplifications were carried out using geneamp® pcr system using standard procedures. the pcr primers are listed in table amplifications were carried out in 35- cycles using the following general temperature profile: initial denaturation for min at 94°c, denaturation for 30- s at 94°c, annealing for 15- s at 40-65°c, and extension for 30- s at 72°c, with a final extension at 72°c for min. specific conditions varied across loci and taxa. the pcr products were subsequently desalted and concentrated using an arrayit® pcr product purification kit on a biomek® laboratory automation workstation, and sequenced in both directions using the abi prism® bigdye™ terminator cycle sequencing reaction kit on abi prism® 3700/3730xl dna analyzers at amnh. sequence fragments were assembled into contigs, and checked for errors and ambiguities against chromatograms using sequencher™. the sequence coverage varied from 2× to 6× depending on the quality of the material. the resultant sequences were checked for potential contamination by blast searches as implemented by the ncbi website http://www.ncbi.nlm.nih.gov. the new sequences have been deposited in genbank under accession numbers hq329247-hq. dna sequence homology the ribosomal dna sequences, the principal source of character data in the present study, vary considerably in length, creating uncertainty in the inference of nucleotide-level homology. to circumvent alignment ambiguity, all potential nucleotide homologies were evaluated in the framework of dynamic homology, as implemented in poy. in this approach nucleotide correspondences and their transformations are simultaneously optimized for each tree topology provided an optimality criterion and a cost function. the integration of alignment and tree searching into a single procedure to produce globally optimal trees allows for evaluating multiple unique homology schemes on a consistent theoretical basis, as both are evaluated using the same criteria. in simultaneous analysis using direct optimization, morphological and molecular characters are co-optimized at the level of sequence alignment in search of a globally optimal solution. the homology scheme corresponding to the topology of the optimal cladogram was used to produce an implied alignment. the implied alignment was used for statistical tests of sequence variation and other analyses requiring fixed character matrices with equal length character strings. dna sequences were partitioned into fragments to reduce computation time and to constrain the assignment of insertions and deletions to unambiguously homologous regions. none of the regions were excluded from the analyses except for the flanking primer sequences. initially, multiple alignments were constructed using clustalx. under default parameter settings. for rdna sequences, the alignments were subsequently partitioned into fragments flanked by highly conservative regions corresponding to elements of secondary structure. the secondary structure features were identified with reference to those inferred for molluscan and other metazoan taxa. lastly, the gaps were removed resulting in a total of matrices of unaligned fragments. all histone h sequences had the same length. a comparison of the translated h open reading frame with the published database of histone sequences confirmed high conservation at the level of amino acids. in this case, the nucleotide homology was unambiguous and the h data set was treated as a single aligned fragment. patterns of nucleotide substitutions individual data partitions were tested for substitution saturation using a non-parametric statistical test based on an information entropy index implemented in dambe.. the actual number of transitions and transversions, and their ratio were obtained by optimizing substitutions on the optimal topology using macclade. departure from homogeneity of base composition across all taxa was assessed using the χ test as implemented in paup* d. the extent of rate variation across sites for individual data partitions and for the entire data set was estimated by the shape parameter α of the gamma distribution. the values of α were estimated by the maximum likelihood method assuming the discrete gamma model with four rate categories as implemented in paup* d. for close inspection of species-level divergence of the problematic pinctada imbricata/fucata/radiata species complex, the genetic distances were obtained for within and among presumed species, and then compared to the inter- and intraspecific levels of sequence divergence inferred for their congeners. in this analysis, the distance matrices for rdna data sets were generated using the minimum evolution objective function under the general time-reversible model of nucleotide substitution using paup* d. the choice of the substitution model was guided by the best-fit model preferred by the akaike information criterion for the combined molecular data set under maximum likelihood optimality criterion. the h data set was not used in this analysis due to missing data for several critical taxa in question. morphological data a previously published data set based on morphological data was used for simultaneous analysis of morphological and molecular data. the morphological characters and methods have been described and illustrated in detail in. the data set included morphological characters scored for exemplar species, representing all valid pterioidean genera. as in the original study, all characters were unweighted and one multistate character was treated as additive based on unequivocal ontogenetic evidence. molecular data were available for all the taxa for which morphological characters were sampled with the exception of malleus anatinus. therefore, for the simultaneous analysis of morphological and molecular data sets, the sequences of a closely related species, m. albus, were used as a surrogate for the missing molecular data for m. anatinus. given the small number of taxa for which morphological data was available, the benefits of minimizing the amount of missing data at the expense of using a composite taxon was justified in the context of the present analysis aimed primarily at resolving supraspecific relationships. for the morphological character matrix, see. character matrices and incongruence-length difference tests several different molecular data sets varying in taxon and character sampling were examined in the present study. the most comprehensive character matrix included sequences for four loci scored for exemplars. the complete data set contained representatives of all supraspecific taxa and approximately 92% of all extant species of the pterioidea. in this data set 43% terminals were represented by sequences for all four loci, 54%-by three loci, 2%-by two loci, and only 1%-by a short fragment of a single locus. thus, character sampling was not satisfactory for 3% of terminal taxa, and their placement must therefore be treated with caution. the second character matrix was restricted to the subset of those taxa for which morphological data was available. the reduced data set was tailored for a combined analysis to have a maximum taxonomic overlap of morphological and molecular character matrices with minimum missing data. this data set contained exemplars, representing all pterioidean genera and families. morphological data was scored for all terminals; of these terminals had a complete complement of molecular data, terminals had sequence data for three loci, and one terminal was represented by a single fragment of the 18s locus. the simultaneous analyses of molecular and morphological data were performed under maximum parsimony optimality criterion with the same cost for transformations of morphological characters as for nucleotide substitutions. in addition to the complete and reduced matrices of combined molecular data, five data sets corresponding to the combined nuclear data and the four individual loci were analyzed in the parsimony framework to explore topological differences, and the extent of resolution and support among the partitions. to evaluate the degree of incongruence between the data sets quantitatively, the incongruence-length difference tests were performed in paup* d. the ild tests were done simultaneously on the four individual locus matrices as well as pairwise on the complete data set. the tests were performed under parsimony using replicates, where each replicate consisted of random-addition sequences, with one tree retained per iteration, followed by tree bisection and reconnection branch swapping to completion. phylogenetic analyses tree searching and optimization phylogenetic analyses of dna sequence data were performed using direct optimization, a heuristic approximation to the optimal tree alignment method, under parsimony optimality criterion. a heuristic extension of generalized optimization, direct optimization includes insertion and deletion events in addition to substitutions as character transformations, allowing for the analysis of sequences of different lengths. additional analyses based on a static alignment and ml optimality criterion were performed to characterize the patterns of nucleotide substitutions and to explore potential sources of systematic error for inferences based on the assumptions of parsimony. these analyses were based on the implied alignment corresponding to the single optimal cladogram obtained from the parsimony analysis of the combined data performed by direct optimization under the equal-cost regime for the combined molecular character data set. parsimony analyses under direct optimization, as implemented in poy, were carried out in parallel on the american museum of natural history parallel computing cluster and a dual-processor personal computer. a stepwise tree searching strategy that combines independent multiple starting points and efficient algorithms designed for escaping local optima was used to insure sufficient sampling of tree space. initially, 250- wagner trees were generated by ras and submitted to branch swapping to completion by alternating rounds of subtree pruning and regrafting and tbr. the resultant optimal cladograms were subjected to iterations of tbr parsimony ratchet, upweighting 20% of nucleotide characters by a factor of and retaining one optimal cladogram per iteration. in the next step, the optimal cladograms were submitted to tree fusing, allowing for up to pairwise clade exchanges, which was followed by another round of tbr. finally, the resultant optimal cladograms were submitted to a final round of tbr under exact three-dimensional direct optimization, a more efficient, but substantially more computationally intensive method. the final equally most-parsimonious trees were rediagnosed to ensure the correctness of tree length calculation and were used to infer implied alignment. all tree searches were performed with unconstrained ingroup and outgroup designation. the ml analyses were performed using raxml. for independent replicates. the results were evaluated for topological congruence and convergence on the optimal likelihood value to ensure that the heuristic searches of tree space were sufficiently thorough. there appears to be no starting-point dependence to successive approximation that could potentially compromise the accuracy of approximation compared to full-optimization of ml searches. the best-fit model of substitution, given the optimal cladogram and the corresponding implied alignment, was selected by modeltest for the combined molecular data set using the akaike information criterion. the model favored by the akaike information criterion corresponded to the general time-reversible model that included estimation of the proportion of invariant sites and assumed a gamma distributed rate variation among sites. however, due to non-independence of p-invar and the shape parameter α of the gamma distribution, the ml analyses were performed under the gtr + Γ model. topological differences between the trees based on the partitioned and combined data do not present a problem for choosing among alternative models of sequence evolution, because tree topology does not strongly affect model estimation unless long internal branches are involved. the ensemble consistency index and the ensemble retention index for the optimal tree were calculated in poy with the dynamic-homology characters transformed into static nucleotide-level characters by static approximation. the robustness of phylogenetic relationships was evaluated using the bremer and jackknife support methods at nucleotide level. relative contribution of individual genes was assessed by partitioned bremer support. the batch files for constrained searches were generated by treerot. the constrained searches were performed by paup* d. for calculating jackknife support values, resampling iterations with 36% of sites removed during each pseudoreplicate were performed using poy. for the partitioned analyses, only jackknife support values were estimated. the resulting trees were visualized using poy and figtree.. the optimization and ancestral state reconstruction of morphological characters were evaluated using macclade. the significance of the length difference between alternative topologies was tested by the kishino-hasegawa test using resampling of estimated log-likelihoods bootstrap with replicates under ml criterion, and the wilcoxon signed-ranks test and the kh test under parsimony criterion as implemented in paup* d. the tests conducted under ml were performed using the gtr + Γ, the same model used for the ml tree searches. the likelihood-based tests are generally considered inapplicable for post-hoc tests of the optimal versus an alternative topology. in the present study, these paired-sites tests were used to assess the significance of alternative outgroup taxon configurations using constrained topologies defined a priori for the reduced molecular data set. indel treatment and sensitivity analysis because direct optimization integrates the implied alignment and the explicit topology, the choice of analytical parameters is critical for phylogenetic inference, as analyses performed under alternative cost matrices can produce widely different hypotheses of relationships. to evaluate the effect of varying alignment parameter values on cladogram topology, independent analyses of the complete data set were performed for a range of parameter combinations under parsimony. sensitivity analysis allows discerning between stable relationships and unstable relationships. twelve parameter combinations were used in the analysis: indel/substitution ratios of and 4; and tv/ts ratios of and ∞. when the tv/ts cost ratio was other than one, the insertion-deletion cost was set relative to the cost of transversions. transversions were never given a higher cost than indels. to account for a likely possibility that insertions and deletions spanning multiple contiguous nucleotide positions could have been accommodated by single events, additional sets of parameters taking into account affine indel costs were carried out. these analyses explored the indel-substitution ratios of and and tv/ts ratios of where the extension gap cost was set at a constant minimum. this corresponded to five ratios of gap opening/gap extension costs: and in addition, a cost regime that maximizes homology of both sequence fragments and individual nucleotide positions in the analysis of sequences of variable length was also examined. this weighting scheme, based on the view of parsimony as a two-taxon analysis, corresponds to the following combination of costs: for substitutions, for gap opening, and for gap extension. in total, parameter combinations were analyzed. the results of the sensitivity analysis are shown as binary plots to designate presence or absence of every clade recovered in the analysis of the complete data set under equal-costs regime. in addition, these results are summarized as a stability tree-a consensus of clades recovered under the entire parameter space explored. despite the fact that character and topological congruence, as measured by the ild and the topological ild respectively, has been proposed as a criterion for choosing among alternative weighting schemes, these measures are partition-dependent and their values are not comparable across different data matrices. because of the lack of such an objective criterion to decide among alternative weighting schemes, results from equal costs were preferred. results under additional parameter value combinations were used to assess clade stability.  dna sequence divergence a detailed description of dna sequence divergence patterns is provided in additional file 1: dna sequence divergence. as expected, the average levels and the range of divergence, based on the estimates of p distance, increased with relative taxonomic rank and, when compared across loci, were the smallest for the 18s and largest for the 16s sequences. with regard to the taxonomically problematic pinctada imbricata/fucata/radiata species complex, even though none of the markers could provide unambiguous cut-off levels for intra- vs. interspecific sequence divergence, in all cases the levels of divergence in pairwise comparisons between p. imbricata, fucata, and radiata exemplars were invariably lower compared to other interspecific comparisons, and were the same or marginally exceeding the limits of the intraspecific sequence divergence range inferred for other species within pinctada. nucleotide base composition overall base composition was homogeneous across all taxa: a χ test rejected the null hypothesis of base-composition stationarity neither for any individual locus, nor for the complete data set. the 16s sequences had a low proportion of cytosine and the prevalence of a + t, that has been previously reported for other molluscan taxa. across all loci, the sequences were to some extent gc-rich. the average and range of base composition for each data partition across taxa. the results of the χ value and the p value of the χ test are listed in the rightmost column. patterns of nucleotide substitution the entropy-based statistical test of the combined molecular data set showed no significant net effect of saturation. it revealed that the 18s and h sequences were unaffected by multiple hits, but suggested the presence of substantial substitution saturation in the 16s sequence. saturation of the 28s sequences can potentially become problematic when extremely asymmetrical topologies are considered, which are probably unrealistic given the results of the phylogenetic analysis. the overall value of the shape parameter α of the gamma distribution was, suggesting extensive heterogeneity in substitution rates across sites. there was a considerable disparity in the levels of α across loci: the highest level of across-site variation was found in 18s, followed by h, 28s, and 16s. it is noteworthy, that all the nuclear loci displayed substantial variation, whereas the mitochondrial locus showed relatively moderate variation in substitution rates among sites. the estimates of α appear to be accurate for such range of α. in h sequences it was likely due to codon bias: virtually no substitutions were inferred for the second position, less than 2% of sites experienced substitutions in the first position, whereas up to roughly 20% of sites in the third position underwent substitution. incongruence-length difference tests the ild tests performed simultaneously on all four molecular data partitions showed that the complete data set was homogeneous, suggesting that combining molecular partitions in a phylogenetic analysis was not likely to reduce phylogenetic accuracy. the pairwise ild tests, revealed, however, that only the rdna partitions were entirely congruent, whereas the h data set was congruent with 16s, but not with the nuclear rdna loci. despite the fact that the appropriate threshold of incongruence is not known, it has been argued that a significance threshold of might be too conservative for the ild test. given the overall homogeneity of the combined data and marginal incongruence of h data, the data-set incongruence was not likely to render an analysis of the combined data misleading. consequently, the phylogenetic inference based on the combined data was the preferred description of pterioidean relationships, but one must be cautious when interpreting the results with regard to the contribution of the h data. the ild tests together with phylogenentic analyses of individual data partitions and the assessment of their relative contribution to support on a combined-data cladogram were used to assess the distribution, nature, and extent of agreement among data sets. combined analyses of molecular data a parsimony analysis of the complete molecular data set performed under equal-costs regime produced a single, well-resolved optimal cladogram. the superfamily pterioidea was recovered as a monophyletic, well-supported clade, which was stable under most combinations of alignment parameters. the few cases of non-monophyly of the superfamily were due entirely to the placement of pulvinites exempla among the outgroup taxa as a sister group to the ostreoidea. this arrangement was recovered only under parameter combinations characterized by the highest indel/substitution ratio, suggesting that the spurious placement was due to missing data, an expected artifact given the fact that p. exempla was the only taxon for which only a small fraction of sequence data was available. excluding p. exempla from the analysis did not affect the relationships of the remaining taxa, suggesting that the effect of the missing data for this taxon on the relationships among remaining ingroup taxa was negligible. the suprageneric relationships were completely resolved and all polytomies were restricted to a single derived subclade within the isognomon clade. all traditionally recognized families-with the obvious exception of the monotypic pulvinitidae-were recovered polyphyletic, whereas all genera were definitively monophyletic. the monophyly of the genera pulvinites and crenatula could not be tested in the present study because these genera were represented by single exemplars. thus, the polyphyly of the families was a result of unconventional relative placement of traditionally recognized genera. the genus-level clades were distributed in a fully pectinate pattern of relationships with the basal-most pulvinites and the others branching off in the following sequence: malleus, pteria, isognomon, pinctada, vulsella, and electroma. the nodes along the backbone had relatively high support and were stable under the majority of alignment parameters with the exception of the vulsella/electroma/crenatula group. despite the high support values obtained under the equal-costs regime, the monophyly of this clade was recovered in only 41% of alignment parameter combinations. the sister-group relationship of crenatula and electroma was robust and insensitive to parameter choice. the monophyly of all the genera was very strongly supported and entirely stable across alignment parameter space. the monophyly of higher taxa among the outgroups was strongly supported in all analyses and was insensitive to the choice of alignment parameters. the ostreoidea was identified as an immediate outgroup to the pterioidea with high support under equal-cost regime, but the pterioidea/ostreoidea sister-group relationship was highly sensitive to alignment ambiguity and was obtained for only four combinations of alignment parameters. effectively, all possible arrangements of the considered superfamilies were recovered. the cladogram generated under the parameter set that maximizes homology of both sequence fragments and base correspondences, had a similar overall structure to the cladogram generated under uniform weighting with several important differences: the placement of isognomon and pteria clades were exchanged; the vulsella/electroma/crenatula clade was not recovered; and the pinnoidea was an immediate pterioidean outgroup. relative contributions of different data partitions, calculated as partitioned bremer support values for each node, are shown in figure each locus provides support for some deep and apical nodes, but not all of them. there also appears to be a slight taxonomic bias in the distribution of the pbrs values: for example, relatively high pbrs values for most nodes within pteria were provided by 28s, whereas comparable in magnitude pbrs values for the majority of the nodes within pinctada were provided by 16s. taken together, 18s and 28s bring close to 75% of the signal and contribute the highest on average support values; h provides the least support and the lowest on average support values. when the support of each gene was normalized by the number of informative characters, 18s appeared as most informative, providing more than twice the support as 28s, whereas 16s and h were the least informative markers in this data set. a ml analysis of the complete molecular data set produced in a single optimal tree. the values of suboptimal trees from % of search replicates were within % of the optimal likelihood value; of these had identical to the optimal likelihood values to the third decimal point and differed only in the placement of several conspecific taxa. the estimated shape parameter α of gamma distribution was; the substitution relative rate parameters for the model were,,,,,. comparison of the optimal parsimony and ml topologies showed non-substantial differences with regard to placement of several species-level taxa within pteria and malleus, and alternate arrangement of conspecific individuals of pinctada imbricata. a parsimony analysis of the reduced molecular data set performed under equal-costs regime resulted in a single well-supported cladogram. in spite of a great difference in taxon sampling between the complete and reduced data sets, the higher-level relationships were nearly identical with only two notable differences: in the analysis of the reduced data set pulvinites exempla emerged as a sister taxon to the outgroup crassostrea virginica and a pinnid, atrina rigida, was recovered as an immediate outgroup to the pterioidea. as suggested by the analyses of the complete data set, missing data and sensitivity to alignment ambiguity coupled with the limited taxon sampling of the reduced data set probably account respectively for the instability in the position of p. exempla and the placement of the immediate outgroup. the phylogenetic analyses of individual molecular data partitions resulted in cladograms that differ in topology, degree of resolution, and relative support measures. however, all the analyses invariably recovered the pterioidea and major genus-level clades monophyletic. combined analyses of morphological and molecular data a parsimony analysis of the combined morphological and the reduced molecular data sets performed under equal-costs regime produced a single, highly supported tree. when compared to the tree generated for the reduced data set analyzed under the same conditions, the addition of morphological data reconstituted monophyletic the pterioidea by placing pulvinites exempla at the base of the pterioidean subtree; the placement of the remaining taxa was identical. similarly, the arrangements of pterioidean genera in the analyses of the combined reduced and complete molecular data sets were identical with the only significant topological difference being the identity of the immediate pterioidean sister-group: the pinnoidea in the former and the ostreoidea in the latter. the paired-sites tests, applied to three topologies constrained to the identical arrangement of pterioidean taxa, but differing solely in the placement of the outgroup taxa, showed that the compared trees were not significantly different. only the monophyly of the ostreoidea-pinnoidea clade could be rejected in one set of tests under parsimony. the significance of the length difference between alternative hypotheses of relationships were evaluated using the kishino-hasegawa test under parsimony and maximum likelihood, and the wilcoxon signed-ranks test. * = significant difference, p <. other abbreviations: ost, ostreoidea; pin, pinnoidea; pte, pterioidea. a parsimony analysis of the combined morphological and the complete molecular data sets produced a well-resolved and strongly supported single optimal cladogram. the topology was essentially the same as of the tree obtained under the same conditions for the complete molecular-only data set: the higher-level relationships along the backbone were identical. the few differences included an alternative placement of several clades within pteria and the rearrangement of several, probably conspecific, representatives of isognomon. the most significant discrepancy between the analyses was the placement of electroma papilionacea. in the combined analysis it was recovered at the base of the ) clade, whereas in the molecular-only study it was nested within it: )). discussion reliability of phylogenetic inference the present study provides the first explicit hypothesis for higher-level phylogeny and species-level relationships for the bivalve superfamily pterioidea based on dna sequence data. nearly complete species-level taxon and character coverage suggest that the random error associated with sampling was not likely to compromise the fidelity of the results inasmuch as supraspecific relationships are concerned. uncertainty in species-level relationships does, however, remain, due to incomplete taxon sampling, differences in character sampling among congeneric individuals, and potentially insufficient phylogenetic signal of chosen markers for resolving recent divergences. sequence alignment was straightforward and unambiguous for the protein-coding gene. for the ribosomal sequences, the effect of the uncertainty in nucleotide positional homology on stability of inferred relationships was evaluated using the sensitivity analysis, which allowed for recognizing specific parts of the phylogeny affected by the uncertainty in nucleotide homology and evaluating the extent of its effect. nucleotide base composition was nearly perfectly stationary across taxa overall as well as for individual loci, without evidence of taxon attraction due to similar nucleotide-composition bias. substitution saturation in dna sequences can potentially underestimate the actual amount of evolutionary change and reduce the phylogenetic signal due to excessive multiple hits in the same nucleotide positions resulting in extensive homoplasy. despite the fact that extensive saturation was detected for transitions in the 16s sequence, its effect on the phylogenetic inference was negligible, as shown by comparison of the results from individual data partitions, incorporating the ts/tv into phylogenetic analysis, character congruence among data partitions, and consistency in the distribution of brs values. furthermore, excluding 16s data from the analysis resulted in a decreased resolution and support throughout the tree. these results are consistent with theoretical and empirical findings suggesting that homoplasy does not necessarily confound phylogenetic inference and that a given site, which is homoplastic in a sub-part of the most-parsimonious tree, can bring phylogenetic information for another sub-part of the same tree. heterogeneity in substitution rates among sites within a nucleotide sequence can have considerable impact on sequence divergence. substantial heterogeneity in substitution rates among lineages resulting in long-branch attraction is another source of misleading inference in molecular data. the agreement between the results of the maximum likelihood analysis under the model allowing for variable base frequencies and parsimony suggested that the predominance of transitions and among-site heterogeneity for most loci, as well as among-lineage heterogeneity in nucleotide substitution rates were not likely to be a source of misleading inference. the absence of long-branch attraction despite considerable differences in branch lengths was consistent with prior observation that extensively sampled phylogenies are not sensitive to this artifact. the doubly uniparental inheritance, which has been reported for a number of bivalve species, can potentially affect the accuracy of phylogenetic reconstruction due to varying degrees of divergence and homogenization between conspecific genomes. within the subclass autolamellibranchiata, the instances of dui were reported only from the most basal extant group, mytiloida. even though the lack of reported instances of dui in a particular group does not necessarily imply that it is not present, the absence of strict dioeciousness in the pterioidea suggests that dui was not a likely to be a problem. all pterioidean species for which reproductive strategies have been elucidated, are rhythmical protandrous hermaphrodites, with occasional instances of protogyny and simultaneous hermaphroditism, and mtdna appears to be maternally inherited. simultaneous analysis of a single data matrix combining independent character sets potentially increases the descriptive efficiency and explanatory power, provides the basis for evaluating the relative support and taking into account the hidden support provided by individual partitions, and is less sensitive to sampling bias. however, the incongruence among data partitions may indicate that they have evolved under different models of evolution and render a combined-data analysis misleading. despite topological differences obtained in the analyses of individual data partitions, the results of ild tests showed that the data sets were largely congruent, implying that a single phylogeny underlies the relationships inferred from individual data partitions. differences in substitution rates across sites, the average pairwise sequence divergence values, the distribution and levels of support, and the extent of resolution obtained in the analyses of individual partitions indicated that the marker loci collectively encompassed sufficiently broad spectrum of dna sequence divergence to produce a nearly completely resolved and well-supported tree. these results also suggested that topological differences among data partitions were probably due to characters that evolved rapidly to provide a reliable signal for sub-terminal nodes of a tree, but produced noise for the resolution of deeper nodes. the overall increased resolution and support of the combined-data trees when compared to the results of partitioned analyses strongly speak for the simultaneous analysis approach. in summary, the present molecular sequence data and the choice of methodology adequately account for sampling biases and potentially confounding effects of various molecular evolutionary processes, providing a strong foundation for accurate, reliable, and robust reconstruction of phylogenetic relationships of the recent pterioidea. phylogeny of the pterioidea the present study corroborates the monophyly of the superfamily pterioidea and definitively rejects that of all nominal pterioidean families. whereas the non-monophyly of some traditional pterioidean groups was anticipated by prior studies, the arrangement of the genera was novel and unexpected, as well as resolved and strongly supported. in terms of taxonomic congruence with morphological studies, the analysis of molecular data offered a considerably different phylogenetic hypothesis. generally assumed close relatedness of commercially important genera pteria and pinctada was rejected by the molecular evidence, which was not entirely surprising given the absence of unambiguously optimized non-homoplastic characters supporting the sister-group relationship of these taxa. the presumed close relatedness of the genera malleus and isognomon, implied by superficial resemblance in shell morphology and similar life habits, was also contradicted by the molecular data. the only morphological character supporting a close relatedness of these two genera-the absence of food grooves on inner demibranchs of the gills-is symplesiomorphic in the context of the dna based topology, as this condition has been also reported from other autolamellibranchiate taxa. contrary to previously proposed hypothesis, this study rejected a close affinity of isognomon and pulvinites, placing the latter at the base of the pterioidean tree. this finding was corroborated by the fact that in the morphology-based analysis the sister-group relationship of the two genera was weakly supported and relied on a unique combination of shared homoplastic characters. the present analysis provided the first molecular evidence for a peculiar clade, recently recognised on morphological grounds, that consists of relatively understudied, low-diversity genera electroma, vulsella, and crenatula, that were traditionally placed in pteriidae, malleidae, and isognomonidae respectively. the arrangement among the genera, however, varied between the two studies. in the morphology-only analysis, where each genus was represented by a single exemplar, electroma was basal to vulsella and crenatula, the sister-group relationship between which was implied by the reduction of the pedo-byssal apparatus. instead, in the molecular analysis, vulsella appeared basal to the electroma clade that nested the sole representative of crenatula, rendering electroma paraphyletic. an unanticipated discovery of polyphyletic electroma requires detailed morphological and taxonomic reappraisal of this poorly known genus. another species of electroma, typically placed in its own subgenus pterelectroma, e. physoides, invariably grouped with the pteria clade, putting into question the taxonomic status of this subgenus. due to high sensitivity to alignment parameters and taxon sampling, uncertainty remains regarding the identity of the sister group of the pterioidea. the paired-sites tests, applied to the alternative topologies differing only in the placement of the outgroup taxa, showed that the trees were not significantly different. consequently, the sister-group of the pterioidea could not be unequivocally inferred in light of the present results. based on the results of the combined-data analyses under uniform weighting and predominant support, the ostreoidea is tentatively proposed here as the most likely pterioidean sister group. implications for homology of morphological characters it has long been acknowledged that extensive homoplasy is a recurrent theme in bivalve phenotypic evolution that can potentially mislead phylogenetic inference based on morphological data alone. optimizing morphological character state transformations onto topology derived from the analysis of molecular data provides means to discern putative homologies and instances of homoplasy. the molecular analyses presented here have confirmed a number of suspected convergences identified on the basis of morphology, but also uncovered a number of hitherto unanticipated ones. the hypotheses of homology of several key characters were evaluated in light of a robust phylogenetic hypothesis based on molecular evidence. a duplicated outer fold of the mantle edge has been considered a derived trait that evolved twice independently in the pterioidea based on the analysis of morphological data and in agreement with presumably "primitive" three-fold arrangement. in contrast, molecular data suggested that the four-fold condition was primitive in the pterioidea, whereas the three-fold condition was a result of a single reversal to the three-fold condition exhibited by the outgroup taxa. this finding was not entirely unexpected, given the fact that analogous transitions in mantle evolution had occurred in other bivalve lineages. in the analysis of morphological data, homorhabdic ctenidia evolved once from heterorhabdic ctenidia. molecular results, on the other hand, implied that the gill grades had evolved iteratively and the ancestral condition could not be unequivocally resolved. the presence of food grooves in the outer demibranchs showed the same pattern. the loss of postlarval byssus and posterior pedo-byssal retractor muscles has been previously considered as evidence for sister-group relationship of vulsella and crenatula. both features are likely to be adaptations for living inside sponges, the habit shared by members of the two genera, where physical stabilization is achieved by the attachment to the host sponge. the relationships inferred from the molecular data invariably grouped crenatula with electroma, a relationship supported by a single unambiguously optimized nonhomoplastic apomorphy: a t-shaped typhlosole, a longitudinal fold of the intestine. molecular results were ambiguous with regard to the immediate relatedness of the vulsella and the crenatula/electroma clades due to alignment uncertainty, and suggested a possibility that the evolution of adaptation for living inside sponges could have been two separate events. this supposition is corroborated by prior ecological observations that the host specificity, and the nature of mutualistic association of vulsella and crenatula species with their demosponge hosts are different: the former is a facultative mutualist with keratosid sponges, whereas the latter is an obligatory mutualist with monaxonid sponges. species of vulsella possess a relatively thick and elongated shell with an extremely abbreviated hinge line and a narrow prismatic margin that provide a rigid "endo-skeleton" for the sponge, hence facilitating the host's growth. species of crenatula, on the other hand, have a very thin shell with a very wide prismatic margin, which appears to be influenced by the sponge growth, and produce highly irregular shell shapes without consistent orientation. the particular form of the multivincular ligament of crenatula allows for curvature of the extensive hinge line, which facilitates attachment along the hinge line to irregular surfaces. the loss of postlarval byssus in distantly related malleus albus also supports the likelihood of convergence. the symmetry and relative placement of abdominal sense organs, paired mechanoreceptors typically found on the surface of the adductor muscle, have been previously shown to be phylogenetically informative within autolamellibranchiata and, specifically, within the pterioidea. in the morphology-based analysis, the clustering of the asos around the anus was inferred as a primitive condition in the pterioidea, whereas the dorsal displacement of the right aso had occurred twice. according to the results of the molecular character analysis, both aspects of the aso position had evolved twice. a set of characters that has historically been used for higher-level systematics of pterioidean and other autolamellibranchiate taxa is associated with the types of ligament and their support structures. according to the molecular-based topology, the multivincular ligament characterizes three distantly related genera: pulvinites, isognomon, and crenatula. morphological data alone also identified an independent origin of multivincular ligament in crenatula, but suggested a sister-group relationship of pulvinites and isognomon on the basis of shared ligament morphology and two other homoplastic characters. it is noteworthy that despite apparent similarity between multivincular systems in pulvinites and isognomon, the pattern of ligostracum deposition and the ontogenetic sequence of origination of ligamental layers along the hinge line differ between the genera. even though not entirely congruent, the results of the molecular and morphological studies reject the long-standing assumption of close relatedness based on a single origin of multivincular ligament in the pterioidea. the present molecular evidence confirmed the taxonomic significance of the shape of the anal funnel, an ear-shaped membranous structure situated at the posterior extremity of the rectum that presumably facilitates expulsion of fecal pellets. the absence of the anal funnel in pulvinites is consistent with its basal placement in the pterioidea as the same condition is displayed by the outgroup taxa. in contrast, in the analysis based on morphological data, the lack of the funnel was inferred as a loss. the molecular and morphological analyses agree that the anal funnel is subtriangular in outline in the more basal pterioidean clades, whereas it is characteristically lanceolate and rounded in two most derived clades, pinctada and vulsella/crenatula/electroma, respectively. given the uncertainty in the identity of pterioidean sister group, it is instructive to explore the effect of alternative combinations of outgroup taxa on the distribution of diagnostic morphological characters using the best available hypothesis of the ingroup relationships. to this end, morphological characters were optimized on three trees constrained to the identical arrangement of the ingroup taxa but differing in the placement of the outgroups. the pinnoidea relationship was defined by five unambiguously optimized, non-homoplastic apomorphic characters that remained unchanged above this node: the descending intestine produced towards the posteroventral side of the posterior adductor muscle; the loss of the anterior adductor muscle; the absence of the pseudonymphae; the presence of a resilifer throughout ontogeny; and the placement of the posterior adductor muscle entirely within the boundary of the inner shell layer. previously, steiner and hammer had mapped a set of morphological characters on to their maximum likelihood tree substantiating their finding of ostreoidean/pterioidean sister-group relationship by two characters: the monomyary and the presence of the anal funnel. the present study confirmed the monomyaryan condition as a synapomorphy of the clade, but not the occurrence of the anal funnel in different superfamilies, interpreted here as a result of independent origins, as was previously suggested on the basis of morphological data. the ostreoidea relationship was defined by three apomorphic characters: the absence of the marginal mantle lobe fusion; the absence of the eulaterofrontal cirri of gill filaments; and the ascending intestine not produced anteriorly passed the stomach. the pterioidea relationship was defined by two apomorphic characters: the absence of the ciliated discs on gill filaments and the clustering of ducts of the digestive diverticula. given comparable levels of homoplasy and the lack of a priori criteria to differentially weigh the characters supporting each of the possibilities, the arrangement supported by the largest number of apomorphies and yielding the shortest tree length, that is ), is the preferred hypothesis of relationship. species-level taxonomy of pinctada despite the fact that the focus of the present study was to establish higher-level relationships within the pterioidea, it provided significant insights into species-level taxonomy of the genus pinctada, a commercially important and the most exhaustively sampled group. the present study suggests that pinctada imbricata, p. fucata, and p. radiata are distinct genealogical units. the taxonomic status of the p. imbricata/fucata/radiata species complex previously remained unsettled due to extensive variation in shell characters within and among populations, extremely wide geographical distribution, and transport and hybridization by humans. traditionally, three distinct species were recognized corresponding to three biogeographic realms: p. imbricata in the western atlantic region, p. radiata in the eastern indian ocean and the red sea regions, and p. fucata in the indo-pacific region. the japanese populations have frequently been regarded as a distinct species, p. martensii, or a subspecies, p. fucata martensii. because of its major role in perliculture, the species complex had become a major research focus in recent years. a large number of comparative genetic analyses using allozyme and dna sequence data indicated significant levels of gene flow over large geographic distances among populations of indian ocean and indo-pacific regions, resulting in low levels of genetic differentiation. the pattern of genetic differentiation among the chinese, japanese, and australian populations using amplified fragment length polymorphism markers was consistent with isolation by distance. shell morphometric studies, mating experiments, and comparative karyotype analyses failed to establish unequivocal species delineations. to date, no diagnostic discrete morphological characters have been identified for populations from the three geographical regions. the present results showed shallow levels of divergence among the members of the p. imbricata/fucata/radiata species complex, comparable to the levels of divergence of conspecific exemplars of other pinctada species. however, the three groups were reciprocally monophyletic showing an apparent clinal geographical structure. therefore, it might be meaningful to treat these individual populations as evolutionary significant units. taking into account the comparatively low divergence values and the lack of diagnostic morphological characters distinguishing these groups, from the taxonomic standpoint it is advocated here to provisionally recognize them at a subspecies level under the senior synonym, p. imbricata. the final word on the taxonomic status of these taxa will depend on the future population-level analysis with dense sampling throughout the entire distribution range and using genetic markers more sensitive for resolving recent divergencies. this study also suggested that pinctada albina and p. nigra, generally distinguished by shell color, represent two colormorphs of the same species. this possibility has been previously suggested based on extremely low levels of genetic divergence between two exemplars of p. albina and one of p. nigra. this finding was consistent with other molecular studies that invariably resolved these species as sister taxa. the present analysis contained three exemplars of each colormorph that were always recovered as a monophyletic group and never formed separate clades diagnosed by color. in addition, based on the fact that p. mazatlanica invariably nested among representatives of p. margaritifera or members of both species formed an unresolved clade, the present study provides preliminary support for conspecificity of these two species. these findings need to be investigated further and cannot be definitively tested given the present sample size.  the most significant conclusion from the phylogenetic analyses presented here is that the state of pterioidean systematics is no longer in chaos. the present results establish a solid-however unconventional- framework for the higher-level taxonomy of the group, focusing future revisionary systematic effort and providing a better grasp of standing alpha-diversity. the phylogeny is crucial for understanding processes and patterns of diversification and morphological evolution of the pterioidea through time. to this end, the major future challenge lies in integrating the knowledge on the recent pterioideans with information on the extensive, well-documented fossil record of the group. as throughout their history, pterioideans remain a globally distributed, relatively common group of bivalves found in most warm, shallow marine environments around the globe, and are likely to display patterns of diversity typical of other marine benthic macrofauna. therefore, the insights gained from the phylogenetically-informed studies of pterioidean bivalves might be broadly applicable across marine invertebrates. supplementary material additional file 1 dna sequence divergence. word doc file containing dna sequence divergence. click here for file additional file 2 partitioned analyses. word doc file displaying partitioned analyses. click here for file acknowledgements i would like to thank friends and colleagues around the world who donated or loaned specimens: r. bieler, a.j. baldinger, p. bouchet, j. gerber, g. giribet, l. groves, m.g. harasewych, d. herbert, k. holmes, japan pearl exporters, m. knapp, k.v. ramana murty, e. kools, p. maestrati, b. marshall, p.g. oliver, c.j.m. lodeiros seijo, r. shearer, l.r. simone, a. tommasi, k. umadevi, Á. valdés, n. voss, and r.c. willan. g. steiner and t. masaoka provided information on specimens, sequences of which were deposited in genbank. i am indebted to e. arizmendi castillo for providing assistance during my visit of guaymas, mexico in r. rorandelli, f. bartolini, m. fusi, m. vannini, and m. gomei have provided invaluable help during fieldwork in the tyrranean sea in summer i would like to thank all participants of the second international marine bivalve workshop and c. printrakoon in particular for partnership during the workshop and continuing collaboration. r. desalle generously made available molecular laboratory resources at the sackler institute for comparative genomics. w.c. wheeler granted access and logistical support for using the parallel computing cluster. m. leslie provided invaluable help for the laboratory component of the project. m. coppolino and s. watson assisted with the american museum of natural history mollusk collection. y. iwasaki had assisted with translations from japanese. f. almeida, a. varón, m. weksler, and x. xia made technical suggestions. p.m. mikkelsen, e. strong, j. hammock, and three anonymous reviewers provided critical suggestions for improving the manuscript. partial funding toward completion of this work was provided by nsf-peet marine bivalve project deb-, nsf dissertation improvement grant deb- amnh annette kade graduate fellowship, the conchologists of america grant, and nsf-atol bivalve tree of life project deb-0732854/0732903/.
this issue is examined here in light of molecular data sets composed of dna sequences for nuclear and mitochondrial loci, and a published character data set of anatomical and shell morphological characters. recent morphological treatment of selected pterioideans and molecular phylogenetic analyses of higher-level relationships in bivalvia have challenged the traditional view that pterioidean families are monophyletic. the present study is the first comprehensive species-level analysis of the pterioidea to produce a well-resolved, robust phylogenetic hypothesis for nearly all extant taxa. four primary conclusions were reached: the results strongly supported the monophyly of the pterioidea; none of the previously defined families were monophyletic; the arrangement of the genera was novel and unanticipated, however strongly supported and robust to changes in alignment parameters; and optimizing key morphological characters onto topologies derived from the analysis of molecular data revealed many instances of homoplasy and uncovered synapomorphies for major nodes. the congruence and contribution of different partitions were quantified, and the sensitivity of clade stability to alignment parameters was explored. the data were analyzed for potential biases due to taxon and character sampling, and idiosyncracies of different molecular evolutionary processes. additionally, a complete species-level sampling of the genus pinctada provided further insights into the on-going controversy regarding the taxonomic identity of major pearl culture species.
 the superfamily pterioidea is a morphologically and ecologically diverse lineage of epifaunal marine bivalves distributed throughout the tropical and subtropical continental shelf regions. this group includes commercially important pearl culture species and model organisms used for medical studies of biomineralization. recent morphological treatment of selected pterioideans and molecular phylogenetic analyses of higher-level relationships in bivalvia have challenged the traditional view that pterioidean families are monophyletic. this issue is examined here in light of molecular data sets composed of dna sequences for nuclear and mitochondrial loci, and a published character data set of anatomical and shell morphological characters.  the present study is the first comprehensive species-level analysis of the pterioidea to produce a well-resolved, robust phylogenetic hypothesis for nearly all extant taxa. the data were analyzed for potential biases due to taxon and character sampling, and idiosyncracies of different molecular evolutionary processes. the congruence and contribution of different partitions were quantified, and the sensitivity of clade stability to alignment parameters was explored.  four primary conclusions were reached: the results strongly supported the monophyly of the pterioidea; none of the previously defined families were monophyletic; the arrangement of the genera was novel and unanticipated, however strongly supported and robust to changes in alignment parameters; and optimizing key morphological characters onto topologies derived from the analysis of molecular data revealed many instances of homoplasy and uncovered synapomorphies for major nodes. additionally, a complete species-level sampling of the genus pinctada provided further insights into the on-going controversy regarding the taxonomic identity of major pearl culture species. 