* Department of Anatomy and Cell Biology
Department of Biological Sciences, Wayne State University
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Key Words: fibrinogen intron intron alignment bird Aves conserved intron sequence avian phylogeny
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Although valuable for resolving phylogenies of recently evolved groups, it is unclear how far back in evolutionary time intron sequences might give phylogenetic resolution. Computer simulations (Moore, Smith, and Prychtiko 1999) suggest that introns do not begin to saturate significantly in avian species until they have diverged for 50 Myr, and it is likely that they retain phylogenetic information of greater antiquity. Diverging introns also accumulate indels of various lengths, which have the potential of being useful characters in establishing evolutionary relationships, but severe alignment problems may occur when more anciently diverged sequences are aligned. Prychitko and Moore (1997, 2000) demonstrated the ability of ß-fibint 7 to provide a well-resolved gene tree for woodpeckers, a subfamily of birds that appears to be approximately 7 to 8 Myr old (Moore, Smith, and Prychitko 1999). Although indels occurred in the ß-fibint 7 sequences of the woodpeckers, they did not pose serious difficulties with alignment.
The phylogenetic relationships among avian orders and among families within most avian orders are poorly understood because of the apparent rapid radiation of modern birds (Sibley and Ahlquist 1990; Hedges 1994; Feduccia 1996; Cooper and Penny 1997; Groth and Barrowclough 1999; Mindell et al. 1999). This explosive period of evolution, which occurred some time between 55 and more than 90 MYA (reviewed by Waddell et al. 1999; Mindell et al. 1999; van Tuinen and Hedges 2001), appears to have generated all of the modern orders of birds in approximately 5 to 10 Myr (Feduccia 1995). Because modern avian orders diverged within a short span of time relatively long ago, the number of synapomorphies that delineate basal relationships is limited. Morphological studies such as those by Cracraft (1988) and Cracraft and Mindell (1989) had difficulty resolving relationships among particular avian orders because of a lack of sufficient synapomorphies. Sibley and Ahlquist's (1990) extensive DNA-DNA hybridization analysis has provided a good beginning for DNA-based inference of avian phylogeny but has been criticized from a number of perspectives (Lanyon 1992; Mindell 1992; Harshman 1994). At the very least, the phylogeny hypothesized by Sibley and Ahlquist (1990) needs to be corroborated.
Several recent studies have explored the value of DNA sequence data from various genes for resolving relationships among major avian lineages (Caspers et al. 1997; Harlid, Janke, and Arnason 1997, 1998; Mindell et al. 1997, 1999; Groth and Barrowclough 1999; Garcia-Moreno and Mindell 2000; van Tuinen, Sibley, and Hedges 2000). All of these should be considered pilot studies because they were based on limited taxic samples and/or genes for which it was unclear whether they would retain phylogenetic signal sufficient to resolve relationships of this antiquity. All of these studies have produced some resolution but have produced some unsatisfactory results as well. Phylogenies based on mtDNA and a small taxic sample root the avian tree in a way that is at odds with traditional concepts of avian classification (Mindell et al. 1997, 1999; van Tuinen, Sibley, and Hedges 2000; Johnson 2001). Increasing the size of the taxic sample appeared to "solve" the rooting problem in one study (van Tuinen, Sibley, and Hedges 2000) but not in another (Johnson 2001), and Braun and Kimball (2002) recovered the traditional root in reanalyses of the mtDNA data if the analysis was restricted to transversion parsimony or maximum likelihood based on complex models. Moreover, the studies of Mindell et al. (1999) and Johnson (2001) inferred other relationships at odds with more traditional classifications, which were strongly supported by morphological characters (e.g., the suboscine and oscine passerines not being sister taxa). The study of van Tuinen et al. (2000), based on the 12S and 16S rRNA and tRNA-Valine genes from the mt-genome and the nuclear 18S rRNA gene, found strong statistical support for five basal nodes (Aves, Paleognathae, Neognathae, Galloanserae, and Neoaves) but insignificant support for any lower level relationships. A phylogenetic analysis of exon sequence from the nuclear encoded RAG-1 gene provides strong statistical support for the traditional root, splitting the ancestor of modern birds into paleognath and neognath lineages, but provides low resolution among some major lineages (Groth and Barrowclough 1999).
Although, all DNA sequencebased studies to date have been based on very limited numbers of taxa, the lack of stable resolution may result from the rates at which the various genes used in the studies are evolving as much as, or more than, the limited taxic sample. Specifically, mtDNA may be evolving too fast and nuclear gene exons evolving too slow to be optimal for resolving relationships among avian orders. A reasonable hypothesis is that nuclear gene introns, which evolve more rapidly than exons but more slowly than mitochondrial-encoded genes may be evolving at rates suitable for resolving relationships among taxa of the antiquity of avian orders and families.
The high substitution rate in introns and the frequent occurrence of indels suggest that their nucleotide sequence is adaptively neutral, although their existence must be perpetuated by natural selection because the intron/exon structure of genes is identical among distantly related taxa (Lewin 1997). ß-fibint 7, for example, has been amplified by PCR from birds (Prychitko and Moore 1997, 2000; Johnson and Clayton 2000), reptiles (Giannasi, Malhotra, and Thorpe 2001), and mammals (C. Krajewski, personal communication) utilizing the same exon-based primers. Although the extent of functional sequences within introns is not known, it is known that introns contain conserved sequences for correct splicing (Keller and Noon 1985; Mount et al. 1992). Conserved sequences that function as alternative splice sites (Leicht et al. 1993), as regulatory elements (Gasch, Hinz, and Renkowitz-Pohl 1989; Schultz et al. 1991), and in maintaining secondary structure of transcripts (Kirby, Muse, and Stephan 1995) have also been identified in some introns.
In this study, ß-fibint 7 was sequenced for 28 species of birds, representing 18 families in nine avian orders, for the following purposes: (1) to assess problems of sequence alignment among these relatively old taxa, (2) to determine if an intron is able to resolve relationships at these taxonomic levels, (3) to examine the conservation of nucleotide sequence and composition (GC content), and (4) to examine patterns of insertions and deletions (indels) over the course of evolutionary history represented by these species. Our taxon sample is very limited, but we include species that probably represent the full temporal extent of neornithine evolution (Mindell et al. 1999; van Tuinen, Sibley, and Hedges 2000). ß-fibint 7 sequences have lower levels of homoplasy and experience fewer multiple substitutions than mitochondrial-encoded protein genes (Prychitko and Moore 1997, 2000); therefore, they may retain sufficient phylogenetic signal to resolve relationships among these most ancient of extant avian lineages.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Methods of DNA isolation, PCR amplification, PCR-product cleaning, and DNA quantitation used in this study are those described by Prychitko and Moore (1997). We had access to different automated DNA sequencing instruments during the course of the study: ABI 377, ABI 3700 (Applied Biosystems), and Pharmacia ALF. For the ABI sequencers, reactions were done using the Taq Big Dye primer cycle sequencing kit (Applied Biosystems) according to the manufacturer's specifications. For the ALF sequencer, reactions were done using the Thermo Sequenase fluorescent-labeled primer cycle sequencing kit (Amersham) according to the protocols described by Prychitko and Moore (1997) and Prychitko, Ries, and Moore (1998).
Sequence Construction and Alignment
The entire ß-fibint 7 sequence for each specimen was reconstructed by linking together overlapping fragments generated using primers FIB-BI7U, FIB-BI7L, FIB-B17U2, FIBL210, FIBU330, and FIB-BI7L2 (Prychtiko and Moore 1997). These primer combinations result in approximately 70% overlap in sequence from the sense and nonsense strands. Two or more DNA isolations were made for each specimen, and in several cases, an isolation was made from a second specimen representing the same species. ß-fibint 7 was amplified and sequenced for these additional isolations and used to verify the sequences reported in this study.
Sequences were aligned with ClustalX (Thompson et al. 1997; Hall 2001). To explore the range of plausible alignments and the consistency of alignment procedures, we initiated several alignments (see Results and Discussion), following a different strategy in each. The alignment presented in the Supplementary Materials and used in further analyses was obtained as follows: Members of the order Piciformes were aligned first, beginning with the infraorder Picides (northern flicker, bar-breasted piculet, Eurasian wryneck, and scaly-throated honeyguide), because we are relatively confident of their relationships and had previously aligned these sequences (Prychitko and Moore 2000). The black-spotted barbet and ivory-billed aracari (infraorder Ramphastides) sequences were aligned separately and then added to the overall Piciformes alignment as a profile (Thompson et al. 1997), because both sequences contain sets of three large, identical deletions totaling 303 nucleotides, and New World barbets and toucans (aracari) are thought to be sister taxa (Lanyon and Hall 1994). Individual sequences and profiles representing other avian orders were subsequently added to the alignment in random order. When an avian order was represented by more than one species, sequences were aligned within the order, before adding that order to the overall alignment. When a new avian order was added to the alignment, a "profile" alignment was performed whereby the sequences to be added created a profile that was aligned to the profile comprising the overall alignment.
Phylogenetic Analysis
Phylogenetic analyses were performed with PAUP* version 4.0b5 (Swofford 1998). Details of each analysis are given in the Results section.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Based on our alignment, we judge that there is a total of 125 indels in the data set.
Base Composition
ß-fibint 7 is AT-rich in woodpeckers (Prychitko and Moore 1997, 2000) and doves (Johnson and Clayton 2000). In this analysis involving avian species representing nine orders, the ß-fibint 7 sequences again consisted of a high percentage of A and T (see table 1). For the 28 species sequenced in this study, the ß-fibint 7 sequences on average consisted of almost 65% A and T, and there is little variance in the base composition among orders (see table 1).
|
|
|
Phylogenetic Analysis
We performed a total of four phylogenetic analyses, three based on nucleotide substitutions only (maximum parsimony [MP], Neighbor-Joining [NJ] and maximum likelihood [ML]) and one based on nucleotide substitutions plus indels (MP).
All MP searches were TBR (tree bisection-reconnection) heuristic searches with 10 replicates of randomly added OTUs conducted with characters unordered and equally weighted. The nucleotide-only MP search found two equally parsimonious topologies (1,155 characters, length = 1,669 steps, CI [consistency index] excluding uninformative characters = 0.61); the two topologies differ with regard to where the ovenbird joins the tree. One of the two MP topologies is shown in figure 3.
|
We used the Tajima-Nei distance estimator to compute the Neighbor-Joining tree because of the apparent ts/tv ratio close to 1 and a base composition that differs substantially from 25% representation by each nucleotide. Tajima-Nei is therefore appropriate because it is the simplest model that includes the essential nucleotide substitution properties (Kumar, Tamura, and Nei 1993). The NJ tree topology computed is identical to that of the MP tree illustrated in figure 3.
The computational demands of ML estimation precluded an exhaustive search based on this many taxa and characters. As an alternative, we did a partial search based on a stepwise strategy simplified from Frati et al. (1997). We estimated the model parameters for the GTR (general time reversible) model, gamma corrected (four rate categories) with a proportion of invariant sites by ML based on the topology of the NJ tree. This gives good approximations because the parameter estimates are not strongly dependent on the topology (Yang, Goldman, and Friday 1995). The estimated parameter values were then fixed, and the ML topology and branch lengths were estimated. We then specified this topology as a user tree to refine the ML estimates of the model parameters and branch lengths. The topology illustrated in figure 3 proved to be the ML tree, as estimated by this procedure, as well as the NJ tree and one of the MP trees.
One thousand bootstrap replicates were generated for the MP analysis with indels; these are presented in figure 3.
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Although programs such as ClustalX facilitate alignment of multiple sequences, there is an element of subjectivity in judging which alignment best reflects the evolutionary history of indels. The several alignment procedures we undertook differed in the order in which taxa were entered into the alignment, in gap-opening penalties (GOP) and gap-extension penalties (GEP), and in the extent to which profiles and the final alignment were fine tuned by ClustalX tools or manually. Nevertheless, all of the resultant alignments were very similar and produced only minor differences in subsequent phylogenetic analysis. With regard to aligning intron sequences such as ß-fibint 7, with large as well as small indels, it was our experience that specifying high GOP (40+) and low GEP (0.05) were important to identify large indels but that subsequent fine tuning of small, restricted regions was facilitated by reducing the GOP (<10).
Alignment of multiple DNA sequences for subsequent phylogenetic analysis is fraught with a potential circularity: Better alignment is achieved if sequences are entered in order of phylogenetic relatedness, but knowledge of relatedness is based on the alignment. To facilitate alignment, we first aligned sequences from species traditionally assigned to the same order (Piciformes, Strigiformes, etc. [see fig. 3]). However, we assumed no knowledge of relationships beyond this point, as sequences and profiles representing distinct orders were entered into the growing alignment in random order. The assignment of species in our study to orders is not contentious; therefore the tree we inferred (fig. 3) should not be an artifact of circular reasoning.
In sum, alignment of this intron among species representing distinct avian families and orders is time consuming, tedious, and uncertain in restricted regions for some species, but it is not an insurmountable or even difficult task. However, adding reptilian or mammalian sequences to the analysis manifested problems both with alignment and with detection of phylogenetic signal. We attempted to amplify and sequence ß-fibint 7 from a crocodilian species (Alligator mississippiensis), intending to use it as an outgroup, but our amplifications were unsuccessful. Although we believe ß-fibint 7 from the alligator could be amplified and sequenced, we abandoned the effort when we obtained snake (Giannasi, Malhotra, and Thorpe 2001) and mammal sequences (C. Krajewski, personal communication). We were unable to align the snake and mammal sequences with our avian sequences using ClustalX, but we were able to align them using the computer program T-Coffee (C. Notredame, http://www.ch.embnet.org/software/TCoffee.html). However, the reptile and mammal sequences were so disparate relative to the avian sequences that they did not provide a strong, statistically supported root for the phylogeny. This appears to establish a limit for the use of ß-fibint 7 in antiquity. The lineages leading to mammals and snakes are thought to have diverged from that leading to birds at least 310 and 260 MYA, respectively (Benton 1990). On the other hand, ß-fibint 7 is readily aligned among the modern avian orders represented in this study, which are thought to have diverged from 55 to more than 90 MYA, and it retains significant information on relationships among the orders and families, as well as the lower taxa. Thus, it appears that ß-fibint 7 would certainly be useful in phylogenetic studies of groups as old as 55 Myr, probably as old as 90 Myr, but not as old as 260 Myr. At the other extreme, ß-fibint 7 has the potential to resolve relations among species of birds that diverged as recently as 2 to 5 MYA (Prychitko and Moore 2000) and perhaps less, although mitochondrial-encoded genes appear more useful for these very recent divergences (Moore, Smith, and Prychitko 1999; Weibel and Moore 2002). The estimate of this window of resolution is only approximate because resolution of specific relationships depends on the length and depth of the internode (Lanyon 1988; Moore and DeFilippis 1997).
Nucleotide Base Composition
It was shown previously that base composition of ß-fibint 7 varied little among woodpeckers (Prychitko and Moore 2000) or among doves (Johnson and Clayton 2000). The present analysis includes representatives from a broad range of avian orders, and base composition remains remarkably constant across lineages. Constancy of base composition is important to phylogenetic inference because variation can result in clustering of unrelated taxa with similar base composition into artificial groups (Lockhart et al. 1994). The genome structure of warm-blooded animals is characterized by isochores, homogeneous regions of DNA, 100 to 300 kb or longer, that are either AT rich or GC rich (Bernardi 1995). Whether the homogenous base composition of isochores is maintained by selection or a biased mutational mechanism is unclear (Li 1997), but presumably ß-fibrinogen is embedded in an AT-rich isochore; the uniform AT richness of ß-fibint 7 over the diversity of modern birds suggests such a condition.
Conserved Sequences Within Introns
It is clear from figure 2 that comparison of aligned sequences from distantly related species is a powerful tool for identifying conserved, and presumably functional, regions within genomes. In our earlier analyses (Prychitko and Moore 1997, 2000), we failed to detect these regions because the woodpeckers are closely related, and too few substitutions have occurred within the family (Picidae) to reveal variation in substitution rate along the length of the intron.
Since the discovery of introns in the 1970s, their nucleotide sequences were thought not to specify function and hence not to be constrained by selection (reviewed by Rodriguez-Trelles, Tarrio, and Ayala 2000; Mattick and Gagen 2001). However, there are now reports of conserved regions within introns where it was inferred that these regions function as alternative splice sites (Leicht et al. 1993, 1995; Croft et al. 2000), as regulatory elements (Gasch, Hinz, and Renkowitz-Pohl 1989; Schultz et al. 1991; Alder et al. 1992; Jackson and Hoffmann 1994; Metz et al. 1998), and as stem regions maintaining secondary structure of transcripts (Kirby, Muse, and Stephan 1995). (Several additional examples of introns affecting gene expression are listed on the ISIS Web site <http://isis.bit.uq.au/papers_function.html>.) Spliceosomal introns characteristically begin with the dinucleotide GT and end with AG (Keller and Noon 1985; Mount et al. 1992). As expected, these four sites are conserved across all avian orders (fig. 2). The three additional conserved sites at the 5' end of the intron (sites 5, 6, and 12) and two at the 3' end (1132 and 1147 [fig. 2]) plausibly function in spicing also. Introns characteristically have a short, conserved region (usually less than 10 nucleotides in length) towards the 3' end, the branch point, which includes the bonding site for the G at the 5' end as it loops back on the intron to form the "lariat" during excision (Mount et al. 1992; Lewin 1997). Typically, there is also a conserved pyrimidinerich region located between the 3' splice site and the branch point, which in mammalian species has been shown to play a role in branch point recognition during splicing (Ruskin et al. 1984; Csank, Taylor, and Martindale 1990; Mount et al. 1992). Although the conserved sites near the 3' end of ß-fibint 7 (fig. 2) do not precisely satisfy these criteria, they are plausibly associated with the branch point.
At least two conspicuous clusters of conserved nucleotides remain beyond the minimal set of functionally conserved nucleotides usually associated with excision: the sequences between sites 139 and 176 and between sites 571 and 683 (fig. 2). Although these sequences are almost certainly constrained by selection, their function is unknown.
Mattick and Gagen (2001) argued that intron transcripts may play an important, but largely overlooked, role in modulating the expression of genes during development. The existence of constrained sequences within ß-fibint 7, which seem to be in excess of those required for excision, is consistent with this hypothesis to the extent of establishing that there are conserved regions and at the same time could explain why intron-embedded regulatory elements have not been detected often in the past. The proportion of conserved sites in a large intron is small (46 of 1,155 sites in ß-fibint 7), and hence introns as a whole evolve at rates approaching those of pseudogenes and synonymous sites of exons (Hughes and Yeager 1997; Li 1997). Moreover, these short, conserved sequences would only be detected if introns were aligned among appropriately distantly related species. This has been done rarely to date.
Avian Phylogeny
Because the mammal and reptile ß-fibint 7 sequences were too distant from the avian sequences to provide a satisfactory outgroup root, we rooted our tree by designating the ostrich as the outgroup. The basal lineage of the avian phylogeny remains somewhat controversial (Mindell et al. 1997, 1999; Garcia-Moreno and Mindell 2000), but we think the preponderance of evidence, especially recent studies based on more slowly evolving nuclear gene sequences, favors the classical basal division of Aves into Paleognathae and Neognathae (Pycraft 1900; Cracraft 1988; Cracraft and Mindell 1989; reviewed by van Tuinen, Sibley, and Hedges 2000). Recent nuclear DNA sequence studies favoring the classical division include Groth and Barrowclough (1999), van Tuinen, Sibley, and Hedges (2000), and Garcia-Moreno and Mindell (2000) (see Braun and Kimball 2002 for a recent review). Some mtDNA studies based on additional sequences and more complex analyses also favor the classical division (Braun and Kimball 2002; Paton, Haddrath, and Baker 2002), as does a novel approach based on nuclear DNA "strings" (Edwards et al. 2002). This is our justification for outgroup rooting with the ostrich designated as the outgroup, but we acknowledge an element of uncertainty regarding the basal divergences depicted in the ß-fibint 7 phylogeny (fig. 3).
Although our study is intended as a pilot study to determine whether ß-fibint 7 has retained phylogenetic signal from the depth in antiquity of avian orders, the limited taxon sample does provide some useful insight into early avian evolution and classification. First, all of the orders represented by two or more species appear monophyletic in the ML tree (Piciformes, Passeriformes, Strigiformes, Columbiformes, Gruiformes, and Anseriformes) (fig. 3). Moreover, the relationships among the avian orders are strongly supported by bootstrap proportions with the exception of the Caprimulgiformes-Columbiformes node. Some of these relationships differ from those inferred by Sibley and Ahlquist (1990) based on DNA-DNA hybridization. For example, the ß-fibint 7 tree has the Piciformes and Passeriformes as sister orders, whereas the DNA-DNA hybridization tree infers Piciformes to be more basal. The relatively close relationship of Piciformes and Passeriformes and the derived position of the clade is consistent with more traditional classifications (e.g., Wetmore 1960; Cracraft 1981; reviewed by Sibley and Ahlquist 1990).
Within orders, relationships among subtaxa are generally strongly supported and consistent with traditional classifications. Relationships within the Piciformes corroborate current hypotheses (Short 1982; Lanyon and Hall 1994). Similarly, the relationships among passeriform taxa correspond exactly to those inferred by Sibley and Ahlquist (1990), when their tree is pruned to the limited number of species included in our study. Of particular interest here is that several studies support the traditional division of Passeriformes into suboscines (wood pewee) and oscines (other passerines) (Ames 1971 [syrinx]; Edwards, Arctander, and Wilson 1991 [mt cyt b]; Mindell et al. 1997 [mt rRNA]; Garcia-Moreno and Mindell 2000 [CHD]; Lovette and Bermingham 2000 [nuclear c mos]; Irestedt et al. 2001 [c-myc and RAG-1]), but some mtDNA-based phylogenies suggest that these traditional groups are not sister taxa (Mindell et al. 1999; Johnson 2001). The ß-fibint 7 data strongly support this traditional perspective of a sister-group relationship of suboscines and oscines (fig. 3). The subsequent split of the oscines into Corvida (blue jay) and Passerida (other oscines) is strongly supported (fig. 3), consistent with several other studies (Bock 1962 [osteology]; Sheldon and Gill 1996 [DNA hybridization]; Mindell et al. 1997 [mt rRNA]; Pasquet et al. 1999 [mt rRNA]; Ericson, Johansson, and Parsons 2000 [c-myc indels]; Garcia-Moreno and Mindell 2000 [CHD]). Two recent studies (Barker, Barrowclough, and Groth 2002; Ericson et al. 2002) indicate that the suborder Tyranni (suboscines) is paraphyletic, as is the parvorder Corvida, but, again, when the trees in those studies are pruned to the taxa in the ß-fibint 7 tree, the topologies are congruent (assuming the wood pewee represents the Tyrannida and the blue jay represents the "core Corvida" of Barker, Barrowclough, and Groth 2002 and the suboscines and the "Corvida," respectively, of Ericson et al. 2002). The four species of owls divide into two strongly supported clades, corresponding to the widely accepted bifurcation of the owls into two families, Tytonidae (bay owl, barn owl) and Strigidae (flammulated owl, hawk owl). Our sample is deficient in species that would have allowed us to assess the uncertain relationships between owls and frogmouths and their allies (Caprimulgiformes [Sibley and Ahlquist 1990; Bleiweiss, Kirsh, and Lapointe 1994]). In our tree, the frogmouth is united with the Columbiformes as opposed to the Strigiformes (Sibley and Ahlquist 1990), but the bootstrap support is weak (43%), and the MP analysis with indels found six equally most-parsimonious trees distinguished by three placements of the frogmouth in combination with two placements of the ovenbird. Thus, the relationship of the frogmouth and the Caprimulgiformes, to which it is perhaps incorrectly assigned, remains uncertain other than to say this lineage diverged early in the radiation of the modern birds. The relationship depicted in figure 2 is supported, however, by an inferred 4-bp deletion (positions 1096 to 1099) shared by the two doves and the frogmouth. The ß-fibint 7 data strongly support the existence of a sister-group relationship between gallinaceous birds and waterfowl, the Galloanserae (Sibley and Ahlquist 1990; Caspers et al. 1997; Groth and Barrowclough 1999; Mindell et al. 1999; Garcia-Moreno and Mindell 2000; van Tuinen et al. 2000).
Finally, the level of support for ordinal and familial relationships among birds provided by the ß-fibint 7 data set is as strong or stronger than that in other studies reported to date. The number of taxa included in all these studies is small and a limited number of taxa are common to all of the studies; thus, conclusions based on comparisons among the studies are tentative. It is important to note, however, that the substitution rate among diverging ß-fibint 7 lineages is greater than that reported for other nuclear genes. Relative to ß-fibint 7 set at a rate of 1, c-mos, RAG-1, and c-myc evolve at 0.36, 0.30, and 0.25, respectively (based on comparison of genetic distances among taxa held in common among the studies: c-mos, [Lovette and Bermingham 2000], RAG-1 [Groth and Barrowclough 1999], c-myc [Irestedt et al. 2001]). The rate of ß-fibint 7 substitution appears well suited for resolving relationships among lineages that diverged 55 to 90 Myr B.P. and perhaps longer, the timeframe of the Cretaceous-Tertiary boundary, which is so important in vertebrate evolution. Moreover, the indels, themselves, harbor significant phylogenetic information. This is apparent in the indel characters displayed on the internodes in figure 3, the increased bootstrap support when the indel matrix is included in the MP analysis and the fact that an analysis based on the indels alone yields a reasonable tree with several significantly supported nodes. The phylogenetic distribution of indels raises a number of intriguing questions. For example, there is a preponderance of deletions over insertions, which suggests the intron is, on average, becoming shorter through evolution, and the indels appear to have been frequent along some lineages and infrequent or absent along others.
![]() |
Supplementary Materials |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
![]() |
Literature Cited |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Alder, H., M. Yoshinouchi, M. B. Prystowsky, P. Appasamy, and R. Baserga. 1992. A conserved region in intron 1 negatively regulates the expression of the PCNA gene. Nucleic Acids Res. 20:1769-1775.[Abstract]
Ames, P. L. 1971. The morphology of the syrinx in passerine birds. Bull. Peabody Mus. Nat. Hist. 37:1-194.
Barker, F. K., G. F. Barrowclough, and J. G. Groth. 2002. A phylogenetic hypothesis for passerine birds: taxonomic and biogeographical implications of an analysis of nuclear DNA sequence data. Proc. R. Soc. Lond. B 269:295-308.[CrossRef][ISI][Medline]
Benton, M. J. 1990. Phylogeny of the major tetrapod groups: morphological data and divergence dates. J. Mol. Evol. 30:409-424.[ISI][Medline]
Bernardi, G. 1995. The human genome: organization and evolutionary history. Annu. Rev. Genet. 29:445-476.[CrossRef][ISI][Medline]
Bleiweiss, R., J. A. W. Kirsch, and F. J. Lapointe. 1994. DNA-DNA hybridization-based phylogeny for "higher" nonpasserines: reevaluating a key portion of the avian family tree. Mol. Phylogenet. Evol. 3:248-255.[CrossRef][Medline]
Bliss, C. I., and R. A. Fisher. 1953. Fitting the negative binomial distribution to biological data. Biometrics 9:1-25.[ISI]
Bock, W. J. 1962. The pneumatic fossa of the humerus in the Passeres. Auk 79:425-443.
Braun, E. L., and R. T. Kimball. 2002. Examining basal avian divergences with mitochondrial sequences: model complexity, taxon sampling, and sequence length. Syst. Biol. 51:614-625.[CrossRef][ISI][Medline]
Caspers, G. J., D. U. De Weerd, J. Wattel, and W. W. De Jong. 1997. Alpha-crystallin sequences support a galliform/anseriform clade. Mol. Phylogenet. Evol. 7:185-188.[CrossRef][ISI][Medline]
Cooper, A., and D. Penny. 1997. Mass survival of birds across the Cretaceous-Tertiary boundary: molecular evidence. Science 275:1109-1113.
Cracraft, J. 1981. Toward a phylogenetic classification of the recent birds of the world (class Aves). Auk 98:681-714.[ISI]
Cracraft, J. 1988. The major clades of birds. Pp. 339361 in M. J. Benton, ed. The phylogeny and classification of the tetrapods, Vol. 1. Systematics Association Special Volume No. 35A. Clarendon, Oxford.
Cracraft, J., and D. P. Mindell. 1989. The early history of modern birds: a comparsion of molecular and morphological evidence. Pp. 389403 in B. Fernholm, K. Bremer, and H. Jornvall, eds. The hierarchy of life. Elsevier Science Publishers, Amsterdam.
Croft, L., S. Schandorff, F. Clark, K. Burrage, P. Arctander, and J. S. Mattick. 2000. ISIS, the intron information system, reveals the high frequency of alternative splicing in the human genome. Nat. Genet. 24:340-341.[CrossRef][ISI][Medline]
Csank, C., F. M. Taylor, and D. W. Martindale. 1990. Nuclear pre-mRNA introns: analysis and comparsion of intron sequences from Tetrahymena thermophila and other eukaryotes. Nucleic Acids Res. 18:5133-5141.[Abstract]
Deinard, A., and D. G. Smith. 2001. Phylogenetic relationships among the macaques: evidence from the nuclear locus NRAMP1. J. Human Evol. 41:45-59.[CrossRef][ISI][Medline]
Edwards, S. V., P. Arctander, and A. C. Wilson. 1991. Mitochondrial resolution of a deep branch in the genealogical tree for perching birds. Proc. R. Soc. Lond. B 243:99-107.[ISI][Medline]
Edwards, S. V., B. Fertil, A. Giron, and P. J. Deschavanne. 2002. A genomic schism in birds revealed by phylogenetic analysis of DNA strings. Syst. Biol. 51:599-613.[CrossRef][ISI][Medline]
Ericson, P. G. P., L. Christidis, A. Cooper, M. Irestedt, J. Jackson, U. S. Johansson, and J. Norman. 2002. A Gondwanan origin of passerine birds supported by DNA sequences of the endemic New Zealand wrens. Proc. R. Soc. London B 269:235-241.[CrossRef][ISI][Medline]
Ericson, P. G. P., U. S. Johansson, and T. J. Parsons. 2000. Major divisions in oscines revealed by insertions in the nuclear gene c-myc: a novel gene in avian phylogenetics. Auk 117:1069-1078.[ISI]
Feduccia, A. 1995. Explosive evolution in Tertiary birds and mammals. Science 267:637-638.[ISI]
Feduccia, A. 1996. The origin and evolution of birds. Yale University Press, New Haven.
Frati, F., C. Simon, J. Sullivan, and D. L. Swofford. 1997. Evolution of the mitochondrial cytochrome oxidase II gene in collembola. J. Mol. Evol. 44:145-158.[ISI][Medline]
Garcia-Moreno, J., and D. P. Mindell. 2000. Rooting a phylogeny with homologous genes on opposite sex chromosomes (gametologs): a case study using avian CHD. Mol. Biol. Evol. 17:1826-1832.
Gasch, A., U. Hinz, and R. Renkowitz-Pohl. 1989. Intron and upstream sequences regulate expression of the Drosophila ß3-tubulin gene in the visceral and somatic musculature, respectively. Proc. Natl. Acad. Sci. USA 86:3215-3218.[Abstract]
Giannasi, N., A. Malhotra, and R. S. Thorpe. 2001. Nuclear and mtDNA phylogenies of the Trimeresurus complex: implications for the gene versus species tree debate. Mol. Phylogenet. Evol. 19:57-66.[CrossRef][ISI][Medline]
Groth, J. G., and G. F. Barrowclough. 1999. Basal divergences in birds and the phylogenetic utility of the nuclear RAG-1 gene. Mol. Phylogenet. Evol. 12:115-123.[CrossRef][ISI][Medline]
Hall, B. G. 2001. Phylogenetic trees made easy: a how-to manual for molecular biologists. Sinuaer Associates, Sunderland, Mass.
Harlid, A., A. Janke, and U. Arnason. 1997. The mtDNA sequence of the ostrich and the divergence between paleognathous and neognathous birds. Mol. Biol. Evol. 14:754-761.[Abstract]
Harlid, A., A. Janke, and U. Arnason. 1998. The complete mitochondrial genome of Rhea americana and early avian divergences. J. Mol. Evol. 46:669-679.[ISI][Medline]
Harshman, J. 1994. Reweaving the tapestry: what can we learn from Sibley and Ahlquist (1990)? Auk 111:377-388.[ISI]
Hedges, S. B. 1994. Molecular evidence for the origin of birds. Proc. Natl. Acad. Sci. USA 91:2621-2624.[Abstract]
Hughes, A. L., and M. Yeager. 1997. Comparative evolutionary rates of introns and exons in murine rodents. J. Mol. Evol. 45:125-130.[ISI][Medline]
Irestedt, M., U. S. Johansson, T. J. Parsons, and P. G. P. Ericson. 2001. Phylogeny of major lineages of suboscines (Passeriformes) analysed by nuclear DNA sequence data. J. Avian Biol. 32:15-25.[CrossRef][ISI][Medline]
Jackson, P. D., and F. M. Hoffmann. 1994. Embryonic expression patterns of the Drosophila decapentaplegic gene: separate regulatory elements control blastoderm expression and lateral ecotodermal expression. Develop. Dynamics 199:28-44.[ISI][Medline]
Johnson, K. P. 2001. Taxon sampling and the phylogenetic position of Passeriformes: evidence from 916 avian cytochrome b sequences. Syst. Biol. 50:128-136.[CrossRef][ISI][Medline]
Johnson, K. P., and D. H. Clayton. 2000. Nuclear and mitochondrial genes contain similar phylogenetic signal for pigeons and doves (Aves: Columbiformes). Mol. Phylogenet. Evol. 14:141-151.[CrossRef][ISI][Medline]
Keller, E. B., and W. A. Noon. 1985. Intron splicing: a conserved internal control signal in introns of Drosophila pre-mRNAs. Nucleic Acids Res. 13:4971-4981.[Abstract]
Kirby, A., S. V. Muse, and W. Stephan. 1995. Maintenance of pre-mRNA secondary structure by epistatic selection. Proc. Natl. Acad. Sci. USA 92:9047-9051.[Abstract]
Koop, B. F., D. A. Tagle, M. Goodman, and J. L. Slightom. 1989. A molecular view of primate phylogeny and important systematic and evolutionary questions. Mol. Biol. Evol. 6:580-612.[Abstract]
Kumar, S., K. Tamura, and M. Nei. 1993. Mega: molecular evolutionary genetics analysis. Version 1.01. University Park, Penn.
Lanyon, S. M. 1988. The stochastic mode of molecular evolution: what consequences for systematic investigations. Auk 105:565-573.[ISI]
Lanyon, S. M. 1992. Review of "Phylogeny and classification of birds. A study in molecular evolution," by C.G. Sibley and J.E. Ahlquist. Condor 94:304-310.
Lanyon, S. M., and J. G. Hall. 1994. Reexamination of barbet monophyly using mitochondrial DNA sequence data. Auk 111:389-397.[ISI]
Leicht, B. G., E. M. S. Lyckegaard, C. M. Benedict, and A. G. Clark. 1993. Conservation of alternative splicing and genomic organization of the myosin alkali light-chain (Mlc1) gene among Drosophila species. Mol. Biol. Evol. 10:769.[Abstract]
Leicht, B. G., S. V. Muse, M. Hanczyc, and A. G. Clark. 1995. Constraints on intron evolution in the gene encoding the myosin alkali light chain in Drosophila. Genetics 139:299-308.
Lewin, B. 1997. Genes VI. Oxford University Press, Oxford, U.K.
Li, W.-H. 1997. Molecular evolution. Sinauer Associates, Sunderland, Mass.
Lockhart, P. J., M. A. Steel, M. D. Hendy, and D. Penny. 1994. Recovering evolutionary trees under a more realistic model of sequence evolution. Mol. Biol. Evol. 11:605-612.
Lovette, I. J., and E. Bermingham. 2000. c-mos variation in songbirds: molecular evolution, phylogenetic implciations, and comparisons with mitochondrial differentiation. Mol. Biol. Evol. 17:1569-1577.
Mattick, J. S., and M. J. Gagen. 2001. The evolution of controlled multitasked gene networks: the role of introns and other noncoding RNAs in the development of complex organisms. Mol. Biol. Evol. 18:1611-1630.
Metz, A., S. Knoechel, P. Buechler, M. Koester, and W. Knoechel. 1998. Structural and functional analysis of the BMP-4 promoter in early embryos of Xenopus laevis. Mech. Develop. 74:29-39.[CrossRef][ISI][Medline]
Miller, K. M., and R. E. Withler. 1996. Sequence analysis of a polymorphic Mhc class II gene in Pacific salmon. Immunogenetics 43:337-351.[CrossRef][ISI][Medline]
Mindell, D. P. 1992. DNA-DNA hybridization and avian phylogeny. Syst. Biol. 41:126-134.
Mindell, D. P., M. D. Sorenson, D. E. Dimcheff, and M. Hasegawa. 1999. Interordinal relationships of birds and other reptiles based on whole mitochondrial genomes. Syst. Biol. 48:138-152.[CrossRef][ISI][Medline]
Mindell, D. P., M. D. Sorenson, C. J. Huddleston, H. C. Miranda, A. Knight, S. J. Sawchuk, and T. Yuri. 1997. Phylogenetic relationships among and within select avian orders based on mitochondrial DNA. Pp. 214250 in D. P. Mindell, ed. Avian molecular evolution and systematics. Academic Press, San Diego.
Moore, W. S., and V. R. DeFilippis. 1997. The window of taxonomic resolution for avian phylogenies based on mitochondrial cytochrome b DNA sequences. Pp. 83119 in D. P. Mindell, ed. Avian molecular evolution and systematics. Academic Press, San Diego.
Moore, W. S., S. M. Smith, and T. Prychitko. 1999. Nuclear gene introns versus mitochondrial genes as molecular clocks. Pp. 745753 in N. Adams and R. Slotow, eds. Proc. 22 Int. Ornithol. Congr. University of Natal, Durbin.
Mount, S. M., C. Burks, G. Hertz, G. D. Stormo, O. White, and C. Fields. 1992. Splicing signals in Drosophila intron size information content and consensus sequences. Nucleic Acids Res. 20:4255-4262.[Abstract]
Pasquet, E., A. Cibois, F. Baillon, and C. Erard. 1999. Relationships between the ant-thrushes Neocossyphus and the flycatcher-thrushes Stizorhina, and their position relative to Myadestes, Entomodestes and some other Turdidae (Passeriformes). J. Zool. Syst. Evol. Res. 37:177-183.[ISI]
Paton, T., O. Haddrath, and A. J. Baker. 2002. Complete mitochondrial DNA genome sequences show that modern birds are not descended from transitional shorebirds. Proc. R. Soc. Biol. Sci. Ser. B 269:839-846.[CrossRef][ISI][Medline]
Prychitko, T. M., and W. S. Moore. 1997. The utility of DNA sequences of an intron from the ß-fibrinogen gene in phylogenetic analysis of woodpeckers (Aves: Picidae). Mol. Phylogenet. Evol. 8:193-204.[CrossRef][ISI][Medline]
Prychitko, T. M., and W. S. Moore. 2000. Comparative evolution of the mitochondrial cytochrome b gene and nuclear ß-fibrinogen intron 7 in woodpeckers. Mol. Biol. Evol. 17:1101-1111.
Prychitko, T. M., E. A. Ries, and W. S. Moore. 1998. Two-step cycle sequencing reduces premature terminations when using primers with high annealing temperatures. Mol. Biotechnology 10:231-236.[ISI][Medline]
Pycraft, W. P. 1900. The morhpology and phylogeny of the Paleognathae (Ratitae and Crypturi) and the Neognathae (Carinatae). Trans. Zool. Soc. Lond. 15:149-290.
Rodriguez-Trelles, F., R. Tarrio, and F. J. Ayala. 2000. Disparate evolution of paralogous introns in the Xdh gene of Drosophila. J. Mol. Evol. 50:123-130.[ISI][Medline]
Ruskin, B., A. R. Krainer, T. Maniatas, and M. R. Green. 1984. Excision of an intact intron as a novel lariat structure during pre-mRNA splicing in vitro. Cell 38:317-331.[ISI][Medline]
Schneider, H., I. Sampaio, M. L. Harada, C. M. L. Barroso, M. P. C. Schneider, J. Czelusniak, and M. Goodman. 1996. Molecular phylogeny of the New World monkeys (Platyrrhini, Primates) based on two unlinked nuclear genes: IRBP intron 1 and epsilon-globin sequences. Am. J. Phys. Anthro. 100:153-179.[CrossRef][ISI][Medline]
Schultz, J. R., T. Tansey, L. Gremke, and R. V. Storti. 1991. A muscle-specific Intron enhancer required for rescue of indirect flight muscle and jump muscle function regulates Drosophila tropomyosin I gene expression. Mol. Cell. Biol. 11:1901-1911.[ISI][Medline]
Shapiro, L. H., and J. P. Dumbacher. 2001. Adenylate kinase intron 5: A new nuclear locus for avian systematics. Auk 118:248-255.[ISI]
Sheldon, F. H., and F. B. Gill. 1996. A reconsideration of songbird phylogeny, with emphasis on the evolution of titmice and their sylvioid relatives. Syst. Biol. 45:473-495.[ISI]
Short, L. L. 1982. Woodpeckers of the world. Delaware Museum of Natural History, Greenville, Del.
Sibley, C. G., and J. E. Ahlquist. 1990. Phylogeny and classification of birds: a study in molecular evolution. Yale University Press, New Haven, Conn.
Sibley, C. G., and J. B. L. Monroe. 1990. Distribution and taxonomy of birds of the world. Yale University Press, New Haven, Conn.
Slade, R. W., C. Moritz, and A. Heideman. 1994. Multiple nuclear-gene phylogenies: application to pinnipeds and comparison with a mitochondrial DNA gene phylogeny. Mol. Biol. Evol. 11:341-356.[Abstract]
Swofford, D. L. 1998. PAUP*: phylogenetic analysis using parsimony (*and other methods). Version 4. Sinauer Associates, Sunderland, Mass.
Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 24:4867-4882.
van Tuinen, M., and S. B. Hedges. 2001. Calibration of avian molecular clocks. Mol. Biol. Evol. 18:206-213.
van Tuinen, M., C. G. Sibley, and B. S. Hedges. 2000. The early history of modern birds inferred from DNA sequences of nuclear and mitochondrial ribosomal genes. Mol. Biol. Evol. 17:451-457.
Waddell, P. J., Y. Cao, M. Hasegawa, and D. P. Mindell. 1999. Assessing the Cretaceous superordinal divergence times within birds and placental mammals by using whole mitochondrial protein sequences and an extended statistical framework. Syst. Biol. 48:119-137.[CrossRef][ISI][Medline]
Walton, A. H., M. A. Nedbal, and R. L. Honeycutt. 2000. Evidence from intron 1 of the nuclear transthyretin (Prealbumin) gene for the phylogeny of African mole-rats (Bathyergidae). Mol. Phylogenet. Evol. 16:467-474.[CrossRef][ISI][Medline]
Weibel, A. C., and W. S. Moore. 2002. A test of a mitochondrial gene-based phylogeny of woodpeckers (genus Picoides) using an independent nuclear gene, ß-fibrinogen intron 7. Mol. Phylogenet. Evol. 22:247-257.[CrossRef][ISI][Medline]
Wetmore, A. 1960. A classification for the birds of the world. Smithson. Misc. Collect. 39:1-37.
Yang, Z. 1994. Estimating the pattern of nucleotide substitution. J. Mol. Evol. 39:105-111.[ISI][Medline]
Yang, Z., N. Goldman, and A. Friday. 1995. Maximum-likelihood trees from DNA sequencesa peculiar statistical estimation problem. Syst. Biol. 44:384-399.[ISI]