Improved bootstrap confidence limits in large-scale phylogenies, with an example from Neo-Astragalus (Leguminosae)

Phylogenetic analyses of large data sets pose special challenges, including the apparent tendency for the bootstrap support for a clade to decline with increased taxon sampling of that clade. We document this decline in data sets with increasing numbers of taxa in Astragalus, the most species-rich angiosperm genus. Support for one subclade, Neo-Astragalus, declined monotonically with increased sampling of taxa inside Neo-Astragalus, irrespective of whether parsimony or neighbor-joining methods were used or of which particular heuristic search algorithm was used (although more stringent algorithms tended to yield higher support). Three possible explanations for this decline were examined, including (1) mistaken assignment of the most recent common ancestor of the taxon sample (and its bootstrap support) with the most recent common ancestor of the clade from which it was sampled; (2) computational limitations of heuristic search strategies; and (3) statistical bias in bootstrap proportions, especially that from random homoplasy distributed among taxa. The best explanation appears to be (3), although computational shortcomings (2) may explain some of the problem. The bootstrap proportion, as currently used in phylogenetic analysis, does not accurately capture the classical notion of confidence assessments on the null hypothesis of nonmonophyly, especially in large data sets. More accurate assessments of confidence as type I error levels (relying on iterated bootstrap methods) remove most of the monotonic decline in confidence with increasing numbers of taxa.

Error, bias, and long-branch attraction in data for two chloroplast photosystem genes in seed plants

Sequences of two chloroplast photosystem genes, psaA and psbB, together comprising about 3,500 bp, were obtained for all five major groups of extant seed plants and several outgroups among other vascular plants. Strongly supported, but significantly conflicting, phylogenetic signals were obtained in parsimony analyses from partitions of the data into first and second codon positions versus third positions. In the former, both genes agreed on a monophyletic gymnosperms, with Gnetales closely related to certain conifers. In the latter, Gnetales are inferred to be the sister group of all other seed plants, with gymnosperms paraphyletic. None of the data supported the modern "anthophyte hypothesis," which places Gnetales as the sister group of flowering plants. A series of simulation studies were undertaken to examine the error rate for parsimony inference. Three kinds of errors were examined: random error, systematic bias (both properties of finite data sets), and statistical inconsistency owing to long-branch attraction (an asymptotic property). Parsimony reconstructions were extremely biased for third-position data for psbB. Regardless of the true underlying tree, a tree in which Gnetales are sister to all other seed plants was likely to be reconstructed for these data. None of the combinations of genes or partitions permits the anthophyte tree to be reconstructed with high probability. Simulations of progressively larger data sets indicate the existence of long-branch attraction (statistical inconsistency) for third-position psbB data if either the anthophyte tree or the gymnosperm tree is correct. This is also true for the anthophyte tree using either psaA third positions or psbB first and second positions. A factor contributing to bias and inconsistency is extremely short branches at the base of the seed plant radiation, coupled with extremely high rates in Gnetales and nonseed plant outgroups.

Phylogenetic systematics of the tribe Millettieae (Leguminosae) based on chloroplast trnK/matK sequences and its implications for evolutionary patterns in Papilionoideae

Phylogenetic relationships in the tribe Millettieae and allies in the subfamily Papilionoideae (Leguminosae) were reconstructed from chloroplast trnK/matK sequences. Sixty-two accessions representing 57 traditionally recognized genera of Papilionoideae were sampled, including 27 samples from Millettieae. Phylogenies were constructed using maximum parsimony and are well resolved and supported by high bootstrap values. A well-supported "core Millettieae" clade is recognized, comprising the four large genera Millettia, Lonchocarpus, Derris, and Tephrosia. Several other small genera of Millettieae are not in the core Millettieae clade. Platycyamus is grouped with Phaseoleae (in part). Ostryocarpus, Austrosteenisia, and Dalbergiella are neither in the core Millettieae or Phaseoleae clade. These taxa, along with core Millettieae and Phaseoleae, form a monophyletic sister group to Indigofereae. Cyclolobium and Poecilanthe are close to Brongniartieae. Callerya and Wisteria belong to a large clade that includes all the legumes that lack the inverted repeat in their chloroplast genome, which confirms previous rbcL and phytochrome gene family phylogenies. The evolutionary history of four characters was examined in Millettieae and allies: the presence of canavanine, inflorescence types, the dehiscence of pods, and the presence of winged pods. trnK/matK sequence analysis suggests that the presence of a pseudoraceme or pseudopanicle and the accumulation of nonprotein amino acids are phylogenetically informative for Millettieae and allies with only a few exceptions.

Testing transposable elements as genetic drive mechanisms using Drosophila P element constructs as a model system

The use of transposable elements (TEs) as genetic drive mechanisms was explored using Drosophila melanogaster as a model system. Alternative strategies, employing autonomous and nonautonomous P element constructs were compared for their efficiency in driving the ry+ allele into populations homozygous for a ry- allele at the genomic rosy locus. Transformed flies were introduced at 1%, 5%, and 10% starting frequencies to establish a series of populations that were monitored over the course of 40 generations, using both phenotypic and molecular assays. The transposon-borne ry+ marker allele spread rapidly in almost all populations when introduced at 5% and 10% seed frequencies, but 1% introductions frequently failed to become established. A similar initial rapid increase in frequency of the ry+ transposon occurred in several control populations lacking a source of transposase. Constructs carrying ry+ markers also increased to moderate frequencies in the absence of selection on the marker. The results of Southern and in situ hybridization studies indicated a strong inverse relationship between the degree of conservation of construct integrity and transposition frequency. These finding have relevance to possible future applications of transposons as genetic drive mechanisms.

Persistent nuclear ribosomal DNA sequence polymorphism in the Amelanchier agamic complex (Rosaceae)

Individual plants of several Amelanchier taxa contain many polymorphic nucleotide sites in the internal transcribed spacers (ITS) of nuclear ribosomal DNA (nrDNA). This polymorphism is unusual because it is not recent in origin and thus has resisted homogenization by concerted evolution. Amelanchier ITS sequence polymorphism is hypothesized to be the result of gene flow between two major North American clades resolved by phylogenetic analysis of ITS sequences. Western North American species plus A. humilis and A. sanguinea of eastern North America form one clade (A), and the remaining eastern North American Amelanchier make up clade B. Five eastern North American taxa are polymorphic at many of the nucleotide sites where clades A and B have diverged and are thought to be of hybrid origin, with A. humilis or A. sanguinea as one parent and various members of clade B as the other parent. Morphological evidence suggests that A. humilis is one of the parents of one of the polymorphic taxa, a microspecies that we refer to informally as A. "erecta." Sequences of 21 cloned copies of the ITS1-5.8S gene-ITS2 region from one A. "erecta" individual are identical to A. humilis sequence or to the clade B consensus sequence, or they are apparent recombinants of A. humilis and clade B ITS repeats. Amelanchier "erecta" and another polymorphic taxon are suspected to be relatively old because both grow several hundred kilometers beyond the range of one of their parents. ITS sequence polymorphisms have apparently persisted in these two taxa perhaps because of polyploidy and/or agamospermy (asexual seed production), which are prevalent in the genus.

Paths and patterns: the biology and physics of swimming bacterial populations

The velocity distribution of swimming micro-organisms depends on directional cues supplied by the environment. Directional swimming within a bounded space results in the accumulation of organisms near one or more surfaces. Gravity, gradients of chemical concentration and illumination affect the motile behaviour of individual swimmers. Concentrated populations of organisms scatter and absorb light or consume molecules, such as oxygen. When supply is one-sided, consumption creates gradients; the presence of the population alters the intensity and the symmetry of the environmental cues. Patterns of cues interact dynamically with patterns of the consumer population. In suspensions, spatial variations in the concentration of organisms are equivalent to variations of mean mass density of the fluid. When organisms accumulate in one region whilst moving away from another region, the force of gravity causes convection that translocates both organisms and dissolved substances. The geometry of the resulting concentration-convection patterns has features that are remarkably reproducible. Of interest for biology are (1) the long-range organisation achieved by organisms that do not communicate, and (2) that the entire system, consisting of fluid, cells, directional supply of consumables, boundaries and gravity, generates a dynamic that improves the organisms' habitat by enhancing transport and mixing. Velocity distributions of the bacterium Bacillus subtilis have been measured within the milieu of the spatially and temporally varying oxygen concentration which they themselves create. These distributions of swimming speed and direction are the fundamental ingredients required for a quantitative mathematical treatment of the patterns. The quantitative measurement of swimming behaviour also contributes to our understanding of aerotaxis of individual cells.

Regulation of the SOS response in Bacillus subtilis: evidence for a LexA repressor homolog

The inducible SOS response for DNA repair and mutagenesis in the bacterium Bacillus subtilis resembles the extensively characterized SOS system of Escherichia coli. In this report, we demonstrate that the cellular repressor of the E. coli SOS system, the LexA protein, is specifically cleaved in B. subtilis following exposure of the cells to DNA-damaging treatments that induce the SOS response. The in vivo cleavage of LexA is dependent upon the functions of the E. coli RecA protein homolog in B. subtilis (B. subtilis RecA) and results in the same two cleavage fragments as produced in E. coli cells following the induction of the SOS response. We also show that a mutant form of the E. coli RecA protein (RecA430) can partially substitute for the nonfunctional cellular RecA protein in the B. subtilis recA4 mutant, in a manner consistent with its known activities and deficiencies in E. coli. RecA430 protein, which has impaired repressor cleaving (LexA, UmuD, and bacteriophage lambda cI) functions in E.coli, partially restores genetic exchange to B. subtilis recA4 strains but, unlike wild-type E. coli RecA protein, is not capable of inducing SOS functions (expression of DNA damage-inducible [din::Tn917-lacZ] operons or RecA synthesis) in B. subtilis in response to DNA-damaging agents or those functions that normally accompany the development of physiological competence. Our results provide support for the existence of a cellular repressor in B. subtilis that is functionally homologous to the E. coli LexA repressor and suggest that the mechanism by which B. subtilis RecA protein (like RecA of E. coli) becomes activated to promote the induction of the SOS response is also conserved.

DNA repair and the evolution of transformation in Bacillus subtilis. II. Role of inducible repair

In Bacillus subtilis, DNA repair and recombination are intimately associated with competence, the physiological state in which the bacterium can bind, take up and recombine exogenous DNA. Previously, we have shown that the homologous DNA transformation rate (ratio of transformants to total cells) increases with increasing UV dosage if cells are transformed after exposure to UV radiation (UV-DNA), whereas the transformation rate decreases if cells are transformed before exposure to UV (DNA-UV). In this report, by using different DNA repair-deficient mutants, we show that the greater increase in transformation rate in UV-DNA experiments than in DNA-UV experiments does not depend upon excision repair or inducible SOS-like repair, although certain quantitative aspects of the response do depend upon these repair systems. We also show that there is no increase in the transformation rate in a UV-DNA experiment when repair and recombination proficient cells are transformed with nonhomologous plasmid DNA, although the results in a DNA-UV experiment are essentially unchanged by using plasmid DNA. We have used din operon fusions as a sensitive means of assaying for the expression of genes under the control of the SOS-like regulon in both competent and noncompetent cell subpopulations as a consequence of competence development and our subsequent experimental treatments. Results indicate that the SOS-like system is induced in both competent and noncompetent subpopulations in our treatments and so should not be a major factor in the differential response in transformation rate observed in UV-DNA and DNA-UV treatments. These results provide further support to the hypothesis that the evolutionary function of competence is to bring DNA into the cell for use as template in the repair of DNA damage.

DNA repair and the evolution of transformation in the bacterium Bacillus subtilis

The purpose of the work reported here is to test the hypothesis that natural genetic transformation in the bacterium Bacillus subtilis has evolved as a DNA repair system. Specifically, tests were made to determine whether transformation functions to provide DNA template for the bacteria] cell to use in recombinational repair. The survivorship and the homologous transformation rate as a function of dose of ultraviolet irradiation (UV) was studied in two experimental treatments, in which cells were either transformed before (DNA-UV), or after (UV-DNA), treatment with UV. The results show that there is a qualitative difference in the relationship between the survival of transformed cells (sexual cells) and total cells (primarily asexual cells) in the two treatments. As predicted by the repair hypothesis, in the UV-DNA treatment, transformed cells had greater average survivorship than total cells, while in the DNA-UV treatment this relationship was reversed. There was also a consistent and qualitative difference between the UV-DNA and DNA-UV treatments in the relationship between the homologous transformation rate (transformed cells/total cells) and UV dosage. As predicted by the repair hypothesis, the homologous transformation rate increases with UV dose in the UV-DNA experiments but decreases with UV dose in the DNA-UV treatments. However, the transformation rate for plasmid DNA does not increase in a UV-DNA treatment. These results support the DNA repair hypothesis for the evolution of transformation in particular, and sex generally.

Inhibition of DNA methyltransferases in vitro by benzo[a]pyrene diol epoxide-modified substrates

Covalent adducts formed from the ultimate carcinogen 7 beta,8 alpha-dihydroxy-9 alpha, 10 alpha-epoxy-7,8,9,10-tetrahydrobenzo[ a]pyrene inhibit the enzyme-catalyzed transfer of methyl groups from S-adenosylmethionine to cytosine residues in DNA. Two DNA methyltransferase enzymes, isolated from the bacterium Haemophilus and mouse spleen nuclei, were tested for their ability to methylate carcinogen-modified substrates in vitro. These model enzymes possess the known methylation activities found in mammalian cells, de novo, and maintenance methylation of CpG-containing nucleotide sequences. The in vitro alkylation of DNA substrates by the carcinogen effectively decreases the methyltransferase reaction of both enzymes in a manner that is directly dependent upon the level of covalent modification of the DNA. Inhibition of de novo methylation activity can be detected at very low levels of carcinogen modification, 1 hydrocarbon residue per 20,000-40,000 nucleotides. Adduct levels in this range are capable of initiating transformation. Both enzymes are inactivated by direct reaction with the carcinogen in the absence of DNA. We also find that carcinogen adducts are capable of inhibiting DNA methylation at CpG sites removed from the primary lesion. These results support the proposal that carcinogen-induced DNA damage can cause alterations in methylation patterns that may eventually lead to heritable changes in gene expression.