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     

Evidence for Involvement of a Zymogen Granule Na+/H+ Exchanger in Enzyme Secretion from Rat Pancreatic Acinar Cells

We have characterized a Na⁺/H⁺ exchanger in the membrane of isolated zymogen granules (ZG) from rat exocrine pancreas and investigated its role in secretagogue-induced enzyme secretion. ZG Na⁺/H⁺ exchanger activity was estimated by measuring Na⁺ or Li⁺ influx and consequent osmotic swelling and lysis of ZG incubated in Na- or Li-acetate. Alternatively, intragranule pH was investigated by measuring absorbance changes in ZG which had been preloaded with the weak base acridine orange. Na⁺- or Li⁺-dependent ZG lysis was enhanced by increasing inward to outward directed H⁺ gradients. Na⁺-dependent ZG lysis was not prevented by an inside-positive K⁺ diffusion potential generated by valinomycin which argues against parallel operation of separate electrogenic Na⁺ and H⁺ permeabilities and for coupled Na⁺/H⁺ exchange through an electroneutral carrier. Na⁺- and Li⁺-dependent ZG lysis was inhibited by EIPA (EC₅₀∼25 μm) and benzamil (EC₅₀∼100 μm), but only weakly by amiloride. Similarly, absorbance changes due to release of acridine orange from acidic granules into the medium were obtained with Na⁺ and Li⁺ salts only, and were inhibited by EIPA, suggesting the presence of a Na⁺/H⁺ exchanger in the membrane. Na⁺ dependent lysis of ZG was inhibited by 0.5 mm MgATP and MgATP-γ-S by about 60% and 35%, respectively. Inhibition by MgATP was prevented by incubation of ZG with alkaline phosphatase (100 U/ml), or by the calmodulin antagonists calmidazolium (0.75 μm), trifluoperazine (100 μm) and W-7 (500 μm), suggesting that the ZG Na⁺/H⁺ exchanger is regulated by a ZG membrane-bound calmodulin-dependent protein kinase. Na⁺ dependence of secretagogue (CCK-OP)-stimulated amylase secretion was investigated in digitonin permeabilized rat pancreatic acini and was higher in acini incubated in Na⁺ containing buffer (30 mm NaCl/105 mm KCl buffer; 6.4 ± 0.4% of total amylase above basal) compared to buffer without Na⁺ (0 mm NaCl/135 mm KCl buffer; 4.7 ± 0.4% of total amylase above basal, P < 0.03). EIPA (50 μm) reduced CCK-OP-induced amylase secretion in Na⁺ containing buffer from 7.5 ± 0.6% to 4.1 ± 0.8% (P < 0.02). In the absence of Na⁺ in the buffer, CCK-OP-stimulated amylase release was not inhibited by 50 μm EIPA. The data suggest that an amiloride insensitive, EIPA inhibitable Na⁺/H⁺ exchanger is present in ZG membranes, which is stimulated by calmodulin antagonists and could be involved in secretagogue-induced enzyme secretion from rat pancreatic acini. Received: 7 December 1995/Revised: 2 April 1996     

Gene Conversion and Evolution of Daphniid Hemoglobins (Crustacea, Cladocera)

The extracellular hemoglobins of cladocerans derive from the aggregation of 12 two-domain globin subunits that are apparently encoded by four genes. This study establishes that at least some of these genes occur as a tandem array in both Daphnia magna and Daphnia exilis. The genes share a uniform structure; a bridge intron separates two globin domains which each include three exons and two introns. Introns are small, averaging just 77 bp, but a longer sequence (2.2–3.2 kb) separates adjacent globin genes. A survey of structural diversity in globin genes from other daphniids revealed three independent cases of intron loss, but exon lengths were identical, excepting a 3-bp insertion in exon 5 of Simocephalus. Heterogeneity in the extent of nucleotide divergence was marked among exons, largely as a result of the pronounced diversification of the terminal exon. This variation reflected, in part, varying exposure to concerted evolution. Conversion events were frequent in exons 1–4 but were absent from exons 5 and 6. Because of this difference, the results of phylogenetic analyses were strongly affected by the sequences employed in this construction. Phylogenies based on total nucleotide divergence in exons 1–4 revealed affinities among all genes isolated from a single species, reflecting the impact of gene conversion events. In contrast, phylogenies based on total nucleotide divergence in exons 5 and 6 revealed affinities among orthologous genes from different taxa. Received: 8 March 1999 / Accepted: 14 July 1999     

Molecular Evolution of Calmodulin and Calmodulin-Like Genes in the Cephalochordate Branchiostoma

Calmodulin is a calcium-binding EF-hand protein that is an activator of many enzymes as well as ion pumps and channels. Due to its multiple targets and its central role in the cell, understanding the evolutionary history of calmodulin genes should provide insights into the origin of genetic complexity in eukaryotes. We have previously isolated and characterized a calmodulin gene from the early-diverging chordate Branchiostoma lanceolatum (CaM1). In this paper, we report the existence of a second calmodulin gene (CaM2) as well as two CaM-like genomic fragments (CaML-2, CaML-3) in B. lanceolatum and a CaM2 and three CaM-like genes (CaML-1, CaML-2, CaML-3) in B. floridae. The CaM-like genes were isolated using low-stringency PCR. Surprisingly, the nucleotide sequences of the B. lanceolatum CaM1 and CaM2 cDNAs differ by 19.3%. Moreover, the CaM2 protein differs at two positions from the amino acid sequence of CaM1; the latter is identical to calmodulins in Drosophila melanogaster, the mollusc Aplysia californica, and the tunicate Halocynthia roretzi. The two B. lanceolatum CaM-like genes are more closely related to the CaM2 than to the CaM1 gene. This relationship is supported by the phylogenetic analyses and the identical exon/intron organization of these three genes, a relationship unique among animal CaM sequences. These data demonstrate the existence of a CaM multigene family in the cephalochordate Branchiostoma, which may have evolved independently from the multigene family in vertebrates. Received: 2 November 1999 / Accepted: 25 April 2000     

Molecular evolution of a portion of the mitochondrial 16S ribosomal gene region in scleractinian corals

Relationships among families and suborders of scleractinian corals are poorly understood because of difficulties 1) in making inferences about the evolution of the morphological characters used in coral taxonomy and 2) in interpreting their 240-million-year fossil record. Here we describe patterns of molecular evolution in a segment of the mitochondrial (mt) 16S ribosomal gene from taxa of 14 families of corals and the use of this gene segment in a phylogenetic analysis of relationships within the order. We show that sequences obtained from scleractinians are homologous to other metazoan 16S ribosomal sequences and fall into two distinct clades defined by size of the amplified gene product. Comparisons of sequences from the two clades demonstrate that both sets of sequences are evolving under similar evolutionary constraints: they do not differ in nucleotide composition, numbers of transition and transversion substitutions, spatial patterns of substitutions, or in rates of divergence. The characteristics and patterns observed in these sequences as well as the secondary structures, are similar to those observed in mt 16S ribosomal DNA sequences from other taxa. Phylogenetic analysis of these sequences shows that they are useful for evaluating relationships within the order. The hypothesis generated from this analysis differs from traditional hypotheses for evolutionary relationships among the Scleractinia and suggests that a reevaluation of evolutionary affinities in the order is needed. Received: 4 September 1996 / Accepted: 7 April 1997     

Rapid evolution of a young L1 (LINE-1) clade in recently speciated rattus taxa

L1 elements are retrotransposons that have been replicating and evolving in mammalian genomes since before the mammalian radiation. Rattus norvegicus shares the young L1_mlvi2 clade only with its sister taxon, Rattus cf moluccarius. Here we compared the L1_mlvi2 clade in these recently diverged species and found that it evolved rapidly into closely related but distinct clades: the L1_mlvi2-rm clade (or subfamily), characterized here from R. cf moluccarius, and the L1_mlvi2-rn clade, originally described in R. norvegicus. In addition to other differences, these clades are distinguished by a cluster of amino acid replacement substitutions in ORF I. Both rat species contain the L1_mlvi2-rm clade, but the L1_mlvi2-rn clade is restricted to R. norvegicus. Therefore, the L1_mlvi2-rm clade arose prior to the divergence of R. norvegicus and R. cf moluccarius, and the L1_mlvi2-rn clade amplified after their divergence. The total number of L1_mlvi2-rm elements in R. cf moluccarius is about the same as the sum of the L1_mlvi2-rm and L1_mlvi2-rn elements in R. norvegicus. The possibility that L1 amplification is in some way limited so that the two clades compete for replicative supremacy as well as the implications of the other distinguishing characteristic of the L1_mlvi2-rn and L1_mlvi2-rm clades are discussed. Received: 6 November 1996 / Accepted: 12 April 1997     

Origin and evolution of GBV-C/hepatitis G virus and relationships with ancient human migrations

The GB virus C/hepatitis G virus (GBV-C/HGV) is a newly identified human RNA virus, belonging to the Flaviviridae family. Persistent infection by GBV-C/HGV is common in humans, and genetically divergent isolates have been identified in different parts of the world. Due to the absence of a real pathogenic role of GBV-C/HGV in liver disease and its extremely low mutation rate, this virus is a potential marker to trace prehistoric links between human populations. In this study, origin and evolution of GBV-C/HGV were examined using a set of fully sequenced strains of worldwide origin. A first phylogenetic analysis, addressed to the short (255 nucleotides) NS5A overlapping coding region by the neighbor-joining method, suggested an ancient African origin of GBV-C/HGV. This notion was confirmed when the same analysis was applied to the genomic regions showing the lowest rate of synonymous substitutions, covering one-fourth (2184 nucleotides) of the total coding potential of the virus genome. By using a multivariate statistical method and extending the analysis to the complete coding region, fine details of the evolutionary history of GBV-C/HGV were further elucidated. By this approach, isolates from Southeast Asia appeared to be the most closely related to those of African origin, consistent with a major route of ancient human migrations from Africa to southeastern parts of the Asian continent. Received: 26 October 2000 / Accepted: 28 February 2001     

Modulation of the Calmodulin-induced Inhibition of Sarcoplasmic Reticulum Calcium Release Channel (Ryanodine Receptor) by Sulfhydryl Oxidation in Single Channel Current Recordings and [3H]Ryanodine Binding

The modulation of the calmodulin-induced inhibition of the calcium release channel (ryanodine receptor) by two sulfhydryl oxidizing compounds, 4-(chloromercuri)phenyl–sulfonic acid (4-CMPS) and 4,4′-dithiodipyridine (4,4′-DTDP) was determined by single channel current recordings with the purified and reconstituted calcium release channel from rabbit skeletal muscle sarcoplasmic reticulum (HSR) and [³H]ryanodine binding to HSR vesicles. 0.1 μm CaM reduced the open probability (P _o ) of the calcium release channel at maximally activating calcium concentrations (50–100 μm) from 0.502 ± 0.02 to 0.137 ± 0.022 (n= 28), with no effect on unitary conductance. 4-CMPS (10–40 μm) and 4,4′-DTDP (0.1–0.3 mm) induced a concentration dependent increase in P _o (> 0.9) and caused the appearance of longer open states. CaM shifted the activation of the calcium release channel by 4-CMPS or 4,4′-DTDP to higher concentrations in single channel recordings and [³H]ryanodine binding. 40 μm 4-CMPS induced a near maximal (P _o > 0.9) and 0.3 mm 4,4′-DTDP a submaximal (P _o = 0.74) channel opening in the presence of CaM, which was reversed by the specific sulfhydryl reducing agent DTT. Neither 4-CMPS nor 4,4′-DTDP affected Ca-[¹²⁵I]calmodulin binding to HSR. 1 mm MgCl₂ reduced P _o from 0.53 to 0.075 and 20–40 μm 4-CMPS induced a near maximal channel activation (P _o > 0.9). These results demonstrate that the inhibitory effect of CaM or magnesium in a physiological concentration is diminished or abolished at high concentrations of 4-CMPS or 4,4′-DTDP through oxidation of activating sulfhydryls on cysteine residues of the calcium release channel. Received: 22 July 1999/Revised: 15 November 1999     

The tryptophan biosynthetic pathway of aphid endosymbionts (Buchnera): Genetics and evolution of plasmid-associated anthranilate synthase (trpEG) within the aphididae

The bacterial endosymbionts (Buchnera) from the aphids Rhopalosiphum padi, R. maidis, Schizaphis graminum, and Acyrthosiphon pisum contain the genes for anthranilate synthase (trpEG) on plasmids made up of one or more 3.6-kb units. Anthranilate synthase is the first as well as the rate-limiting enzyme in the tryptophan biosynthetic pathway. The amplification of trpEG on plasmids may result in an increase of enzyme protein and overproduction of this essential amino acid, which is required by the aphid host. The nucleotide sequence of trpEG from endosymbionts of different species of aphids is highly conserved, as is an approximately 500-bp upstream DNA segment which has the characteristics of an origin of replication. Phylogenetic analyses were performed using trpE and trpG from the endosymbionts of these four aphids as well as from the endosymbiont of Schlechtendalia chinensis, in which trpEG occurs on the chromosome. The resulting phylogeny was congruent with trees derived from sequences of two chromosome-located bacterial genes (part of trpB and 16S ribosomal DNA). In turn, trees obtained from plasmid-borne and bacterial chromosome-borne sequences were congruent with the tree resulting from phylogenetic analysis of three aphid mitochondrial regions (portions of the small and large ribosomal DNA subunits, as well as cytochrome oxidase II). Congruence of trees based on genes from host mitochondria and from bacteria adds to previous support for exclusively vertical transmission of the endosymbionts within aphid lineages. Congruence with trees based on plasmid-borne genes supports the origin of the plasmid-borne trpEG from the chromosomal genes of the same lineage and the absence of subsequent plasmid exchange among endosymbionts of different species of aphids. Received: 22 August 1995 / Accepted: 6 September 1995     

The mitochondrial DNA molecule of the hagfish (myxine glutinosa) and vertebrate phylogeny

The vertebrates are traditionally classified into two distinct groups, Agnatha (jawless vertebrates) and Gnathostomata (jawed vertebrates). Extant agnathans are represented by hagfishes (Myxiniformes) and lampreys (Petromyzontiformes), frequently grouped together within the Cyclostomata. Whereas the recognition of the Gnathostomata as a clade is commonly acknowledged, a consensus has not been reached regarding whether or not Cyclostomata represents a clade. In the present study we have used newly established sequences of the protein-coding genes of the mitochondrial DNA molecule of the hagfish to explore agnathan and gnathostome relationships. The phylogenetic analysis of Pisces, using echinoderms as outgroup, placed the hagfish as a sister group of Vertebrata sensu stricto, i.e., the lamprey and the gnathostomes. The phylogenetic analysis of the Gnathostomata identified a basal divergence between gnathostome fishes and a branch leading to birds and mammals, i.e., between ``Anamnia' and Amniota. The lungfish has a basal position among gnathostome fishes with the teleosts as the most recently evolving lineage. The findings portray a hitherto unrecognized polarity in the evolution of bony fishes. The presently established relationships are incompatible with previous molecular studies. Received: 15 August 1997 / Accepted: 1 October 1997     

Regulation of Ca2+-activated K+ Channel from Rabbit Distal Colon Epithelium by Phosphorylation and Dephosphorylation

The Ca²⁺-activated maxi K⁺ channel is predominant in the basolateral membrane of the surface cells in the distal colon. It may play a role in the regulation of the aldosterone-stimulated Na⁺ reabsorption from the intestinal lumen. Previous measurements of these basolateral K⁺ channels in planar lipid bilayers and in plasma membrane vesicles have shown a very high sensitivity to Ca²⁺ with a K _0.5 ranging from 20 nm to 300 nm, whereas other studies have a much lower sensitivity to Ca²⁺. To investigate whether this difference could be due to modulation by second messenger systems, the effect of phosphorylation and dephosphorylation was examined. After addition of phosphatase, the K⁺ channels lost their high sensitivity to Ca²⁺, yet they could still be activated by high concentrations of Ca²⁺ (10 μm). Furthermore, the high sensitivity to Ca²⁺ could be restored after phosphorylation catalyzed by a cAMP dependent protein kinase. There was no effect of addition of protein kinase C. In agreement with the involvement of enzymatic processes, lag periods of 30–120 sec for dephosphorylation and of 10–280 sec for phosphorylation were observed. The phosphorylation state of the channel did not influence the single channel conductance. The results demonstrate that the high sensitivity to Ca²⁺ of the maxi K⁺ channel from rabbit distal colon is a property of the phosphorylated form of the channel protein, and that the difference in Ca²⁺ sensitivity between the dephosphorylated and phosphorylated forms of the channel protein is more than one order of magnitude. The variety in Ca²⁺ sensitivities for maxi K⁺ channels from tissue to tissue and from different studies on the same tissue could be due to modification by second messenger systems. Received: 28 February 1995/Revised: 22 December 1995     

Common features of analogous replacement histone H3 genes in animals and plants

Phylogenetic analysis of histone H3 protein sequences demonstrates the independent origin of the replacement histone H3 genes in animals and in plants. Multiple introns in the replacement histone H3 genes of animals in a pattern distinct from that in plant replacement H3 genes supports this conclusion. It is suggested that replacement H3 genes arose at the same time that, independently, multicellular forms of animals and of plants evolved. Judged by the degree of invariant and functionally constrained amino acid positions, histones H3 and H4, which form together the tetramer kernel of the nucleosome, have co-evolved with equal rates of sequence divergence. Residues 31 and 87 in histone H3 are the only residues that consistently changed across each gene duplication event that created functional replacement histone H3 variant forms. Once changed, these residues have remained invariant across divergent speciation. This suggests that they are required to allow replacement histone H3 to participate in the assembly of nucleosomes in non–S-phase cells. The abundant occurrence of polypyrimidine sequences in the introns of all replacement H3 genes, and the replacement of an intron by a polypyrimidine motif upstream of the alfalfa replacement H3 gene, suggests a function. It is speculated that they may contribute to the characteristic cell-cycle-independent pattern of replacement histone H3 genes by binding nucleosome-excluding proteins.     

Evolution of trophic types in emperor fishes ( Lethrinus,Lethrinidae, Percoidei) based on cytochrome B gene sequence variation

Three trophic categories exist within emperor fishes, genus Lethrinus, relating to body form and dentition type. One group contains low-bodied, high speed, stalking predators with conical teeth. Another group comprises high-bodied, slow speed carnivores with molariform teeth capable of crushing hard-shelled benthic prey. A third group is also high-bodied but with conical teeth feeding mostly on small or soft-shelled benthic prey. Inferring the evolution of these trophic types within Lethrinus using morphology is problematic since these characters are typically correlated with feeding mode and are potentially homoplasious. We use mitochondrial DNA sequences, to independently determine a phylogenetic hypothesis for Lethrinus, which are not dependent on morphological characters relating to trophic categories. We analyzed complete cytochrome b gene sequences (1140 bp) for 20 species of Lethrinus, representing the three trophic types, and for 13 outgroup species, including four other representatives of the Lethrinidae. A monophyletic Lethrinidae did not resolve, but the monophyly of Lethrinus is well supported. In addition, two major clades within Lethrinus are well supported. One of these clades exclusively contains low-bodied species with conical teeth while the other clade only comprises the high-bodied species with molariform teeth. A high-bodied species with conical teeth, Lethrinus miniatus, appears most ancestral and sister to all other Lethrinus species. We hypothesize that this generalist trophic type was the evolutionary precursor to both of the other primary trophic types.     

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⁻⁴.     

Molecular Data from the Chloroplast rpoC1 Gene Suggest a Deep and Distinct Dichotomy of Contemporary Spermatophytes into Two Monophyla: Gymnosperms (Including Gnetales) and Angiosperms

Partial sequences of the rpoC1 gene from two species of angiosperms and three species of gymnosperms (8330 base pairs) were determined and compared. The data obtained support the hypothesis that angiosperms and gymnosperms are monophyletic and none of the recent groups of the latter is sister to angiosperms. Received: 20 November 1998 / Accepted: 26 April 1999     

The Sperm Nuclear Basic Proteins (SNBPs) of the Sponge Neofibularia nolitangere: Implications for the Molecular Evolution of SNBPs

We have isolated and characterized for the first time, the SNBPs from an organism (Neofibularia nolitangere) of the phylum Porifera (Sponges). We have shown that these proteins consist of histones which, as expected, exhibit an amino acid composition very similar to that of other eukaryotic histones. The finding of histones in the sperm of these primitive organisms provides support to the notion that histones (SNBPs of the histone, H, type) were the proteins present at the onset of SNBP evolution. In contrast, a discrete number of alternative SNBP types (protamine-like, PL; protamine, P, types) seem to have appeared later on in the course of evolution and are found in both protostomes and deuterostomes, most likely as a result of processes of parallel evolution. Received: 5 March 1997 / Accepted: 6 March 1997     

Rapid Diversification of Marine Picophytoplankton with Dissimilar Light-Harvesting Structures Inferred from Sequences of Prochlorococcus and Synechococcus (Cyanobacteria)

Cultured isolates of the unicellular planktonic cyanobacteria Prochlorococcus and marine Synechococcus belong to a single marine picophytoplankton clade. Within this clade, two deeply branching lineages of Prochlorococcus, two lineages of marine A Synechococcus and one lineage of marine B Synechococcus exhibit closely spaced divergence points with low bootstrap support. This pattern is consistent with a near-simultaneous diversification of marine lineages with divinyl chlorophyll b and phycobilisomes as photosynthetic antennae. Inferences from 16S ribosomal RNA sequences including data for 18 marine picophytoplankton clade members were congruent with results of psbB and petB and D sequence analyses focusing on five strains of Prochlorococcus and one strain of marine A Synechococcus. Third codon position and intergenic region nucleotide frequencies vary widely among members of the marine picophytoplankton group, suggesting that substitution biases differ among the lineages. Nonetheless, standard phylogenetic methods and newer algorithms insensitive to such biases did not recover different branching patterns within the group, and failed to cluster Prochlorococcus with chloroplasts or other chlorophyll b-containing prokaryotes. Prochlorococcus isolated from surface waters of stratified, oligotrophic ocean provinces predominate in a lineage exhibiting low G + C nucleotide frequencies at highly variable positions. Received: 18 January 1997 / Accepted: 18 May 1997