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     

Replacement of nucleosomal histones by histone H1-like proteins during spermiogenesis in Cnidaria: Evolutionary implications

We have analyzed the chromosomal protein composition of the sperm from several species belonging to three different classes (Hydrozoa, Scyphozoa, Anthozoa) of the phylum Cnidaria. In every instance, the sperm nuclear basic proteins (SNBPs) were found to consist of one to two major protein fractions that belong to the histone H1 family, as can be deduced from their amino acid composition and solubility in dilute perchloric acid, and the presence of a trypsin-resistant core. In those species where mature spawned sperm could be obtained, we were able to show that these proteins completely replace the somatic histones from the stem cells that are present at the onset of spermatogenesis. The presence of a highly specialized histone H1 molecule in the sperm of this phylum provides support for the idea that the protamine-like proteins (PL) from higher groups in the phylogenetic tree (and possibly protamines as well) may all have evolved from a primitive histone H1 ancestor.     

Evolution of Chitin-Binding Proteins in Invertebrates

Analysis of a group of invertebrate proteins, including chitinases and peritrophic matrix proteins, reveals the presence of chitin-binding domains that share significant amino acid sequence similarity. The data suggest that these domains evolved from a common ancestor which may be a protein containing a single chitin-binding domain. The duplication and transposition of this chitin-binding domain may have contributed to the functional diversification of chitin-binding proteins. Sequence comparisons indicated that invertebrate and plant chitin binding domains do not share significant amino acid sequence similarity, suggesting that they are not coancestral. However, both the invertebrate and the plant chitin-binding domains are cysteine-rich and have several highly conserved aromatic residues. In plants, cysteines have been elucidated in maintaining protein folding and aromatic amino acids in interacting with saccharides [Wright HT, Sanddrasegaram G, Wright CS (1991) J Mol Evol 33:283–294]. It is likely that these residues perform similar functions in invertebrates. We propose that the invertebrate and the plant chitin-binding domains share similar mechanisms for folding and saccharide binding and that they evolved by convergent evolution. Furthermore, we propose that the disulfide bonds and aromatic residues are hallmarks for saccharide-binding proteins. Received: 2 March 1998 / Accepted: 17 July 1998     

Evolution of the Larval Cuticle Proteins Coded by the Secondary Sex Chromosome Pair: X2 and Neo-Y of Drosophila miranda: II. Comparison at the Amino Acid Sequence Level

The larval cuticle proteins (LCPs) are encoded by a multigene family, Lcp1–4, located at the right arm of the metacentric autosome 2 (2R) in Drosophila melanogaster. Due to a chromosome fusion the Lcp locus of Drosophila miranda is situated on a pair of secondary sex chromosomes, the X2 and neo-Y chromosomes. Comparing the deduced amino acid sequences of the autosomal D. melanogaster loci with the sex-chromosomal loci of D. miranda, we were able to trace the evolution of the Lcp loci with respect to their different chromosomal inheritance. The length of the signal peptide is conserved in all four LCPs, while the size of the mature LCPs varies. Conserved protein motifs became obvious from the alignment, indicating regions of structural and functional importance. Analyzing intra- and interspecific sequence similarities of the Lcp gene families allowed us to reconstruct the phylogeny of the gene cluster. Alignment with cuticle amino acid sequences originating from divergent insect species reveals motifs already present in the primordial insect LCPs. These motifs indicate different levels of constraint acting during the evolution of the LCPs. Received: 27 December 1995 / Accepted: 30 April 1996     

Rapid Rate of Control-Region Evolution in Pacific Butterflyfishes (Chaetodontidae)

Sequence differences in the tRNA-proline (tRNA^pro) end of the mitochondrial control-region of three species of Pacific butterflyfishes accumulated 33–43 times more rapidly than did changes within the mitochondrial cytochrome b gene (cytb). Rapid evolution in this region was accompanied by strong transition/transversion bias and large variation in the probability of a DNA substitution among sites. These substitution constraints placed an absolute ceiling on the magnitude of sequence divergence that could be detected between individuals. This divergence ``ceiling' was reached rapidly and led to a decay in the relative rate of control-region/cytb b evolution. A high rate of evolution in this section of the control-region of butterflyfishes stands in marked contrast to the patterns reported in some other fish lineages. Although the mechanism underlying rate variation remains unclear, all taxa with rapid evolution in the 5′-end of the control-region showed extreme transition biases. By contrast, in taxa with slower control-region evolution, transitions accumulated at nearly the same rate as transversions. More information is needed to understand the relationship between nucleotide bias and the rate of evolution in the 5′-end of the control-region. Despite strong constraints on sequence change, phylogenetic information was preserved in the group of recently differentiated species and supported the clustering of sequences into three major mtDNA groupings. Within these groups, very similar control-region sequences were widely distributed across the Pacific Ocean and were shared between recognized species, indicating a lack of mitochondrial sequence monophyly among species. Received: 30 June 1996 / Accepted: 15 May 1997     

The Evolution of the Conserved ATPase Domain (CAD): Reconstructing the History of an Ancient Protein Module

The AAA proteins (ATPases Associated with a variety of cellular Activities) are found in eubacterial, archaebacterial, and eukaryotic species and participate in a large number of cellular processes, including protein degradation, vesicle fusion, cell cycle control, and cellular secretory processes. The AAA proteins are characterized by the presence of a 230 to 250-amino acid ATPase domain referred to as the Conserved ATPase Domain or CAD. Phylogenetic analysis of 133 CAD sequences from 38 species reveal that AAA CADs are organized into discrete groups that are related not only in structure but in cellular function. Evolutionary analyses also indicate that the CAD was present in the last common ancestor of eubacteria, archaebacteria, and eukaryotes. The eubacterial CADs are found in metalloproteases, while CAD-containing proteins in the archaebacterial and eukaryotic lineages appear to have diversified by a series of gene duplication events that lead to the establishment of different functional AAA proteins, including proteasomal regulatory, NSF/Sec, and Pas proteins. The phylogeny of the CADs provides the basis for establishing the patterns of evolutionary change that characterize the AAA proteins. Received: 28 January 1997 / Accepted: 8 May 1997     

On the Evolution of Dopa decarboxylase (Ddc) and Drosophila Systematics

We have sequenced most of the coding region of the gene Dopa decarboxylase (Ddc) in 24 fruitfly species. The Ddc gene is quite informative about Drosophila phylogeny. Several outstanding issues in Drosophila phylogeny are resolved by analysis of the Ddc sequences alone or in combination with three other genes, Sod, Adh, and Gpdh. The three species groups, melanogaster, obscura, and willistoni, are each monophyletic and all three combined form a monophyletic group, which corresponds to the subgenus Sophophora. The Sophophora subgenus is the sister group to all other Drosophila subgenera (including some named genera, previously considered outside the Drosophila genus, namely, Scaptomyza and Zaprionus, which are therefore downgraded to the category of subgenus). The Hawaiian Drosophila and Scaptomyza are a monophyletic group, which is the sister clade to the virilis and repleta groups of the subgenus Drosophila. The subgenus Drosophila appears to be paraphyletic, although this is not definitely resolved. The two genera Scaptodrosophila and Chymomyza are older than the genus Drosophila. The data favor the hypothesis that Chymomyza is older than Scaptodrosophila, although this issue is not definitely resolved. Molecular evolution is erratic. The rates of nucleotide substitution in 3rd codon position relative to positions 1 + 2 vary from one species lineage to another and from gene to gene. Received: 2 June 1998 / Accepted: 3 September 1998     

Getting In or Out: Early Segregation Between Importers and Exporters in the Evolution of ATP-Binding Cassette (ABC) Transporters

ATP-binding cassette (ABC) systems, also called traffic ATPases, are found in eukaryotes and prokaryotes and almost all participate in the transport of a wide variety of molecules. ABC systems are characterized by a highly conserved ATPase module called here the ABC module, involved in coupling transport to ATP hydrolysis. We have used the sequence of one of the first representatives of bacterial ABC transporters, the MalK protein, to collect 250 closely related sequences from a nonredundant protein sequence database. The sequences collected by this objective method are all known or putative ABC transporters. After having eliminated short protein sequences and duplicates, the 197 remaining sequences were subjected to a phylogenetic analysis based on a mutational similarity matrix. An unrooted tree for these modules was found to display two major branches, one grouping all collected uptake systems and the other all collected export systems. This remarkable disposition strongly suggests that the divergence between these two functionally different types of ABC systems occurred once in the history of these systems and probably before the differentiation of prokaryotes and eukaryotes. We discuss the implications of this finding and we propose a model accounting for the generation and the diversification of ABC systems. Received: 23 February 1997 / Accepted: 7 April 1998     

Nuclear and Nucleomorph SSU rDNA Phylogeny in the Cryptophyta and the Evolution of Cryptophyte Diversity

The plastid-bearing members of the Cryptophyta contain two functional eukaryotic genomes of different phylogenetic origin, residing in the nucleus and in the nucleomorph, respectively. These widespread and diverse protists thus offer a unique opportunity to study the coevolution of two different eukaryotic genomes within one group of organisms. In this study, the SSU rRNA genes of both genomes were PCR-amplified with specific primers and phylogenetic analyses were performed on different data sets using different evolutionary models. The results show that the composition of the principal clades obtained from the phylogenetic analyses of both genes was largely congruent, but striking differences in evolutionary rates were observed. These affected the topologies of the nuclear and nucleomorph phylogenies differently, resulting in long-branch attraction artifacts when simple evolutionary models were applied. Deletion of long-branch taxa stabilized the internal branching order in both phylogenies and resulted in a completely resolved topology in the nucleomorph phylogeny. A comparison of the tree topologies derived from SSU rDNA sequences with characters previously used in cryptophyte systematics revealed that the biliprotein type was congruent, but the type of inner periplast component incongruent, with the molecular trees. The latter is indicative of a hidden cellular dimorphism (cells with two periplast types present in a single clonal strain) of presumably widespread occurrence throughout cryptophyte diversity, which, in consequence, has far-reaching implications for cryptophyte systematics as it is practiced today.     

The Family of Major Royal Jelly Proteins and Its Evolution

A cDNA encoding a new member of the gene family of major royal jelly proteins (MRJPs) from the honeybee, Apis mellifera, was isolated and sequenced. Royal jelly (RJ) is a secretion of the cephalic glands of nurse bees. The origin and biological function of the protein component (12.5%, w/w) of RJ is unknown. We show that the MRJP gene family encodes a group of closely related proteins that share a common evolutionary origin with the yellow protein of Drosophila melanogaster. Yellow protein functions in cuticle pigmentation in D. melanogaster. The MRJPs appear to have evolved a novel nutritional function in the honeybee. Received: 26 September 1998 / Accepted: 28 February 1999     

Evolution of Histone H4 and H3 Genes in Different Ciliate Lineages

The histones H4 are known as highly conserved proteins. However, in ciliates a high degree of variation was found compared both to other eukaryotes and between the ciliate species. To date, only H4 histones of species belonging to two distantly related classes have been investigated. In order to obtain more detailed information on histone H4 variation in ciliates we undertook a comprehensive sequence analysis of PCR-amplified internal H4 fragments from 12 species belonging to seven out of the nine currently recognized ciliate classes. In addition, we used PCR primers to amplify longer fragments of H3 and H4 genes including the intergenic region. The encoded amino acid sequences reveal a high number of differences when compared with those of other eukaryotes and the ciliate species investigated. Furthermore, in some species H4 gene variants were detected, which result in amino acid differences. The greatest number of substitutions and insertions found was in the amino terminal region of the H4 histones. However, all sequences possess a conserved region corresponding to those of all other eukaryotic H4 histones. The histone gene variations were used to reconstruct phylogenetic relationships. The tree from our data matches perfectly with the ribosomal RNA data: The heterotrichs, which were considered as a late branching lineage, diverge at the base of the ciliate tree and groups formerly thought to represent ancestral lineages now appear as highly derived ciliates. Received: 4 April 1997 / Accepted: 1 August 1997     

Evolution of DNA Polymerase Families: Evidences for Multiple Gene Exchange Between Cellular and Viral Proteins

A phylogenetic analysis of the five major families of DNA polymerase is presented. Viral and plasmid sequences are included in this compilation along with cellular enzymes. The classification by Ito and Braithwaite (Ito and Braithwaite 1991) of the A, B, C, D, and X families has been extended to accommodate the ``Y family' of DNA polymerases that are related to the eukaryotic RAD30 and the bacterial UmuC gene products. After analysis, our data suggest that no DNA polymerase family was universally conserved among the three biological domains and no simple evolutionary scenario could explain that observation. Furthermore, viruses and plasmids carry a remarkably diverse set of DNA polymerase genes, suggesting that lateral gene transfer is frequent and includes non-orthologous gene displacements between cells and viruses. The relationships between viral and host genes appear very complex. We propose that the gamma DNA polymerase of the mitochondrion replication apparatus is of phage origin and that this gene replaced the one in the bacterial ancestor. Often there was no obvious relation between the viral and the host DNA polymerase, but an interesting exception concerned the family B enzymes: in which ancient gene exchange can be detected between the viruses and their hosts. Additional evidence for horizontal gene transfers between cells and viruses comes from an analysis of the small damage-inducible DNA polymerases. Taken together, these findings suggest a complex evolutionary history of the DNA replication apparatus that involved significant exchanges between viruses, plasmids, and their hosts.     

The Molecular Evolution of Visual Pigments of Freshwater Crayfishes (Decapoda: Cambaridae)

This study examines the diverse maximum wavelength absorption (λ_max) found in crayfishes (Decapoda: Cambaridae and Parastacidae) and the associated genetic variation in their opsin locus. We measured the wavelength absorption in the photoreceptors of six species that inhabit environments of different light intensities (i.e., burrows, streams, standing waters, and subterranean waters). Our results indicate that there is relatively little variation in λ_max (522–530 nm) among species from different genera and families. The existing variation did not correlate with the habitat differences of the crayfishes studied. We simultaneously sequenced the rhodopsin gene to identify the amino acid replacements that affect shifts in maximum wavelength absorption. We then related these to changes that correlated with shifts in λ_max by reconstructing ancestral character states using a maximum-likelihood approach. Using amino acid sequences obtained from five species (all were 301 amino acids in length), we identified a number of candidates for producing shifts of 4 to 8 nm in λ_max. These amino acid replacements occurred in similar regions to those involved in spectral shifts in vertebrates. Received: 12 March 1997 / Accepted: 3 June 1997     

Sociality and the Rate of rDNA Sequence Evolution in Wasps (Vespidae) and Honeybees (Apis)

Sequence data of mitochondrial 16S ribosomal DNA (mt-rDNA) and nuclear 28S ribosomal DNA (nuc-rDNA) were compared in two honeybee species (Apis mellifera and Apis dorsata) and a selection of 22 wasp species (Vespidae) with different levels of sociality. The averge substitution rates in mt-rDNA and nuc-rDNA were almost-equal in solitary species. In species with larger nests, however, the difference between the nuclear and the mitochondrial substitution rate significantly increased. The average substitution ratio, ψ (nucleotide substitutions in mt-rDNA/nucleotide substitutions in nuc-rDNA) was 1.48 ± 0.12 (SE) among the solitary Eumeninae, 3.70 ± 0.15 among five primitive social Stenogastrinae species, 3.24 ± 0.20 among five Polistinae species, 5.76 ± 0.33 among nine highly eusocial Vespinae, and 12.7 in the two Apis species. The high egg-laying rate and the effective population size skew between the sexes may contribute to the rise of the substitution ratio in the highly eusocial species. Drift and bottleneck effects in the mitochondrial DNA pool during speciation events as well as polyandry may further enhance this phenomenon. Received: 12 January 1998 / Accepted: 28 April 1998     

The Evolution of the MAP Kinase Pathways: Coduplication of Interacting Proteins Leads to New Signaling Cascades

The MAP-kinase pathways are intracellular signaling modules that are likely to exist in all eukaryotes. We provide an evolutionary model for these signaling pathways by focusing on the gene duplications that have occurred since the divergence of animals from yeast. Construction of evolutionary trees with confidence assessed by bootstrap clearly shows that the mammalian JNK and p38 pathways arose from an ancestral hyperosmolarity pathway after the split from yeast and before the split from C. elegans. These coduplications of interacting proteins at the MAPK and MEK levels have since evolved toward substrate specificity, thus giving distinct pathways. Mammalian duplications since the split from C. elegans are often associated with divergent tissue distribution but do not appear to confer detectable substrate specificity. The yeast kinase cascades have undergone similar fundamental functional changes since the split from mammals, with duplications giving rise to central signaling components of the filamentous and hypoosmolarity pathways. Experimentally defined cross-talk between yeast pheromone and hyperosmolarity pathways is mirrored with corresponding cross-talk in mammalian pathways, suggesting the existence of ancient orthologous cross-talk; our analysis of gene duplications at all levels of the cascade is consistent with this model but does not always provide significant bootstrap support. Our data also provide insights at different levels of the cascade where conflicting experimental evidence exists. Received: 2 December 1998 / Accepted: 9 June 1999     

Determination of the evolutionary relationships in Rattus sensu lato (Rodentia : Muridae) using L1 (LINE-1) amplification events

We determined ∼215 bp of DNA sequence from the 3′-untranslated region (UTR) of 240 cloned L1 (LINE-1) elements isolated from 22 species of Rattus sensu lato and Rattus sensu stricto murine rodents. The sequences were sorted into different L1 subfamilies, and oligonucleotides cognate to them were hybridized to genomic DNA of various taxa. From the distribution of the L1 subfamilies in the various species, we inferred the partial phylogeny of Rattus sensu lato. The four Maxomys species comprise a well-defined clade separate from a monophyletic cluster that contains the two Leopoldamys and four Niviventer species. The Niviventer/Leopoldamys clade, in turn, shares a node with the clade that contains Berylmys, Sundamys, Bandicota, and Rattus sensu stricto. The evolutionary relationships that we deduced agree with and significantly extend the phylogeny of Rattus sensu lato established by other molecular criteria. Furthermore, the L1 amplification events scored here produced a unique phylogenetic tree, that is, in no case did a character (a given L1 amplification event) appear on more than one branch. The lack of homoplasy found in this study supports the robustness of L1 amplification events as phylogenetic markers for the study of mammalian evolution. Received: 8 November 1996 / Accepted: 11 April 1997     

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     

Characterization and Evolution of the Mitochondrial DNA Control Region in Hornbills (Bucerotiformes)

We determined the mitochondrial DNA control region sequences of six Bucerotiformes. Hornbills have the typical avian gene order and their control region is similar to other avian control regions in that it is partitioned into three domains: two variable domains that flank a central conserved domain. Two characteristics of the hornbill control region sequence differ from that of other birds. First, domain I is AT rich as opposed to AC rich, and second, the control region is approximately 500 bp longer than that of other birds. Both these deviations from typical avian control region sequence are explainable on the basis of repeat motifs in domain I of the hornbill control region. The repeat motifs probably originated from a duplication of CSB-1 as has been determined in chicken, quail, and snowgoose. Furthermore, the hornbill repeat motifs probably arose before the divergence of hornbills from each other but after the divergence of hornbills from other avian taxa. The mitochondrial control region of hornbills is suitable for both phylogenetic and population studies, with domains I and II probably more suited to population and phylogenetic analyses, respectively.     

Evolution of Microsatellite Alleles in Four Species of Mice (Genus Apodemus)

Microsatellite length variation was investigated at a highly variable microsatellite locus in four species of Apodemus. Information obtained from microsatellite allele sequences was contrasted with allele sizes, which included 18 electromorphs. Additional analysis of a 400-bp unique sequence in the flanking region identified 26 different haplotype sequences or ``true' alleles in the sample. Three molecular mechanisms, namely, (1) addition/deletion of repeats, (2) substitutions and indels in the flanking region, and (3) mutations interrupting the repeat, contributed to the generation of allelic variation. Size homoplasy can be inferred for alleles within populations, from different populations of the same species, and from different species. We propose that microsatellite flanking sequences may be informative markers for investigating mutation processes in microsatellite repeats as well as phylogenetic relationships among alleles, populations, and species. Received: 3 November 1999 / Accepted: 2 May 2000