Evolution of influenza virus genes

The nucleotide sequences of the eight different influenza A virus segments (genes) were compared among 14 different subtypes. These comparisons demonstrate the presence of molecular clocks in the viral genes; they accumulated both silent and amino acid-changing substitutions at approximately constant rates with respect to time during evolution. In addition, comparison of the rates of evolution among the eight viral genes, excluding the P2 gene, revealed a rapid and roughly equal rate of silent substitution for different genes. The P2 gene exception is explained as the result of recombination (reassortment) between distantly related strains. The rate of amino acid-changing substitution differs greatly from gene to gene. The rate of silent substitution was estimated to be 1.1 X 10(-2)/site/year on the average--that is, about 2 X 10(6) times higher than eukaryotic gene equivalents, which is remarkable. Strain A/USSR/90/77 was shown to evolve with a rate that is similar to those of other strains but to behave as if replication was frozen during a certain period (Nakajima et al. 1978). The frozen period was estimated to be 25 yr on the basis of the molecular clock. A similar analysis revealed another example of frozen replication--in this case, apparently for a period of about 9 yr- -in a duck strain, A/duck/Ontario/77.      

A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes

A new method is proposed for estimating the number of synonymous and nonsynonymous nucleotide substitutions between homologous genes. In this method, a nucleotide site is classified as nondegenerate, twofold degenerate, or fourfold degenerate, depending on how often nucleotide substitutions will result in amino acid replacement; nucleotide changes are classified as either transitional or transversional, and changes between codons are assumed to occur with different probabilities, which are determined by their relative frequencies among more than 3,000 changes in mammalian genes. The method is applied to a large number of mammalian genes. The rate of nonsynonymous substitution is extremely variable among genes; it ranges from 0.004 X 10(-9) (histone H4) to 2.80 X 10(-9) (interferon gamma), with a mean of 0.88 X 10(-9) substitutions per nonsynonymous site per year. The rate of synonymous substitution is also variable among genes; the highest rate is three to four times higher than the lowest one, with a mean of 4.7 X 10(-9) substitutions per synonymous site per year. The rate of nucleotide substitution is lowest at nondegenerate sites (the average being 0.94 X 10(-9), intermediate at twofold degenerate sites (2.26 X 10(-9)). and highest at fourfold degenerate sites (4.2 X 10(-9)). The implication of our results for the mechanisms of DNA evolution and that of the relative likelihood of codon interchanges in parsimonious phylogenetic reconstruction are discussed.     

Evolution of glucagon genes

Statistical analyses of DNA sequences of the preproglucagon genes from bovine, human, hamster, and anglerfish suggest that a gene duplication creating two anglerfish genes (AF I and II) occurred about 160 Myr ago, long after the separation of fish and mammals. The analyses further suggest that the internal duplication producing the glucagon and glucagon-like peptide II (GLP-II) regions occurred about 1.2 billion years ago, which would indicate that the GLP-II region was present in the ancestral anglerfish sequence but was silenced or deleted before the gene duplication separating AF I and II. The glucagon-like peptide I (GLP-I) was derived from a duplication of the ancestral glucagon region about 800 Myr ago. The rate of synonymous substitution in these genes is approximately 4.3 x 10(-9) substitutions per year per synonymous site. The rate of nonsynonymous substitution in the signal peptide region is about 1.1 x 10(-9) substitutions per year per nonsynonymous site, a high rate comparable to that in the C-peptide region of preproinsulin. The rate of nonsynonymous substitution in the glicentin-related pancreatic polypeptide (GRPP) region is 0.63 x 10(-9) for the comparisons between mammalian species and 1.8 x 10(-9) for the comparisons between fish and mammals; the moderate rate in mammals suggests a physiological role for GRPP. The glucagon region is extremely conservative; no nonsynonymous substitution is observed in the mammalian genes, and a nonsynonymous rate of 0.18 x 10(-9) was obtained from the comparisons between fish and mammals. In the GLP-I region, the rate of nonsynonymous substitution was estimated to be 0.08 x 10(-9) for the comparisons between mammalian species and 0.30 x 10(- 9) for the comparisons between fish and mammals. In the GLP-II region, the rate was estimated to be 0.25 x 10(-9) for the comparisons between mammalian species. Thus, GLP-I and II are also very conservative, which suggests an important physiological role for these peptides.      

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

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

Comparative analysis of human masculinity

To study rapidly evolving male specific Y (MSY) genes we retrieved and analyzed nine such genes. VCY, HSFY and RBMY were found to have functional X gametologs, but the rest did not. Using chimpanzee orthologs for XKRY, CDY, HSFY, PRY, and TSPY, the average silent substitution is estimated as 0.017 +/- 0.006/site and the substitution rate is 1.42 x 10(-9)/site/year. Except for VCY, all other loci possess two or more pseudogenes on the Y chromosome. Sequence differences from functional genes show that BPY2, DAZ, XKRY, and RBMY each have one pseudogene for each one that is human specific, while others were generated well before the human-chimpanzee split, by means of duplication, retro-transposition or translocation. Some functional MSY gene duplication of VCY, CDY and HSFY, as well as X-linked VCX and HSFX duplication, occurred in the lineage leading to humans; these duplicates have accumulated nucleotide substitutions that permit their identification.     

Evolution in bacteria: Evidence for a universal substitution rate in cellular genomes

Summary This paper constructs a temporal scale for bacterial evolution by tying ecological events that took place at known times in the geological past to specific branch points in the genealogical tree relating the 16S ribosomal RNAs of eubacteria, mitochondria, and chloroplasts. One thus obtains a relationship between time and bacterial RNA divergence which can be used to estimate times of divergence between other branches in the bacterial tree. According to this approach,Salmonella typhimurium andEscherichia coli diverged between 120 and 160 million years (Myr) ago, a date which fits with evidence that the chief habitats occupied now by these two enteric species became available that long ago.The median extent of divergence betweenS. typhimurium andE. coli at synonymous sites for 21 kilobases of protein-coding DNA is 100%. This implies a silent substitution rate of 0.7–0.8%/Myr—a rate remarkably similar to that observed in the nuclear genes of mammals, invertebrates, and flowering plants. Similarities in the substitution rates of eucaryotes and procaryotes are not limited to silent substitutions in protein-coding regions. The average substitution rate for 16S rRNA in eubacteria is about 1%/50 Myr, similar to the average rate for 18S rRNA in vertebrates and flowering plants. Likewise, we estimate a mean rate of roughly 1%/25 Myr for 5S rRNA in both eubacteria and eucaryotes.For a few protein-coding genes of these enteric bacteria, the extent of silent substitution since the divergence ofS. typhimurium andE. coli is much lower than 100%, owing to extreme bias in the usage of synonymous codons. Furthermore, in these bacteria, rates of amino acid replacement were about 20 times lower, on average, than the silent rate. By contranst, for the mammalian genes studied to date, the average replacement rate is only four to five times lower than the rate of silent substitution.     

Neuroglobins from the zebrafish Danio rerio and the pufferfish Tetraodon nigroviridis

Neuroglobin is a recently discovered respiratory, porphyrin-containing protein that is expressed in the brain of mouse and man. Here we show that neuroglobin is also present in the teleost fish. Complete cDNA sequences are reported from the pufferfish Tetraodon nigroviridis and the zebrafish Danio rerio. In addition, the neuroglobin gene of T. nigroviridis was sequenced, demonstrating the conservation of the B12.2, E11.0 and G7.0 introns plus the presence of an additional intron in the 5' noncoding region. The fish neuroglobins each comprise 159 amino acids and are 84.3% identical. Phylogenetic analyses show a basal position of the neuroglobins within the metazoan globin tree. An enhanced amino acid substitution rate was estimated for the fish neuroglobins ( approximately 0.93 x 10(-9) amino acid substitutions per site and year) compared with their mammalian proteins ( approximately 0.39 x 10(-9) replacements per site and year).     

Rates of Synonymous Substitution and Base Composition of Nuclear Genes in Drosophila

We compared the rates of synonymous (silent) substitution among various genes in a number of species of Drosophila. First, we found that even for a particular gene, the rate of synonymous substitution varied considerably with Drosophila lineages. Second, we showed a large variation in synonymous substitution rates among nuclear genes in Drosophila. These rates of synonymous substitution were correlated negatively with C content and positively with A content at the third codon positions. Nucleotide sequences were also compared between pseudogenes and their functional homologs. The C content of the pseudogenes was lower than that of the functional genes and the A content of the former was higher than that of the latter. Because the synonymous substitution for functional genes and the nucleotide substitution for pseudogenes are exempted from any selective constraint at the protein level, these observations could be explained by a biased pattern of mutation in the Drosophila nuclear genome. Such a bias in the mutation pattern may affect the molecular clock (local clock) of each nuclear gene of each species. Finally, we obtained the average rates of synonymous substitution for three gene groups in Drosophila; 11.0 x 10(-9), 17.5 x 10(-9) and 27.1 x 10(-9)/site/year.     

Molecular Evolution of the Small Subunit of Ribulose Bisphosphate Carboxylase: Nucleotide Substitution and Gene Conversion

The nucleotide sequences encoding the mature portion of 31 ribulose 1.5-bisphosphate carboxylase small subunit (SSU) genes from 17 genera of plants, green algae and cyanobacteria were examined. Among the 465 pairwise sequence comparisons, SSU multigene family members within the same species were more similar to each other in nonsynonymous or replacement nucleotide substitutions (RNS) than they were to SSU sequences in any other organism. The concerted evolution of independent SSU gene lineages within closely related plant species suggests that homogenization of RNS positions has occurred at least once in the life of each genus. The rate of expected RNS among mature SSU sequences was calculated to be 1.25 X 10(-9)/site/yr for the first 70 million years (MY) of divergence with a significant slowing to 0.13 X 10(-9)/site/yr for the next 1,400 MY. The data suggest that mature SSU sequences do not accumulate more than 20% differences in the RNS positions without compensatory changes in other components of this enzyme system. During the first 70 MY of divergence between species, the rate of expected synonymous or silent nucleotide substitutions (SNS) is approximately 6.6 X 10(-9)/site/yr. This is five times the RNS rate and is similar to the silent rate observed in animals. In striking contrast, SNS and RNS do not show this correlation among SSU gene family members within a species. A mechanism involving gene conversion within the exons followed by selection for biased gene conversion products with conservation of RNS positions and divergence of SNS positions is discussed. A SSU gene tree based on corrected RNS for 31 SSU sequences is presented and agrees well with a species tree based on morphological and cytogenetic traits for the 17 genera examined. SSU gene comparisons may be useful in predicting phylogenetic relationships and in some cases divergence times of various plant, algal and cyanobacterial species.     

GnRH-associated peptide (GAP) is cosecreted with GnRH into the hypophyseal portal blood of ovariectomized sheep

The human amyloid beta protein is a major component of brain amyloid found in patients with Alzheimer's disease. As an initial step to understand the biological function of its precursor protein, we have isolated cDNA for the mouse homolog of the human beta protein precursor. Comparison of the predicted amino acid sequence with that of human revealed a quite high degree of homology (96.8%), and the calculated evolutionary rate of the mRNA at amino acid substitution site was relatively low (0.1 x 10(-9)/site/year). The mRNA was abundant in brain and kidney, and also detected in other tissues at low level. These results indicated that this protein is highly conserved through mammalian evolution and may be involved in a basic biological process(es).     

Xanthine Dehydrogenase (XDH): Episodic Evolution of a ``Neutral' Protein

We investigated the evolution of xanthine dehydrogenase (Xdh) in 34 species from the three multicellular kingdoms, including one plant, two fungi, and three animal phyla, two classes of vertebrates, four orders of mammals, and two orders of insects. We adopted a model-based maximum-likelihood framework of inference. After accounting for among-site rate variation and heterogeneous nucleotide composition of the sequences using the discrete gamma distribution, and using nonhomogeneous nonstationary representations of the substitution process, the rate of amino acid replacement is 30.4 x 10(-10)/site/year when Drosophila species are compared but only approximately 18 x 10(-10)/site/year when comparisons are made between mammal orders, between insect orders, or between different animal phyla and approximately 11 x 10(-10)/site/year when comparisons are made between birds and mammals, between fungi, or between the three multicellular kingdoms. To account for these observations, the rate of amino acid replacement must have been eight or more times higher in some lineages and at some times than in others. Spastic evolution of Xdh appears to be related to the particularities of the genomes in which the locus is embedded.     

Evolution of nucleic acids coding for ribonucleases: the mRNA sequence of mouse pancreatic ribonuclease

The cDNA of mouse pancreatic mRNA has been cloned. After the library was screened with a rat ribonuclease cDNA probe, the positive clones were isolated and sequenced. There were no differences from the previously determined protein sequence. The mRNA codes for a preribonuclease of 149 amino acid residues including a signal peptide of 25 amino acids. The 3' noncoding region has a length of 260 bp, and the total mRNA length is approximately 940 bp. Comparison with the rat pancreatic ribonuclease sequence showed a high rate of nucleotide substitution. Within the coding region, nonsynonymous and synonymous substitution rates are 4.3 X 10(-9) and 15 X 10(-9) nucleotide substitutions/site/year, respectively. The latter value is one of the highest rates observed in the molecular evolution of mammalian nuclear genes. In the signal sequences the synonymous substitution rate is much lower and about the same as the nonsynonymous rate. Signal sequences of other mouse and rat proteins also exhibit little difference between synonymous and nonsynonymous rates. The sequences of rat and mouse pancreatic ribonuclease messengers were compared with those of bovine pancreatic, seminal, and brain ribonuclease. While the 3' noncoding regions of rat and mouse are very similar, as are those of the three bovine messengers, there is no significant similarity between both rodent and the three bovine messengers for the greater part of these regions. There is a duplication of approximately 50 nucleotides in the 3' noncoding region of the bovine messengers, with a region rich in A and C in between. The presence of this structural feature may be correlated with recent gene duplications that have occurred in the bovine genome.     

Human alpha-fetoprotein. Fluorescence studies on binding and proximity relationships for fatty acids and bilirubin.

The binding of bilirubin and the polyene fatty acids cis-parinaric acid and cis-eleostearic acid to human alpha-fetoprotein was studied using fluorescence quenching and fluorescence enhancement techniques. alpha-Fetoprotein has three fatty acid binding sites of decreasing affinity (association constants 2.1 x 10(7) M-1 9.1 X 10(5) M-1, and 1.4 x 10(5) M-1) and one relatively strong and one relatively weak bilirubin binding site (association constants 1.1 x 10(7) M-1 and 1.8 x 10(5) M-1). These association constants are slightly weaker than the corresponding association constants for binding to human albumin. Competition experiments failed to show preferential binding of polyunsaturated fatty acids. Fluorescence quenching was used to determine 11 ligand-ligand and ligand-tryptophanyl residue distances. Each of these 11 calculated distances (ranging from 19 A to 32 A) was within 5 A of the corresponding distances measured previously for human albumin (Berde, C.B., Hudson, B.S., Simoni, R.D., and Sklar, L.A. 1979, J. Biol. Chem. 254, 391-400). Thus, in addition to previously described sequence homology, immunologic cross-reactivity, and other similarities, human albumin and human alpha-fetoprotein have functional and geometric homologies.     

Hormonal regulation during secretion of alpha-fetoprotein in hepatoma cells grown in synthetic medium

The variant cell line of H4-II-E-C3 cells derived from the Reuber H-35 hepatoma cells has been established using protein- and lipid-free synthetic medium. This H4-II-E-C3-V line can synthesize and secrete considerable amounts of alpha-fetoprotein (AFP) and albumin. The addition of 5 X 10(-7) M dexamethasone to the medium stimulated the excretion of AFP without increasing total AFP synthesis, whereas 8.7 X 10(-8) M insulin inhibited the excretion of AFP without a significant inhibition of intracellular AFP synthesis. However, neither dexamethasone nor insulin altered either the cellular or secreted levels of albumin. Cells were pulse labeled with [35S]methionine and then chased after addition of excess unlabeled methionine. AFP appeared in the medium after 10 min, and 50% of the protein was secreted after 110 min. The rate of secretion of AFP was much slower than that of albumin, 50% of which was secreted after 25 min. Dexamethasone, 5 X 10(-7) M, caused a marked enhancement in the rate of AFP secretion, with 50% released after 75 min. Insulin, 8.7 X 10(-8) M, by contrast, caused a marked delay in AFP secretion with only 20% released after 180 min and then a plateau was approached. Since the intracellular AFP was excreted 55% after 180 min the remaining 25% of newly made AFP was suggested to be degraded during secretion. The kinetics of movement of AFP during secretion and endoglycosidase H treatment of intracellular and secreted AFP suggested that insulin impeded the transport of AFP from the rough endoplasmic reticulum to the Golgi apparatus.     

Alpha-fetoprotein receptors in a human breast cancer cell line

Evidence is presented for the existence of specific receptors for alpha-fetoprotein on the surface of MCF-7 human breast cancer cells. At 4 degrees C, the binding of alpha-fetoprotein to these cells displayed a biphasic saturation curve. Scatchard analysis revealed the presence of at least two binding sites with dissociation constants of 4.5 X 10(-9) M (2,000 sites/cell) and 1.3 X 10(-8) M (135,000 sites/cell), respectively. Binding was inhibited by 85% in the presence of a 5,000-fold excess of unlabeled alpha-fetoprotein and by 50% with the same excess of serum albumin. Competition by other serum proteins was not significant. At 37 degrees C, alpha-fetoprotein was endocytosed and the uptake curve reached a plateau after 3-4 hours of incubation.     

Synthesis and characterization of a recombinant fragment of human alpha-fetoprotein with antigenic selectivity versus albumin

A DNA sequence coding for human alpha-fetoprotein amino acid sequence 38-119 was synthesized and cloned in a bacterial expression vector. The alpha-fetoprotein sequence was selected as the least homologous to albumin, since the two proteins have an overall amino acid identity of approximately 38%. A chimeric protein was obtained which was purified by preparative electrophoresis and characterized in its primary structure by fast atom bombardment mass spectometry. About 70% of the alpha-fetoprotein sequence was physically mapped and found to correspond to the amino acids encoded in the synthetic gene. The use of this recombinant protein allowed the selection of monoclonal antibodies recognizing both the recombinant fragment and native alpha-fetoprotein. These antibodies should allow the development of an immunoassay for alpha-fetoprotein with absolute selectivity versus albumin. This might result in more sensitive clinical determinations, avoiding the possibility of cross-reactions.     

Amino acid sequence of the ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit from Euglena

Summary The amino acid sequence of the ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) small subunit (SSU) from Euglena has been established by alignment of the sequence of peptides obtained by cleavage with chymotrypsin, trypsin, Staphylococcus aureus protease or formic acid. The Euglena SSU has 138 amino acids and thus represents longest SSU sequence described so far. Homology is only 41% with cyanobacteria SSU and about 51% with higher plant SSU, whereas it is around 75% between higher plants. The largest homologous portion between all the known SSU sequences is localized in the second half and covers about 20 amino acids. The phylogenetic tree based on known SSU sequences has been established and the rate of amino acid substitution for SSU is estimated to be about 1.35×10^-9 per year and per site. Despite heterogeneity in amino acid sequence, we found that the overall secondary structure is fairly well conserved.Abbreviations DABITC Dimethyl amino azobenzene isothiocyanate - HPLC high pressure liquid chromatography - Kd Kilo daltons - LSU large subunit - PITC phenyl isothiocyanate - RuBisCO ribulose-1,5-bisphosphate carboxylase/oxygenase - SDS sodium dodecyl sulfate - SSU small subunit - TFA trifluoric acetic acid     

Molecular drift of the bride of sevenless (boss) gene in Drosophila

DNA sequences were determined for three to five alleles of the bride-of- sevenless (boss) gene in each of four species of Drosophila. The product of boss is a transmembrane receptor for a ligand coded by the sevenless gene that triggers differentiation of the R7 photoreceptor cell in the compound eye. Population parameters affecting the rate and pattern of molecular evolution of boss were estimated from the multinomial configurations of nucleotide polymorphisms of synonymous codons. The time of divergence between D. melanogaster and D. simulans was estimated as approximately 1 Myr, that between D. teissieri and D. yakuba as approximately 0.75 Myr, and that between the two pairs of sibling species as approximately 2 Myr. (The boss genes themselves have estimated divergence times approximately 50% greater than the species divergence times.) The effective size of the species was estimated as approximately 5 x 10(6), and the average mutation rate was estimated as 1-2 x 10(-9)/nucleotide/generation. The ratio of amino acid polymorphisms within species to fixed differences between species suggests that approximately 25% of all possible single-step amino acid replacements in the boss gene product may be selectively neutral or nearly neutral. The data also imply that random genetic drift has been responsible for virtually all of the observed differences in the portion of the boss gene analyzed among the four species.