THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

Molluscan memory of injury: evolutionary insights into chronic pain and neurological disorders.

Molluscan preparations have yielded seminal discoveries in neuroscience, but the experimental advantages of this group have not, until now, been complemented by adequate molecular or genomic information for comparisons to genetically defined model organisms in other phyla. The recent sequencing of the transcriptome and genome of Aplysia californica, however, will enable extensive comparative studies at the molecular level. Among other benefits, this will bring the power of individually identifiable and manipulable neurons to bear upon questions of cellular function for evolutionarily conserved genes associated with clinically important neural dysfunction. Because of the slower rate of gene evolution in this molluscan lineage, more homologs of genes associated with human disease are present in Aplysia than in leading model organisms from Arthropoda (Drosophila) or Nematoda (Caenorhabditis elegans). Research has hardly begun in molluscs on the cellular functions of gene products that in humans are associated with neurological diseases. On the other hand, much is known about molecular and cellular mechanisms of long-term neuronal plasticity. Persistent nociceptive sensitization of nociceptors in Aplysia displays many functional similarities to alterations in mammalian nociceptors associated with the clinical problem of chronic pain. Moreover, in Aplysia and mammals the same cell signaling pathways trigger persistent enhancement of excitability and synaptic transmission following noxious stimulation, and these highly conserved pathways are also used to induce memory traces in neural circuits of diverse species. This functional and molecular overlap in distantly related lineages and neuronal types supports the proposal that fundamental plasticity mechanisms important for memory, chronic pain, and other lasting alterations evolved from adaptive responses to peripheral injury in the earliest neurons. Molluscan preparations should become increasingly useful for comparative studies across phyla that can provide insight into cellular functions of clinically important genes.

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^