Skip to main content

A multi-gene phylogeny of Cephalopoda supports convergent morphological evolution in association with multiple habitat shifts in the marine environment

Abstract

Background

The marine environment is comprised of numerous divergent organisms living under similar selective pressures, often resulting in the evolution of convergent structures such as the fusiform body shape of pelagic squids, fishes, and some marine mammals. However, little is known about the frequency of, and circumstances leading to, convergent evolution in the open ocean. Here, we present a comparative study of the molluscan class Cephalopoda, a marine group known to occupy habitats from the intertidal to the deep sea. Several lineages bear features that may coincide with a benthic or pelagic existence, making this a valuable group for testing hypotheses of correlated evolution. To test for convergence and correlation, we generate the most taxonomically comprehensive multi-gene phylogeny of cephalopods to date. We then create a character matrix of habitat type and morphological characters, which we use to infer ancestral character states and test for correlation between habitat and morphology.

Results

Our study utilizes a taxonomically well-sampled phylogeny to show convergent evolution in all six morphological characters we analyzed. Three of these characters also correlate with habitat. The presence of an autogenic photophore (those relying upon autonomous enzymatic light reactions) is correlated with a pelagic habitat, while the cornea and accessory nidamental gland correlate with a benthic lifestyle. Here, we present the first statistical tests for correlation between convergent traits and habitat in cephalopods to better understand the evolutionary history of characters that are adaptive in benthic or pelagic environments, respectively.

Discussion

Our study supports the hypothesis that habitat has influenced convergent evolution in the marine environment: benthic organisms tend to exhibit similar characteristics that confer protection from invasion by other benthic taxa, while pelagic organisms possess features that facilitate crypsis and communication in an environment lacking physical refuges. Features that have originated multiple times in distantly related lineages are likely adaptive for the organisms inhabiting a particular environment: studying the frequency and evolutionary history of such convergent characters can increase understanding of the underlying forces driving ecological and evolutionary transitions in the marine environment.

Background

Large-scale shifts in habitat during evolution require populations to respond to new selective pressures, often resulting in a cascade of novel morphologies. For example, the initial land invasion by marine tetrapods required significant changes in physiology, morphology, and life history, which led to the origin of amniotic eggs, modified limbs, scales to reduce water loss, and modified lungs (e.g. [1]). Distantly related lineages that transition to similar habitats are under the same kinds of selective pressures, often resulting in convergent structures. Canonical examples include spines in plants inhabiting arid environments and flight morphologies in insects, birds, and bats. Convergence remains largely understudied in many marine invertebrates, where continued habitat shifts and diversification have led to a number of diverse lineages that possess multiple potentially convergent characters.

Potentially adaptive morphological features shared between distantly related organisms inhabiting similar environments provide an opportunity to test for correlation between morphology and habitat to provide new hypotheses about the factors influencing evolution and diversification. For example, stickleback fish have long been studied for their numerous instances of convergent characters such as body shape, pelvis morphology, skeletal armor, and mating preferences, associated with transitions from marine to freshwater benthic and limnetic habitats (e.g. [25]). In sticklebacks, recent work by Jones et al. [6] identified a core set of genome-wide loci that are associated with the transition from marine to freshwater habitats, shedding light on the underlying genes that are responding to and driving adaptive evolution. In other marine taxa, numerous cases of convergent evolution have been found, including: very similar antifreeze glycoproteins that occur in phylogenetically distinct Arctic and Antarctic fishes [7], the presence of marine bioluminescence in almost all marine invertebrate phyla (e.g. [810]), and the thunniform body shape seen in whales, ichthyosaurs, and several species of large pelagic fishes [11]. However, in cephalopods, studies that examine the potential for correlation between habitat and morphology have been largely limited to physiology. Childress [12] argued that many of the observed relationships between metabolism and depth in oceanic squids were the result of adaptation rather than phylogenetic history, because similar metabolic rates consistently coincide with depth, even across distantly related taxa. Seibel and Carlini [13] tested this hypothesis in a phylogenetic context using a molecular-based phylogeny, and found that correlation between changes in metabolism and depth existed in the pelagic squids.

The marine environment can be broadly divided into two macrohabitats, each promoting a suite of adaptations in the inhabitants. The largely two-dimensional demersal habitats include both the benthic and benthopelagic habitat, where organisms live on or just above the bottom, respectively. Demersal habitats are dominated by animals with defensive techniques that do not require extended periods of rapid movement, such as releasing ink and then jetting away, burrowing in sediment or emitting toxic substances. Many demersal fishes and cephalopods commonly utilize cryptic behaviors that involve changing body coloration, texture, and/or shape to avoid predation (e.g. [14, 15]). In contrast, the three-dimensional pelagic environment is dominated by animals that are capable of rapid movement, that often have body plans capable of responding to stimuli from several directions, and that use cryptic techniques such as transparency, reflection, and ventral bioluminescence (e.g. [14, 16]) to allow for improved predation, communication, and predator avoidance. Transitions between demersal and pelagic habitats during evolution should promote shifts in suites of morphological features, a hypothesis that can be examined with phylogenetic comparative analyses on lineages living in both habitats.

Modern cephalopods occupy nearly every habitat of the marine environment from the rocky intertidal to the deep sea, making them an excellent group in which to test for convergent evolution. The two most speciose lineages, Decapodiformes (pelagic squids, bobtail squids, and cuttlefishes) and Octopodiformes (octopuses and Vampyroteuthis infernalis) are each comprised of demersal and pelagic taxa. Although the two cephalopod lineages are morphologically very distinct, some features seem to appear consistently in species inhabiting similar environments. Demersal taxa generally tend to have a more bulbous body, corneas covering the eye, and exhibit some of the most advanced crypsis seen in the marine environment (e.g. [14, 1719]). Pelagic cephalopods often possess photophores for communication and crypsis, fins to help with propulsion, the ability to become almost entirely transparent, and a flexible internal shell for stability. However, there are exceptions to these trends. For example, the squid family Loliginidae, the most common group in the lineage Myopsida (Decapodiformes) is an important fisheries group that lives in coastal waters and lays its eggs in the benthic environment, but feeds in the water column (e.g. [20]). Loliginids have a streamlined, muscular body morphology and internal shell similar to many pelagic squids, while bearing corneas reminiscent of benthic squids, bacteriogenic photophores (in some species), and accessory reproductive structures found largely in benthic species.

To study the phylogenetic and environmental forces driving character evolution requires a well-sampled phylogeny. In cephalopods, molecular data sets have provided many new phylogenetic hypotheses, but have either had limited taxon sampling (e.g. [13, 21]) or recovered limited resolution for some major extant lineages [2225]. Here, we assemble previously published data from up to 10 genes for over 400 cephalopod taxa (representing 42 of 47 known families) to create the most comprehensive molecular analysis of cephalopod phylogeny to date, and subsequently use the phylogeny to test the hypothesis that habitat shifts drive convergent evolution in cephalopods. We first use the phylogeny to reconstruct ancestral character states for several characters whose evolutionary history remain elusive: the accessory nidamental gland, a specialized structure that provides antimicrobial protection to eggs laid on the ocean floor; the branchial canal, a structure often thought to be a synapomorphy uniting bobtail squids, cuttlefishes, and loliginids; the cornea, which covers the eye in many demersal taxa and whose status as a synapomorphy remains unclear; the presence of the right oviduct, another important potential synapomorphy for major clades; and photophores, which occur in many lineages of Cephalopoda (e.g. [22, 26, 27]). We then are able to test whether any of these characters are correlated with either demersal or pelagic lifestyles, thus providing new information regarding convergent evolution in the marine environment.

Results

Based on extensive new phylogenetic analyses including most available cephalopod molecular data, we find that transitions between demersal and pelagic habitats were common during the evolutionary history of cephalopods. We report strong evidence for convergent evolution of multiple phenotypic traits and correlated evolution between habitat and some of these traits, including accessory nidamental glands (ANG), bioluminescent photophores, and corneas. Below, we detail our results relating to new phylogenetic analyses, convergent trait evolution, habitat shifts, and correlations between habitat and morphology. These results are indicative of convergent adaptations to demersal or pelagic marine habitats in cephalopods.

A well-sampled phylogeny

Our phylogenetic analyses include the highest degree of taxonomic diversity of any published study on cephalopods to date (42 of 47 coleoid families). This phylogeny allows explicit testing of hypotheses of convergent adaptations of cephalopods to different marine environments. Results from phylogenetic analysis of the reduced 188-taxon dataset (Figure1) confirms monophyly for almost all major orders and for all families studied except Octopodidae, which is a highly morphologically divergent group with numerous species complexes [28]. The sepioid orders at the base of the decapodiform clade did not form a monophyletic group, and the position of Loliginidae was ambiguous. Such ambiguity may be due to factors such as substitution saturation and/or rate heterogeneity: we found saturation in several mitochondrial genes, more pronounced in decapodiforms than octopods (Additional file 1 Figures S5-S14), which could hinder attempts to resolve ancient, rapid radiations [29]. We accounted for this phylogenetic uncertainty in tests of character evolution by analyzing all hypotheses on multiple phylogenies estimated from bootstrapped, pseudoreplicated data sets. Despite phylogenetic uncertainty, we uncovered clear signals of convergent and correlated character evolution.

Figure 1
figure 1

Maximum-likelihood topology of the dataset containing taxa represented by more than three loci. Support values of nodes were generated via 1000 bootstrap replicates and nodes are labeled as follows: black circle (100%), grey circle (90-99%), grey square (80-89%), white square (70-79%), white triangle (51-69%). Nodes lacking symbols represent bipartitions recovered in less than 50% of the bootstrapped topologies. The presence of sequence data for a particular locus is indicated by a black square following the terminals.

Multiple transitions between demersal and pelagic habitats

Our phylogenetic analyses support multiple habitat transitions during the history of cephalopods, which allows us to test whether convergent phenotypes are associated with such shifts. Habitat reconstruction across the 1000 bootstraps of our 188-taxon tree (Figure2, condensed to family, see also Additional file 1 Figure S2), indicates that transitions to the pelagic realm are relatively common, while reciprocal transitions to demersal habitats have been rare. The lack of viable outgroup taxa, coupled with rapid evolution, which erases deep history, led to ancestral state analyses uninformative about the habitat state of the ancestor to Cephalopoda. Nevertheless, the deep ancestor of octopuses, true squids, loliginids, bobtail squids and cuttlefish (Coleoidea) received more support for a demersal state, with 19% of bootstrap trees recovering significantly greater proportional likelihood at this node. In contrast, we recover a pelagic existence for the ancestor of the octopuses and Vampyroteuthis (Octopodiformes), with 14% of bootstrapped topologies yielding significantly higher proportional likelihoods for this outcome compared to zero in favor of a demersal ancestor of Octopodiformes, a finding contradictory to previous hypotheses [30]. Young et al. [30] proposed that the immediate ancestor of Vampyroteuthis had an oral orientation that pre-adapted it for settling on the ocean floor, indicating a benthic ancestor. The subsequent transition to a demersal lifestyle in the cirrate/incirrate ancestor (33% of bootstrapped trees favoring the demersal state compared to zero favoring a pelagic state) represents the sole case of a pelagic-to- demersal transition that we recovered. The ancestral cirrate lineage also underwent a secondary transition to the pelagos, a scenario consistent with previous hypotheses based on morphology [30]. The ancestral decapodiform squid was likely demersal in nature, with recent transitions to the pelagic realm in some sepiolid lineages as well as an ancient transition to the pelagic environment occurring in the ancestor of the pelagic orders Oegopsida and Bathyteuthoidea, the latter finding reinforced by fossil data [31].

Figure 2
figure 2

Character state analyses on the 188-taxon tree summarized at the family-level. Ancestral states for habitat type and six morphological traits were reconstructed and tested for correlation on each of 1000 bootstrap replicates of our 188-taxon dataset. Extant morphological state presence for the six characters examined is shown in tip boxes as follows: accessory nidamental gland (yellow); branchial canal (pink); corneas (violet); right oviduct (teal); autogenic photophore (green); and bacteriogenic photophore (orange). The taxa found in pelagic habitats are shaded in grey while demersal taxa are white in the tree at right. More detail on character states and coding can be found in Additional file 2. Internal branches are shaded with: trait color, when ASR tests significantly favored the presence of the trait; white, when ASR tests significantly favored an absent state; black, when neither state fit significantly better.

Convergent evolution in Cephalopoda

Ancestral state reconstruction of six morphological characters on the 188-taxon dataset suggests that convergent evolution is a common phenomenon across cephalopods, as all six morphological characters we examined showed evidence of convergence (collapsed to family in Figure 2, see Figure 1 for complete topology). However, these characters were chosen based on a hypothesis that habitat has shaped their evolution, and therefore may not represent overall rates of convergent evolution in cephalopods.

ANG

An accessory nidamental gland (ANG) is present in squid families Idiosepiidae, Sepiolidae, Sepiidae, Spirulidae, Loliginidae, and Chtenopterygidae (e.g. [20, 32]). The ANG is present solely in females where it is thought to provide an additional protective coating of mucous on eggs. Female nautiloids and octopodiforms lack the accessory nidamental gland. Several octopods, as well as some oegopsids, lack a primary nidamental gland entirely, instead relying on oviducal glands on the oviducts that aid in egg secretion. Our reconstruction analysis finds the ancestral female decapodiform squid likely possessed an ANG (76% of bootstraps yielding significantly higher likelihoods for the present state). Although gain/loss rates for this trait are relatively low (not shown), losses likely occurred following transitions to pelagic lifestyles in the squid lineage.

Corneas

Our results are consistent with the hypothesis that the cornea is convergent between squids and octopods, a finding also supported by developmental data [33]. The decapodiform ancestor likely possessed a one-part cornea, which was lost with the transition to a pelagic lifestyle. In the octopod lineage, the basal pelagic groups Vampyroteuthis and Cirrata both lack a cornea. In the incirrate octopods, two pelagic clades possess reduced corneas or have lost corneas altogether, while all benthic octopods possess a fully formed two-part cornea.

Autogenic and bacteriogenic photophores

Bioluminescent organs, or photophores, have been acquired numerous times throughout Cephalopoda and rates of gain and loss suggest a very dynamic history (Figure2). Of all characters tested, autogenic photophores are the sole trait for which asymmetrical rates of gain and loss significantly improve the fit of the character model. Photophores appear to be lost at similarly moderate rates regardless of habitat, but show significant rates of gain in pelagic lineages, which is echoed by the structural diversity seen in these organs. Autogenic organs were likely present in the ancestor of the bathyteuthoids and oegopsids. If Spirulidae is sister to Bathyteuthoidea + Oegopsida, then autogenic bioluminescence may be synapomorphic among all three lineages. However, if Spirulidae is more closely related to sepioids, then we predict that autogenic photophores evolved separately in Spirulidae and in Bathyteuthoidea + Oegopsida. Autogenic bioluminescence in the Octopodiformes likely originated separately in Vampyroteuthis, the cirrate octopod genus Stauroteuthis, and in both genera (Bolitaena and Japatella) of the incirrate octopod family Bolitaenidae. While both forms of bioluminescence were included in subsequent habitat correlation analysis, bacteriogenic organs are not known from pelagic taxa.

Bacteriogenic photophores are known only from some demersal species from both sepiolid and loliginid lineages. Bacteriogenic photophores (those harboring symbiotic bacteria from the family Vibrionaceae, e.g. [3436]) occur in the demersal Loliginidae and Sepiolidae clades and appear to have evolved separately. Ancestral reconstruction supports the hypothesis of convergence for this character, as the MRCA of sepiolids and loliginids is predicted to have lacked a photophore.

Branchial canal

The branchial canal is present at the base of the gill lamellae of Vampyroteuthis, the incirrate octopods, Loliginidae, Bathyteuthoidea, and Oegopsida. Our findings suggest that this character is ancestral in the Octopodiformes (lost in the cirrates) and evolved convergently in the decapodiforms.

Right oviduct

Statistical analyses indicate that the right oviduct was likely present in the octopodiform ancestor and subsequently lost in the ancestral cirrate and gained independently within the decapodiform lineage. Ancestral state reconstruction also suggests that the right oviduct may have evolved convergently in the oegopsids and in Idiosepiidae, where is it highly reduced [27].

Correlation between phenotypes and marine habitats

Of the six characters we examined, three show correlation with marine habitat, supporting hypotheses that these characters are adaptive. A demersal lifestyle is correlated with the presence of accessory nidamental glands in females (p = 0.0147) and marginally correlated with the presence of corneas (p = 0.0747). In addition, we find a significant correlation between pelagic lifestyle and autogenic photophore presence (p = 0.0089). In contrast, branchial canal, right oviduct, and bacteriogenic photophores show no correlation with habitat type (Table1).

Table 1 Results of Pagel’s correlation analysis

Discussion

Convergence occurs frequently in the marine environment. In cases where distantly related taxa inhabit similar, or sympatric habitats, the presence of similar characteristics may be indicative of adaptation: in the open ocean, features such as a fusiform (or thunniform) body and silvering are present in a breadth of organisms to aid in hunting and predator avoidance (e.g. [16]). However, in cases where taxa are more closely related, it can be difficult to distinguish between phylogenetic and environmental influences. Comparative methods enable scientists to infer adaptations in the context of phylogeny, which provides insight into the interplay between phenotype and environment (e.g. [3741]) although some argue that these methods only take into account current usage and therefore cannot accurately represent historical processes [42, 43]. Martins [39] summarized the two general approaches that incorporate phylogeny when attempting to classify characters as adaptations: reconstructing ancestral character states (e.g. [42, 43] and testing for statistical correlation between habitat and morphology in a phylogenetic context (e.g. [38]). Both approaches have pitfalls: a large number of taxa are required for calculating statistical correlation, and a thorough understanding of single traits is required to reconstruct their ancestral states. Here, we overcome both of these challenges with a new, much larger, phylogenetic data set and a thorough examination of a small number of morphological characters, illustrating that cephalopods are an attractive group for testing hypotheses of convergence and adaptation in the marine environment. Modern cephalopods, not unlike extinct ammonoids [44] likely have also undergone several major habitat shifts and presently occupy a plethora of habitats, ranging from the intertidal, to the open ocean, to the hydrothermal vents. Transitions to such a diverse array of niches are echoed in the presence of a myriad of morphological characteristics, many of which appear to be convergent. Here, we explicitly tested hypotheses of morphological convergence in taxa living in similar marine environments, which required large-scale reconstruction of cephalopod phylogeny.

Phylogenetic relationships

Historically, the high rates of evolution evident in cephalopods coupled with their radiation into numerous habitats has confounded our ability to recover consistent phylogenetic relationships using either morphological [22, 27, 45, 46] or molecular [2224, 4753] data. Our analyses include the most comprehensive taxon sampling to date and provide a more rigorous assessment of phylogenetic relationships. Here, we recover the two major coleoid lineages with higher support values than in any previous study: Octopodiformes (all octopuses and Vampyroteuthis) and Decapodiformes (bobtail squids, cuttlefishes, pygmy squids, and true squids), as well as several lower-level relationships, discussed below.

Within Octopodiformes, we recovered monophyly for the higher-level groups Cirrata (the pelagic, finned octopods) and Incirrata (the benthic and pelagic octopods that lack fins), the pelagic octopod super-family Argonautoidea, and the pelagic octopod families Vitrellidonellidae and Bolitaenidae. The benthic family Octopodidae, which contains the majority of octopod species, was rendered non-monophyletic due to the internal placement of the pelagic vitrelledonellid + bolitaenid clade. The lack of monophyly for the family Octopodidae is not surprising. Although some morphological studies have suggested the Octopodidae to be monophyletic (e.g. [46]), current molecular and combined analyses have demonstrated that neither the family Octopodidae, nor the genus Octopus, are monophyletic (e.g. [48, 54]), suggesting a very dynamic evolutionary history. Our study finds support for the division of Octopodidae into several separate families, although more rigorous morphological and molecular examination is warranted before taxonomic changes can be made.

Within Decapodiformes, we recovered several important relationships, such as the monophyly of the pelagic lineages, Bathyteuthoidea, Oegopsida, Spirulida, and Myopsida, while the benthic lineages formed several distinct clades. While support values for the clade containing Bathyteuthoidea and Oegopsida are higher than in any previous study, the position of Spirulida (the ram’s horn squid) and Myopsida has low bootstrap support (56%). The benthic pygmy squids (Idiosepiidae) was the sister group to the rest of the Decapodiformes (with low support), while the benthic bobtail squids (Sepiolidae and Sepiadariidae) and cuttlefishes (Sepiidae) formed separate clades. Previous morphological analyses suggest that the sepioid orders are likely monophyletic (e.g. [26, 27, 55, 56]), but molecular data has been unable to recover consistent monophyly (e.g. [2224, 47, 49, 51, 53, 57]). The lack of resolution at the base of the decapodiform clade could also be due to the lack of a close living outgroup. The two coleoid lineages, Octopodiformes and Decapodiformes, are known to exhibit a substantial degree of genetic differentiation in their ribosomal genes [58, 59], which likely causes difficulties in rooting the Decapodiformes clade. The ambiguity in terms of clade placement among the benthic squids and Myopsida is not a new phenomenon. Naef [26] discussed the difficulties in placing myopsids relative to other squids, eventually placing loliginids (Myopsida) as the sister group to the pelagic oegopsids in the order Teuthoidea, citing characters such as the branchial canal, and gladius (the modified internal shell) as potential synapomorphies. Previous analyses have been ambiguous; some have placed myopsids near the sepioids (e.g. [22, 55]) and others as sister to the oegopsids/bathyteuthoids [24, 25]. It is possible that increased gene or taxon sampling could remedy the ambiguity over the placement of the loliginids. In particular, the inclusion of a second myopsid family, Australiteuthidae, which has yet remained uncollected for molecular data, could provide additional insight.

Character correlation and habitat

Habitat states were assigned based on the literature descriptions, which generally reference taxa as occupying a coastal benthic or an oceanic pelagic habitat (see the Methods section for details). Using this conservative binary coding scheme has several benefits. First, since habitat data are often limited to geographic coordinates and depth, we are able to confidently assign habitat states to more taxa in the dataset than would have been possible using multi-state habitat categories. For example, in many cases for nearshore coastal taxa the substantial data available in terms of their habitat and depth makes coding discrete states more straightforward (e.g. [60]). However, other groups were more challenging to classify, such as some cirrate octopods, which spend time in both the benthos and pelagos. Here, we condense the wide array of marine habitats into two states: neritic, bottom dwelling (demersal) and those inhabiting the water column in the open ocean (pelagic) in order to classify organisms as conservatively as possible to generate new hypotheses on character evolution and test for correlation.

Cephalopods possess numerous features that have been considered at one time or another to be convergent such as accessory nidamental glands, branchial canals, corneas, fins, internal shells, photophores, the ability to become transparent, egg brooding, and buoyancy. Several of these characteristics are not presently quantifiable due to a lack of data for living specimens or material: for example, we know very little about egg brooding in the true squids, although it has been known to occur in some pelagic taxa [61]. Transparency is of great interest to a multitude of researchers (e.g. [62]) but new insights are hindered by the technical difficulties in acquiring fresh specimens for most species. The presence of fins and internal shell morphology are also characters that may also be influenced by habitat, but were not quantified here, largely due to our inability to identify putative homologous structures. The fins of the deep-sea cirrate octopods contain “fin supports” which are one form of internal shell, while the fins of squids are muscular and are completely independent of the internal shell, making it difficult to determine independence and test for correlation, particularly in octopods. Furthermore, it has been extremely difficult to divide the internal shell into distinctive character states as many different morphologies exist: oegopsid, bathyteuthoid, and loliginid squid all have an internal shell in the form of a gladius, but its evolutionary history is unknown. Cuttlefishes have a cuttlebone, and sepiolids have what may be a highly reduced gladius in some species, but again, more work is needed to determine if separate character states are needed. The internal shell and fins represent a common problem with linking morphological and molecular data in cephalopods: one has to have a reasonable hypothesis of homology a priori. Here, we chose to focus on six characters whose morphology could be divided into distinct states with a substantial degree of confidence based on the authors’ examination of specimens and the literature. Of the six characters we examined, three correlate to habitat: accessory nidamental glands, corneas, and photophores. These characters are correlated when analyzed under alternate, pseudoreplicated topologies, indicating the results are robust to ambiguities in our phylogeny. Each character is discussed in detail below.

Accessory nidamental glands (ANG)

The importance to an individual’s fitness of protecting its reproductive investment is evident in the number of behavioral, chemical, and life history strategies that have evolved. Common solutions to avoiding predation and fouling on egg masses laid in benthic habitats include maternal care (e.g. [63, 64]), and chemical defenses (e.g. [65]). While female incirrate octopods typically guard and aerate egg masses they have deposited in dens, the demersal bobtail squids and cuttlefishes typically attach clutches, or festoons, of 50–400 eggs to a rocky or sandy substrate and provide no post-spawning care. The accessory nidamental gland, present in all demersal decapodiform females, secretes an additional layer of mucous on eggs prior to their attachment to the benthic substrate (e.g. [32, 6668]). The dense and diverse microbial community housed within the accessory nidamental glands and transmitted to the eggs may contain strains with antimicrobial properties, thus providing some degree of protection for vulnerable embryos against fouling organisms or egg predators [69, 70]. An ANG is also present the pelagic families Chtenopterygidae and Spirulidae and although they differ slightly in position when compared to the more benthic decapodiforms, Nesis [20] indicated that this structure is likely homologous. Adult chtenopterygids are meso- to bathypelagic, but it is unclear whether their eggs are deposited on the sediment, released into the midwater, or brooded as seen in some other pelagic squid groups [61]. We find evidence that the presence of an ANG correlates with a benthic lifestyle and was present in the ancestral decapodiform squid. In contrast, female octopods (both benthic and pelagic) lack a nidamental gland and an ANG, but provide substantial post-spawning care: females lay single eggs or strings in a den or shell, and tend to them constantly by maintaining a fresh, oxygenated, water supply and staving off predators until shortly before hatching, at which point the senescent female dies quickly (e.g. [64]).

Corneas

Our analyses indicate that the cephalopod cornea has evolved separately in the decapodiform and octopod lineages: all demersal bobtail squids, cuttlefishes, loliginids, and octopods possess a cornea, with reductions and/or complete loss occurring in pelagic forms. Young et al. [30] argued that the cornea was convergent and likely arose in response to similar selective pressures in the benthos, a finding further supported by developmental data [33]. Corneas are thought to serve several functions, ranging from keeping the eye socket free of sediment, to providing protection against ultraviolet light for taxa living in the photic zone (e.g. [71]). While it is possible that protection from ultraviolet light may play a minor role in corneal function, many benthic taxa never occupy the photic zone but still possess a fully formed cornea, indicating that a primary protective function is more likely.

Our results indicate that the ancestor of the decapodiforms likely lived in the demersal zone and possessed a cornea, which was lost as species shifted to a pelagic lifestyle. The association between habitat and corneal morphology is less clear in the octopod lineage because members of the pelagic incirrate superfamily Argonautoidea possess a cornea. The results of the discrete analysis suggest the origins of the octopod cornea correspond to the transition to a benthic habitat in the lineage ancestral to incirrates. The second pelagic octopod clade contains two families (Bolitaenidae and Vitrellidonellidae) wherein the cornea appears to be reduced or absent, which further supports the correlation between cornea and benthic habitat. One issue with evaluating corneal morphology in this group is that the pelagic octopods are extremely gelatinous animals, who deteriorate quickly upon fixation, making it difficult to determine if the cornea is simply reduced or completely absent.

Photophores

Bioluminescent organs throughout animals are extremely morphologically diverse (for reviews, see [10, 36, 72]) and can produce light via various endogenous biochemistries, so-called ‘autogenic’ organs, or via symbiotic bacteria harbored within specialized chambers, or ‘bacteriogenic’ organs [73, 74]. Similar to other organisms, both forms of photophores in cephalopods can be morphologically complex, relying on optical features such as reflectors, lenses and chromatophore-based light filters [75, 76]. We find that, in general, the ability to bioluminesce autogenically is highly correlated with pelagic existence in cephalopods (p = 0.0089).

Autogenic organs are found in many pelagic species of squids, octopods and in Vampyroteuthis (e.g. [36]). From study of two oegopsid species, the mechanism of light production appears to involve an intrinsic luciferase photoprotein reacting with dietary coelenterazine [7780]. The high variability in expression and morphology of autogenic photophores across and within families historically has been interpreted as evidence of separate origins [81, 82]. While experimentally establishing the adaptive value of these organs remains difficult owing to the challenges of sampling and maintaining deep-sea organisms, hypotheses abound. For instance, photophores frequently are located ventral to the eyes and opaque internal organs such as the ink sac, or dispersed across the surface of the skin, presumably serving to counter-illuminate the shadow cast by organs or body and thus provide a means of camouflage against deeper-dwelling predators [36, 83]. However, photophores also are found along the arms or encircling the mouth, allowing the animal to escape predation, attract or startle prey, or communicate with conspecifics [8487]. We find that autogenic photophores were likely present in some form in the ancestor of pelagic squid, and have undergone considerable diversification in terms of anatomical positioning, morphological complexity, and emission properties. Bioluminescence in octopod lineages appears to have evolved independently in three separate lineages: in Vampyroteuthis, within the cirrates, and in the incirrate family Bolitaenidae. The rapid transition rates of this trait suggest that any number of evolutionary mechanisms could be in play such as: selection acting either directly or indirectly on photophore expression and/or stochastic processes varying substantially through time and across lineages. Future work investigating to what extent cephalopods might rely on shared molecular mechanisms for bioluminescence may provide clues as to how such a complex trait could exhibit rapid rates of gain and loss.

Branchial canal, right oviduct, bacteriogenic photophore

Of the six characters examined, the branchial canal, right oviduct and bacteriogenic photophores did not correlate with habitat, although each has evolved multiple times within cephalopods. A branchial canal is present at the base of the gill lamellae in both benthic and pelagic cephalopods (absent from Cirrata) and is used to allow for passage of seawater into the gills [27]. Traditionally, the branchial canal was thought to be a synapomorphy uniting myopsids with the pelagic squids [27], but our ancestral state reconstructions suggest this state is convergent between Decapodiformes and Octopodiformes and within Decapodiformes. A right oviduct is present in both benthic and pelagic taxa, and ancestral state analysis predicts that this trait may be convergent between squids and octopods due to a loss early in decapodiform evolution. Also, if our prediction accurately reflects evolution, then the right oviduct found in Idiosepiidae and Oegopsida is the result of convergent evolution. While previous characterizations of the right oviduct in Idiosepiidae have suggested a highly modified morphology is due to a loss of function (e.g. [88]), our reconstruction analysis offers an alternative hypothesis of homoplasy and possible gain. Bacteriogenic photophores are far less common among cephalopods than compared to autogenic organs, and are found only in some members of the more benthic families Sepiolidae and Loliginidae. Although homology among bacteriogenic organs has been debated [89], our phylogenetic analyses indicate that bacteriogenic organs evolved separately in sepiolid and loliginid clades. Bacteriogenic photophores in both lineages typically develop on the ink sac as a bilobed structure, feature prominent crystalline lenses, and can recruit and culture related species of luminescent Vibrionaceae bacteria on the organ’s epithelium lining [34, 35, 90]. Evidence from the sepiolid Euprymna indicates that bacteriogenic organs may be used for counter-illumination [91]. Understanding both the molecular and ecological factors that have given rise to such a remarkable parallel will yield tremendous insight into how recurrent phenotypes evolve.

Conclusions

Convergent evolution occurs frequently in the marine environment, and cephalopods are no exception. Here, we identify convergent morphological characteristics in Cephalopoda, using both phylogenetic and statistical approaches and find correlation between habitat and morphology. The diverse nature of cephalopod habitat and morphology has interested researchers for hundreds of years, and yet the evolutionary history of many key characteristics has remained elusive. We provide conclusive evidence that accessory nidamental glands, corneas, and photophores are convergent and correlated with habitat. The accessory nidamental gland likely protects developing embryos from fouling microorganisms and is generally present only in decapodiforms known to attach eggs to the substrate. A cornea, one of the main components of the camera-type eye found in cephalopods and vertebrates, is present only in taxa that spend time associated with the benthos, and has evolved separately in the squid and octopod lineage albeit to serve similar purposes of protection against sediment and in some cases, ultraviolet light. Photophores, present in numerous marine phyla, are gained and lost many times throughout the cephalopod lineage, but are most commonly found in pelagic taxa.

The findings presented here in phylogenetics and character evolution create a foundation for understanding convergent evolution in the marine environment. The numerous biotic and abiotic factors that influence spatial and temporal expression for a given trait make studying convergence challenging, particularly in cephalopods where experimental studies linking function and habitat are often intractable. However, the advent of high-throughput transcriptome sequencing provides new means for generating hypotheses concerning the molecular basis for morphologically convergent evolution in experimentally challenging organisms like cephalopods. Future work utilizing these tools to study how convergent characters originate will provide us with new insight into how molecular processes such as gene duplication, gene sharing, and co-option can result in the evolution of similar morphological novelties.

Methods

Data assembly and alignment

Of the roughly 750 valid cephalopods species, sensu[92], many are rarely observed and poorly known to science, which has limited taxon sampling for molecular data. Here, all sequence data available for six nuclear genes (Histone H3a, octopine dehydrogenase, pax6, opsin, 18 S rRNA, and several regions of 28 S rRNA) and four mitochondrial genes (cytochrome c oxidase subunit I, cytochrome B, 12 S rRNA, and 16 S rRNA) were downloaded from GenBank and compiled into a local MySQL database yielding 409 taxonomic units with varying degrees of gene coverage (Additional file 3). Each locus was vetted for possible contamination: individual genes were aligned in Muscle under default parameters (see [93]) and analyzed in a maximum likelihood framework in RAxML v.7.0.4 [94] by generating a single tree under the GTR + GAMMA distribution model. Putative contaminants were identified as terminals with long branches relative to their sister taxa and also fell into incorrect clades (e.g. a squid in an octopus clade). These potential “contaminants” were then subjected to additional blast searches in GenBank, and the top 10 “hits” were evaluated. Any sequence identified as highly similar to another cephalopod was not considered to be a contaminant and was included in the final data matrix. Only sequences identified as a non-cephalopod (most commonly bacteria) were excluded from the analyses (see Additional file 4 for a list of identified contaminants).

Only one representative sequence for a given taxon was included for each gene in the final data matrix. When multiple copies of a gene were available on GenBank, as was often the case with COI, we included either the longest read (if length variation existed in GenBank files) or else the first read in the database. Sequences were aligned using default parameters in MAFFT [95], except for COI, which was aligned using MAFFT(−einsi) for improved accuracy. To further improve alignment accuracy, protein-coding sequences were then reverse transcribed using dna2pep.py, [96] and codon-aligned using TranslatorX [97] and RevTrans [96] to generate final DNA alignments that reflected codon positions. Lastly, we visually inspected each alignment for frame shift errors.

Molecular phylogenetic analysis

Taxon versus data coverage

To test for correlation between morphology and habitat, we required a well-sampled, robust phylogeny. The largest dataset with 408 taxa may not be optimal in this case because many terminals are represented by only a few genes and may be placed incorrectly due to the heterogeneity of the data matrix rather than to phylogenetic history [98, 99]. To identify the dataset that was both taxon-rich and relatively robust, nine datasets were generated and subjected to phylogenetic analyses. The largest dataset included all taxa with one or more genes present (408 taxa), while the smallest included only taxa that had at least nine loci present (17 taxa). The dataset we chose as the best topology to test for subsequent character evolution and correlation is comprised of 188 taxa represented by at least four genes (see Additional file 3 for taxa list). This dataset allows the inclusion of the greatest taxonomic diversity while maintaining robust support for all higher-level clades (greater than 75%). Furthermore, bootstrap support values for these clades evaluated on datasets with greater data matrix density (i.e., more genes, fewer taxa) were never less than 75% (Table2).

Table 2 Bootstrap support for key nodes across increasingly refined datasets

Prior to formal analyses, ModelTest [100] was used to determine which model best fit the data for each partition. Each dataset was analyzed in RAxML 7.2.8 [94] under the GTR + GAMMA model (−N 10), and support was calculated with 1,000 replicates of bootstrap resampling. To identify the dataset to use for tests of character correlation, we measured clade support and compared values (Table2) for eight nodes commonly examined in taxonomic and molecular studies: Octopodiformes (V. infernalis and all octopods), Octopoda, Incirrata, Cirrata, Decapodiformes, Loliginidae, Bathyteuthoidea + Oegopsida, and Oegopsida. The data set that contained at least three representative taxa in each node and support values above 80% was used for all subsequent analyses.

Data partitions and analyses

Given the diverse nature of the genes in our data matrix, using the same likelihood parameter values for the entire dataset is unlikely to provide a good fit (e.g. [101, 102]). To reduce the risks associated with both over- and under-parameterization, we tested partitioning schemes most likely to represent groupings of genes that may have evolved under similar conditions (Table1). We computed log-likelihoods for each of the sixteen partition strategies using a 100 replicate maximum likelihood search (RAxML v. 7.2.8) and identified the partitioning strategy that best fit the data by comparing the log-likelihoods from each of the lesser-partitioned analyses to the most heavily partitioned (sixteen) scheme using a modified BIC test [103] as implemented by McGuire et al. [104]. This modified version of the BIC test holds the number of branch lengths constant across alternative partitioning strategies, therefore returning identical estimates to the standard BIC calculations [104].

We analyzed the final data matrix (188 taxonomic units, each represented by four or more loci) under a likelihood setting in RAxML 7.2.8 with 300 replicate tree searches under GTR + GAMMA. Each gene and the 3rd position codon sites were treated as separate partitions (16-partition scheme, Table1). All species of Nautilus were listed as outgroup taxa. We then evaluated support using 1000 nonparametric bootstrap replicates (−f b) in RAxML. For comparison, the above tree search was also conducted on the complete data set (Additional file 1, Figure1).

The partitioning strategy that best fit the data (16-part, Additional file 1 Table S1) allowed for rate parameters for the sites of each ribosomal gene as well as the non-synonymous (1st/2nd codon positions) and synonymous (3rd codon positions) sites for each coding gene to be estimated independently [105, 106]. It is worth noting that the 6-partition strategy (in which nuclear, mitochondrial genes, ribosomal, 1st/2nd coding positions and 3rd codon positions are differentiated) yielded likelihoods that were close to, but still less likely than the fully partitioned scheme and the BIC evaluation confirmed that the additional parameters in the 16-part scheme significantly improved the fit of this data set.

While additional phylogenetic analyses (not shown) were carried out in PhyloBayes (v2.3) using a heterogeneous-site model (GTRCAT), our dataset failed to converge after surpassing the recommended maximum of 10,000 MCMC cycles on each of 6 chains and running for over 60 days [106]. Due to the lack of convergence in the posterior distribution, we carried out likelihood-based character analyses using our best ML tree (Figure1) and ML bootstrap pseudoreplicates.

Character evolution and correlation with habitat

To attempt to identify characters that have arisen via convergent evolution in response to similar selective pressures, we performed ancestral states reconstructions using our 188-taxon topology. We created a binary character matrix for habitat (benthic or pelagic) and for six features of presumed taxonomic importance (e.g. [26, 27, 30]) that can be putatively associated with either demersal or pelagic lifestyles [30]: accessory nidamental gland, branchial canal, cornea, right oviduct, autogenic and bacteriogenic photophores. With the exception of the accessory nidamental glands and bacteriogenic photophores, the remaining characters are present in both squid and octopod lineages, thus increasing our ability to test for correlation. Character coding was done primarily through literature, although all pelagic octopod taxa were re-examined at the Santa Barbara Museum of Natural History to confirm corneal morphology.

For the habitat data, we chose to cluster habitats into two primary categories (demersal/pelagic) to be able to conduct statistical tests of correlation, and to identify broad patterns and hypotheses that can then be further examined in on a much finer scale. All habitat data for Decapodiformes was coded from the literature [60, 107]: all Octopodiformes data was coded by FGH. There is some inherent risk in subdividing the expanse of the marine environment into two discrete habitats, as each macro habitat encompasses a number of smaller ecosystems. There are numerous habitats that can fall within this spectrum, and there is some inherent risk as each macro habitat encompasses a number of smaller ecosystems. Although there are many potential subdivisions within the demersal/pelagic realms, there are still some general trends seen across the entire habitat in terms of selective pressures. For example, all taxa living in demersal environments, must to some degree, compete for space. Alternately, taxa inhabiting the pelagic realm are more impacted by the three-dimensional nature of their surroundings and must be able to respond to stimuli from numerous directions.

Likelihood-based ancestral state reconstructions were carried out in BayesTraits [108] on our 188-taxon topology for habitat type (benthic/pelagic), and for characters (present/absent): accessory nidamental gland, cornea, bacteriogenic photophore, autogenic photophore, branchial canal, and right oviduct. Prior to ASR analysis, each trait was fit to a single-transition rate and asymmetrical rate model. The best-fitting model was determined with the likelihood-ratio test.

We calculated the proportional likelihoods of either ancestral state on each of the 1000 bootstrap ML trees generated from our dataset to incorporate topological uncertainty. We specifically evaluated character state likelihoods for the high-level nodes of the best ML tree recovered in at least 90% of bootstrap replicates. For nodes that received low branch support (e.g., basal Decapodiformes), we evaluated ancestral state likelihoods for their most recent common ancestor (MRCA), such as for Oegopsida + Loliginidae. As it remains unclear whether the MRCA of Oegopsida + Loliginidae lived earlier or later than the MCRA of Sepiolida + Loliginidae, we estimated the character state likelihoods for both ancestors (Figure2, Additional file 1 Figure S2). To determine the extent to which one character state received significantly more support over the other state at a particular node, we calculated the number of bootstrap replicates for which the ratio between the two proportional likelihood scores (e.g., L present/L absent) was greater than 7.4 [109]. In addition, for each node of interest, we calculated the likelihoods across 1000 bootstrap trees fixed for either state using BayesTraits’ ‘fossil’ command. State models which yielded log-likelihood scores having two or more units higher than the other state model are regarded as a significantly better fit [110, 111].

To test whether certain cephalopod characteristics are correlated with demersal or pelagic existence, we used Pagel’s correlations as implemented in BayesTraits [108] on our 188-taxon topology. We tested for the correlation between habitat and each binary trait by comparing the fit of a dependent and independent character model on bootstrapped ML trees. Under the independent model, the number and timing of gains and losses for each trait (e.g., habitat and the morphological trait) occur independently, whereas the dependent model provides additional parameters that affect the probability of a character’s transitioning on a branch depending on the state of the other character. To account for the phylogenetic uncertainty in our dataset, we fit these models for each character set on 1000 bootstrap ML topologies. We then calculated median log-likelihoods for each model and used the log-ratio test to determine which characters are significantly correlated.

Authors’ information

Annie R Lindgren and Molly S Pankey are co-first authors.

References

  1. Carrol R: Vertebrate Paleontology and Evolution. 1988, W.H. Freeman and Company, New York

    Google Scholar 

  2. Bell MA: Lateral plate evolution in the threespine stickleback: getting nowhere fast. Genetica. 2001, 112–113: 445-461.

    PubMed  Google Scholar 

  3. Bell MA, Orti G, Walker JA, Koenings JP: Evolution of Pelvic Reduction in Threespine Stickleback Fish: A Test of Competing Hypotheses. Evolution. 1993, 47: 906-914. 10.2307/2410193.

    Google Scholar 

  4. Shapiro MD, Summers BR, Balabhadra S, Aldenhoven JT, Miller AL, Cunningham CB, Bell MA, Kingsley DM: The Genetic Architecture of Skeletal Convergence and Sex Determination in Ninespine Sticklebacks. Curr Biol. 2009, 19: 1140-1145. 10.1016/j.cub.2009.05.029.

    PubMed  CAS  PubMed Central  Google Scholar 

  5. Mobley KB, Lussetti D, Johansson F, Englund G, Bokma F: Morphological and genetic divergence in Swedish postglacial stickleback (Pungitius pungitius) populations. BMC Evol Biol. 2011, 11: 287-10.1186/1471-2148-11-287.

    PubMed  PubMed Central  Google Scholar 

  6. Jones FC, Grabherr MG, Chan YF, Russell P, Mauceli E, Johnson J, Swofford R, Pirun M, Zody MC, White S, Birney E, Searle S, Schmutz J, Grimwood J, Dickson MC, Myers RM, Miller CT, Summers BR, Knecht AK, Brady SD, Zhang H, Pollen AA, Howes T, Amemiya C, Lander ES, Palma FD, Lindblad-Toh K, Kingsley DM, Team BIGSP& WGA: The genomic basis of adaptive evolution in threespine sticklebacks. Nature. 2012, 484: 55-61. 10.1038/nature10944.

    PubMed  CAS  PubMed Central  Google Scholar 

  7. Chen L, DeVries AL, Cheng CC: Evolution of antifreeze glycoprotein gene from a trypsinogen gene in Antarctic notothenioid fish. PNAS. 1997, 94: 3811-3816. 10.1073/pnas.94.8.3811.

    PubMed  CAS  PubMed Central  Google Scholar 

  8. Herring PJ: Bioluminescent signals and the role of reflectors. J. Optics A-Pure and App. Optics. 2000, 2: R29-R38. 10.1088/1464-4258/2/6/202.

    Google Scholar 

  9. Widder EA: Bioluminesence. In Adaptive Mechanisms in the Ecology of Vision. 1999, Kluwer Academic Publishers, Dordrecht, the Netherlands, 555-581.

    Google Scholar 

  10. Haddock SHD, Moline MA, Case JF: Bioluminescence in the Sea. Ann. Rev. Mar. Sci. 2010, 2: 443-493. 10.1146/annurev-marine-120308-081028.

    PubMed  Google Scholar 

  11. Donley JM, Sepulveda CA, Konstantinidis P, Gemballa S, Shadwick RE: Convergent evolution in mechanical design of lamnid sharks and tunas. Nature. 2004, 429: 61-65. 10.1038/nature02435.

    PubMed  CAS  Google Scholar 

  12. Childress JJ: Are there physiological and biochemical adaptations of metabolism in deep-sea animals? Trends Ecol. Evol. 1995, 10: 30-36.

    CAS  Google Scholar 

  13. Seibel BA, Carlini DB: Metabolism of pelagic cephalopods as a function of habitat depth: a reanalysis using phylogenetically independent contrasts. Biol Bull. 2001, 201: 1-5. 10.2307/1543519.

    PubMed  CAS  Google Scholar 

  14. McFall-Ngai MJ: Crypsis in the pelagic environment. Am. Zool. 1990, 30: 175-188.

    Google Scholar 

  15. Hanlon RT, Messenger JB: Cephalopod Behavior. 1996, Cambridge University Press, Cambridge

    Google Scholar 

  16. Packard A: Cephalopods and fish: The limits of convergence. Biol Rev. 1972, 47: 241-307. 10.1111/j.1469-185X.1972.tb00975.x.

    CAS  Google Scholar 

  17. Finn JK, Tregenza T, Norman MK: Defensive tool use in a coconut-carrying octopus. Current Biol. 2009, 19: R1069-R1070. 10.1016/j.cub.2009.10.052.

    CAS  Google Scholar 

  18. Hanlon RT, Chiao CC, Mäthger LM, Barbosa A, Buresch KC, Chubb C: Cephalopod dynamic camouflage: bridging the continuum between background matching and disruptive coloration. Philos. Trans. R. Soc. London B. 2009, 364: 429-437. 10.1098/rstb.2008.0270.

    CAS  Google Scholar 

  19. Norman M, Finn JK, Tregenza T: Dynamic mimicry in an Indo-Malayan octopus. Proc. R. Soc. London B. 2001, 268: 1755-1758. 10.1098/rspb.2001.1708.

    CAS  Google Scholar 

  20. Nesis A: Cephalopods of the World. 1982, T.F.H Publications, New Jersey

    Google Scholar 

  21. Strugnell JM, Nishiguchi MK: Molecular phylogeny of coleoid cephalopods (Mollusca: Cephalopoda) inferred from three mitochondrial and six nuclear loci: a comparison of alignment, implied alignment and analysis methods. J. Moll. Stud. 2007, 73: 399-410. 10.1093/mollus/eym038.

    Google Scholar 

  22. Lindgren AR, Giribet G, Nishiguchi MK: A combined approach to the phylogeny of Cephalopoda (Mollusca). Cladistics. 2004, 20: 454-486. 10.1111/j.1096-0031.2004.00032.x.

    Google Scholar 

  23. Strugnell JM, Norman M, Jackson J, Drummond AJ, Cooper A: Molecular phylogeny of coleoid cephalopods (Mollusca: Cephalopoda) using a multigene approach; the effect of data partitioning on resolving phylogenies in a Bayesian framework. Mol. Phylogenet. Evol. 2005, 37: 426-441. 10.1016/j.ympev.2005.03.020.

    PubMed  CAS  Google Scholar 

  24. Lindgren AR: Molecular inference of phylogenetic relationships among Decapodiformes (Mollusca: Cephalopoda) with special focus on the squid order Oegopsida. Mol. Phylogenet. Evol. 2010, 56: 77-90. 10.1016/j.ympev.2010.03.025.

    PubMed  Google Scholar 

  25. Allcock L, Cooke IR, Strugnell JM: What can the mitochondrial genome reveal about higher-level phylogeny of the molluscan class Cephalopoda?. Zool. J. Linnean Soc. 2011, 161: 563-586.

    Google Scholar 

  26. Naef A: Die Cephalopoden. Fauna e flora del Golfo di Napoli. 1923, Israel Programs for Scientific Translations, Jerusalem, 1-

    Google Scholar 

  27. Young RE, Vecchione M: Analysis of morphology to determine primary sister taxon relationships within coleoid cephalopods. Am. Malacological Bull. 1996, 12: 91-112.

    Google Scholar 

  28. Norman MK, Hochberg FG: The current state of octopus taxonomy. Phuket Mar. Biol. Cent. Res. Bull. 2005, 66: 127-154.

    Google Scholar 

  29. Whitfield JB, Lockhart PJ: Deciphering ancient rapid radiations. Trends Ecol. Evol. 2007, 22: 258-265. 10.1016/j.tree.2007.01.012.

    PubMed  Google Scholar 

  30. Young RE, Vecchione M, Donovan D: The evolution of coleoid cephalopods and their present biodiversity and ecology. South African Journal of Marine Science. 1998, 20: 393-420. 10.2989/025776198784126287.

    Google Scholar 

  31. Kröger B, Servais T, Zhang Y: The origin and initial rise of pelagic cephalopods in the Ordovician. PLoS One. 2009, 4: e7262-10.1371/journal.pone.0007262.

    PubMed  PubMed Central  Google Scholar 

  32. Bloodgood R: Squid accessory nidamental gland – ultrastructure and association with bacteria. Tissue and Cell. 1977, 9: 197-208. 10.1016/0040-8166(77)90016-7.

    PubMed  CAS  Google Scholar 

  33. Naef A: Die Cephalopoden. Embryologie. 1928, Amerind Publishing Company, New Dehli, 2-

    Google Scholar 

  34. Fidiopiastis P, Boletzky SV, Ruby EG: A new niche for Vibrio logei, the predominant light organ symbiont of squids in the genus Sepiola. J Bacteriol. 1998, 180: 59-64.

    Google Scholar 

  35. McFall-Ngai MJ, Montgomery M: The anatomy and morphology of the adult bacterial light organ of Euprymna scolopes Berry (Cephalopoda: Sepiolidae). Biol Bull. 1990, 179: 332-339. 10.2307/1542325.

    Google Scholar 

  36. Herring PJ: Luminescent organs. The Mollusca. Edited by: Wilbur KM. 1988, Academic, San Diego, 449-485. 11

    Google Scholar 

  37. Coddington JA: Review of: The Explanation of Organic Diversity: The Comparative Method and Adaptations for Mating. Cladistics. 1985, 1: 102-107.

    Google Scholar 

  38. Harvey PH, Pagel MD: The Comparative Method in Evolutionary Biology. 1991, Oxford University Press, London

    Google Scholar 

  39. Martins EP: Adaptation and the comparative method. Trends Ecol. Evol. 2000, 15: 296-299. 10.1016/S0169-5347(00)01880-2.

    PubMed  Google Scholar 

  40. Westoby M, Leishman MR, Lord JM: On misinterpreting the phylogenetic correction. J Ecol. 2007, 83: 531-553.

    Google Scholar 

  41. Price T: Correlated evolution and independent contrasts. Philos. Trans. R. Soc. London Ser. B. 1997, 352: 519-552. 10.1098/rstb.1997.0036.

    CAS  Google Scholar 

  42. Baum DA, Larson A: A phylogenetic methodology for studying character macroevolution. Syst. Zool. 1991, 40: 1-28. 10.2307/2992218.

    Google Scholar 

  43. Leroi AM, Rose MR, Lauder GV: What does the comparative method reveal about adaptation? Am. Nat. 1994, 143: 381-402. 10.1086/285609.

    Google Scholar 

  44. Monnet C, Baets K, Klug C: Parallel evolution controlled by adaptation and covariation in ammonoid cephalopods. BMC Evol Biol. 2011, 11: 115-10.1186/1471-2148-11-115.

    PubMed  PubMed Central  Google Scholar 

  45. Voss GL: Evolution and phylogenetic relationships of deep-sea octopods (Cirrata and Incirrata). In The Mollusca: Paleontology and neontology of cephalopods. 1988, Academ, London, 12-

    Google Scholar 

  46. Voight JR: Cladistic analysis of the octopods based on morphological characters. J. Moll. Stud. 1997, 63: 311-325. 10.1093/mollus/63.3.311.

    Google Scholar 

  47. Carlini DB, Graves JE: Phylogenetic analysis of cytochrome c oxidase subunit I sequences to determine higher-level relationships within the coleoid cephalopods. Bull Mar Sci. 1999, 64: 57-76.

    Google Scholar 

  48. Carlini DB, Young RE, Vecchione M: A molecular phylogeny of Octopoda (Mollusca: Cephalopoda) evaluated in light of morphological evidence. Molec. Phylogenet. Evol. 2001, 21: 388-397. 10.1006/mpev.2001.1022.

    PubMed  CAS  Google Scholar 

  49. Bonnaud L, Boucher-Rodoni R, Monnerott M: Relationship of some coleoid cephalopods established by 3’ end of the 16 S rDNA and cytochrome oxidase III gene sequence comparison. Am. Malacological Bull. 1996, 12: 87-90.

    Google Scholar 

  50. Bonnaud L, Boucher-Rodoni R, Monnerott M: Phylogeny of cephalopods inferred from mitochondrial DNA sequences. Mol. Phylogenet. Evol. 1997, 7: 44-54. 10.1006/mpev.1996.0366.

    PubMed  CAS  Google Scholar 

  51. Akasaki T, Nikaido K, Tsuchiya K, Segawa S, Hasegawa M, Okada N: Extensive mitochondrial gene arrangements in coleoid Cephalopoda and their phylogenetic implications. Mol. Phylogenet. Evol. 2006, 38: 648-658. 10.1016/j.ympev.2005.10.018.

    PubMed  CAS  Google Scholar 

  52. Yokobori S, Fukuda N, Nakamura M, Aoyama T, Oshima T: Long-term conservation of six duplicated structural genes in cephalopod mitochondrial genomes. Mol Biol Evol. 2004, 21: 2034-2046. 10.1093/molbev/msh227.

    PubMed  CAS  Google Scholar 

  53. Yokobori S, Lindsay M, Yoshida M, Tsuchiya K, Yamagashi A, Maruyama T, Oshima T: Mitochondrial genome structure and evolution in the living fossil vampire squid, Vampyroteuthis infernalis, and extant cephalopods. Mol. Phylogenet. Evol. 2007, 44: 898-910. 10.1016/j.ympev.2007.05.009.

    PubMed  CAS  Google Scholar 

  54. Strugnell JM, Rogers AD, Prodöhl PA, Collins MA, Allcock AL: The thermohaline expressway: the Southern Ocean as a centre of origin for deep‒sea octopuses. Cladistics. 2008, 24: 853-860. 10.1111/j.1096-0031.2008.00234.x.

    Google Scholar 

  55. Berthold T, Engeser T: Phylogenetic analysis and systematization of the Cephalopoda (Mollusca). Verh. Naturwiss. Ver. Hamburg. 1987, 29: 187-220.

    Google Scholar 

  56. Vecchione M, Young RE, Donovan DT, Rodhouse PG: Reevaluation of coleoid cephalopod relationships based on modified arms in the Jurassic coleoid Mastigophora. Lethaia. 1999, 32: 113-118.

    Google Scholar 

  57. Carlini DB, Reece KB, Graves JE: Actin gene family evolution and the phylogeny of coleoid cephalopods (Mollusca: Cephalopoda). Mol Biol Evol. 2000, 17 (9): 1353-1370. 10.1093/oxfordjournals.molbev.a026419.

    PubMed  CAS  Google Scholar 

  58. Bonnaud L, Saihi A, Boucher-Rodoni R: Are 28s rDNA and 18s rDNA informative for cephalopod phylogeny?. Bull Mar Sci. 2002, 71 (1): 197-208.

    Google Scholar 

  59. Lindgren AR, Daly M: The impact of length-variable data and alignment criterion on the phylogeny of Decapodiformes (Mollusca: Cephalopoda). Cladistics. 2007, 23: 464-476. 10.1111/j.1096-0031.2007.00160.x.

    Google Scholar 

  60. Jereb P, Roper CFE: Cephalopods of the World. An annotated and illustrated catalogue of cephalopod species known to date: Chambered nautiluses and sepioids. In FAO Species Catalogue for Fisheries Purposes. 2005, Food and Agriculture Organization of the United Nations, Rome, 1-

    Google Scholar 

  61. Seibel BA, Hochberg FG, Carlini DB: Life history of Gonatus onyx (Cephalopoda: Teuthoidea): deep-sea spawning and post-spawning egg care. Mar Biol. 2000, 137: 519-526. 10.1007/s002270000359.

    Google Scholar 

  62. Zylinski S, Johnsen S: Mesopelagic cephalopods switch between transparency and pigmentation to optimize camouflage in the deep. Current Biol. 2011, 21: 1937-1941. 10.1016/j.cub.2011.10.014.

    CAS  Google Scholar 

  63. Boletzky S: Embryonic development of cephalopods at low temperatures. Antarct Sci. 1994, 6: 139-142.

    Google Scholar 

  64. Voight JR, Grehan AJ: Egg brooding by deep-sea octopuses in the North Pacific Ocean. Biol Bull. 2000, 198: 94-100. 10.2307/1542807.

    PubMed  CAS  Google Scholar 

  65. Bekendorff K, Davis AR, Brenner JB: Chemical defense in the egg massess of benthic invertebrates: an assessment of antibacterial activity in 39 mollusks and 4 polychaetes. J Invertebr Pathol. 2001, 78: 109-118. 10.1006/jipa.2001.5047.

    Google Scholar 

  66. Barbieri E, Barry K, Child A, Wainwright N: Antimicrobial activity in the microbial community of the accessory nidamental gland and egg cases of Loligo pealei (Cephalopoda: Loliginidae). Biol Bull. 1997, 193: 275-276.

    Google Scholar 

  67. Barbieri E, Gulledge J, Moser D, Chieng CC: New evidence for bacterial diversity in the accessory nidamental gland of the squid Loligo pealei. Biol Bull. 1996, 191: 316-317.

    Google Scholar 

  68. Kaufman MY, Ikeda Y, Patton C, Van Dykhuizen G, Epel D: Bacterial symbionts colonize the accessory nidamental gland of the squid Loligo opalescens via horizontal transmission. Biol Bull. 1998, 194: 36-43. 10.2307/1542511.

    Google Scholar 

  69. Pichon D, Gaia V, Norman M, Boucher R: Phylogenetic diversity of epibiotic bacteria in the accessory nidamental glands of squids (Cephalopoda: Loliginidae and Idiosepiidae). Mar Biol. 2005, 147: 1323-1332. 10.1007/s00227-005-0014-5.

    Google Scholar 

  70. Nair P, Sherief PM: Antibacterial activity in the accessory nidamental gland extracts of the Indian squid, Loligo duvauceli Orbigny. Indian J. Mar. Sci. 2010, 39: 100-104.

    Google Scholar 

  71. Piatigorsky J: Enigma of the abundant water-soluble cytoplasmic proteins of the cornea: the “refraction” hypothesis. Cornea. 2001, 20: 853-858. 10.1097/00003226-200111000-00015.

    PubMed  CAS  Google Scholar 

  72. Harvey EN: Bioluminescence. 1952, Academic, New York

    Google Scholar 

  73. Kishitani T: Preliminary report on the luminous symbiosis in Sepiola birostrata Sasaki. Proc. Imp. Acad. Jpn. 1928, 4: 393-396.

    Google Scholar 

  74. Kishitani T: L’tétude de I’organe photognène du Loligo edulis Hoyle (Notes préliminaires). Proc. Imp. Acad. Jpn. 1928, 4: 609-612.

    Google Scholar 

  75. Arnold J, Young R: Ultrastructure of a cephalopod photophore. I. Structure of the photogenic tissue. Biol Bull. 1974, 174: 507-521.

    Google Scholar 

  76. Young RE, Arnold J: The functional morphology of a ventral photophore from the mesopelagic squid Abralia trigonura. Malacologia. 1982, 23: 135-163.

    Google Scholar 

  77. Young R, Roper CF, Mangold K, Leisman G, Hochberg FG: Luminescence from non-bioluminescent tissues in oceanic cephalopods. Mar Biol. 1979, 53: 69-77. 10.1007/BF00386530.

    Google Scholar 

  78. Takahashi H, Isobe M: Photoprotein of luminous squid, Symplectoteuthis oualaniensis and reconstruction of the luminous system. Chem Lett. 1994,  : 843-846.

  79. Isobe M, Kuse M, Tanji N, Matsuda T: Cysteine-390 is the binding site of luminous substance with symplectin, a photoprotein from Okinawan squid, Symplectoteuthis oualaniensis. Proc. Jpn. Acad. Ser. B Phys. Biol. Sci. 2008, 84: 386-932. 10.2183/pjab.84.386.

    PubMed  CAS  PubMed Central  Google Scholar 

  80. Tsujii F: Bioluminescence reaction catalyzed by membrane-bound luciferase in the “firefly squid” Watasenia scintillans. BBA Biomembranes. 2002, 1564: 189-917. 10.1016/S0005-2736(02)00447-9.

    Google Scholar 

  81. Herring PJ: Luminescence in cephalopods and fish. In Biology of Cephalopods. Symposia of the Zoological Society of London. 1977, Academic, London, 449-489.

    Google Scholar 

  82. Herring PJ, Dilly PN, Cope C: Different types of photophore in the oceanic squids Octopoteuthis and Taningia (Cephalopoda: Octopoteuthidae). J. Zool. 1992, 227: 479-491. 10.1111/j.1469-7998.1992.tb04408.x. Lond

    Google Scholar 

  83. Young RE: Oceanic bioluminescence - an overview of general functions. Bull Mar Sci. 1983, 33: 829-845.

    Google Scholar 

  84. Robison BH, Young RE: Bioluminescence in pelagic octopods. Pac Sci. 1981, 35: 39-44.

    Google Scholar 

  85. Young RE, Seapy R, Mangold KM, Hochberg FG: Luminescent flashing in the midwater squids Pterygioteuthis microlampas and P. giardi. Mar Biol. 1982, 69: 299-308. 10.1007/BF00397495.

    Google Scholar 

  86. Robison BH, Reisenbichler KR, Hunt JC, Haddock SHD: Light production by the arm tips of the deep-sea cephalopod Vampyroteuthis infernalis. Biol Bull. 2003, 205: 102-109. 10.2307/1543231.

    PubMed  Google Scholar 

  87. Bush SL, Robison BH, Caldwell RH: Behaving in the dark: loctomotor, chromatic, postural, and bioluminescent behaviors of the deep-sea squid Octopoteuthis deletron Young 1972. Biol Bull. 2009, 216: 7-22.

    PubMed  Google Scholar 

  88. Young RE, Vecchione M, Sweeney MJ: Bathyteuthoidea:. The Tree of Life Web Project. -[tolweb.org/Bathyteuthoidea]

  89. Nishiguchi MK, Lopez J, Boletzky SV: Enlightenment of old ideas from new investigations: more questions regarding the evlution of bacteriogenic light organs in squids. Evol. Devel. 2004, 61: 41-49.

    Google Scholar 

  90. Guerrero-Ferreira R, Nishiguchi MK: Biodiversity among luminescent symbionts from squid of the genera Uroteuthis, Loliolus and Euprymna (Mollusca: Cephalopoda). Cladistics. 2007, 23: 497-506. 10.1111/j.1096-0031.2007.00155.x.

    PubMed  PubMed Central  Google Scholar 

  91. Jones B, Nishiguchi MK: Counterillumination in the Hawaiian bobtail squid, Euprymna scolopes Berry (Mollusca: Cephalopoda). Mar Biol. 2004, 144: 1151-1155. 10.1007/s00227-003-1285-3.

    Google Scholar 

  92. Sweeney MJ, Roper CF: Classification, type localities and type repositories of recent Cephalopoda. Smithsonian Cont. Zool. 1998, 586: 561-599.

    Google Scholar 

  93. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nuc. Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.

    CAS  Google Scholar 

  94. Stamatakis A: RAxML-VI-HPC: Maximum Likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.

    PubMed  CAS  Google Scholar 

  95. Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nuc. Acids Res. 2002, 30: 3059-3066. 10.1093/nar/gkf436.

    CAS  Google Scholar 

  96. Wernersson R: Virtual Ribosome - a comprehensive translation tool wth support for sequence feature integration. Nuc. Acids Res. 2006, 34: W385-W388. 10.1093/nar/gkl252.

    CAS  Google Scholar 

  97. Abscal F, Zardoya R, Telford MJ: TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translations. Nuc. Acids Res. 2010, 38: W7-W13. 10.1093/nar/gkq291.

    Google Scholar 

  98. Sanderson M, Shaffer HB: Troubleshooting molecular phylogenetic analyses. Ann. Rev. Ecol. Syst. 2002, 33: 49-72. 10.1146/annurev.ecolsys.33.010802.150509.

    Google Scholar 

  99. Thomson RC, Shaffer HB: Sparse supermatrices for phylogenetic inference: Taxonomy, alignment, rogue taxa, and the phylogeny of living turtles. Syst Biol. 2010, 59: 42-58. 10.1093/sysbio/syp075.

    PubMed  Google Scholar 

  100. Posada D, Crandall K: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.

    PubMed  CAS  Google Scholar 

  101. Nylander JA, Ronquist AF, Huelsenbeck JP, Nieves Aldrey JL: Bayesian phylogenetic analysis of combined data. Syst Biol. 2004, 53: 47-67. 10.1080/10635150490264699.

    PubMed  Google Scholar 

  102. Brandley MC, Schmitz A, Reeder TW: Partitioned Bayesian analyses, partition choice, and the phylogenetic relationships of scinid lizards. Syst Biol. 2005, 54: 373-390. 10.1080/10635150590946808.

    PubMed  Google Scholar 

  103. Shwarz G: Estimating the dimensions of a model. Ann Stat. 1978, 6: 461-464. 10.1214/aos/1176344136.

    Google Scholar 

  104. McGuire JA, Witt CC, Altshuler DL, Remsen JV: Phylogenetic systematics and biogeography of hummingbirds: Bayesian and Maximum Likelihood analyses of partitioned data and selection of an appropriate partitioning strategy. Syst Biol. 2007, 56: 837-856. 10.1080/10635150701656360.

    PubMed  CAS  Google Scholar 

  105. Lartillot N, Phillipe H: A Bayesian mixsture model for across-site heterogeneities in the amino-acid replacement process. Mol Biol Evol. 2004, 21: 1095-1109. 10.1093/molbev/msh112.

    PubMed  CAS  Google Scholar 

  106. Lartillot N, Phillipe H: Computing Bayes factors using thermodynamic integration. Syst Biol. 2006, 55: 195-207. 10.1080/10635150500433722.

    PubMed  Google Scholar 

  107. Jereb P, Roper CFE: Cephalopods of the World. An annotated and illustrated catalogue of cephalopod species known to date: Myopsid and Oegopsid Squids. 2010, Food and Agriculture Organization of the United Nations, Rome, 2-

    Google Scholar 

  108. Pagel M, Meade A: Bayesian analysis of correlated evolution of discrete characters by reversible-jump Markov chain Monte Carlo. Am Nat. 2006, 167: 808-825. 10.1086/503444.

    PubMed  Google Scholar 

  109. Schluter D, Price T, Mooers AØ, Ludwig D: Likelihood of ancestor states in adaptive radiation. Evolution. 1997, 51: 1699-1711. 10.2307/2410994.

    Google Scholar 

  110. Edwards AFW: Likelihood. 1972, Cambridge University Press, Cambridge

    Google Scholar 

  111. Pagel M: The maximum likelihood approach to reconstructing ancestral character states of disccrete characters on phylogenies. Syst Biol. 1999, 48: 612-622. 10.1080/106351599260184.

    Google Scholar 

  112. Lindgren AR: Molecular inference of phylogenetic relationships among Decapo diformes (Mollusca: Cephalopoda) with special focus on the squid Order Oegopsida. Mol. Phylogenet. Evol. 2010, 56: 1-14. 10.1016/j.ympev.2010.03.029.

    Google Scholar 

  113. Mindell DP, Honeycutt RL: Ribosomal RNA in vertebrates: evolution and phylogenetic applications. Annu Rev Ecol Evol Syst. 1990, 21: 541-566. 10.1146/annurev.es.21.110190.002545.

    Google Scholar 

  114. Phillips MJ, Delsuc F, Penny D: Genome-scale phylogeny and the detection of systematic biases. Mol. Biol. Evolution. 2004, 21: 1455-1458. 10.1093/molbev/msh137.

    CAS  Google Scholar 

  115. Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.

    PubMed  CAS  Google Scholar 

  116. Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. 2002, Sinauer Associates, Sunderland, Massachusetts

    Google Scholar 

  117. Xia X, Xie Z: DAMBE: Data analysis in molecular biology and evolution. J Hered. 2001, 92: 371-373. 10.1093/jhered/92.4.371.

    PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We are very grateful to two anonymous reviewers, and to Markus Pfenninger, who provided a very rapid, helpful, and comprehensive review of this manuscript. Thanks to the Santa Barbara Museum of Natural History for granting the authors access to its vast cephalopod collection. The authors are grateful for the assistance of Patrizia Jereb, who supplied a substantial amount of information for our habitat and depth-coding matrix. Thanks also to Dr. Michael Vecchione for providing editorial comments to an earlier draft of this manuscript. Support for this research was provided by DEB-0906033 to A.R.L. and DEB-0643840 to T.H.O.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Annie R Lindgren.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ARL, MSP, and THO conceived of the study, participated in its design and coordination, and helped to draft the manuscript. ARL and MSP generated the molecular database from GenBank and vetted it for possible contaminants. MSP and THO wrote the scripts used in the data analysis pipeline. ARL, MSP and FGH generated the character matrix; MSP tested character correlations and assessed ancestral reconstruction. FGH provided access to museum collections, oversaw habitat/character coding, and helped to edit the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12862_2012_2069_MOESM1_ESM.doc

Additional file 1: Supplementary Material: contents [112-117]. (DOC 3 MB)

12862_2012_2069_MOESM2_ESM.pdf

Additional file 2: Habitat and morphological character data used in analyses. Abbreviations for Depth range coded from the literature are: “ukn” = data is not available for a particular species, “?” = a depth could not be defined as no species name was designated, “D” = depth in day, “N” = depth at night, “juv” = depth for juveniles, “OM” = ontogenetic migratory, “VM” = vertical migrator. Habitat states are 0 = pelagic, 1 = benthic/demersal. Character states are (0 = absent, 1 = present) for ANG, Branchial canal, Autogenic photophore, Bacteriogenic Photophore, and Right Oviduct. States for Cornea are: 0 = absent, 1 = 1-part cornea present, 2 = 2-part cornea present. (PDF 150 KB)

12862_2012_2069_MOESM3_ESM.pdf

Additional file 3: Appendix 2.GenBank sequence identifiers (GI numbers) for loci used in analyses listed alphapbetically. Changes to the species name used in this manuscript are listed in bold. Checkmarks denote taxa used in primary analysis (188-taxa), shown in Figure1. (PDF 51 KB)

Additional file 4: Appendix 3.GI numbers for all contaminant sequences. (PDF 185 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lindgren, A.R., Pankey, M.S., Hochberg, F.G. et al. A multi-gene phylogeny of Cephalopoda supports convergent morphological evolution in association with multiple habitat shifts in the marine environment. BMC Evol Biol 12, 129 (2012). https://doi.org/10.1186/1471-2148-12-129

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-12-129

Keywords