Multiple DNA and protein sequence alignment based on segment-to-segment comparison.

In this paper, a new way to think about, and to construct, pairwise as well as multiple alignments of DNA and protein sequences is proposed. Rather than forcing alignments to either align single residues or to introduce gaps by defining an alignment as a path running right from the source up to the sink in the associated dot-matrix diagram, we propose to consider alignments as consistent equivalence relations defined on the set of all positions occurring in all sequences under consideration. We also propose constructing alignments from whole segments exhibiting highly significant overall similarity rather than by aligning individual residues. Consequently, we present an alignment algorithm that (i) is based on segment-to-segment comparison instead of the commonly used residue-to-residue comparison and which (ii) avoids the well-known difficulties concerning the choice of appropriate gap penalties: gaps are not treated explicity, but remain as those parts of the sequences that do not belong to any of the aligned segments. Finally, we discuss the application of our algorithm to two test examples and compare it with commonly used alignment methods. As a first example, we aligned a set of 11 DNA sequences coding for functional helix-loop-helix proteins. Though the sequences show only low overall similarity, our program correctly aligned all of the 11 functional sites, which was a unique result among the methods tested. As a by-product, the reading frames of the sequences were identified. Next, we aligned a set of ribonuclease H proteins and compared our results with alignments produced by other programs as reported by McClure et al. [McClure, M. A., Vasi, T. K. & Fitch, W. M. (1994) Mol. Biol. Evol. 11, 571-592]. Our program was one of the best scoring programs. However, in contrast to other methods, our protein alignments are independent of user-defined parameters.

Development of a bovine X chromosome linkage group and painting probes to assess cattle, sheep, and goat X chromosome segment homologies.

The X chromosome linkage group is conserved in placental mammals. However, X chromosome morphological differences, due to internal chromosome rearrangements, exist among mammalian species. We have developed bovine chromosome painting probes for Xp and Xq to assess segment homologies between the submetacentric bovine X chromosome and the acrocentric sheep and goat X chromosomes. These painting probes and their corresponding DNA libraries were developed by chromosome micromanipulation, DNA micropurification, microcloning, and PCR amplification. The bovine Xp painting probe identified an interstitially located homologous segment in the sheep and goat Xq region, most probably resulting from chromosome inversion. Ten type II (microsatellite) markers obtained from the bovine Xq library and five other X chromosome assigned, but unlinked, markers were used to generate a linkage map for Xq spanning 89.4 centimorgans. The chromosome painting probes and molecular markers generated in this study would be useful for comparative mapping and tracing of internal X chromosome rearrangements in all ruminant species and would contribute to the understanding of mammalian sex chromosome evolution.

A de novo missense mutation of the beta subunit of the epithelial sodium channel causes hypertension and Liddle syndrome, identifying a proline-rich segment critical for regulation of channel activity.

Liddle syndrome is a mendelian form of hypertension characterized by constitutively elevated renal Na reabsorption that can result from activating mutations in the beta or gamma subunit of the epithelial Na channel. All reported mutations have deleted the last 45-76 normal amino acids from the cytoplasmic C terminus of one of these channel subunits. While these findings implicate these terminal segments in the normal negative regulation of channel activity, they do not identify the amino acid residues that are critical targets for these mutations. Potential targets include the short highly conserved Pro-rich segments present in the C terminus of beta and gamma subunits; these segments are similar to SH3-binding domains that mediate protein-protein interaction. We now report a kindred with Liddle syndrome in which affected patients have a mutation in codon 616 of the beta subunit resulting in substitution of a Leu for one of these highly conserved Pro residues. The functional significance of this mutation is demonstrated both by the finding that this is a de novo mutation appearing concordantly with the appearance of Liddle syndrome in the kindred and also by the marked activation of amiloride-sensitive Na channel activity seen in Xenopus oocytes expressing channels containing this mutant subunit (8.8-fold increase compared with control oocytes expressing normal channel subunits; P = 0.003). These findings demonstrate a de novo missense mutation causing Liddle syndrome and identify a critical channel residue important for the normal regulation of Na reabsorption in humans.

Invariant chain made in Escherichia coli has an exposed N-terminal segment that blocks antigen binding to HLA-DR1 and a trimeric C-terminal segment that binds empty HLA-DR1.

Invariant chain (Ii), a membrane glycoprotein, binds class II major histocompatibility complex (MHC) glycoproteins, probably via its class II-associated Ii peptide (CLIP) segment, and escorts them toward antigen-containing endosomal compartments. We find that a soluble, trimeric ectodomain of Ii expressed and purified from Escherichia coli blocks peptide binding to soluble HLA-DR1. Proteolysis indicates that Ii contains two structural domains. The C-terminal two-thirds forms an alpha-helical domain that trimerizes and interacts with empty HLA-DR1 molecules, augmenting rather than blocking peptide binding. The N-terminal one-third, which inhibits peptide binding, is proteolytically susceptible over its entire length. In the trimer, the N-terminal domains act independently with each CLIP segment exposed and free to bind an MHC class II molecule, while the C-terminal domains act as a trimeric unit.

Expression of PTPH1, a rat protein tyrosine phosphatase, is restricted to the derivatives of a specific diencephalic segment.

Studies to date have identified only a few proteins that are expressed in a segment-specific manner within the mammalian brain. Here we report that a nonreceptor protein tyrosine phosphatase, PTPH1, is selectively expressed in the adult thalamus. Expression of PTPH1 mRNA is detected in most, but not all, thalamic nuclei. Nuclei that are derived embryonically from the dorsal thalamus and project to the neocortex express this gene, whereas those derived from the ventral thalamus do not. PTPH1 mRNA expression is also restricted to the dorsal thalamus during development and, thus, can serve as a specific marker for the dorsal thalamic nuclei. Since the subcellular localization of PTPH1 protein is not known, its functional role is not clear. However, the restriction of its expression to the thalamic nuclei that have thalamocortical connections suggests that PTPH1 may play a role in the maintenance of these connections or in determining the physiological properties of thalamic relay nuclei.

Binding of the transcription activator NRI (NTRC) to a supercoiled DNA segment imitates association with the natural enhancer: an electron microscopic investigation.

Electron microscopic visualization indicates that the transcription activator NRI (NTRC) binds with exceptional selectivity and efficiency to a sequence-induced superhelical (spiral) segment inserted upstream of the glnA promoter, accounting for its observed ability to substitute for the natural glnA enhancer. The cooperative binding of NRI to the spiral insert leads to protein oligomerization which, at higher concentration, promotes selective coating of the entire superhelical segment with protein. Localization of NRI at apical loops is observed with negatively supercoiled plasmid DNA. With a linear plasmid, bending of DNA is observed. We confirm that NRI is a DNA-bending protein, consistent with its high affinity for spiral DNA. These results prove that spiral DNA without any homology to the NRI-binding sequence site can substitute for the glnA enhancer by promoting cooperative activator binding to DNA and facilitating protein oligomerization. Similar mechanisms might apply to other prokaryotic and eukaryotic activator proteins that share the ability to bend DNA and act efficiently as multimers.

Relatedness threshold for the production of female sexuals in colonies of a polygynous ant, Myrmica tahoensis, as revealed by microsatellite DNA analysis.

The genetic relationships of colony members in the ant Myrmica tahoensis were determined on the basis of highly polymorphic microsatellite DNA loci. These analyses show that colonies fall into one of two classes. In roughly half of the sampled colonies, workers and female offspring appear to be full sisters. The remaining colonies contain offspring produced by two or more queens. Colonies that produce female sexuals are always composed of highly related females, while colonies that produce males often show low levels of nestmate relatedness. These results support theoretical predictions that workers should skew sex allocation in response to relatedness asymmetries found within colonies. The existence of a relatedness threshold below which female sexuals are not produced suggests a possible mechanism for worker perception of relatedness. Two results indicate that workers use genetic cues, not queen number, in making sex-allocation decisions. (i) The number of queens in a colony was not significantly correlated with either the level of relatedness asymmetry or the sex ratio. (ii) Sex-ratio shifts consistent with a genetically based mechanism of relatedness assessment were seen in an experiment involving transfers of larvae among unrelated nests. Thus workers appear to make sex-allocation decisions on the basis of larval cues and appear to be able to adjust sex ratios long after egg laying.