*Department of Biology, Trinity University;
and
Department of Biology, Institute of Molecular Evolutionary Genetics, and Life Sciences Consortium, Pennsylvania State University
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Parasitic plants provide evolutionary genetic mutants for investigation of gene function. The cessation of photosynthetic activity that has occurred in some parasites provides us with a valuable natural experimentan opportunity to examine the genetic consequences of the shutdown of a major (perhaps the major) physiological process of most green plants (dePamphilis 1995
; dePamphilis, Young, and Wolfe 1997
). Genes not known to have an essential function outside of photosynthesis can be examined to see if they are still maintained in nonphotosynthetic plants (Wolfe, Morden, and Palmer 1992
; dePamphilis 1995
). Indeed, some of the plastid (chloroplast) genes of these plants have been lost or have become pseudogenes (dePamphilis 1995
; Nickrent et al. 1998
). If the genes have remained intact, their sequences and gene products can be examined for insight into whether they still have value to the plant.
In these nonphotosynthetic plants, genes on the plastid (chloroplast) genome may be much less constrained than are the same genes in photosynthetic plants. Nearly all of the plastid genome is composed of genes either directly involved in photosynthesis or involved in transcription and translation of plastid genes. Many plastid genes, including all those encoding photosynthetic proteins, are known to be completely expendable in at least one nonphotosynthetic plant, for they have been lost or have become pseudogenes in Epifagus virginiana (dePamphilis and Palmer 1990
; Wolfe, Morden, and Palmer 1992
). The monospecific Epifagus is holoparasitic, meaning that it has lost photosynthetic ability and derives its nutritional needs from other plants. It is a member of the family Orobanchaceae in the order Lamiales, and it has lost over half of the plastid DNA relative to Nicotiana, a photosynthetic member of the closely related order Solanales. The phylogeny of Epifagus and its relatives (dePamphilis, Young, and Wolfe 1997
; Young, Steiner, and dePamphilis 1999
) reveals several separate losses of photosynthesis in Orobanchaceae, which has recently been defined to include both photosynthetic and holoparasitic plants (Young, Steiner, and dePamphilis 1999
). However, some genes are present as intact open reading frames (ORFs) in all plants of this group, even Epifagus. Although present in the holoparasites, these genes might experience less translational demand than in photosynthetic plants, due to the diminished role of the plastid in the holoparasites (Morden et al. 1991
). Thus, purifying selection may be relaxed in these genes.
One of these persistent genes, matK, is found in every species we have examined. Despite this, its sequence is one of the least conserved of plastid genes (Olmstead and Palmer 1994
; Soltis and Soltis 1998
). matK is a popular choice for plant systematics studies, from the interspecific level (Johnson et al. 1996
) to the ordinal level (Johnson and Soltis 1995
; Hilu and Liang 1997
). The coding region of matK (
1,500 bp) and its flanking regions of a few hundred base pairs of DNA on each side together constitute a group II (lariat-forming) intron residing in the trnK gene (fig. 1).
In Epifagus, the matK gene is a full-length ORF. However, it is no longer inside an intron, because the trnK exons have been deleted, leaving matK a freestanding gene (Wolfe, Morden, and Palmer 1992
). The plastid genome of Epifagus has experienced many other deletions as well. Only 42 plastid genes remain, compared with 112 genes in a typical photosynthetic plant such as Nicotiana tabacum (Wolfe, Morden, and Palmer 1992
). In addition to the photosynthetic genes, a number of genes involved in gene expression have also been lost, including about half of the tRNAs, trnK among them. matK, however, appears to have been retained to aid in the splicing of other group II introns of the plastid genome (Ems et al. 1995
). The persistence of matK, even in Epifagus, whose plastid genome has undergone such drastic change, raises the possibility that its function may be very important, even in holoparasites. If so, we might expect that relaxation of purifying selection has not occurred. Natural selection may continue to impose constraints on the evolution of this molecule, and these constraints may be similar to those found in related photosynthetic plants.
In photosynthetic plants, and in most holoparasites as well (except for Epifagus), the DNA of the trnK intron codes for two functions. It codes for a splicing enzyme, maturase K (MatK), and because the intron is a "self-splicing" ribozyme, the regions flanking the matK gene have a folded RNA structure thought to be required for splicing (Mohr, Perlman, and Lambowitz 1993
). Although some group II introns have been observed to self-splice in vitro, they probably all require a chaperoning maturase to splice in vivo (Michel, Umesono, and Ozeki 1989
; Mohr, Perlman, and Lambowitz 1993
). MatK is thought to be distantly related to other maturases contained in other group II introns, but although many of these have reverse transcriptase (RT) domains and a domain X, the MatK protein has only the domain X and a small remnant of the RT domain. Thus, although the ancestor of MatK probably had reverse transcriptase ability, MatK has lost it. The domain X occupies a location similar to that of the "thumb" and "connection" domains of the HIV-1 RT, indicating a probable role in the binding of the intron RNA for splicing. In 34 published group II intron sequences from all kingdoms of organisms, both domain X and the remnant (regions V, VI, and VII) of the RT domain were more conserved than the rest of MatK (Mohr, Perlman, and Lambowitz 1993
).
If purifying selection is conserving the function of the MatK protein, we expect the nonsynonymous substitution rate to be lower than the synonymous substitution rate in the protein-coding region. Likewise, if substitution rates in paired-stem regions of the RNA structure are lower than substitution rates in loops, it would indicate constraint due to RNA function.
Several methods exist for comparing amounts of purifying selection. Some are based on two-taxon comparisons (Zhang, Kumar, and Nei 1997
). A more complete use of phylogenetic data is made in tree-based tests (Takezaki, Rzhetsky, and Nei 1995
; Yang 1998
; Suzuki and Gojobori 1999
). Some of these tests (Yang 1998
; Suzuki and Gojobori 1999
) are designed for protein-coding regions; the other is a distance tree-based test for sister groups (Takezaki, Rzhetsky, and Nei 1995
). Only part of the region we wish to analyze is protein-coding. We introduce a parsimony tree-based method, "subtree contrasts," which can be used for coding and noncoding regions and is not limited to sister group comparisons.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Alignment
In several previous studies (e.g., Soltis et al. 1996
; Plunkett, Soltis, and Soltis 1997
; Xiang, Soltis, and Soltis 1998
), indels within the matK coding region were easy to align by eye, but in one (Gadek, Wilson, and Quinn 1996
), indel-containing regions were removed from the phylogenetic analysis.
In this study, the data set for the trnK intron (which includes the matK gene) includes at least 96 indels, making sequence alignment a complex task. Because coding and flanking regions are expected to have different evolutionary dynamics, the matK gene was aligned separately from the flanking regions. Many alignments were generated for each region using CLUSTAL W (Thompson, Higgins, and Gibson 1994
). The one best alignment was chosen for each region using a phylogenetic optimality criterion (Wheeler 1995
). Specifically, we chose the alignment that produced parsimony trees with the least homoplasy, as judged by the rescaled consistency index (RC; Farris 1989
). For the coding-region alignment, we explored gap weights (G) of 517 and gap extension weights (E) of 17, with and without transition/transversion (ti/tv) weighting. In the flanking regions, because the lack of a reading frame was assumed to permit more frequent indel events, we explored G = 210 and E = 15, with and without ti/tv weighting. In these analyses, gaps were coded as indel characters (see Phylogenetic Analysis, below). Alignments were also refined manually and compared against the CLUSTAL alignments.
Phylogenetic Analysis
Maximum-parsimony analyses of nearly the entire intron were conducted, beginning 19 bp after trnK exon 1 and ending 26 bp before trnK exon 2 (fig. 1
). This region was approximately 2,420 bp long. Separate analyses of the coding region and of the combined flanking regions were conducted. Additional indel characters were included (23 characters for the coding region and 73 for the flanking regions). Simple gaps (those that did not overlap with other gaps) were each coded as single characters, regardless of length. Complex regions of overlapping gaps were found in both coding and flanking regions. Each of these was coded as a single multistate, unordered character (Baum, Sytsma, and Hoch 1994
). The full data set was also analyzed without these indels to determine their effect. PAUP* was used for both analyses, utilizing 1,000 random-addition replicates and tree bisection-reconnection (TBR) branch swapping. Bootstrap support (Felsenstein 1985
) was calculated for each node on the phylogeny of the entire data set.
|
|
Analysis of the Protein-Coding Region
To examine constraint on protein evolution, we used the coding region only and compared changes at third codon positions, which were mostly (70%) synonymous, with changes at first and second positions combined, which were mostly (
98%) nonsynonymous. Assuming that third codon positions were evolving in an unconstrained manner, tests of constraint were conducted by testing whether the mean number of steps at third codon positions differed from the mean at first and second positions combined. Because the number of steps per character was not normally distributed, but followed a Poisson distribution, t-tests and linear regressions could not be used. A more suitable test, Wilcoxon's nonparametric test, was used. The data take the form of a list of DNA sequence positions, each with two variables. The first tells whether the position belonged to codon position type "first and second codon positions" or "third codon positions." The second contains the number of steps. The test, in effect, asks whether there are two distributions for base changes relative to the two codon position types. The following questions were addressed: Are the bases at third codon positions changing more frequently than those at first and second positions for the photosynthetic taxa? For holoparasites? For the broomrapes alone? A "yes" to any of these questions was taken as evidence of constraint in the group examined.
We also tested whether molecular evolution in the holoparasites was less constrained than that in the photosynthetic taxa. Subtree contrasts were made using the ordinal logistic regression model (Sokal and Rohlf 1995
, pp. 767778; Sall and Lehman 1996
, pp. 224234). The data set was constructed with a list of positions. One variable was called GH and had a value of G for the data from the photosynthetic subtree and H for the data from the holoparasitic subtrees. The third variable (POS) told whether the position belonged to the class "first and second codon positions" or "third codon positions." The dependent variable was ordinal and contained the number of steps. Because so few characters had more than three steps, the "number of steps" variable was divided into four categories: 0, 1, 2, and >2. It is the significance of the interaction effect that determines whether the degrees of constraint differ, for example, among photosynthetic and holoparasitic plants of this data set. This method allows comparison of different subtrees that may have different ages, numbers of species, or amounts of molecular evolution.
The distribution of amino acid changes along the gene was generated for photosynthetic and holoparasitic subtrees. Each amino acid was coded as variant (1) or invariant (0), and a correlation was sought between the distributions based on the two subtrees.
Analysis of Flanking Regions
The same subtree contrast method was used to look for constraint on the RNA secondary structure. The flanking region was divided into bases involved in paired stem structures (stem-pair bases) and others (nonstem bases). For all taxa in this study, stem-pairing was assumed to occur as in Nicotiana (Michel, Umesono, and Ozeki 1989
). If RNA secondary structures are not conserved, the statistical tests will be weakened but may still be revealing. We also assumed (for the purposes of the test) that nonstem bases were not constrained in their evolution. If this assumption is false, and they are constrained, then the chance of detecting significant levels of constraint in the stem-pair bases is weakened but may still show whether levels of constraint differ in the two regions. Once again, the phylogeny was based on the entire intron, but the five taxa that were not sequenced for the flanking regions were removed prior to identifying subtrees. The same set of constraint tests was performed, this time on stem-pair bases versus nonstem bases.
Coding Versus Flanking Regions
To compare substitution rates in the coding region with those of the flanking regions, we chose the subtree of photosynthetic taxa for which all regions were sequenced (fig. 3D
). We also eliminated regions containing indels, because including these might underestimate rates. Again, we used the subtree contrasts method to test whether coding regions and flanking regions differed in their average substitution rates.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
When the entire Epifagus plastid sequence was published, matK was reported to be 200 bp shorter at the 5' end than in other angiosperms (Wolfe, Morden, and Palmer 1992
). However, Mohr, Perlman, and Lambowitz (1993) noted that a sequencing error resulting in a frameshift may have occurred. We have sequenced matK from an individual from the same Michigan population of Epifagus (dePamphilis and Palmer 1990
) from which the published sequence derived, as well as an individual from a population from Indiana, and have found both to contain full-length ORFs. With the exception of one nucleotide, our sequence from the Michigan individual agrees with the published sequence, so we use the corrected sequence in our analysis.
The best alignment parameters for the matK gene (under the stated optimality criterion of maximizing RC) were gap weight (G) = 15 and gap extension weight (E) = 5, with a ti/tv weight of ½. For the 700 -bp 5'-flanking region and the
230-bp 3'-flanking region, which were aligned in a combined analysis, the best parameters were G = 3 and E = 2, with a ti/tv weight of 1/1. Additional efforts to improve the alignment by eye resulted in trees with equal or lower RCs. This alignment shows that length differences are many times more common in the flanking regions than in the coding region (table 1
).
|
Protein-Coding Region
The phylogeny (fig. 2
) was then used to test for significant differences in substitution rates among codon positions as an indication of the level of constraint on the matK gene. Wilcoxon tests of evolutionary changes on photosynthetic and holoparasitic subtrees (fig. 3
) were used to examine constraint on the gene, first among photosynthetic plants and then among holoparasites. The photosynthetic subtree had a greater rate of change at third codon positions than at the other positions (2 = 50.1, df = 1, P < 0.0001; table 2
), suggesting that the gene is evolving under constraint. Third-position change averaged 1.7 (±0.4) times the (per-base-pair) change at first and second positions. We then analyzed the holoparasites. Because the number of steps on the Harveya-Hyobanche subtree was too small to analyze separately, the steps on the two holoparasitic subtrees were added together and used to represent holoparasitic evolution. In the holoparasitic subtrees, constraint on matK was also observed (
2 = 15.2, df = 1, P < 0.0001). In the broomrapes alone, constraint was evident as well (
2 = 13.2, df = 1, P = 0.0003). Using the interaction effect tests of the ordinal logistic regression model, matK in the holoparasites was significantly less constrained than in photosynthetic plants (
2 = 4.96, df = 1, P = 0.0259). Likewise, matK in the broomrapes alone was less constrained than in photosynthetic plants (
2 = 5.85, df = 1, P = 0.0156).
|
Although the protein domains tested did not show differences in constraint, when amino acid sites were examined individually, there was an excellent correlation (r2 = 0.205, P < 0.0001) between sites that were invariant on the photosynthetic subtree (fig. 3A ) and those that were invariant over the holoparasitic subtrees (fig. 3B and C ).
Intron RNA Secondary Structure
We then looked for evidence of constraint on the ribozyme structure, comparing the stem-pair regions with nonstem regions. Paired stem regions show less change than nonstem regions in photosynthetic plants (2 = 54.9, df = 1, P < 0.0001; table 2
), holoparasites combined (
2 = 53.4, df = 1, P < 0.0001), and the broomrapes alone (
2 = 51.0, df = 1, P < 0.0001). The combined holoparasites do not differ significantly in degree of constraint when compared with photosynthetic plants (
2 = 0.121, df = 1, P = 0.738). Likewise, the subtree of broomrapes alone does not differ from that of photosynthetic plants in degree of constraint (
2 = 0.228, df = 1, P = 0.633). For photosynthetic plants, the nonstem regions show 2.3 times the rate of change of the stem-pair regions.
Coding Versus Flanking Regions
We also estimated the relative substitution rates for the various regions using the subtree of photosynthetic plants for which both regions were sequenced (fig. 3D
). For these comparisons, we excluded indel-containing regions. In terms of steps per character over the photosynthetic subtree, which is a relative indication of evolutionary rate, the flanking regions averaged 0.365 and the coding region averaged 0.486. These two values were significantly different (P = 0.0005). Within the flanking regions, the nonstem regions are expected to be relatively unconstrained, as are the third codon positions of the coding region. Nonstem regions averaged 0.474 steps per character, and third codon positions averaged 0.673. These two values were significantly different (P = 0.0003).
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Parsimony-Based Constraint Analyses
Although we refer to taxa as "photosynthetic" and "holoparasitic" for the sake of convenience, we compared rates of change only in Orobanchaceae and their close relatives. All hypotheses were tested based on a phylogeny generated from the intron sequences themselves. The accuracy of this plastid DNA phylogeny does not depend on these sequences alone, but is confirmed by analyses of the rps2 gene (Young, Steiner, and dePamphilis 1999
). Trees derived from rps2 are less resolved than those based on matK, but they are congruent.
Protein Constraint
These analyses allowed us to see that both the matK gene and its flanking regions continue to evolve constrained by purifying selection, even after loss of photosynthesis. Our tests show that there is significant constraint on MatK evolution in the photosynthetic plants we examined, but the level of constraint is low, as has been reported for most angiosperms examined (table 3 ). The ratio of third-position changes to the average of the other positions is only 1.7 (±0.4). This is similar to the ratio of 1.4 for matK in the Saxifragaceae (Johnson and Soltis 1995
) and is a much lower level of constraint than is found in rbcL, which has a ratio of 5.5 in the Saxifragaceae (Johnson and Soltis 1995
). The level of constraint is not significantly reduced in domain X or domains V, VI, VII, and X. However, similar sites are invariant in the photosynthetic plants and in the holoparasites, suggesting a similarity in the constraint they have experienced.
|
Ribozyme Constraint
Both the photosynthetic plants and the holoparasites in this study show constraint on the evolution of the putative paired stem regions. The nonstem regions change 2.3 times as fast as the paired stem regions. The paired stem regions are constrained in the holoparasites at nearly the same level as in their photosynthetic relatives. This suggests that despite the loss of photosynthesis and the deletions that have occurred in the plastid of the holoparasites, the self-splicing ribozyme function of the trnK intron is still important, perhaps because the trnK gene is still necessary. Both matK and the rest of the trnK intron remain selectively constrained in the holoparasites that have them. This fact deepens the mystery behind the loss of trnK in Epifagus (Wolfe, Morden, and Palmer 1992
).
The comparison of coding and flanking regions shows an intriguing result. Among photosynthetic plants, the substitution rate in the flanking regions is only 54% of that at third codon positions in the coding region. Even when the stem-pair regions are excluded, the remaining nonstem regions have only 70% of the substitution rate of third positions, suggesting that parts of the nonstem regions are also under constraint of some kind. However, indels are much more frequent in the flanking regions (at least 61 per 1,000 bp, as compared with at least 14 per 1,000 bp for the coding region). Selection is constraining substitutions in the flanking region, while still allowing a large number of indels.
The subtree contrasts method presented here can be used to test for differences and two-way interactions whenever a phylogeny can be divided into subtrees. The subtrees can be monophyletic or paraphyletic. For example, warm-blooded and cold-blooded vertebrate subtrees could be tested for differing substitution rates and for differing degrees of constraint on protein evolution. The method could be adapted for the use of synonymous and nonsynonymous sites, rather than the first, second, and third positions used here. It also could be adapted for the use of likelihood-based estimates of branch lengths. When only one subtree in each category is analyzed, the conclusions hold only for those subtrees examined, but if multiple subtrees are available for each category, the capacity for generalization will increase (dePamphilis 1995
).
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
1 Keywords: group II intron
purifying selection
ribozyme
matK
Epifagus virginiana,
parasitic plant
subtree contrasts
2 Address for correspondence and reprints: Nelson D. Young, Department of Biology, Trinity University, San Antonio, Texas 78212. E-mail: nyoung{at}trinity.edu
![]() |
literature cited |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Baum, D. A., K. J. Sytsma, and P. C. Hoch. 1994. A phylogenetic analysis of Epilobium (Onagraceae) based on nuclear ribosomal DNA sequences. Syst. Bot. 19:363388[ISI]
Cronquist, A. 1981. An integrated system of classification of flowering plants. Columbia University Press, New York
dePamphilis, C. W. 1995. Genes and genomes. Pp. 177205 in M. C. Press and J. D. Graves, eds. Parasitic plants. Chapman and Hall, London
dePamphilis, C. W., and J. D. Palmer. 1990. Loss of photosynthetic and chlororespiratory genes from the plastid genome of a parasitic flowering plant. Nature 348:337339
dePamphilis, C. W., N. D. Young, and A. D. Wolfe. 1997. Evolution of plastid gene rps2 in a lineage of hemiparasitic and holoparasitic plants: many losses of photosynthesis and complex patterns of rate variation. Proc. Natl. Acad. Sci. USA 94:73677372
Ems, S., C. W. Morden, C. Dixon, K. H. Wolfe, C. W. dePamphilis, and J. D. Palmer. 1995. Transcription, splicing and editing of plastid RNAs in the non-photosynthetic plant Epifagus virginiana. Plant Mol. Biol. 29:721733
Farris, J. S. 1989. The retention index and the rescaled consistency index. Cladistics 5:417419
Felsenstein, J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783791
Gadek, P. A., P. G. Wilson, and C. J. Quinn. 1996. Phylogenetic reconstruction in Myrtaceae using matK, with particular reference to the position of Psiloxylon and Heteropyxis. Aust. Syst. Bot. 9:283290
Hilu, K. W., and H. Liang. 1997. The matK gene: sequence variation and application in plant systematics. Am. J. Bot. 84:830839[Abstract]
Jarrell, D. C., and M. T. Clegg. 1995. Systematic implications of the chloroplast-encoded matK gene in the tribe Vandeae (Orchidaceae) (abstract). Am. J. Bot. 82(Suppl.):137
Johnson, L. A., J. L. Schultz, D. E. Soltis, and P. S. Soltis. 1996. Monophyly and generic relationships of Polemoniaceae based on matK sequences. Am. J. Bot. 83:12071224[ISI]
Johnson, L. A., and D. E. Soltis. 1995. Phylogenetic inference in Saxifragaceae sensu stricto and Gilia (Polemoniaceae) using matK sequences. Ann. Mo. Bot. Gard. 82:149175
Li, W.-H., and D. Graur. 1991. Fundamentals of molecular evolution. Sinauer, Sunderland, Mass
Liang, H., and K. W. Hilu. 1996. Application of the matK gene sequences to grass systematics. Can. J. Bot. 74:125134[ISI]
Michel, F., K. Umesono, and H. Ozeki. 1989. Comparative and functional anatomy of group II catalytic intronsa review. Gene 82:530
Mohr, G., P. S. Perlman, and A. M. Lambowitz. 1993. Evolutionary relationships among group II intron-encoded proteins and identification of a conserved domain that may be related to maturase function. Nucleic Acids Res. 21:49914997[Abstract]
Morden, C. W., K. H. Wolfe, C. W. dePamphilis, and J. D. Palmer. 1991. Plastid translation and transcription genes in a non-photosynthetic plant: intact, missing and pseudo genes. EMBO J. 10:32813288[Abstract]
Nickrent, D. L., C. W. dePamphilis, A. D. Wolfe, A. E. Colwell, N. D. Young, and R. J. Duff. 1998. Molecular phylogenetic and evolutionary studies of parasitic plants. Pp. 211241 in D. Soltis, P. Soltis, and J. Doyle, eds. Molecular systematics of plants II: DNA sequencing. Kluwer Academic Publishers, Boston
Olmstead, R. G., and J. D. Palmer. 1994. Chloroplast DNA systematics: a review of methods and data analysis. Am. J. Bot. 81:12051224[ISI]
Plunkett, G. M., D. E. Soltis, and P. S. Soltis. 1997. Clarification of the relationship between Apiaceae and Araliaceae based on matK and rbcL sequence data. Am. J. Bot. 84:565580[Abstract]
Sall, J., and A. Lehman. 1996. JMP start statistics. Duxbury Press, Belmont, Calif
SAS Institute. 1994. JMP. SAS Institute. Cary, N.C
Sokal, R. R., and F. J. Rohlf. 1995. Biometry. 3rd edition. W. H. Freeman, New York
Soltis, D. E., R. K. Kuzoff, E. Conti, R. Gornall, and K. Ferguson. 1996. matK and rbcL gene sequence data indicate that saxifraga (Saxifragaceae) is polyphyletic. Am. J. Bot. 83:371382[ISI]
Soltis, D. E., and P. S. Soltis. 1998. Choosing an approach and an appropriate gene for phylogenetic analysis. Pp. 142 in D. Soltis, P. Soltis, and J. Doyle, eds. Molecular systematics of plants II: DNA sequencing. Kluwer Academic Publishers, Boston
Steele, K. P., and R. Vilgalys. 1994. Phylogenetic analyses of Polemoniaceae using nucleotide sequences of the plastid gene matK. Syst. Bot. 19:126142
Suzuki, Y., and T. Gojobori. 1999. A method for detecting positive selection at single amino acid sites. Mol. Biol. Evol. 16:13151328[Abstract]
Takezaki, N., A. Rzhetsky, and M. Nei. 1995. Phylogenetic tests of molecular clock and linearized trees. Proc. Natl. Acad. Sci. USA 12:823833
Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:46734680[Abstract]
Wheeler, W. C. 1995. Sequence alignment, parameter sensitivity, and the phylogenetic analysis of molecular data. Syst. Biol. 44:321331[ISI]
Wolfe, K. H., C. W. Morden, and J. D. Palmer. 1992. Function and evolution of a minimal plastid genome from a nonphotosynthetic parasitic plant. Proc. Natl. Acad. Sci. USA 89:1064810652
Xiang, Q.-Y., D. E. Soltis, and P. S. Soltis. 1998. Phylogenetic relationships of Cornaceae and close relatives inferred from matK and rbcL sequences. Am. J. Bot. 85:285297[Abstract]
Yang, Z. 1998. Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol. Biol. Evol. 15:568573[Abstract]
Young, N. D., K. E. Steiner, and C. W. dePamphilis. 1999. The evolution of parasitism in Scrophulariaceae/Orobanchaceae: plastid gene sequences refute an evolutionary transition series. Ann. Mo. Bot. Gard. 86:876893
Zhang, J., S. Kumar, and M. Nei. 1997. Small-sample tests of episodic adaptive evolution: a case study of primate lysozymes. Mol. Biol. Evol. 14:13351338