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.

Matchings and phylogenetic trees

This paper presents a natural coordinate system for phylogenetic trees using a correspondence with the set of perfect matchings in the complete graph. This correspondence produces a distance between phylogenetic trees, and a way of enumerating all trees in a minimal step order. It is useful in randomized algorithms because it enables moves on the space of trees that make random optimization strategies “mix” quickly. It also promises a generalization to intermediary trees when data are not decisive as to their choice of tree, and a new way of constructing Bayesian priors on tree space.

Directed seed dispersal by bellbirds in a tropical cloud forest

A fundamental goal of plant population ecology is to understand the consequences for plant fitness of seed dispersal by animals. Theories of seed dispersal and tropical forest regeneration suggest that the advantages of seed dispersal for most plants are escape from seed predation near the parent tree and colonization of vacant sites, the locations of which are unpredictable in space and time. Some plants may gain in fitness as a fortuitous consequence of disperser behavior if certain species of dispersers nonrandomly place seeds in sites predictably favorable for seedling establishment. Such patterns of directed dispersal by vertebrates long have been suggested but never demonstrated for tropical forest trees. Here we report the pattern of seed distribution and 1-year seedling survival generated by five species of birds for a neotropical, shade-tolerant tree. Four of the species dispersed seeds to sites near the parent trees with microhabitat characteristics similar to those at random locations, whereas the fifth species, a bellbird, predictably dispersed seeds under song perches in canopy gaps. The pattern of seedling recruitment was bimodal, with a peak near parent trees and a second peak, corresponding to bellbird song perches, far (>40 m) from parent trees. Seedling survival was higher for seeds dispersed by bellbirds than by the other species, because of a reduction in seedling mortality by fungal pathogens in gaps. Thus, bellbirds play a significant role in seed dispersal by providing directed dispersal to favorable sites and therefore may influence plant recruitment patterns and species diversity in Neotropical forests.

Phytochrome D Acts in the Shade-Avoidance Syndrome in Arabidopsis by Controlling Elongation Growth and Flowering Time1

Shade avoidance in higher plants is regulated by the action of multiple phytochrome (phy) species that detect changes in the red/far-red ratio (R/FR) of incident light and initiate a redirection of growth and an acceleration of flowering. The phyB mutant of Arabidopsis is constitutively elongated and early flowering and displays attenuated responses to both reduced R/FR and end-of-day far-red light, conditions that induce strong shade-avoidance reactions in wild-type plants. This indicates that phyB plays an important role in the control of shade avoidance. In Arabidopsis phyB and phyD are the products of a recently duplicated gene and share approximately 80% identity. We investigated the role played by phyD in shade avoidance by analyzing the responses of phyD-deficient mutants. Compared with the monogenic phyB mutant, the phyB-phyD double mutant flowers early and has a smaller leaf area, phenotypes that are characteristic of shade avoidance. Furthermore, compared with the monogenic phyB mutant, the phyB-phyD double mutant shows a more attenuated response to a reduced R/FR for these responses. Compared with the phyA-phyB double mutant, the phyA-phyB-phyD triple mutant has elongated petioles and displays an enhanced elongation of internodes in response to end-of-day far-red light. These characteristics indicate that phyD acts in the shade-avoidance syndrome by controlling flowering time and leaf area and that phyC and/or phyE also play a role.

Visualization of Cavitated Vessels in Winter and Refilled Vessels in Spring in Diffuse-Porous Trees by Cryo-Scanning Electron Microscopy1

Xylem cavitation in winter and recovery from cavitation in the spring were visualized in two species of diffuse-porous trees, Betula platyphylla var. japonica Hara and Salix sachalinensis Fr. Schm., by cryo-scanning electron microscopy after freeze-fixation of living twigs. Water in the vessel lumina of the outer three annual rings of twigs of B. platyphylla var. japonica and of S. sachalinensis gradually disappeared during the period from January to March, an indication that cavitation occurs gradually in these species during the winter. In April, when no leaves had yet expanded, the lumina of most of the vessels of both species were filled with water. Many vessel lumina in twigs of both species were filled with water during the period from the subsequent growth season to the beginning of the next winter. These observations indicate that recovery in spring occurs before the onset of transpiration and that water transport through twigs occurs during the subsequent growing season. We found, moreover, that vessels repeat an annual cycle of winter cavitation and spring recovery from cavitation for several years until irreversible cavitation occurs.

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.

Twilight-zone and canopy shade induction of the Athb-2 homeobox gene in green plants.

We present evidence that a novel phytochrome (other than phytochromes A and B, PHYA and PHYB) operative in green plants regulates the "twilight-inducible" expression of a plant homeobox gene (Athb-2). Light regulation of the Athb-2 gene is unique in that it is not induced by red (R)-rich daylight or by the light-dark transition but is instead induced by changes in the ratio of R to far-red (FR) light. These changes, which normally occur at dawn and dusk (end-of-day FR), also occur during the daytime under the canopy (shade avoidance). By using pure light sources and phyA/phyB null mutants, we demonstrated that the induction of Athb-2 by changes in the R/FR ratio is mediated for the most part by a novel phytochrome operative in green plants. Furthermore, PHYB plays a negative role in repressing the accumulation of Athb-2 mRNA in the dark and a minor role in the FR response. The strict correlation of Athb-2 expression with FR-induced growth phenomena suggests a role for the Athb-2 gene in mediating cell elongation. This interpretation is supported by the finding that the Athb-2 gene is expressed at high levels in rapidly elongating etiolated seedlings. Furthermore, as either R or FR light inhibits cell elongation in etiolated tissues, they also down-regulate the expression of Athb-2 mRNA. Thus, these data support the notion that changes in light quality perceived by a novel phytochrome regulate plant development through the action of the Athb-2 homeobox gene.

Hyperplane arrangements, interval orders, and trees.

A hyperplane arrangement is a finite set of hyperplanes in a real affine space. An especially important arrangement is the braid arrangement, which is the set of all hyperplanes xi - xj = 1, 1 trees. For instance, the number of labeled interval orders that can be obtained from n intervals I1,..., In of generic lengths is counted. There is also discussed an arrangement due to N. Linial whose number of regions is the number of alternating (or intransitive) trees, as defined by Gelfand, Graev, and Postnikov [Gelfand, I. M., Graev, M. I., and Postnikov, A. (1995), preprint]. Finally, a refinement is given, related to counting labeled trees by number of inversions, of a result of Shi [Shi, J.-Y. (1986), Lecture Notes in Mathematics, no. 1179, Springer-Verlag] that a certain deformation of the braid arrangement has (n + 1)n-1 regions.

Calculating the probability of multitaxon evolutionary trees: bootstrappers Gambit.

The reconstruction of multitaxon trees from molecular sequences is confounded by the variety of algorithms and criteria used to evaluate trees, making it difficult to compare the results of different analyses. A global method of multitaxon phylogenetic reconstruction described here, Bootstrappers Gambit, can be used with any four-taxon algorithm, including distance, maximum likelihood, and parsimony methods. It incorporates a Bayesian-Jeffreys'-bootstrap analysis to provide a uniform probability-based criterion for comparing the results from diverse algorithms. To examine the usefulness of the method, the origin of the eukaryotes has been investigated by the analysis of ribosomal small subunit RNA sequences. Three common algorithms (paralinear distances, Jukes-Cantor distances, and Kimura distances) support the eocyte topology, whereas one (maximum parsimony) supports the archaebacterial topology, suggesting that the eocyte prokaryotes are the closest prokaryotic relatives of the eukaryotes.