Inferring phylogenies from DNA sequences of unequal base compositions.

A new method for computing evolutionary distances between DNA sequences is proposed. Contrasting with classical methods, the underlying model does not assume that sequence base compositions (A, C, G, and T contents) are at equilibrium, thus allowing unequal base compositions among compared sequences. This makes the method more efficient than the usual ones in recovering phylogenetic trees from sequence data when base composition is heterogeneous within the data set, as we show by using both simulated and empirical data. When applied to small-subunit ribosomal RNA sequences from several prokaryotic or eukaryotic organisms, this method provides evidence for an early divergence of the microsporidian Vairimorpha necatrix in the eukaryotic lineage.

Signaling by N- and C-terminal sequences of parathyroid hormone-related protein in hippocampal neurons.

Parathyroid hormone-related protein (PTHrP) is synthesized in the brain, and a single type of cloned receptor for the N-terminal portion of PTHrP and PTH is present in the central nervous system. Nothing is known about the physiological actions or signaling pathways used by PTHrP in the brain. Using cultured rat hippocampal neurons, we demonstrate that N-terminal PTHrP[1-34] and PTH[1-34] signal via cAMP and cytosolic calcium transients. The cAMP response showed strong acute (< or = 6 h) homologous and heterologous desensitization after preincubation with PTHrP or PTH. In contrast, the acute calcium response did not desensitize after preincubation with PTHrP; in fact, preincubation dramatically recruited additional responsive neurons. Unexpectedly, C-terminal PTHrP[107-139], which does not bind or activate the cloned PTH/PTHrP receptor, signaled in neurons via cytosolic calcium but not cAMP. Although some neurons responded to both PTHrP[1-34] and PTHrP[107-139], others responded only to PTHrP[1-34]. We conclude that certain hippocampal neurons exhibit dual signaling in response to PTHrP[1-34] and that some neurons have a receptor for C-terminal PTHrP that signals only via cytosolic calcium.

[Leu5]enkephalin-encoding sequences are targets for a specific DNA-binding factor.

A DNA-binding factor with high affinity and specificity for the [Leu5]enkephalin-encoding sequences in the prodynorphin and proenkephalin genes has been characterized. The factor has the highest affinity for the [Leu5]-enkephalin-encoding sequence in the dynorphin B-encoding region of the prodynorphin gene, has relatively high affinity for other [Leu5]enkephalin-encoding sequences in the prodynorphin and proenkephalin genes, but has no apparent affinity for similar DNA sequences coding for [Met5]-enkephalin in the prodynorphin or proopiomelanocortin genes. The factor has been named [Leu5]enkephalin-encoding sequence DNA-binding factor (LEF). LEF has a nuclear localization and is composed of three subunits of about 60, 70, and 95 kDa, respectively. The highest levels were observed in rat testis, cerebellum, and spleen and were generally higher in late embryonal compared to newborn or adult animals. LEF activity was also recorded in human clonal tumor cell lines. LEF inhibited the transcription of reporter genes in artificial gene constructs where a [Leu5]enkephalin-encoding DNA fragment had been inserted between the transcription initiation site and the coding region of the reporter genes. These observations suggest that the [Leu5]enkephalin-encoding sequences in the prodynorphin and proenkephalin genes also have regulatory functions realized through interaction with a specific DNA-binding factor.

Regulation of glycolipid synthesis in HL-60 cells by antisense oligodeoxynucleotides to glycosyltransferase sequences: effect on cellular differentiation.

Treatment of the human promyelocytic leukemia cell line HL-60 with antisense oligodeoxynucleotides to UDP-N-acetylgalactosamine:beta-1,4-N-acetylgalactosaminyl-transferase (GM2-synthase; EC 2.4.1.92) and CMP-sialic acid:alpha-2,8-sialyltransferase (GD3-synthase; EC 2.4.99.8) sequences effectively down-regulated the synthesis of more complex gangliosides in the ganglioside synthetic pathways after GM3, resulting in a remarkable increase in endogenous GM3 with concomitant decreases in more complex gangliosides. The treated cells underwent monocytic differentiation as judged by morphological changes, adherent ability, and nitroblue tetrazolium staining. These data provide evidence that the increased endogenous ganglioside GM3 may play an important role in regulating cellular differentiation and that the antisense DNA technique proves to be a powerful tool in manipulating glycolipid synthesis in the cell.

Toward the therapeutic editing of mutated RNA sequences.

If RNA editing could be rationally directed to mutated RNA sequences, genetic diseases caused by certain base substitutions could be treated. Here we use a synthetic complementary RNA oligonucleotide to direct the correction of a premature stop codon mutation in dystrophin RNA. The complementary RNA oligonucleotide was hybridized to a premature stop codon and the hybrid was treated with nuclear extracts containing the cellular enzyme double-stranded RNA adenosine deaminase. When the treated RNAs were translated in vitro, a dramatic increase in expression of a downstream luciferase coding region was observed. The cDNA sequence data are consistent with deamination of the adenosine in the UAG stop codon to inosine by double-stranded RNA adenosine deaminase. Injection of oligonucleotide-mRNA hybrids into Xenopus embryos also resulted in an increase in luciferase expression. These experiments demonstrate the principle of therapeutic RNA editing.

Conformation, energy, and folding ability of selected amino acid sequences.

Evolutionary selection of sequences is studied with a knowledge-based Hamiltonian to find the design principle for folding to a model protein structure. With sequences selected by naive energy minimization, the model structure tends to be unstable and the folding ability is low. Sequences with high folding ability have only the low-lying energy minimum but also an energy landscape which is similar to that found for the native sequence over a wide region of the conformation space. Though there is a large fluctuation in foldable sequences, the hydrophobicity pattern and the glycine locations are preserved among them. Implications of the design principle for the molecular mechanism of folding are discussed.

Identification of planarian homeobox sequences indicates the antiquity of most Hox/homeotic gene subclasses.

The homeotic gene complex (HOM-C) is a cluster of genes involved in the anteroposterior axial patterning of animal embryos. It is composed of homeobox genes belonging to the Hox/HOM superclass. Originally discovered in Drosophila, Hox/HOM genes have been identified in organisms as distantly related as arthropods, vertebrates, nematodes, and cnidarians. Data obtained in parallel from the organization of the complex, the domains of gene expression during embryogenesis, and phylogenetic relationships allow the subdivision of the Hox/HOM superclass into five classes (lab, pb/Hox3, Dfd, Antp, and Abd-B) that appeared early during metazoan evolution. We describe a search for homologues of these genes in platyhelminths, triploblast metazoans emerging as an outgroup to the great coelomate ensemble. A degenerate PCR screening for Hox/HOM homeoboxes in three species of triclad planarians has revealed 10 types of Antennapedia-like genes. The homeobox-containing sequences of these PCR fragments allowed the amplification of the homeobox-coding exons for five of these genes in the species Polycelis nigra. A phylogenetic analysis shows that two genes are clear orthologues of Drosophila labial, four others are members of a Dfd/Antp superclass, and a seventh gene, although more difficult to classify with certainty, may be related to the pb/Hox3 class. Together with previously identified Hox/HOM genes in other flatworms, our analyses demonstrate the existence of an elaborate family of Hox/HOM genes in the ancestor of all triploblast animals.

A human RNase E-like activity that cleaves RNA sequences involved in mRNA stability control.

We have detected an endoribonucleolytic activity in human cell extracts that processes the Escherichia coli 9S RNA and outer membrane protein A (ompA) mRNA with the same specificity as RNase E from E. coli. The human enzyme was partially purified by ion-exchange chromatography, and the active fractions contained a protein that was detected with antibodies shown to recognize E. coli RNase E. RNA containing four repeats of the destabilizing motif AUUUA and RNA from the 3' untranslated region of human c-myc mRNA were also found to be cleaved by E. coli RNase E and its human counterpart in a fashion that may suggest a role of this activity in mammalian mRNA decay. It was also found that RNA containing more than one AUUUA motif was cleaved more efficiently than RNA with only one or a mutated motif. This finding of a eukaryotic endoribonucleolytic activity corresponding to RNase E indicates an evolutionary conservation of the components of mRNA degradation systems.

Documentation of reticulate evolution in peonies (Paeonia) using internal transcribed spacer sequences of nuclear ribosomal DNA: implications for biogeography and concerted evolution.

The internal transcribed spacers (ITS) of nuclear ribosomal DNA of 33 species of genus Paeonia (Paeoniaceae) were sequenced. In section Paeonia, different patterns of nucleotide additivity were detected in 14 diploid and tetraploid species at sites that are variable in the other 12 species of the section, suggesting that reticulate evolution has occurred. Phylogenetic relationships of species that do not show additivity, and thus ostensibly were not derived through hybridization, were reconstructed by parsimony analysis. The taxa presumably derived through reticulate evolution were then added to the phylogenetic tree according to additivity from putative parents. The study provides an example of successfully using ITS sequences to reconstruct reticulate evolution in plants and further demonstrates that the sequence data could be highly informative and accurate for detecting hybridization. Maintenance of parental sequences in the species of hybrid origin is likely due to slowing of concerted evolution caused by the long generation time of peonies. The partial and uneven homogenization of parental sequences displayed in nine species of putative hybrid origin may have resulted from gradients of gene conversion. The documented hybridizations may have occurred since the Pleistocene glaciations. The species of hybrid origin and their putative parents are now distantly allopatric. Reconstruction of reticulate evolution with sequence data, therefore, provides gene records for distributional histories of some of the parental species.

Isolation of multiple sequences from the Plasmodium falciparum genome that encode conserved domains homologous to those in erythrocyte-binding proteins.

Open reading frames in the Plasmodium falciparum genome encode domains homologous to the adhesive domains of the P. falciparum EBA-175 erythrocyte-binding protein (eba-175 gene product) and those of the Plasmodium vivax and Plasmodium knowlesi Duffy antigen-binding proteins. These domains are referred to as Duffy binding-like (DBL), after the receptor that determines P. vivax invasion of Duffy blood group-positive human erythrocytes. Using oligonucleotide primers derived from short regions of conserved sequence, we have developed a reverse transcription-PCR method that amplifies sequences encoding the DBL domains of expressed genes. Products of these reverse transcription-PCR amplifications include sequences of single-copy genes (including eba-175) and variably transcribed genes that cross-hybridize to multiple regions of the genome. Restriction patterns of the multicopy genes show a high degree of polymorphism among different parasite lines, whereas single-copy genes are generally conserved. Characterization of the single-copy genes has identified a gene (ebl-1) that is related to eba-175 and is likely to be involved in erythrocyte invasion.

Molecular coevolution of mammalian ribosomal gene terminator sequences and the transcription termination factor TTF-I.

Both the DNA elements and the nuclear factors that direct termination of ribosomal gene transcription exhibit species-specific differences. Even between mammals--e.g., human and mouse--the termination signals are not identical and the respective transcription termination factors (TTFs) which bind to the terminator sequence are not fully interchangeable. To elucidate the molecular basis for this species-specificity, we have cloned TTF-I from human and mouse cells and compared their structural and functional properties. Recombinant TTF-I exhibits species-specific DNA binding and terminates transcription both in cell-free transcription assays and in transfection experiments. Chimeric constructs of mouse TTF-I and human TTF-I reveal that the major determinant for species-specific DNA binding resides within the C terminus of TTF-I. Replacing 31 C-terminal amino acids of mouse TTF-I with the homologous human sequences relaxes the DNA-binding specificity and, as a consequence, allows the chimeric factor to bind the human terminator sequence and to specifically stop rDNA transcription.

Coiled bodies contain U7 small nuclear RNA and associate with specific DNA sequences in interphase human cells.

Coiled bodies (CBs) are nuclear organelles whose structures appear to be highly conserved in evolution. In rapidly cycling cells, they are typically located in the nucleoplasm but are often found in contact with the nucleolus. The CBs in human cells contain a unique protein, called p80-coilin. Studies on amphibian oocyte nuclei have revealed a protein within the "sphere" organelle that shares significant structural similarity to p80-coilin. Spheres and CBs are also highly enriched in small nuclear ribonucleoproteins and other RNA-processing components. We present evidence that, like spheres, CBs contain U7 small nuclear RNA (snRNA) and associate with specific chromosomal loci. Using biotinylated 2'-O-methyl oligonucleotides complementary to the 5' end of U7 snRNA and fluorescence in situ hybridization, we show that U7 is distributed throughout the nucleoplasm, excluding nucleoli, and is concentrated in CBs. Interestingly, we found that CBs often associate with subsets of the histone, U1, and U2 snRNA gene loci in interphase HeLa-ATCC and HEp-2 monolayer cells. However, in a strain of suspension-grown HeLa cells, called HeLa-JS1000, we found a much lower rate of association between CBs and snRNA genes. Possible roles for CBs in the metabolism of these various histone and snRNAs are discussed.

The penicillin gene cluster is amplified in tandem repeats linked by conserved hexanucleotide sequences.

The penicillin biosynthetic genes (pcbAB, pcbC, penDE) of Penicillium chrysogenum AS-P-78 were located in a 106.5-kb DNA region that is amplified in tandem repeats (five or six copies) linked by conserved TTTACA sequences. The wild-type strains P. chrysogenum NRRL 1951 and Penicillium notatum ATCC 9478 (Fleming's isolate) contain a single copy of the 106.5-kb region. This region was bordered by the same TTTACA hexanucleotide found between tandem repeats in strain AS-P-78. A penicillin overproducer strain, P. chrysogenum E1, contains a large number of copies in tandem of a 57.9-kb DNA fragment, linked by the same hexanucleotide or its reverse complementary TGTAAA sequence. The deletion mutant P. chrysogenum npe10 showed a deletion of 57.9 kb that corresponds exactly to the DNA fragment that is amplified in E1. The conserved hexanucleotide sequence was reconstituted at the deletion site. The amplification has occurred within a single chromosome (chromosome I). The tandem reiteration and deletion appear to arise by mutation-induced site-specific recombination at the conserved hexanucleotide sequences.

Binding of purified multiple antibiotic-resistance repressor protein (MarR) to mar operator sequences.

Elevated expression of the marORAB multiple antibiotic-resistance operon enhances the resistance of Escherichia coli to various medically significant antibiotics. Transcription of the operon is repressed in vivo by the marR-encoded protein, MarR, and derepressed by salicylate and certain antibiotics. The possibility that repression results from MarR interacting with the marO operator-promoter region was studied in vitro using purified MarR and a DNA fragment containing marO. MarR formed at least two complexes with marO DNA, bound > 30-fold more tightly to it than to salmon sperm DNA, and protected two separate 21-bp sites within marO from digestion by DNase I. Site I abuts the downstream side of the putative -35 transcription-start signal and includes 4 bp of the -10 signal. Site II begins 13 bp downstream of site I, ending immediately before the first base pair of marR. Site II, approximately 80% homologous to site I, is not required for repression since a site II-deleted mutant (marO133) was repressed in trans by wild-type MarR. The absence of site II did not prevent MarR from complexing with the site I of marO133. Salicylate bound to MarR (Kd approximately 0.5 mM) and weakened the interaction of MarR with sites I and II. Thus, repression of the mar operon, which curbs the antibiotic resistance of E. coli, correlates with the formation of MarR-site I complexes. Salicylate appears to induce the mar operon by binding to MarR and inhibiting complex formation, whereas tetracycline and chloramphenicol, which neither bind MarR nor inhibit complex formation, must induce by an indirect mechanism.

Evidence for five divergent thioredoxin h sequences in Arabidopsis thaliana.

Five different clones encoding thioredoxin homologues were isolated from Arabidopsis thaliana cDNA libraries. On the basis of the sequences they encode divergent proteins, but all belong to the cytoplasmic thioredoxins h previously described in higher plants. The five proteins obtained by overexpressing the coding sequences in Escherichia coli present typical thioredoxin activities (NADP(+)-malate dehydrogenase activation and reduction by Arabidopsis thioredoxin reductase) despite the presence of a variant active site, Trp-Cys-Pro-Pro-Cys, in three proteins in place of the canonical Trp-Cys-Gly-Pro-Cys sequence described for thioredoxins in prokaryotes and eukaryotes. Southern blots show that each cDNA is encoded by a single gene but suggest the presence of additional related sequences in the Arabidopsis genome. This very complex diversity of thioredoxins h is probably common to all higher plants, since the Arabidopsis sequences appear to have diverged very early, at the beginning of plant speciation. This diversity allows the transduction of a redox signal into multiple pathways.