Predicting the Impact of Describing New Species on Phylogenetic Patterns

Although our inventory of Earth’s biodiversity remains incomplete, we still require analyses using the Tree of Life to understand evolutionary and ecological patterns. Because incomplete sampling may bias our inferences, we must evaluate how future additions of newly discovered species might impact analyses performed today. We describe an approach that uses taxonomic history and phylogenetic trees to characterize the impact of past species discoveries on phylogenetic knowledge using patterns of branch-length variation, tree shape, and phylogenetic diversity. This provides a framework for assessing the relative completeness of taxonomic knowledge of lineages within a phylogeny. To demonstrate this approach, we use recent large phylogenies for amphibians, reptiles, flowering plants, and invertebrates. Well-known clades exhibit a decline in the mean and range of branch lengths that are added each year as new species are described. With increased taxonomic knowledge over time, deep lineages of well-known clades become known such that most recently described new species are added close to the tips of the tree, reflecting changing tree shape over the course of taxonomic history. The same analyses reveal other clades to be candidates for future discoveries that could dramatically impact our phylogenetic knowledge. Our work reveals that species are often added non-randomly to the phylogeny over multiyear time-scales in a predictable pattern of taxonomic maturation. Our results suggest that we can make informed predictions about how new species will be added across the phylogeny of a given clade, thus providing a framework for accommodating unsampled undescribed species in evolutionary analyses.

Accelerated diversification correlated with functional traits shapes extant diversity of the early divergent angiosperm family Annonaceae

A major goal of phylogenetic systematics is to understand both the patterns of diversification and the processes by which these patterns are formed. Few studies have focused on the ancient, species-rich Magnoliales clade and its diversification pattern. Within Magnoliales, the pantropically distributed Annonaceae are by far the most genus-rich and species-rich family-level clade, with c. 110 genera and c. 2,400 species. We investigated the diversification patterns across Annonaceae and identified traits that show varied associations with diversification rates using a time-calibrated phylogeny of 835 species (34.6% sampling) and 11,211 aligned bases from eight regions of the plastid genome (rbcL, matK, ndhF, psbA-trnH, trnL-F, atpB-rbcL, trnS-G, and ycf1). Twelve rate shifts were identified using BAMM: in Annona, Artabotrys, Asimina, Drepananthus, Duguetia, Goniothalamus, Guatteria, Uvaria, Xylopia, the tribes Miliuseae and Malmeeae, and the Desmos-Dasymaschalon-Friesodielsia-Monanthotaxis clade. TurboMEDUSA and method-of-moments estimator analyses showed largely congruent results. A positive relationship between species richness and diversification rate is revealed using PGLS. Our results show that the high species richness in Annonaceae is likely the result of recent increased diversification rather than the steady accumulation of species via the 'museum model'. We further explore the possible role of selected traits (habit, pollinator trapping, floral sex expression, pollen dispersal unit, anther septation, and seed dispersal unit) in shaping diversification patterns, based on inferences of BiSSE, MuSSE, HiSSE, and FiSSE analyses. Our results suggest that the liana habit, the presence of circadian pollinator trapping, androdioecy, and the dispersal of seeds as single-seeded monocarp fragments are closely correlated with higher diversification rates; pollen aggregation and anther septation, in contrast, are associated with lower diversification rates.

The Evolution of HD2 Proteins in Green Plants

In eukaryotes, protein deacetylation is carried out by two well-conserved histone deacetylase (HDAC) families: RPD3/HDA1 and SIR2. Intriguingly, model plants such as Arabidopsis express an additional plant-specific HDAC family, termed type-2 HDACs (HD2s). Transcriptomic analyses from more than 1300 green plants generated by the 1000 plants (1KP) consortium showed that HD2s appeared early in green plant evolution, the first members being detected in several streptophyte green alga. The HD2 family has expanded via several rounds of successive duplication; members are expressed in all major green plant clades. Interestingly, angiosperm species express new HD2 genes devoid of a zinc-finger domain, one of the main structural features of HD2s. These variants may have been associated with the origin and/or the biology of the ovule/seed.

Patterns of chromosomal variation in natural populations of the neoallotetraploid Tragopogon mirus (Asteraceae)

Cytological studies have shown many newly formed allopolyploids (neoallopolyploids) exhibit chromosomal variation as a result of meiotic irregularities, but few naturally occurring neoallopolyploids have been examined. Little is known about how long chromosomal variation may persist and how it might influence the establishment and evolution of allopolyploids in nature. In this study we assess chromosomal composition in a natural neoallotetraploid, Tragopogon mirus, and compare it with T. miscellus, which is an allotetraploid of similar age (~40 generations old). We also assess whether parental gene losses in T. mirus correlate with entire or partial chromosome losses. Of 37 T. mirus individuals that were karyotyped, 23 (62%) were chromosomally additive of the parents, whereas the remaining 14 individuals (38%) had aneuploid compositions. The proportion of additive versus aneuploid individuals differed from that found previously in T. miscellus, in which aneuploidy was more common (69%; Fisher's exact test, P=0.0033). Deviations from parental chromosome additivity within T. mirus individuals also did not reach the levels observed in T. miscellus, but similar compensated changes were observed. The loss of T. dubius-derived genes in two T. mirus individuals did not correlate with any chromosomal changes, indicating a role for smaller-scale genetic alterations. Overall, these data for T. mirus provide a second example of prolonged chromosomal instability in natural neoallopolyploid populations.

Genetic evidence suggests that homosporous ferns with high chromosome numbers are diploid

Homosporous ferns have usually been considered highly polyploid because they have high chromosome numbers (average n = 57.05). In angiosperms, species with chromosome numbers higher than n = 14 generally have more isozymes than those with lower numbers, consistent with their polyploidy. By extrapolation, homosporous ferns would be expected to have many isozymes. However, ongoing surveys indicate that within fern genera, species having the lowest chromosome numbers have the number of isozymes considered typical of diploid seed plants. Only species above these lowest numbers have additional isozymes. Therefore, homosporous ferns either have gone through repeated cycles of polyploidy and gene silencing or were initiated with relatively high chromosome numbers. The latter possibility represents a radical departure from currently advocated hypotheses of fern evolution and suggests that there may be fundamental differences between the genomes of homosporous ferns and those of higher plants. These hypotheses can be tested by genetic, karyological, and molecular techniques.

On the origins of species: does evolution repeat itself in polyploid populations of independent origin?

Multiple origins of the same polyploid species pose the question: Does evolution repeat itself in these independently formed lineages? Tragopogon is a unique evolutionary model for the study of recent and recurrent allopolyploidy. The allotetraploids T. mirus (T. dubius x T. porrifolius) and T. miscellus (T. dubius x T. pratensis) formed repeatedly following the introduction of three diploids to the United States. Concerted evolution has consistently occurred in the same direction (resulting in loss of T. dubius rDNA copies). Both allotetraploids exhibit homeolog loss, with the same genes consistently showing loss, and homeologs of T. dubius preferentially lost in both allotetraploids. We have also documented repeated patterns of tissue-specific silencing in multiple populations of T. miscellus. Hence, some aspects of genome evolution may be "hardwired," although the general pattern of loss is stochastic within any given population. On the basis of the study of F(1) hybrids and synthetics, duplicate gene loss and silencing do not occur immediately following hybridization or polyploidization, but gradually and haphazardly. Genomic approaches permit analysis of hundreds of loci to assess the frequency of homeolog loss and changes in gene expression. This methodology is particularly promising for groups such as Tragopogon for which limited genetic and genomic resources are available.

Tomato spotted wilt virus Found in Five Species of the genus Tragopogon in a Florida Greenhouse

An obviously unhealthy plant identified as Tragopogon mirus Ownbey (remarkable goatsbeard) was sent for diagnosis to the Division of Plant Industry (DPI), Gainesville, FL in May of 2008. T. mirus is a recently formed allotetraploid that has T. dubius Scop. and T. porrifolius L. (goatsbeard or salsify) as parents. The parents (family Asteraceae) are diploid and originate from Eurasia. They were introduced to the northwest United States in the early 1900s. The allotetraploid T. mirus, which does not occur in Eurasia, was discovered in 1949 and named in 1950. It has been found in the northwest states of Washington and Idaho. It has also been found in Arizona (4). The plant sent to the DPI was grown in a greenhouse for research purposes at the Botany Department of the University of Florida (Alachua County). Symptoms exhibited on the leaves included mottling, chlorotic and necrotic spots, and mild distortion. Epidermal leaf strips from a mottled leaf were stained with the Orange-Green protein stain and Azure A nucleic acid stain (1). With a light microscope, granular inclusions typical for Tomato spotted wilt virus (TSWV) (1) were seen in leaf strips from both stains. The remainder of the leaf was ground in buffer and tested serologically for TSWV by TSWV-specific ImmunoStrips (Agdia, Elkhart, IN). The ImmunoStrip was positive for the presence of TSWV. This test was confirmed by double-antibody sandwich-ELISA using antiserum and conjugate for TSWV (Agdia). Further serological testing of other Tragopogon species with similar symptoms growing in the same greenhouse revealed that T. miscellus (another recently formed allotetraploid found in the northwestern United States; parents T. dubius and T. pratensis), T. dubius, T. porrifolius, and T. pratensis were also infected with TSWV. Total RNA was extracted from symptomatic leaves of T. mirus, T. dubius, T. porrifolius, and T. miscellus. Reverse transcription-PCR was performed with universal tospovirus primers BR60 and BR65 that amplify part of the nucleocapsid protein gene (2). Target amplicons of 454 bp were produced for all four samples. The PCR product from T. porrifolius was cloned and sequenced. The resulting sequence (GenBank Accession No. FJ655913) shows high homology, 98%, to several isolates of the Tomato spotted wilt virus deposited in the GenBank (Accession Nos. AY870391, AY744477, and AF020659). T. porrifolius has been reported to be naturally infected with TSWV in Italy (3); however, to our knowledge, this is the first report of this virus in the allotetraploids T. mirus and T. miscellus and in the diploids T. dubius and T. pratensis. This report adds five new Asteraceae weeds to the list of possible reservoirs of TSWV in the United States. References: (1) J. R. Edwardson and R. G. Christie. Univ. Fla. Inst. Food Agric. Sci. Bull. 894. 1996. (2) M. Eiras et al. Fitopatol. Bras. 26:170, 2001. (3) G. Parrella et al. J. Plant Pathol. 85:227. 2003. (4) D. E. Soltis et al. Biol. J. Linn. Soc. 82:2004.

Gene loss and silencing in Tragopogon miscellus (Asteraceae): comparison of natural and synthetic allotetraploids

Whole-genome duplication (polyploidisation) is a widespread mechanism of speciation in plants. Over time, polyploid genomes tend towards a more diploid-like state, through downsizing and loss of duplicated genes (homoeologues), but relatively little is known about the timing of gene loss during polyploid formation and stabilisation. Several studies have also shown gene transcription to be affected by polyploidisation. Here, we examine patterns of gene loss in 10 sets of homoeologues in five natural populations of the allotetraploid Tragopogon miscellus that arose within the past 80 years following independent whole-genome duplication events. We also examine 44 first-generation synthetic allopolyploids of the same species. No cases of homoeologue loss arose in the first allopolyploid generation, but after 80 years, 1.6% of homoeologues were lost in natural populations. For seven homoeologue sets we also examined transcription, finding that 3.4% of retained homoeologues had been silenced in the natural populations, but none in the synthetic plants. The homoeologue losses and silencing events found were not fixed within natural populations and did not form a predictable pattern among populations. We therefore show haphazard loss and silencing of homoeologues, occurring within decades of polyploid formation in T. miscellus, but not in the initial generation.

Rapid concerted evolution of nuclear ribosomal DNA in two Tragopogon allopolyploids of recent and recurrent origin

We investigated concerted evolution of rRNA genes in multiple populations of Tragopogon mirus and T. miscellus, two allotetraploids that formed recurrently within the last 80 years following the introduction of three diploids (T. dubius, T. pratensis, and T. porrifolius) from Europe to North America. Using the earliest herbarium specimens of the allotetraploids (1949 and 1953) to represent the genomic condition near the time of polyploidization, we found that the parental rDNA repeats were inherited in roughly equal numbers. In contrast, in most present-day populations of both tetraploids, the rDNA of T. dubius origin is reduced and may occupy as little as 5% of total rDNA in some individuals. However, in two populations of T. mirus the repeats of T. dubius origin outnumber the repeats of the second diploid parent (T. porrifolius), indicating bidirectional concerted evolution within a single species. In plants of T. miscellus having a low rDNA contribution from T. dubius, the rDNA of T. dubius was nonetheless expressed. We have apparently caught homogenization of rDNA repeats (concerted evolution) in the act, although it has not proceeded to completion in any allopolyploid population yet examined.

Evolution of DNA amounts across land plants (embryophyta)

BACKGROUND AND AIMS

DNA C-values in land plants (comprising bryophytes, lycophytes, monilophytes, gymnosperms and angiosperms) vary approximately 1000-fold from approx. 0.11 to 127.4 pg. To understand the evolutionary significance of this huge variation it is essential to evaluate the phylogenetic component. Recent increases in C-value data (e.g. Plant DNA C-values database; release 2.0, January 2003; http://www.rbgkew.org.uk/cval/homepage.html) together with improved consensus of relationships between and within land plant groups makes such an analysis timely.

METHODS

Insights into the distribution of C-values in each group of land plants were gained by superimposing available C-value data (4119 angiosperms, 181 gymnosperms, 63 monilophytes, 4 lycophytes and 171 bryophytes) onto phylogenetic trees. To enable ancestral C-values to be reconstructed for clades within land plants, character-state mapping with parsimony and MacClade was also applied.

KEY RESULTS AND CONCLUSIONS

Different land plant groups are characterized by different C-value profiles, distribution of C-values and ancestral C-values. For example, the large ( approximately 1000-fold) range yet strongly skewed distribution of C-values in angiosperms contrasts with the very narrow 12-fold range in bryophytes. Further, character-state mapping showed that the ancestral genome sizes of both angiosperms and bryophytes were reconstructed as very small (i.e. < or =1.4 pg) whereas gymnosperms and most branches of monilophytes were reconstructed with intermediate C-values (i.e. >3.5, <14.0 pg). More in-depth analyses provided evidence for several independent increases and decreases in C-values; for example, decreases in Gnetaceae (Gymnosperms) and heterosperous water ferns (monilophytes); increases in Santalales and some monocots (both angiosperms), Pinaceae, Sciadopityaceae and Cephalotaxaceae (Gymnosperms) and possibly in the Psilotaceae + Ophioglossaceae clade (monilophytes). Thus, in agreement with several focused studies within angiosperm families and genera showing that C-values may both increase and decrease, it is apparent that this dynamic pattern of genome size evolution is repeated on a broad scale across land plants.

Phylogeny of Saxifragales (angiosperms, eudicots): analysis of a rapid, ancient radiation

Rapid, ancient radiations pose one of the most difficult challenges for phylogenetic estimation. We used DNA sequence data of 9,006 aligned base pairs from five genes (chloroplast atpB, matK, rbcL, and 18S and 26S nrDNA) to elucidate relationships among major lineages of Saxifragales (angiosperms, eudicots). These relationships were poorly supported in previous studies, apparently because the lineages originated in rapid succession. Using an array of methods that explicitly incorporate assumptions about evolutionary process (weighted maximum parsimony, maximum likelihood, LogDet/paralinear transformed distances), we show that the initial diversification of Saxifragales was indeed rapid. We suggest that the poor resolution of our best phylogenetic estimate is not due to violations of assumptions or to combining data partitions having conflicting histories or processes. We show that estimated branch lengths during the initial diversification are exceedingly short, and we estimate that acquiring sufficient sequence data to resolve these relationships would require an extraordinary effort (approximately 10(7) bp), assuming a linear increase in branch support with branch length. However, our simulation of much larger data sets containing a distribution of phylogenetic signal similar to that of the five sampled gene sequences suggests a limit to achievable branch support. Using statistical tests of differences in the likelihoods of topologies, we evaluated whether the initial radiation of Saxifragales involved the simultaneous origin of major lineages. Our results are consistent with predictions that resolving the branching order of rapid, ancient radiations requires sampling characters that evolved rapidly at the time of the radiation but have since experienced a slower evolutionary rate.

Phylogenetic relationships and evolution in Chrysosplenium (Saxifragaceae) based on matK sequence data

Chrysosplenium (Saxifragaceae) consists of 57 species widely distributed in temperate and arctic regions of the Northern Hemisphere, with two species restricted to the southern part of South America. Species relationships within the genus are highly problematic. The genus has traditionally been divided into two groups, sometimes recognized as sections (Oppositifolia and Alternifolia), based on leaf arrangement, or, alternatively, into 17 series. Based on morphological features, Hara suggested that the genus originated in South America and then subsequently migrated to the Northern Hemisphere. We conducted phylogenetic analyses of DNA sequences of the chloroplast gene matK for species of Chrysosplenium to elucidate relationships, test Hara's biogeographic hypothesis for the genus, and examine chromosomal and gynoecial diversification. These analyses revealed that both sections Oppositifolia and Alternifolia are monophyletic and form two large sister clades. Hence, leaf arrangement is a good indicator of relationships within this genus. Hara's series Pilosa and Macrostemon are each also monophyletic; however, series Oppositifolia, Alternifolia, and Nepalensia are clearly not monophyletic. MacClade reconstructions suggest that the genus arose in Eastern Asia, rather than in South America, with several independent migration events from Asia to the New World. In one well-defined subclade, species from eastern and western North America form a discrete clade, with Old World species as their sister group, suggesting that the eastern and western North American taxa diverged following migration to that continent. The South American species forms a clade with species from eastern Asia; this disjunction may be the result of ancient long-distance dispersal. Character mapping demonstrated that gynoecial diversification is dynamic, with reversals from inferior to half-inferior ovaries, as well as to ovaries that appear superior. Chromosomal evolution also appears to be labile with several independent origins of n = 12 (from an original number of n = 11) and multiple episodes of aneuploidy.

Phylogenetic relationships of Loasaceae subfamily Gronovioideae inferred from matK and ITS sequence data

Members of subfamily Gronovioideae are distinctive among Loasaceae in their androecial and gynoecial simplicity. The four genera of the subfamily differ, however, in chromosome number, floral novelties, and pollen exine sculpturing, which led to suggestions that the Gronovioideae were polyphyletic. Phylogenetic analyses based on sequences of the chloroplast gene matK and the internal transcribed spacer region (ITS) of nuclear rDNA have been conducted using parsimony and maximum likelihood methods to assess the monophyly of Gronovioideae and to determine the sister group relationships of gronovioid genera. The results show Gronovioideae are monophyletic and placed as the sister to Mentzelia. Within Gronovioideae, Petalonyx is sister to a clade consisting of Cevallia, Gronovia, and Fuertesia. Among the remaining Loasaceae, subfamily Mentzelioideae, as originally circumscribed, is paraphyletic. Subfamily Loasoideae is placed as the sister to the Gronovioideae-Mentzelia clade.

Structural homology and developmental transformations associated with ovary diversification in Lithophragma (Saxifragaceae)

Lithophragma, comprising only ten species, encompasses a remarkable diversity of ovary positions, reported to range from inferior to superior. The structural homology of the gynoecium and developmental transformations associated with ovary diversification are investigated for Lithophragma. Scanning electron and light microscopy indicate that all species of Lithophragma have epigynous flowers. Lithophragma campanulatum, L. glabrum, and L. heterophyllum have ovaries that externally appear nearly superior, but are actually shallowly inferior or "pseudosuperior." The inferior ovaries of Lithophragma species can be conceptually divided into superior and inferior regions that meet at the point of perianth and androecial insertion. Static and ontogenetic allometry reveal that across the species of Lithophragma the lengths of these two ovary regions are coordinated. Ovary regions in mature flowers display an approximately linear relationship that can be expressed through the allometric equation SL = -0.5314 IL + 2.0348 (where SL and IL are the lengths of the superior and inferior regions of the ovary, respectively; r = 0.7683, df = 35, P = 2.45 × 10). Mapping ontogenetic allometries onto a recent phylogeny for Lithophragma shows that ovary position evolution is bidirectional and has shifted toward greater superiority in some species and greater inferiority in others.

The role of genetic and genomic attributes in the success of polyploids

In 1950, G. Ledyard Stebbins devoted two chapters of his book Variation and Evolution in Plants (Columbia Univ. Press, New York) to polyploidy, one on occurrence and nature and one on distribution and significance. Fifty years later, many of the questions Stebbins posed have not been answered, and many new questions have arisen. In this paper, we review some of the genetic attributes of polyploids that have been suggested to account for the tremendous success of polyploid plants. Based on a limited number of studies, we conclude: (i) Polyploids, both individuals and populations, generally maintain higher levels of heterozygosity than do their diploid progenitors. (ii) Polyploids exhibit less inbreeding depression than do their diploid parents and can therefore tolerate higher levels of selfing; polyploid ferns indeed have higher levels of selfing than do their diploid parents, but polyploid angiosperms do not differ in outcrossing rates from their diploid parents. (iii) Most polyploid species are polyphyletic, having formed recurrently from genetically different diploid parents. This mode of formation incorporates genetic diversity from multiple progenitor populations into the polyploid "species"; thus, genetic diversity in polyploid species is much higher than expected by models of polyploid formation involving a single origin. (iv) Genome rearrangement may be a common attribute of polyploids, based on evidence from genome in situ hybridization (GISH), restriction fragment length polymorphism (RFLP) analysis, and chromosome mapping. (v) Several groups of plants may be ancient polyploids, with large regions of homologous DNA. These duplicated genes and genomes can undergo divergent evolution and evolve new functions. These genetic and genomic attributes of polyploids may have both biochemical and ecological benefits that contribute to the success of polyploids in nature.

Phylogenetics of flowering plants based on combined analysis of plastid atpB and rbcL gene sequences

Following (1) the large-scale molecular phylogeny of seed plants based on plastid rbcL gene sequences (published in 1993 by Chase et al., Ann. Missouri Bot. Gard. 80:528-580) and (2) the 18S nuclear phylogeny of flowering plants (published in 1997 by Soltis et al., Ann. Missouri Bot. Gard. 84:1-49), we present a phylogenetic analysis of flowering plants based on a second plastid gene, atpB, analyzed separately and in combination with rbcL sequences for 357 taxa. Despite some discrepancies, the atpB-based phylogenetic trees were highly congruent with those derived from the analysis of rbcL and 18S rDNA, and the combination of atpB and rbcL DNA sequences (comprising approximately 3000 base pairs) produced increased bootstrap support for many major sets of taxa. The angiosperms are divided into two major groups: noneudicots with inaperturate or uniaperturate pollen (monocots plus Laurales, Magnoliales, Piperales, Ceratophyllales, and Amborellaceae-Nymphaeaceae-Illiciaceae) and the eudicots with triaperturate pollen (particularly asterids and rosids). Based on rbcL alone and atpB/rbcL combined, the noneudicots (excluding Ceratophyllum) are monophyletic, whereas in the atpB trees they form a grade. Ceratophyllum is sister to the rest of angiosperms with rbcL alone and in the combined atpB/rbcL analysis, whereas with atpB alone, Amborellaceae, Nymphaeaceae, and Illiciaceae/Schisandraceae form a grade at the base of the angiosperms. The phylogenetic information at each codon position and the different types of substitutions (observed transitions and transversions in the trees vs. pairwise comparisons) were examined; taking into account their respective consistency and retention indices, we demonstrate that third-codon positions and transitions are the most useful characters in these phylogenetic reconstructions. This study further demonstrates that phylogenetic analysis of large matrices is feasible.

Timing the eastern Asian-eastern North American floristic disjunction: molecular clock corroborates paleontological estimates

Sequence data of the chloroplast gene rbcL were used to estimate the time of the well-known eastern Asian-eastern North American floristic disjunction. Sequence divergence of rbcL was examined for 22 species of 11 genera (Campsis, Caulophyllum, Cornus, Decumaria, Liriodendron, Menispermum, Mitchella, Pachysandra, Penthorum, Podophyllum, and Phryma) representing a diverse array of flowering plants occurring disjunctly in eastern Asia and eastern North America. Divergence times of putative disjunct species pairs were estimated from synonymous substitutions, using rbcL molecular clocks calibrated for Cornus. Relative rate tests were performed to assess rate constancy of rbcL evolution among lineages. Corrections of estimates of divergence times for each species pair were made based on rate differences of rbcL between Cornus and other species pairs. Results of these analyses indicate that the time of divergence of species pairs examined ranges from 12.56 +/- 4.30 million years to recent (<0.31 million years), with most within the last 10 million years (in the late Miocene and Pliocene). These results suggest that the isolation of most morphologically similar disjunct species in eastern Asia and eastern North America occurred during the global climatic cooling period that took place throughout the late Tertiary and Quaternary. This estimate is closely correlated with paleontological evidence and in agreement with the hypothesis that considers the eastern Asian-eastern North American floristic disjunction to be the result of the range restriction of a once more or less continuously distributed mixed mesophytic forest of the Northern Hemisphere that occurred during the late Tertiary and Quaternary. This implies that in most taxa the disjunction may have resulted from vicariance events. However, long-distance dispersal may explain the disjunct distribution of taxa with low divergence, such as Menispermum.

Diversification of the North American shrub genus Ceanothus (Rhamnaceae): conflicting phylogenies from nuclear ribosomal DNA and chloroplast DNA

Ceanothus comprises ∼55 morphologically and ecologically diverse species of woody perennials endemic to North America. Interpretations of the natural history of Ceanothus have served as a general model of evolution for woody perennials with simple entomophilous pollination systems, but these interpretations lacked explicit phylogenetic context. We used cladistic analysis of sequences of the chloroplast-encoded matK and the internal transcribed spacers (ITS) and 5.8S coding region of nuclear ribosomal DNA (nrDNA) to reconstruct the phylogeny of Ceanothus. The nuclear and organellar phylogenies exhibited very low levels of both topological and character congruence. Subgenera Ceanothus and Cerastes are monophyletic sister taxa in both phylogenies, but both data sets suffer from a lack of resolution below the level of subgenus. Lack of taxonomic congruence between the two data sets may be a result of introgression and/or lineage sorting. The ITS tree was accepted as the better estimate of a species phylogeny for Ceanothus, on the assumption that nuclear markers are less prone to introgression. Three of five polytypic species in the ITS data set were paraphyletic, and four of six polytypic species in the matK data set were paraphyletic. This study demonstrates the degree to which matched independent data sets can produce conflicting summaries of evolutionary history.

Phylogenetic relationships of the silver saxifrages (Saxifraga, sect. Ligulatae haworth): implications for the evolution of substrate specificity, life histories, and biogeography

The silver saxifrages (Saxifraga sect. Ligulatae Haworth; Saxifragaceae) exhibit remarkable variation of substrate specialization, with strictly calcicole to calcifuge species, as well as life histories which range from semelparity to iteroparity. They occur almost exclusively in the European mountain ranges and display high levels of endemism. Sequences from chloroplast and nuclear ribosomal DNA were obtained to resolve phylogenetic relationships among the silver saxifrages and related taxa and to gain insight into the evolution of substrate specificity, life history, and biogeography. The resulting phylogenies suggested that (1) Saxifraga sect. Ligulatae, as traditionally defined, does not constitute a monophyletic group; (2) lime-secreting hydathodes in calcifuge species apparently represent a secondary nonaptation; (3) semelparity evolved independently two or three times in the silver saxifrages and allied sections, possibly in response to climatic changes that occured during the Pleistocene; and (4) narrow endemics, for example S. cochlearis, likely evolved from the fragmentation of the widespread S. paniculata into refugial populations that became isolated during the glacial maxima of the Pleistocene.

The phylogeny of land plants inferred from 18S rDNA sequences: pushing the limits of rDNA signal?

Previous studies of the phylogeny of land plants based on analysis of 18S ribosomal DNA (rDNA) sequences have generally found weak support for the relationships recovered and at least some obviously spurious relationships, resulting in equivocal inferences of land plant phylogeny. We hypothesized that greater sampling of both characters and taxa would improve inferences of land plant phylogeny based on 18S rDNA sequences. We therefore conducted a phylogenetic analysis of complete (or nearly complete) 18S rDNA sequences for 93 species of land plants and 7 green algal relatives. Parsimony analyses with equal weighting of characters and characters state changes and parsimony analyses weighting (1) stem bases half as much as loop bases and (2) transitions half as much as transversions did not produce substantially different topologies. Although the general structure of the shortest trees is consistent with most hypotheses of land plant phylogeny, several relationships, particularly among major groups of land plants, appear spurious. Increased character and taxon sampling did not substantially improve the performance of 18S rDNA in phylogenetic analyses of land plants, nor did analyses designed to accommodate variation in evolutionary rates among sites. The rate and pattern of 18S rDNA evolution across land plants may limit the usefulness of this gene for phylogeny reconstruction at deep levels of plant phylogeny. We conclude that the mosaic structure of 18S rDNA, consisting of highly conserved and highly variable regions, may contain historical signal at two levels. Rapidly evolving regions are informative for relatively recent divergences (e.g., within angiosperms, seed plants, and ferns), but homoplasy at these sites makes it difficult to resolve relationships among these groups. At deeper levels, changes in the highly conserved regions of small-subunit rDNAs provide signal across all of life. Because constraints imposed by the secondary structure of the rRNA may affect the phylogenetic information content of 18S rDNA, we suggest that 18S rDNA sequences be combined with other data and that methods of analysis be employed to accommodate these differences in evolutionary patterns, particularly across deep divergences in the tree of life.

The earliest angiosperms: evidence from mitochondrial, plastid and nuclear genomes

Angiosperms have dominated the Earth's vegetation since the mid-Cretaceous (90 million years ago), providing much of our food, fibre, medicine and timber, yet their origin and early evolution have remained enigmatic for over a century. One part of the enigma lies in the difficulty of identifying the earliest angiosperms; the other involves the uncertainty regarding the sister group of angiosperms among extant and fossil gymnosperms. Here we report a phylogenetic analysis of DNA sequences of five mitochondrial, plastid and nuclear genes (total aligned length 8,733 base pairs), from all basal angiosperm and gymnosperm lineages (105 species, 103 genera and 63 families). Our study demonstrates that Amborella, Nymphaeales and Illiciales-Trimeniaceae-Austrobaileya represent the first stage of angiosperm evolution, with Amborella being sister to all other angiosperms. We also show that Gnetales are related to the conifers and are not sister to the angiosperms, thus refuting the Anthophyte Hypothesis. These results have far-reaching implications for our understanding of diversification, adaptation, genome evolution and development of the angiosperms.

Angiosperm phylogeny inferred from multiple genes as a tool for comparative biology

Comparative biology requires a firm phylogenetic foundation to uncover and understand patterns of diversification and evaluate hypotheses of the processes responsible for these patterns. In the angiosperms, studies of diversification in floral form, stamen organization, reproductive biology, photosynthetic pathway, nitrogen-fixing symbioses and life histories have relied on either explicit or implied phylogenetic trees. Furthermore, to understand the evolution of specific genes and gene families, evaluate the extent of conservation of plant genomes and make proper sense of the huge volume of molecular genetic data available for model organisms such as Arabidopsis, Antirrhinum, maize, rice and wheat, a phylogenetic perspective is necessary. Here we report the results of parsimony analyses of DNA sequences of the plastid genes rbcL and atpB and the nuclear 18S rDNA for 560 species of angiosperms and seven non-flowering seed plants and show a well-resolved and well-supported phylogenetic tree for the angiosperms for use in comparative biology.

Pollen tube Growth and Self-Incompatibility in Heuchera micrantha var. diversifolia (Saxifragaceae)

A crossing study and an analysis of pollen tube growth were conducted in diploid and autotetraploid Heuchera micrantha var. diversifolia to distinguish between the possible mechanisms that could promote the high outcrossing rates observed and that could also result in the absence of fruit set following self-fertilization. The crossing study indicated that no fruit set occurred after self-pollination, whereas fruit set occurred in all of the hand-pollinated outcrosses. After 4 d, the self-pollinated flowers shriveled and abscised. Pollen tube growth following hand pollination was assessed in selfed and outcrossed flowers using fluorescence microscopy. The self-pollinated flowers exhibited far fewer pollen tubes than did the outcrossed flowers. Furthermore, in self-pollinated flowers, some of the pollen tubes extended into the style; fewer than one-half of the pollen tubes reached the base of the style and still fewer reached the ovules. The variable length of pollen tube growth, the uniform timing of floral abortion after self-pollination, and the absence of variability among individuals in the level of fruit set following self-pollination are all consistent with a system that lies somewhere between classic gametophytic self-incompatibility and late-acting self-incompatibility as the mechanism that is most likely operating in H. micrantha var. diversifolia. A similar "nonstandard" system may be present in other Saxifragaceae, such as Tolmiea and Lithophragma, as well as in Ribes, the sister group of Saxifragaceae. Our data also indicate that ploidal level (diploid vs. autotetraploid) has no influence on the extent or mechanism of self-incompatibiltiy in autopolyploid H. micrantha var. diversifolia.

Phylogenetic relationships of the enigmatic angiosperm family Podostemaceae inferred from 18S rDNA and rbcL sequence data

The phylogenetic relationships of some angiosperm families have remained enigmatic despite broad phylogenetic analyses of rbcL sequences. One example is the aquatic family Podostemaceae, the relationships of which have long been controversial because of major morphological modifications associated with their aquatic habit. Podostemaceae have variously been associated with Piperaceae, Nepenthaceae, Polygonaceae, Caryophyllaceae, Scrophulariaceae, Rosaceae, Crassulaceae, and Saxifragaceae. Two recent analyses of rbcL sequences suggest a possible sister-group relationship of Podostemaceae to Crassulaceae (Saxifragales). However, the branch leading to Podostemaceae was long, and use of different outgroups resulted in alternative placements. We explored the phylogenetic relationships of Podostemaceae using 18S rDNA sequences and a combined rbcL + 18S rDNA matrix representing over 250 angiosperms. In analyses based on 18S rDNA data, Podostemaceae are not characterized by a long branch; the family consistently appears as part of a Malpighiales clade that also includes Malpighiaceae, Turneraceae, Passifloraceae, Salicaceae, Euphorbiaceae, Violaceae, Linaceae, Chrysobalanaceae, Trigoniaceae, Humiriaceae, and Ochnaceae. Phylogenetic analyses based on a combined 18S rDNA + rbcL data set (223 ingroup taxa) with basal angiosperms as the outgroup also suggest that Podostemaceae are part of a Malpighiales clade. These searches swapped to completion, and the shortest trees showed enhanced resolution and increased internal support compared to those based on 18S rDNA or rbcL alone. However, when Gnetales are used as the outgroup, Podostemaceae appear with members of the nitrogen fixing clade (e.g., Elaeagnaceae, Ulmaceae, Rhamnaceae, Cannabaceae, Moraceae, and Urticaceae). None of the relationships suggested here for Podostemaceae receives strong bootstrap support. Our analyses indicate that Podostemaceae are not closely allied with Crassulaceae or with other members of the Saxifragales clade; their closest relatives, although still uncertain, appear to lie elsewhere in the rosids.

The eastern Asian and eastern and western North American floristic disjunction: congruent phylogenetic patterns in seven diverse genera

One of the most remarkable examples of intercontinental disjunction of the North Temperate Flora involves eastern Asia and eastern and western North America. Although there has been considerable interest in this phytogeographic pattern for over 150 years (e.g., Gray, 1859; Li, 1952; Graham, 1972; Boufford and Spongberg, 1983; Wu, 1983; Tiffney, 1985a, 1985b), relationships among taxa displaying the disjunction remain obscure. Understanding phylogenetic relationships is, however, a prerequisite for historical biogeographic analyses of this distributional pattern. To understand better the relationships of taxa displaying this intercontinental disjunction, phylogenetic analyses were conducted using a variety of DNA data sets for species of four genera (Cornus, Boykinia, Tiarella, and Trautvetteria) that occur in eastern Asia, eastern North America, and western North America. An area cladogram was constructed for each of the four genera, all of which show a similar pattern of relationship: the eastern Asian species are sister to all North American species. An identical phylogenetic pattern is also found in three other taxa exhibiting this disjunction (Aralia sect. Aralia, Calycanthus, and Adiantum pedatum). The congruent phylogenetic pattern found in these seven diverse genera raises the possibility of a common origin of the eastern Asia, eastern and western North America disjunction. The data are in agreement with the long-standing hypothesis that this well-known floristic disjunction represents the fragmentation of a once continuous Mixed Mesophytic forest community and suggest that the disjunction may have involved only two major vicariance events: an initial split between Eurasia and North America, followed by the isolation of floras between eastern and western North America. However, congruence between phylogenies and geographic distributions does not necessarily indicate an identical phytogeographic history. Taxa exhibiting the same phylogenetic pattern may have originated at different geological times. Analysis of divergence times using the molecular clock indicates that species of Cornus, Boykinia, and Calycanthus may have diverged at different geological times, suggesting that the floristic disjunction involving eastern Asia and North America may not be simple; it may have involved multiple historical events at very different geological times in different genera.

The phylogenetic potential of entire 26S rDNA sequences in plants

18S ribosomal RNA genes are the most widely used nuclear sequences for phylogeny reconstruction at higher taxonomic levels in plants. However, due to a conservative rate of evolution, 18S rDNA alone sometimes provides too few phylogenetically informative characters to resolve relationships adequately. Previous studies using partial sequences have suggested the potential of 26S or large-subunit (LSU) rDNA for phylogeny retrieval at taxonomic levels comparable to those investigated with 18S rDNA. Here we explore the patterns of molecular evolution of entire 26S rDNA sequences and their impact on phylogeny retrieval. We present a protocol for PCR amplification and sequencing of entire (approximately 3.4 kb) 26S rDNA sequences as single amplicons, as well as primers that can be used for amplification and sequencing. These primers proved useful in angiosperms and Gnetales and likely have broader applicability. With these protocols and primers, entire 26S rDNA sequences were generated for a diverse array of 15 seed plants, including basal eudicots, monocots, and higher eudicots, plus two representatives of Gnetales. Comparisons of sequence dissimilarity indicate that expansion segments (or divergence domains) evolve 6.4 to 10.2 times as fast as conserved core regions of 26S rDNA sequences in plants. Additional comparisons indicate that 26S rDNA evolves 1.6 to 2.2 times as fast as and provides 3.3 times as many phylogenetically informative characters as 18S rDNA; compared to the chloroplast gene rbcL, 26S rDNA evolves at 0.44 to 1.0 times its rate and provides 2.0 times as many phylogenetically informative characters. Expansion segment sequences analyzed here evolve 1.2 to 3.0 times faster than rbcL, providing 1.5 times the number of informative characters. Plant expansion segments have a pattern of evolution distinct from that found in animals, exhibiting less cryptic sequence simplicity, a lower frequency of insertion and deletion, and greater phylogenetic potential.

Inferring complex phylogenies using parsimony: an empirical approach using three large DNA data sets for angiosperms

To explore the feasibility of parsimony analysis for large data sets, we conducted heuristic parsimony searches and bootstrap analyses on separate and combined DNA data sets for 190 angiosperms and three outgroups. Separate data sets of 18S rDNA (1,855 bp), rbcL (1,428 bp), and atpB (1,450 bp) sequences were combined into a single matrix 4,733 bp in length. Analyses of the combined data set show great improvements in computer run times compared to those of the separate data sets and of the data sets combined in pairs. Six searches of the 18S rDNA + rbcL + atpB data set were conducted; in all cases TBR branch swapping was completed, generally within a few days. In contrast, TBR branch swapping was not completed for any of the three separate data sets, or for the pairwise combined data sets. These results illustrate that it is possible to conduct a thorough search of tree space with large data sets, given sufficient signal. In this case, and probably most others, sufficient signal for a large number of taxa can only be obtained by combining data sets. The combined data sets also have higher internal support for clades than the separate data sets, and more clades receive bootstrap support of > or = 50% in the combined analysis than in analyses of the separate data sets. These data suggest that one solution to the computational and analytical dilemmas posed by large data sets is the addition of nucleotides, as well as taxa.

Chloroplast gene sequence data suggest a single origin of the predisposition for symbiotic nitrogen fixation in angiosperms

Of the approximately 380 families of angiosperms, representatives of only 10 are known to form symbiotic associations with nitrogen-fixing bacteria in root nodules. The morphologically based classification schemes proposed by taxonomists suggest that many of these 10 families of plants are only distantly related, engendering the hypothesis that the capacity to fix nitrogen evolved independently several, if not many, times. This has in turn influenced attitudes toward the likelihood of transferring genes responsible for symbiotic nitrogen fixation to crop species lacking this ability. Phylogenetic analysis of DNA sequences for the chloroplast gene rbcL indicates, however, that representatives of all 10 families with nitrogen-fixing symbioses occur together, with several families lacking this association, in a single clade. This study therefore indicates that only one lineage of closely related taxa achieved the underlying genetic architecture necessary for symbiotic nitrogen fixation in root nodules.

Phylogenetic relationships of dennstaedtioid ferns: evidence from rbcL sequences

The ferns are an ancient group of vascular plants that have yielded a staggering array of systematic problems. Among fern classifications, the number of genera in some families has ranged over 10-fold, and some genera have been treated in up to five different families. Three main groups of leptosporangiate ferns have been recognized; the adiantoid, polypodioid, and dennstaedtioid lines. To clarify relationships among genera and families of dennstaedtioid ferns, we sequenced 1320 bp of the chloroplast gene rbcL from 45 species representing 13 families. Sequence divergence for rbcL averaged 0.9% among species within genera, 10.3% among genera within families, and 14.8% among families, suggesting that the data are appropriate for phylogenetic analysis at the generic and familial levels in ferns. Maximum parsimony analysis resulted in four shortest trees of equal length. The strict consensus tree supported many aspects of previously published hypotheses of relationship based on morphological and cytological variations. For example, the tree ferns (which form a single clade) and Hymenophyllaceae appear as sister groups to the dennstaedtioid ferns on all shortest trees. However, Polypodiaceae and adiantoid ferns, groups traditionally considered separate from the dennstaedtioid families, emerged within the dennstaedtioid clade. This analysis also suggests relationships of some problematical genera, such as Monachosorum, Calochlaena, and Lonchitis. Examination of additional DNA sequences of nuclear genes and developmental studies are needed to evaluate further the relationships suggested by phylogenetic analysis of rbcL sequence data.

Genetic variation and random drift in autotetraploid populations

The rate of decay of genetic variation is determined for randomly mating autotetraploid populations of finite size, and the equilibrium homozygosity under mutation and random drift is calculated. It is shown that heterozygosity is lost at a slower rate than in diploid populations, and that the equilibrium heterozygosity with mutation and random drift is higher than for diploids. Outcrossing populations as well as populations that randomly self are analyzed. A method of comparing genetic variation between autotetraploid and diploid populations is proposed. Our treatment suggests that the "gametic homozygosity" provides a unified approach for comparing genotypes within a population as well as comparing genetic variation between populations with different levels of ploidy.

An rbcL sequence from a Miocene Taxodium (bald cypress)

During the past decade, ancient DNAs from both animals and plants have been successfully extracted and analyzed. Recently, the age of DNA that can be recovered and sequenced was increased manyfold by the amplification and sequencing of a DNA fragment from a Magnolia fossil obtained from the Miocene Clarkia deposit (17-20 million yr old). However, the validity of this report has been questioned based on models predicting that DNA should be completely degraded after 4 million yr. We report here the successful amplification, sequencing, and analysis of a 1320-base-pair portion of the chloroplast gene rbcL from a Miocene Taxodium specimen, also from the Clarkia site. These data not only validate the earlier report of sequence data for a Magnolia species from the same site but also suggest that it may be possible to isolate and sequence DNAs routinely from the Clarkia deposit. The ability to recover and sequence DNAs of such age offers enormous research possibilities in the areas of molecular evolution, biogeography, and systematics.

rbcL sequence divergence and phylogenetic relationships in Saxifragaceae sensu lato

Phylogenetic relationships are often poorly understood at higher taxonomic levels (family and above) despite intensive morphological analysis. An excellent example is Saxifragaceae sensu lato, which represents one of the major phylogenetic problems in angiosperms at higher taxonomic levels. As originally defined, the family is a heterogeneous assemblage of herbaceous and woody taxa comprising 15 subfamilies. Although more recent classifications fundamentally modified this scheme, little agreement exists regarding the circumscription, taxonomic rank, or relationships of these subfamilies. The recurrent discrepancies in taxonomic treatments of the Saxifragaceae prompted an investigation of the power of chloroplast gene sequences to resolve phylogenetic relationships within this family and between the Saxifragaceae and other major plant lineages. Sequence data from the gene rbcL (ribulose-1,5-bisphosphate carboxylase, large subunit) reveal that (i) Saxifragaceae sensu lato is at least paraphyletic, and probably polyphyletic, (ii) the genera Parnassia and Brexia are only distantly related to other members of Saxifragaceae, and (iii) representatives of the Solanaceae (subclass Asteridae) appear more closely related to Saxifragaceae (subclass Rosidase) than traditionally maintained. These data illustrate the value of chloroplast gene sequence data in resolving genetic, and hence phylogenetic, relationships among members of the most taxonomically complex groups.

Chloroplast DNA variation in a wild plant, tolmiea menziesii

Few studies of cpDNA have provided evolutionary and/or phylogenetic information at the intraspecific level. We analyzed restriction site variation using 19 endonucleases in 37 populations representing both diploid (2n = 14) and autotetraploid (2n = 28) Tolmiea menziesii. Seven restriction site mutations and five length mutations were observed. Although diploid and tetraploid Tolmiea have been intensively studied using nuclear markers, cpDNA variation provided additional evolutionary insights not revealed previously. The chloroplast genomes of diploid and tetraploid Tolmiea are as distinct as those of many pairs of congeneric species of angiosperms. Based on outgroup comparisons, the primitive chloroplast genome is present in tetraploid rather than diploid Tolmiea. These findings suggest that either: (1) diploid and tetraploid Tolmiea may have diverged since the origin of the autotetraploid, (2) the original diploid donor of the cytoplasm present in the tetraploid subsequently became extinct, or (3) the diploid was actually derived from the tetraploid via polyhaploidy. cpDNA variation also revealed that despite their close geographic proximity, diploid and tetraploid Tolmiea do not experience cytoplasmic gene flow. Last, three cytoplasmically distinct groups of diploid populations exist, two of which occupy distinct geographic areas. These findings demonstrate that, at least in some plant species, restriction fragment analysis of cpDNA can provide important evolutionary and phylogenetic information at low taxonomic levels.