BACKGROUND
we previously used mitochondrial dna analyses to describe the history of population expansion and divergence of the bananaquit , the most abundant and widely distributed songbird in the west indies, which also occurs in the humid tropics of continental america  <cit> .
the geographic distribution and relationship of mtdna rflp haplotypes suggested that populations of this species were periodically invasive, and then quiescent.
here we increase our sample of populations and genes to dissect more finely the population history of bananaquits.
the bananaquit is a small nectivorous and frugivorous emberizine bird  that is an abundant resident throughout the west indies, except for cuba  <cit> .
on the continent, it is widely distributed from southern mexico through much of south america.
geographic variation in plumage coloration led to the recognition of  <dig> subspecies throughout the range of the bananaquit  <cit> , but mtdna rflp analysis failed to support many of these distinctions  <cit> .
the mtdna analyses also revealed levels of genetic divergence between some populations within the range of genetic differences observed between species  <cit> .
building on our earlier studies of the bananaquit  <cit> , we have now included representatives of more populations from the west indies and the continent, expanded our nucleotide sample of the mitochondrial genome of each individual to  <dig> base pairs, and have added at least  <dig> base pairs representing three nuclear gene regions.
our phylogenetic analysis is rooted by the addition of five individuals of three species of emberizid finches representing three genera.
our principal objectives were to determine the ancestral area of c. flaveola, to infer the number and directions of colonization events within the west indies and between the islands and the continent, and to assess the power of nuclear versus mitochondrial genes to inform analyses of colonization history and subsequent population dynamics.
our phylogenetic results locate the geographic ancestral area of the bananaquit in the region of the greater antilles and bahamas, and strongly support the hypothesis that continental bananaquits were derived from west indian ancestors.
in addition, comparisons of mitochondrial and nuclear gene genealogies exposed a complex and dynamic history of bananaquits in the greater and lesser antilles that would have remained masked without recourse to both maternal and biparental genealogical markers.
RESULTS
sequencing, nucleotide composition and saturation
we obtained clean sequences for all  <dig> c. flaveola distributed within the west indies and south-central america  and for six outgroup samples, with a few exceptions, for six mitochondrial genes and three nuclear genes .
all sequences were submitted to genbank .
cloning was performed on nine samples for bfib <dig>   <dig> samples for rag <dig>  and two samples for chdz.
cloning error rate was estimated to be about  <dig> % and cloning errors were eliminated from cloned alleles.
several properties of the data were examined to ascertain the utility of the genes for phylogenetic analysis .
none of the mitochondrial genes exhibited stop codons when translated into amino acids.
a test of homogeneity indicated that sequences had similar base frequencies across taxa .
as expected, mtdna sequences showed an anti-guanine bias, except for nd <dig>  which showed an anti-cytosine bias and is the only one of the six genes encoded by the l strand; all other protein-coding mitochondrial genes are encoded by the h strand.
the majority of variable sites in the mitochondrial datasets were third-position substitutions, as expected for protein-coding genes.
the mitochondrial sequence dataset included more variation than the nuclear dataset, with a mean percentage of variable sites ranging from  <dig>  to  <dig>  for mitochondrial genes compared to  <dig>   to  <dig>  for nuclear genes.
we calculated the relative substitution rate between mitochondrial and nuclear data by dividing the net average mitochondrial distance and the net average nuclear distance between the bahamas individuals and the rest of the c. flaveola phylogeny.
the mean substitution rate was approximately  <dig> times slower for nuclear genes than for mitochondrial genes.
plots of the transition/transversion ratio against uncorrected distances indicated that nucleotide substitution in the mitochondrial genome was not saturated .
all plots were linear whether plotted by position or by gene.
thus, all nucleotides were employed in the ensuing phylogenetic analyses of c. flaveola.
phylogenetic analyses: outgroups
the closest relatives to c. flaveola  <cit>  chosen here as outgroups to root our phylogenies  were divergent by at least  <dig> % from a mitochondrial perspective and at least 2% from a nuclear perspective.
phylogenetic analyses: congruence among tree topologies using different approaches and different genes
phylogenetic trees generated using the three approaches – mp, ml and bayesian – were congruent for both mitochondrial and nuclear datasets.
we also failed to reject homogeneity of phylogenetic trees based on each mitochondrial gene region independently  or each nuclear gene region independently .
therefore, for the subsequent phylogenetic analyses, we combined all six mitochondrial genes for the mitochondrial dataset and all three nuclear genes for the nuclear dataset.
concatenation of the different genes for nuclear data can be problematic as independent sorting of alleles in different parts of the genome and stochasticity of lineage sorting may lead to erroneous species trees  <cit> .
however, our results should be unaffected by this issue because the nuclear and mitochondrial trees are congruent, with the exception of the mtdna of hispaniolan birds .
the gtr+i+g model of evolution was selected using the aic criterion in mr modeltest for both the mitochondrial and the nuclear datasets.
phylogenetic analysis using mtdna
when describing the phylogeographic results, we use the following area name codes: la = lesser antilles; pr = puerto rico; ga = greater antilles; sca = south and central america; bh/qr = bahamas/quintana roo; dr = dominican republic ; ja/cy = jamaica/cayman islands).
a consensus phylogenetic tree based on a bayesian analysis combining all six mitochondrial genes identified four major mitochondrial clades representing bh and qr, ga , la , and sca .
relationships were strongly supported with marginal posterior probability values ≥  <dig>  for all  <dig> major nodes.
the continental  clade of bananaquits was sister to the pr and la clades, which were imbedded within the ga  clades.
this topology indicates derivation of the continental clade of bananaquits from within the ga. within the pr/la clade, sequences from pr had a sister relationship to the rest of the clade, and bananaquits on grenada and saint vincent form a distinct phylogenetic group sister to populations in the northern la.
within the sca clade, samples from central america  were grouped and distinct from samples from south america.
within south america, two groups were distinguished: birds from trinidad and venezuela and birds from peru and bolivia.
a sample from mexico  sequenced with nd <dig>  fell within the central america group  when added to our mitochondrial data.
phylogenetic analysis using nuclear dna
bayesian analysis combining all three nuclear genes  distinguished four clades:  bh,  ja and cy,  dr together with isolated haplotypes from montserrat, barbados, saint lucia and pr, and  pr, la and sca.
all three genes showed the affiliation of continental haplotypes with pr/la haplotypes and the distinctiveness of dr and ja/cy haplotypes.
nevertheless each gene shows a specific pattern.
for example, bfib <dig>  the most variable gene, unites the four dr and seven la/pr haplotypes, a group that is separated from the other sequences by five nucleotide changes.
chdz, the least variable gene, separates with only one nucleotide change the ja/cy samples on one side and dr samples on the other side from the rest of the phylogeny .
mitochondrial versus nuclear phylogenetic analyses
in comparison with the mitochondrial tree, the tree based on nuclear data resolved the more recent divergences poorly.
for example, the nuclear data grouped the continental samples with those from the la.
nonetheless, the nuclear phylogeny was consistent with the mitochondrial tree in the basal branching pattern for c. flaveola within the west indies.
the distinctiveness of dr birds in the nuclear tree was surprising given their relative proximity to ja/cy birds in the mitochondrial tree.
because the discrepancy between mitochondrial and nuclear data suggests a distinctive phylogenetic history for the dr population, we sequenced eight additional dr individuals with both mitochondrial  and nuclear  genes.
those additional sequences confirmed that dr birds grouped with cy and ja birds from a mitochondrial perspective and formed a separate clade from the other populations from a nuclear perspective .
estimates of population splitting times
the likelihood-ratio test based on the six mitochondrial genes failed to reject the null hypothesis of rate constancy .
therefore, to estimate divergence times, we used linear mtdna substitution rates of  <dig>  and 3% ma- <dig>  which span the range of most calibrations .
the average genetic distance between bh/qr birds and all other populations was  <dig>  ±  <dig> , placing the ingroup root between  <dig>  and  <dig>  mya.
the average divergence between the ga and la/pr/sca groups was  <dig>  ±  <dig> , indicating that the split between those two clades occurred at approximately the same time, i.e.,  <dig>  to  <dig>  mya.
the divergence between the la/pr clade and the sca clade was  <dig>  ±  <dig> , implying a split about  <dig>  to  <dig>  mya.
finally, pr diverged from the la approximately  <dig> – <dig>  mya.
discussion
in this paper, we have described phylogeographic relationships of bananaquit populations to determine the history of colonization within the caribbean basin.
additional analyses using coalescent models to infer population parameters require a different sampling regime  and will be reported separately .
we provide a well-resolved phylogeny of the bananaquit in the west indies, based on multiple genes.
aspects of the bananaquit phylogeny highlight points of general interest that are not commonly appreciated for birds.
first, we infer that populations sampled in the greater antilles and bahamas were the source of mainland colonization for c. flaveola, which is the reverse of the usual continent-to-island direction of colonization.
second, our study demonstrated the dynamic nature of the bananaquit distribution, resulting from several phases of overlapping expansions, especially through the lesser antilles.
third, a discrepancy between mitochondrial and nuclear haplotypes among bananaquits on hispaniola  suggested the recent replacement of the mitochondrial genome by introgression from jamaica.
mitochondrial versus nuclear data: relative coalescence time, sorting of ancestral polymorphism and inferences concerning phylogenetic history
dna sequences of mitochondrial genes provide a single estimate of a species phylogeny based on one linkage group  <cit> .
in our analysis of bananaquit phylogeny, we adopted a multi-gene approach based on three nuclear genes as well as dna sequences for six mitochondrial genes.
the more rapid coalescence of mtdna represents an advantage as it is more likely to resolve nodes with short internode distances, and the probability of incomplete lineage sorting is reduced compared to nuclear dna  <cit> .
however, an evolutionary reconstruction based on multiple independent genes increases confidence that the true species tree has been recovered, and improves phylogenetic resolution by averaging out potentially misleading effects of ancestral polymorphisms  <cit> .
furthermore, analytical problems resulting from the stochastic nature of lineage sorting and nucleotide substitution are reduced with an increased sample of genes and nucleotides.
our aim, here, was to explore the evolutionary history of the bananaquit populations within the west indies from the vantage points offered by the different genes independently and collectively, considering their different mutation and substitution rates, coalescence times, modes of inheritance and effective population sizes.
we obtained a well-resolved phylogeny of c. flaveola based on analysis of the six mitochondrial genes.
the substitution rate of bananaquit mtdna was high relative to the observed rate of lineage splitting, as evidenced by abundant synapomorphies and complete lineage sorting  of the major clades.
the data exhibited no saturation and approximated clock-like nucleotide substitution.
thus, we are confident that the tree presented in figure  <dig> represents a strongly supported hypothesis of matrilineal relationship.
in contrast, nuclear genes demonstrated considerably lower resolution regarding the phylogenetic relationships of bananaquit populations.
the substitution rate of the nuclear genes was, on average,  <dig> times slower than the mitochondrial genes, and thus closely related clades strongly supported by the mitochondrial data  were not distinguished by the nuclear data.
nonetheless, the allelic divergence at nuclear loci and the relative frequency of nuclear alleles across populations  clearly identifies the genetic distinctiveness of the bananaquit clades we have named "bahamas", "jamaica" and "puerto rico", and lends strong support to the hypothesis that these clades separated early in the diversification of coereba.
the slow nuclear sorting rate relative to population diversification resulted in many examples of incomplete lineage sorting compared to the mitochondrial perspective.
alternatively, the nuclear data could be interpreted as revealing abundant gene flow mediated by male bananaquits, biparental gene flow followed by localized episodes of selective sweeps acting on the mitochondrial genome, or drift associated with reduction in female population size.
we suspect incomplete lineage sorting because deep in the bananaquit phylogeny both the nuclear and mitochondrial data consistently identify the same clades, with the single exception that individuals on hispaniola  have mitochondrial dna closely related to ja/cy bananaquits, whereas the nuclear markers are allied with those on pr and the north central la.
phylogenetic history of the dominican republic population on hispaniola revealed by the discrepancy between nuclear and mitochondrial data
nuclear markers separate dr from ja/cy birds, whereas those populations are grouped in the same clade from a mitochondrial perspective.
this surprising result was confirmed by sequencing eight additional samples from this population for nuclear  and mitochondrial  genes .
recent introgression of ja/cy mitochondrial haplotypes into the dr population is the most likely explanation of conflicting positions of the dr samples in the nuclear and mitochondrial trees .
although ja/cy mtdna completely replaced the local dr mtdna, we found no evidence of introgression in the nuclear genome.
the challenge is to understand the complete introgression of the ja mitotype against a predominately dr nuclear background.
mitochondrial genomes are particularly susceptible to introgression, although this phenomenon is not clearly understood .
the ca.
four-fold smaller effective population size of the mtdna genome compared to the nuclear genome implies a higher probability for introgression by mtdna.
however, this outcome would be relatively unlikely if the number of ja immigrants was small relative to the population of resident dr bananaquits.
of course, if immigration from ja coincided with a population crash of hispaniolan bananaquits, the likelihood of stochastic fixation of the ja mtdna clade would have increased accordingly.
alternatively, mitochondrial introgression could have been favoured by selection, while the nuclear genomes of the two populations were incompatible.
this type of selection-driven introgression has been demonstrated in drosophila from mtdna microinjection studies  <cit>  and more recently in wild goats from a molecular phylogenetic study  <cit> .
the authors of the latter study suggested that proto-hemitragus mtdna could have invaded the ancestral population of capra because the former were better adapted to high altitudes.
however, we cannot imagine differences in the selective environments of basal bananaquit populations that might have led to the selective advantage of the ja mitochondrion, considering that the major islands of the greater antilles have a similar range of environments.
a third hypothesis is that sexual selection could account for the observation, much in the same manner as hypothesized for dendroica occidentalis warblers along the west coast and islands of british colombia  <cit> .
in this case, hybrid ja/dr males would have a mating disadvantage, slowing the introgression of nuclear genes relative to the mitochrondrial genomes passed on only by the females.
each backcross generation would dilute the contribution of ja nuclear alleles, while ja mitotypes might become fixed by chance.
however, unless the ja mitochondrion experienced a selective advantage, this scenario is likely only in the case of greatly reduced dr populations.
other biological processes that can create topological incongruence between gene trees include lineage sorting, paralogy, and lateral gene transfer  <cit> .
these alternatives can easily be discounted in the case of dr bananaquits.
first, because several independent nuclear gene sequences support the complete divergence of the nuclear genomes of ja/cy and dr birds, the combination of dr nuclear and ja/cy mitochondrial genomes is not likely the result of incomplete lineage sorting of an ancestral polymorphism in nuclear alleles.
although scenarios implying gene duplications and gene deletions could be inferred for one nuclear gene only, it is highly unlikely that the same scenario may have occurred independently and identically in the three unlinked nuclear genes .
indeed, the absence of substantial divergence between ja/cy and dr mitochondrial genomes indicates recent derivation.
second, phylogenetic trees using paralogous sequences may be misinterpreted if one assumes that the nuclear sequences are orthologous.
orthologous genes derive from the same locus whereas paralogous genes derive from different loci that originated by gene duplication  <cit> .
however, reconstructing an erroneous gene tree based on paralogy or gene duplication/deletion is highly improbable because it would imply one multiple-gene or genome-wide duplication event at the root of the tree, a panmictic dr/ja/cy population until recently , and several recent nuclear gene deletions.
finally, lateral gene transfer other than by hybridization/introgression can happen when an organism incorporates genetic material from another distantly related organism, and this can create topological incongruence between gene trees.
however, animals seem to be largely unaffected by this phenomenon  <cit>  and it has never been reported for birds.
the evolutionary history of bananaquits
the mitochondrial phylogeny pictured in figure  <dig> confirms the phylogeographic structure of bananaquits revealed by the earlier rflp analysis  <cit> .
regional mtdna haplotype clades displayed levels of sequence divergence typically in excess of 2%, with the largest distances exceeding 5% .
the nuclear dna data confirm the mtdna-based hypothesis of early divergence amongst bananaquit populations today revealed by three highly differentiated regional assemblages: bh, ja/cy, and pr/la/sca.
bananaquit lineages representing the bh, ja and pr were probably isolated soon after the initial geographic spread of the group in the ga. subsequent expansions from the bh and ja were geographically restricted.
bananaquits from the bh probably spread to the yucatán peninsula of mexico early in the history of that clade based on the mtdna sequence of a single bird collected from qr  and on the island itself at gibara , but breeding has not been established  <cit> .
to summarize, bananaquits from coastal yucatan deserve additional study, but data in hand suggest that the matriline allied to the bh is narrowly restricted in comparison to the predominant mtdna haplotype clade in middle america.
the lineage of bananaquits found today on ja/cy is also narrowly distributed from a joint nuclear and mtdna perspective.
given the predominance of the ja/cy matriline in the dr, however, it is clear that at least female bananaquits have moved between the islands relatively recently.
thus, migration or range expansion of the ja clade group has been recent in comparison to that of the bh mtdna haplotype clade.
nuclear and mitochondrial data indicate that most of the contemporary distribution of bananaquits results from the spread of birds whose descendents are now found in pr, la and sca.
the geographic expansion of this clade is relatively recent in comparison to the original diversification of bananaquits, and the first split of this expanding lineage  separated island from continental birds.
subsequent diversification on the mainland led to reciprocally monophyletic mtdna haplotype clades, one predominant in south america and the other in middle america.
in the islands, pr populations became isolated from la populations, and then bananaquits in the north central la became separated from the populations on grenada and st. vincent.
the pattern of evolutionary quiescence followed by geographic expansion is also repeated in the mtdna record of bananaquits in the la.
this history was clearly revealed in our earlier research  <cit>  based on a mtdna rflp analysis of a larger sample of birds per island, and demonstrated that bananaquits had recently expanded their distribution through the northern islands in the chain.
however, we did not anticipate the northward expansion of the north-central lesser antilles bananaquit clade all the way to the british virgin islands.
it is particularly noteworthy that bananaquits crossed the anegada gap, which serves as a significant biogeographic divide in the west indies, to colonize islands that would have been connected to pr, or nearly so, during pleistocene low sea level stands.
this seems improbable if the british virgin islands had been occupied by the puerto rican population of bananaquits, and suggests that a catastrophic event might have caused the extirpation of populations in the northern la and the british virgin islands.
the alternative scenario of north-central la bananaquits displacing pr like birds in the british virgin islands as a result of competitive superiority is lessened by the evident failure of this lineage to invade puerto rico.
our sampling in pr included coastal areas on the eastern end of the island close to the virgin islands.
at the other extreme of the la, grenada and saint vincent have maintained their evolutionary separation from islands to the north, but no significant separation from one another as measured genetically or morphologically.
both islands harbor melanistic forms of the bananaquit, as well as shared mtdna haplotypes, that are not found on islands to the north.
the evolutionary connection between grenada and st. vincent might reflect their separation by a shallow bank of islands  that would have connected the islands during pleistocene low sea level stands  <cit> .
this is not the case for other islands in the la, which are mostly  separated one from the other by deep-water channels.
phylogenetic support for taxonomic distinctions among subspecies
our molecular systematic appraisal of bananaquits permits us to assess phylogenetic support for taxonomic distinctions among the many named subspecies.
two subspecies groups were recognized by paynter  <cit> : bahamensis  and flaveola , but most populations in the west indies have received specific or subspecific epithets at some time in the past, bond  <cit>  recognizing  <dig> subspecies.
our data confirm the genetic distinctiveness of the bahamensis birds; however, the flaveola group encompasses mtdna haplotype clades that are nearly as divergent as flaveola and bahamensis groups.
for example the ja/cy/dr clade is approximately 5% diverged in mitochondrial sequence from the remainder of the populations recognized as the flaveola subspecies by paynter  <cit> .
additionally, subspecies such as portoricensis on puerto rico represent distinct evolutionary lineages.
however, more often than not closely related or identical mtdna haplotypes encompass several subspecies.
for instance, the northern la mtdna haplotype clade includes four subspecies , the two panama bananaquits comprise different subspecies , and birds from bolivia and peru belong to three different subspecies .
although finer genetic distinctions consistent with morphological differences between populations may exist, this study adds evidence that presently recognized bird subspecies often do not represent historically and phylogenetically equivalent evolutionary lineages .
revision of earlier phylogenetic and biogeographic history of the bananaquit
our analysis also allows us to review the phylogenetic and biogeographic history of the lineage.
bond  <cit>  interpreted the distribution of c. flaveola in the west indies as resulting from two invasions, one from south america spreading north through the la and west to ja, and the other from central america spreading north and east to the bh.
our initial mitochondrial rflp-based survey of bananaquits  <cit>  indicated a different, more complex history.
however lacking samples from bananaquit populations on dominican republic and the bahamas, as well as a suitable outgroup, we could not satisfactorily address bond's  <cit>  hypothesis of a continental origin for the species  <cit> .
results reported here confirm and extend the earlier molecular phylogenetic analyses.
the sister-group relationship between sca bananaquit populations rejects bond's hypothesis, which predicted that the two continental regions should be most closely related to the island populations founded from south and central america, respectively.
islands as the ancestral area of the bananaquit and "reverse colonisation" of the continent
because species diversity on islands decreases with distance from continental source areas  <cit> , biologists have assumed that colonization of islands or archipelagos is a one-way process, and that "reverse colonization" to the mainland against a diversity gradient rarely, if ever, occurs.
our study confirms previous indications that island species also can colonize continental areas.
although the depauperate biotas of remote islands and archipelagos have little chance of recolonizing continental landmasses, near islands, such as the west indies, have larger biotas with fewer impediments to dispersal, considerably improving the probability of back migration, or onward migration.
colonization from the islands to the mainland has similarly been reported within the west indies for several bird species .
island to mainland colonization has also been demonstrated in monarch flycatchers from the solomon islands to new guinea/australia  <cit>  .
in the case of the bananaquit, we have provided both mitochondrial and nuclear evidence that permits strong inference that birds from the islands have colonized the mainland.
the nesting of the south american to panama samples within the west indian clades and the deep paraphyly of bananaquits from the bh and ja clade with respect to all other bananaquits, including all mainland birds , represent the strongest lines of evidence that the ancestral node of the extant lineages is most parsimoniously placed in the greater antilles.
furthermore, ancestral area analysis  based on the bremer method  <cit>  clearly excludes the continent as the ancestral area of coereba.
finally, near outgroups for the bananaquit are principally emberizid finches from the west indies  <cit> .
we cannot place the origin  of the species itself, but this issue is not relevant to the present study.
to conclude, islands, including the west indies, might be significant sources of biodiversity for continents, and we emphasize the importance of considering reverse colonisation for interpreting biogeographic patterns.
diversification and range expansion of coereba flaveola in the west indies: alternative scenarios
alternative scenarios of diversification and range expansion, from the greater antilles/bahamas, can explain the phylogenetic patterns observed in both the mitochondrial and nuclear gene trees.
they differ mainly in whether the continent was invaded directly from the greater antilles or through the lesser antilles.
in the simplest scenario, bananaquits expanded from the ga to independently colonize the la  and the continent  approximately  <dig>  to  <dig> mya, through an eastern route.
considering the affinity of south american genotypes to those on pr rather than dr or ja, a western route can be excluded.
this scenario is consistent from a mitochondrial point of view.
it does not exclude the possibility that several phases of expansion might have occurred within the la .
however, direct colonization of south america from ga seems unlikely for several reasons.
first, bananaquits, which do not flock and are not strong fliers, would have difficulty making a long over-water flight.
second, virtually all colonization events in the west indies occur in stepping-stone fashion from one island to the next, rather than over long distances  <cit> .
third, remnants of the nuclear genome of the south american birds in pr and la  favour dispersal through the islands rather than directly from the ga.
we envision colonization from ga to sca through la as occurring in the following way .
first, bananaquits spread from the dr to pr and through the la approximately  <dig>  to  <dig> mya.
evidence of this expansion is seen in the clade of haplotypes in the nuclear phylogeny recovered from dr, pr, and la.
this expansion did not reach the continent of south america.
a second expansion would have occurred from pr through the la to sca about  <dig>  to  <dig>  mya, as seen in the phylogenetic link between sca, pr, and la populations, exclusive of other ga populations, particularly dr .
this is also corroborated by the fact that the combined nuclear genes  show the affiliation of continental haplotypes with la and pr haplotypes, to the exclusion of dr haplotypes.
it is also possible that the dr/pr/la and pr/la/sca clades represent the same expansion, followed by lineage sorting.
a third, more recent expansion from pr through the la, but not to the continent, is indicated by the shared mitochondrial haplotypes of those two populations about  <dig>  to  <dig>  mya.
after the first and second expansions, pr birds would have had sufficient time to diverge from dr and continental populations from a mitochondrial perspective, but incompletely from a nuclear perspective.
this would explain why the la clade does not include dr mitochondrial sequences but does retain remnants of the dr nuclear genome.
CONCLUSIONS
analysis of genetic diversity in the bananaquit shows that the different rates of nuclear and mitochondrial genome processes, and incongruities between nuclear and mitochondrial phylogenies, can be used to dissect the often complex phylogeographic history of populations.
we emphasize the importance of including multiple genes with distinct evolutionary histories.
our study is the first, to our knowledge, to demonstrate such biogeographic complexity for a single species.
first, we demonstrated that the ancestral node of the extant bananaquit lineages is most parsimoniously placed in the greater antilles, indicating that islands, including the west indies, might be significant sources of biodiversity for continents.
we suggest that the "reverse colonisation" phenomenon deserves more systematic attention in phylogenetic analyses, and assessing its significance is of prime importance in studies of biogeography and ecology.
second, several phases of colonization and mixing among islands occurred before the bananaquit populations became isolated long enough to develop reciprocal monophyly in mtdna.
third, the mixed nuclear and mitochondrial genomes in the dominican republic and the remnants of old dominican republic nuclear lineages in the lesser antilles suggest that expanding bananaquit populations mix with resident populations during expansion phases.
thus, genetic compatibility is maintained over periods exceeding three million years and perhaps much more, despite morphological and song differentiation.
this does not exclude the possibility that expansions can follow extinctions of local populations or that colonists can competitively exclude local populations, as suggested by non-overlapping distributions of recently expanded mtdna lineages in several species of lesser antillean birds  <cit> .

 we previously used mitochondrial dna analyses to describe the history of population expansion and divergence of the bananaquit, the most abundant and widely distributed songbird in the west indies, which also occurs in the humid tropics of continental america. the geographic distribution and relationship of mtdna rflp haplotypes suggested that populations of this species were periodically invasive, and then quiescent. here we increase our sample of populations and genes to dissect more finely the population history of bananaquits. the bananaquit is a small nectivorous and frugivorous emberizine bird that is an abundant resident throughout the west indies, except for cuba. on the continent, it is widely distributed from southern mexico through much of south america. geographic variation in plumage coloration led to the recognition of subspecies throughout the range of the bananaquit, but mtdna rflp analysis failed to support many of these distinctions. the mtdna analyses also revealed levels of genetic divergence between some populations within the range of genetic differences observed between species. building on our earlier studies of the bananaquit, we have now included representatives of more populations from the west indies and the continent, expanded our nucleotide sample of the mitochondrial genome of each individual to base pairs, and have added at least base pairs representing three nuclear gene regions. our phylogenetic analysis is rooted by the addition of five individuals of three species of emberizid finches representing three genera. our principal objectives were to determine the ancestral area of c. flaveola, to infer the number and directions of colonization events within the west indies and between the islands and the continent, and to assess the power of nuclear versus mitochondrial genes to inform analyses of colonization history and subsequent population dynamics. our phylogenetic results locate the geographic ancestral area of the bananaquit in the region of the greater antilles and bahamas, and strongly support the hypothesis that continental bananaquits were derived from west indian ancestors. in addition, comparisons of mitochondrial and nuclear gene genealogies exposed a complex and dynamic history of bananaquits in the greater and lesser antilles that would have remained masked without recourse to both maternal and biparental genealogical markers.  sequencing, nucleotide composition and saturation we obtained clean sequences for all c. flaveola distributed within the west indies and south-central america and for six outgroup samples, with a few exceptions, for six mitochondrial genes and three nuclear genes. all sequences were submitted to genbank. cloning was performed on nine samples for bfib samples for rag and two samples for chdz. cloning error rate was estimated to be about % and cloning errors were eliminated from cloned alleles. several properties of the data were examined to ascertain the utility of the genes for phylogenetic analysis. none of the mitochondrial genes exhibited stop codons when translated into amino acids. a test of homogeneity indicated that sequences had similar base frequencies across taxa. as expected, mtdna sequences showed an anti-guanine bias, except for nd which showed an anti-cytosine bias and is the only one of the six genes encoded by the l strand; all other protein-coding mitochondrial genes are encoded by the h strand. the majority of variable sites in the mitochondrial datasets were third-position substitutions, as expected for protein-coding genes. the mitochondrial sequence dataset included more variation than the nuclear dataset, with a mean percentage of variable sites ranging from to for mitochondrial genes compared to to for nuclear genes. we calculated the relative substitution rate between mitochondrial and nuclear data by dividing the net average mitochondrial distance and the net average nuclear distance between the bahamas individuals and the rest of the c. flaveola phylogeny. the mean substitution rate was approximately times slower for nuclear genes than for mitochondrial genes. plots of the transition/transversion ratio against uncorrected distances indicated that nucleotide substitution in the mitochondrial genome was not saturated. all plots were linear whether plotted by position or by gene. thus, all nucleotides were employed in the ensuing phylogenetic analyses of c. flaveola. phylogenetic analyses: outgroups the closest relatives to c. flaveola chosen here as outgroups to root our phylogenies were divergent by at least % from a mitochondrial perspective and at least 2% from a nuclear perspective. phylogenetic analyses: congruence among tree topologies using different approaches and different genes phylogenetic trees generated using the three approaches – mp, ml and bayesian – were congruent for both mitochondrial and nuclear datasets. we also failed to reject homogeneity of phylogenetic trees based on each mitochondrial gene region independently or each nuclear gene region independently. therefore, for the subsequent phylogenetic analyses, we combined all six mitochondrial genes for the mitochondrial dataset and all three nuclear genes for the nuclear dataset. concatenation of the different genes for nuclear data can be problematic as independent sorting of alleles in different parts of the genome and stochasticity of lineage sorting may lead to erroneous species trees. however, our results should be unaffected by this issue because the nuclear and mitochondrial trees are congruent, with the exception of the mtdna of hispaniolan birds. the gtr+i+g model of evolution was selected using the aic criterion in mr modeltest for both the mitochondrial and the nuclear datasets. phylogenetic analysis using mtdna when describing the phylogeographic results, we use the following area name codes: la = lesser antilles; pr = puerto rico; ga = greater antilles; sca = south and central america; bh/qr = bahamas/quintana roo; dr = dominican republic; ja/cy = jamaica/cayman islands). a consensus phylogenetic tree based on a bayesian analysis combining all six mitochondrial genes identified four major mitochondrial clades representing bh and qr, ga, la, and sca. relationships were strongly supported with marginal posterior probability values ≥ for all major nodes. the continental clade of bananaquits was sister to the pr and la clades, which were imbedded within the ga clades. this topology indicates derivation of the continental clade of bananaquits from within the ga. within the pr/la clade, sequences from pr had a sister relationship to the rest of the clade, and bananaquits on grenada and saint vincent form a distinct phylogenetic group sister to populations in the northern la. within the sca clade, samples from central america were grouped and distinct from samples from south america. within south america, two groups were distinguished: birds from trinidad and venezuela and birds from peru and bolivia. a sample from mexico sequenced with nd fell within the central america group when added to our mitochondrial data. phylogenetic analysis using nuclear dna bayesian analysis combining all three nuclear genes distinguished four clades: bh, ja and cy, dr together with isolated haplotypes from montserrat, barbados, saint lucia and pr, and pr, la and sca. all three genes showed the affiliation of continental haplotypes with pr/la haplotypes and the distinctiveness of dr and ja/cy haplotypes. nevertheless each gene shows a specific pattern. for example, bfib the most variable gene, unites the four dr and seven la/pr haplotypes, a group that is separated from the other sequences by five nucleotide changes. chdz, the least variable gene, separates with only one nucleotide change the ja/cy samples on one side and dr samples on the other side from the rest of the phylogeny. mitochondrial versus nuclear phylogenetic analyses in comparison with the mitochondrial tree, the tree based on nuclear data resolved the more recent divergences poorly. for example, the nuclear data grouped the continental samples with those from the la. nonetheless, the nuclear phylogeny was consistent with the mitochondrial tree in the basal branching pattern for c. flaveola within the west indies. the distinctiveness of dr birds in the nuclear tree was surprising given their relative proximity to ja/cy birds in the mitochondrial tree. because the discrepancy between mitochondrial and nuclear data suggests a distinctive phylogenetic history for the dr population, we sequenced eight additional dr individuals with both mitochondrial and nuclear genes. those additional sequences confirmed that dr birds grouped with cy and ja birds from a mitochondrial perspective and formed a separate clade from the other populations from a nuclear perspective. estimates of population splitting times the likelihood-ratio test based on the six mitochondrial genes failed to reject the null hypothesis of rate constancy. therefore, to estimate divergence times, we used linear mtdna substitution rates of and 3% ma- which span the range of most calibrations. the average genetic distance between bh/qr birds and all other populations was ±, placing the ingroup root between and mya. the average divergence between the ga and la/pr/sca groups was ±, indicating that the split between those two clades occurred at approximately the same time, i.e., to mya. the divergence between the la/pr clade and the sca clade was ±, implying a split about to mya. finally, pr diverged from the la approximately – mya. discussion in this paper, we have described phylogeographic relationships of bananaquit populations to determine the history of colonization within the caribbean basin. additional analyses using coalescent models to infer population parameters require a different sampling regime and will be reported separately. we provide a well-resolved phylogeny of the bananaquit in the west indies, based on multiple genes. aspects of the bananaquit phylogeny highlight points of general interest that are not commonly appreciated for birds. first, we infer that populations sampled in the greater antilles and bahamas were the source of mainland colonization for c. flaveola, which is the reverse of the usual continent-to-island direction of colonization. second, our study demonstrated the dynamic nature of the bananaquit distribution, resulting from several phases of overlapping expansions, especially through the lesser antilles. third, a discrepancy between mitochondrial and nuclear haplotypes among bananaquits on hispaniola suggested the recent replacement of the mitochondrial genome by introgression from jamaica. mitochondrial versus nuclear data: relative coalescence time, sorting of ancestral polymorphism and inferences concerning phylogenetic history dna sequences of mitochondrial genes provide a single estimate of a species phylogeny based on one linkage group. in our analysis of bananaquit phylogeny, we adopted a multi-gene approach based on three nuclear genes as well as dna sequences for six mitochondrial genes. the more rapid coalescence of mtdna represents an advantage as it is more likely to resolve nodes with short internode distances, and the probability of incomplete lineage sorting is reduced compared to nuclear dna. however, an evolutionary reconstruction based on multiple independent genes increases confidence that the true species tree has been recovered, and improves phylogenetic resolution by averaging out potentially misleading effects of ancestral polymorphisms. furthermore, analytical problems resulting from the stochastic nature of lineage sorting and nucleotide substitution are reduced with an increased sample of genes and nucleotides. our aim, here, was to explore the evolutionary history of the bananaquit populations within the west indies from the vantage points offered by the different genes independently and collectively, considering their different mutation and substitution rates, coalescence times, modes of inheritance and effective population sizes. we obtained a well-resolved phylogeny of c. flaveola based on analysis of the six mitochondrial genes. the substitution rate of bananaquit mtdna was high relative to the observed rate of lineage splitting, as evidenced by abundant synapomorphies and complete lineage sorting of the major clades. the data exhibited no saturation and approximated clock-like nucleotide substitution. thus, we are confident that the tree presented in figure represents a strongly supported hypothesis of matrilineal relationship. in contrast, nuclear genes demonstrated considerably lower resolution regarding the phylogenetic relationships of bananaquit populations. the substitution rate of the nuclear genes was, on average, times slower than the mitochondrial genes, and thus closely related clades strongly supported by the mitochondrial data were not distinguished by the nuclear data. nonetheless, the allelic divergence at nuclear loci and the relative frequency of nuclear alleles across populations clearly identifies the genetic distinctiveness of the bananaquit clades we have named "bahamas", "jamaica" and "puerto rico", and lends strong support to the hypothesis that these clades separated early in the diversification of coereba. the slow nuclear sorting rate relative to population diversification resulted in many examples of incomplete lineage sorting compared to the mitochondrial perspective. alternatively, the nuclear data could be interpreted as revealing abundant gene flow mediated by male bananaquits, biparental gene flow followed by localized episodes of selective sweeps acting on the mitochondrial genome, or drift associated with reduction in female population size. we suspect incomplete lineage sorting because deep in the bananaquit phylogeny both the nuclear and mitochondrial data consistently identify the same clades, with the single exception that individuals on hispaniola have mitochondrial dna closely related to ja/cy bananaquits, whereas the nuclear markers are allied with those on pr and the north central la. phylogenetic history of the dominican republic population on hispaniola revealed by the discrepancy between nuclear and mitochondrial data nuclear markers separate dr from ja/cy birds, whereas those populations are grouped in the same clade from a mitochondrial perspective. this surprising result was confirmed by sequencing eight additional samples from this population for nuclear and mitochondrial genes. recent introgression of ja/cy mitochondrial haplotypes into the dr population is the most likely explanation of conflicting positions of the dr samples in the nuclear and mitochondrial trees. although ja/cy mtdna completely replaced the local dr mtdna, we found no evidence of introgression in the nuclear genome. the challenge is to understand the complete introgression of the ja mitotype against a predominately dr nuclear background. mitochondrial genomes are particularly susceptible to introgression, although this phenomenon is not clearly understood. the ca. four-fold smaller effective population size of the mtdna genome compared to the nuclear genome implies a higher probability for introgression by mtdna. however, this outcome would be relatively unlikely if the number of ja immigrants was small relative to the population of resident dr bananaquits. of course, if immigration from ja coincided with a population crash of hispaniolan bananaquits, the likelihood of stochastic fixation of the ja mtdna clade would have increased accordingly. alternatively, mitochondrial introgression could have been favoured by selection, while the nuclear genomes of the two populations were incompatible. this type of selection-driven introgression has been demonstrated in drosophila from mtdna microinjection studies and more recently in wild goats from a molecular phylogenetic study. the authors of the latter study suggested that proto-hemitragus mtdna could have invaded the ancestral population of capra because the former were better adapted to high altitudes. however, we cannot imagine differences in the selective environments of basal bananaquit populations that might have led to the selective advantage of the ja mitochondrion, considering that the major islands of the greater antilles have a similar range of environments. a third hypothesis is that sexual selection could account for the observation, much in the same manner as hypothesized for dendroica occidentalis warblers along the west coast and islands of british colombia. in this case, hybrid ja/dr males would have a mating disadvantage, slowing the introgression of nuclear genes relative to the mitochrondrial genomes passed on only by the females. each backcross generation would dilute the contribution of ja nuclear alleles, while ja mitotypes might become fixed by chance. however, unless the ja mitochondrion experienced a selective advantage, this scenario is likely only in the case of greatly reduced dr populations. other biological processes that can create topological incongruence between gene trees include lineage sorting, paralogy, and lateral gene transfer. these alternatives can easily be discounted in the case of dr bananaquits. first, because several independent nuclear gene sequences support the complete divergence of the nuclear genomes of ja/cy and dr birds, the combination of dr nuclear and ja/cy mitochondrial genomes is not likely the result of incomplete lineage sorting of an ancestral polymorphism in nuclear alleles. although scenarios implying gene duplications and gene deletions could be inferred for one nuclear gene only, it is highly unlikely that the same scenario may have occurred independently and identically in the three unlinked nuclear genes. indeed, the absence of substantial divergence between ja/cy and dr mitochondrial genomes indicates recent derivation. second, phylogenetic trees using paralogous sequences may be misinterpreted if one assumes that the nuclear sequences are orthologous. orthologous genes derive from the same locus whereas paralogous genes derive from different loci that originated by gene duplication. however, reconstructing an erroneous gene tree based on paralogy or gene duplication/deletion is highly improbable because it would imply one multiple-gene or genome-wide duplication event at the root of the tree, a panmictic dr/ja/cy population until recently, and several recent nuclear gene deletions. finally, lateral gene transfer other than by hybridization/introgression can happen when an organism incorporates genetic material from another distantly related organism, and this can create topological incongruence between gene trees. however, animals seem to be largely unaffected by this phenomenon and it has never been reported for birds. the evolutionary history of bananaquits the mitochondrial phylogeny pictured in figure confirms the phylogeographic structure of bananaquits revealed by the earlier rflp analysis. regional mtdna haplotype clades displayed levels of sequence divergence typically in excess of 2%, with the largest distances exceeding 5%. the nuclear dna data confirm the mtdna-based hypothesis of early divergence amongst bananaquit populations today revealed by three highly differentiated regional assemblages: bh, ja/cy, and pr/la/sca. bananaquit lineages representing the bh, ja and pr were probably isolated soon after the initial geographic spread of the group in the ga. subsequent expansions from the bh and ja were geographically restricted. bananaquits from the bh probably spread to the yucatán peninsula of mexico early in the history of that clade based on the mtdna sequence of a single bird collected from qr and on the island itself at gibara, but breeding has not been established. to summarize, bananaquits from coastal yucatan deserve additional study, but data in hand suggest that the matriline allied to the bh is narrowly restricted in comparison to the predominant mtdna haplotype clade in middle america. the lineage of bananaquits found today on ja/cy is also narrowly distributed from a joint nuclear and mtdna perspective. given the predominance of the ja/cy matriline in the dr, however, it is clear that at least female bananaquits have moved between the islands relatively recently. thus, migration or range expansion of the ja clade group has been recent in comparison to that of the bh mtdna haplotype clade. nuclear and mitochondrial data indicate that most of the contemporary distribution of bananaquits results from the spread of birds whose descendents are now found in pr, la and sca. the geographic expansion of this clade is relatively recent in comparison to the original diversification of bananaquits, and the first split of this expanding lineage separated island from continental birds. subsequent diversification on the mainland led to reciprocally monophyletic mtdna haplotype clades, one predominant in south america and the other in middle america. in the islands, pr populations became isolated from la populations, and then bananaquits in the north central la became separated from the populations on grenada and st. vincent. the pattern of evolutionary quiescence followed by geographic expansion is also repeated in the mtdna record of bananaquits in the la. this history was clearly revealed in our earlier research based on a mtdna rflp analysis of a larger sample of birds per island, and demonstrated that bananaquits had recently expanded their distribution through the northern islands in the chain. however, we did not anticipate the northward expansion of the north-central lesser antilles bananaquit clade all the way to the british virgin islands. it is particularly noteworthy that bananaquits crossed the anegada gap, which serves as a significant biogeographic divide in the west indies, to colonize islands that would have been connected to pr, or nearly so, during pleistocene low sea level stands. this seems improbable if the british virgin islands had been occupied by the puerto rican population of bananaquits, and suggests that a catastrophic event might have caused the extirpation of populations in the northern la and the british virgin islands. the alternative scenario of north-central la bananaquits displacing pr like birds in the british virgin islands as a result of competitive superiority is lessened by the evident failure of this lineage to invade puerto rico. our sampling in pr included coastal areas on the eastern end of the island close to the virgin islands. at the other extreme of the la, grenada and saint vincent have maintained their evolutionary separation from islands to the north, but no significant separation from one another as measured genetically or morphologically. both islands harbor melanistic forms of the bananaquit, as well as shared mtdna haplotypes, that are not found on islands to the north. the evolutionary connection between grenada and st. vincent might reflect their separation by a shallow bank of islands that would have connected the islands during pleistocene low sea level stands. this is not the case for other islands in the la, which are mostly separated one from the other by deep-water channels. phylogenetic support for taxonomic distinctions among subspecies our molecular systematic appraisal of bananaquits permits us to assess phylogenetic support for taxonomic distinctions among the many named subspecies. two subspecies groups were recognized by paynter: bahamensis and flaveola, but most populations in the west indies have received specific or subspecific epithets at some time in the past, bond recognizing subspecies. our data confirm the genetic distinctiveness of the bahamensis birds; however, the flaveola group encompasses mtdna haplotype clades that are nearly as divergent as flaveola and bahamensis groups. for example the ja/cy/dr clade is approximately 5% diverged in mitochondrial sequence from the remainder of the populations recognized as the flaveola subspecies by paynter. additionally, subspecies such as portoricensis on puerto rico represent distinct evolutionary lineages. however, more often than not closely related or identical mtdna haplotypes encompass several subspecies. for instance, the northern la mtdna haplotype clade includes four subspecies, the two panama bananaquits comprise different subspecies, and birds from bolivia and peru belong to three different subspecies. although finer genetic distinctions consistent with morphological differences between populations may exist, this study adds evidence that presently recognized bird subspecies often do not represent historically and phylogenetically equivalent evolutionary lineages. revision of earlier phylogenetic and biogeographic history of the bananaquit our analysis also allows us to review the phylogenetic and biogeographic history of the lineage. bond interpreted the distribution of c. flaveola in the west indies as resulting from two invasions, one from south america spreading north through the la and west to ja, and the other from central america spreading north and east to the bh. our initial mitochondrial rflp-based survey of bananaquits indicated a different, more complex history. however lacking samples from bananaquit populations on dominican republic and the bahamas, as well as a suitable outgroup, we could not satisfactorily address bond's hypothesis of a continental origin for the species. results reported here confirm and extend the earlier molecular phylogenetic analyses. the sister-group relationship between sca bananaquit populations rejects bond's hypothesis, which predicted that the two continental regions should be most closely related to the island populations founded from south and central america, respectively. islands as the ancestral area of the bananaquit and "reverse colonisation" of the continent because species diversity on islands decreases with distance from continental source areas, biologists have assumed that colonization of islands or archipelagos is a one-way process, and that "reverse colonization" to the mainland against a diversity gradient rarely, if ever, occurs. our study confirms previous indications that island species also can colonize continental areas. although the depauperate biotas of remote islands and archipelagos have little chance of recolonizing continental landmasses, near islands, such as the west indies, have larger biotas with fewer impediments to dispersal, considerably improving the probability of back migration, or onward migration. colonization from the islands to the mainland has similarly been reported within the west indies for several bird species. island to mainland colonization has also been demonstrated in monarch flycatchers from the solomon islands to new guinea/australia. in the case of the bananaquit, we have provided both mitochondrial and nuclear evidence that permits strong inference that birds from the islands have colonized the mainland. the nesting of the south american to panama samples within the west indian clades and the deep paraphyly of bananaquits from the bh and ja clade with respect to all other bananaquits, including all mainland birds, represent the strongest lines of evidence that the ancestral node of the extant lineages is most parsimoniously placed in the greater antilles. furthermore, ancestral area analysis based on the bremer method clearly excludes the continent as the ancestral area of coereba. finally, near outgroups for the bananaquit are principally emberizid finches from the west indies. we cannot place the origin of the species itself, but this issue is not relevant to the present study. to conclude, islands, including the west indies, might be significant sources of biodiversity for continents, and we emphasize the importance of considering reverse colonisation for interpreting biogeographic patterns. diversification and range expansion of coereba flaveola in the west indies: alternative scenarios alternative scenarios of diversification and range expansion, from the greater antilles/bahamas, can explain the phylogenetic patterns observed in both the mitochondrial and nuclear gene trees. they differ mainly in whether the continent was invaded directly from the greater antilles or through the lesser antilles. in the simplest scenario, bananaquits expanded from the ga to independently colonize the la and the continent approximately to mya, through an eastern route. considering the affinity of south american genotypes to those on pr rather than dr or ja, a western route can be excluded. this scenario is consistent from a mitochondrial point of view. it does not exclude the possibility that several phases of expansion might have occurred within the la. however, direct colonization of south america from ga seems unlikely for several reasons. first, bananaquits, which do not flock and are not strong fliers, would have difficulty making a long over-water flight. second, virtually all colonization events in the west indies occur in stepping-stone fashion from one island to the next, rather than over long distances. third, remnants of the nuclear genome of the south american birds in pr and la favour dispersal through the islands rather than directly from the ga. we envision colonization from ga to sca through la as occurring in the following way. first, bananaquits spread from the dr to pr and through the la approximately to mya. evidence of this expansion is seen in the clade of haplotypes in the nuclear phylogeny recovered from dr, pr, and la. this expansion did not reach the continent of south america. a second expansion would have occurred from pr through the la to sca about to mya, as seen in the phylogenetic link between sca, pr, and la populations, exclusive of other ga populations, particularly dr. this is also corroborated by the fact that the combined nuclear genes show the affiliation of continental haplotypes with la and pr haplotypes, to the exclusion of dr haplotypes. it is also possible that the dr/pr/la and pr/la/sca clades represent the same expansion, followed by lineage sorting. a third, more recent expansion from pr through the la, but not to the continent, is indicated by the shared mitochondrial haplotypes of those two populations about to mya. after the first and second expansions, pr birds would have had sufficient time to diverge from dr and continental populations from a mitochondrial perspective, but incompletely from a nuclear perspective. this would explain why the la clade does not include dr mitochondrial sequences but does retain remnants of the dr nuclear genome.  analysis of genetic diversity in the bananaquit shows that the different rates of nuclear and mitochondrial genome processes, and incongruities between nuclear and mitochondrial phylogenies, can be used to dissect the often complex phylogeographic history of populations. we emphasize the importance of including multiple genes with distinct evolutionary histories. our study is the first, to our knowledge, to demonstrate such biogeographic complexity for a single species. first, we demonstrated that the ancestral node of the extant bananaquit lineages is most parsimoniously placed in the greater antilles, indicating that islands, including the west indies, might be significant sources of biodiversity for continents. we suggest that the "reverse colonisation" phenomenon deserves more systematic attention in phylogenetic analyses, and assessing its significance is of prime importance in studies of biogeography and ecology. second, several phases of colonization and mixing among islands occurred before the bananaquit populations became isolated long enough to develop reciprocal monophyly in mtdna. third, the mixed nuclear and mitochondrial genomes in the dominican republic and the remnants of old dominican republic nuclear lineages in the lesser antilles suggest that expanding bananaquit populations mix with resident populations during expansion phases. thus, genetic compatibility is maintained over periods exceeding three million years and perhaps much more, despite morphological and song differentiation. this does not exclude the possibility that expansions can follow extinctions of local populations or that colonists can competitively exclude local populations, as suggested by non-overlapping distributions of recently expanded mtdna lineages in several species of lesser antillean birds. 
our results suggest that multi-gene phylogenies will permit improved comparative analysis of the evolutionary histories of different lineages in the same geographical setting, which provide replicated "natural experiments" for testing evolutionary hypotheses. phylogenetic concordance between mitochondrial and nuclear genes at the base of the bananaquit phylogeny supported a west indian origin for continental populations. although nuclear genes provided insufficient information at the tips of the tree to further evaluate relationships of closely allied but strongly supported mitochondrial dna clades, the discrepancy between mitochondrial and nuclear data in the population of dominican republic suggested that the mitochondrial genome was recently acquired by introgression from jamaica. this study represents one of the most complete phylogeographic analyses of its kind and reveals three patterns that are not commonly appreciated in birds: island to mainland colonization, multiple expansion phases, and mitochondrial genome replacement. we used multi-gene analyses to investigate the evolutionary history of this species throughout its distribution in the west indies and in south and middle america. we infer the ancestral area of the present-day bananaquit populations, report on the species' phylogenetic, biogeographic and evolutionary history, and propose scenarios for its diversification and range expansion. multi-gene analysis showing genetic remnants of successive colonization events in the lesser antilles reinforced earlier research demonstrating that bananaquits alternate periods of invasiveness and colonization with biogeographic quiescence. we sequenced six mitochondrial genes and three nuclear genes for forty-four bananaquits and three outgroup species. the bananaquit is a small nectivorous and frugivorous emberizine bird that is an abundant resident throughout the caribbean region.
 the bananaquit is a small nectivorous and frugivorous emberizine bird that is an abundant resident throughout the caribbean region. we used multi-gene analyses to investigate the evolutionary history of this species throughout its distribution in the west indies and in south and middle america. we sequenced six mitochondrial genes and three nuclear genes for forty-four bananaquits and three outgroup species. we infer the ancestral area of the present-day bananaquit populations, report on the species' phylogenetic, biogeographic and evolutionary history, and propose scenarios for its diversification and range expansion.  phylogenetic concordance between mitochondrial and nuclear genes at the base of the bananaquit phylogeny supported a west indian origin for continental populations. multi-gene analysis showing genetic remnants of successive colonization events in the lesser antilles reinforced earlier research demonstrating that bananaquits alternate periods of invasiveness and colonization with biogeographic quiescence. although nuclear genes provided insufficient information at the tips of the tree to further evaluate relationships of closely allied but strongly supported mitochondrial dna clades, the discrepancy between mitochondrial and nuclear data in the population of dominican republic suggested that the mitochondrial genome was recently acquired by introgression from jamaica.  this study represents one of the most complete phylogeographic analyses of its kind and reveals three patterns that are not commonly appreciated in birds: island to mainland colonization, multiple expansion phases, and mitochondrial genome replacement. the detail revealed by this analysis will guide evolutionary analyses of populations in archipelagos such as the west indies, which include islands varying in size, age, and geological history. our results suggest that multi-gene phylogenies will permit improved comparative analysis of the evolutionary histories of different lineages in the same geographical setting, which provide replicated "natural experiments" for testing evolutionary hypotheses. 