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.

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.

Female preference for swords in Xiphophorus helleri reflects a bias for large apparent size

Swordtail fish (Poeciliidae: genus Xiphophorus) are a paradigmatic case of sexual selection by sensory exploitation. Female preference for males with a conspicuous “sword” ornament is ancestral, suggesting that male morphology has evolved in response to a preexisting bias. The perceptual mechanisms underlying female mate choice have not been identified, complicating efforts to understand the selection pressures acting on ornament design. We consider two alternative models of receiver behavior, each consistent with previous results. Females could respond either to specific characteristics of the sword or to more general cues, such as the apparent size of potential mates. We showed female swordtails a series of computer-altered video sequences depicting a courting male. Footage of an intact male was preferred strongly to otherwise identical sequences in which portions of the sword had been deleted selectively, but a disembodied courting sword was less attractive than an intact male. There was no difference between responses to an isolated sword and to a swordless male of comparable length, or between an isolated sword and a homogenous background. Female preference for a sworded male was abolished by enlarging the image of a swordless male to compensate for the reduction in length caused by removing the ornament. This pattern of results is consistent with mate choice being mediated by a general preference for large males rather than by specific characters. Similar processes may account for the evolution of exaggerated traits in other systems.

Crossover isomer bias is the primary sequence-dependent property of immobilized Holliday junctions

Recombination of genes is essential to the evolution of genetic diversity, the segregation of chromosomes during cell division, and certain DNA repair processes. The Holliday junction, a four-arm, four-strand branched DNA crossover structure, is formed as a transient intermediate during genetic recombination and repair processes in the cell. The recognition and subsequent resolution of Holliday junctions into parental or recombined products appear to be critically dependent on their three-dimensional structure. Complementary NMR and time-resolved fluorescence resonance energy transfer experiments on immobilized four-arm DNA junctions reported here indicate that the Holliday junction cannot be viewed as a static structure but rather as an equilibrium mixture of two conformational isomers. Furthermore, the distribution between the two possible crossover isomers was found to depend on the sequence in a manner that was not anticipated on the basis of previous low-resolution experiments.

Protein folding kinetics exhibit an Arrhenius temperature dependence when corrected for the temperature dependence of protein stability

The anomalous temperature dependence of protein folding has received considerable attention. Here we show that the temperature dependence of the folding of protein L becomes extremely simple when the effects of temperature on protein stability are corrected for; the logarithm of the folding rate is a linear function of 1/T on constant stability contours in the temperature–denaturant plane. This convincingly demonstrates that the anomalous temperature dependence of folding derives from the temperature dependence of the interactions that stabilize proteins, rather than from the super Arrhenius temperature dependence predicted for the configurational diffusion constant on a rough energy landscape. However, because of the limited temperature range accessible to experiment, the results do not rule out models with higher order temperature dependences. The significance of the slope of the stability-corrected Arrhenius plots is discussed.

The atypical codon usage of the plant psbA gene may be the remnant of an ancestral bias

The psbA gene of the chloroplast genome has a codon usage that is unusual for plant chloroplast genes. In the present study the evolutionary status of this codon usage is tested by reconstructing putative ancestral psbA sequences to determine the pattern of change in codon bias during angiosperm divergence. It is shown that the codon biases of the ancestral genes are much stronger than all extant flowering plant psbA genes. This is related to previous work that demonstrated a significant increase in synonymous substitution in psbA relative to other chloroplast genes. It is suggested, based on the two lines of evidence, that the codon bias of this gene currently is not being maintained by selection. Rather, the atypical codon bias simply may be a remnant of an ancestral codon bias that now is being degraded by the mutation bias of the chloroplast genome, in other words, that the psbA gene is not at equilibrium. A model for the evolution of selective pressure on the codon usage of plant chloroplast genes is discussed.

Prolonged production of NADPH oxidase-corrected granulocytes after gene therapy of chronic granulomatous disease

Little is known about the potential for engraftment of autologous hematopoietic stem cells in human adults not subjected to myeloablative conditioning regimens. Five adult patients with the p47phox deficiency form of chronic granulomatous disease received intravenous infusions of autologous CD34+ peripheral blood stem cells (PBSCs) that had been transduced ex vivo with a recombinant retrovirus encoding normal p47phox. Although marrow conditioning was not given, functionally corrected granulocytes were detectable in peripheral blood of all five patients. Peak correction occurred 3–6 weeks after infusion and ranged from 0.004 to 0.05% of total peripheral blood granulocytes. Corrected cells were detectable for as long as 6 months after infusion in some individuals. Thus, prolonged engraftment of autologous PBSCs and continued expression of the transduced gene can occur in adults without conditioning. This trial also piloted the use of animal protein-free medium and a blood-bank-compatible closed system of gas-permeable plastic containers for culture and transduction of the PBSCs. These features enhance the safety of PBSCs directed gene therapy.

Sources of selection bias in evaluating social programs: An interpretation of conventional measures and evidence on the effectiveness of matching as a program evaluation method

This paper decomposes the conventional measure of selection bias in observational studies into three components. The first two components are due to differences in the distributions of characteristics between participant and nonparticipant (comparison) group members: the first arises from differences in the supports, and the second from differences in densities over the region of common support. The third component arises from selection bias precisely defined. Using data from a recent social experiment, we find that the component due to selection bias, precisely defined, is smaller than the first two components. However, selection bias still represents a substantial fraction of the experimental impact estimate. The empirical performance of matching methods of program evaluation is also examined. We find that matching based on the propensity score eliminates some but not all of the measured selection bias, with the remaining bias still a substantial fraction of the estimated impact. We find that the support of the distribution of propensity scores for the comparison group is typically only a small portion of the support for the participant group. For values outside the common support, it is impossible to reliably estimate the effect of program participation using matching methods. If the impact of participation depends on the propensity score, as we find in our data, the failure of the common support condition severely limits matching compared with random assignment as an evaluation estimator.

How Escherichia coli can bias the results of molecular cloning: Preferential selection of defective genomes of hepatitis C virus during the cloning procedure

Cloned PCR products containing hepatitis C virus (HCV) genomic fragments have been used for analyses of HCV genomic heterogeneity and protein expression. These studies assume that the clones derived are representative of the entire virus population and that subsets are not inadvertently selected. The aim of the present study was to express HCV structural proteins. However, we found that there was a strong cloning selection for defective genomes and that most clones generated initially were incapable of expressing the HCV proteins. The HCV structural region (C-E1-E2-p7) was directly amplified by long reverse transcription–PCR from the plasma of an HCV-infected patient or from a control plasmid containing a viable full-length cDNA of HCV derived from the same patient but cloned in a different vector. The PCR products were cloned into a mammalian expression vector, amplified in Escherichia coli, and tested for their ability to produce HCV structural proteins. Twenty randomly picked clones derived from the HCV-infected patient all contained nucleotide mutations leading to absence or truncation of the expected HCV products. Of 25 clones derived from the control plasmid, only 8% were fully functional for polyprotein synthesis. The insertion of extra nucleotides in the region just upstream of the start codon of the HCV insert led to a statistically significant increase in the number of fully functional clones derived from the patient (42%) and from the control plasmid (72–92%). Nonrandom selection of clones during the cloning procedure has enormous implications for the study of viral heterogeneity, because it can produce a false spectrum of genomic diversity. It can also be an impediment to the construction of infectious viral clones.

Evolution of codon usage bias in Drosophila

We first review what is known about patterns of codon usage bias in Drosophila and make the following points: (i) Drosophila genes are as biased or more biased than those in microorganisms. (ii) The level of bias of genes and even the particular pattern of codon bias can remain phylogenetically invariant for very long periods of evolution. (iii) However, some genes, even very tightly linked genes, can change very greatly in codon bias across species. (iv) Generally G and especially C are favored at synonymous sites in biased genes. (v) With the exception of aspartic acid, all amino acids contribute significantly and about equally to the codon usage bias of a gene. (vi) While most individual amino acids that can use G or C at synonymous sites display a preference for C, there are exceptions: valine and leucine, which prefer G. (vii) Finally, smaller genes tend to be more biased than longer genes. We then examine possible causes of these patterns and discount mutation bias on three bases: there is little evidence of regional mutation bias in Drosophila, mutation bias is likely toward A+T (the opposite of codon usage bias), and not all amino acids display the preference for the same nucleotide in the wobble position. Two lines of evidence support a selection hypothesis based on tRNA pools: highly biased genes tend to be highly and/or rapidly expressed, and the preferred codons in highly biased genes optimally bind the most abundant isoaccepting tRNAs. Finally, we examine the effect of bias on DNA evolution and confirm that genes with high codon usage bias have lower rates of synonymous substitution between species than do genes with low codon usage bias. Surprisingly, we find that genes with higher codon usage bias display higher levels of intraspecific synonymous polymorphism. This may be due to opposing effects of recombination.

Effects of passage history and sampling bias on phylogenetic reconstruction of human influenza A evolution

In this paper we determine the extent to which host-mediated mutations and a known sampling bias affect evolutionary studies of human influenza A. Previous phylogenetic reconstruction of influenza A (H3N2) evolution using the hemagglutinin gene revealed an excess of nonsilent substitutions assigned to the terminal branches of the tree. We investigate two hypotheses to explain this observation. The first hypothesis is that the excess reflects mutations that were either not present or were at low frequency in the viral sample isolated from its human host, and that these mutations increased in frequency during passage of the virus in embryonated eggs. A set of 22 codons known to undergo such “host-mediated” mutations showed a significant excess of mutations assigned to branches attaching sequences from egg-cultured (as opposed to cell-cultured) isolates to the tree. Our second hypothesis is that the remaining excess results from sampling bias. Influenza surveillance is purposefully biased toward sequencing antigenically dissimilar strains in an effort to identify new variants that may signal the need to update the vaccine. This bias produces an excess of mutations assigned to terminal branches simply because an isolate with no close relatives is by definition attached to the tree by a relatively long branch. Simulations show that the magnitude of excess mutations we observed in the hemagglutinin tree is consistent with expectations based on our sampling protocol. Sampling bias does not affect inferences about evolution drawn from phylogenetic analyses. However, if possible, the excess caused by host-mediated mutations should be removed from studies of the evolution of influenza viruses as they replicate in their human hosts.

An essential role for Bruton's [corrected] tyrosine kinase in the regulation of B-cell apoptosis.

Mutations of the Bruton's tyrosine kinase (btk) gene cause X-linked agammaglobulinemia (XLA) in humans and X-linked immune deficiency (Xid) in mice. To establish the BTK role in B-cell activation we examined the responses of wild-type and Xid B cells to stimulation through surface IgM and CD40, the transducers of thymus independent-type 2 and thymus-dependent activation, respectively. Wild-type BTK was necessary for proliferation induced by soluble anti-IgM (a prototype for thymus independent-type 2 antigen), but not for responses to soluble CD40 ligand (CD40L, the B-cell activating ligand expressed on T-helper cells). In the absence of wild-type BTK, B cells underwent apoptotic death after stimulation with anti-IgM. In the presence of wild-type but not mutated BTK, anti-IgM stimulation reduced apoptotic cell death. In contrast, CD40L increased viability of both wild-type and Xid B cells. Importantly, viability after stimulation correlated with the induced expression of bcl-XL. In fresh ex vivo small resting B cells from wild-type mice there was only barely detectable bcl-XL protein, but there was more in the larger, low-density ("activated") splenic B cells and peritoneal B cells. In vitro Bcl-XL induction following ligation of sIgM-required BTK, was cyclosporin A (CsA)-sensitive and dependent on extracellular Ca2+. CD40-mediated induction of bcl-x required neither wild-type BTK nor extracellular Ca2+ and was insensitive to CsA. These results indicate that BTK lies upstream of bcl-XL in the sIgM but not the CD40 activation pathway. bcl-XL is the first induced protein to be placed downstream of BTK.

An Escherichia coli chromosomal "addiction module" regulated by guanosine [corrected] 3',5'-bispyrophosphate: a model for programmed bacterial cell death.

"Addiction modules" consist of two genes. In most of them the product of one is long lived and toxic while the product of the second is short lived and antagonizes the toxic effect; so far, they have been described mainly in a number of prokaryotic extrachromosomal elements responsible for the postsegregational killing effect. Here we show that the chromosomal genes mazE and mazF, located in the Escherichia coli rel operon, have all of the properties required for an addiction module. Furthermore, the expression of mazEF is regulated by the cellular level of guanosine [corrected] 3',5'-bispyrophosphate, the product of the RelA protein under amino acid starvation. These properties suggest that the mazEF system may be responsible for programmed cell death in E. coli and thus may have a role in the physiology of starvation.

Neighboring base composition and transversion/transition bias in a comparison of rice and maize chloroplast noncoding regions.

The correspondence between the transversion/transition ratio and the neighboring base composition in chloroplast DNA is examined. For 18 noncoding regions of the chloroplast genome, alignments between rice (Oryza sativa) and maize (Zea mays) were generated by two different methods. Difficulties of aligning noncoding DNA are discussed, and the alignments are analyzed in a manner that reduces alignment artifacts. Sequence divergence is < 10%, so multiple substitutions at a site are assumed to be rare. Observed substitutions were analyzed with respect to the A+T content of the two immediately flanking bases. It is shown that as this content increases, the proportion of transversions also increases. When both the 5'- and 3'-flanking nucleotides are G or C (A+T content of 0), only 25% of the observed substitutions are transversions. However, when both the 5'- and 3'-flanking nucleotides are A or T (A+T content of 2), 57% of the observed substitutions are transversions. Therefore, the influence of flanking base composition on substitutions, previously reported for a single noncoding region, is a general feature of the chloroplast genome.