| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
,1
* Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana 70803
1 Corresponding author: Museum of Natural Science, 119 Foster Hall, Baton Rouge, LA 70803.
E-mail: mcarli1{at}lsu.edu
 |
ABSTRACT |
|---|
 TOP ABSTRACT MATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: SAMPLING LOCALITIES...ACKNOWLEDGEMENTSLITERATURE CITED |
|---|
Over most of its taxonomic history, the genus Passerina (Aves: Cardinalidae) has been composed of six species of small (13–20 g), sexually dichromatic songbirds. The collective breeding ranges of these species encompass most of Mexico, the United States, and southern Canada (Figure 1 and supplemental Figure 1 at http://www.genetics.org/supplemental/). Within the genus, Passerina cyanea and P. amoena have typically been considered sister species because they are morphologically similar and form a well-known hybrid zone where their breeding distributions overlap in the western Great Plains and eastern foothills of the Rocky Mountains in North America (SIBLEY and SHORT 1959; EMLEN et al. 1975; KROODSMA 1975; BAKER and BAKER 1990; BAKER and BOYLAN 1999). While the monophyly of Passerina has never been questioned (HELLACK and SCHNELL 1977; TAMPLIN et al. 1993), there have been attempts to subsume Guiraca caerulea and some species of the genus Cyanocompsa into the genus (PHILLIPS et al. 1964; BLAKE 1969; MAYR and SHORT 1970; PAYNTER 1970).
|
The branching structure of the mitochondrial tree and the demography of Passerina buntings suggest that the unexpected P. amoena/P. caerulea sister relationship could reflect discordance between the mitochondrial gene tree and the species tree. The internodes in the mitochondrial phylogeny are short, possibly indicating rapid divergences, and the population sizes for some of the species, including P. cyanea, P. amoena, and P. caerulea, are large—on the order of tens of millions of individuals (PAYNE 1992; GREENE et al. 1996). This combination of short internodes and large population sizes is among the most difficult scenarios for inferring phylogenies from single locus data sets given the elevated risk of discordance between a species tree and its underlying gene trees (Figure 2). Using the calculation of HEIN et al. (2005) for the probability of species tree/gene tree discordance and realistic parameter estimates for Passerina buntings, the probability of discordance is estimated to be as high as 55% (see Comparing mtDNA and nuclear DNA signals in DISCUSSION).
|
|
MATERIALS AND METHODS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: SAMPLING LOCALITIES...ACKNOWLEDGEMENTSLITERATURE CITED |
|---|
25 mg of pectoral muscle using a DNeasy tissue kit (QIAGEN, Valencia, CA), following the manufacturer's recommended protocol. Each individual was amplified at 10 nuclear loci (Table 1) using the following PCR conditions in a 25 µl amplification reaction: 50 ng template DNA (2.5 µl of QIAGEN DNA extracts), 1 µl of 10 mM dNTPs (2.5 mM each dATP, dTTP, dCTP, and dGTP), 1 µl of each primer (10 mM, Table 1), 2.5 µl 10x buffer with MgCl2 (15 mM), 0.1 µl Taq (5 units/µl AmpliTaq DNA polymerase, Applied Biosystems, Foster City, CA), and 16.9 µl sterile dH2O. The thermocycling profile was as follows: an initial 95° denaturation for 2 min followed by 35 cycles consisting of a 30-sec, 95°-denaturation step, a 30-sec, locus-specific temperature primer annealing step (Table 1), and a 2-min, 72°-extension step, and a final extension of 5 min at 72°. We checked for amplification by electrophoresing 2.5 µl of each PCR amplicon on a 1% agarose gel and then used 20% polyethylene glycol (PEG) to purify the remaining 22.5 µl of the PCR product from successful amplifications.
|
Phylogenetic analyses:
Sequence alignments of all 10 nuclear loci were unambiguous and checked visually. In all phylogenetic analyses, gaps were alternately classified as a fifth state or as missing data. Phylogenetic analyses were conducted on a concatenated data set since partition homogeneity tests revealed no significant conflicts between the data from each of the individual loci. Tests for recombination (see Population genetic analyses below) indicated that three loci (βact3, MPLPR, and TGFβ2) had undergone recombination. Therefore, phylogenetic analyses were conducted on both the entire concatenated data set (5232 bp) and on a reduced concatenated data set containing only the largest independently segregating block of sequence for βact3, MPLPR, and TGFβ2 (4668 bp). Using PAUP* v 4.0b10 (SWOFFORD 1998) we estimated a neighbor-joining tree from the p-distance matrix and calculated likelihood scores on a series of nested substitution models. These were then evaluated using ModelTest v 3.8 (POSADA and CRANDALL 1998). Due to its advantages over hierarchical likelihood-ratio tests (POSADA and BUCKLEY 2004), we used the Bayesian Information Criterion (BIC) to determine the best-fit model of sequence evolution for the concatenated data set as well as for each locus independently (Table 2).
|
Using parameters estimated on a neighbor-joining tree (SAITOU and NEI 1987), we used PAUP* v 4.0b10 (SWOFFORD 1998) to perform heuristic searches [tree bisection–reconnection (TBR) branch swapping; 10 random-addition replicates] under the maximum-likelihood (ML) optimality criterion for both the 8- and 16-terminal-tip concatenated data sets. We assessed nodal support using 100 ML heuristic bootstrap replicates with the search parameters described above. We also performed a variety of constrained ML searches (Table 3) on the 8- and 16-terminal-tip data sets using the sequence evolution model and search parameters, except without bootstrap replicates, as described above. Using a Shimodaira–Hasegawa test (SH) (SHIMODAIRA and HASEGAWA 1999) with full optimization and 1000 bootstrap replicates, we assessed whether the likelihoods of the unconstrained and constrained trees were significantly different.
|
Minimize-deep-coalescences analyses:
Because well-supported species trees inferred from concatenated data sets can be discordant with the true species tree, especially if the species divergence was recent (DEGNAN and ROSENBERG 2006; EDWARDS et al. 2007; KUBATKO and DEGNAN 2007), we inferred a species tree using a method that considers the individual gene histories. Specifically, we assessed the fit of the gene trees onto three possible rooted species trees for P. cyanea, P. amoena, and P. caerulea, using the minimize-deep-coalescences method as implemented in the program Mesquite (MADDISON 1997; MADDISON and KNOWLES 2006). This method attempts to reconcile gene trees contained within species trees by minimizing the number of deep coalescences, those cases where the coalescence of two gene copies predates the speciation event. Using a parsimony optimality criterion, the "best" species tree is the one that requires the fewest deep-coalescent events in contained gene genealogies. The number of deep-coalescent events, or costs, can be summed across multiple markers to accommodate multilocus data sets. Using Mesquite (MADDISON and MADDISON 2004), we compared the fit of the genealogies from each locus into the possible species trees using methods slightly modified from the recently described method of MADDISON and KNOWLES (2006).
We used MrBayes v 3.1.2 (HUELSENBECK and RONQUIST 2001; RONQUIST and HUELSENBECK 2003) to sample the tree space of the gene trees for each locus independently, using completely phased haplotypes and the same model parameters with uninformative priors as in the ML tree searches and the same MrBayes run parameters as described above. After discarding the first 250,000 generations as burn-in, we used PAUP* to calculate a majority-rule consensus gene tree from each locus-specific search. We then overlaid each unrooted locus-specific gene tree onto the three possible species trees (Table 3) and tallied the number of deep- coalescent events required to reconcile the gene trees with the species trees. The total cost of each species tree was determined by summing the locus-specific tallies (Table 4), with the species tree having the lowest cost being inferred as the correct species tree.
|
 |
In light of the current methodological limitations, we focused on calculating the conditional probability of a gene tree given a species tree. Locus-specific gene trees were inferred in three ways. First, we randomly chose a single individual to represent P. amoena, P. cyanea, and P. caerulea, which along with one individual from each outgroup taxon (P. rositae and C. cyanoides), resulted in a data set with five terminal tips. For each locus, we estimated a gene tree under the maximum-likelihood search criterion as described above and then used COAL to calculate the probability of each ML gene tree given three possible species trees (Table 3) using COAL. ML inferences for Tropo6 produced >100 equally likely gene trees and so we eliminated Tropo6 from the COAL analysis.
We also inferred locus-specific gene trees from the 8-terminal-tip data set. Since COAL requires fully resolved bifurcating gene trees and since locus-specific ML searches of this 8-terminal-tip data set produced numerous equally likely topologies with unresolved polytomies, we inferred locus-specific neighbor-joining gene trees, using uncorrected p-distances. In the Tropo6 data set, some of the sequences were identical across individuals. Accordingly, neighbor-joining analyses failed to produce a fully resolved tree for Tropo6 and we excluded it from COAL calculations.
In an attempt to incorporate a larger sample of the possible gene topologies and branch lengths, we sampled a collection of gene trees from the locus-specific Bayesian tree searches described above (see Minimize-deep-coalescences analyses). Our approach was to calculate the probability of each of the three possible species trees (Table 3) given each of the last 1000 gene trees from the Bayesian searches, using COAL, and average the probabilities of those gene trees to generate a more robust estimate of the likelihoods of the three species trees. Unfortunately, in each locus-specific Bayesian tree search some gene trees contained unresolved polytomies, which COAL cannot handle; therefore, we abandoned this method of inferring gene trees.
We replicated the COAL calculations, for both the five- and eight-terminal tip data sets, using three different internal:terminal branch-length ratios (1:1, 1:2, and 1:10). All internal:terminal branch-length ratios, with the exception of 1:1, correspond to topologies with internode lengths that are short, relative to the terminal branches, which appears to be the general pattern observed in Passerina buntings. To assess the difference in probabilities of the gene trees if the species tree united P. cyanea and P. amoena as sisters (hypothesis 1) with one in which P. caerulea and P. amoena were sisters (hypothesis 2), we conducted a test of proportions (SOKAL and ROHLF 1995) whereby we tallied the number of gene trees that were more probable under hypothesis 1 and compared that with a tally of the number of gene trees that were more probable under hypothesis 2.
Population genetic analyses:
Both phylogenetic methods described above assume that discordance between gene trees and the species tree is due to genetic drift (i.e., lineage sorting stochasticity). Because of the well-documented hybrid zone between P. amoena and P. cyanea, it would be ideal if a gene flow parameter could be incorporated into the phylogenetic analysis. This is not currently possible. We thus examined the utility of divergence population genetics for inferring species relationships indirectly. The Isolation with Migration (IM) software package (HEY and NIELSEN 2004; HEY 2005) can estimate pairwise divergence times between P. cyanea, P. amoena, and P. caerulea, while jointly estimating pairwise migration rates. A sister relationship is inferred to be between the species with the most recent divergence estimate.
IM combines coalescent theory with Bayesian methodologies to simultaneously estimate multiple population genetic parameters for two diverging populations. Briefly, using a Markov chain Monte Carlo (MCMC) approach, IM generates posterior probability distributions of six parameters (
population1, population2, ancestral, m1, m2, t). These parameters can then be used to calculate population demographic estimates, i.e., effective population sizes of the two populations (Ne1 and Ne2), the ancestral population size (i.e., the population size prior to initial divergence = NeA), migration rates between the two populations (m1 and m2, which allows for asymmetrical migration rates), and the time since divergence of the two populations (t). A more recent version of IM (HEY 2005) also allows for the inclusion of a population splitting parameter, s, which allows for exponential population growth or decline in each descendant population.
IM assumes selective neutrality of each marker and no intralocus recombination. We tested for neutrality by calculating Tajima's D (TAJIMA 1983) for each locus using Arlequin v 3.0 (EXCOFFIER et al. 2005). None of the 10 loci showed a significant departure from neutrality (each P > 0.25). To investigate intralocus recombination, we used TOPALi v 2 (MCGUIRE and WRIGHT 1998, 2000), which uses a sliding-window approach to search for putative recombination breakpoints along an alignment of DNA sequences. Using the difference in sum-of-squares (DSS) method, with a 100-bp window and a 10-bp increment, we looked for signals of recombination, manifested by a significant DSS peak along the sequence alignment. The statistical significance of DSS peaks was assessed using 100 parametric bootstraps. Significant DSS peaks were detected in 3 loci, βact3, MPLPR, and TGFβ2 (Table 1). For each of those 3 loci, we only included the longest independently segregating block of sequence from the phased haplotype data in the IM runs along with the total available phased haplotype sequences from each of the other seven markers.
Using IM, we estimated three pairwise divergence times, one between P. cyanea and P. amoena, one between P. caerulea and P. amoena, and one between P. cyanea and P. caerulea. Initial IM runs, conducted using the infinite sites (IS) model (KIMURA 1969) and employing unrealistically large maximum priors (t = 20, q1 = q2 = qa = 500, m1 = m2 = 50, or m1 = m2 = 0 for parameter estimation without gene flow) were used to establish appropriate parameter priors for all subsequent analyses (WON and HEY 2005). Once appropriate priors were identified, we ran IM three times with identical starting conditions, with the exception of the random number seed, to assess convergence. All runs began with a burn-in period of 100,000 steps and were allowed to continue until the minimum effective sample size (ESS) was >100 (HEY 2005). Because runs with the same starting conditions produced qualitatively similar results, we only present data from the longest run (>3 x 109 steps). For each pairwise estimate of divergence time, we performed a variety of IM analyses with varying numbers of parameters. The simplest model included four parameters (1, 2, A, and t); other models included five (1, 2, A, t, and m—equal migration between populations or 1, 2, A, t, and s—allowing for population growth) or six parameters (1, 2, A, t, m1, and m2—allowing for different, directional migration between populations). All model parameter estimates in IM are scaled to the mutation rate, µ. To translate scaled divergence time estimates to real time, we included mutation rates for each locus calculated as the sequence length multiplied by 3.6 x 10–9 substitutions/site/year (AXELSSON et al. 2004).
|
RESULTS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: SAMPLING LOCALITIES...ACKNOWLEDGEMENTSLITERATURE CITED |
|---|
|
The topologies inferred from the tree searches of the concatenated data sets did not differ when gaps were coded as missing data or as a fifth state (topologies not shown).
Minimize-deep-coalescences analyses:
Overlaying the majority-rule consensus gene trees from the Bayesian tree searches on the three possible species trees resulted in three different deep coalescent costs for each gene—one cost for each species tree (Table 4). For four locus-specific gene trees (AETC, Rho1, TGFβ2, and Tropo6) the fewest number of deep coalescent events were observed when the species topology united P. amoena and P. cyanea as sister species. The fewest deep coalescent events for one locus (MC1R) occurred when P. amoena and P. caerulea were sister species, and the P. cyanea/P. caerulea sister species tree contained the fewest deep coalescences when overlaid with the ODC and MPLPR gene trees.
Summing the number of deep coalescent events needed to fit the gene trees into the species tree across loci produced a total cost of 42 for the P. amoena/P. cyanea sister species tree (Table 4). Fitting the gene trees into the other two species trees required more deep coalescences: 46 for the P. amoena/P. caerulea tree and 45 for the P. cyanea/P. caerulea tree.
Calculating probabilities of gene trees given possible species trees:
After randomly choosing one individual to represent P. amoena, P. cyanea, and P. caerulea, and eliminating Tropo6 from further analyses, 14 ML gene topologies were inferred from phylogenetic tree searches conducted on each locus separately (Table 5). Multiple gene trees with equal likelihoods were inferred from three loci (AETC—three gene trees, βact3—two, and MPLPR—three). Replicate COAL calculations that differed only in the ratio of internal to terminal branch lengths (1:1, 1:2, and 1:10) all gave qualitatively similar results. When terminal branch lengths were longer than internal branch lengths (1:2 and 1:10), 3 of the 14 gene topologies were equally probable whether the given species tree had P. cyanea and P. amoena as sister species or a sister species relationship between P. amoena and P. cyanea (Table 5). Those 3 gene topologies (one of two βact3 gene trees and two of three MPLPR gene trees) were excluded from further analyses, because with respect to these relationships, these topologies are uninformative. Of the remaining 9 ML gene topologies (one of the three equally likely AETC gene trees was chosen at random), 8 were more probable given a species tree uniting P. cyanea and P. amoena as sisters, and 1 was more probable given a P. amoena–P. caerulea species tree (Table 5, test of proportions, nobs = 9, P = 0.05). When the internal:terminal branch length ratio was 1:1, 4 of the 14 gene topologies were equally probable given either species tree and were subsequently eliminated. Seven of the remaining 8 (1 AETC gene tree chosen at random) ML gene topologies were more probable with P. cyanea/P. amoena as sister species (Table 4, test of proportions, nobs = 8, P = 0.1).
|
|
1, 2 and A were inconclusive; the distributions of the posterior probabilities were nearly flat in all comparisons (data not shown). Allowing for growth or decline in the descendant populations did not qualitatively alter the divergence time estimates for any of the pairwise comparisons; the most recent divergence was between P. amoena and P. cyanea, with divergence times later for both the P. caerulea/P. amoena and P. caerulea/P. cyanea comparisons (Table 7).
|
|
Estimates of migration rates between the species varied. The posterior probability distributions of m1 and m2 between P. cyanea and P. amoena were flat, but nonzero (Figure 4, C and D). In contrast, migration into P. caerulea from both P. cyanea and P. amoena was effectively zero (m2 in both comparisons, Figure 4D). However, the parameter distributions for migration into P. cyanea and P. amoena from P. caerulea both show evidence for nonzero gene flow (Figure 4C). The posterior probability distributions for m1 were also similar for both comparisons involving P. caerulea, one peak at zero and a secondary, nonzero peak indicating gene flow into either P. cyanea or P. amoena from P. caerulea (Figure 4C).
While the divergence time parameter estimates did change with the inclusion of the gene flow parameter (Figure 4 and Table 7), the data were not significantly more likely under the more complex model in any of the pairwise comparisons. In all three comparisons twice the difference between the likelihoods of the data under the complex model (allowing for directional migration) and the simple model (no migration) were less than the critical value (
= 5.991, d.f. = 2, P = 0.05)
|
DISCUSSION |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: SAMPLING LOCALITIES...ACKNOWLEDGEMENTSLITERATURE CITED |
|---|
). Mutations in MC1R are known to correlate with variations in melanin-based plumage in a variety of taxa, including birds (MUNDY 2005; HOEKSTRA 2006). Within Passerina, P. amoena and P. caerulea are the only species with putatively melanin-based plumages (HILL and MCGRAW 2006) and selection acting on MC1R might explain the differences we found between inferences based on MC1R and those based on other loci. There is a nonsynonymous change at amino acid residue 102 (Val 
Leu) that unites P. caerulea and P. amoena (Val) to the exclusion of P. cyanea (Leu), which does not have any melanistic plumage. However, P. rositae and C. cyanoides (neither of which have melanin-based plumage) share the Val allele with P. amoena and P. caerulea indicating that the allele (Leu) found in P. cyanea is likely derived and probably not related to the differences in plumage patterns between P. caerulea, P. amoena, and P. cyanea. Therefore, while we cannot completely rule out the possibility that selection on the melanin-based plumages of P. amoena and P. caerulea has resulted in the observed patterns, the available data suggest selection is an unlikely cause. Below, we discuss the benefits of using the two types of analyses for resolving sister relationships in Passerina and then with respect to closely related taxa in general.
Comparing mtDNA and nuclear DNA signals:
The sister relationship between P. cyanea and P. amoena supported by the multilocus analyses conflicts with a well-supported P. caerulea/P. amoena sister relationship inferred from mitochondrial data (KLICKA et al. 2001). A possible explanation for the discrepancies between the mtDNA and nuclear results is differential lineage sorting of the mtDNA with respect to the true species tree (BRAUN and BRUMFIELD 1998). Assuming ((P. cyanea, P. amoena), P. caerulea) is the correct topology, one can estimate the probability of incongruence between a gene tree and the species tree (FELSENSTEIN 2004; HEIN et al. 2005). Incorporating a range of internode values based on the divergence time estimates from IM and a range of potential effective ancestral population sizes, the chance a gene tree differs from the true species tree in the Passerina system ranges from 0% (Ne = 10,000 individuals, internode length = 100,000 years) to 65% (ancestral Ne = 1,000,000, internode length = 100,000 years, Figure 2). Using an internode length of 350,000 years, which is the difference in the maximum-likelihood estimates of the divergence time for the P. cyanea x P. amoena (138,000 years ago) and the P. caerulea x P. amoena (486,000 years ago) comparisons (Table 7), and assuming no postspeciation gene flow, the probabilities of incongruence range from 0% (ancestral Ne = 10,000) to 55% (ancestral Ne = 1,000,000) (Figure 2). While an accurate estimate of the ancestral population size is unavailable, an effective size of 10,000 or 20,000 may be unrealistically small. Current population census estimates for P. cyanea are on the order of 20–80 million individuals (PAYNE 1992). Furthermore, while some coalescent-based ancestral population-size estimates tend to be relatively small (HEY 2005; WON and HEY 2005), an estimate in birds based on sequence data from multiple nuclear markers suggested an ancestral population size on the order of hundreds of thousands of individuals (JENNINGS and EDWARDS 2005). Large population sizes and short internode lengths combine to increase the chance of incomplete and incorrect lineage sorting of gene copies, confusing the phylogenetic relationships of the focal taxa. The analyses of a multilocus data set presented herein suggest that the stochastic nature of lineage sorting of the mitochondrial genome resulted in a gene tree that is discordant with the species tree.
In closely related taxa, gene flow, in addition to lineage sorting problems, may limit the ability to correctly infer the species tree from gene sequences. P. cyanea and P. amoena currently hybridize where their breeding ranges overlap, although the amount of introgression is currently unknown, and hybrid specimens between other members of the genus have been reported. The samples of P. amoena and P. cyanea included in this study were collected far from the Great Plains contact zone (APPENDIX) to minimize the chance that recent introgression between the two species might be obscuring the true evolutionary relationships between P. cyanea, P. amoena, and P. caerulea. Ongoing studies of the genetic structure of the hybrid zone between P. cyanea and P. amoena suggest little mitochondrial or nuclear introgression is occurring (M. CARLING, unpublished data). The data from the hybrid zone, along with the nonsignificant difference in the likelihood of the data to the more complex model in IM, indicate that gene flow between P. cyanea and P. amoena has probably not had a strong influence on the speciational history inferred in this study.
Implications for uncovering phylogenetic signals in closely related taxa:
Integrating top-down phylogenetic methods with bottom-up population genetic methods was successful in that the same sister relationship was supported by both types of analyses. The data set of 10 loci sampled from either one (P. rositae and C. cyanoides) or two individuals per species (P. cyanea, P. amoena, and P. caerulea) generated a relatively well-supported species tree despite the low level of variation at each locus. One possible solution to increase the resolution of individual gene trees is to increase the length of the sampled sequence. Theoretically, lengthening the sampled sequence should increase the chance of detecting variable sites, which may then offer additional support for a particular topology (DOYLE 1992). In practice however, lengthening loci may not be a straightforward solution to a problem of poor resolution in phylogeny reconstruction. As the length of the sampled locus increases, so does the chance of sampling sites that have undergone recombination. We detected recombination in 3 of the 10 loci, a ratio similar to a recent multilocus study investigating the speciational history of Australian grass finches (JENNINGS and EDWARDS 2005). Both studies sampled loci between 200 and 900 bp, which demonstrates that potential effects of recombination should be considered in any population genetic or phylogenetic study employing sequence data from nuclear loci. Here, the results of the phylogenetic analyses were similar regardless of whether we analyzed the entire sequence length or just the longest independently segregating block of sequence from the 3 loci (βact3, MPLPR, and TGFβ2) that showed evidence of recombination (Table 3, Figure 3).
The sampling design used here closely mimicked the design that resulted in the greatest accuracy of species tree inference in MADDISON and KNOWLES'S 2006 investigation into the utility of using the minimize-deep-coalescences method for phylogenetic reconstruction. Using simulated data sets that differed both in the number of individuals sampled and the number of loci sampled, they assessed how often the minimize-deep-coalescences method recovered the true species tree, that is, the tree on which the sequences were simulated. A sampling design of nine loci and three individuals per species proved to be quite successful in recovering the correct species tree. The results of the modified minimize-deep-coalescences method used here provide empirical support for their best sampling design and serve as a useful reference point for researchers designing sampling protocols with the purpose of inferring phylogenies of closely related taxa.
Conclusions:
Using DNA sequence data from a suite of nuclear loci, phylogenetic, and population genetic analyses supported the traditionally accepted evolutionary relationships between three North American passerines. In contrast to the findings of an earlier study that used data from the mitochondrial gene cytochrome b to infer the phylogeny of the genus (KLICKA et al. 2001), we found, using both phylogenetic and population genetic-based methods, P. cyanea and P. amoena to be more closely related to each other than either is to P. caerulea. These results underscore the importance of sampling multiple individuals and multiple loci when estimating multiple demographic parameters simultaneously under complicated evolutionary scenarios, e.g., when allowing for postdivergence gene flow. Unfortunately, IM and other coalescent-based methods available for inferring population genetic parameters are often limited in some ways (e.g., how many populations they can handle, the sequence evolution models employed, etc). The continued improvement of these methods will facilitate the combined analytical approach we advocate in studies with a greater number of focal taxa.
|
APPENDIX: SAMPLING LOCALITIES FOR INDIVIDUALS INCLUDED IN THE STUDY TESTING SPECIES DIVERGENCE HYPOTHESES IN PASSERINA BUNTINGS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: SAMPLING LOCALITIES...ACKNOWLEDGEMENTSLITERATURE CITED |
|---|
|
a Tissue obtained through loan from the University of Washington Burke Museum. All other tissue samples were obtained from the Genetic Resources Collection at the Louisiana State University Museum of Natural Science.
|
ACKNOWLEDGEMENTS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: SAMPLING LOCALITIES...ACKNOWLEDGEMENTSLITERATURE CITED |
|---|
|
FOOTNOTES |
|---|
|
LITERATURE CITED |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: SAMPLING LOCALITIES...ACKNOWLEDGEMENTSLITERATURE CITED |
|---|
AXELSSON, E., N. G. C. SMITH, H. SUNDSTROM, S. BERLIN and H. ELLEGREN, 2004 Male-biased mutation rate and divergence in autosomal, Z-linked and W-linked introns of chicken and turkey. Mol. Biol. Evol. 21: 1538–1547.
BAKER, M. C., and A. E. M. BAKER, 1990 Reproductive-behavior of female buntings - isolating mechanisms in a hybridizing pair of species. Evolution 44: 332–338.[CrossRef]
BAKER, M. C., and J. T. BOYLAN, 1999 Singing behavior, mating associations and reproductive success in a population of hybridizing Lazuli and Indigo Buntings. Condor 101: 493–504.
BLAKE, C. H., 1969 Notes of the Indigo Bunting. Bird Banding 40: 133–139.
BRAUN, M. J., and R. T. BRUMFIELD, 1998 Enigmatic phylogeny of skuas: an alternative hypothesis. Proc. R. Soc. Lond. Ser. B 265: 995–999.
CHEVIRON, Z. A., S. J. HACKETT and R. T. BRUMFIELD, 2006 Sequence variation in the coding region of the melanocortin-1 receptor gene (MC1R) is not associated with plumage variation in the Blue-crowned Manakin (Lepidothrix coronata). Proc. R. Soc. Lond. Ser. B 273: 1613–1618.[Medline]
DEGNAN, J. H., and N. A. ROSENBERG, 2006 Discordance of species trees with their most likely gene trees. PLoS Genet. 2: 0762–0768.
DEGNAN, J. H., and L. A. SALTER, 2005 Gene tree distributions under the coalescent process. Evolution 59: 24–37.[Medline]
DOYLE, J. J., 1992 Gene trees and species trees—molecular systematics as one-character taxonomy. Syst. Bot. 17: 144–163.[CrossRef]
EDWARDS, S. V., L. LIU and D. K. PEARL, 2007 High-resolution species trees without concatenation. Proc. Natl. Acad. Sci. USA 104: 5936–5941.
EMLEN, S. T., J. D. RISING and W. L. THOMPSON, 1975 A behavioral and morphological study of sympatry in the Indigo and Lazuli buntings of the Great Plains. Wilson Bull. 87: 145–179.
EXCOFFIER, L., G. LAVAL and S. SCHNEIDER, 2005 Arlequin ver 3.0: an integrated software package for population genetics data analysis. Evol. Bioinformatics Online 1: 47–50.
FELSENSTEIN, J., 2004 Inferring Phylogenies. Sinauer Associates, Sunderland, MA.
GILLESPIE, J. H., and C. H. LANGLEY, 1979 Are evolutionary rates really variable? J. Mol. Evol. 13: 27–34.[CrossRef][Medline]
GREENE, E., V. R. MUEHTER and W. DAVISON, 1996 Lazuli Bunting (Passerina amoena), pp. 1–24 in The Birds of North America, No. 232, edited by A. POOLE and F. GILL. The Academy of Natural Sciences, Washington, DC/The American Ornithologists' Union, Philadelphia.
HEIN, J., M. K. SCHIERUP and C. WIUF, 2005 Gene Genealogies, Variation and Evolution. Oxford University Press, New York.
HELLACK, J. J., and G. SCHNELL, 1977 Phenetic analysis of the subfamily Cardinalinae using external and skeletal characters. Wilson Bull. 89: 130–148.
HESLEWOOD, M. M., M. S. ELPHINSTONE, S. C. TIDEMANN and P. R. BAVERSTOCK, 1998 Myoglobin intron variation in the Gouldian Finch Erythrura gouldiae assessed by temperature gradient gel electrophoreses. Electrophoresis 19: 142–151.[CrossRef][Medline]
HEY, J., 2005 On the number of New World founders: a population genetic portrait of the peopling of the Americas. PLoS Biol. 3: 965–975.
HEY, J., and R. NIELSEN, 2004 Multilocus methods for estimating population sizes, migration rates, and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persimilis. Genetics 167: 747–760.
HILL, G. E., and K. J. MCGRAW, 2006 Bird Coloration, Vol. I. Harvard University Press, Cambridge, MA.
HOEKSTRA, H. E., 2006 Genetics, development and the evolution of adaptive pigmentation in vertebrates. Heredity 97: 222–234.[CrossRef][Medline]
HUELSENBECK, J. P., and F. RONQUIST, 2001 MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics 17: 754–755.
JENNINGS, W. B., and S. V. EDWARDS, 2005 Speciational history of Australian Grass Finches (Poephila) inferred from thirty gene trees. Evolution 59: 2033–2047.[Medline]
KIMURA, M., 1969 The number of heterozygous nucleotide sites maintained in a finite population due to steady flux of mutations. Genetics 61: 893–903.
KLICKA, J., A. J. FRY, R. M. ZINK and C. W. THOMPSON, 2001 A cytochrome-b perspective on Passerina bunting relationships. Auk 118: 611–623.
KROODSMA, R. L., 1975 Hybridization in buntings (Passerina) in North Dakota and eastern Montana. Auk 92: 66–80.
KUBATKO, L. S., and J. H. DEGNAN, 2007 Inconsistency of phylogenetic estimates from concatenated data under coalescence. Syst. Biol. 56: 17–24.[CrossRef][Medline]
MADDISON, W. P., 1997 Gene trees in species trees. Syst. Biol. 46: 523–536.[CrossRef]
MADDISON, W. P., and L. L. KNOWLES, 2006 Inferring phylogeny despite incomplete lineage sorting. Syst. Biol. 55: 21–30.[CrossRef][Medline]
MADDISON, W. P., and D. R. MADDISON, 2004 Mesquite: a modular system for evolutionary analysis, Version 1.01. http://mesquiteproject.org.
MAYR, E., and L. L. J. SHORT, 1970 Species taxa of North American birds, a contribution to comparative systematics. Pub. 9, Nuttall Ornithological Club, Cambridge, MA.
MCGUIRE, G., and F. WRIGHT, 1998 TOPAL: recombination detection in DNA and protein sequences. Bioinformatics 14: 219–220.
MCGUIRE, G., and F. WRIGHT, 2000 TOPAL 2.0: improved detection of mosaic sequences within multiple alignments. Bioinformatics 16: 130–134.
MUNDY, N. I., 2005 A window on the genetics of evolution: MC1R and plumage colouration in birds. Proc. R. Soc. Lond. Ser. B 272: 1633–1640.[Medline]
NATURESERVE, 2006 NatureServe Explorer: An Online Encyclopedia of Life, Version 6.1. http://www.natureserve.org/explorer.
NEI, M., and W. H. LI, 1979 Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc. Natl. Acad. Sci. USA 76: 5269–5273.
PAMILO, P., and M. NEI, 1988 Relationships between gene trees and species trees. Mol. Biol. Evol. 5: 568–583.[Abstract]
PAYNE, R. B., 1992 Indigo Bunting, pp. 1–24 in The Birds of North America, No. 4, edited by A. POOLE, P. STETTENHEIM and F. GILL. The Academy of Natural Sciences, Washington, DC/The American Ornithologists' Union, Philadelphia.
PAYNTER, R. A., 1970 Subfamily Cardinalidae, pp. 216–245 in Check-list of Birds of the World, edited by R. A. PAYNTER and R. W. STORER. Museum of Comparative Zoology/Heffernan Press, Worcester, MA.
PHILLIPS, A. R., J. T. MARSHALL and G. MONSON, 1964 The Birds of Arizona. University of Arizona Press, Tucson, AZ.
POSADA, D., and T. R. BUCKLEY, 2004 Model selection and model averaging in phylogenetics: advantages of the AIC and Bayesian approaches over likelihood ratio tests. Syst. Biol. 53: 793–808.[CrossRef][Medline]
POSADA, D., and K. A. CRANDALL, 1998 MODELTEST: testing the model of DNA substitution. Bioinformatics 14: 817–818.
PRIMMER, C. R., T. BORGE and G. P. SAETRE, 2002 Single-nucleotide polymorphism characterization in species with limited available sequence information: high nucleotide diversity revealed in the avian genome. Mol. Ecol. 11: 603–612.[CrossRef][Medline]
RIDGELY, R. S., T. F. ALLNUTT, T. BROOKS, D. K. MCNICOL, D. W. MEHLMAN et al., 2003 Digital distribution maps of the birds of the Western hemisphere, version 1.0. NatureServe, Arlington, VA.
RONQUIST, F., and J. P. HUELSENBECK, 2003 MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.
ROSENBERG, N. A., 2002 The probability of topological concordance of gene trees and species trees. Theor. Popul. Biol. 61: 225–247.[CrossRef][Medline]
ROSENBERG, N. A., 2003 The shapes of neutral gene genealogies in two species: probabilities of monophyly, paraphyly, and polyphyly in a coalescent model. Evolution 57: 1465–1477.[CrossRef][Medline]
SAITOU, N., and M. NEI, 1987 The neighbor-joining method—a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4: 406–425.[Abstract]
SHIMODAIRA, H., and M. HASEGAWA, 1999 Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol. Biol. Evol. 16: 1114–1116.
SIBLEY, C. G., and L. L. J. SHORT, 1959 Hybridization in the buntings (Passerina) of the Great Plains. Auk 76: 443–463.
SLADE, R. W., C. MORITZ, A. HEIDEMAN and P. T. HALE, 1993 Rapid assessment of single-copy nuclear DNA variation in diverse species. Mol. Ecol. 2: 359–373.[Medline]
SOKAL, R. R., and F. J. ROHLF, 1995 Biometry. W. H. Freeman, New York.
SWOFFORD, D. L., 1998 PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods), Version 4. Sinauer Associates, Sunderland, MA.
TAJIMA, F., 1983 Evolutonary relationships of DNA sequences in finite populations. Genetics 105: 437–460.
TAKAHATA, N., 1989 Gene genealogy in three related populations—consistency probability between gene and population trees. Genetics 122: 957–966.
TAKAHATA, N., 1995 A genetic perspective on the origin and history of humans. Annu. Rev. Ecol. Syst. 26: 343–372.[CrossRef]
TAKAHATA, N., and M. NEI, 1985 Gene genealogy and variance of interpopulational nucleotide differences. Genetics 110: 325–344.
TAKAHATA, N., and M. SLATKIN, 1990 Genealogy of neutral genes in two partially isolated populations. Theor. Popul. Biol. 38: 331–350.[CrossRef][Medline]
TAMPLIN, J. W., J. W. DEMASTES and J. V. REMSEN, 1993 Biochemical and morphometric relationships among some members of the Cardinalinae. Wilson Bull. 105: 93–113.
WAKELEY, J., 1996 Pairwise differences under a general model of population subdivision. J. Genet. 75: 81–89.
WAKELEY, J., and J. HEY, 1998 Testing speciation models with DNA sequence data, pp. 157–167 in Molecular Approaches to Ecology and Evolution, edited by R. DESALLE and R. SCHIERWATER. Birkhauser Verlag, Basel, Switzerland.
WON, Y. J., and J. HEY, 2005 Divergence population genetics of chimpanzees. Mol. Biol. Evol. 22: 297–307.
| ||||||||
