Analysis of Orthologous Retrovirus-Like Elements in the White-Footed Mouse,Peromyscus leucopus

Three loci in the genome of the white-footed mouse, Peromyscus leucopus, were examined for the presence or absence of orthologous copies of the retrovirus-like element mys using polymerase chain reaction. We examined these loci in 28 mice collected throughout the P. leucopus species range. Mys insertions were present in only one of the individuals examined at the mys-1 and mys-7 loci. Conversely, the mys-6 element was found in several individuals, but the presence of this element was limited to northern latitudes. Because the long terminal repeats (LTRs) of a given element are expected to be identical at the time of retrotransposition into the genome, and to accumulate changes over evolutionary time, within-element LTR sequence comparisons can be used to estimate the relative age of insertions. Within-element LTR differences are greater in mys-6 than in mys-1 or mys-7. The LTRs from orthologous mys-6 elements of six mice were sequenced. The alignment revealed 13 of the 22 differences between the right and left LTRs that were shared by all orthologous mys-6 sites, suggesting that relative to its time of insertion into the genome, mys-6 has only recently spread across the northern part of the species range. Received: 23 January 1996 / Accepted: 24 April 1996     

Isolation and Characterization of Five Actin cDNAs from the Cestode Diphyllobothrium dendriticum: A Phylogenetic Study of the Multigene Family

Five cDNAs (pDidact2–pDidact6), representing different actin genes, were isolated from a Diphyllobothrium dendriticum cDNA library, and the DNA as well as the putative amino acid sequences were determined. The corresponding Didact2 and Didact4 genes code for peptides 376 amino acids long, with molecular weights 41,772 and 41,744 Da, respectively, while the deduced Didact3 protein is 377 amino acids long and weighs 41,912 Da. The pDidact5 and -6 cDNAs lack nucleotides corresponding to three to six amino acids at the amino-terminus. Two of the five cDNAs contain the conventional AATAAA as the putative polyadenylation signal, one has the common variant ATTAAA, whereas the hexanucleotide AATAGA is found 15 and 18 nucleotides, respectively, upstream of the poly(A) site in two of the cDNAs. Phylogenetic studies including 102 actin protein sequences revealed that there are at least four different types of cestode actins. In this study three of these types were found to be expressed in the adult D. dendriticum tapeworm. Structurally the cestode actin groupings differ from each other to an extent seen only among the metazoan actins between the vertebrate muscle and cytoplasmic isoforms. In the phylogenetic trees constructed, cestode actins were seen to map to two different regions, one on the border of the metazoan actins and the other within this group. It is, however, difficult to say whether the cestode actins branched off early in the metazoan evolution or if this position in the phylogenetic tree only reflects upon differences in evolutionary rate. Received: 19 June 1996 / Accepted: 20 August 1996     

How did glycogen structure evolve to satisfy the requirement for rapid mobilization of glucose? A problem of physical constraints in structure building

Optimization of molecular design in cellular metabolism is a necessary condition for guaranteeing a good structure–function relationship. We have studied this feature in the design of glycogen by means of the mathematical model previously presented that describes glycogen structure and its optimization function [Meléndez-Hevia et al. (1993), Biochem J 295: 477–483]. Our results demonstrate that the structure of cellular glycogen is in good agreement with these principles. Because the stored glucose in glycogen must be ready to be used at any phase of its synthesis or degradation, the full optimization of glycogen structure must also imply the optimization of every intermediate stage in its formation. This case can be viewed as a molecular instance of the eye problem, a classical paradigm of natural selection which states that every step in the evolutionary formation of a functional structure must be functional. The glycogen molecule has a highly optimized structure for its metabolic function, but the optimization of the full molecule has meaning and can be understood only by taking into account the optimization of each intermediate stage in its formation. Received: 23 October 1996 / Accepted: 21 April 1997     

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     

Evolution of 28S Ribosomal DNA in Chaetognaths: Duplicate Genes and Molecular Phylogeny

The chaetognaths are an extraordinarily homogeneous phylum of animals at the morphological level, with a bauplan that can be traced back to the Cambrian. Despite the attention of zoologists for over two centuries, there is little agreement on classification within the phylum. We have used a molecular biological approach to investigate the phylogeny of extant chaetognaths. A rapidly evolving expansion segment toward the 5′ end of 28S ribosomal DNA (rDNA) was amplified using the polymerase chain reaction (PCR), cloned, and sequenced from 26 chaetognath samples representing 18 species. An unusual finding was the presence of two distinct classes of 28S rDNA gene in chaetognaths; our analyses suggest these arose by a gene (or gene cluster) duplication in a common ancestor of extant chaetognaths. The two classes of chaetognath 28S rDNA have been subject to different rates of molecular evolution; we present evidence that both are expressed and functional. In phylogenetic reconstructions, the two classes of 28S rDNA yield trees that root each other; these clearly demonstrate that the Aphragmophora and Phragmophora are natural groups. Within the Aphragmophora, we find good support for the groupings denoted Solidosagitta, Parasagitta, and Pseudosagitta. The relationships between several well-supported groups within the Aphragmophora are uncertain; we suggest this reflects rapid, recent radiation during chaetognath evolution. Received: 19 March 1996 / Accepted: 5 August 1996     

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     

An Estimate of Divergence Time of Parazoa and Eumetazoa and That of Cephalochordata and Vertebrata by Aldolase and Triose Phosphate Isomerase Clocks

Previously we suggested that four proteins including aldolase and triose phosphate isomerase (TPI) evolved with approximately constant rates over long periods covering the whole animal phyla. The constant rates of aldolase and TPI evolution were reexamined based on three different models for estimating evolutionary distances. It was shown that the evolutionary rates remain essentially unchanged in comparisons not only between different classes of vertebrates but also between vertebrates and arthropods and even between animals and plants, irrespective of the models used. Thus these enzymes might be useful molecular clocks for inferring divergence times of animal phyla. To know the divergence time of Parazoa and Eumetazoa and that of Cephalochordata and Vertebrata, the aldolase cDNAs from Ephydatia fluviatilis, a freshwater sponge, and the TPI cDNAs from Ephydatia fluviatilis and Branchiostoma belcheri, an amphioxus, have been cloned and sequenced. Comparisons of the deduced amino acid sequences of aldolase and TPI from the freshwater sponge with known sequences revealed that the Parazoa–Eumetazoa split occurred about 940 million years ago (Ma) as determined by the average of two proteins and three models. Similarly, the aldolase and TPI clocks suggest that vertebrates and amphioxus last shared a common ancestor around 700 Ma and they possibly diverged shortly after the divergence of deuterostomes and protostomes.     

Glycosylation of Acetylcholinesterase Forms in Microsomal Membranes from Normal and Dystrophic Lama2dy Mouse Muscle

Abstract: The distribution and glycosylation of acetylcholinesterase (AChE) forms in vesicles derived from sarcoplasmic reticulum of normal muscle (NMV) were investigated and compared with those from dystrophic muscle vesicles (DMV). AChE activity was similar in NMV and DMV. Most of the AChE in NMV and half in DMV were released with Triton X-100. Asymmetric (A₁₂) and globular hydrophilic and amphiphilic (G^H₄, G^A₄, G^A₂, and G^A₁) AChE species occurred in NMV and DMV, the lighter forms being predominant. The percentage of G^H₄ and G^A₄ decreased in DMV. A fraction of the AChE that could not be extracted with detergent was detached with collagenase. Most of the detergent-released A₁₂ AChE from NMV and nearly half in DMV failed to bind to Ricinus communis agglutinin (RCA-I). Conversely, the collagenase-detached isoforms bound to RCA, revealing that asymmetric AChE associated with internal membranes or basal lamina differed in glycosylation. Moreover, nearly half of G^A₄ AChE in DMV and a few in NMV bound to RCA. Most of the RCA-unreactive G^A₄ forms in NMV come from sarcolemma. The results indicate that dystrophy induces minor changes in the distribution and glycosylation of AChE forms in internal membranes of muscle.     

Combined use of 13C chemical shift and 1Hα−13Cα heteronuclear NOE data in monitoring a protein NMR structure refinement

Summary A large portion of the ¹³C resonance assignments for murine epidermal growth factor (mEGF) at pH 3.1 and 28°C has been determined at natural isotope abundance. Sequence-specific ¹³C assignments are reported for 100% of the assignable C, 96% of the C, 86% of the aromatic and 70% of the remaining peripheral aliphatic resonances of mEGF. A good correlation was observed between experimental and back-calculated C chemical shifts for regions of regular -sheet structure. These assignments also provide the basis for interpreting ¹H–¹³C heteronuclear NOE (HNOE) values in mEGF at natural isotope abundance. Some of the backbone polypeptide segments with high internal mobility, indicated by these ¹H–¹³C HNOE measurements, correlate with locations of residues involved in the putative mEGF-receptor binding site. Using four families of mEGF structures obtained over the last few years, we demonstrate that standard deviations between experimental and back-calculated C values can be used to monitor the refinement of this protein's structure, particularly for -sheet regions. Improved agreement between calculated and observed values of C is correlated with other measures of structure quality, including lowered values of residual constraint violations and more negative values of conformational energy. These results support the view that experimental conformation-dependent chemical shifts, C, can provide a reliable source of information for monitoring the process of protein structure refinement and are potentially useful restraints for driving the refinement.Abbreviations HSQC heteronuclear single-quantum coherence spectroscopy - PFG pulsed-field gradient - TOCSY ¹H-¹H total correlation spectroscopy - EGF epidermal growth factor - mEGF murine EGF - hEGF human EGF - hTGF human type- transforming growth factor - DIPSI spm-locking pulse sequence - NOE nuclear Overhauser effect - HNOE heteronuclear Overhauser effect     

Dynamic modelling of a helical peptide in solution using NMR data: Multiple conformations and multi-spin effects

Summary Nuclear Overhauser effect (NOE) measurements on molecules in solution provide information about only the ensemble-averaged properties of these molecules. An algorithm is presented that uses a list of NOEs to produce an ensemble of molecules that on average agrees with these NOEs, taking into account the effect of surrounding spins on the buildup of each NOE (spin diffusion). A simplified molecular dynamics simulation on several copies of the molecule in parallel is restrained by forces that are derived directly from differences between calculated and measured NOEs. The algorithm is tested on experimental NOE data of a helical peptide derived from bovine pancreatic trypsin inhibitor.