The Effect of Branch Length Variation on the Selection of Models of Molecular Evolution

Models of sequence evolution play an important role in molecular evolutionary studies. The use of inappropriate models of evolution may bias the results of the analysis and lead to erroneous conclusions. Several procedures for selecting the best-fit model of evolution for the data at hand have been proposed, like the likelihood ratio test (LRT) and the Akaike (AIC) and Bayesian (BIC) information criteria. The relative performance of these model-selecting algorithms has not yet been studied under a range of different model trees. In this study, the influence of branch length variation upon model selection is characterized. This is done by simulating sequence alignments under a known model of nucleotide substitution, and recording how often this true model is recovered by different model-fitting strategies. Results of this study agree with previous simulations and suggest that model selection is reasonably accurate. However, different model selection methods showed distinct levels of accuracy. Some LRT approaches showed better performance than the AIC or BIC information criteria. Within the LRTs, model selection is affected by the complexity of the initial model selected for the comparisons, and only slightly by the order in which different parameters are added to the model. A specific hierarchy of LRTs, which starts from a simple model of evolution, performed overall better than other possible LRT hierarchies, or than the AIC or BIC. Received: 2 October 2000 / Accepted: 4 January 2001     

Phylogenetic Analysis of Reptilian Hemoglobins: Trees, Rates, and Divergences

Phylogenetic relationships among reptiles were examined using previously published and newly determined hemoglobin sequences. Trees reconstructed from these sequences using maximum-parsimony, neighbor-joining, and maximum-likelihood algorithms were compared with a phylogenetic tree of Amniota, which was assembled on the basis of published morphological data. All analyses differentiated α chains into α^A and α^D types, which are present in all reptiles except crocodiles, where only α^A chains are expressed. The occurrence of the α^D chain in squamates (lizards and snakes only in this study) appears to be a general characteristic of these species. Lizards and snakes also express two types of β chains (βI and βII), while only one type of β chain is present in birds and crocodiles. Reconstructed hemoglobin trees for both α and β sequences did not yield the monophyletic Archosauria (i.e., crocodilians + birds) and Lepidosauria (i.e., Sphenodon+ squamates) groups defined by the morphology tree. This discrepancy, as well as some other poorly resolved nodes, might be due to substantial heterogeneity in evolutionary rates among single hemoglobin lineages. Estimation of branch lengths based on uncorrected amino acid substitutions and on distances corrected for multiple substitutions (PAM distances) revealed that relative rates for squamate α^A and α^D chains and crocodilian β chains are at least twice as high as those of the rest of the chains considered. In contrast to these rate inequalities between reptilian orders, little variation was found within squamates, which allowed determination of absolute evolutionary rates for this subset of hemoglobins. Rate estimates for hemoglobins of lizards and snakes yielded 1.7 (α^A) and 3.3 (β) million years/PAM when calibrated with published divergence time vs. PAM distance correlates for several speciation events within snakes and for the squamate ↔ sphenodontid split. This suggests that hemoglobin chains of squamate reptiles evolved ∼3.5 (α^A) or ∼1.7 times (β) faster than their mammalian equivalents. These data also were used to obtain a first estimate of some intrasquamate divergence times. Received: 15 September 1997 / Accepted: 4 February 1998     

Implications of Fluctuations in Substitution Rates: Impact on the Uncertainty of Branch Lengths and on Relative-Rate Tests

Many tests of the lineage dependence of substitution rates, computations of the error of evolutionary distances, and simulations of molecular evolution assume that the rate of evolution is constant in time within each lineage descended from a common ancestor. However, estimates of the index of dispersion of numbers of mammalian substitutions suggest that the rate has time-dependent variations consistent with a fractal-Gaussian-rate Poisson process, which assumes common descent without assuming rate constancy. While this model does not affect certain relative-rate tests, it substantially increases the uncertainty of branch lengths. Thus, fluctuations in the rate of substitution cannot be neglected in calculations that rely on evolutionary distances, such as the confidence intervals of divergence times and certain phylogenetic reconstructions. The fractal-Gaussian-rate Poisson process is compared and contrasted with previous models of molecular evolution, including other Poisson processes, the fractal renewal process, a Lévy-stable process, a fractional-difference process, and a log-Brownian process. The fractal models are more compatible with mammalian data than the nonfractal models considered, and they may also be better supported by Darwinian theory. Although the fractal-Gaussian-rate Poisson process has not been proven to have better agreement with data or theory than the other fractal models, its Gaussian nature simplifies the exploration of its impact on evolutionary distance errors and relative-rate tests. Received: 29 September 1999 / Accepted: 20 January 2000     

Rates of Molecular Evolution in RNA Viruses: A Quantitative Phylogenetic Analysis

The study of rates of nucleotide substitution in RNA viruses is central to our understanding of their evolution. Herein we report a comprehensive analysis of substitution rates in 50 RNA viruses using a recently developed maximum likelihood phylogenetic method. This analysis revealed a significant relationship between genetic divergence and isolation time for an extensive array of RNA viruses, although more rate variation was usually present among lineages than would be expected under the constraints of a molecular clock. Despite the lack of a molecular clock, the range of statistically significant variation in overall substitution rates was surprisingly narrow for those viruses where a significant relationship between genetic divergence and time was found, as was the case when synonymous sites were considered alone, where the molecular clock was rejected less frequently. An analysis of the ecological and genetic factors that might explain this rate variation revealed some evidence of significantly lower substitution rates in vector-borne viruses, as well as a weak correlation between rate and genome length. Finally, a simulation study revealed that our maximum likelihood estimates of substitution rates are valid, even if the molecular clock is rejected, provided that sufficiently large data sets are analyzed. Received: 23 February 2001 / Accepted: 3 July 2001     

Dynamics and Phylogenetic Implications of MtDNA Control Region Sequences in New World Jays (Aves: Corvidae)

To study the evolution of mtDNA and the intergeneric relationships of New World Jays (Aves: Corvidae), we sequenced the entire mitochondrial DNA control region (CR) from 21 species representing all genera of New World jays, an Old World jay, crows, and a magpie. Using maximum likelihood methods, we found that both the transition/transversion ratio (κ) and among site rate variation (α) were higher in flanking domains I and II than in the conserved central domain and that the frequency of indels was highest in domain II. Estimates of κ and α were much more influenced by the density of taxon sampling than by alternative optimal tree topologies. We implemented a successive approximation method incorporating these parameters into phylogenetic analysis. In addition we compared our study in detail to a previous study using cytochrome b and morphology to examine the effect of taxon sampling, evolutionary rates of genes, and combined data on tree resolution. We found that the particular weighting scheme used had no effect on tree topology and little effect on tree robustness. Taxon sampling had a significant effect on tree robustness but little effect on the topology of the best tree. The CR data set differed nonsignificantly from the tree derived from the cytochrome b/morphological data set primarily in the placement of the genus Gymnorhinus, which is near the base of the CR tree. However, contrary to conventional taxonomy, the CR data set suggested that blue and black jays (Cyanocorax sensu lato) might be paraphyletic and that the brown jay Psilorhinus (=Cyanocorax) morio is the sister group to magpie jays (Calocitta), a phylogenetic hypothesis that is likely as parsimonious with regard to nonmolecular characters as monophyly of Cyanocorax. The CR tree also suggests that the common ancestor of NWJs was likely a cooperative breeder. Consistent with recent systematic theory, our data suggest that DNA sequences with high substitution rates such as the CR may nonetheless be useful in reconstructing relatively deep phylogenetic nodes in avian groups. Received: 10 November 1999 / Accepted: 16 March 2000     

Molecular Phylogenetic Inference of the Woolly Mammoth Mammuthus primigenius, Based on Complete Sequences of Mitochondrial Cytochrome b and 12S Ribosomal RNA Genes

Complete sequences of cytochrome b (1,137 bases) and 12S ribosomal RNA (961 bases) genes in mitochondrial DNA were successfully determined from the woolly mammoth (Mammuthus primigenius), African elephant (Loxodonta africana), and Asian elephant (Elephas maximus). From these sequence data, phylogenetic relationships among three genera were examined. Molecular phylogenetic trees reconstructed by the neighbor-joining and the maximum parsimony methods provided an identical topology both for cytochrome b and 12S rRNA genes. These results support the ``Mammuthus-Loxodonta' clade, which is contrary to some previous morphological reports that Mammuthus is more closely related to Elephas than to Loxodonta. Received: 8 April 1997 / Accepted: 23 July 1997     

The Evolution of Codon Preferences in Drosophila: A Maximum-Likelihood Approach to Parameter Estimation and Hypothesis Testing

Synonymous codon usage in related species may differ as a result of variation in mutation biases, differences in the overall strength and efficiency of selection, and shifts in codon preference—the selective hierarchy of codons within and between amino acids. We have developed a maximum-likelihood method to employ explicit population genetic models to analyze the evolution of parameters determining codon usage. The method is applied to twofold degenerate amino acids in 50 orthologous genes from D. melanogaster and D. virilis. We find that D. virilis has significantly reduced selection on codon usage for all amino acids, but the data are incompatible with a simple model in which there is a single difference in the long-term N _e, or overall strength of selection, between the two species, indicating shifts in codon preference. The strength of selection acting on codon usage in D. melanogaster is estimated to be |N _e s|≈ 0.4 for most CT-ending twofold degenerate amino acids, but 1.7 times greater for cysteine and 1.4 times greater for AG-ending codons. In D. virilis, the strength of selection acting on codon usage for most amino acids is only half that acting in D. melanogaster but is considerably greater than half for cysteine, perhaps indicating the dual selection pressures of translational efficiency and accuracy. Selection coefficients in orthologues are highly correlated (ρ= 0.46), but a number of genes deviate significantly from this relationship. Received: 20 December 1998 / Accepted: 17 February 1999     

Evolutionary Rate and Genetic Heterogeneity of Human T-Cell Lymphotropic Virus Type II (HTLV-II) Using Isolates from European Injecting Drug Users

Seven new Italian and two new British HTLV-II isolates were obtained from injecting drug users and the entire long terminal repeat (LTR) region was sequenced. Restriction analysis showed that all the Italian isolates are of the IIb subtype, whereas the British isolates are of the IIa subtype. To understand whether the further differentiation of each two principal HTLV-II subtypes in several subgroups could be statistically supported by phylogenetic analysis, the neighbor-joining, parsimony, and maximum likelihood methods were used. The separation between IIa and IIb is very well supported by all three methods. At least two phylogenetic subgroups exist within the HTLV-IIa and at least three within the HTLV-IIb subtype. In the present analysis, no statistical support was obtained for additional phylogroups. Two particular subgroups seem interesting because they include all European and North American injecting drug user strains within the IIa and IIb subtypes, respectively. These data confirm that European HTLV-II infection among drug users is probably derived from North America. They also suggest that though a certain differentiation by restriction analysis in different subgroups is possible, carefully interpreted phylogenetic analyses remain necessary. Using the likelihood ratio test, a molecular clock for the drug user strains was calibrated. A fixation rate between 1.08 × 10⁻⁴ and 2.7 × 10⁻⁵ nucleotide substitutions per site per year was calculated for the IIa and IIb injecting drug user strains. This is the lowest fixation rate so far reported for RNA viruses, including for HIV, which typically range between 10⁻² and 10⁻⁴.     

Efficiency of the Neighbor-Joining Method in Reconstructing Deep and Shallow Evolutionary Relationships in Large Phylogenies

The neighbor-joining (NJ) method is widely used in reconstructing large phylogenies because of its computational speed and the high accuracy in phylogenetic inference as revealed in computer simulation studies. However, most computer simulation studies have quantified the overall performance of the NJ method in terms of the percentage of branches inferred correctly or the percentage of replications in which the correct tree is recovered. We have examined other aspects of its performance, such as the relative efficiency in correctly reconstructing shallow (close to the external branches of the tree) and deep branches in large phylogenies; the contribution of zero-length branches to topological errors in the inferred trees; and the influence of increasing the tree size (number of sequences), evolutionary rate, and sequence length on the efficiency of the NJ method. Results show that the correct reconstruction of deep branches is no more difficult than that of shallower branches. The presence of zero-length branches in realized trees contributes significantly to the overall error observed in the NJ tree, especially in large phylogenies or slowly evolving genes. Furthermore, the tree size does not influence the efficiency of NJ in reconstructing shallow and deep branches in our simulation study, in which the evolutionary process is assumed to be homogeneous in all lineages. Received: 7 March 2000 / Accepted: 2 August 2000     

Phylogenetic Characterization of the MIP Family of Transmembrane Channel Proteins

The ubiquitous major intrinsic protein (MIP) family includes several transmembrane channel proteins known to exhibit specificity for water and/or neutral solutes. We have identified 84 fully or partially sequenced members of this family, have multiply aligned over 50 representative, divergent, fully sequenced members, have used the resultant multiple alignment to derive current MIP family-specific signature sequences, and have constructed a phylogenetic tree. The tree reveals novel features relevant to the evolutionary history of this protein family. These features plus an evaluation of functional studies lead to the postulates: (i) that all current MIP family proteins derived from two divergent bacterial paralogues, one a glycerol facilitator, the other an aquaporin, and (ii) that most or all current members of the family have retained these or closely related physiological functions. Received: 19 April 1996/Revised: 3 June 1996     

Phylogenetic Analysis of the Friedreich Ataxia GAA Trinucleotide Repeat

Friedreich ataxia is an autosomal recessive neurodegenerative disorder associated with a GAA repeat expansion in the first intron of the gene (FRDA) encoding a novel, highly conserved, 210 amino acid protein known as frataxin. Normal variation in repeat size was determined by analysis of more than 600 DNA samples from seven human populations. This analysis showed that the most frequent allele had nine GAA repeats, and no alleles with fewer than five GAA repeats were found. The European and Syrian populations had the highest percentage of alleles with 10 or more GAA repeats, while the Papua New Guinea population did not have any alleles carrying more than 10 GAA repeats. The distributions of repeat sizes in the European, Syrian, and African American populations were significantly different from those in the Asian and Papua New Guinea populations (p < 0.001). The GAA repeat size was also determined in five nonhuman primates. Samples from 10 chimpanzees, 3 orangutans, 1 gorilla, 1 rhesus macaque, 1 mangabey, and 1 tamarin were analyzed. Among those primates belonging to the Pongidae family, the chimpanzees were found to carry three or four GAA repeats, the orangutans had four or five GAA repeats, and the gorilla carried three GAA repeats. In primates belonging to the Cercopithecidae family, three GAA repeats were found in the mangabey and two in the rhesus macaque. However, an AluY subfamily member inserted in the poly(A) tract preceding the GAA repeat region in the rhesus macaque, making the amplified sequence approximately 300 bp longer. The GAA repeat was also found in the tamarin, suggesting that it arose at least 40 million years ago and remained relatively small throughout the majority of primate evolution, with a punctuated expansion in the human genome. Received: 18 August 2000 / Accepted: 10 November 2000     

Variation in the Pattern of Nucleotide Substitution Across Sites

A model of nucleotide substitution that allows the transition/transversion rate bias to vary across sites was constructed. We examined the fit of this model using likelihood-ratio tests by analyzing 13 protein coding genes and 1 pseudogene. Likelihood-ratio testing indicated that a model that allows variation in the transition/transversion rate bias across sites provided a significant improvement in fit for most protein coding genes but not for the pseudogene. When the analysis was repeated with parameters estimated separately for first, second, and third codon positions, strong heterogeneity was uncovered for the first and second codon positions; the variation in the transition/transversion rate was generally weaker at the third codon position. The transition rate bias and branch lengths are underestimated when variation in the transition/transversion rate was not accommodated, suggesting that it may be important to accommodate variation in the pattern of nucleotide substitution for accurate estimation of evolutionary parameters. Received: 4 November 1997 / Accepted: 19 May 1998     

The Effect of Recombination on the Accuracy of Phylogeny Estimation

Phylogenetic studies based on DNA sequences typically ignore the potential occurrence of recombination, which may produce different alignment regions with different evolutionary histories. Traditional phylogenetic methods assume that a single history underlies the data. If recombination is present, can we expect the inferred phylogeny to represent any of the underlying evolutionary histories? We examined this question by applying traditional phylogenetic reconstruction methods to simulated recombinant sequence alignments. The effect of recombination on phylogeny estimation depended on the relatedness of the sequences involved in the recombinational event and on the extent of the different regions with different phylogenetic histories. Given the topologies examined here, when the recombinational event was ancient, or when recombination occurred between closely related taxa, one of the two phylogenies underlying the data was generally inferred. In this scenario, the evolutionary history corresponding to the majority of the positions in the alignment was generally recovered. Very different results were obtained when recombination occurred recently among divergent taxa. In this case, when the recombinational breakpoint divided the alignment in two regions of similar length, a phylogeny that was different from any of the true phylogenies underlying the data was inferred.     

A Phylogenetic Assessment of the Eukaryotic Light-Harvesting Antenna Proteins, with Implications for Plastid Evolution

The light-harvesting complexes (LHCs) are a superfamily of chlorophyll-binding proteins present in all photosynthetic eukaryotes. The Lhc genes are nuclear-encoded, yet the pigment–protein complexes are localized to the thylakoid membrane and provide a marker to follow the evolutionary paths of plastids with different pigmentation. The LHCs are divided into the chlorophyll a/b-binding proteins of the green algae, euglenoids, and higher plants and the chlorophyll a/c-binding proteins of various algal taxa. This work examines the phylogenetic position of the LHCs from three additional taxa: the rhodophytes, the cryptophytes, and the chlorarachniophytes. Phylogenetic analysis of the LHC sequences provides strong statistical support for the clustering of the rhodophyte and cryptomonad LHC sequences within the chlorophyll a/c-binding protein lineage, which includes the fucoxanthin–chlorophyll proteins (FCP) of the heterokonts and the intrinsic peridinin–chlorophyll proteins (iPCP) of the dinoflagellates. These associations suggest that plastids from the heterokonts, haptophytes, cryptomonads, and the dinoflagellate, Amphidinium, evolved from a red algal-like ancestor. The Chlorarachnion LHC is part of the chlorophyll a/b-binding protein assemblage, consistent with pigmentation, providing further evidence that its plastid evolved from a green algal secondary endosymbiosis. The Chlorarachnion LHC sequences cluster with the green algal LHCs that are predominantly associated with photosystem II (LHCII). This suggests that the green algal endosymbiont that evolved into the Chlorarachnion plastid was acquired following the emergence of distinct LHCI and LHCII complexes. Received: 25 February 1998 / Accepted: 13 May 1998     

Substitution Model of Sequence Evolution for the Human Immunodeficiency Virus Type 1 Subtype B gp120 Gene over the C2-V5 Region

Phylogenetic analyses frequently rely on models of sequence evolution that detail nucleotide substitution rates, nucleotide frequencies, and site-to-site rate heterogeneity. These models can influence hypothesis testing and can affect the accuracy of phylogenetic inferences. Maximum likelihood methods of simultaneously constructing phylogenetic tree topologies and estimating model parameters are computationally intensive, and are not feasible for sample sizes of 25 or greater using personal computers. Techniques that initially construct a tree topology and then use this non-maximized topology to estimate ML substitution rates, however, can quickly arrive at a model of sequence evolution. The accuracy of this two-step estimation technique was tested using simulated data sets with known model parameters. The results showed that for a star-like topology, as is often seen in human immunodeficiency virus type 1 (HIV-1) subtype B sequences, a random starting topology could produce nucleotide substitution rates that were not statistically different than the true rates. Samples were isolated from 100 HIV-1 subtype B infected individuals from the United States and a 620 nt region of the env gene was sequenced for each sample. The sequence data were used to obtain a substitution model of sequence evolution specific for HIV-1 subtype B env by estimating nucleotide substitution rates and the site-to-site heterogeneity in 100 individuals from the United States. The method of estimating the model should provide users of large data sets with a way to quickly compute a model of sequence evolution, while the nucleotide substitution model we identified should prove useful in the phylogenetic analysis of HIV-1 subtype B env sequences. Received: 4 October 2000 / Accepted: 1 March 2001     

Phylogenetic Studies of Complete Mitochondrial DNA Molecules Place Cartilaginous Fishes Within the Tree of Bony Fishes

It is commonly acknowledged that cartilaginous fishes, Chondrichthyes, have a basal position among the Gnathostomata (jawed vertebrates). In order to explore this relationship we have sequenced the complete mitochondrial genome of the spiny dogfish, Squalus acanthias, and included it in a phylogenetic analysis together with a number of bony fishes and amniotes. The phylogenetic reconstructions placed the dogfish among the bony fishes. Thus, and contrary to the common view, the analyses have shown that the position of the sharks is not basal among the gnathostomes. The presently recognized phylogenetic position of the dogfish was identified irrespective of the outgroup used, echinoderms or agnathan fishes. The lungfish was the most basal gnathostome fish, while the teleosteans had an apical position in the piscine tree. A basal position of the dogfish among the gnathostomes was statistically rejected, but the phylogenetic relationship among the coelacanth, spiny dogfish, and teleosts was not conclusively resolved. The findings challenge the current theory that sharks and other chondrichthyans, if monophyletic, are the sister group to all other extant gnathostomes. The results open to question the status of several morphological characters commonly used in piscine phylogenetic reconstruction, most notably the presence versus absence of endochondral bone in the endoskeleton, the macromeric versus micromeric structure of the exoskeleton, and the presence/absence of swimbladder and/or lung. The study also confirmed recent findings demonstrating that the origin of the amniotes is deeper than the diversification of extant bony fishes. Received: 12 March 1998 / Accepted: 12 June 1998     

Toward Assigning Helical Regions in Alignments of Ribosomal RNA and Testing the Appropriateness of Evolutionary Models

We suggest a nucleotide substitution model that takes correlation between base-paired nucleotides into account. The model includes the estimation of the transition–transversion ratio and allows inference of the shape parameter of a discrete gamma distribution to include rate heterogeneity. A Cox-test statistic, applied to a diatom ribosomal RNA alignment, shows that the suggested correlation model explains evolution of the stem region better than usual independence models. Moreover, the Cox-test procedure is extended to shed some light upon the problem of assigning helical regions in a secondary structure based alignment. This approach provides an estimate of the percentage of stem positions that do not appear to be correlated. Received: 4 March 1999 / Accepted: 10 May 1999     

Positive Darwinian Selection Promotes Heterogeneity Among Members of the Antifreeze Protein Multigene Family

A variety of organisms have independently evolved proteins exhibiting antifreeze activity that allows survival at subfreezing temperatures. The antifreeze proteins (AFPs) bind ice nuclei and depress the freezing point by a noncolligative absorption–inhibition mechanism. Many organisms have a heterogeneous suite of AFPs with variation in primary sequence between paralogous loci. Here, we demonstrate that the diversification of the AFP paralogues is promoted by positive Darwinian selection in two independently evolved AFPs from fish and beetle. First, we demonstrate an elevated rate of nonsynonymous substitutions compared to synonymous substitutions in the mature protein coding region. Second, we perform phylogeny-based tests of selection to demonstrate a subset of codons is subjected to positive selection. When mapped onto the three-dimensional structure of the fish antifreeze type III antifreeze structure, these codons correspond to amino acid positions that surround but do not interrupt the putative ice-binding surface. The selective agent may be related to efficient binding to diverse ice surfaces or some other aspect of AFP function. Received: 27 February 2001 / Accepted: 12 September 2001     

Muscarinic Activation of BK Channels Induces Membrane Oscillations in Glioma Cells and Leads to Inhibition of Cell Migration

Patients with cerebral tumors often present with elevated levels of acetylcholine (ACh) in their cerebrospinal fluid. This motivated us to investigate physiological effects of ACh on cultured human astrocytoma cells (U373) using a combination of videomicroscopy, calcium microspectrofluorimetry and perforated patch-clamp recording. Astrocytoma cells exhibited the typical morphological changes associated with cell migration; polarized cells displayed prominent lamellipodia and associated membrane ruffling at the anterior of the cell, and a long tail region that periodically contracted into the cell body as the cell moved forward. Bath application of the ACh receptor agonist, muscarine, reversibly inhibited cell migration. In conjunction with this inhibition, ACh induced a dose-dependent, biphasic increase in resting intracellular free calcium concentration ([Ca²⁺] _i ) associated with periodic Ca²⁺ oscillations during prolonged ACh applications. The early transient rise in [Ca²⁺] _i was abolished by ionomycin and thapsigargin but was insensitive to caffeine and ryanodine while the plateau phase was strictly dependent on external calcium. The Ca²⁺ response to ACh was mimicked by muscarine and abolished by the muscarinic antagonists, atropine or 4-DAMP, but not by pirenzepine. Using perforated patch-clamp recordings combined with fluorescent imaging, we demonstrated that ACh-induced [Ca²⁺] _i oscillations triggered membrane voltage oscillations that were due to the activation of voltage-dependent, Ca²⁺-sensitive K⁺ currents. These K⁺ currents were blocked by intracellular injection of EGTA, or by extracellular application of TEA, quinine, or charybdotoxin, but not by apamin. These studies suggest that activation of muscarinic receptors on glioma cells induce the release of Ca²⁺ from intracellular stores which in turn activate Ca²⁺-dependent (BK-type) K⁺ channels. Furthermore, this effect was associated with inhibition of cell migration, suggesting an interaction of this pathway with glioma cell migration. Received: 17 December/Revised: 17 March 2000     

Phylogenetic Analysis of Invertebrate Lysozymes and the Evolution of Lysozyme Function

We isolated and sequenced the cDNAs coding for lysozymes of six bivalve species. Alignment and phylogenetic analysis showed that, together with recently described bivalve lysozymes, the leech destabilase, and a number of putative proteins from extensive genomic and cDNA analyses, they belong to the invertebrate type of lysozymes (i type), first described by Jollès and Jollès (1975). We determined the genomic structure of the gene encoding the lysozyme of Mytilus edulis, the common mussel. We provide evidence that the central exon of this gene is homologous to the second exon of the chicken lysozyme gene, belonging to the c type. We propose that the origin of this domain can be traced back in evolution to the origin of bilaterian animals. Phylogenetic analysis suggests that i-type proteins form a monophyletic family. Received: 21 May 2001 / Accepted: 22 October 2001