Bootstrap confidence levels for phylogenetic trees

Evolutionary trees are often estimated from DNA or RNA sequence data. How much confidence should we have in the estimated trees? In 1985, Felsenstein [Felsenstein, J. (1985) Evolution 39, 783–791] suggested the use of the bootstrap to answer this question. Felsenstein’s method, which in concept is a straightforward application of the bootstrap, is widely used, but has been criticized as biased in the genetics literature. This paper concerns the use of the bootstrap in the tree problem. We show that Felsenstein’s method is not biased, but that it can be corrected to better agree with standard ideas of confidence levels and hypothesis testing. These corrections can be made by using the more elaborate bootstrap method presented here, at the expense of considerably more computation.

Molecular and Enzymatic Characterization of Three Phosphoinositide-Specific Phospholipase C Isoforms from Potato1

Many cellular responses to stimulation of cell-surface receptors by extracellular signals are transmitted across the plasma membrane by hydrolysis of phosphatidylinositol-4,5-bisphosphate (PIP2), which is cleaved into diacylglycerol and inositol-1,4,5-tris-phosphate by phosphoinositide-specific phospholipase C (PI-PLC). We present structural, biochemical, and RNA expression data for three distinct PI-PLC isoforms, StPLC1, StPLC2, and StPLC3, which were cloned from a guard cell-enriched tissue preparation of potato (Solanum tuberosum) leaves. All three enzymes contain the catalytic X and Y domains, as well as C2-like domains also present in all PI-PLCs. Analysis of the reaction products obtained from PIP2 hydrolysis unequivocally identified these enzymes as genuine PI-PLC isoforms. Recombinant StPLCs showed an optimal PIP2-hydrolyzing activity at 10 μm Ca2+ and were inhibited by Al3+ in equimolar amounts. In contrast to PI-PLC activity in plant plasma membranes, however, recombinant enzymes could not be activated by Mg2+. All three stplc genes are expressed in various tissues of potato, including leaves, flowers, tubers, and roots, and are affected by drought stress in a gene-specific manner.

Bootstrap confidence levels for phylogenetic trees.

Evolutionary trees are often estimated from DNA or RNA sequence data. How much confidence should we have in the estimated trees? In 1985, Felsenstein [Felsenstein, J. (1985) Evolution 39, 783-791] suggested the use of the bootstrap to answer this question. Felsenstein's method, which in concept is a straightforward application of the bootstrap, is widely used, but has been criticized as biased in the genetics literature. This paper concerns the use of the bootstrap in the tree problem. We show that Felsenstein's method is not biased, but that it can be corrected to better agree with standard ideas of confidence levels and hypothesis testing. These corrections can be made by using the more elaborate bootstrap method presented here, at the expense of considerably more computation.

The human U5-200kD DEXH-box protein unwinds U4/U6 RNA duplices in vitro

Splicing of nuclear precursors of mRNA (pre-mRNA) involves dynamic interactions between the RNA constituents of the spliceosome. The rearrangement of RNA–RNA interactions, such as the unwinding of the U4/U6 duplex, is believed to be driven by ATP-dependent RNA helicases. We recently have shown that spliceosomal U5 small nuclear ribonucleoproteins (snRNPs) from HeLa cells contain two proteins, U5–200kD and U5–100kD, which share homology with the DEAD/DEXH-box families of RNA helicases. Here we demonstrate that purified U5 snRNPs exhibit ATP-dependent unwinding of U4/U6 RNA duplices in vitro. To identify the protein responsible for this activity, U5 snRNPs were depleted of a subset of proteins under high salt concentrations and assayed for RNA unwinding. The activity was retained in U5 snRNPs that contain the U5–200kD protein but lack U5–100kD, suggesting that the U5–200kD protein could mediate U4/U6 duplex unwinding. Finally, U5–200kD was purified to homogeneity by glycerol gradient centrifugation of U5 snRNP proteins in the presence of sodium thiocyanate, followed by ion exchange chromatography. The RNA unwinding activity was found to reside exclusively with the U5–200kD DEXH-box protein. Our data raise the interesting possibility that this RNA helicase catalyzes unwinding of the U4/U6 RNA duplex in the spliceosome.

Polyadenylation of stable RNA precursors in vivo

Polyadenylation at the 3′ terminus has long been considered a specific feature of mRNA and a few other unstable RNA species. Here we show that stable RNAs in Escherichia coli can be polyadenylated as well. RNA molecules with poly(A) tails are the major products that accumulate for essentially all stable RNA precursors when RNA maturation is slowed because of the absence of processing exoribonucleases; poly(A) tails vary from one to seven residues in length. The polyadenylation process depends on the presence of poly(A) polymerase I. A stochastic competition between the exoribonucleases and poly(A) polymerase is proposed to explain the accumulation of polyadenylated RNAs. These data indicate that polyadenylation is not unique to mRNA, and its widespread occurrence suggests that it serves a more general function in RNA metabolism.

AtGRP7, a nuclear RNA-binding protein as a component of a circadian-regulated negative feedback loop in Arabidopsis thaliana

The endogenous clock that drives circadian rhythms is thought to communicate temporal information within the cell via cycling downstream transcripts. A transcript encoding a glycine-rich RNA-binding protein, Atgrp7, in Arabidopsis thaliana undergoes circadian oscillations with peak levels in the evening. The AtGRP7 protein also cycles with a time delay so that Atgrp7 transcript levels decline when the AtGRP7 protein accumulates to high levels. After AtGRP7 protein concentration has fallen to trough levels, Atgrp7 transcript starts to reaccumulate. Overexpression of AtGRP7 in transgenic Arabidopsis plants severely depresses cycling of the endogenous Atgrp7 transcript. These data establish both transcript and protein as components of a negative feedback circuit capable of generating a stable oscillation. AtGRP7 overexpression also depresses the oscillation of the circadian-regulated transcript encoding the related RNA-binding protein AtGRP8 but does not affect the oscillation of transcripts such as cab or catalase mRNAs. We propose that the AtGRP7 autoregulatory loop represents a “slave” oscillator in Arabidopsis that receives temporal information from a central “master” oscillator, conserves the rhythmicity by negative feedback, and transduces it to the output pathway by regulating a subset of clock-controlled transcripts.

Heterogeneous nuclear ribonucleoprotein A1 binds to the transcription-regulatory region of mouse hepatitis virus RNA

A cellular protein, previously described as p35/38, binds to the complementary (−)-strand of the leader RNA and intergenic (IG) sequence of mouse hepatitis virus (MHV) RNA. The extent of the binding of this protein to IG sites correlates with the efficiency of the subgenomic mRNA transcription from that IG site, suggesting that it is a requisite transcription factor. We have purified this protein and determined by partial peptide sequencing that it is heterogeneous nuclear ribonucleoprotein (hnRNP) A1, an abundant, primarily nuclear protein. hnRNP A1 shuttles between the nucleus and cytoplasm and plays a role in the regulation of alternative RNA splicing. The MHV(−)-strand leader and IG sequences conform to the consensus binding motifs of hnRNP A1. Recombinant hnRNP A1 bound to these two RNA regions in vitro in a sequence-specific manner. During MHV infection, hnRNP A1 relocalizes from the nucleus to the cytoplasm, where viral replication occurs. These data suggest that hnRNP A1 is a cellular factor that regulates the RNA-dependent RNA transcription of the virus.

Specific atomic groups and RNA helix geometry in acceptor stem recognition by a tRNA synthetase

Oligonucleotides that recapitulate the acceptor stems of tRNAs are substrates for aminoacylation by many tRNA synthetases in vitro, even though these substrates are missing the anticodon trinucleotides of the genetic code. In the case of tRNAAla a single acceptor stem G⋅U base pair at position 3·70 is essential, based on experiments where the wobble pair has been replaced by alternatives such as I⋅U, G⋅C, and A⋅U, among others. These experiments led to the conclusion that the minor-groove free 2-amino group (of guanosine) of the G⋅U wobble pair is essential for charging. Moreover, alanine-inserting tRNAs (amber suppressors) that replace G⋅U with mismatches such as G⋅A and C⋅A are partially active in vivo and can support growth of an Escherichia coli tRNAAla knockout strain, leading to the hypothesis that a helix irregularity and nucleotide functionalities are important for recognition. Herein we investigate the charging in vitro of oligonucleotide and full-length tRNA substrates that contain mismatches at the position of the G⋅U pair. Although most of these substrates have undetectable activity, G⋅A and C⋅A variants retain some activity, which is, nevertheless, reduced by at least 100-fold. Thus, the in vivo assays are much less sensitive to large changes in aminoacylation kinetic efficiency of 3·70 variants than is the in vitro assay system. Although these functional data do not clarify all of the details, it is now clear that specific atomic groups are substantially more important in determining kinetic efficiency than is a helical distortion. By implication, the activity of mutant tRNAs measured in the in vivo assays appears to be more dependent on factors other than aminoacylation kinetic efficiency.

Cell type-specific loss of atp6 RNA editing in cytoplasmic male sterile Sorghum bicolor

RNA editing and cytoplasmic male sterility are two important phenomena in higher plant mitochondria. To determine whether correlations might exist between the two, RNA editing in different tissues of Sorghum bicolor was compared employing reverse transcription–PCR and subsequent sequence analysis. In etiolated shoots, RNA editing of transcripts of plant mitochondrial atp6, atp9, nad3, nad4, and rps12 genes was identical among fertile or cytoplasmic male sterile plants. We then established a protocol for mitochondrial RNA isolation from plant anthers and pollen to include in these studies. Whereas RNA editing of atp9, nad3, nad4, and rps12 transcripts in anthers was similar to etiolated shoots, mitochondrial atp6 RNA editing was strongly reduced in anthers of the A3Tx398 male sterile line of S. bicolor. atp6 transcripts of wheat and selected plastid transcripts in S. bicolor showed normal RNA editing, indicating that loss of atp6 RNA editing is specific for cytoplasmic male sterility S. bicolor mitochondria. Restoration of fertility in F1 and F2 lines correlated with an increase in RNA editing of atp6 transcripts. Our data suggest that loss of atp6 RNA editing contributes to or causes cytoplasmic male sterility in S. bicolor. Further analysis of the mechanism of cell type-specific loss of atp6 RNA editing activity may advance our understanding of the mechanism of RNA editing.

Extensive purification of a putative RNA polymerase I holoenzyme from plants that accurately initiates rRNA gene transcription in vitro

RNA polymerase I (pol I) is a nuclear enzyme whose function is to transcribe the duplicated genes encoding the precursor of the three largest ribosomal RNAs. We report a cell-free system from broccoli (Brassica oleracea) inflorescence that supports promoter-dependent RNA pol I transcription in vitro. The transcription system was purified extensively by DEAE-Sepharose, Biorex 70, Sephacryl S300, and Mono Q chromatography. Activities required for pre-rRNA transcription copurified with the polymerase on all four columns, suggesting their association as a complex. Purified fractions programmed transcription initiation from the in vivo start site and utilized the same core promoter sequences required in vivo. The complex was not dissociated in 800 mM KCl and had a molecular mass of nearly 2 MDa based on gel filtration chromatography. The most highly purified fractions contain ≈30 polypeptides, two of which were identified immunologically as RNA polymerase subunits. These data suggest that the occurrence of a holoenzyme complex is probably not unique to the pol II system but may be a general feature of eukaryotic nuclear polymerases.

Nonequilibrium mechanism of transcription termination from observations of single RNA polymerase molecules

Cessation of transcription at specific terminator DNA sequences is used by viruses, bacteria, and eukaryotes to regulate the expression of downstream genes, but the mechanisms of transcription termination are poorly characterized. To elucidate the kinetic mechanism of termination at the intrinsic terminators of enteric bacteria, we observed, by using single-molecule light microscopy techniques, the behavior of surface-immobilized Escherichia coli RNA polymerase (RNAP) molecules in vitro. An RNAP molecule remains at a canonical intrinsic terminator for ≈64 s before releasing DNA, implying the formation of an elongation-incompetent (paused) intermediate by transcription complexes that terminate but not by those that read through the terminator. Analysis of pause lifetimes establishes a complete minimal mechanism of termination in which paused intermediate formation is both necessary and sufficient to induce release of RNAP at the terminator. The data suggest that intrinsic terminators function by a nonequilibrium process in which terminator effectiveness is determined by the relative rates of nucleotide addition and paused state entry by the transcription complex.

A curved RNA helix incorporating an internal loop with G·A and A·A non-Watson–Crick base pairing

The crystal structure of the RNA dodecamer 5′-GGCC(GAAA)GGCC-3′ has been determined from x-ray diffraction data to 2.3-Å resolution. In the crystal, these oligomers form double helices around twofold symmetry axes. Four consecutive non-Watson–Crick base pairs make up an internal loop in the middle of the duplex, including sheared G·A pairs and novel asymmetric A·A pairs. This internal loop sequence produces a significant curvature and narrowing of the double helix. The helix is curved by 34° from end to end and the diameter is narrowed by 24% in the internal loop. A Mn2+ ion is bound directly to the N7 of the first guanine in the Watson–Crick region following the internal loop and the phosphate of the preceding residue. This Mn2+ location corresponds to a metal binding site observed in the hammerhead catalytic RNA.

Mutation rates among RNA viruses

The rate of spontaneous mutation is a key parameter in modeling the genetic structure and evolution of populations. The impact of the accumulated load of mutations and the consequences of increasing the mutation rate are important in assessing the genetic health of populations. Mutation frequencies are among the more directly measurable population parameters, although the information needed to convert them into mutation rates is often lacking. A previous analysis of mutation rates in RNA viruses (specifically in riboviruses rather than retroviruses) was constrained by the quality and quantity of available measurements and by the lack of a specific theoretical framework for converting mutation frequencies into mutation rates in this group of organisms. Here, we describe a simple relation between ribovirus mutation frequencies and mutation rates, apply it to the best (albeit far from satisfactory) available data, and observe a central value for the mutation rate per genome per replication of μg ≈ 0.76. (The rate per round of cell infection is twice this value or about 1.5.) This value is so large, and ribovirus genomes are so informationally dense, that even a modest increase extinguishes the population.

Altered phosphorylation and intracellular distribution of a (CUG)n triplet repeat RNA-binding protein in patients with myotonic dystrophy and in myotonin protein kinase knockout mice

Myotonic dystrophy (DM) is associated with expansion of CTG repeats in the 3′-untranslated region of the myotonin protein kinase (DMPK) gene. The molecular mechanism whereby expansion of the (CUG)n repeats in the 3′-untranslated region of DMPK gene induces DM is unknown. We previously isolated a protein with specific binding to CUG repeat sequences (CUG-BP/hNab50) that possibly plays a role in mRNA processing and/or transport. Here we present evidence that the phosphorylation status and intracellular distribution of the RNA CUG-binding protein, identical to hNab50 protein (CUG-BP/hNab50), are altered in homozygous DM patient and that CUG-BP/hNab50 is a substrate for DMPK both in vivo and in vitro. Data from two biological systems with reduced levels of DMPK, homozygous DM patient and DMPK knockout mice, show that DMPK regulates both phosphorylation and intracellular localization of the CUG-BP/hNab50 protein. Decreased levels of DMPK observed in DM patients and DMPK knockout mice led to the elevation of the hypophosphorylated form of CUG-BP/hNab50. Nuclear concentration of the hypophosphorylated CUG-BP/hNab50 isoform is increased in DMPK knockout mice and in homozygous DM patient. DMPK also interacts with and phosphorylates CUG-BP/hNab50 protein in vitro. DMPK-mediated phosphorylation of CUG-BP/hNab50 results in dramatic reduction of the CUG-BP2, hypophosphorylated isoform, accumulation of which was observed in the nuclei of DMPK knockout mice. These data suggest a feedback mechanism whereby decreased levels of DMPK could alter phosphorylation status of CUG-BP/hNab50, thus facilitating nuclear localization of CUG-BP/hNab50. Our results suggest that DM pathophysiology could be, in part, a result of sequestration of CUG-BP/hNab50 and, in part, of lowered DMPK levels, which, in turn, affect processing and transport of specific subclass of mRNAs.

Functional topography of nascent RNA in elongation intermediates of RNA polymerase

To determine the dynamics of transcript extrusion from Escherichia coli RNA polymerase (RNAP), we used degradation of the RNA by RNases T1 and A in a series of consecutive elongation complexes (ECs). In intact ECs, even extremely high doses of the RNases were unable to cut the RNA closer than 14–16 nt from the 3′ end. Our results prove that all of the cuts detected within the 14-nt zone are derived from the EC that is denatured during inactivation of the RNases. The protected zone monotonously translocates along the RNA after addition of new nucleotides to the transcript. The upstream region of the RNA heading toward the 5′ end is cleaved and dissociated from the EC, with no effect on the stability and activity of the EC. Most of the current data suggest that an 8- to 10-nt RNA⋅DNA hybrid is formed in the EC. Here, we show that an 8- to 10-nt RNA obtained by truncating the RNase-generated products further with either GreB or pyrophosphate is sufficient for the high stability and activity of the EC. This result suggests that the transcript–RNAP interaction that is required for holding the EC together can be limited to the RNA region involved in the 8- to 10-nt RNA⋅DNA hybrid.