Isolation, characterization, and expression analysis of three actin genes in the New Zealand black-footed abalone, Haliotis iris

Three actin genes -- H. irisA1, H. irisA2, and H. irisA3 -- from the mollusc Haliotis iris (New Zealand black-footed abalone) were isolated by polymerase chain reaction (PCR). The genes were similar to molluscan (84.1% to 94.9%) and vertebrate (84.5% to 86.6%) actins. The sequence similarity between the genes ranged from 88.5% to 93.2%. The greatest disparity, 32.3%, was found over a 99-nt region located at nt 808-906 of H. irisA1, corresponding to amino acids 212-244 of the three actins. The H. iris actin gene family contains at least eight members. Reverse transcription (RT)-PCR analysis of the three genes showed H. irisA1 and H. irisA2 were expressed at low levels in fertilized eggs and blastula stages and at high levels in trochophore and veliger larvae. H. irisA3 was detected in fertilized eggs; it was not detected in the blastula stages and at high levels in the trochophore and veliger larvae. The structure and expression of the three actin genes are discussed.

Rapid evolution of reproductive proteins in abalone and Drosophila

Observations from different taxa, including plants, protozoa, insects and mammals, indicate that proteins involved in reproduction evolve rapidly. Several models of adaptive evolution have been proposed to explain this phenomenon, such as sexual conflict, sexual selection, self versus non-self recognition and pathogen resistance. Here we discuss the potential role of sexual conflict in the rapid evolution of reproductive genes in two different animal systems, abalone (Haliotis) and Drosophila. In abalone, we reveal how specific interacting sperm-egg proteins were identified and discuss this identification in the light of models for rapid protein evolution and speciation. For Drosophila, we describe the genomic approaches taken to identify male accessory gland proteins and female reproductive tract proteins. Patterns of protein evolution from both abalone and Drosophila support the predicted patterns of rapid protein evolution driven by sexual conflict. We stress however that other selective pressures may contribute to the rapid evolution that is observed. We conclude that the key to distinguishing between sexual conflict and other mechanisms of protein evolution will be an integration of genetic, experimental and theoretical data.

Genes expressed in a turrid venom duct: divergence and similarity to conotoxins

The toxoglossate mollusks are a large group of venomous animals (>10,000 species) conventionally divided into three groups, the cone snails, the auger snails, and the turrid snails; turrids account for >90% of the biodiversity of toxoglossans. Only the venoms of cone snails have been intensively investigated, with little work focused on turrids. We report the first broad characterization of genes expressed in venom ducts of any turrid species. Twenty-three different cDNA clones encoding putative toxins were characterized from the venom duct of the turrine species Lophiotoma olangoensis Olivera 2002 and belong to 16 different gene families. Of the 16 different Lophiotoma olangoensis gene families that encode putative toxins, for only 1 was there clear evidence of sequence similarity with any conotoxin gene family. The I-like gene family of Lophiotoma olangoensis was found to be related to the K channel-targeted I(2) conotoxin superfamily. Most putative Lophiotoma toxins are cysteine-rich polypeptides, with a significant fraction much larger (>80 amino acids) than the toxins from cone snails. A small number were not cysteine-rich but had hydrophobic amino acid clusters interspersed with arginine residues. This is only 1 of >10,000 different turrid venoms that needs to be characterized. From this study, a common origin with Conus for one family of putative turrid toxins is indicated.

Molecular Evolution and Functionally Important Structures of Molluscan Dermatopontin: Implications for the Origins of Molluscan Shell Matrix Proteins

A major shell matrix protein originally obtained from a freshwater snail is a molluscan homologue of Dermatopontins, a group of Metazoan proteins also called TRAMP (tyrosine-rich acidic matrix protein). We sequenced and identified 14 molluscan homologues of Dermatopontin from eight snail species belonging to the order Basommatophora and Stylommatophora. The bassommatophoran Dermatopontins fell into three types, one is suggested to be a shell matrix protein and the others are proteins having more general functions based on gene expression analyses. N-glycosylation is inferred to be important for the function involved in shell calcification, because potential N-glycosylation sites were found exclusively in the Dermatopontins considered as shell matrix proteins. The stylommatophoran Dermatopontins fell into two types, also suggested to comprise a shell matrix protein and a protein having a more general function. Phylogenetic analyses using maximum likelihood and Bayesian methods revealed that gene duplication events occurred independently in both basommatophoran and stylommatophoran lineages. These results suggest that the dermatopontin genes were co-opted for molluscan calcification at least twice independently after the divergence of basommatophoran and stylommatophoran lineages, or more recently than we have expected.

Evolutionary responses of natives to introduced species: what do introductions tell us about natural communities?

Biological invasions dramatically affect the distribution, abundance and reproduction of many native species. Because of these ecological effects, exotic species can also influence the evolution of natives exposed to novel interactions with invaders. Evolutionary changes in natives in response to selection from exotics are usually overlooked, yet common responses include altered anti-predator defenses, changes in the spectrum of resources and habitats used, and other adaptations that allow native populations to persist in invaded areas. Whether a native population is capable of responding evolutionarily to selection from invaders will depend on the demographic impact of the invader, the genetic architecture and genetic variability of the native population and potentially the history of previous invasions. In some cases, natives will fail to evolve or otherwise adapt, and local or global extinction will result. In other cases, adaptive change in natives may diminish impacts of invaders and potentially promote coexistence between invaders and natives. Here, we review the evidence for evolutionary responses of native species to novel community members. We also discuss how the effects of introduced species may differ from those caused by natural range expansions of native species. Notably, introduced species may come from remote biotas with no previous evolutionary history with the native community. In addition, the rate of addition of introduced species into communities is much greater than all but the most extreme cases of historical biotic exchange. Understanding the evolutionary component of exotic/native species interactions is critical to recognizing the long-term impacts of biological invasions, and to understanding the role of evolutionary processes in the assembly and dynamics of natural communities.

Rangewide phylogeography of a terrestrial slug in Europe: evidence for Alpine refugia and rapid colonization after the Pleistocene glaciations

Intraspecific phylogeographical patterns largely depend on the life history traits of a species. Especially species with a high degree of cold tolerance, limited requirements towards habitat preferences, and relatively low active dispersal capacities may have responded in a different way to the Pleistocene climatological fluctuations than the majority of taxa studied so far. To evaluate this possibility, we studied Arion fuscus (Muller, 1774), a common and widespread European terrestrial slug, from 88 locations (N = 964). Sequence variation was assessed for fragments of the mitochondrial 16S rDNA and COI genes by means of single-strand conformation polymorphisms (SSCP) and subsequent DNA sequencing. Additionally, eight allozyme loci were scored in 843 individuals. Phylogenetic analysis revealed the presence of two major evolutionary lineages, one in the Balkan region and another in the Alps and the rest of Europe. The sequence divergence between the two lineages was limited (3.3%), but gene flow between the regions was absent, suggesting that the two regions have been isolated since the late Pliocene or early Pleistocene. Allozyme differentiation among geographical regions and mitochondrial DNA (mtDNA) lineages was low. The geographical patterns observed in our data showed that (i) haplotype and nucleotide diversities are very low in northern Europe, suggesting that single haplotypes rapidly colonized large areas; (ii) recently expanded haplotype clades have restricted distribution ranges, suggesting that current gene flow is low; and (iii) genetic diversity in the Alps is much higher than in other regions and estimated past gene flow from the Eastern Alps to other regions was high, suggesting that this was a refugial zone during the Pleistocene. This full-range phylogeography suggests the existence of an alternative refugial zone, situated north of the refugial areas currently recognized in most other taxa.

Endopolyploidy as a morphogenetic factor of development

This paper summarizes the works published by author and his co-workers in the Russian journal Tsitologiya concerning endopolyploidy in mollusks and appraises this phenomenon in general. Both ontogenetic and phylogenetic aspects of endopolyploidy have been studied. In the snail Succinea lauta, a complex examination of endomitosis has been performed. A regular replacement of the normal (complete) proliferative mitosis by abnormal (incomplete) restitutional mitosis, and then by Geitler's classic endomitosis has been demonstrated. We examined 29 bivalve and 82 gastropod species for the presence of polyploid cells in glandular tissues and ganglia. In the bivalve species, ordinary diploid cells form various tissues, while in the gastropods, the role of polyploidy in tissue development appears to increase in phylogenesis. The rise of endopolyploidy and cell giantism in histogeneses of a variety of animal and plant species is widely known. It is believed to be a regular event in the evolution of certain groups. To give a universal interpretation of endopolyploidy, we proposed that a single polyploid cell be better considered as an endoclone. In this case, evolutionary transformation of diploid cell clones into polyploid endoclones may be viewed as Dogel's oligomerization applied to cell-tissue level. From this viewpoint, major properties of an oligomerized system (intensification of function, functional efficiency (ergonomy), increased genomes reliability, simplification of the intra- and supersystem regulations, and acceleration of development) can be considered as principal peculiarities of polyploid growth strategy. The above peculiarities allow one to consider endopolyploidy as an additional means of integrative onto(histo)genetic regulations and correlations and as an important evolutionary factor (coordinations) acting through natural selection. Thus, in general, endopolyploidy is an adaptive morphogenetic factor, but its concrete role may differ in different tissues and organisms depending on cell specialization and histogenetic particularities.

Hypotaurocyamine kinase evolved from a gene for arginine kinase

Hypotaurocyamine kinase (HTK) is a member of the highly conserved family of phosphagen kinases that includes creatine kinase (CK) and arginine kinase (AK). HTK is found only in sipunculid worms, and it shows activities for both the substrates hypotaurocyamine and taurocyamine. Determining how HTK evolved in sipunculids is particularly insightful because all sipunculid-allied animals have AK and only some sipunculids have HTK. We determined the cDNA sequence of HTK from the sipunculid worm Siphonosoma cumanense for the first time, cloned it in pMAL plasmid and expressed it in E. coli as a fusion protein with maltose-binding protein. The cDNAderived amino acid sequence of Siphonosoma HTK showed high amino acid identity with molluscan AKs. Nevertheless, the recombinant enzyme of Siphonosoma HTK showed no activity for the substrate arginine, but showed activity for taurocyamine. Comparison of the amino acid sequences of HTK and AK indicated that the amino acid residues necessary for the binding of the substrate arginine in AK have been completely lost in Siphonosoma HTK sequence. The phylogenetic analysis indicated that the HTK amino acid sequence was placed just outside the molluscan AK cluster, which formed a sister group with the arthropod and nematode AKs. These results suggest that Siphonosoma HTK evolved from a gene for molluscan AK. Moreover, to confirm this assertion, we determined by PCR that the gene for Siphonosoma HTK has a 5-exon/4-intron structure, which is homologous with that of the molluscan AK genes. Further, the positions of splice junctions were conserved exactly between the two genes. Thus, we conclude that Siphonosoma HTK has evolved from a primordial gene for molluscan AK.

Flow cytometric DNA content analysis of neoplastic cells in haemolymph of the cockle Cerastoderma edule

Epizootiological outbreaks of disseminated neoplasia (DN) have been reported in association with mass mortalities in various bivalve species including the cockle Cerastoderma edule. A flow cytometric (FCM) procedure to study DNA content was successfully adapted and tested in haemolymph cells (haemocytes and neoplastic cells) of the cockle. The FCM results were similar to those obtained by histological analysis (DN diagnosis and haemolymph cell features). FCM analysis revealed differences in DNA content among normal haemocytes (diploid) and neoplastic cells. Four types of cells with abnormal DNA content were found in the haemolymph of affected animals: hypodiploid, hyperdiploid, triploid-sesploid and pentaploid. Our results suggest that the flow cytometric DNA content analysis can be applied to identify neoplastic cell types and to study the association between different cell types and the DN progression or remission in this edible and commercially important bivalve species.

Molecular cloning of a molluscan gonadotropin-releasing hormone receptor orthologue specifically expressed in the gonad

Despite their economic importance, only very little information is available regarding (neuro)endocrine mechanisms of reproduction in bivalve molluscs. To gain insights into the molecular control of gonadic development of these animals, G protein-coupled receptors (GPCR) expressed in the gonad of the pacific oyster Crassostrea gigas were investigated. One such receptor was cloned by RT-PCR using oligonucleotide primers derived from consensus sequences of various vertebrate (neuro)peptide receptors. This receptor named Cg-GnRH-related receptor (Cg GnRH-R) exhibits a high degree of amino acid sequence identity with both vertebrate GnRH receptors and insect AKH receptors. Quantitative RT-PCR shows a specific expression of Cg-GnRH-R in both male and female gonads during the reproductive cycle. This demonstrates for the first time the plausible involvement of a GnRH receptor orthologue in the control of reproduction in a protostomian invertebrate.

Sequence divergence and conservation in organizationally distinct subfamilies of Donax trunculus satellite DNA

Characterization of a low-copy number DTF1 satellite DNA detected in the bivalve mollusk Donax trunculus revealed extensive grouping of monomer sequence variants into subfamilies identified by distinctive combinations of diagnostic nucleotides. It can be anticipated that a large number of subfamilies exists in the genome. In addition to the tandem organization of 169 bp long monomers, at least one subfamily was created through amplification of adjacent repeats in a higher order register. This complex satellite unit consists of two distinctive monomer variants that differ both in specific nucleotide changes and in a deleted segment partially substituted with a short unrelated sequence element. Most of the nucleotide substitutions differing between subfamilies are highly homogenized within a corresponding group of monomer variants, and intra-subfamily variability in general is low. Nucleotide diversity analysis of all sequenced variants of DTF1 satellite revealed the presence of two conserved segments, while the rest of the monomer sequence shows uniform and considerably higher level of variability. The persistence of conserved segments stands in contrast to the sequence and organizational divergence of monomer variant groups, and may indicate constraints in the evolution of DTF1 satellite repeats.

Developmental expression of two Haliotis asinina hemocyanin isoforms

Hemocyanins are large copper-containing respiratory proteins that play a role in oxygen transport in many molluscs. In some species only one hemocyanin isoform is present while in others two are expressed. The physiological relevance of these isoforms is unclear and the developmental and tissue-specific expression of hemocyanin genes is largely unknown. Here we show that two hemocyanin genes in the gastropod Haliotis asinina, which encode H. asinina hemocyanin (HaH1) and HaH2 isoforms, are developmentally expressed. These genes initially are expressed in a small number of mesenchyme cells at trochophore and pre-torsional veliger stages, with HaH1 expression slightly preceding HaH2. These cells largely are localized to the visceral mass, although a small number of cells are present in head and foot regions. Following metamorphosis the isoforms show overlapping as well as isoform-specific expression profiles, suggesting some degree of isoform-specific function.

Low genetic diversity in a marine nature reserve: re-evaluating diversity criteria in reserve design

Little consideration has been given to the genetic composition of populations associated with marine reserves, as reserve designation is generally to protect specific species, communities or habitats. Nevertheless, it is important to conserve genetic diversity since it provides the raw material for the maintenance of species diversity over longer, evolutionary time-scales and may also confer the basis for adaptation to environmental change. Many current marine reserves are small in size and isolated to some degree (e.g. sea loughs and offshore islands). While such features enable easier management, they may have important implications for the genetic structure of protected populations, the ability of populations to recover from local catastrophes and the potential for marine reserves to act as sources of propagules for surrounding areas. Here, we present a case study demonstrating genetic differentiation, isolation, inbreeding and reduced genetic diversity in populations of the dogwhelk Nucella lapillus in Lough Hyne Marine Nature Reserve (an isolated sea lough in southern Ireland), compared with populations on the local adjacent open coast and populations in England, Wales and France. Our study demonstrates that this sea lough is isolated from open coast populations, and highlights that there may be long-term genetic consequences of selecting reserves on the basis of isolation and ease of protection.

Evidence for biased use of sperm sources in wild female giant cuttlefish (Sepia apama)

In species where females store sperm from their mates prior to fertilization, sperm competition is particularly probable. Female Sepia apama are polyandrous and have access to sperm from packages (spermatangia) deposited by males onto their buccal area during mating and to sperm stored in internal sperm-storage organs (receptacles) located below the beak. Here, we describe the structure of the sperm stores in the female's buccal area, use microsatellite DNA analyses to determine the genetic diversity of stored sperm and combine these data with offspring genotypes to determine the storage location of paternal sperm. The number of male genotypes represented in the sperm receptacles was significantly lower than that found among the spermatangia. Estimation of the volumes of sperm contained in the receptacles and the spermatangia were statistically comparable; however, paternal sperm were more likely to have come from spermatangia than from the sperm receptacles. These results confirm a genetic polyandrous mating system in this species and suggest that fertilization pattern with respect to the sperm stores used is not random.

Molecular cloning and functional characterization of the Aplysia FMRFamide-gated Na+ channel

FMRFamide-gated Na+ channel (FaNaC) is the only known peptide-gated ion channel, which belongs to the epithelial Na+ channel/degenerin (ENaC/DEG) family. We have cloned a putative FaNaC from the Aplysia kurodai CNS library using PCR, and examined its characteristics in Xenopus oocytes. A. kurodai FaNaC (AkFaNaC) comprised with 653 amino acids, and the sequence predicts two putative membrane domains and a large extracellular domain as in other members of the ENaC/DEG family. In oocytes expressing AkFaNaC, FMRFamide evoked amiloride-sensitive Na+ current. Different from the known FaNaCs (Helix and Helisoma FaNaCs), AkFaNaC was blocked by external Ca2+ but not by Mg2+. Also, desensitization of the current was enhanced by Mg2+ but not by Ca2+. The FMRFamide-gated current was depressed in both low and high pH. These results indicate that AkFaNaC is an FaNaC of Aplysia, and that the channel has Aplysia specific functional domains.

Mutational and structural studies of the diisopropylfluorophosphatase from Loligo vulgaris shed new light on the catalytic mechanism of the enzyme

The active site, the substrate binding site, and the metal binding sites of the diisopropylfluorophosphatase (DFPase) from Loligo vulgaris have been modified by means of site-directed mutagenesis to improve our understanding of the reaction mechanism. Enzymatic characterization of mutants located in the major groove of the substrate binding pocket indicates that large hydrophobic side chains at these positions are favorable for substrate turnover. Moreover, the active site residue His287 proved to be beneficial, but not essential, for DFP hydrolysis. In most cases, hydrophobic side chains at position 287 led to significant catalytic activities although reduced relative to the wild-type enzyme. With respect to the Ca-1 binding site, where catalysis occurs, various mutants indicated that the net charge at this calcium-binding site as well as the relative positions of the charged calcium ligands is crucial for catalytic activity. The importance of the electrostatic potential at the active site was furthermore revealed by various mutations of residues lining the interior of the central water-filled tunnel, which traverses the entire protein structure. In this respect, the structural features of residue His181, which is located at the opposite end of the DFPase tunnel relative to the active site, were characterized extensively. It was concluded that a tunnel-spanning hydrogen bond network, which includes a large number of apparently slow exchanging water molecules, relays any modifications in the electrostatics of the system to the active site, thus affecting the catalytic reactivity of the enzyme.

The HOX Gene Cluster in the Bivalve Mollusc Mytilus galloprovincialis

The clustered Hox genes play a central role in the regulation of development in bilaterian animals. In this study, we analyzed the homeobox-containing genes in a bivalve mollusc, the mussel Mytilus galloprovincialis, an unsegmented spiralian lophotrochozoan. We isolated and characterized four Hox cluster genes using the polymerase chain reaction with specific primers. Molecular alignments and phylogenetic analysis indicate that these mussel genes are homologs of the anterior group (pb ortholog), paralog group 3, and central group (PG4/Dfd and PG5/Scr) genes. The putative homeodomain sequences were designated Mgox1, Mgox2, Mgox3, and Mgox4.

Precursors of Novel Gla-Containing Conotoxins Contain a Carboxy-Terminal Recognition Site That Directs γ-Carboxylation,

Vitamin K-dependent gamma-glutamyl carboxylase catalyzes the conversion of glutamyl residues to gamma-carboxyglutamate. Its substrates include vertebrate proteins involved in blood coagulation, bone mineralization, and signal transduction and invertebrate ion channel blockers known as conotoxins. Substrate recognition involves a recognition element, the gamma-carboxylation recognition site, typically located within a cleavable propeptide preceding the targeted glutamyl residues. We have purified two novel gamma-carboxyglutamate-containing conotoxins, Gla-TxX and Gla-TxXI, from the venom of Conus textile. Their cDNA-deduced precursors have a signal peptide but no apparent propeptide. Instead, they contain a C-terminal extension that directs gamma-carboxylation but is not found on the mature conotoxin. A synthetic 13-residue "postpeptide" from the Gla-TxXI precursor reduced the K(m) for the reaction of the Conus gamma-carboxylase with peptide substrates, including FLEEL and conantokin-G, by up to 440-fold, regardless of whether it was positioned at the N- or C-terminal end of the mature toxin. Comparison of the postpeptides to propeptides from other conotoxins suggested some common elements, and amino acid substitutions of these residues perturbed gamma-carboxylation of the Gla-TxXI peptide. The demonstration of a functional and transferable C-terminal postpeptide in these conotoxins indicates the presence of the gamma-carboxylation recognition site within the postpeptide and defines a novel precursor structure for vitamin K-dependent polypeptides. It also provides the first formal evidence to prove that gamma-carboxylation occurs as a post-translational rather than a cotranslational process.

Cerastodermaglaucum5S ribosomal DNA: characterization of the repeat unit, divergence with respect toCerastoderma edule, and PCR–RFLPs for the identification of both cockles

The 5S rDNA repeat unit of the cockle Cerastoderma glaucum from the Mediterranean and Baltic coasts was PCR amplified and sequenced. The length of the units was 539–568 bp, of which 120 bp were assigned to the 5S rRNA gene and 419–448 bp to the spacer region, and the G/C content was 46%–49%, 54%, and 44%–47%, respectively. Two types of units (A and B), differing in the spacer, were distinguished based on the percentage of differences and clustering in phylogenetic trees. A PCR assay with specific primers for each unit type indicated that the occurrence of both units is not restricted to the sequenced individuals. The 5S rDNA units of C. glaucum were compared with new and previously reported sequences of Cerastoderma edule. The degree of variation observed in C. edule was lower than that in C. glaucum and evidence for the existence of units A and B in C. edule was not found. The two cockles have the same coding region but displayed numerous fixed differences in the spacer region and group separately in the phylogenetic trees. Digestion of the 5S rDNA PCR product with the restriction enzymes HaeIII and EcoRV revealed two RFLPs useful for cockle identification.Key words: Cerastoderma, cockle identification, 5S ribosomal DNA, nontranscribed spacer variation, PCR-RFLP.

Oxidative folding of conotoxins sharing an identical disulfide bridging framework

Conotoxins are short, disulfide-rich peptide neurotoxins produced in the venom of predatory marine cone snails. It is generally accepted that an estimated 100,000 unique conotoxins fall into only a handful of structural groups, based on their disulfide bridging frameworks. This unique molecular diversity poses a protein folding problem of relationships between hypervariability of amino acid sequences and mechanism(s) of oxidative folding. In this study, we present a comparative analysis of the folding properties of four conotoxins sharing an identical pattern of cysteine residues forming three disulfide bridges, but otherwise differing significantly in their primary amino acid sequence. Oxidative folding properties of M-superfamily conotoxins GIIIA, PIIIA, SmIIIA and RIIIK varied with respect to kinetics and thermodynamics. Based on rates for establishing the steady-state distribution of the folding species, two distinct folding mechanisms could be distinguished: first, rapid-collapse folding characterized by very fast, but low-yield accumulation of the correctly folded form; and second, slow-rearrangement folding resulting in higher accumulation of the properly folded form via the reshuffling of disulfide bonds within folding intermediates. Effects of changing the folding conditions indicated that the rapid-collapse and the slow-rearrangement mechanisms were mainly determined by either repulsive electrostatic or productive noncovalent interactions, respectively. The differences in folding kinetics for these two mechanisms were minimized in the presence of protein disulfide isomerase. Taken together, folding properties of conotoxins from the M-superfamily presented in this work and from the O-superfamily published previously suggest that conotoxin sequence diversity is also reflected in their folding properties, and that sequence information rather than a cysteine pattern determines the in vitro folding mechanisms of conotoxins.

Isolation of Hox and ParaHox genes in the bivalve Pecten maximus

The Hox cluster genes encode a set of transcription factors that have been shown to control spatial patterning mechanisms in bilaterian organism development. The ParaHox cluster is the evolutionary sister of the Hox cluster. The two are believed to descend from an ancestral ProtoHox cluster of four genes from which the three ParaHox genes (Gsx, Xlox, and Cdx) and the four Hox gene classes are believed to have originated. Although molluscs are among the most successful lophotrochozoan groups, very little work has been devoted to the characteristics of their homeotic genes. Using polymerase chain reaction-based approaches, we isolated 13 different Pecten maximus (Bivalvia: Pteriomorphia) sequences corresponding to all the genes of the four Hox cluster classes and to genes Xlox and Cdx of the ParaHox cluster. Comparison of results with those obtained in other lophotrochozoans seems to confirm the considerable homogeneity of the Hox and ParaHox genes in these taxa both as regards the presence of nearly all the genes of the two clusters and the marked sequence resemblance among orthologous genes.

Phylogenetic relationships among major species of japanese coleoid cephalopods (Mollusca: Cephalopoda) using three mitochondrial DNA sequences

Phylogenetic relationships among 36 species of major coleoid cephalopods from Japanese waters were studied using partial sequences of three mitochondrial genes, 16S rDNA, 12S rDNA, and cytochrome c oxidase subunit I gene. Octopoda and Decapoda were monophylic groups. Within Sepioidea, Sepiadariidae and Sepiolidae were not closely related to Sepiidae, but rather related to Teuthoidea. Sepiidae with a distinct calcareous shell formed a single cluster. Myopsida was closely related to Oegopsida. Within Octopoda, Opisthoteuthis depressa and Argonauta argo diverged earlier than Octopodiidae. The common octopuses in Japanese waters were separated into three clusters. The first cluster occupied a basal position, and includes large-sized octopuses, such as Enteroctopus dofleini and Octopus (Paroctopus) conispadiceus from the continental shelf and upper slope. The second cluster consisted of long-armed octopuses, such as O. ornatus, O. minor, and O. sasakii. The third cluster contained small- to medium-sized octopus, such as Amphioctopus fangsiao, A. areolatus, O. cyaneus, and O. vulgaris, in which several species possess ocelli on the web. The second cluster formed the sister group to the third cluster.