The Geographic Distribution of Monoamine Oxidase Haplotypes Supports a Bottleneck During the Dispersion of Modern Humans from Africa

Every genetic locus mingles the information about the evolutionary history of the human species with the history of its own evolution. Therefore, to address the question of the origin of humans from a genetic point of view, evolutionary histories from many genetic loci have to be gathered and compared. We have studied two genes residing on the X chromosome encoding monoamine oxidases A and B (MAOA and MAOB). Both genes have been suggested to play a role in psychiatric and/or behavioral traits. To search for DNA variants of the MAO genes, the sequences of exonic and flanking intronic regions of these two genes were determined in a group of Swedish males. The sequence analysis revealed several novel polymorphisms in the MAO genes. Haplotypes containing high-frequency MAOA polymorphisms were constructed, and their frequencies were determined in additional samples from Caucasian, Asian, and African populations. We found two common haplotypes with similar frequencies in Caucasian and Asian populations. However, only one of them was also the most frequent haplotype in Africans, while the other haplotype was present in only one Kenyan male. This profound change in haplotype frequencies from Africans to non-Africans supports a possible bottleneck during the dispersion of modern humans from Africa. Received: 19 April 2000 / Accepted: 9 October 2000     

Sequential Loss of Two Neighboring Exons of the Tropoelastin Gene During Primate Evolution

Previous evidence has demonstrated the absence of exons 34 and 35 within the 3′ end of the human tropoelastin (ELN) gene. These exons encode conserved polypeptide domains within tropoelastin and are found in the ELN gene in vertebrate species ranging from chickens to rats to cows. We have analyzed the ELN gene in a variety of primate species to determine whether the absence of exons 34 and 35 in humans either is due to allelic variation within the human population or is a general characteristic of the Primates order. An analysis of the 3′ end of the ELN gene in several nonhuman primates and in 546 chromosomes from humans of varying ethnic background demonstrated a sequential loss of exons 34 and 35 during primate evolution. The loss of exon 35 occurred at least 35–45 million years ago, when Catarrhines diverged from Platyrrhines (New World monkeys). Exon 34 loss, in contrast, occurred only about 6–8 million years ago, when Homo separated from the common ancestor shared with chimpanzees and gorillas. Loss of both exons was probably facilitated by Alu-mediated recombination events and possibly conferred a functional evolutionary advantage in elastic tissue. Received: 6 July 1998 / Accepted: 18 February 1999     

Predominance of one T-cell antigen receptor BV haplotype in African populations

 The human T-cell antigen receptor (TCR) is the counter-receptor for the HLA/peptide complex displayed on the surface of antigen-presenting cells. It confers antigen specificity on T lymphocytes and therefore plays a central role in pathogen recognition and host response. The most frequently used form of the TCR is a heterodimer composed of variable α and β chains. We investigated allele frequencies for four variable-region gene segments of the β chain (2S1, 3S1, 8S3, and 15S1) in 146 Caucasians and 165 Africans. The results reveal significant unexpected differences between the two populations for allele frequencies, phenotypes, genotypes, and haplotypes. Among Caucasians, there are 43 phenotypes, whereas there are 31 among the Africans studied. There are 17 haplotypes in the Caucasian sample but only 10 in Africans. This loss of diversity is largely due to the high frequency of one haplotype in the African sample which represents 65% of the informative chromosomes. At least one copy of this haplotype is present in 90% of informative individuals. As a result, 29% of Africans are homozygous for the common haplotype. Less genetic diversity at TCRBV is unexpected, since Africans usually show greater genetic diversity than other ethnic groups. For example, there are approximately twice as many HLA haplotypes in Africans compared to Caucasians. Homozygosity is also unexpected because it reduces the number of TCR variants available to recognize HLA pathogen-derived peptide complexes. Received: 23 September 1999 / Revised: 2 November 1999     

Mhc Allelic Diversity and Modern Human Origins

Thirty complete coding sequences of human major histocompatibility complex (Mhc) class II DRB alleles, spanning 237 codons, were analyzed for phylogenetic information using distance, parsimony, and likelihood approaches. Allelic genealogies derived from different parts of the coding sequence (exon 2, the 5′ and 3′ ends of exon 2, respectively, and exons 3–6) were compared. Contrary to prior assertions, a rigorous analysis of allelic genealogies in this gene family cannot be used to justify the claim that the lineage leading to modern humans contained on average at least 100,000 individuals. Phylogenetic inferences based upon the exon 2 region of the DRB loci are complicated by selection and recombination, so this part of the gene does not provide a complete and accurate view of allelic relationships. Attempts to reconstruct human history from genetic data must use realistic models which consider the complicating factors of nonequilibrium populations, recombination, and different patterns of selection. Received: 19 February 1997 / Accepted: 12 June 1997     

Distribution of corticotropin-releasing hormone promoter polymorphism in different ethnic groups: evidence for natural selection in human populations

 The regulatory region of the corticotropin-releasing hormone (CRH) is highly conserved across species and plays a crucial role in the response of the organism to stress. Release of CRH initiates a cascade of events leading to the release of cortisol and the regulation of inflammatory and immune events. In this report we describe polymorphisms in the 5′ regulatory region of the CRH gene in humans. We studied the distribution of CRH alleles in three different African populations, in white UK Caucasoids, and in a Chinese population. In the African and UK populations we found three new polymorphisms which cosegregated, resulting in two alleles, A1 and A2. Gene frequencies for A1 and A2 were extremely divergent between the African and the UK populations. The African A1 frequency ranged from 0.27–0.3, while the UK Caucasoid frequency was 0.9. Compound alleles could be assigned by taking into account the previously described biallelic polymorphism at position 225 in the CRH promoter. The A2B1 compound allele is the commonest in contemporary African human populations (allele frequency range 0.44–0.61) and was the only allele observed in a population of chimpanzees from Sierra Leone. Wright's F_ST for the A2B1 allele over the four sampled populations was 0.612, a value exceeded in human populations only by loci which have apparently been subject to natural selection. Taken together, these findings support A2B1 as the ancestral allele and suggest that the CRH genomic region may have been subject to strong disruptive selection throughout human evolution. Received: 29 October 1998 / Revised: 24 March 1999     

Complex and Tandem Repeat Structure of Subtelomeric Regions in the Taiwan Cricket, <Emphasis Type="BoldItalic">Teleogryllus taiwanemma</Emphasis>

Telomeres of most insects are composed of simple (TTAGG) _n repeats that are synthesized by telomerase. However, in some dipteran insects such as Drosophila melanogaster, (TTAGG) _n repeats or telomerase activity has not been detected. Although telomere structure is well documented in Diptera and Lepidoptera, very limited information is available on lower insect groups. To understand general aspects of telomere function and evolution in insects, we endeavored to characterize structures of the telomeric and subtelomeric regions in a lower insect, the Taiwan cricket, Teleogryllus taiwanemma. FISH analysis of this insect's chromosomes demonstrated (TTAGG) _n repeat elements in all distal ends. Just proximal to the telomeric repeats, the highly conserved 9-kb long terminal unit (LTU) sequences are tandemly repeated. These were observed in four of six chromosomes, three autosomal ends, and one X-chromosomal end. LTU sequences represent about 0.2% of the T. taiwanemma genome. Each LTU contains a core (TTAGG)₈-like sequence (TRLS) and five types of conserved sequences—ST (short telomere associated), J (joint), X, SR (satellite sequence rich), and Y—which vary in length from about 150 bp to 2.7 kb. The LTU sequence is defined as ST–J–TRLS–SR–X–Y–X–Y–X. Most LTU regions may be derived from the ancestral common sequence, which is observed in ST regions six times and at many other LTU sites. We could not find the LTU-like sequence in three other crickets including the closest species, T. emma, suggesting that the LTU in T. taiwanemma has been rapidly amplified in subtelomeric regions through recent evolutional events. It is also suggested that the highly conserved structure of the LTU is maintained by recombination and may contribute to telomere elongation, as seen in dipteran insects. Received: 6 August 2001/Accepted: 10 October 2001     

Rates of Nucleotide Substitution in Angiosperm Mitochondrial DNA Sequences and Dates of Divergence Between Brassica and Other Angiosperm Lineages

We obtained 16 nucleotide sequences (∼1400 bp each) of the first intron of the mitochondrial (mt) gene for NADH subunit 4 (nad4) from 10 species of Brassicaceae. Using these new sequences and five published sequences from GenBank, we constructed a phylogenetic tree of the Brassicaceae species under study and showed that the rate of nucleotide substitution in the first intron of nad4 is very low, about 0.16–0.23 × 10⁻⁹ substitution per site per year, which is about half of the silent rate in exons of nad4. The ratios of substitution rates in this intron, ITS, and IGS are approximately 1:23:73, where ITS is the nuclear intergenic spacer between 18S and 25S rRNA genes and IGS is the intergenic spacer of 5S rRNA genes. A segment (335 bp) in the first intron of nad4 in Brassicaceae species that is absent in wheat was considered as a nonfunctional sequence and used to estimate the neutral rate (the rate of mutation) in mtDNA to be 0.5–0.7 × 10⁻⁹ substitution per site per year, which is about three times higher than the substitution rate in the rest of the first intron of nad4. We estimated that the dates of divergence are 170–235 million years (Myr) for the monocot–dicot split, 112–156 Myr for the Brassicaceae–Lettuce split, 14.5–20.4 Myr for the Brassica–Arabidopsis split, and 14.5–20.4 Myr for the Arabidopsis–Arabideae split. Received: 14 July 1998 / Accepted: 1 October 1998     

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     

Evaluation of Sequence Variation and Selection in the Bindin Locus of the Red Sea Urchin, Strongylocentrotus franciscanus

Recent evidence suggests that gamete recognition proteins may be subjected to directed evolutionary pressure that enhances sequence variability. We evaluated whether diversity enhancing selection is operating on a marine invertebrate fertilization protein by examining the intraspecific DNA sequence variation of a 273-base pair region located at the 5′ end of the sperm bindin locus in 134 adult red sea urchins (Strongylocentrotus franciscanus). Bindin is a sperm recognition protein that mediates species-specific gamete interactions in sea urchins. The region of the bindin locus examined was found to be polymorphic with 14 alleles. Mean pairwise comparison of the 14 alleles indicates moderate sequence diversity (p-distance = 1.06). No evidence of diversity enhancing selection was found. It was not possible to reject the null hypothesis that the sequence variation observed in S. franciscanus bindin is a result of neutral evolution. Statistical evaluation of expected proportions of replacement and silent nucleotide substitutions, observed versus expected proportions of radical replacement substitutions, and conformance to the McDonald and Kreitman test of neutral evolution all indicate that random mutation followed by genetic drift created the polymorphisms observed in bindin. Observed frequencies were also highly similar to results expected for a neutrally evolving locus, suggesting that the polymorphism observed in the 5′ region of S. franciscanus bindin is a result of neutral evolution. Received: 19 June 1998 / Accepted: 2 August 2000     

Phylogeographic Patterns and Evolution of the Mitochondrial DNA Control Region in Two Neotropical Cats (Mammalia, Felidae)

The ocelot (Leopardus pardalis) and margay (L. wiedii) are sister-species of Neotropical cats which evolved from a lineage that migrated into South America during the formation of the Panamanian land bridge 3–5 million years ago. Patterns of population genetic divergence of each species were studied by phylogenetic analyses of mitochondrial DNA (mtDNA) control region sequences in individuals sampled across the distribution of these taxa. Abundant genetic diversity and remarkably concordant phylogeographic partitions for both species were observed, identifying parallel geographic regions which likely reflect historical faunal barriers. Inferred aspects of phylogeography, population genetic structure, and demographic history were used to formulate conservation recommendations for these species. In addition, observed patterns of sequence variation provided insight into the molecular evolution of the mtDNA control region in closely related felids. Received: 26 January 1998 / Accepted: 14 May 1998     

Comparison of Substitution Rates in ZFX and ZFY Introns of Sheep and Goat Related Species Supports the Hypothesis of Male-Biased Mutation Rates

There is a growing body of evidence that males serve as the major generators of mutations, due to the larger number of cell divisions involved in sperm compared to egg production. In mammals, this hypothesis (referred to as ``male-driven evolution') has been tested by comparison of nucleotide substitution rates on the X and Y sex chromosomes in a limited number of taxa, predominantly primates and rodents. This study asks whether male-driven evolution is a more general phenomenon among mammals, by comparison of paralogous ZFX and ZFY intron sequences in sheep and goat species (the tribe Caprini). The male-to-female mutation ratio, α_m, was estimated to be between 2.93 (95% CI, 1.51–8.61) and 3.94 (95% CI, 1.25–32.29) when calculated using pairwise distance and branch length, respectively, suggesting that the Caprini are subject to weak, male-driven evolution. Comparison to published values for primates, felids, and rodents implies that there may be some correlation with reproductive life span. However, this is difficult to test with current data because confidence intervals are large and overlapping. Nonindependent evolution of paralogous sequences and/or the presence of selective constraints could lead to inaccurate estimates of α_m. No evidence for gene conversion between the ZFX and the ZFY introns was found, and this suggests that they have evolved independently during the radiation of the Caprini. Finally, there was no apparent evidence that these introns are subject to selective constraints, although low levels of intraspecific polymorphism reduce the power of neutrality tests. Received: 13 February 2001 / Accepted: 23 May 2001     

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     

Recent Selection on Synonymous Codon Usage in Drosophila

Evidence from a variety of sources indicates that selection has influenced synonymous codon usage in Drosophila. It has generally been difficult, however, to distinguish selection that acted in the distant past from ongoing selection. However, under a neutral model, polymorphisms usually reflect more recent mutations than fixed differences between species and may, therefore, be useful for inferring recent selection. If the ancestral state is preferred, selection should shift the frequency distribution of derived states/site toward lower values; if the ancestral is unpreferred, selection should increase the number of derived states/site. Polymorphisms were classified as ancestrally preferred or unpreferred for several genes of D. simulans and D. melanogaster. A computer simulation of coalescence was employed to derive the expected frequency distributions of derived states/site under various modifications of the Wright–Fisher neutral model, and distributions of test statistics (t and Mann–Whitney U) were derived by appropriate sampling. One-tailed tests were applied to transformed frequency data to assess whether the two frequency distributions deviated from neutral expectations in the direction predicted by selection on codon usage. Several genes from D. simulans appear to be subject to recent selection on synonymous codons, including one gene with low codon bias, esterase-6. Selection may also be acting in D. melanogaster. Received: 15 April 1998 / Accepted: 13 May 1999     

Mineral Catalysis of a Potentially Prebiotic Aldol Condensation

Minerals may have played a significant role in chemical evolution. In the course of investigating the chemistry of phosphonoacetaldehyde (PAL), an analogue of glycolaldehyde phosphate, we have observed a striking case of catalysis by the layered hydroxide mineral hydrotalcite ([Mg₂Al(OH)₆][Cl.nH₂O]). In neutral or moderately basic aqueous solutions, PAL is unreactive even at a concentration of 0.1 M. In the presence of a large excess of NaOH (2 M), the compound undergoes aldol condensation to produce a dimer containing a C3–C4 double-bond. In dilute neutral solutions and in the presence of the mineral, however, condensation takes place rapidly, to produce a dimer which is almost exclusively the C2–C3 unsaturated product. Received: 11 February 1998 / Accepted: 12 May 1998     

Sequence of PRAT Satellite DNA ``Frozen' in Some Coleopteran Species

The intriguing diversity of highly abundant satellite repeats found even among closely related species can result from processes leading to dramatic changes in copy number of a particular sequence in the genome and not from rapid accumulation of mutations. To test this hypothesis, we investigated the distribution of the PRAT satellite DNA family, a highly abundant major satellite in the coleopteran species Palorus ratzeburgii, in eight species belonging to the related genera (Tribolium, Tenebrio, Latheticus), the subfamily (Pimeliinae), and the family (Chrysomelidae). Dot blot analysis and PCR assay followed by Southern hybridization revealed that the PRAT satellite, in the form of low-copy number repeats, was present in all tested species. The PRAT satellite detected in the species Pimelia elevata has been sequenced, and compared with previously cloned PRAT monomers from Palorus ratzeburgii and Palorus subdepressus. Although the two Palorus species diverged at least 7 Myr ago, and the subfamily Pimeliinae separated from the genus Palorus 50–60 Myr ago, all PRAT clones exhibit high mutual homology, with average variability relative to the common consensus sequence of 1.3%. The presence of ancestral mutations found in PRAT clones from all three species as well as the absence of species diagnostic mutations illustrate extremely slow sequence evolution. This unexpectedly high conservation of PRAT satellite DNA sequence might be induced by a small bias of turnover mechanisms favoring the ancestral sequence in the process of molecular drive.     

The Nonrandom Location of Synonymous Codons Suggests That Reading Frame-Independent Forces Have Patterned Codon Preferences

Biased codon usage is common in eukaryotic and prokaryotic genes. Evidence from Escherichia, Saccharomyces, and Drosophila indicates that it favors translational efficiency and accuracy. However, to date no functional advantages have been identified in the codon–anticodon interactions involving the most frequently used (preferred) codons. Here we present evidence that forces not related to the individual codon–anticodon interaction may be involved in determining which synonymous codons are preferred or avoided. We show that the ``off-frame' trinucleotide motif preferences inferrable from Drosophila coding regions are often in the same direction as Drosophila's ``in-frame' codon preferences, i.e., its codon usage. The off-frame preferences were inferred from the nonrandomness of the location of confamilial synonymous codons along coding regions—a pattern often described as a context dependence of nucleotide choice at synonymous positions or as codon-pair bias. We relied on randomizations of the location of confamilial codons that do not alter, and cannot be influenced by, the encoded amino acid sequences, codon usage, or base composition of the genes examined. The statistically significant congruency of in-frame and off-frame trinucleotide preferences suggests that the same kind of reading-frame-independent force(s) may also influence synonymous codon choice. These forces may have produced biases in codon usage that then led to the evolution of the translational advantages of these motifs as preferred codons. Under this scenario, tRNA pool size differences between preferred and nonpreferred codons initially were evolved to track the default overrepresentation of codons with preferred motifs. The motif preference hypothesis can explain the structuring of codon preferences and the similarities in the codon usages of distantly related organisms. Received: 10 November 1998 / Accepted: 23 February 1999     

Evaluating a Putative Bottleneck in a Population of Bowhead Whales from Patterns of Microsatellite Diversity and Genetic Disequilibria

A size-selected Balaena mysticetus genomic library was screened for clones containing simple sequence repeat, or microsatellite, loci. A total of 11 novel loci was identified. These loci were combined with a set of 9 published loci, for a total of 20 markers, and were scored across a sample of 108 bowhead whales from the Bering–Chukchi–Beaufort Seas population of bowhead whales. Genetic variability was measured in terms of polymorphism information content values and unbiased heterozygosity. From the latter, estimates of long-term effective population size were obtained. In addition, gametic phase disequilibrium among loci was investigated. Moderate to high levels of polymorphism were found overall, and the long-term effective size estimates were large relative to total population size. Tests of heterozygosity excess (Cornuet and Luikart 1996) and allele frequency distribution (Luikart et al. 1998) indicated that the possibility of a recent genetic bottleneck in the Bering–Chukchi–Beaufort Seas population of bowhead whales is highly unlikely. However, the fact that five loci displayed a statistically significant heterozygote deficiency remains to be explained. Received: 3 November 1998 / Accepted: 28 April 1999     

Interelement Selection in the Regulatory Region of the copia Retrotransposon

We report the results of an analysis of naturally occurring cis-regulatory variation within and between two families of the copia Drosophila long terminal repeat (LTR) retrotransposon. The copia 5′ LTR and adjacent untranslated leader region (ULR) consists of a number of well-characterized sequence motifs which play a role in regulating expression of the element. In order to understand the evolutionary forces which may be responsible for generating and maintaining copia regulatory sequence variation, we have quantified levels of naturally occurring copia LTR-ULR nucleotide variation and subjected the data to a series of tests of neutrality. Our analysis indicates that the copia LTR-ULR has been subject to negative purifying selection within families and positive adaptive selection between families. We discuss these findings with respect to the regulatory evolution of retrotransposons and the phenomenon of interelement selection. Received: 5 February 1998 / Accepted: 14 May 1998     

Molecular Evolution of Cytochrome b of Subterranean Mole Rats, Spalax ehrenbergi Superspecies, in Israel

We describe the molecular evolution of cytochrome b of blind subterranean mole rats. We examined 12 individuals for nucleotide differences in the region of 402 base pairs of mitochondrial cytochrome b. Each individual represents a different population from the entire ecological and speciational range of the four chromosomal species in Israel (2n= 52, 54, 58, and 60) belonging to the Spalax ehrenbergi superspecies. Our results indicate the following. (i) There are seven first-position transitional differences, compared to 34 variable third positions, with no observed second-position substitutions. (ii) A maximum of four amino acids differences occurs across the range. (iii) Within-species diversity increases southward. Only 1 autoapomorphic substitution characterizes either 2n= 52 or 2n= 54, but 6–11 substitutions characterize 2n= 58, and 9–13 substitutions characterize 2n= 60. (iv) Both parsimony and maximum-likelihood trees suggest two monophyletic groups: (a) 2n= 52 and 54, and (b) 2n= 58 and 60, as identified earlier by other protein and DNA markers. (v) Mitochondrial cytochrome b heterogeneity is significantly correlated with climatic factors (rainfall) and biotic factors (body size and allozymes). We hypothesize that two selective regimes direct cytochrome b evolution in the S. ehrenbergi superspecies: (i) purifying selection in the flooded, mesic, hypoxic northern range of 2n= 52 and 54 and (ii) diversifying selection in the climatically spatiotemporal, xeric, and variable southern range of 2n= 58 and 60. Thus, the molecular evolution of mitochondrial cytochrome b in S. ehrenbergi is explicable by opposite selective stresses across the range of S. ehrenbergi in Israel, associated with the ecological adaptive radiation of the complex. Received: 23 October 1998 / Accepted: 2 May 1999