Ascidian and Amphioxus Adh Genes Correlate Functional and Molecular Features of the ADH Family Expansion During Vertebrate Evolution

The alcohol dehydrogenase (ADH) family has evolved into at least eight ADH classes during vertebrate evolution. We have characterized three prevertebrate forms of the parent enzyme of this family, including one from an urochordate (Ciona intestinalis) and two from cephalochordates (Branchiostoma floridae and Branchiostoma lanceolatum). An evolutionary analysis of the family was performed gathering data from protein and gene structures, exon–intron distribution, and functional features through chordate lines. Our data strongly support that the ADH family expansion occurred 500 million years ago, after the cephalochordate/vertebrate split, probably in the gnathostome subphylum line of the vertebrates. Evolutionary rates differ between the ancestral, ADH3 (glutathione-dependent formaldehyde dehydrogenase), and the emerging forms, including the classical alcohol dehydrogenase, ADH1, which has an evolutionary rate 3.6-fold that of the ADH3 form. Phylogenetic analysis and chromosomal mapping of the vertebrate Adh gene cluster suggest that family expansion took place by tandem duplications, probably concurrent with the extensive isoform burst observed before the fish/tetrapode split, rather than through the large-scale genome duplications also postulated in early vertebrate evolution. The absence of multifunctionality in lower chordate ADHs and the structures compared argue in favor of the acquisition of new functions in vertebrate ADH classes. Finally, comparison between B. floridae and B. lanceolatum Adhs provides the first estimate for a cephalochordate speciation, 190 million years ago, probably concomitant with the beginning of the drifting of major land masses from the Pangea. Received: 10 April 2001 / Accepted: 23 May 2001     

Deuterostomic Actin Genes and the Definition of the Chordates: cDNA Cloning and Gene Organization for Cephalochordates and Hemichordates

The evolutionary relationship of muscle and nonmuscle actin isoforms in deuterostomia was studied by the isolation and characterization of two actin genes from the cephalochordate Branchiostoma lanceolatum and two from the hemichordate Saccoglossus kowalevskii The Branchiostoma genes specify a muscle and a nonmuscle actin type, respectively. Together with earlier results on muscle actins from vertebrates and urochordates, a N-terminal sequence signature is defined for chordate muscle actins. These diagnostic amino acid residues separate the chordates from the echinoderms and other metazoa. Although the two Saccoglossus actins characterized so far lack the diagnostic residues, in line with the presumptive phylogenetic position of hemichordates outside the chordates, a definitive conclusion can only be expected once the full complement of actin genes of Saccoglossus is established. Comparison of the intron patterns of the various deuterostomic actin genes shows that intron 330-3, which is present in all vertebrate genes, is conspicuously absent from nonvertebrate genes. The possible origin of this intron is discussed. Received: 4 July 1997 / Accepted: 29 August 1997     

Characterization of the Hydra Lamin and Its Gene: A Molecular Phylogeny of Metazoan Lamins

We report sequences for nuclear lamins from the teleost fish Danio and six invertebrates. These include two cnidarians (Hydra and Tealia), one priapulid, two echinoderms, and the cephalochordate Branchiostoma. Combining these results with earlier data on Drosophila, Caenorhabditis elegans, and various vertebrates, the following conclusions on lamin evolution can be drawn. First, all invertebrate lamins resemble in size the vertebrate B-type lamin. Second, all lamins described previously for amphibia, birds and mammals as well as the first lamin of a fish, characterized here, show a cluster of 7 to 12 acidic residues in the tail domain. Since this acidic cluster is absent from all invertebrate lamins including that of the cephalochordate Branchiostoma, it was acquired with the vertebrate lineage. The larger A-type lamin of differentiated cells must have arisen subsequently by gene duplication and insertion of an extra exon. This extra exon of the vertebrate A-lamins is the only major change in domain organization in metazoan lamin evolution. Third, the three introns of the Hydra and Priapulus genes correspond in position to the last three introns of vertebrate B-type lamin genes. Thus the entirely different gene organization of the C. elegans and Drosophila Dmo genes seems to reflect evolutionary drift, which probably also accounts for the fact that C. elegans has the most diverse lamin sequence. Finally we discuss the possibility that two lamin types, a constitutively expressed one and a developmentally regulated one, arose independently on the arthropod and vertebrate lineages. Received: 4 February 1999 / Accepted: 1 April 1999     

Cloning,Recombinant Expression and Characterization of Wild Type-105-Trp-Calmodulin of the Green Alga Mougeotia scalaris

The single calmodulin gene (CaM) of the green alga Mougeotia scalaris (Hassall) was cloned, sequenced and the CDNA inserted into the prokaryotic expression vector pGEX-2T. The recombinant calmodulin protein (CAM) was expressed as a fusion product together with glutathione S-transferase and isolated on glutathione sepharose. After cleavage and purification, the CaM was characterized by Ca²⁺-dependent shift in SDS-PAGE, by activation of cyclic 3′,5′nucleotide phosphodiesterase (PDE) and sensitivity to the inhibitors trifluoperazine and calmidazolium, with native Mougeotia CaM as control. Using Ca²⁺ buffers in the PDE test, affinity to Ca²⁺ of Mougeotia CaM was found to be diminished fivefold compared to maize or bovine brain CaMs. There was also a 20-fold increase of half maximal activation (K_act) in the PDE test for Mougeotia CaM relative to maize CaM, while the K_act of maize CaM to that of bovine brain CaM was almost the same. The derived amino acid sequences of CaM from Mougeotia and Zea mays revealed three major conservative amino acid exchanges, including unique 105-Trp (Mougeotia) → Leu (maize). In Mougeotia CaM the 105-Trp, including the neighbouring side chains of 92-Phe and 141-Phe, putatively form a hydrophobic ring interaction, as revealed by molecular modelling.     

Evolution of Orthologous Intronless and Intron-Bearing Globin Genes in Two Insect Species

While globin genes ctt-2β and ctt-9.1 in Chironomus thummi thummi each have a single intron, all of the other insect globin genes reported so far are intronless. We analyzed four globin genes linked to the two intron-bearing genes in C. th. thummi. Three have a single intron at the same position as ctt-2β and ctt-9.1; the fourth is intronless and lies between intron bearing genes. Finally, in addition to its intron, one gene (ctt-13RT) was recently interrupted by retrotransposition. Phylogenetic analyses show that the six genes in C. th. thummi share common ancestry with five globin genes in the distantly related species C. tentans, and that a 5-gene ancestral cluster predates the divergence of the two species. One gene in the ancestral cluster gave rise to ctn-ORFB in C. tentans, and duplicated in C. th. thummi to create ctt-11 and ctt-12. From parsimonious calculations of evolutionary distances since speciation, ctt-11, ctt-12, and ctn-ORFB evolved rapidly, while ctn-ORFE in C. tentans evolved slowly compared to other globin genes in the clusters. While these four globins are under selective pressure, we suggest that most chironomid globin genes were not selected for their unique function. Instead, we propose that high gene copy number itself was selected because conditions favored organisms that could synthesize more hemoglobin. High gene copy number selection to produce more of a useful product may be the basis of forming multigene families, all of whose members initially accumulate neutral substitutions while retaining essential function. Maintenance of a large family of globin genes not only ensured high levels of hemoglobin production, but may have facilitated the extensive divergence of chironomids into as many as 5000 species. Received: 31 December 1996 / Accepted: 16 May 1997     

Slow Gating of Gap Junction Channels and Calmodulin

Certain COOH-terminus mutants of connexin32 (Cx32) were previously shown to form channels with unusual transjuctional voltage (V _j ) sensitivity when tested heterotypically in oocytes against Cx32 wild type. Junctional conductance (G _j ) slowly increased by severalfold or decreases to nearly zero with V _j positive or negative, respectively, at mutant side, and V _j positive at mutant side reversed CO₂-induced uncoupling. This suggested that the CO₂-sensitive gate might be a V _j -sensitive slow gate. Based on previous data for calmodulin (CaM) involvement in gap junction function, we have hypothesized that the slow gate could be a CaM-like pore plugging molecule (cork gating model). This study describes a similar behavior in heterotypic channels between Cx32 and each of four new Cx32 mutants modified in cytoplasmic-loop and/or COOH-terminus residues. The mutants are: ML/NN+3R/N, 3R/N, ML/NN and ML/EE; in these mutants, N or E replace M105 and L106, and N replace R215, R219 and R220. This study also reports that inhibition of CaM expression strongly reduces V _j and CO₂ sensitivities of two of the most effective mutants, suggesting a CaM role in slow and chemical gating. Received: 19 April 2000/Revised: 11 August 2000     

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     

Selection and Methionine Accumulation in the Fat Body Protein 2 Gene (FBP2), a Duplicate of the Drosophila Alcohol Dehydrogenase (ADH) Gene

The Drosophila fat body protein 2 gene (Fbp2) is an ancient duplication of the alcohol dehydrogenase gene (Adh) which encodes a protein that differs substantially from ADH in its methionine content. In D. melanogaster, there is one methionine in ADH, while there are 51 (20% of all amino acids) in FBP2. Methionine is involved in 46% of amino acid replacements when Fbp2 DNA sequences are compared between D. melanogaster and D. pseudoobscura. Methionine accumulation does not affect conserved residues of the ADH-ADH^r-FBP2 multigene family. The multigene family has evolved by replacement of mildly hydrophobic amino acids by methionine with no apparent reversion. Its short-term evolution was compared between two Drosophila species, while its long-term evolution was compared between two genera belonging respectively to acalyptrate and calyptrate Diptera, Drosophila and Sarcophaga. The pattern of nucleotide substitution was consistent with an independent accumulation of methionines at the Fbp2 locus in each lineage. Under a steady-state model, the rate of methionine accumulation was constant in the lineage leading to Drosophila, and was twice as fast as that in the calyptrate lineage. Substitution rates were consistent with a slight positive selective advantage for each methionine change in about one-half of amino acid sites in Drosophila. This shows that selection can potentially account for a large proportion of amino acid replacements in the molecular evolution of proteins. Received: 12 December 1994 / Accepted: 15 April 1996     

The single calmodulin gene of the cephalochordate Branchiostoma

The cDNA and gene for calmodulin (CaM) from the cephalochordate Branchiostoma were isolated and characterized. The nucleotide sequence of the Branchiostoma CaM cDNA is about 80% identical to the CaM of Drosophila and Aplysia. However, all nucleotide substitutions are silent, therefore the amino acid sequences of all these CaMs are identical. Branchiostoma and Aplysia CaM genes have the same exon/intron organization. PCR, Northern and genomic Southern analyses showed that Branchiostoma CaM is encoded by a single copy gene, while fish are known to have at least four CaM genes. These results fit the hypothesis that major gene duplication events occurred close to the origin of vertebrates, i.e., after the divergence of the cephalochordate lineage.     

Sequence Diversity in the 16S–23S Intergenic Spacer Region (ISR) of the rRNA Operons in Representatives of the Escherichia coli ECOR Collection

The ribosomal RNA multigene family in Escherichia coli comprises seven rrn operons of similar, but not identical, sequence. Four operons (rrnC, B, G, and E) contain genes in the 16S–23S intergenic spacer region (ISR) for tRNA^Glu-2 and three (rrnA, D, and H) contain genes for tRNA^Ile-1 and tRNA^Ala-1B. To increase our understanding of their molecular evolution, we have determined the ISR sequence of the seven operons in a set of 12 strains from the ECOR collection. Each operon was specifically amplified using polymerase chain reaction primers designed from genes or open reading frames located upstream of the 16S rRNA genes in E. coli K12. With a single exception (ECOR 40), ISRs containing one or two tRNA genes were found at the same respective loci as those of strain K12. Intercistronic heterogeneity already found in K12 was representative of most variation among the strains studied and the location of polymorphic sites was the same. Dispersed nucleotide substitutions were very few but 21 variable sites were found grouped in a stem-loop, although the secondary structure was conserved. Some regions were found in which a stretch of nucleotides was substituted in block by one alternative, apparently unrelated, sequence (as illustrated by the known putative insertion of rsl in K12). Except for substitutions of different sizes and insertions/deletions found in the ISR, the pattern of nucleotide variation is very similar to that found for the 16S rRNA gene in E. coli. Strains K12 and ECOR 40 showed the highest intercistronic heterogeneity. Most strains showed a strong tendency to homogenization. Concerted evolution could explain the notorious conservation of this region that is supposed to have low functional restrictions. Received: 31 July 1997 / Accepted: 17 October 1997     

Updating Carbamoylphosphate Synthase (CPS) Phylogenies: Occurrence and Phylogenetic Identity of Archaeal CPS Genes

Among Bacteria the carA and carB genes encoding the small (CarA) and large (CarB) subunits of carbamoylphosphate synthase (CPS) have been lost in certain symbionts (Haemophylus influenzae) and in most obligate intracellular parasites (Chlamydiae, Spirochaetes, Mycoplasmatales, Rickettsiae) having genome sizes in the 0.7- to 1.1-Mb range. Compared to Bacteria, Archaea exhibit a more varied pattern of CPS gene losses and an unusual propensity to incorporate CPS genes derived from both Bacteria and other Archaea. Schematically they fall into three groups. Group 1 taxa (the crenarchaeon Aeropyrum pernix and the euryarchaea Pyrococcus horikoshi and Pyrococcus abyssii) lack CPS genes altogether. Group 2 taxa (comprising Halobacteriales, Thermoplasmales, Methanococcales, Methanomicrobiales, Archaeoglobales) harbor CPS genes whose encoded CarB and CarA subunit proteins are ostensibly bacterial in origin; that is, they are intermixed with bacterial homologues on a phylogeny of concatenated CarA and CarB sequences and are not distinguishable from bacterial sequences after searching for domain-specific amino acid residue positions. Group 3 taxa (the crenarchaea Pyrobaculum aerophilum, Sulfolobus solfataricus, and Sulfolobus tokodaii and the euryarchaeon Pyrococcus furiosus) harbor CPS genes whose encoded proteins appear to be archaeal: consistent with an archaeal origin, the CarA and CarB sequences in this group possess both unique signatures and signatures affiliating them to Eukarya. Based on the topology of the clade comprising the four Group 3 taxa, we argue that CPS genes of P. furiosus (a euryarchaeon) and those of the crenarchaea P. aerophilum, S. solfataricus, and S. tokodaii are of a single type, resulting from the two genes being laterally transferred from a crenarchaeon to P. furiosus.     

Evolution of MADS-box gene induction by FLO/LFY genes

Some MADS-box genes function as floral homeotic genes. The Arabidopsis LFY gene is a positive regulator of floral homeotic genes, and homologs of the FLO/LFY gene family in other angiosperms and gymnosperms are likely to have a similar function. To investigate the origin of the floral homeotic gene regulatory cascade involving the FLO/LFY gene, FLO/LFY homologs were cloned from a leptosporangiate fern (Ceratopteris richardii), two eusporangiate ferns (Angiopteris lygodiifolia and Botrychium multifidum var. robustum), three fern allies (Psilotum nudum, Equisetum arvense, and Isoetes asiatica), and a moss (Physcomitrella patens). The FLO/LFY gene phylogenetic tree indicates that both duplication and loss of FLO/LFY homologs occurred during the course of vascular plant evolution. The expression patterns of the Ceratopteris LFY genes (CrLFY1 and 2) were assessed. CrLFY1 expression was prominent in tissues including shoot tips and circinate reproductive leaves, but very weak in other tissues examined. Expression of CrLFY2 was also prominent in tissues, including shoot tips and circinate reproductive leaves. These patterns of expression are dissimilar to that of any Ceratopteris MADS-box gene previously reported, suggesting that the induction of MADS-box genes by FLO/LFY is not established at the stage of ferns. Received: 4 January 2001 / Accepted: 28 February 2001     

Sucrose Phosphate Synthase Genes in Plants Belong to Three Different Families

We present phylogenetic analyses to demonstrate that there are three families of sucrose phosphate synthase (SPS) genes present in higher plants. Two data sets were examined, one consisting of full-length proteins and a second larger set that covered a highly conserved region including the 14-3-3 binding region and the UDPGlu active site. Analysis of both datasets showed a well supported separation of known genes into three families, designated A, B, and C. The genomic sequences of Arabidopsis thaliana include a member in each family: two genes on chromosome 5 belong to Family A, one gene on chromosome 1 to Family B, and one gene on chromosome 4 to Family C. Each of three Citrus genes belong to one of the three families. Intron/exon organization of the four Arabidopsis genes differed according to phylogenetic analysis, with members of the same family from different species having similar genomic organization of their SPS genes. The two Family A genes on Arabidopsis chromosome 5 appear to be due to a recent duplication. Analysis of published literature and ESTs indicated that functional differentiation of the families was not obvious, although B family members appear not to be expressed in roots. B family genes were cloned from two Actinidia species and southern analysis indicated the presence of a single gene family, which contrasts to the multiple members of Family A in Actinidia. Only two family C genes have been reported to date. Received: 17 April 2001 / Accepted: 27 August 2001     

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

The larval cuticle protein genes (Lcps) represent a multigene family 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 chromosome. Comparing the DNA sequences from D. miranda and D. melanogaster organization and the gene arrangement of Lcp1–Lcp4 are similar, although the intergene distances vary considerably. The greatest difference between Lcp1 and Lcp2 is due to the occurrence of a pseudogene in D. melanogaster which is not present in D. miranda. Thus the cluster of the four Lcp genes existed already before the separation of the melanogaster and obscura group. Intraspecific homogenizations of different cluster units must have occurred repeatedly between the Lcp1/Lcp2 and Lcp3/Lcp4 sequence types. The most obvious example is exon 2 of the Lcp3 gene in D. miranda, which has been substituted by the corresponding section of the Lcp4 gene rather recently. The homogenization must have occurred before the translocation which generated the neo-Y chromosome. Lcp3 of D. melanogaster has therefore no orthologous partner in D. miranda. Rearrangements in the promoter regions of the D. miranda Lcp genes have generated new, potentially functional CAAT-box motifs. Since three of the Lcp alleles on the neo-Y are not expressed and Lcp3 is expressed only at a reduced level, it is suggestive to speculate that the rearrangements might be involved as cis-regulatory elements in the up-regulation of the X2-chromosomal Lcp alleles, in Drosophila an essential process for dosage compensation. The Lcp genes on the neo-Y chromosome have accumulated more base substitutions than the corresponding alleles on the X2. Received: 27 December 1995 / Accepted: 30 April 1996     

Evolution of Duplicated <Emphasis Type="BoldItalic">reggie</Emphasis> Genes in Zebrafish and Goldfish

Invertebrates, tetrapod vertebrates, and fish might be expected to differ in their number of gene copies, possibly due the occurrence of genome duplication events during animal evolution. Reggie (flotillin) genes code for membrane-associated proteins involved in growth signaling in developing and regenerating axons. Until now, there appeared to be only two reggie genes in fruitflies, mammals, and fish. The aim of this research was to search for additional copies of reggie genes in fishes, since a genome duplication might have increased the gene copy number in this group. We report the presence of up to four distinct reggie genes (two reggie-1 and two reggie-2 genes) in the genomes of zebrafish and goldfish. Phylogenetic analyses show that the zebrafish and goldfish sequence pairs are orthologous, and that the additional copies could have arisen through a genome duplication in a common ancestor of bony fish. The presence of novel reggie mRNAs in fish embryos indicates that the newly discovered gene copies are transcribed and possibly expressed in the developing and regenerating nervous system. The intron/exon boundaries of the new fish genes characterized here correspond with those of human genes, both in location and phase. An evolutionary scenario for the evolution of reggie intron-exon structure, where loss of introns appears to be a distinctive trait in invertebrate reggie genes, is presented. Received: 24 January 2001 / Accepted: 27 July 2001     

Recognition of β-calcineurin by the domains of calmodulin: thermodynamic and structural evidence for distinct roles

Calcineurin (CaN, PP2B, PPP3), a heterodimeric Ca²⁺‐calmodulin‐dependent Ser/Thr phosphatase, regulates swimming in Paramecia, stress responses in yeast, and T‐cell activation and cardiac hypertrophy in humans. Calcium binding to CaN_B (the regulatory subunit) triggers conformational change in CaN_A (the catalytic subunit). Two isoforms of CaN_A (α, β) are both abundant in brain and heart and activated by calcium‐saturated calmodulin (CaM). The individual contribution of each domain of CaM to regulation of calcineurin is not known. Hydrodynamic analyses of (Ca²⁺)₄‐CaM_1–148 bound to βCaNp, a peptide representing its CaM‐binding domain, indicated a 1:1 stoichiometry. βCaNp binding to CaM increased the affinity of calcium for the N‐ and C‐domains equally, thus preserving intrinsic domain differences, and the preference of calcium for sites III and IV. The equilibrium constants for individual calcium‐saturated CaM domains dissociating from βCaNp were ～1 μM. A limiting K_d ≤ 1 nM was measured directly for full‐length CaM, while thermodynamic linkage analysis indicated that it was approximately 1 pM. βCaNp binding to ¹⁵N‐(Ca²⁺)₄‐CaM_1–148 monitored by ¹⁵N/¹HN HSQC NMR showed that association perturbed the N‐domain of CaM more than its C‐domain. NMR resonance assignments of CaM and βCaNp, and interpretation of intermolecular NOEs observed in the ¹³C‐edited and ¹²C‐¹⁴N‐filtered 3D NOESY spectrum indicated anti‐parallel binding. The sole aromatic residue (Phe) located near the βCaNp C‐terminus was in close contact with several residues of the N‐domain of CaM outside the hydrophobic cleft. These structural and thermodynamic properties would permit the domains of CaM to have distinct physiological roles in regulating activation of βCaN. Proteins 2011. © 2010 Wiley‐Liss, Inc.     

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     

The Quality of merC, a Module of the mer Mosaic

We examined a region of high variability in the mosaic mercury resistance (mer) operon of natural bacterial isolates from the primate intestinal microbiota. The region between the merP and merA genes of nine mer loci was sequenced and either the merC, the merF, or no gene was present. Two novel merC genes were identified. Overall nucleotide diversity, π (per 100 sites), of the merC gene was greater (49.63) than adjacent merP (35.82) and merA (32.58) genes. However, the consequences of this variability for the predicted structure of the MerC protein are limited and putative functional elements (metal-binding ligands and transmembrane domains) are strongly conserved. Comparison of codon usage of the merTP, merC, and merA genes suggests that several merC genes are not coeval with their flanking sequences. Although evidence of homologous recombination within the very variable merC genes is not apparent, the flanking regions have higher homologies than merC, and recombination appears to be driving their overall sequence identities higher. The synonymous codon usage bias (EN_C) values suggest greater variability in expression of the merC gene than in flanking genes in six different bacterial hosts. We propose a model for the evolution of MerC as a host-dependent, adventitious module of the mer operon. Received: 2 June 2000 / Accepted: 23 October 2000     

The Role of Yeast VDAC Genes on the Permeability of the Mitochondrial Outer Membrane

In addition to the POR1 gene, which encodes the well-characterized voltage dependent anion-selective channel (YVDAC1) of the mitochondrial outer membrane, the yeast Saccharomyces cerevisiae contains a second gene (POR2) encoding a protein (YVDAC2) with 50% sequence identity to YVDAC1. Mitochondria isolated from yeast cells deleted for the POR1 gene (Δpor1) had a profoundly reduced outer membrane permeability as measured by the ability of an intermembrane space dehydrogenase to oxidize exogenously added NADH. Mitochondria missing either YVDAC1 or both YVDAC1 and YVDAC2 showed a 2-fold increase in the rate of NADH oxidation when the outer membrane was deliberately damaged. Mitochondria from parental cells showed only a 10% increase indicating that the outer membrane is highly permeable to NADH. In the absence of YVDAC1, we calculate that the outer membrane permeability to NADH is reduced 20-fold. The low NADH permeability in the presence of YVDAC2 was not due to the low levels of YVDAC2 expression as mitochondria from cells expressing levels of YVDAC2 comparable to those of YVDAC1 in parental cells showed no substantial increase in NADH permeability, indicating a minimal role of YVDAC2 in this permeability. The residual permeability may be due to other pathways because cells missing both genes can still grow on nonfermentable carbon sources. However, YVDAC1 is clearly the major pathway for NADH flux through the outer membrane in these mitochondria. Received: 23 May 1997/Revised: 3 October 1997