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.

Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy

It has long been assumed that HIV-1 evolution is best described by deterministic evolutionary models because of the large population size. Recently, however, it was suggested that the effective population size (Ne) may be rather small, thereby allowing chance to influence evolution, a situation best described by a stochastic evolutionary model. To gain experimental evidence supporting one of the evolutionary models, we investigated whether the development of resistance to the protease inhibitor ritonavir affected the evolution of the env gene. Sequential serum samples from five patients treated with ritonavir were used for analysis of the protease gene and the V3 domain of the env gene. Multiple reverse transcription–PCR products were cloned, sequenced, and used to construct phylogenetic trees and to calculate the genetic variation and Ne. Genotypic resistance to ritonavir developed in all five patients, but each patient displayed a unique combination of mutations, indicating a stochastic element in the development of ritonavir resistance. Furthermore, development of resistance induced clear bottleneck effects in the env gene. The mean intrasample genetic variation, which ranged from 1.2% to 5.7% before treatment, decreased significantly (P < 0.025) during treatment. In agreement with these findings, Ne was estimated to be very small (500–15,000) compared with the total HIV-1 RNA copy number. This study combines three independent observations, strong population bottlenecking, small Ne, and selection of different combinations of protease-resistance mutations, all of which indicate that HIV-1 evolution is best described by a stochastic evolutionary model.

Phylogeny of rice genomes with emphasis on origins of allotetraploid species

The rice genus, Oryza, which comprises 23 species and 9 recognized genome types, represents an enormous gene pool for genetic improvement of rice cultivars. Clarification of phylogenetic relationships of rice genomes is critical for effective utilization of the wild rice germ plasm. By generating and comparing two nuclear gene (Adh1 and Adh2) trees and a chloroplast gene (matK) tree of all rice species, phylogenetic relationships among the rice genomes were inferred. Origins of the allotetraploid species, which constitute more than one-third of rice species diversity, were reconstructed based on the Adh gene phylogenies. Genome types of the maternal parents of allotetraploid species were determined based on the matK gene tree. The phylogenetic reconstruction largely supports the previous recognition of rice genomes. It further revealed that the EE genome species is most closely related to the DD genome progenitor that gave rise to the CCDD genome. Three species of the CCDD genome may have originated through a single hybridization event, and their maternal parent had the CC genome. The BBCC genome species had different origins, and their maternal parents had either a BB or CC genome. An additional genome type, HHKK, was recognized for Oryza schlechteri and Porteresia coarctata, suggesting that P. coarctata is an Oryza species. The AA genome lineage, which contains cultivated rice, is a recently diverged and rapidly radiated lineage within the rice genus.

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.

The RDP-II (Ribosomal Database Project)

The Ribosomal Database Project (RDP-II), previously described by Maidak et al. [Nucleic Acids Res. (2000), 28, 173–174], continued during the past year to add new rRNA sequences to the aligned data and to improve the analysis commands. Release 8.0 (June 1, 2000) consisted of 16 277 aligned prokaryotic small subunit (SSU) rRNA sequences while the number of eukaryotic and mitochondrial SSU rRNA sequences in aligned form remained at 2055 and 1503, respectively. The number of prokaryotic SSU rRNA sequences more than doubled from the previous release 14 months earlier, and ~75% are longer than 899 bp. An RDP-II mirror site in Japan is now available (http://wdcm.nig.ac.jp/RDP/html/index.html). RDP-II provides aligned and annotated rRNA sequences, derived phylogenetic trees and taxonomic hierarchies, and analysis services through its WWW server (http://rdp.cme.msu.edu/). Analysis services include rRNA probe checking, approximate phylogenetic placement of user sequences, screening user sequences for possible chimeric rRNA sequences, automated alignment, production of similarity matrices and services to plan and analyze terminal restriction fragment polymorphism experiments. The RDP-II email address for questions and comments has been changed from curator@cme.msu.edu to rdpstaff@msu.edu.

PlantsP: a functional genomics database for plant phosphorylation

The PlantsP database is a curated database that combines information derived from sequences with experimental functional genomics information. PlantsP focuses on plant protein kinases and protein phosphatases. The database will specifically provide a resource for information on a collection of T-DNA insertion mutants (knockouts) in each protein kinase and phosphatase in Arabidopsis thaliana. PlantsP also provides a curated view of each protein that includes a comprehensive annotation of functionally related sequence motifs, sequence family definitions, alignments and phylogenetic trees, and descriptive information drawn directly from the literature. PlantsP is available at http://PlantsP.sdsc.edu.

Recursive partitioning for tumor classification with gene expression microarray data

Precise classification of tumors is critically important for cancer diagnosis and treatment. It is also a scientifically challenging task. Recently, efforts have been made to use gene expression profiles to improve the precision of classification, with limited success. Using a published data set for purposes of comparison, we introduce a methodology based on classification trees and demonstrate that it is significantly more accurate for discriminating among distinct colon cancer tissues than other statistical approaches used heretofore. In addition, competing classification trees are displayed, which suggest that different genes may coregulate colon cancers.

Detection of a major gene for resistance to fusiform rust disease in loblolly pine by genomic mapping.

Genomic mapping has been used to identify a region of the host genome that determines resistance to fusiform rust disease in loblolly pine where no discrete, simply inherited resistance factors had been previously found by conventional genetic analysis over four decades. A resistance locus, behaving as a single dominant gene, was mapped by association with genetic markers, even though the disease phenotype deviated from the expected Mendelian ratio. The complexity of forest pathosystems and the limitations of genetic analysis, based solely on phenotype, had led to an assumption that effective long-term disease resistance in trees should be polygenic. However, our data show that effective long-term resistance can be obtained from a single qualitative resistance gene, despite the presence of virulence in the pathogen population. Therefore, disease resistance in this endemic coevolved forest pathosystem is not exclusively polygenic. Genomic mapping now provides a powerful tool for characterizing the genetic basis of host pathogen interactions in forest trees and other undomesticated, organisms, where conventional genetic analysis often is limited or not feasible.

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.

Island colonization and evolution of the insular woody habit in Echium L. (Boraginaceae).

Numerous island-inhabiting species of predominantly herbaceous angiosperm genera are woody shrubs or trees. Such "insular woodiness" is strongly manifested in the genus Echium, in which the continental species of circummediterranean distribution are herbaceous, whereas endemic species of islands along the Atlantic coast of north Africa are woody perennial shrubs. The history of 37 Echium species was traced with 70 kb of noncoding DNA determined from both chloroplast and nuclear genomes. In all, 239 polymorphic positions with 137 informative sites, in addition to 27 informative indels, were found. Island-dwelling Echium species are shown to descend from herbaceous continental ancestors via a single island colonization event that occurred < 20 million years ago. Founding colonization appears to have taken place on the Canary Islands, from which the Madeira and Cape Verde archipelagos were invaded. Colonization of island habitats correlates with a recent origin of perennial woodiness from herbaceous habit and was furthermore accompanied by intense speciation, which brought forth remarkable diversity of forms among contemporary island endemics. We argue that the origin of insular woodiness involved response to counter-selection of inbreeding depression in founding island colonies.

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 and with the enumeration of 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.