Evolutionary conservation of residues in vertebrate DNA polymerase N conferring low fidelity and bypass activity.

POLN is a nuclear A-family DNA polymerase encoded in vertebrate genomes. POLN has unusual fidelity and DNA lesion bypass properties, including strong strand displacement activity, low fidelity favoring incorporation of T for template G and accurate translesion synthesis past a 5S-thymine glycol (5S-Tg). We searched for conserved features of the polymerase domain that distinguish it from prokaryotic pol I-type DNA polymerases. A Lys residue (679 in human POLN) of particular interest was identified in the conserved 'O-helix' of motif 4 in the fingers sub-domain. The corresponding residue is one of the most important for controlling fidelity of prokaryotic pol I and is a nonpolar Ala or Thr in those enzymes. Kinetic measurements show that K679A or K679T POLN mutant DNA polymerases have full activity on nondamaged templates, but poorly incorporate T opposite template G and do not bypass 5S-Tg efficiently. We also found that a conserved Tyr residue in the same motif not only affects sensitivity to dideoxynucleotides, but also greatly influences enzyme activity, fidelity and bypass. Protein sequence alignment reveals that POLN has three specific insertions in the DNA polymerase domain. The results demonstrate that residues have been strictly retained during evolution that confer unique bypass and fidelity properties on POLN.

Crystal structures of progressive Ca2+ binding states of the Ca2+ sensor Ca2+ binding domain 1 (CBD1) from the CALX Na+/Ca2+ exchanger reveal incremental conformational transitions.

Na(+)/Ca(2+) exchangers (NCX) constitute a major Ca(2+) export system that facilitates the re-establishment of cytosolic Ca(2+) levels in many tissues. Ca(2+) interactions at its Ca(2+) binding domains (CBD1 and CBD2) are essential for the allosteric regulation of Na(+)/Ca(2+) exchange activity. The structure of the Ca(2+)-bound form of CBD1, the primary Ca(2+) sensor from canine NCX1, but not the Ca(2+)-free form, has been reported, although the molecular mechanism of Ca(2+) regulation remains unclear. Here, we report crystal structures for three distinct Ca(2+) binding states of CBD1 from CALX, a Na(+)/Ca(2+) exchanger found in Drosophila sensory neurons. The fully Ca(2+)-bound CALX-CBD1 structure shows that four Ca(2+) atoms bind at identical Ca(2+) binding sites as those found in NCX1 and that the partial Ca(2+) occupancy and apoform structures exhibit progressive conformational transitions, indicating incremental regulation of CALX exchange by successive Ca(2+) binding at CBD1. The structures also predict that the primary Ca(2+) pair plays the main role in triggering functional conformational changes. Confirming this prediction, mutagenesis of Glu(455), which coordinates the primary Ca(2+) pair, produces dramatic reductions of the regulatory Ca(2+) affinity for exchange current, whereas mutagenesis of Glu(520), which coordinates the secondary Ca(2+) pair, has much smaller effects. Furthermore, our structures indicate that Ca(2+) binding only enhances the stability of the Ca(2+) binding site of CBD1 near the hinge region while the overall structure of CBD1 remains largely unaffected, implying that the Ca(2+) regulatory function of CBD1, and possibly that for the entire NCX family, is mediated through domain interactions between CBD1 and the adjacent CBD2 at this hinge.

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.

Review: Families & Change: Coping With Stressful Events and Transitions

A review of Families & Change: Coping With Stressful Events and Transitions-4th edition by Sharon J. Price, Christine A. Price and Patrick C. McKenry.

Chromosomal organization and evolutionary conservation of the human chromosome 19q linkage group containing three DNA repair genes

A series of human-rodent somatic cell hybrids were investigated by Southern blot analysis for the presence or absence of twenty-six molecular markers and three isozyme loci from human chromosome 19. Based on the co-retention of these markers in the various independent hybrid clones containing portions of human chromosome 19 and on pulsed field mapping, chromosome 19 is divided into twenty ordered regions. The most likely marker order for the chromosome is: (LDLR, C3)-(cen-MANNB)-D19S7-PEPD-D19S9-GPI-TGF$ \beta$-(CYP2A, NCA, CGM2, BCKAD)-PSG1a-(D19S8, XRCC1)-(D19S19, ATP1A3)-(D19S37, APOC2)-CKMM-ERCC2-ERCC1-(D19S62, D19S51)-D19S6-D19S50-D19S22-(CGB, FTL)-qter.^ The region of 19q between the proximal marker D19S7 and the distal gene coding for the beta subunit of chorionic gonadotropin (CGB) is about 37 Mb in size and covers about 37 cM genetic distance. The ration of genetic to physical distance on 19q is therefore very close to the genomic average OF 1 cM/Mb. Estimates of physical distances for intervals between chromosome 19 markers were calculated using a mapping function which estimates distances based on the number of breaks in hybrid clone panels. The consensus genetic distances between individual markers (established at HBM10) were compared to these estimates of physical distances. The close agreement between the two estimates suggested that spontaneously broken hybrids are as appropriate for this type of study as radiation hybrids.^ All three DNA repair genes located on chromosome 19 were found to have homologues on Chinese hamster chromosome 9, which is hemizygous in CHO cells, providing an explanation for the apparent ease with which mutations at these loci were identified in CHO cells. Homologues of CKMM and TGF$\beta$ (from human chromosome 19q) and a mini-satellite DNA specific to the distal region of human chromosome 19q were also mapped to Chinese hamster 9. Markers from 19p did not map to this hamster chromosome. Thus the q-arm of chromosome 19, at least between the genes PEPD and ERCC1, appears to be a linkage group which is conserved intact between humans and Chinese hamsters. ^

Population genetics of VNTR loci and their applications in evolutionary studies

Variable number of tandem repeats (VNTR) are genetic loci at which short sequence motifs are found repeated different numbers of times among chromosomes. To explore the potential utility of VNTR loci in evolutionary studies, I have conducted a series of studies to address the following questions: (1) What are the population genetic properties of these loci? (2) What are the mutational mechanisms of repeat number change at these loci? (3) Can DNA profiles be used to measure the relatedness between a pair of individuals? (4) Can DNA fingerprint be used to measure the relatedness between populations in evolutionary studies? (5) Can microsatellite and short tandem repeat (STR) loci which mutate stepwisely be used in evolutionary analyses?^ A large number of VNTR loci typed in many populations were studied by means of statistical methods developed recently. The results of this work indicate that there is no significant departure from Hardy-Weinberg expectation (HWE) at VNTR loci in most of the human populations examined, and the departure from HWE in some VNTR loci are not solely caused by the presence of population sub-structure.^ A statistical procedure is developed to investigate the mutational mechanisms of VNTR loci by studying the allele frequency distributions of these loci. Comparisons of frequency distribution data on several hundreds VNTR loci with the predictions of two mutation models demonstrated that there are differences among VNTR loci grouped by repeat unit sizes.^ By extending the ITO method, I derived the distribution of the number of shared bands between individuals with any kinship relationship. A maximum likelihood estimation procedure is proposed to estimate the relatedness between individuals from the observed number of shared bands between them.^ It was believed that classical measures of genetic distance are not applicable to analysis of DNA fingerprints which reveal many minisatellite loci simultaneously in the genome, because the information regarding underlying alleles and loci is not available. I proposed a new measure of genetic distance based on band sharing between individuals that is applicable to DNA fingerprint data.^ To address the concern that microsatellite and STR loci may not be useful for evolutionary studies because of the convergent nature of their mutation mechanisms, by a theoretical study as well as by computer simulation, I conclude that the possible bias caused by the convergent mutations can be corrected, and a novel measure of genetic distance that makes the correction is suggested. In summary, I conclude that hypervariable VNTR loci are useful in evolutionary studies of closely related populations or species, especially in the study of human evolution and the history of geographic dispersal of Homo sapiens. (Abstract shortened by UMI.) ^

CHROMOSOMAL MAPPING, LINKAGE RELATIONSHIPS, AND EVOLUTIONARY STUDIES OF THE HUMAN ANONYMOUS DNA CLONE D1S1 (IN SITU HYBRIDIZATION, PRIMATES, GENE MAPPING, AUTOSOMAL DOMINANT RETINITIS PIGMENTOSA, GENE DUPLICATION)

D1S1, an anonymous human DNA clone originally called (lamda)Ch4-H3 or (lamda)H3, was the first single copy mapped to a human chromosome (1p36) by in situ hybridization. The chromosomal assignment has been confirmed in other laboratories by repeating the in situ hybridization but not by another method. In the present study, hybridization to a panel of hamster-human somatic cell hybrids revealed copies of D1S1 on both chromosomes 1 and 3. Subcloning D1S1 showed that the D1S1 clone itself is from chromosome 3, and the sequence detected by in situ hybridization is at least two copies of part of the chromosome 3 copy. This finding demonstrates the importance of verifying gene mapping with two methods and questions the accuracy of in situ hybridization mapping.^ Non-human mammals have only one copy of D1S1, and the non-human primate D1S1 map closely resembles the human chromosome 3 copy. Thus, the human chromosome 1 copies appear to be part of a very recent duplication that occurred after the divergence between humans and the other great apes.^ A moderately informative HindIII D1S1 RFLP was mapped to chromosome 3. This marker and 12 protein markers were applied to a linkage study of autosomal dominant retinitis pigmentosa (ADRP). None of the markers proved linkage, but adding the three families examined to previously published data raises the ADRP:Rh lod score to 1.92 at (THETA) = 0.30. ^

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. ^

MATHEMATICAL STUDIES ON THE EVOLUTIONARY CHANGE OF DNA SEQUENCES

With the aim of understanding the mechanism of molecular evolution, mathematical problems on the evolutionary change of DNA sequences are studied. The problems studied and the results obtained are as follows: (1) Estimation of evolutionary distance between nucleotide sequences. Studying the pattern of nucleotide substitution for the case of unequal substitution rates, a new mathematical formula for estimating the average number of nucleotide substitutions per site between two homologous DNA sequences is developed. It is shown that this formula has a wider applicability than currently available formulae. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. (2) Biases of the estimates of nucleotide substitutions obtained by the restriction enzyme method. The deviation of the estimate of nucleotide substitutions obtained by the restriction enzyme method from the true value is investigated theoretically. It is shown that the amount of the deviation depends on the nucleotides in the recognition sequence of the restriction enzyme used, unequal rates of substitution among different nucleotides, and nucleotide frequences, but the primary factor is the unequal rates of nucleotide substitution. When many different kinds of enzymes are used, however, the amount of average deviation is generally small. (3) Distribution of restriction fragment lengths. To see the effect of undetectable restriction fragments and fragment differences on the estimate of nucleotide differences, the theoretical distribution of fragment lengths is studied. This distribution depends on the type of restriction enzymes used as well as on the relative frequencies of four nucleotides. It is shown that undetectability of small fragments or fragment differences gives a serious underestimate of nucleotide substitutions when the length-difference method of estimation is used, but the extent of underestimation is small when the site-difference method is used. (4) Evolutionary relationships of DNA sequences in finite populations. A mathematical theory on the expected evolutionary relationships among DNA sequences (nucleons) randomly chosen from the same or different populations is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author). UMI ^

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.) ^

Molecular and evolutionary genetics of the X -linked visual pigment genes in humans and New World monkeys

Normal humans have one red and at least one green visual pigment genes. These genes are tightly linked as tandem repeats on the X chromosome and each of them has six exons. There is only one X-linked visual pigment gene in New World monkeys (NWMs) but the locus has three polymorphic alleles encoding red, yellow and green visual pigments, respectively. The spectral properties of the squirrel monkey and the marmoset (both NWMs) have been studied and partial sequences of the three alleles are available. To study the evolutionary history of these X-linked opsin genes in humans and NWMs, coding and intron sequences of the three squirrel monkey alleles and the three marmoset alleles were amplified by PCR followed by subcloning and sequencing. Introns 2 and 4 of the human red and green pigment genes were also sequenced. The results obtained are as follows: (1) The sequences of introns 2 and 4 of the human red and green opsin genes are significantly more similar between the two genes than are coding sequences, contrary to the usual situation where coding regions are better conserved in evolution than are introns. The high similarities in the two introns are probably due to recent gene conversion events during evolution of the human lineage. (2) Phylogenetic analysis of both intron and exon sequences indicates that the phylogenetic tree of the available primate opsin genes is the same as the species tree. The two human genes were derived from a gene duplication event after the divergence of the human and NWM lineages. The three alleles in each of the two NWM species diverged after the split of the two NWMs but have persisted in the population for at least 5 million years. (3) Allelic gene conversion might have occurred between the three squirrel monkey alleles. (4) A model of additive effect of hydroxyl-bearing amino acids on spectral tuning is proposed by treating some unknown variables as groups. Under the assumption that some residues have no effect, it is found that at least five amino acid residues, at positions 178 (3 nm), 180 (5 nm), 230 ($-$4 nm), 277 (9 nm) and 285 (13 nm), have linear spectral tuning effects. (5) Adaptive evolution of the opsin genes to different spectral peaks was observed at four residues that are important for spectral tuning. ^