Reconstructing the Diversification of α-Esterases: Comparing the Gene Clusters of Drosophila buzzatii and D. melanogaster

A cluster composed of 10 active α-esterase genes and a pseudogene is distributed over 60 kb in the Drosophila melanogaster genome. This paper describes the corresponding cluster in Drosophila buzzatii, whose lineage diverged from that of D. melanogaster when the subgenera Drosophila and Sophophora diverged about 50 Mya. With three exceptions we find that the composition of the cluster is conserved in the two lineages. The location of αE1 in D. melanogaster differs from that of its nearest relative in D. buzzatii, and αE4 has duplicated independently in the two lineages. The nature of these differences indicates that a mechanism exists whereby copies of genes can be placed in opposite orientation and nonadjacent positions within a gene cluster, although this does not seem to be a feature of earlier events in the cluster's evolution. The rates of amino acid change are not significantly different between orthologs, but the rates differ sevenfold among paralogs, indicating that very different selective forces are acting on the genes of the cluster. Mapping of sequence differences onto a model of the tertiary structure of the enzymes indicates that motifs contributing to substrate binding and catalysis have changed radically in the αE4s and suggest that this subgroup of α-esterases may be evolving into a substantially different functional niche. Received: 4 January 2000 / Accepted: 18 April 2000     

An operon for histidine biosynthesis in Streptomyces coelicolor

Summary In order to demonstrate that a cluster of five his genes (eight cistrons) on the circular chromosome of S. coelicolor is an operon, a constitutive mutant was characterized biochemically, and some aspects of enzyme repression were studied.The specific activities of three enzymes, two of which coded by genes of the his cluster and one specified by a his gene located far from the his cluster, were tested under different repression and derepression conditions and at various times of grwoth in a constitutive his mutant, in two leaky his mutants and in the wild type strains of S. coelicolor.The results of such investigations demonstrate that the constitutive mutant is derepressed exclusively for the synthesis of enzymes coded by genes of the his cluster; moreover only the synthesis of such enzymes takes place in a strictly coordinate way, suggesting that the his cluster behaves as a single unit of expression and regulation.     

Characterization of the ectoine biosynthesis genes of haloalkalotolerant obligate methanotroph “Methylomicrobium alcaliphilum 20Z”

The genes involved in biosynthesis of the major compatible solute ectoine (1,4,5,6-tetrahydro-2-methylpyrimidine carboxylic acid) in halotolerant obligate methanotroph “Methylomicrobium alcaliphilum 20Z” were studied. The complete nucleotide sequences of the structural genes encoding l-aspartokinase (Ask), l-2,4-diaminobutyric acid transaminase (EctB), l-2,4-diaminobutyric acid acetyltransferase (EctA), and l-ectoine synthase (EctC) were defined and shown to be transcribed as a single operon ectABCask. Phylogenetic analysis revealed high sequence identities (34–63%) of the Ect proteins to those from halophilic heterotrophs with the highest amino acid identities being to Vibrio cholerae enzymes. The chromosomal DNA fragment from “M. alcaliphilum 20Z” containing ectABC genes and putative promoter region was expressed in Escherichia coli. Recombinant cells could grow in the presence of 4% NaCl and synthesize ectoine. The data obtained suggested that despite the ectoine biosynthesis pathway being evolutionary well conserved with respect to the genes and enzymes involved, some differences in their organization and regulation could occur in various halophilic bacteria.Dedicated to the 70th birthday of Professor Gerhard Gottschalk who inspired our studies on methylotrophic haloalkaliphiles.     

Molecular Evolution of Nitrogen Fixation: The Evolutionary History of the nifD, nifK, nifE, and nifN Genes

The pairs of nitrogen fixation genes nifDK and nifEN encode for the α and β subunits of nitrogenase and for the two subunits of the NifNE protein complex, involved in the biosynthesis of the FeMo cofactor, respectively. Comparative analysis of the amino acid sequences of the four NifD, NifK, NifE, and NifN in several archaeal and bacterial diazotrophs showed extensive sequence similarity between them, suggesting that their encoding genes constitute a novel paralogous gene family. We propose a two-step model to reconstruct the possible evolutionary history of the four genes. Accordingly, an ancestor gene gave rise, by an in-tandem paralogous duplication event followed by divergence, to an ancestral bicistronic operon; the latter, in turn, underwent a paralogous operon duplication event followed by evolutionary divergence leading to the ancestors of the present-day nifDK and nifEN operons. Both these paralogous duplication events very likely predated the appearance of the last universal common ancestor. The possible role of the ancestral gene and operon in nitrogen fixation is also discussed. Received: 21 June 1999 / Accepted: 1 March 2000     

Reversion of anE. coli strain carrying an IS1-activatedbgl operon under nonselective conditions is predominantly due to deletions within the structural genes

Thebgl operon ofEscherichia coli, which encodes the genes necessary for transport and catabolism of β-glucosides, is silent in wild-type cells and is activated by the transposition of IS elements. The silent form of the operon appears to be the stable state. We isolated Bgl^- revertants of an activated strain after growth under nonselective conditions to understand whether activation of the cryptic operon by IS elements is reversible. Genetic and molecular analyses revealed that a majority of revertants carry deletions of thebgl structural genes, indicating that an irreversible alteration has occurred in the operon. Implications of these results for the evolution and maintenance of cryptic genes are discussed. [Yakkundi A., Moorthy S. and Mahadevan S. 1998 Reversion of anE. coli strain carrying an IS1-activatedbgl operon under nonselective conditions is predominantly due to deletions within the structural genes.J. Genet. 77, 21–26]     

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     

Duplication and divergence of the genes of the α-esterase cluster ofDrosophila melanogaster

The α-esterase cluster of D. melanogaster contains 11 esterase genes dispersed over 60 kb. Embedded in the cluster are two unrelated open reading frames that have sequence similarity with genes encoding ubiquitin-conjugating enzyme and tropomyosin. The esterase amino acid sequences show 37–66% identity with one another and all but one have all the motifs characteristic of functional members of the carboxyl/cholinesterase multigene family. The exception has several frameshift mutations and appears to be a pseudogene. Patterns of amino acid differences among cluster members in relation to generic models of carboxyl/cholinesterase protein structure are broadly similar to those among other carboxyl/cholinesterases sequenced to date. However the α-esterases differ from most other members of the family in: their lack of a signal peptide; the lack of conservation in cysteines involved in disulfide bridges; and in four indels, two of which occur in or adjacent to regions that align with proposed substrate-binding sites of other carboxyl/cholinesterases. Phylogenetic analyses clearly identify three simple gene duplication events within the cluster. The most recent event involved the pseudogene which is located in an intron of another esterase gene. However, relative rate tests suggest that the pseudogene remained functional after the duplication event and has become inactive relatively recently. The distribution of indels also suggests a deeper node in the gene phylogeny that separates six genes at the two ends of the cluster from a block of five in the middle. Received: 18 January 1996 / Accepted: 12 March 1996     

Phylogenetic Depth of the Bacterial Genera Aquifex and Thermotoga Inferred from Analysis of Ribosomal Protein, Elongation Factor, and RNA Polymerase Subunit Sequences

The phylogenetic placement of the Aquifex and Thermotoga lineages has been inferred from (i) the concatenated ribosomal proteins S10, L3, L4, L23, L2, S19, L22, and S3 encoded in the S10 operon (833 aa positions); (ii) the joint sequences of the elongation factors Tu(1α) and G(2) coded by the str operon tuf and fus genes (733 aa positions); and (iii) the joint RNA polymerase β- and β′-type subunits encoded in the rpoBC operon (1130 aa positions). Phylogenies of r-protein and EF sequences support with moderate (r-proteins) to high statistical confidence (EFs) the placement of the two hyperthermophiles at the base of the bacterial clade in agreement with phylogenies of rRNA sequences. In the more robust EF-based phylogenies, the branching of Aquifex and Thermotoga below the successive bacterial lineages is given at bootstrap proportions of 82% (maximum likelihood; ML) and 85% (maximum parsimony; MP), in contrast to the trees inferred from the separate EF-Tu(1α) and EF-G(2) data sets, which lack both resolution and statistical robustness. In the EF analysis MP outperforms ML in discriminating (at the 0.05 level) trees having A. pyrophilus and T. maritima as the most basal lineages from competing alternatives that have (i) mesophiles, or the Thermus genus, as the deepest bacterial radiation and (ii) a monophyletic A. pyrophilus–T. maritima cluster situated at the base of the bacterial clade. RNAP-based phylogenies are equivocal with respect to the Aquifex and Thermotoga placements. The two hyperthermophiles fall basal to all other bacterial phyla when potential artifacts contributed by the compositionally biased and fast-evolving Mycoplasma genitalium and Mycoplasma pneumoniae sequences are eschewed. However, the branching order of the phyla is tenuously supported in ML trees inferred by the exhaustive search method and is unresolved in ML trees inferred by the quartet puzzling algorithm. A rooting of the RNA polymerase-subunit tree at the mycoplasma level seen in both the MP trees and the ML trees reconstructed with suboptimal amino acid substitution models is not supported by the EF-based phylogenies which robustly affiliate mycoplasmas with low-G+C gram-positives and, most probably, reflects a ``long branch attraction' artifact. Received: 22 September 1999 / Accepted: 11 January 2000     

Extraordinary Conservation of Cysteines Among Homologous Chironomus Silk Proteins sp185 and sp220

Aquatic larvae of the midge, Chironomus tentans, synthesize a 185-kDa silk protein (sp185) with the cysteine-containing motif Cys-X-Cys-X-Cys (where X is any residue) every 20–28 residues. We report here the cloning and full-length sequence of cDNAs encoding homologous silk proteins from Chironomus pallidivittatus (sp185) and Chironomus thummi (sp220). Deduced amino acid sequences reveal proteins of nearly identical mass composed of 72 blocks of 20–28 residues, 61% of which can be described by the motif X_5–8-Cys-X₅-(Trp/Phe/Tyr)-X₄-Cys-X-Cys-X-Cys. Spatial arrangement of these residues is preserved more than surrounding sequences. cDNA clones enabled us to map the genes on polytene chromosomes and identify for the first time the homolog of the Camptochironomus Balbiani ring 3 locus in Chironomus thummi. The apparent molecular weight difference between these proteins (185 vs 220 kDa) is not attributable to primary structure and may be due to differential N-linked glycosylation. DNA distances and codon substitutions indicate that the C. tentans and C. pallidivittatus genes are more related to each other than either is to C. thummi; however, substitution rates for the 5′- and 3′-halves of these genes are different. Blockwise sequence comparisons suggest intragenic variation in that some regions evolved slower or faster than the mean and may have been subjected to different selective pressures. Received: 30 August 1996 / Accepted: 6 November 1996     

Heat Shock Protein 70 Family: Multiple Sequence Comparisons, Function, and Evolution

The heat shock protein 70 kDa sequences (HSP70) are of great importance as molecular chaperones in protein folding and transport. They are abundant under conditions of cellular stress. They are highly conserved in all domains of life: Archaea, eubacteria, eukaryotes, and organelles (mitochondria, chloroplasts). A multiple alignment of a large collection of these sequences was obtained employing our symmetric-iterative ITERALIGN program (Brocchieri and Karlin 1998). Assessments of conservation are interpreted in evolutionary terms and with respect to functional implications. Many archaeal sequences (methanogens and halophiles) tend to align best with the Gram-positive sequences. These two groups also miss a signature segment [about 25 amino acids (aa) long] present in all other HSP70 species (Gupta and Golding 1993). We observed a second signature sequence of about 4 aa absent from all eukaryotic homologues, significantly aligned in all prokaryotic sequences. Consensus sequences were developed for eight groups [Archaea, Gram-positive, proteobacterial Gram-negative, singular bacteria, mitochondria, plastids, eukaryotic endoplasmic reticulum (ER) isoforms, eukaryotic cytoplasmic isoforms]. All group consensus comparisons tend to summarize better the alignments than do the individual sequence comparisons. The global individual consensus ``matches' 87% with the consensus of consensuses sequence. A functional analysis of the global consensus identifies a (new) highly significant mixed charge cluster proximal to the carboxyl terminus of the sequence highlighting the hypercharge run EEDKKRRER (one-letter aa code used). The individual Archaea and Gram-positive sequences contain a corresponding significant mixed charge cluster in the location of the charge cluster of the consensus sequence. In contrast, the four Gram-negative proteobacterial sequences of the alignment do not have a charge cluster (even at the 5% significance level). All eukaryotic HSP70 sequences have the analogous charge cluster. Strikingly, several of the eukaryotic isoforms show multiple mixed charged clusters. These clusters were interpreted with supporting data related to HSP70 activity in facilitating chaperone, transport, and secretion function. We observed that the consensus contains only a single tryptophan residue and a single conserved cysteine. This is interpreted with respect to the target rule for disaggregating misfolded proteins. The mitochondrial HSP70 connections to bacterial HSP70 are analyzed, suggesting a polyphyletic split of Trypanosoma and Leishmania protist mitochondrial (Mt) homologues separated from Mt-animal/fungal/plant homologues. Moreover, the HSP70 sequences from the amitochondrial Entamoeba histolytica and Trichomonas vaginalis species were analyzed. The E. histolytica HSP70 is most similar to the higher eukaryotic cytoplasmic sequences, with significantly weaker alignments to ER sequences and much diminished matching to all eubacterial, mitochondrial, and chloroplast sequences. This appears to be at variance with the hypothesis that E. histolytica rather recently lost its mitochondrial organelle. T. vaginalis contains two HSP70 sequences, one Mt-like and the second similar to eukaryotic cytoplasmic sequences suggesting two diverse origins. Received: 29 January 1998 / Accepted: 14 May 1998     

Anopheles stephensi Dox-A2 Shares Common Ancestry with Genes from Distant Groups of Eukaryotes Encoding a 26S Proteasome Subunit and Is in a Conserved Gene Cluster

The sequence of a cloned Anopheles stephensi gene showed 72% inferred amino acid identity with Drosophila melanogaster Dox-A2 and 93% with its putative ortholog in Anopheles gambiae. Dox-A2 is the reported but herein disputed structural locus for diphenol oxidase A2. Database searches identified Dox-A2 related gene sequences from 15 non-insect species from diverse groups. Phylogenetic trees based on alignments of inferred protein sequences, DNA, and protein motif searches and protein secondary structure predictions produced results consistent with expectations for genes that are orthologous. The only inconsistency was that the C-terminus appears to be more primitive in the yeasts than in plants. In mammals, plants, and yeast these genes have been shown to code for a non-ATPase subunit of the PA700 (19S) regulatory complex of 26S proteasome. The analyses indicated that the insect genes contain no divergent structural features, which taken within an appraisal of all available data, makes the reported alternative function highly improbable. A plausible additional role, in which the 26S proteasome is implicated in regulation of phenol oxidase, would also apply to at least the mammalian genes. No function has yet been reported for the other included sequences. These were from genome projects and included Caenorhabiditus elegans, Arabidopsis thaliana, Fugu rubripes, and Toxoplasma gondii. A consensus of the results predicts a protein containing exceptionally long stretches of helix with a hydrophilic C-terminus. Phosphorylation site motifs were identified at two conserved positions. Possible SRY and GATA-1 binding motifs were found at conserved positions upstream of the mosquito genes. The location of A. stephensi Dox-A2 was determined by in situ hybridization at 34D on chromosome arm 3R. It is in a conserved gene cluster with respect to the other insects. However, the A. stephensi cluster contains a gene showing significant sequence identity to human and pigeon carnitine acetyltransferase genes, therefore showing divergence with the distal end of the D. melanogaster cluster. Received: 3 July 1998 / Accepted: 22 December 1999     

The complete mitochondrial DNA (mtDNA) of the donkey and mtDNA comparisons among four closely related mammalian species-pairs

The nucleotide sequence of the complete mitochondrial genome of the donkey, Equus asinus, was determined. The length of the molecule is 16,670 bp. The length, however, is not absolute due to pronounced heteroplasmy caused by variable numbers of two types of repetitive motifs in the control region. The sequence of the repeats is (a) 5′-CACACCCA and (b) 5′-TGCGCGCA, respectively. The order of (a) and (b) can be expressed as {n[2(a)+(b)]+m(a)}. In 32 different clones analyzed the number of n and m ranged from 0 to 9 and 1 to 7. The two rRNA genes, the 13 peptide-coding genes, and the 22 tRNA genes of the donkey and the horse, Equus caballus, were compared in detail. Total nucleotide difference outside the control region was 6.9%. Nucleotide difference between peptide-coding genes ranged from 6.4% to 9.4% with a mean of 8.0%. In the inferred protein sequences of the 13 peptide-coding genes the amino acid difference was 0.2–8.8%, and the mean for the 13 concatenated amino acid sequences was 1.9%. In the 22 tRNA genes, the mean difference was 3.5%, and that in the two rRNA genes was 4.1%. The mtDNA differences between the donkey and the horse suggest that the evolutionary separation of the two species occurred ≈9 million years ago. Analyses of differences among the mtDNAs of three other species-pairs, harbor seal/grey seal, fin whale/blue whale, and Homo/common chimpanzee, showed that the relative evolutionary rate of individual peptide-coding genes varies among different species-pairs and modes of comparison. The findings show that the superimposition of sequence data of one lineage for resolving and dating evolutionary divergences of other lineages should be performed with caution unless based on comprehensive data. Received: 15 October 1995 / Accepted: 15 April 1996     

Organization, Structure, and Evolution of the Nonadult Rat β-Globin Gene Cluster

The β-globin gene cluster of Wistar rat was extensively cloned and the embryonic genes were mapped and sequenced. Four overlapping λ Dash recombinant clones cover about 31 kb and contain four nonadult β-globin genes, 5′–ε1–γ1–γ2–ψγ3–3′. The ε1 and γ2 are active genes, since their protein products were detected in the fetal stage of the rat (Iwahara et al., J Biochem 119:360–366, 1996). The γ1 locus might be a pseudogene, since the ATA box in the promoter region is mutated to GTA; however, no other defect is observed. The ψγ3 locus is a truncated pseudogene because a 19-base deletion, which causes a shift of the reading frame, is observed between the second nucleotide of the putative codon 68 and codon 76. A sequence comparison suggests that the ψγ3 might be produced by a gene conversion event of the proto-γ-globin gene set. Possible histories of the evolution of rat nonadult β-globin genes are discussed. Received: 6 August 1998 / Accepted: 12 February 1999     

The sequences of heat shock protein 40 (DnaJ) homologs provide evidence for a close evolutionary relationship between the Deinococcus- Thermus group and cyanobacteria

The genes encoding for heat shock protein 40 (Hsp40 or DnaJ) homologs were cloned and sequenced from the archaebacterium Halobacterium cutirubrum and the eubacterium Deinococcus proteolyticus to add to sequences from the gene banks. These genes were identified downstream of the Hsp70 (or DnaK) genes in genomic fragments spanning this region and, as in other prokaryotic species, Hsp70-Hsp40 genes are likely part of the same operon. The Hsp40 homolog from D. proteolyticus was found to be lacking a central 204 base pair region present in H. cutirubrum that encodes for the four cysteine-rich domains of the repeat consensus sequence CxxCxGxG (where x is any amino acid), present in most Hsp40 homologs. The available sequences from various archaebacteria, eubacteria, and eukaryotes show that the same deletion is also present in the homologs from Thermus aquaticus and two cyanobacteria, but in no other species tested. This unique deletion and the clustering of homologs from the Deinococcus–Thermus group and cyanobacterial species in the Hsp40 phylogenetic trees suggest a close evolutionary relationship between these groups as was also shown recently for Hsp70 sequences (R.S. Gupta et al., J Bacteriol 179:345–357, 1997). Sequence comparisons indicate that the Hsp40 homologs are not as conserved as the Hsp70 sequences. Phylogenetic analysis provides no reliable information concerning evolutionary relationship between prokaryotes and eukaryotes and their usefulness in this regard is limited. However, in phylogenetic trees based on Hsp40 sequences, the two archaebacterial homologs showed a polyphyletic branching within Gram-positive bacteria, similar to that seen with Hsp70 sequences. Received: 30 January 1997 / Accepted: 22 March 1997     

Sub-epitopic dissection of HCV E1<Subscript>315–328</Subscript>HRMAWDMMMNWSPT sequence by similarity analysis

Summary. Our labs are focused on identifying amino acid sequences having the ability to react specifically with the functional binding site of a complementary antibody. Our epitopic definition is based on the analysis of the similarity level of antigenic amino acid sequences to the host proteome. Here, the similarity profile to the human proteome of an HCV E1 immunodominant epitope, i.e. the HCV E1_315–328HRMAWDMMMNWSPT sequence, led to i) characterizing the immunoreactive HCV E1 315–328 region as a sequence endowed with a low level of similarity to human proteins; ii) defining 2 contiguous immunodominant linear determinants respectively located at the NH₂ and COOH terminus of the conserved viral antigenic sequence. This study supports the hypothesis that low sequence similarity to the host’s proteome modulates the pool of epitopic amino acid sequences in a viral antigen, and appears of potential value in defining immunogenic viral peptide sequences to be used in immunotherapeutic approaches for HCV treatment. Authors’ address: D. Kanduc, Department of Biochemistry and Molecular Biology, University of Bari, Via Orabona 4, Bari 70125, Italy