Molecular Evolution of the Metazoan Extracellular Matrix: Cloning and Expression of Structural Proteins from the Demosponges Suberites domuncula and Geodia cydonium

One crucial event during evolution to multicellularity was the development of either direct cell–cell contact or indirect interaction via extracellular matrix (ECM) molecules. The identification of those polypeptides provides conclusive data on the phylogenetic relationship of metazoan phyla and helps us to understand the position of the Metazoa among the other kingdoms. Recently it became evident that the ECM of sponges is amazingly complex; it is composed of fibrous molecules, e.g., collagen, and their corresponding receptors, which are highly similar to those existing in other metazoan phyla. While these data already support the view of monophyly of Metazoa, additional studies are required to understand whether these molecules, which are similar in their primary sequence, also have the same function throughout the metazoan kingdom. In the present study we identified the ligand for one of the autopomorphic characters of Metazoa, the single-transmembrane receptor protein with the receptor tyrosine kinase (RTK) from G. cydonium, as an example: the putative mucus-like protein from G. cydonium. This protein was upregulated during autograft fusion in the homologous system with kinetics similar to those of the RTK. Additionally, a cDNA was isolated from S. domuncula whose deduced polypeptide displays a high sequence similarity to dermatopontin, an ECM molecule found exclusively in Metazoa. Furthermore, it is documented that expression of the fibrous ECM molecule collagen is regulated by the characteristic metazoan morphogens myotrophin and endothelial monocyte-activating polypeptide. These data indicate that the ECM of sponges is not an unstructured ground substance but provides the basis for integrated cell communication. Received: 26 October 2000 / Accepted: 1 February 2001     

Molecular evolution of the Metazoan protein kinase C multigene family

Protein kinases C (PKCs) comprise closely related Ser/Thr kinases, ubiquitously present in animal tissues; they respond to second messengers, e.g., Ca²⁺ and/or diacylglycerol, to express their activities. Two PKCs have been sequenced from Geodia cydonium, a member of the lowest multicellular animals, the sponges (Porifera). One sponge G. cydonium PKC, GCPKC1, belongs to the ``novel' (Ca<sup>2+</sup>-independent) PKC (nPKC) subfamily while the second one, GCPKC2, has the hallmarks of the ``conventional' (Ca²⁺-dependent) PKC (cPKC) subfamily. The alignment of the Ser/Thr catalytic kinase domains, of the predicted aa sequences for these cDNAs with respective segments from previously reported sequences, revealed highest homology to PKCs from animals but also distant relationships to Ser/Thr kinases from protozoa, plants, and bacteria. However, a comparison of the complete structures of the sponge PKCs, which are—already—identical to those of nPKCs and cPKCs from higher metazoa, with the structures of protozoan, plant, and bacterial Ser/Thr kinases indicates that the metazoan PKCs have to be distinguished from the nonmetazoan enzymes. These data indicate that metazoan PKCs have a universal common ancestor which they share with the nonmetazoan Ser/Thr kinases with respect to the kinase domain, but they differ from them in overall structural composition. Received: 10 January 1996 / Accepted: 12 March 1996     

Structural Evidence for the Evolution of Pyrogenic Toxin Superantigens

Cells from metazoan organisms are eliminated in a variety of physiological and pathophysiological processes by apoptosis. In this report, we describe the cloning and characterization of molecules from the marine sponges Geodia cydonium and Suberites domuncula, whose domains show a high similarity to those that are found in molecules of the vertebrate Bcl-2 superfamily and of the death receptors. The Bcl-2 proteins contain up to four Bcl-2 homology regions (BH). Two Bcl-2-related molecules have been identified from sponges that are provided with two of those regions, BH1 and BH2, and are termed Bcl-2 homology proteins (BHP). The G. cydonium molecule, BHP1_GC, has a putative size of 28,164, while the related sequence from S. domuncula, BHP1_SD, has a M _r of 24,187. Phylogenetic analyses of the entire two sponge BHPs revealed a high similarity to members of the mammalian Bcl-2 superfamilies and to the Caenorhabditis elegans Ced-9. When the two domains, BH1 and BH2, are analyzed separately, again the highest similarity was found to the members of the Bcl-2 superfamily, but a clearly lower relationship to the C. elegans BH1 and BH2 domains in Ced-9. In unrooted phylogenetic trees the sponge BH1 and BH2 are grouped among the mammalian sequences and are only distantly related to the C. elegans BH domains. The analysis of the gene structure of the G. cydonium BHP showed that the single intron present is located within the BH2 domain at the same position as in C. elegans and rat Bcl-x_L. In addition, a sponge molecule comprising two death domains has been characterized from G. cydonium. The two death domains of the potential proapoptotic molecule GC_DD2, M _r 24,970, share a high similarity with the Fas-FADD/MORT1 domains. A death domain-containing molecule has not been identified in the C. elegans genome. The phylogenetic analysis revealed that the sponge domain originated from an ankyrin building block from which the mammalian Fas-FADD/MORT1 evolved. It is suggested that the apoptotic pathways that involve members of the Bcl-2 superfamily and of the death receptors are already present in the lowest metazoan phylum, the Porifera. Received: 27 July 1999 / Accepted: 28 December 1999     

Extensive Gene Duplication in the Early Evolution of Animals Before the Parazoan–Eumetazoan Split Demonstrated by G Proteins and Protein Tyrosine Kinases from Sponge and Hydra

To know whether genes involved in cell–cell communication typical of multicellular animals dramatically increased in concert with the Cambrian explosion, the rapid evolutionary burst in the major groups of animals, and whether these genes exist in the sponge lacking cell cohesiveness and coordination typical of eumetazoans, we have carried out cloning of the G-protein α subunit (Gα) and the protein tyrosine kinase (PTK) cDNAs from Ephydatia fluviatilis (freshwater sponge) and Hydra magnipapillata strain 105 (hydra). We obtained 13 Gα and 20 PTK cDNAs. Generally animal gene families diverged first by gene duplication (subtype duplication) that gave rise to diverse subtypes with different primary functions, followed by further gene duplication in the same subtype (isoform duplication) that gave rise to isoform genes with virtually identical function. Phylogenetic trees of Gα and PTK families including cDNAs from sponge and hydra revealed that most of the present-day subtypes had been established in the very early evolution of animals before the parazoan–eumetazoan split, the earliest branching among the extant animal phyla, by extensive subtype duplication: for PTK and Gα families, 23 and 9 subtype duplications were observed in the early stage before the parazoan–eumetazoan split, respectively, and after that split, only 2 and 1 subtype duplications were found, respectively. After the separation from arthropods, vertebrates underwent frequent isoform duplications before the fish–tetrapod split. Furthermore, rapid amino acid changes appear to have occurred in concert with the extensive subtype duplication and isoform duplication. Thus the pattern of gene diversification during animal evolution might be characterized by bursts of gene duplication interrupted by considerably long periods of silence, instead of proceeding gradually, and there might be no direct link between the Cambrian explosion and the extensive gene duplication that generated diverse functions (subtypes) of these families. Received: 4 November 1998 / Accepted: 17 November 1998     

Origin and Evolution of Paralogous rRNA Gene Clusters Within the Flatworm Family Dugesiidae (Platyhelminthes, Tricladida)

Analysis of the 18S rDNA sequences of five species of the family Dugesiidae (phylum Platyhelminthes, suborder Tricladida, infraorder Paludicola) and eight species belonging to families Dendrocoelidae and Planaridae and to the infraorder Maricola showed that members of the family Dugesiidae have two types of 18S rDNA genes, while the rest of the species have only one. The duplication event also affected the ITS-1, 5.8S, ITS-2 region and probably the 28S gene. The mean divergence value between the type I and the type II sequences is 9% and type II 18S rDNA genes are evolving 2.3 times more rapidly than type I. The evolutionary rates of type I and type II genes were calibrated from biogeographical data, and an approximate date for the duplication event of 80–120 million years ago was calculated. The type II gene was shown, by RT-PCR, to be transcribed in adult individuals of Schmidtea polychroa, though at very low levels. This result, together with the fact that most of the functionally important positions for small-subunit rRNA in prokaryotes have been conserved, indicates that the type II gene is probably functional. Received: 24 March 1998 / Accepted: 17 March 1999     

Multiple Protein Tyrosine Phosphatases in Sponges and Explosive Gene Duplication in the Early Evolution of Animals Before the Parazoan–Eumetazoan Split

Protein tyrosine phosphatases (PTPs) regulate various physiological events in animal cells. They comprise a diverse family which are classified into two categories, receptor type and nonreceptor type. From the domain organization and phylogenetic tree, we have classified known PTPs into 17 subtypes (9 receptor-type and 8 nonreceptor-type PTPs) which are characterized by different organization of functional domain and independent cluster in tree. The receptor type PTPs are thought to be implicated in cell–cell adhesion by association of cell adhesion molecules. Since sponges are the most primitive multicellular animals and are thought to be lacking cell cohesiveness and coordination typical of eumetazoans, cloning and sequencing of PTP cDNAs of Ephydatia fluviatilis (freshwater sponge) have been conducted by RT-PCR to determine whether or not sponges have PTP genes in their genomes. We have isolated nine PTPs, of which five are possibly receptor type. A phylogenetic tree including the sponge PTPs revealed that most of the gene duplications that gave rise to the 17 subtypes had been completed in the very early evolution of animals before the parazoan–eumetazoan split, the earliest branching among extant animal phyla. The family tree also revealed the rapid evolutionary rate of PTP subtypes in the early stage of animal evolution. Received: 22 October 1998 / Accepted: 27 November 1998     

Gene Order Breakpoint Evidence in Animal Mitochondrial Phylogeny

Multiple genome rearrangement methodology facilitates the inference of animal phylogeny from gene orders on the mitochondrial genome. The breakpoint distance is preferable to other, highly correlated but computationally more difficult, genomic distances when applied to these data. A number of theories of metazoan evolution are compared to phylogenies reconstructed by ancestral genome optimization, using a minimal total breakpoints criterion. The notion of unambiguously reconstructed segments is introduced as a way of extracting the invariant aspects of multiple solutions for a given ancestral genome; this enables a detailed reconstruction of the evolution of non-tRNA mitochondrial gene order. Received: 15 July 1998 / Accepted: 5 March 1999     

Limitations of Metazoan 18S rRNA Sequence Data: Implications for Reconstructing a Phylogeny of the Animal Kingdom and Inferring the Reality of the Cambrian Explosion

We document the phylogenetic behavior of the 18S rRNA molecule in 67 taxa from 28 metazoan phyla and assess the effects of among-site rate variation on reconstructing phylogenies of the animal kingdom. This empirical assessment was undertaken to clarify further the limits of resolution of the 18S rRNA gene as a phylogenetic marker and to address the question of whether 18S rRNA phylogenies can be used as a source of evidence to infer the reality of a Cambrian explosion. A notable degree of among-site rate variation exists between different regions of the 18S rRNA molecule, as well as within all classes of secondary structure. There is a significant negative correlation between inferred number of nucleotide substitutions and phylogenetic information, as well as with the degree of substitutional saturation within the molecule. Base compositional differences both within and between taxa exist and, in certain lineages, may be associated with long branches and phylogenetic position. Importantly, excluding sites with different degrees of nucleotide substitution significantly influences the topology and degree of resolution of maximum-parsimony phylogenies as well as neighbor-joining phylogenies (corrected and uncorrected for among-site rate variation) reconstructed at the metazoan scale. Together, these data indicate that the 18S rRNA molecule is an unsuitable candidate for reconstructing the evolutionary history of all metazoan phyla, and that the polytomies, i.e., unresolved nodes within 18S rRNA phylogenies, cannot be used as a single or reliable source of evidence to support the hypothesis of a Cambrian explosion. Received: 9 December 1997 / Accepted: 23 March 1998     

Evolution of the Integrin α and β Protein Families

A phylogenetic analysis of vertebrate and invertebrate α integrins supported the hypothesis that two major families of vertebrate α integrins originated prior to the divergence of deuterostomes and protostomes. These two families include, respectively, the αPS1 and αPS2 integrins of Drosophila melanogaster, and each family has duplicated repeatedly in vertebrates but not in Drosophila. In contrast, a third family (including αPS3) has duplicated in Drosophila but is absent from vertebrates. Vertebrate αPS1 and αPS2 family members are found on human chromosomes 2, 12, and 17. Linkage of these family members may have been conserved since prior to the origin of vertebrates, and the two genes duplicated simultaneously. A phylogenetic analysis of β integrins did not clearly resolve whether vertebrate β integrin genes duplicated prior to the origin of vertebrates, although it suggested that at least the gene encoding vertebrate β4 may have done so. In general, the phylogeny of neither α nor β integrins showed a close correspondence with patterns of α–β heterodimer formation or other functional characteristics. One major exception to this trend involved αL, αM, αX, and αD, a monophyletic group of immune cell-expressed α integrins, which share a number of common functional characteristics and have evolved in coordinated fashion with their β integrin partners. Received: 22 June 2000 / Accepted: 11 September 2000     

Structure of the Transmembrane Cysteine Residues in Phospholamban

Phospholamban, a 52-residue membrane protein, associates to form a pentameric complex of five long α-helices traversing the sarcoplasmic reticulum membrane of cardiac muscle cells. The transmembrane domain of the protein is largely hydrophobic, with only three cysteine residues having polar side chains, yet it functions as a Ca²⁺-selective ion channel. In this report, infrared spectroscopy is used to probe the conformation of the three cysteine side chains and to establish whether the free S-H groups form intrahelical hydrogen bonds in the pentameric complex. Vibrational spectra of a transmembrane peptide were obtained which corresponded to the transmembrane domain of wild-type phospholamban and three peptides each containing a cysteine ⇒ alanine substitution. The observed S-H frequencies argue that each of the sulfhydryl groups is hydrogen-bonded to an i-4 backbone carbonyl oxygen. Electrostatic calculations on a model of phospholamban based on molecular dynamics and mutagenesis studies, show that the sulfhydryl groups may significantly contribute to the electrostatic potential field of the protein. Received: 22 July 1996/Revised: 10 October 1996     

Ligand-binding domains in vitellogenin receptors and other LDL-receptor family members share a common ancestral ordering of cysteine-rich repeats

Insect vitellogenin and yolk protein receptors (VgR/YPR) are newly discovered members of the low-density lipoprotein receptor (LDLR) family, which is characterized by a highly conserved arrangement of repetitive modular elements homologous to functionally unrelated proteins. The insect VgR/YPRs are unique in having two clusters of complement-type cysteine-rich (class A) repeats or modules, with five modules in the first cluster and seven in the second cluster, unlike classical LDLRs which have a single seven-module cluster, vertebrate VgRs and very low density lipoprotein receptors (VLDLR) which have a single eight-module cluster, and LDLR-related proteins (LRPs) and megalins which have four clusters of 2–7, 8, 10, and 11 modules. Alignment of clusters across subfamilies by conventional alignment programs is problematic because of the repetitive nature of the component modules which may have undergone rearrangements, duplications, and deletions during evolution. To circumvent this problem, we ``fingerprinted' each class A module in the different clusters by identifying those amino acids that are both relatively conserved and relatively unique within the cluster. Intercluster reciprocal comparisons of fingerprints and aligned sequences allowed us to distinguish four cohorts of modules reflecting shared recent ancestry. All but two of the 57 modules examined could be assigned to one of these four cohorts designated A, B, C, and D. Alignment of clusters based on modular cohorts revealed that all clusters are derived from a single primordial cluster of at least seven modules with a consensus arrangement of CDCADBC. All extant clusters examined are consistent with this consensus, though none matches it perfectly. This analysis also revealed that the eight-module clusters in vertebrate VgRs, insect VgR/YPRs, and LRP/megalins are not directly homologous with one another. Assignment of modules to cohorts permitted us to properly align 32 class A clusters from all four LDLR subfamilies for phylogenetic analysis. The results revealed that smaller one-cluster and two-cluster members of the family did not originate from the breakup of a large two-cluster or four-cluster receptor. Similarly, the LRP/megalins did not arise from the duplication of a two-cluster insect VgR/YPR-like progenitor. Rather, it appears that the multicluster receptors were independently constructed from the same single-cluster ancestor. Received: 16 January 1997 / Accepted: 21 August 1997     

Role and Source of ATP for Activation of Nonselective Cation Channels by AlF Complex in Guinea Pig Chromaffin Cells

Intracellular dialysis with the solution containing the G protein activator, AlF complex, induced an inward nonselective cation current (I _NS) at −55 mV in chromaffin cells. Amplitudes of I _NS induced by dialysis with ATP-free AlF solutions progressively diminished as cells were pretreated with cyanide, a mitochondrial inhibitor. After a 10-min pretreatment, generation of I _NS by the AlF complex depended on exogenous ATP delivered from pipette solution. The relationship between amplitudes of I _NS and concentrations of MgATP was well expressed by a rectangular hyperbola with an EC₅₀ of 0.265 mm. This result suggests that the cyanide treatment almost depleted ATP near the plasma membrane. On the other hand, a similar cyanide treatment of adrenal medullary preparations did not induce a marked decrease in cellular ATP content. GTP, ITP, or UTP could not substitute for ATP in generation of I _NS by the AlF complex. Similarly, the substitution of ATP with non- or poorly hydrolyzable ATP analogues did not aid in generating I _NS. Bath application of the kinase inhibitor, H-7 (100 μm), suppressed AlF-induced I _NS in a manner depending on intracellular Mg²⁺. We conclude that ATP is a prerequisite for generation of I _NS as a phosphoryl donor and that mitochondria is the main source of ATP. Received: 17 April 1996/Revised: 26 July 1996     

The Dentin Matrix Protein 1 Gene of Prototherian and Metatherian Mammals

Mineralization of tooth dentin (the deposition of hydroxyapatite crystals in and around collagen type I fibers of the extracellular matrix) requires the involvement of several genes, among them the gene coding for the dentin matrix protein 1, DMP1. We determined the exon–intron organization of the cattle DMP1 gene and used this information to amplify by the polymerase chain reaction homologous gene fragments from the genomic DNA of two species of metatherian (marsupial) mammals and one prototherian (monotreme) species. The translated proto- and metatherian protein sequences are highly divergent from the eutherian sequences but retain the general characteristics of the DMP1 (high acidity, serine-richness, multiple glycosylation sites, and the presence of the RGD cell attachment tripeptide). They therefore appear to be functional even though, evolutionarily, teeth are in a regression phase in prototherians. It is possible, therefore, that DMP1 is also involved in other functions besides dentinogenesis. The DMP1 gene appears to evolve rapidly and apparently tolerates non-frame-shifting insertions/deletions throughout the coding sequence. Received: 22 February 1998 / Accepted: 11 July 1998