The sensitivity of the commonly used progressive multiple sequence alignment method has been greatly improved for the alignment of divergent protein sequences. Firstly, individual weights are assigned to each sequence in a partial alignment in order to down-weight near-duplicate sequences and up-weight the most divergent ones. Secondly, amino acid substitution matrices are varied at different alignment stages according to the divergence of the sequences to be aligned. Thirdly, residue-specific gap penalties and locally reduced gap penalties in hydrophilic regions encourage new gaps in potential loop regions rather than regular secondary structure. Fourthly, positions in early alignments where gaps have been opened receive locally reduced gap penalties to encourage the opening up of new gaps at these positions. These modifications are incorporated into a new program, CLUSTAL W which is freely available.

/pdf/clustal-w-improving-the-sensitivity-of-progressive-multiple-2na2s9sq2h.pdf

Clustal w: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

Some simple formulae were obtained which enable us to estimate evolutionary distances in terms of the number of nucleotide substitutions (and, also, the evolutionary rates when the divergence times are known). In comparing a pair of nucleotide sequences, we distinguish two types of differences; if homologous sites are occupied by different nucleotide bases but both are purines or both pyrimidines, the difference is called type I (or “transition” type), while, if one of the two is a purine and the other is a pyrimidine, the difference is called type II (or “transversion” type). Letting P and Q be respectively the fractions of nucleotide sites showing type I and type II differences between two sequences compared, then the evolutionary distance per site is K = — (1/2) ln {(1 — 2P — Q) }. The evolutionary rate per year is then given by k = K/(2T), where T is the time since the divergence of the two sequences. If only the third codon positions are compared, the synonymous component of the evolutionary base substitutions per site is estimated by K'S = — (1/2) ln (1 — 2P — Q). Also, formulae for standard errors were obtained. Some examples were worked out using reported globin sequences to show that synonymous substitutions occur at much higher rates than amino acid-altering substitutions in evolution.

A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences.

The University of Wisconsin Genetics Computer Group (UWGCG) has been organized to develop computational tools for the analysis and publication of biological sequence data. A group of programs that will interact with each other has been developed for the Digital Equipment Corporation VAX computer using the VMS operating system. The programs available and the conditions for transfer are described.

/pdf/a-comprehensive-set-of-sequence-analysis-programs-for-the-luj6w12mfm.pdf

A comprehensive set of sequence analysis programs for the VAX

The abbreviated name,‘mfold web server’,describes a number of closely related software applications available on the World Wide Web (WWW) for the prediction of the secondary structure of single stranded nucleic acids. The objective of this web server is to provide easy access to RNA and DNA folding and hybridization software to the scientific community at large. By making use of universally available web GUIs (Graphical User Interfaces),the server circumvents the problem of portability of this software. Detailed output,in the form of structure plots with or without reliability information,single strand frequency plots and ‘energy dot plots’, are available for the folding of single sequences. A variety of ‘bulk’ servers give less information,but in a shorter time and for up to hundreds of sequences at once. The portal for the mfold web server is http://www.bioinfo.rpi.edu/applications/ mfold. This URL will be referred to as ‘MFOLDROOT’.

/pdf/mfold-web-server-for-nucleic-acid-folding-and-hybridization-3pzoe1h05m.pdf

Mfold web server for nucleic acid folding and hybridization prediction

https://njc.rockefeller.edu/pdf4/Class2-LeeAmbros1993.pdf

The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14

Molecular genetics of the apolipoprotein B gene in pigs in relation to atherosclerosis.

Two strains of mice, DBA and C3H, have a fetal globin polypeptide chain which differs in electrophoretic mobility from the corresponding fetal chain of the C57B1 strain. Mice of the DBA and C3H strains also differ from those of the C57B1 in adult hemoglobin type. Results of backcrossing the (DBA x C57B1) hybrid to the C57B1 suggest that the fetal chain locus and the adult β-chain locus are closely linked.

http://science.sciencemag.org/content/sci/160/3830/885.full.pdf

Fetal Hemoglobin Variants in Mice

We describe two novel arrangements of the human fetal globin gene region: one chromosome with two linked A gamma genes (A gamma-A gamma) and two chromosomes with two linked G gamma genes (G gamma-G gamma). The gamma genes of these three chromosomes were cloned and the unusual 5' A gamma gene and one of the unusual 3' G gamma genes were partially sequenced. Both of these unusual genes differ from the genes normally found at their respective locations by a nucleotide substitution at the site of the single coding region difference between normal G gamma and A gamma genes. In both cases, the substitution is identical to the nucleotide found at that position in the normal neighboring gene. The unusual 3' G gamma gene also differs from normal A gamma genes at two other nucleotide positions, but both differences appear to be "private" or exclusive to this particular gene. These unusual fetal globin gene arrangements could have arisen from point mutations or from gene conversions of limited extent, the boundaries of which have been determined for all three chromosomes.

Two novel arrangements of the human fetal globin genes: Gγ-Gγ and Aγ-Aγ

Earlier studies of protein polymorphisms led to the description of 13 linked loci thought to encode the human salivary proline-rich proteins (PRPs). However, more recent studies at the DNA level have shown that there are only six genes which encode PRPs. The present study was undertaken in order to reconcile these observations. Nucleotide and decoded amino acid sequences from each of the six genes were compared with the available protein sequence data for PRPs. This analysis allowed assignment of the PmF, PmS and Pe proteins to the PRB1 locus, the G1 protein to the PRB3 locus, the Po protein to the PRB4 locus, the Ps protein to the PRB2 locus, and the CON1 and CON2 proteins to the PRB4 locus. Correlations between insertion/deletion RFLPs and PRP protein phenotypes were observed for the PmF, PmS, Gl and CON2 proteins. Our overall analysis indicates that in many instances several proteins previously considered to be the products of separate loci are actually proteolytic cleavage products of a large precursor specified by one or other of the six genes identified at the DNA level. Our analysis also demonstrates that some of the "null" alleles proposed to occur at 11 of the 13 loci in the earlier genetic studies, are actually productive alleles having alterations at proteolytic cleavage sites within the relevant precursor protein. The absence of cleavage leads to the persistence of longer precursor peptides not resolved electrophoretically, concurrently with an absence of the smaller PRPs seen when cleavage occurs.

https://www.genetics.org/content/genetics/120/1/255.full.pdf

Oliver Smithies

Papers

Molecular genetics of the apolipoprotein B gene in pigs in relation to atherosclerosis.

Fetal Hemoglobin Variants in Mice

Two novel arrangements of the human fetal globin genes: Gγ-Gγ and Aγ-Aγ

Many protein products from a few loci: assignment of human salivary proline-rich proteins to specific loci.

Gene action in the human haptoglobins: I. Dissociation into constituent polypeptide chains