# Tracing the colonization history of the Indian Ocean scops-owls (Strigiformes: Otus) with further insight into the spatio-temporal origin of the Malagasy avifauna

Jérôme Fuchs, Jean-Marc Pons, Steven M. Goodman, Vincent Bretagnolle, Martin Melo, Rauri Ck Bowie, David Currie, Roger Safford, Munir Z. Virani, Simon Thomsett, et al.

► To cite this version:

Jérôme Fuchs, Jean-Marc Pons, Steven M. Goodman, Vincent Bretagnolle, Martin Melo, et al. Tracing the colonization history of the Indian Ocean scops-owls (Strigiformes: Otus) with further insight into the spatio-temporal origin of the Malagasy avifauna. BMC Evolutionary Biology, 2008, 8, pp.197-111. 10.1186/1471-2148-8-197. doi:10.1186/1471-2148-8-197

HAL Id: hal-00319403

https://hal.science/hal-00319403v1

Submitted on 25 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
Research article

Open Access

# Tracing the colonization history of the Indian Ocean scops-owls (Strigiformes: Otus) with further insight into the spatio-temporal origin of the Malagasy avifauna

Jérôme Fuchs*¹,²,³,⁴, Jean-Marc Pons¹,², Steven M Goodman⁵,⁶, Vincent Bretagnolle⁷, Martim Melo⁸, Rauri CK Bowie⁴, David Currie⁹, Roger Safford¹⁰, Munir Z Virani¹¹,¹²,¹³, Simon Thomsett¹¹,¹³, Alawi Hija¹⁴, Corinne Cruaud¹⁵ and Eric Pasquet¹,²

Address: ¹UMR5202 «Origine, Structure et Evolution de la Biodiversité», Département Systématique et Evolution, Muséum National d'Histoire Naturelle, 55 Rue Buffon, 75005 Paris, France, ²Service Commun de Systématique Moléculaire, IFR CNRS 101, Muséum National d'Histoire Naturelle, 43, rue Cuvier, 75005 Paris, France, ³DST-NRF Centre of Excellence at the Percy FitzPatrick Institute, University of Cape Town, Rondebosch 7701, Cape Town, Republic of South Africa, ⁴Museum of Vertebrate Zoology and Department of Integrative Biology, 3101 Valley Life Science Building, University of California, Berkeley, CA 94720-3160, USA, ⁵Field Museum of Natural History, 1400 South Lake Shore Drive, Chicago, IL 60605, USA, ⁶Vahatra, BP3972, Antananarivo (101), Madagascar, ⁷CEBC-CNRS, Beauvoir sur Niort, 79360, France, ⁸UMR5175 Centre d'Ecologie Fonctionnelle et Evolutive, 1919 Route de Mende, F-34293 Montpellier Cedex 5, France, ⁹Nature Seychelles, PO Box 1310, Victoria, Mahé, Republic of Seychelles, ¹⁰BirdLife International, Wellbrook Court, Girton Road, Cambridge CB3 0NA, UK, ¹¹The Peregrine Fund, 5668 West Flying Hawk Lane, Boise Idaho 83709, USA, ¹²Department of Biology, Leicester University, LE1 7RH, UK, ¹³Department of Ornithology, National Museums of Kenya, P.O. Box 40658, Nairobi, Republic of Kenya, ¹⁴Department of Environment, Ministry of Agriculture, Natural Resources, Environment and Cooperatives, Zanzibar Revolutionary Government, P.O. Box 811, Zanzibar, United Republic of Tanzania and ¹⁵Genoscope, Centre National de Séquençage, 2 rue Gaston Crémieux, CP5706, 91057 Evry Cedex, France

Email: Jérôme Fuchs* - jerome.fuchs@berkeley.edu; Jean-Marc Pons - pons@mnhn.fr; Steven M Goodman - SGoodman@wwf.mg; Vincent Bretagnolle - breta@cebc.cnrs.fr; Martim Melo - Martim.PINHEIRO-DE-MELO@cefe.cnrs.fr; Rauri CK Bowie - bowie@berkeley.edu; David Currie - kirtlands@batelnet.bs; Roger Safford - Roger.Safford@birdlife.org; Munir Z Virani - tpf@africaonline.co.ke; Simon Thomsett - sthomsett@gmail.com; Alawi Hija - alawihija99@hotmail.com; Corinne Cruaud - craud@genoscope.cns.fr; Eric Pasquet - pasquet@mnhn.fr

Published: 9 July 2008

Received: 14 January 2008

BMC Evolutionary Biology 2008, 8:197 doi:10.1186/1471-2148-8-197

Accepted: 9 July 2008

This article is available from: http://www.biomedcentral.com/1471-2148/8/197

© 2008 Fuchs et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## Abstract

**Background:** The island of Madagascar and surrounding volcanic and coralline islands are considered to form a biodiversity hotspot with large numbers of unique taxa. The origin of this endemic fauna can be explained by two different factors: vicariance or over-water-dispersal. Deciphering which factor explains the current distributional pattern of a given taxonomic group requires robust phylogenies as well as estimates of divergence times. The lineage of Indian Ocean scops-owls (*Otus*: Strigidae) includes six or seven species that are endemic to Madagascar and portions of the Comoros and Seychelles archipelagos; little is known about the species limits, biogeographic affinities and relationships to each other. In the present study, using DNA sequence data gathered from six loci, we examine the biogeographic history of the Indian Ocean scops-owls. We also compare the pattern and timing of colonization of the Indian Ocean islands by scops-owls with divergence times already proposed for other bird taxa.

**Results:** Our analyses revealed that Indian Ocean islands scops-owls do not form a monophyletic assemblage: the Seychelles *Otus insularis* is genetically closer to the South-East Asian endemic *O*.

<footnote>* Corresponding author</footnote>
sunia than to species from the Comoros and Madagascar. The Pemba Scops-owls *O. pembaensis*, often considered closely related to, if not conspecific with *O. rutilus* of Madagascar, is instead closely related to the African mainland *O. senegalensis*. Relationships among the Indian Ocean taxa from the Comoros and Madagascar are unresolved, despite the analysis of over 4000 bp, suggesting a diversification burst after the initial colonization event. We also highlight one case of putative back-colonization to the Asian mainland from an island ancestor (*O. sunia*). Our divergence date estimates, using a Bayesian relaxed clock method, suggest that all these events occurred during the last 3.6 myr; albeit colonization of the Indian Ocean islands were not synchronous, *O. pembaensis* diverged from *O. senegalensis* about 1.7 mya while species from Madagascar and the Comoro diverged from their continental sister-group about 3.6 mya. We highlight that our estimates coincide with estimates of diversification from other bird lineages.

**Conclusion:** Our analyses revealed the occurrence of multiple synchronous colonization events of the Indian Ocean islands by scops-owls, at a time when faunistic exchanges involving Madagascar was common as a result of lowered sea-level that would have allowed the formation of stepping-stone islands. Patterns of diversification that emerged from the scops-owls data are: 1) a star-like pattern concerning the order of colonization of the Indian Ocean islands and 2) the high genetic distinctiveness among all Indian Ocean taxa, reinforcing their recognition as distinct species.

**Background**

The island of Madagascar is considered a biodiversity hotspot with an intriguing endemic fauna [1]. The origin of the island's peculiar and highly unique fauna can be explained by two contrasting processes: vicariance when Madagascar became separated from the African landmass approximately 165 mya and from India approximately 88 mya [2,3], or over-water dispersal from the African, Australian and Eurasian landmasses. In near proximity to Madagascar are several archipelagos, which have different geological histories. The eastern islands in the Seychelles archipelago are of granitic origin and likely broke off from India and reached their current position some 55–75 mya when India drifted northwards [4], while the western portion of the Seychelles is comprised of recent atolls. In contrast, the volcanic Comoros archipelago is of relatively recent age (0–11 mya), hence the only plausible explanation for the colonization of its biota is by over-ocean dispersal from Africa, Australia, Madagascar or Eurasia.

The origin, timing and modes (vicariance versus over-water dispersal) of colonization by animals of Madagascar and in some cases the surrounding Comoro islands have been the focus of several research programs [5-10]. Given the antiquity of the break-up of Indo-Madagascar from Africa, which considerably pre-dates the first known fossil records of most modern families or genera of animals, in particular vertebrates, it is now assumed that the majority of extant groups arrived on Madagascar via over-water dispersal [11]. From recent phylogenetic studies, as well as traditional taxonomy, it is inferred that most of the Malagasy avifauna originated from African ancestors, although colonization events from Eurasia and Australia have also been documented for flying organisms such as bats and birds [12-16].

Interpreting the origin of certain components of the Malagasy avifauna is in many cases difficult when based on current taxonomical classifications alone as several genera of birds are shared between Madagascar, Africa and Eurasia. Among these are scops-owls of the genus *Otus*. As currently defined, *Otus* is present in five biogeographic areas (Indo-Malaya, Afrotropics, Nearctic, Neotropics, Palearctic) [17]. However, monophyly of *Otus* (sensu [17]) is uncertain; some molecular studies suggested that the New World *Otus* (*Megascops*-including *Otus flammeolus*, see [18]), which differ from Old World *Otus* by song type, are genetically more closely related to the widespread owl genera *Strix* and *Bubo* [19,20] whereas the African White-faced Owl (*O. leucotis*) has closer affinities with the genus *Asio* [20]. The above-mentioned taxa excluded, *Otus* species have their center of diversity in Eurasia (26 species); secondary radiations occur in the Indian Ocean (six or seven species) and Africa (four species). To our knowledge, the earliest known *Otus* fossil, a distal end of the right humerus, is from western Kenya and dates from the Miocene (16.5–18.5 mya; [21]). This partial fossil is morphologically close to *O. senegalensis*, but its relationship among members of the genus is not clear. The western Indian Ocean taxa are often thought to constitute a superfamily (*Rutilus* group) of five or six species [22]: *O. capnodes* from Anjouan, *O. mayottensis* from Mayotte, *O. moheliensis* from Mohéli, *O. pauliani* from Grande Comore, and *O. rutilus*/*O. madagascariensis* from Madagascar (see [23] for a discussion about phylogeography and taxonomic status of the two Malagasy forms), and *O. insularis* from the granitic Seychelles. The latter species shows apparent affinities, based on vocalization data, to the Indonesian *O. magicus* [24]. Most of the western Indian Ocean taxa are poorly known: specific status has been proposed only within the last twenty years for *Otus*
pauliani, O. capnodes, O. madagascariensis and O. mayottensis using both biometric and vocalization data [25-28], while O. moheliensis was first described in 1998 [29].

The evolutionary and biogeographic history of the Indian Ocean *Otus* taxa has not been the focus of a phylogenetic study. In this paper, we use two nuclear introns (myoglobin intron-2 and TGFB2 intron-5) and four mitochondrial protein coding-genes (ND2, ND3, ATP6, cytochrome-*b*) in order to propose a first multi-locus phylogeny of scops-owls and to track their colonization history of the western Indian Ocean islands in space and time. We additionally compare the biogeographic affinities of Malagasy scops-owl species to those of other Malagasy avian lineages which have been the focus of recent genetic studies [12,13,16].

## Results

### Sequence properties

We obtained between 680 and 726 bp (*Otus letitia ussuriensis* and *Aegolius acadicus*, respectively) for the myoglobin intron-2 resulting in a final alignment of 749 bp. Among the 749 bp, 151 were variable (20%) and 56 were parsimony informative (7.5%). Maximum Likelihood (ML) analyses yielded one tree (-ln = 2198.81) that slightly differs topologically from the 50% majority consensus rule tree obtained from the Bayesian analyses (-ln = 2442.98). We obtained between 561 and 593 bp (O. leucotis and O. pembanesis/O. senegalensis Allele 1, respectively) for the TGFB2 intron-5 resulting in a final alignment of 605 bp. Three individuals were found to be length-variant heterozygotes: O. leucotis possesses a CCT duplication in a region with a CCT pattern in all other species, O. rutilus (FMNH 431150) possesses a one base pair deletion (G) in position 171 of the alignment, and O. senegalensis possesses a one base pair insertion (A in position 559 of our alignment; this insertion was also found in the two O. pembanensis individuals sequenced). The two O. senegalensis alleles also differ by two further mutations; these two alleles clustered together as the sister-group to O. pembaensis in a ML analysis (tree not shown). Therefore, we use the consensus sequence (the two single nucleotide polymorphisms were coded using the appropriate IUPAC code) from the two O. senegalensis alleles for further phylogenetic analyses. Only the alleles without the insertion/deletion were included in the phylogenetic analyses for O. leucotis and O. rutilus, as the insertion/deletion events were autapomorphic in both cases. Among the 601 base pairs retained for the analyses, 163 were variable (27%) and 71 were parsimony informative (11.8%). ML analyses yielded one tree (-ln = 2131.71) that slightly differs from the 50% majority consensus rule tree obtained from the Bayesian analyses (-ln = 2179.36).

The topologies obtained from the nuclear loci were very similar to each other, delineating the primary clades without achieving resolution at the tips (see Additional Files 1 and 2). The 50% majority consensus rule tree obtained from the Bayesian analyses (Figure 1, -ln = 4638.29) and Maximum Parsimony strict consensus tree (105675 equally parsimonious trees of 412 steps, CI = 0.82, RI = 0.86) of the two concatenated nuclear loci provided a well-resolved topology for inter-generic relationships as well as some resolution of relationships among the primary *Otus* lineages but failed to provide resolution among members of the Indian Ocean radiation.

The concatenated mitochondrial sequences retained for analyses were 2983 bp long (1047 bp for ND2, 684 bp for ATP6, 351 bp for ND3 and 901 bp for cytochrome-b) and correspond to the positions 5246 to 6281 (ND2), 9240 to 9923 (ATP6), 10776 to 11120 (ND3), and 15011 to 15911 (cytochrome-b) of the Gallus gallus mitochondrial genome sequence [30]; GenBank accession number X52392). The ATP6 and cytochrome-b sequences contained no insertions, deletions and stop-codons in the reading frame. The ND2 sequence of Aegolius acadicus exhibits a two-codon insertion (CAA ACC) just before the stop codon. All the ND3 sequences exhibited the pyrimidine insertion (T for O. capnodes, C for all other species analyzed) previously reported for several clades of birds [31]; this extra-nucleotide was removed before phylogenetic analyses. Partitioning the gene by codon positions significantly improved the fit of models to the data for all four mitochondrial loci, as inferred from the Bayes Factor (BF) values (BFND2 = 826.3, BFATP6 = 681.7, BFND3 = 319.1, BFcytb = 1048.8). Mitochondrial gene trees were very similar to each other (ND2: ML -ln = 10556.45, BI partitioned by codon position = 10161.23 – Additional File 3; ATP6: ML -ln = 6550.30, BI partitioned by codon position = 6232.87 – Additional File 4; ND3: ML -ln = 3266.29, BI partitioned by codon position -ln = 3146.47 – Additional File 5; cytochrome-b, ML -ln = 8623.79, BI partitioned by codon position -ln = 8128.06 – Additional File 6), albeit levels or resolution varied among genes. As expected, no conflict was detected between the individual evolutionary histories of the mitochondrial gene trees (as inferred from posterior probabilities). The 50% majority-rule tree obtained from the concatenated analyses of the mitochondrial genes (partitioned by gene and codon position: -ln = 27725.63, Figure 1) was very similar in terms of topology and number of supported nodes to the MP strict consensus tree (two equally most parsimonious trees of 5922 steps, CI = 0.42, RI = 0.60).

### Phylogenetic results

The individual trees obtained from the two nuclear introns and mitochondrial data sets were very similar to each other and no incongruence was detected between the
<figure><img src="image_2.png" /><figcaption>**Figure 1**</figcaption></figure>

Fifty percent majority-rule consensus tree resulting from the Bayesian mixed-model analyses of the nuclear (left, arithmetic mean, -ln = 4638.29) and mitochondrial (right, arithmetic mean, -ln = 27725.63) data sets. Values close to nodes represent MP bootstrap percentages and BI posterior probabilities. Grey blocks highlight Indian Ocean taxa. Species between quotes indicate samples for which geographic origin is unknown (captive individuals). Branch lengths of the outgroup (*Tyto alba*) were reduced by a scale of two for graphical purpose.

nuclear and mitochondrial topologies (Figure 1), accord-
ing to the criteria defined in the Material and Methods sec-
tion. Further, there was usually strong congruence for
nodal support among the different methods. Most of the
nodes present in the 50% majority-rule consensus tree
resulting from the Bayesian analyses performed on the
concatenated data set (partitioned by gene and codon
position- 14 partitions-, -ln = 32195.91, Figure 2) were
very well supported in the parsimony analyses too (two
equally most parsimonious trees, 6340 steps, CI = 0.44, RI
= 0.62).

The African species *O. ireneae* was recovered as the sister-
taxon of all remaining *Otus sensu stricto* scops-owls. All the
western Indian Ocean taxa clustered in a clade that also
contains the Eurasian *O. scops*, the São Tomé endemic *O.
hartlaubi*, the Pemba Island endemic *O. pembaensis*, the

African mainland *O. senegalensis*, the Philippine taxa *O. mirus* and *O. longicornis*, as well as the Indo-Malayan *O. sunia*. The western Indian Ocean taxa were not recovered as a monophyletic lineage since the Seychelles *O. insularis* was more closely related to the Indo-Malayan *O. sunia* than to any other taxon occurring in the western Indian Ocean region. Uncorrected-p mitochondrial distances among members of the Indian Ocean taxa/*O. sunia* clade range between 4.6% (between *O. capnodes* and *O. rutilus*) and 7.0% (*O. sunia* and *O. rutilus*), with a mean of 5.3% (s.d. = 0.7%). Most of the relationships between the Comorian and Malagasy taxa (*capnodes*, *mayottensis*, *mohe-liensis*, *pauliani*, *rutilus*) did not receive statistical support and short inter-nodes characterized most of the branches among these lineages. We attribute this lack of resolution to rapid speciation events ('hard polytomy') rather than a lack of sufficient character sampling ('soft polytomy')
<figure><img src="image_3.png" /><figcaption>**Figure 2**</figcaption></figure>

Fifty percent majority-rule consensus tree resulting from the mixed-model analyses of the concatenated data set (14 partitions, arithmetic mean, -ln = 32195.91). Values close to nodes represent MP bootstrap percentages and BI posterior probabilities. Note the occurrence of color morphs in most of the *Otus sensu stricto* species. Colors for *Otus* taxa names refer to geographic distribution (green: South-East Asia, red: Africa and Blue: Indian Ocean Islands). The genera *Ptilopsis* (African White-faced Owl) and *Megascops* (New World Screech Owls) refer to taxa that were previously included in *Otus*. Species between quotes indicate samples for which geographic origin is unknown (captive individuals). The branch length of the outgroup (*Tyto alba*) was reduced by a scale of two for graphical purpose. Pictures were modified from [32].
because 1) we sampled several genes with different evolutionary properties resulting in a final alignment of more than 4300 bp and, 2) the nodes above and below the polytomy received bootstrap percentages of 70% or more or posterior probabilities of 0.95 or greater. The remaining *Otus* species clustered in a second large clade. Within this latter clade, *O. lettia* *ussuriensis* (eastern Russia) did not cluster with another *O. lettia* sample from Laos, but with the Philippine taxon *O. megalotis*, suggesting that further work with more complete geographic sampling is needed to address the evolutionary history of the *O. lettia/O. megalotis* species complex.

## Dating analyses

The biogeographic history inferred from the topology of the concatenated analyses is intriguing as several faunistic exchanges involving the Indian Ocean islands and Indo-Malaya region occurred (split *O. mirus*-*O. longicornis* from the Indian Ocean taxa and split between *O. insularis* and *O. sunia*). These faunistic exchanges imply either multiple colonization events of the western Indian Ocean islands or one re-colonization of the mainland by *O. sunia*, about 0.25–0.30 mya. The main radiation of western Indian Ocean island taxa (2.5 mya, 95% HPD = 1.2–4.0 Table 1) occurred soon after the initial colonization event (3.6 mya, 95% HPD = 1.8–6.0). Our analyses revealed that *O. hartlaubi*, endemic to the African Atlantic island of São Tomé, and *O. pembaensis*, restricted to Pemba Island off the east African coast, have strong affinities with the African mainland species *O. senegalensis*. A very close relationship between *O. senegalensis* and *O. pembaensis* is further supported by the fact that these two taxa share a one nucleotide insertion in the TGFB2 locus. The position of *O. harlaubi* in the mitochondrial and concatenated tree renders this insertion paraphyletic. Yet, it is worth noting that the two alleles of *O. senegalensis* differ in length at this site. Considering that the first colonizers of São Tomé and Pemba were probably in small numbers when compared to the continental and widely distributed *O. senegalensis*, we regard the discrepancy between the mitochondrial and nuclear trees as being due to the different effective popu-

<table><caption>Table 1: Posterior distribution of divergence times for some selected nodes.</caption><thead><tr><th>Calibration point</th><th>mya (95% HPD)</th></tr></thead><tbody><tr><td><i>O. sunia</i>/<i>O. insularis</i></td><td>1.1 (0.3–2.0)</td></tr><tr><td><i>Otus pauliani</i>/<i>O. moheliensis</i> *</td><td>0.7 (0.5–0.9)</td></tr><tr><td><i>Otus</i> Indian Ocean radiation</td><td>2.5 (1.2–4.0)</td></tr><tr><td><i>Otus</i> Indian Ocean colonization</td><td>3.6 (1.8–6.0)</td></tr><tr><td>Pemba island colonization</td><td>1.7 (0.7–3.3)</td></tr><tr><td>Clade 1/Clade 2</td><td>9.3 (4.7–15.1)</td></tr><tr><td><i>O. ireneae</i>/remaining <i>Otus</i></td><td>11.7 (6.0–19.0)</td></tr><tr><td>Asioninae/sister-group *</td><td>17.3 (16.7–19.3)</td></tr></tbody></table>

mya: million years ago. * denote the two calibration points.

lation size of the markers (the mitochondrial genome has an effective population size that is one fourth of the nuclear genome) and random processes (coalescence). The basal split within *Otus sensu stricto* occurred about 11.7 mya (95% HPD = 6.0–19.0) and the divergence between the two primary clades about 9.3 mya (95% HPD = 4.7–15.1, Figure 2).

## Discussion

### Phylogeny of *Otus sensu stricto*

Previous studies [18,20] highlighted that several species usually assigned to *Otus* (*Ptilopsis leucotis*, *Megascops*) are not directly related to *Otus*. Our analyses, using independent samples and additional genes, confirmed the lack of direct relationships among these three lineages. Given that these results were previously suggested, we do not discuss them further.

Our analyses based on nuclear and mitochondrial sequence data provide the first phylogenetic hypothesis on the diversification of *Otus sensu stricto*. We suggest that *Otus* consists of at least three primary lineages. The first lineage consists of the 'relictual' Sokoke Scops-owl (*O. ireneae*), endemic to the coastal forests of Kenya and Tanzania. This species is sometimes considered to form a superspecies with the patchily distributed Sandy Scops-owl (*O. icterorhynchus*) [22], the only African *Otus* species we were not able to sample.

The second major lineage is restricted to South-East Asia and consists of *O. spilocephalus*, the *O. bakkamoena*-*O. lettia*-*O. lempiji* superspecies, as well as *Otus lettia* *ussuriensis* and *O. megalotis*. As inferred from our results, the current taxonomy within this group appears to be problematic as *O. lettia* *lettia* and *O. l. ussuriensis* are not recovered as sister-taxa. *O. l. ussuriensis*, restricted to Sakhalin, Ussuriland and North-East China, is sometimes considered to be related to the Japanese Scops-owl (*O. semitorques*), based on voice, plumage and iris color data [22]. We could not include the latter species in the present work, but highlight here that this relationship needs to be further tested using molecular data.

The third major lineage includes all the remaining *Otus* species we sampled and is divided in two subclades that consist of 1) *O. scops*, *O. senegalensis*, *O. hartlaubi* and *O. pembaensis*, and 2) *O. longicornis*, *O. mirus*, *O. mayottensis*, *O. rutilus*, *O. capnodes*, *O. moheliensis*, *O. pauliani*, *O. insularis* and *O. sunia*. Relationships among these taxa are discussed below.

### Phylogeny and origin of the western Indian Ocean *Otus*

Data accumulated in recent years on the biology and distribution of the western Indian Ocean scops-owls [26-29] together with our analyses based on DNA sequences of six
loci, shed new light on the evolutionary history of these
birds. The phylogeny we propose here is in full agreement
with recent work based on phenotypic characters suggest-
ing that O. pauliani, O. capnodes, O. moheliensis and O.
mayottensis represent distinct evolutionary lineages
[25,28,29]. In addition we find good support for Asian
biogeographic affinities of the western Indian Ocean
islands Otus spp., and highlight that they may constitute a
paraphyletic assemblage.

The Indo-Malayan *O. sunia* formed a well-supported clade
with *O. insularis* from the granitic Seychelles, which was
very closely related to the Indian Ocean lineage, the latter
forming a sequentially paraphyletic assemblage. The
nested position of *O. sunia/O. insularis* within the Mala-
gasy-Comorian clade, supported by a posterior probabil-
ity of 0.98 (but not by MP bootstrap percentage), suggests
a recent re-colonization of the mainland from an island-
distributed ancestor. Yet, we also acknowledge that poste-
rior probabilities could be misleading when short inter-
nodes/polytomies are involved [33], which is the case
here, and we await the implementations of reversible-
jump Markov chain Monte Carlo algorithms to explore
these aspects [33].

It is also worth noting that in a similar biogeographic
comparison of a different group of birds, non-monophyly
of western Indian Ocean taxa was also retrieved by Warren
et al. [13]; the Indo-Malayan bulbul *Hypsipetes madagas-
cariensis* was nested within the western Indian Ocean taxa,
suggesting a similar scenario to the one we present for
*Otus*. Colonization of the mainland from an island dis-
tributed ancestor is generally regarded as unlikely
because: 1) mainland taxa are considered more competi-
tive than island taxa and, 2) insular populations are
smaller and produce fewer emigrants compared to those
on continents [34,35]. However, empirical cases of conti-
nental re-colonization from islands are accumulating
[16,36-38], including cases from Madagascar [39]; this
suggests that islands can also act as colonization sources
for continental faunas. Data on the phylogeographic
structure within the Indo-Malaya *O. sunia* complex may
help to decipher which hypothesis (multiple coloniza-
tions of the Indian Ocean islands or re-colonization of the
mainland) best explains the current pattern. Indeed, if
strong and ancient phylogeographic structure occurs
within *O. sunia*, we would expect two independent colo-
nizations of the western Indian Ocean islands. In contrast,
if genetic data indicate weak differentiation among *O.
sunia* populations and patterns of population expansion
from western to eastern Indo-Malaya, this would support
the hypothesis of re-colonization of the mainland from
the Seychelles Islands.

Vocalisations as a tool to infer evolutionary relationships
among scops-owls

The genus Otus is fairly homogeneous in plumage relative to many bird genera. However, the Indian Ocean taxa show strong differences from each other in structural, plumage and vocal characters [26,28]. Considering, for example, only the songs of the Comorian taxa (which inhabit islands as close as 50 km from each other), O. pauliani gives a very long series of chaw notes repeated at about 2/sec; O. moheliensis a sequence of hisses; O. capnodes a high-pitched whistled peeooee; and O. mayottensis a series usually of 3–11 deep, single hoots. Intuitively, these differences may be used to argue against close evolutionary relationships; however, our data indicate that the relationship is indeed close in all cases, as could be predicted from the islands' proximity to one another. This study confirms that the vocal and morphological differences are indeed associated with distinct evolutionary lineages, but suggests that they are not related in any simple or obvious way to the evolutionary distance between these lineages, and therefore must be used with caution in identifying affinities between taxa (or lineages). Therefore, we highlight here that the close relationship suggested by Marshall [24] between O. insularis and O. magicus should be further tested using molecular data, especially if we consider the considerable distance between the two areas (over 6000 km).

Comparison with the biogeographic history of other avian
lineages that colonised the Indian Ocean islands

The geographic and temporal origins of certain western
Indian Ocean island bird taxa have received attention in
recent years [e.g. [5,9,10,13,16,40-45] this study] (Table
2).

Our estimate (3.6 mya, 95% HPD: 1.8–6.0) for the timing
of colonization of the Indian Ocean islands by *Otus* coin-
cides with estimates of at least seven other lineages of
birds (Table 2), suggesting that the Indian Ocean islands
avifauna was highly enriched at that time. The period
associated with these multiple independent colonizations
corresponds with the emergence of the volcanic islands of
the Comoros archipelago (Mohéli 5 mya, Anjouan 11.5–
3.9 mya, Mayotte 11.5–7.7 mya and Grande Comore 0.5
mya; [50,51], thus providing possible stepping stones for
dispersal between Africa and Madagascar. This possibility
is also highlighted by the fact that several species or pop-
ulations with African biogeographic affinities (*Streptopelia
capicola*, *Turtur tympanistria*, *Turdus bewsheri*) [53-55], col-
onized the Comoros islands but the colonization of
Madagascar has not yet been achieved.

All faunal exchanges that unambiguously involve Mada-
gascar and Indo-Malaya, or Madagascar and the Seychelles
occurred during the last 3.5 mya [[13,49], this study]
<table><caption>Table 2: Summary of divergence dates and geographic origins involving Comorian and/or Malagasy taxa (Note that dating methods and calibration points vary among the studies).</caption><thead><tr><td>Taxa</td><td>Geographic origin</td><td>Date estimate</td><td>Reference</td></tr></thead><tbody><tr><td>Philepittinae</td><td>Africa or Indo-Malaya</td><td>41.2 ± 2.5 mya</td><td>[40,41]</td></tr><tr><td>Vangidae</td><td>Africa or Indo-Malaya</td><td>(1) 19.7 mya, (16.8–27.0); (2) 28 ± 4.0 mya</td><td>[10,42]</td></tr><tr><td>Bernieridae</td><td>Africa</td><td>25.2 mya (21.4–31.7)</td><td>[42]</td></tr><tr><td>Streptopelia picturata/ Nesoenas mayeri</td><td>Africa</td><td>20.8 ± 4.1 mya</td><td>[43]</td></tr><tr><td>Coracina cinerea</td><td>Australasia</td><td>18.8 ± 2.2 mya</td><td>[16,43]</td></tr><tr><td>Hartlaubius auratus</td><td>Africa or Indo-Malaya</td><td>12.9–17 mya</td><td>[44]</td></tr><tr><td>Ispidina madagascariensis</td><td>Africa</td><td>13.5 ± 2.6 mya</td><td>[9,45]</td></tr><tr><td>Corythornis cristata</td><td>Africa</td><td>5.5 ± 1.1 mya</td><td>[9,45]</td></tr><tr><td>Indian Ocean Dicrurus</td><td>Africa</td><td>4.7 mya (95% CI: 2.7–7.4)</td><td>[46]</td></tr><tr><td>Motacilla flaviventris</td><td>Africa</td><td>4.5 ± 0.3 mya</td><td>[47]</td></tr><tr><td>Otus</td><td>Indo-Malaya</td><td>3.6 mya (95% HPD: 1.8–6.0)</td><td>This study</td></tr><tr><td>Nectarinia souimanga clade</td><td>Africa</td><td>1.9–3.9 mya</td><td>[48]</td></tr><tr><td>Nectarinia notata clade</td><td>Africa</td><td>1.5–3.5 mya</td><td>[48]</td></tr><tr><td>Zosterops borbonicus lineage</td><td>Indo-Malaya</td><td>1.8 mya</td><td>[49]</td></tr><tr><td>Indian Ocean Hypsipetes</td><td>Indo-Malaya</td><td>0.6–2.6 mya</td><td>[13]</td></tr><tr><td>Zosterops maderaspatanus lineage</td><td>Africa</td><td>1.2 mya</td><td>[49]</td></tr><tr><td>Anas</td><td>Australasia</td><td>No divergence date</td><td>[50]</td></tr></tbody></table>

mya: million years ago

(note that the biogeographic history of the Philepittinae, Vangidae and the Sturnidae genus *Hartlaubia* are still uncertain or ambiguous). Warren et al. [48] hypothesised for members of the genus *Zosterops* that the colonizations could have been favored by dramatic sea-level shifts that occurred during the last 2.5 mya, which would have allowed the emergence of currently submerged land-masses between the Seychelles and Madagascar [56-58], implying a 'stepping-stone' model. Our divergence dates estimates are slightly older than those of Warren et al. [49], although confidence intervals are largely overlapping. This biogeographic hypothesis fits the three unambiguous described cases involving Madagascar and Indo-

Malaya, or Madagascar and the Seychelles (*Hypsipetes*, *Otus* and the *Zosterops borbonicus* lineage).

Subfossil remains of three extinct small owl species have been described from the Mascarene Islands (La Réunion, Mauritius, Rodrigues). Based on certain osteological features, these three species have been included in their own genus, *Mascarenotus*, which has been suggested to be derived from *Otus* [59]. The relationships of *Mascarenotus* with respect to the other owl lineages still needs clarification, as this genus could represent a recent and derived off-shot of the western Indian Ocean lineage, possibly a first off-shot of the colonization process from the Sey-
chelles to more westerly regional islands or even an unrelated lineage of owls.

Whereas the continental biogeographic affinities of certain western Indian Ocean island bird taxa seem largely resolved, the timing of colonizations or faunal affinities among volcanic islands are less well understood and the only aspect that is emerging is the absence of a common pattern, whatever the initial geographic origin. Indeed, all the studies that have been conducted so far indicate explosive diversification and a lack of resolution among the Comoros islands species. These data indicate that once the initial colonization was successful on any of the islands in this archipelago, dispersion and then diversification between nearby islands occurred quickly and randomly. The lack of biogeographic structure at the archipelago scale may thus be partly explained by the geographical arrangement of islands and the short inter-islands distances (maximum distance between two islands in the Comoros archipelago is 90 km) that probably favored colonization by alternative routes.

One final factor that could hide common patterns of diversification is unequal rates of extinction and recolonization across lineages. Indeed, not all lineages would face the same risk of extinction on islands as their characteristics (ecological requirements, population size) often considerably differ. For example, middle-sized birds, like owls (70–120 g), may be more prone to extinction on islands than small birds, like sunbirds (12–15 g) [60]. Yet, even if uneven extinctions rates occurred amongst these lineages, it can be concluded that the overall diversification pattern in the Indian Ocean islands is star-like.

## Conclusion

Our analyses revealed the occurrence of multiple synchronous colonization events of the Indian Ocean islands by scops-owls, at a time when faunistic exchanges involving Madagascar was common as a result of lowered sea-level that would have allowed the formation of stepping-stone islands. Patterns of diversification that emerged from the scops-owls data are: 1) a star-like pattern concerning the order of colonization of the Indian Ocean islands and 2) the high genetic distinctiveness among all Indian Ocean taxa, reinforcing their recognition as distinct species.

## Methods

### Taxonomic sampling

We obtained tissue samples from all western Indian Ocean *Otus* taxa, as well as samples from several Indo-Malayan and Afrotropical species (Table 3), focusing on as many super-species complexes as possible (*sensu* [22]). We obtained tissues for seven of these super-species complexes. We were unable to obtain samples of the distinctive *O. rufescens* and *O. sagittatus*, as well as representatives

of four super-species with localised and distant distributions relative to the western Indian Ocean (*brooki/angeliinae* from Sumatra/Java/Borneo; *mantananensis/magicus* from the Lesser Sundas/Philippines/Mollucas; *collari/manadensis/beccarii* from Sangihe/Sulawesi/Biak and *enganensis/alius/umbra* from islands off Sumatra). Since the monograph of Marks et al. [22], one further *Otus* species, *O. thilohoffmani* has been described from Sri Lanka [61]. This species is only known in museum collections by the type specimen, deposited in the National Museum Colombo (Sri Lanka), and based on morphology, has been suggested to be related to either *O. rufescens* or *O. spi locephalus* [61]. As a consequence, we did not have access to a tissue sample of this newly described species. With the exception of *O. magicus*, considered by some authors to include *O. insularis* because of similarities in vocalizations [24], none of the species we were unable to sample have been considered closely related to the Indian Ocean taxa [22]. We included two individuals per color morph for *O. rutilus* (*sensu* [23]). Representatives of the Strigidae genera *Aegolius*, *Athene* and *Glaucidium* (Surniinae), *Asio* (Asioninae), and *Bubo*, *Strix*, New World *'Otus'* (*Megascops*) and the African White-faced Owl ('*Otus' Ptilopsis leucotis*) (Striginae) were included as proximate outgroups. We rooted our trees using sequences from a representative of the Tytonidae (Barn Owl *Tyto alba*), which has been recovered as the sister-group of the Strigidae in molecular and morphological analyses [62,63].

### Laboratory procedure

Total DNA was extracted from frozen, EDTA or alcohol preserved tissues (liver, blood, feathers, muscle) using a CTAB-based protocol [64] with an overnight Proteinase K (0,1 mg.ml-1) digestion. ND2, ATP6 and ND3 were amplified and sequenced using primer pairs L5219/ H6313 [65], L9245/H9947 [66] and L10755/H11151 [67], respectively. A 900 bp portion of cytochrome-b was amplified with the primer pairs L14967-H15487 and L15424-H15916 [23]). Myoglobin intron-2 and TGFB2 intron-5 were amplified with primers Myo2/Myo3F [68,69] and tgf5/tgf6 [70], respectively. The amplification and sequencing protocol were standard [23].

### Phylogenetic analyses

Molecular phylogenies were estimated using parsimony (P) and model-based approaches (maximum likelihood [ML], and Bayesian inferences [BI]), as implemented in *PHYML* v2.4 [71] and *MRBAYES* 3.1 [72-74]. Parsimony analyses were conducted with *PAUP* v4.0b10 [75] using the heuristic tree bisection and reconnection branch-swapping (TBR) algorithm with 100 random addition replicates. Likelihood models were estimated with *MRMODELTEST* 2.0 [76] using the Akaike Information Criterion [77]. The selected models are listed in Table 4. Clade support in the individual gene trees for the ML and
<table><caption>Table 3: List of samples used and GenBank accession numbers for the six loci analysed.</caption><thead><tr><th>Species</th><th>Voucher/Tissue number</th><th>Geographic</th><th>Myoglobin</th><th>TGFB2</th><th>Cytochrome-b</th><th>ND2</th><th>ATP6</th><th>ND3</th></tr></thead><tbody><tr><td>Aegolius acadicus</td><td>MVZ 118707 (T)</td><td>USA</td><td>EU601093</td><td>EU600970</td><td>U89172</td><td>EU601051</td><td>EU601160</td><td>EU601013</td></tr><tr><td>Asio otus</td><td>MVZ 180184 (T)</td><td>USA</td><td>EU601097</td><td>EU600975</td><td>AF082067</td><td>EU601055</td><td>EU601165</td><td>EU601018</td></tr><tr><td>Athene noctua</td><td>MNHN 1995-99 (T)</td><td>France</td><td>EU601089</td><td>EU600966</td><td>AJ003948</td><td>No sequence</td><td>EU601156</td><td>EU601009</td></tr><tr><td>Bubo bubo</td><td>MNHN 24-55 (T)</td><td>France</td><td>EU601069</td><td>EU600949</td><td>AJ003969</td><td>EU601029</td><td>EU601137</td><td>EU600992</td></tr><tr><td>Bubo virginianus</td><td>MVZ 179340 (T)</td><td>USA</td><td>EU601092</td><td>EU600969</td><td>AF168106</td><td>EU601050</td><td>EU601159</td><td>EU601012</td></tr><tr><td>Glaucidium cuculoides</td><td>MNHN 33-9C (JF150, B)</td><td>Laos</td><td>EU601088</td><td>EU600982</td><td>No sequence</td><td>EU601047</td><td>EU601155</td><td>No sequence</td></tr><tr><td>Glaucidium gnoma</td><td>MVZ 179345 (T)</td><td>USA</td><td>EU601094</td><td>EU600972</td><td>AJ003994</td><td>No sequence</td><td>EU601162</td><td>EU601015</td></tr><tr><td>Otus asio</td><td>MVZ 179828 (T)</td><td>USA</td><td>EU601096</td><td>EU600974</td><td>DQ190845</td><td>EU601054</td><td>EU601164</td><td>EU601017</td></tr><tr><td>Otus 'bakkamoena'</td><td>UWBM 67511 (T)</td><td>Captive</td><td>EU601074</td><td>EU600954</td><td>EU601110</td><td>EU601034</td><td>EU601141</td><td>EU600997</td></tr><tr><td>Otus capnodes</td><td>MNHN 30-10J (B)</td><td>Anjouan</td><td>EU601078</td><td>EU600957</td><td>EU601114</td><td>EU601038</td><td>EU601145</td><td>EU601000</td></tr><tr><td>Otus hartlaubi</td><td>MNHH 32-04G (B)</td><td>São Tomé</td><td>EU601072</td><td>EU600952</td><td>EU601108</td><td>EU601032</td><td>EU601139</td><td>EU600995</td></tr><tr><td>Otus hoyi</td><td>ZMUC 114834 (B)</td><td>Bolivia</td><td>EU601061</td><td>EU600942</td><td>EU601103</td><td>EU601024</td><td>EU601130</td><td>EU600985</td></tr><tr><td>Otus insularis</td><td>D. Currie 5H21863 (F)</td><td>Mahé, Seychelles</td><td>EU601059</td><td>EU600940</td><td>EU601101</td><td>EU601022</td><td>EU601128</td><td>EU600983</td></tr><tr><td>Otus insularis</td><td>D. Currie 5H21866 (F)</td><td>Mahé, Seychelles</td><td>EU601060</td><td>EU600941</td><td>EU601102</td><td>EU601023</td><td>EU601129</td><td>EU600984</td></tr><tr><td>Otus ireneae</td><td>MNHN 32-06J (M. Virani C30577) (B)</td><td>Kenya</td><td>EU601077</td><td>EU600956</td><td>EU601113</td><td>EU601037</td><td>EU601144</td><td>EU600999</td></tr><tr><td>Otus kenicottii</td><td>MVZ 182896 (T)</td><td>USA</td><td>EU601095</td><td>EU600973</td><td>DQ190850</td><td>EU601053</td><td>EU601163</td><td>EU601016</td></tr><tr><td>Otus koepckeii</td><td>ZMUC 115283 (B)</td><td>Peru</td><td>EU601062</td><td>EU600943</td><td>EU601104</td><td>EU601025</td><td>EU601131</td><td>EU600986</td></tr><tr><td>Otus 'lempiji'</td><td>UWBM 73860 (T)</td><td>Captive</td><td>EU601076</td><td>EU600981</td><td>EU601112</td><td>EU601036</td><td>EU601143</td><td>No sequence</td></tr><tr><td>Otus lettia lettia</td><td>MNHN 33-4C (JF142, B)</td><td>Laos</td><td>EU601073</td><td>EU600953</td><td>EU601109</td><td>EU601033</td><td>EU601140</td><td>EU600996</td></tr><tr><td>Otus lettia ussuriensis</td><td>UWBM 75379 (T)</td><td>Russia</td><td>EU601075</td><td>EU600955</td><td>EU601111</td><td>EU601035</td><td>EU601142</td><td>EU600998</td></tr><tr><td>Otus leucotis</td><td>FMNH 429716 (T)</td><td>Congo RD</td><td>EU601085</td><td>EU600963</td><td>EU601120</td><td>EU601044</td><td>EU601152</td><td>EU601006</td></tr><tr><td>Otus longicornis</td><td>FMNH 433020 (T)</td><td>Luzon, Philippines</td><td>EU601084</td><td>EU600962</td><td>EU601119</td><td>EU601043</td><td>EU601151</td><td>EU601005</td></tr><tr><td>Otus longicornis</td><td>ZMUC 114206 (B)</td><td>Isabela, Philippines</td><td>EU601063</td><td>No sequence</td><td>No sequence</td><td>EU601026</td><td>EU601132</td><td>EU600987</td></tr><tr><td>Otus mayottensis</td><td>MNHN R22 (F)</td><td>Mayotte</td><td>EU601087</td><td>EU600965</td><td>EU601122</td><td>EU601046</td><td>EU601154</td><td>EU601008</td></tr><tr><td>Otus megalotis</td><td>FMNH 433019 (T)</td><td>Luzon, Philippines</td><td>EU601083</td><td>EU600961</td><td>EU601118</td><td>EU601041</td><td>EU601150</td><td>EU601004</td></tr><tr><td>Otus megalotis</td><td>ZMUC 114208 (B)</td><td>Isabela, Philippines</td><td>EU601064</td><td>EU600944</td><td>EU601105</td><td>EU601027</td><td>EU601133</td><td>EU600988</td></tr><tr><td>Otus mirus</td><td>FMNH 357429 (T)</td><td>Mindanao, Philippines</td><td>EU601099</td><td>EU600978</td><td>EU601126</td><td>EU601057</td><td>No sequence</td><td>EU601020</td></tr><tr><td>Otus moheliensis</td><td>MNHN E-135 (F)</td><td>Mohéli</td><td>EU601086</td><td>EU600964</td><td>EU601121</td><td>EU601045</td><td>EU601153</td><td>EU601007</td></tr><tr><td>Otus pauliani</td><td>MNHN R24 (F)</td><td>Grande Comore</td><td>EU601100</td><td>EU600979</td><td>EU601125</td><td>EU601058</td><td>No sequence</td><td>EU601021</td></tr><tr><td>Otus pembaensis</td><td>MNHN uncatalogued (B)</td><td>Pemba Island</td><td>EU601090</td><td>EU600967</td><td>EU601123</td><td>EU601048</td><td>EU601157</td><td>EU601010</td></tr></tbody></table>
<table><caption>Table 3: List of samples used and GenBank accession numbers for the six loci analysed. (Continued)</caption><tr><td>Otus pembaensis</td><td>MNHN uncatalogued (B)</td><td>Pemba Island</td><td>EU601091</td><td>EU600968</td><td>EU601124</td><td>EU601049</td><td>EU601158</td><td>EU601011</td></tr><tr><td>Otus roboratus</td><td>ZMUC 114634 (B)</td><td>Ecuador</td><td>EU601065</td><td>EU600945</td><td>EU601106</td><td>EU601028</td><td>EU601134</td><td>No sequence</td></tr><tr><td>Otus rutilus</td><td>FMNH 393149 (T)</td><td>Madagascar</td><td>EU601082</td><td>EU600960</td><td>EF198256</td><td>EF198290</td><td>EU601149</td><td>EU601003</td></tr><tr><td>Otus rutilus</td><td>FMNH 431150/ 431152 (T)</td><td>Madagascar</td><td>EU601068</td><td>EU600948</td><td>EF198273</td><td>EF198307</td><td>EU601136</td><td>EU600991</td></tr><tr><td>Otus rutilus</td><td>FMNH 396240 (T)</td><td>Madagascar</td><td>EU601066</td><td>EU600946</td><td>EF198270</td><td>EF198304</td><td>EU601135</td><td>EU600989</td></tr><tr><td>Otus rutilus</td><td>FMNH 427395 (T)</td><td>Madagascar</td><td>EU601067</td><td>EU600947</td><td>EF198296</td><td>EF198262</td><td>No sequence</td><td>EU600990</td></tr><tr><td>Otus scops</td><td>MNHN 23-5F (T)</td><td>France</td><td>EU601079</td><td>EU600958</td><td>EU601115</td><td>EU601039</td><td>EU601146</td><td>EU601001</td></tr><tr><td>Otus senegalensis</td><td>MVZ uncatalogued (B)</td><td>South Africa</td><td>EU601098</td><td>EU600976, EU600977</td><td>EU601127</td><td>EU601053</td><td>EU601166</td><td>EU601019</td></tr><tr><td>Otus spilocephalus</td><td>MNHN 15-58 (B)</td><td>China</td><td>EU601080</td><td>EU600980</td><td>EU601116</td><td>EU601040</td><td>EU601147</td><td>No sequence</td></tr><tr><td>Otus sunia</td><td>MNHN 6-98 (B)</td><td>Thailand</td><td>EU601081</td><td>EU600959</td><td>EU601117</td><td>EU601041</td><td>EU601148</td><td>EU601002</td></tr><tr><td>Strix aluco</td><td>MNHN CG 1996-114 (T)</td><td>France</td><td>EU601070</td><td>EU600950</td><td>EU601107</td><td>EU601030</td><td>EU601138</td><td>EU600993</td></tr><tr><td>Strix woodfordi</td><td>MNHN 32-01H</td><td>Kenya</td><td>EU601071</td><td>EU600951</td><td>AJ004066</td><td>EU601031</td><td>No sequence</td><td>EU600994</td></tr><tr><td>Tyto alba</td><td>MVZ 180644</td><td>USA</td><td>DQ881879</td><td>EU600971</td><td>AJ004073</td><td>EU601052</td><td>EU601161</td><td>EU601014</td></tr></table>

B refers to blood, F to feather and T to tissue. Taxonomy follows [22]. Geographic origins of the taxa between quotes are unknown. Acronyms for specimen collections include: FMNH – Field Museum of Natural History, Chicago; MNHN – Muséum National d'Histoire Naturelle, Paris; MVZ, Museum of Vertebrate Zoology, Berkeley; UWBM – University of Washington, Burke Museum, Seattle; ZMUC – Zoological Museum, University of Copenhagen.

MP analyses was assessed by 1000 non-parametric bootstrap replicates [78]. The six gene regions sequenced differ considerably in their properties and substitution dynamics, as inferred from the parameters of the models (Table 4). Consequently, analyses of concatenated data set were only performed using a mixed-model strategy. The relevance of partitioning the data set by gene and/or codon position was assessed with the Bayes Factor ($B_F$) [79,80]. Fourteen partitions (myoglobin intron-2, TGFB2 intron-5, first, second and third codon position of each of the four mitochondrial genes) were considered for the Bayesian concatenated analyses according to the functional properties of the markers. Bayesian analyses for the concatenated data set were performed allowing base frequencies, rate matrix, shape parameter and proportion of invariable sites to vary between the partitions (using the unlink and preset commands). Between four and six incrementally heated Metropolis-coupled MCMC chains were run for 15 million generations with trees sampled every 100 generations. The first $2 \times 10^6$ generations (20000 trees) were discarded ('burn-in' period) and the posterior probabilities were estimated from the remaining sampled generations. The default temperature for chain heating (T = 0.2) resulted in not satisfactorily mixing among chains for the concatenated data set; we therefore lowered the temperature to T = 0.05, which resulted in swap frequencies between chains within the 20-70% interval. Two independent Bayesian runs initiated from random starting

trees were performed for each data set, and the log-likelihood values and posterior probabilities were checked to ascertain that the chains had reached convergence. We also checked that the Potential Scale Reduction Factor (PSRF) approached 1.0 for all parameters and that the average standard deviation of split frequencies converged towards zero. We detected significant incongruence between the individual gene trees by comparing the topologies and nodal support obtained under different analytical methods (ML, BI). Criteria for incongruence were set at 70% for the bootstrap values [81], and at 0.95 for posterior probabilities [72].

## Molecular dating analyses

Owls have a very rich fossil record [22,82], yet, the taxonomic history of most of these taxa remains controversial and, to our knowledge, no cladistic analyses including fossil and modern taxa has been conducted to date. This hampers the use of most of these fossils as calibration points in our analyses. The family Strigidae is usually divided into three subfamilies or tribes: Striginae, Surniinae and Asioninae [21]. These three subfamilies are defined by a combination of shared characters [[82] fide [83]]. The least inclusive of these three subfamilies, the Asioninae, consists of three genera and nine species (two of the genera, Pseudoscops and Nesasio are endemic to Jamaica and the Solomon Islands, respectively, whereas Asio is widely distributed). The Asionae differ osteologi-
<table><caption>Table 4: Model selected and parameters values with their 95% credibility intervals when applicable (obtained with MRBAYES).</caption><thead><tr><th></th><th>Myoglobin</th><th>TGFB2</th><th>ND2</th><th>ATP6</th><th>ND3</th><th>Cytochrome-b</th></tr></thead><tbody><tr><td>Model</td><td>K80 + Γ</td><td>HKY + Γ</td><td>GTR + Γ + I</td><td>GTR + Γ + I</td><td>GTR + Γ + I</td><td>GTR + Γ + I</td></tr><tr><td>1<sup>st</sup> position</td><td>NA</td><td>NA</td><td>GTR + Γ + I</td><td>GTR + Γ + I</td><td>SYM + Γ</td><td>GTR + Γ + I</td></tr><tr><td>2<sup>nd</sup> position</td><td>NA</td><td>NA</td><td>GTR + Γ + I</td><td>GTR + Γ + I</td><td>GTR + Γ</td><td>GTR + Γ + I</td></tr><tr><td>3<sup>rd</sup> position</td><td>NA</td><td>NA</td><td>GTR + Γ + I</td><td>GTR + Γ + I</td><td>HKY + Γ + I</td><td>GTR + Γ</td></tr><tr><td>Freq A</td><td>0.25</td><td>0.25 (0.22–0.28)</td><td>0.37 (0.34–0.39)</td><td>0.33 (0.30–0.36)</td><td>0.32 (0.28–0.36)</td><td>0.32 (0.30–0.35)</td></tr><tr><td>Freq C</td><td>0.25</td><td>0.21 (0.18–0.24)</td><td>0.41 (0.39–0.43)</td><td>0.44 (0.41–0.46)</td><td>0.41 (0.37–0.45)</td><td>0.44 (0.42–0.47)</td></tr><tr><td>Freq G</td><td>0.25</td><td>0.24 (0.21–0.27)</td><td>0.06 (0.05–0.07)</td><td>0.07 (0.06–0.08)</td><td>0.09 (0.07–0.11)</td><td>0.08 (0.07–0.09)</td></tr><tr><td>Freq T</td><td>0.25</td><td>0.30 (0.27–0.33)</td><td>0.16 (0.15–0.18)</td><td>0.16 (0.14–0.18)</td><td>0.18 (0.16–0.21)</td><td>0.16 (0.14–0.17)</td></tr><tr><td>A-C</td><td>NA</td><td>NA</td><td>0.016<br>(0.012–0.020)</td><td>0.020<br>(0.013–0.027)</td><td>0.022<br>(0.013–0.032)</td><td>0.014<br>(0.009–0.020)</td></tr><tr><td>A-G</td><td>NA</td><td>NA</td><td>0.663<br>(0.607–0.718)</td><td>0.564<br>(0.484–0.642)</td><td>0.535<br>(0.429–0.637)</td><td>0.549<br>(0.474–0.621)</td></tr><tr><td>A-T</td><td>NA</td><td>NA</td><td>0.021<br>(0.014–0.029)</td><td>0.035<br>(0.022–0.050)</td><td>0.015<br>(0.005–0.029)</td><td>0.032<br>(0.019–0.047)</td></tr><tr><td>C-G</td><td>NA</td><td>NA</td><td>0.005<br>(0.0002–0.014)</td><td>0.024<br>(0.012–0.039)</td><td>0.007<br>(0.0002–0.218)</td><td>0.010<br>(0.0002–0.020)</td></tr><tr><td>C-T</td><td>NA</td><td>NA</td><td>0.232<br>(0.188–0.281)</td><td>0.323<br>(0.252–0.399)</td><td>0.327<br>(0.240–0.425)</td><td>0.321<br>(0.259–0.390)</td></tr><tr><td>G-T</td><td>NA</td><td>NA</td><td>0.062<br>(0.035–0.094)</td><td>0.035<br>(0.010–0.069)</td><td>0.094<br>(0.051–0.146)</td><td>0.074<br>(0.041–0.114)</td></tr><tr><td>Γ</td><td>0.10 (0.09–0.12)</td><td>0.586<br>(0.382–0.936)</td><td>0.997<br>(0.826–1.195)</td><td>0.901<br>(0.676–1.158)</td><td>0.910<br>(0.604–1.295)</td><td>0.869<br>(0.702–1.069)</td></tr><tr><td>I</td><td>NA</td><td>NA</td><td>0.338<br>(0.302–0.374)</td><td>0.400<br>(0.345–0.447)</td><td>0.368<br>(0.287–0.433)</td><td>0.465<br>(0.426–0.501)</td></tr><tr><td>Ts/Tv</td><td>4.24 (3.08–5.75)</td><td>5.90 (4.32–7.795)</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>-ln (ML)</td><td>2198.81</td><td>2131.71</td><td>10556.45</td><td>6550.30</td><td>3266.29</td><td>8623.79</td></tr><tr><td>-ln (BI)</td><td>2442.98</td><td>2179.36</td><td>10582.97</td><td>6585.41</td><td>3305.87</td><td>8658.36</td></tr><tr><td>-ln (BI partitioned)</td><td>NA</td><td>NA</td><td>10161.23</td><td>6232.87</td><td>3146.47</td><td>8128.06</td></tr></tbody></table>

NA means not applicable. Reported likelihood scores for the BI refer to the arithmetic mean.

cally from other owls by having: 1) the anterior rim of the internal trochlea not protruding more anteriorly than the anterior rim of the external trochlea, 2) the external calca- neal ridge bent posteriorly, 3) bony loop broad and 4) tubercle for Musculus tibialis antiquus displaced exter- nally [83]. The most ancient fossil having this combina- tion of characters Intulula tinnipara, has been dated from the Early Miocene (23.7–16.4 mya) [84]. We therefore used this date as a minimum age for the split between Asio otus (the member of the Asioninae we sampled) and its closest relative. As a second calibration point, we used the split between O. pauliani (endemic to Grande Comore, a volcanic island) and its closest relative at 0.5 mya. This date corresponds to the emergence of Grande Comore [51], and is thus the oldest possible age for the coloniza- tion of the O. pauliani lineage. We used these two calibra- tion points in combination.

We used BEAST v1.4.6 [85-87] to estimate the divergence dates within the genus Otus. We assigned the best fitting model, as estimated by MRMODELTEST2, to each of the six loci. We used an exponential distribution for the fossil cal- ibration bound [88]. We set the lower bound of the expo- nential distribution to 16.4 mya, which correspond to the lowest bound of the Early Miocene epoch and the expo-

nential mean to 2.5, so that the 95% distribution proba- bility fell within the 23.7–16.4 mya interval, corresponding to the Early Miocene, which is the epoch for the first Asioninae fossil. For the geological calibration point (emergence of Grande Comore), we used a normal distribution with the mean and standard deviation set to 0.5 mya and 0.1 mya, respectively. We assumed a Yule Speciation Process for the tree prior and an Uncorrelated Lognormal distribution for the molecular clock model [86]. We used default prior distributions for all other parameters and ran MCMC chains for 75 million genera- tions, as the effective sample size for some parameter esti- mates was not large enough using the default length (10 million generations).

## Authors' contributions

JF did part of the laboratory work, edited and aligned the sequences, performed the phylogenetic and dating analy- ses and wrote the first draft of the manuscript. JMP did part of the laboratory work and helped to draft the manu- script. SMG provided tissue samples for the study, wrote several sections of the manuscript and commented on multiple versions. VB collected tissue samples in the field and commented on the manuscript. MM assisted with sample collection and with the draft of the manuscript.
RCKB helped with tissue sample acquisition, provided funding and commented on multiple versions of the manuscript. DC acquired field samples. RS acquired field samples, advised on taxon sampling in relation to morphology and vocalizations, and helped to prepare and finalize the manuscript. MZV collected blood samples and commented on the manuscript. ST helped with collecting samples. AH helped with capturing Otus pembaensis. CC did part of the laboratory work. EP collected tissue samples in the field and commented on the manuscript. All authors have read and approved the final version of the manuscript.

## Additional material

### Additional File 1

Fifty percent majority-rule consensus tree (arithmetic mean -ln = 2442.98) obtained from the Bayesian Inference analyses of myoglobin intron-2 (749 bp). Values next to branches represent MP/ML bootstrap percentages (below) and BI posterior probabilities (above). Gray blocks represent the Indian Ocean taxa. Species between quotes indicate samples for which geographic origin is unknown (captive individuals). The phylogram represents the relationships among owls as inferred from myoglobin intron-2 sequence data.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2148-8-197-S1.pdf]

### Additional File 2

Fifty percent majority-rule consensus tree (arithmetic mean -ln = 3427.03) obtained from the Bayesian Inference analyses TGFB2 intron-5 (602 bp). Values next to branches represent MP/ML bootstrap percentages (below) and BI posterior probabilities (above). Gray blocks represent the Indian Ocean taxa. Species between quotes indicate samples for which geographic origin is unknown (captive individuals). The phylogram represents the relationships among owls as inferred from TGFB2 intron-5 sequence data.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2148-8-197-S2.pdf]

### Additional File 3

Fifty percent majority-rule consensus tree (arithmetic mean -ln = 10161.23) obtained from the Bayesian Inference analyses of the mitochondrial ND2 gene (1041 bp) under a mixed-model strategy (partitioned by codon position). Values next to branches represent MP/ML bootstrap percentages (below) and BI posterior probabilities (above). Gray blocks represent the Indian Ocean taxa. Species between quotes indicate samples for which geographic origin is unknown (captive individuals). The phylogram represents the relationships among owls as inferred from ND2 sequence data.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2148-8-197-S3.pdf]

### Additional File 4

Fifty percent majority-rule consensus tree (arithmetic mean -ln = 6232.87) obtained from the Bayesian Inference analyses of the mitochondrial ATP6 gene (684 bp) under a mixed-model strategy (partitioned by codon position). Values next to branches represent MP/ML bootstrap percentages (below) and BI posterior probabilities (above). Gray blocks represent the Indian Ocean taxa. Species between quotes indicate samples for which geographic origin is unknown (captive individuals). The phylogram represents the relationships among owls as inferred from ATP6 sequence data.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2148-8-197-S4.pdf]

### Additional File 5

Fifty percent majority-rule consensus tree (arithmetic mean -ln = 3146.47) obtained from the Bayesian Inference analyses of the mitochondrial ND3 gene (351 bp) under a mixed-model strategy (partitioned by codon position). Values next to branches represent MP/ML bootstrap percentages (below) and BI posterior probabilities (above). Gray blocks represent the Indian Ocean taxa. Species between quotes indicate samples for which geographic origin is unknown (captive individuals). The phylogram represents the relationships among owls as inferred from ND3 sequence data.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2148-8-197-S5.pdf]

### Additional File 6

Fifty percent majority-rule consensus tree (arithmetic mean -ln = 8128.06) obtained from the Bayesian Inference analyses of the mitochondrial Cytochrome-b gene (1041 bp) under a mixed-model strategy (partitioned by codon position). Values next to branches represent MP/ML bootstrap percentages (below) and BI posterior probabilities (above). Gray blocks represent the Indian Ocean taxa. Species between quotes indicate samples for which geographic origin is unknown (captive individuals). The phylogram represents the relationships among owls as inferred from Cytochrome-b sequence data.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2148-8-197-S6.pdf]

## Acknowledgements

We are very grateful to D. Harper (Leicester University), B. Smith and A. McKechnie (University of Pretoria), S. Birks (University of Washington, Burke Museum), J. Bates, D. Willard and S. Hackett (Field Museum, Chicago), J. Fjeldså and J.-B. Kristensen (Zoological Museum, University of Copenhagen) and all the field collectors for kindly loaning us tissue samples and to G. Mayr (Forschungsinstitut Senckenberg), S. Olson (Smithsonian Institution) and D. Johnson (Global Owl Project) for help with the bibliography of fossil owls. A. Tillier, C. Bonillo and J. Lambourdière provided invaluable help in the laboratory. This work was supported by the Plan Pluriformation 'Etat et structure phylogénétique de la biodiversité actuelle et fossile', the 'Consortium National de Recherche en Génomique', and the 'Service de Systématique Moléculaire' of the Muséum National d'Histoire Naturelle (IFR101). It is part of the agreement n°2005/67 between the Genoscope and the Muséum National d'Histoire Naturelle on the project 'Macrophylogeny of life' directed by Guillaume Lecointre. Fieldwork in the Comoros Islands was funded by a research grant from Ministère de l'Écologie et du Développement Durable ('programme Ecosystèmes Tropicaux');
C. Attié, R. Bruckert, R.-M. Lafontaine and R. Zotier are kindly thanked by EP and VB for assistance in the Comoros. The Phongsaly Forest Conservation and Rural Development Project, a Lao-European cooperation and its staff, P. Rousseau, C. Hatten, Y. Varelides, R. Humphrey and Y. Tipavanh provided invaluable help and assistance during fieldwork in Laos, which was conducted with A. Cibois, M. Ruedi and R. Kirsch. The 'Société des Amis du Muséum' provided fundings to JF for his fieldwork and is sincerely acknowledged. This research was supported by a postdoctoral fellowship to JF from the DST/NRF Centre of Excellence at the Percy FitzPatrick Institute. MM was supported by a post-doctoral grant from the Ministère de la Recherche and the Université de Montpellier II (France).

## References

1. Goodman SM, Benstead JP, [eds]: The natural history of *Madagascar*. Chicago The University of Chicago Press; 2003.

2. Coffin MF, Rabinowitz PD: *The Mesozoic East African and Madagascan conjugate continental margins*. In *Geology and geophysics of continental margins* Edited by: Watkins JS, Zhiqiang F, McMillen KJ. Tulsa, Oklahoma American Association of Petroleum Geologists; 1992:207-246.

3. Storey M, Mahoney JJ, Saunders AD, Duncan RA, Kelley SP, Coffin MC: *Timing of hot-spot related volcanism and the breakup of Madagascar and India*. *Science* 1995, **267**:852-855.

4. Norton IO, Sclater JG: *A model for the evolution of the Indian Ocean and the break-up of Gondwanaland*. *J Geophys Res* 1979, **84**:6803-6830.

5. Cibois A, Schulenberg TS, Pasquet E: *An endemic radiation of Malagasy songbirds is revealed by mitochondrial DNA sequence data*. *Evolution* 2001, **55**:1198-1206.

6. Nágy ZT, Joger U, Wink M, Glaw F, Vences M: *Multiple colonization of Madagascar and Socotra by colubrid snakes: Evidence from nuclear and mitochondrial gene phylogenies*. *Proc R Soc Ser B-Bio* 2003, **270**:2613-2621.

7. Vences M, Vieites DR, Glaw F, Brinkmann H, Kosuch J, Veith M, Meyer A: *Multiple overseas dispersal in amphibians*. *Proc R Soc Ser B-Bio* 2003, **270**:2435-2442.

8. Yoder AD, Burns MM, Zehr S, Delefosse T, Veron G, Goodman SM, Flynn JJ: *Single origin of Malagasy Carnivora from an African ancestor*. *Nature* 2003, **421**:734-737.

9. Marks BD, Willard DE: *Phylogenetic relationships of the Madagascar Pygmy Kingfisher (Ispidina madagascariensis)*. *Auk* 2005, **122**:1271-1280.

10. Fuchs J, Fjeldså J, Pasquet E: *An ancient African radiation of corvoid birds (Aves: Passeriformes) detected by mitochondrial and nuclear sequence data*. *Zool Scr* 2006, **35**:375-385.

11. Yoder AD, Nowak MD: *Has vicariance or dispersal been the predominant biogeographic force in Madagascar? Only time will tell*. *Annu Rev Ecol Evol Syst* 2006, **37**:405-431.

12. Shapiro B, Sibthorpe D, Rambaut A, Austin J, Wragg GM, Bininda-Emons ORP, Lee PLM, Cooper A: *Flight of the dodo*. *Science* 2002, **295**:1683.

13. Warren BH, Bermingham E, Prys-Jones RP, Thébaud C: *Tracking island colonization history and phenotypic shifts in Indian Ocean bulbuls (Hypsipetes: Pycnonotidae)*. *Biol J Linn Soc* 2005, **85**:271-287.

14. Bates PJJ, Ratrimomanarivo FH, Harrison D, Goodman SM: *A review of pipistrelles and serotines (Chiroptera: Vespertilionidae) from Madagascar, including the description of a new species of Pipistrellus*. *Acta Chiropterol* 2006, **8**:299-324.

15. Goodman SM, Cardiff SG, Ranivo J, Russell AL, Yoder AD: *A new species of *Emballonura* (Emballonuridae: Chiroptera) from the dry regions of Madagascar*. *Am Mus Novit* 2006, **3538**:1-24.

16. Fuchs J, Cruaud C, Couloux A, Pasquet E: *Complex biogeographic history of the cuckoo-shrikes and allies (Passeriformes: Campephagidae) revealed by mitochondrial and nuclear sequence data*. *Mol Phylogenet Evol* 2007, **44**:138-153.

17. Dickinson EC, [eds]: *The Howard and Moore complete checklist of the birds of the world* 3rd edition. Princeton New Jersey, Princeton University Press; 2003.

18. Proudfoot GA, Gehlbach FR, Honeycutt RL: *Mitochondrial DNA variation and phylogeography of the Eastern and Western Screech Owls*. *Condor* 2007, **109**:617-627.

19. Miranda HC Jr, Kennedy RS, Mindell DP: *Phylogenetic placement of *Mimizuku gurneyi* (Aves: Strigidae) inferred from mitochondrial DNA*. *Auk* 1997, **114**:315-323.

20. Wink M, Heidrich P: *Molecular systematics and evolution of the Owls*. In *Owls: An identification guide to the owls of the World* Edited by: König C, Becking JH, Weick F, Sussex, Pica Press; 1999:39-58.

21. Walker CA, Dyke GK: *New records of fossil birds of prey from the Miocene of Kenya*. *Historical Biology* 2006, **18**:91-94.

22. Marks JS, Cannings RJ, Mikkola H: *Family Strigidae (typical owls)*. In *Handbook of the birds of the World. Volume 5 Barn Owls to hummingbirds* Edited by: del Hoyo J, Elliott A, Sargatal J. Barcelona, Lynx Editions; 1999:76-243.

23. Fuchs J, Pons J-M, Pasquet E, Raherilalao MJ, Goodman SM: *Geographical structure of the genetic variation in the Malagasy Scops-owl inferred from mitochondrial sequence data*. *Condor* 2007, **109**:408-418.

24. Marshall JT: *Systematics of smaller Asian night birds based on voice*. Ornithological Monographs AOU; 1978.

25. Herremans M, Louette M, Stevens J: *Conservation status and vocal and morphological description of the Grand Comoro Scops Owl Otus pauliani Benson 1960*. *Bird Conserv Int* 1991, **1**:123-133.

26. Safford RJ: *Rediscovery, taxonomy and conservation of the Anjouan Scops Owl Otus capnodes* (Gurney, 1889). *Bird Conserv Int* 1993, **3**:57-74.

27. Lewis A: *Mayotte Scops Owl Otus rutilus mayottensis*. *Bull Afr Bird Club* 1998, **5**:33-34.

28. Rasmussen PC, Schulenberg TS, Hawkins F, Voninavoko R: *Geographic variation in the Malagasy Scops-owl (Otus rutilus auct.)*: The existence of an unrecognized species on Madagascar and the taxonomy of other Indian Ocean taxa. *Bull Br Ornithol Club* 2000, **120**:75-102.

29. Lafontaine R-M, Moulaert N: *Une nouvelle espèce de petit-duc (Otus, Aves) aux Comores: taxonomie et statut de conservation*. *J Afr Zool* 1998, **112**:163-169.

30. Desjardins P, Morais R: *Sequence and gene organization of the chicken mitochondrial genome a novel gene order in higher vertebrates*. *J Mol Biol* 1990, **212**:599-634.

31. Mindell DP, Sorenson MD, Dimcheff DE: *An extra nucleotide is not translated in mitochondrial ND3 of some birds and turkeys*. *Mol Biol Evol* 1998, **15**:1568-1571.

32. del Hoyo J, Elliott A, Sargatal J: *Handbook of the birds of the World. Volume 5 Barn Owls to hummingbirds* Barcelona, Lynx Editions; 1999.

33. Lewis PO, Holder MT, Holsinger KE: *Polytomies and Bayesian phylogenetic inference*. *Syst Biol* 2005, **54**:241-253.

34. MacArthur RH, Wilson EO: *The theory of island biogeography* Princeton New Jersey, Princeton University Press; 1967.

35. Mayr E, Diamond J: *The birds of northern Melanesia: Speciation, ecology, and biogeography* Oxford, Oxford University Press; 2001.

36. Carine MA, Russell SJ, Santos-Guerra A, Francisco-Ortega J: *Relationships of the Macaronesian and Mediterranean floras: Molecular evidence for multiple colonizations into Macaronesia and back colonization of the continent in Convolvulus (Convolvulaceae)*. *Am J Bot* 2004, **91**:1070-1095.

37. Nicholson KE, Glor RE, Kolbe JJ, Larson A, Hedges BS, Lobos JB: *Mainland colonization by island lizards*. *J Biogeogr* 2005, **32**:929-938.

38. Filardi C, Moyle RG: *Single origin of a pan-Pacific bird group and upstream colonization of Australasia*. *Nature* 2006, **438**:216-219.

39. Raxworthy CJ, Forstner MRJ, Nussbaum RA: *Chameleon radiation by oceanic dispersal*. *Nature* 2002, **415**:784-787.

40. Moyle RG, Chesser RT, Prum RO, Schikler P, Cracraft J: *Phylogeny and evolutionary history of Old World Suboscine birds (Aves: Eurylaimides)*. *Am Mus Novit* 2006, **3544**:1-22.

41. Irestedt M, Ohlson JI, Zuccon D, Källersjö M, Ericson PGP: *Nuclear DNA from old collections of avian study skins reveals the evolutionary history of the Old World suboscines (Aves, Passeriformes)*. *Zool Scr* 2006, **35**:567-580.

42. Beresford P, Barker FK, Ryan PG, Crowe TM: *African endemics span the tree of songbirds (Passeri): Molecular systematics of several evolutionary 'enigmas'*. *Proc R Soc Ser B-Bio* 2005, **272**:849-858.
43. Fuchs J: Biogéographie comparée des oiseaux forestiers de l'Ancien Monde: phylogénie et datation moléculaire. In PhD thesis Université Pierre et Marie Curie, Paris; 2006.

44. Zuccon D, Cibois A, Pasquet E, Ericson PGP: Nuclear and mitochondrial sequence data reveal the major lineages of starlings, mynas and related taxa. Mol Phylogenet Evol 2006, 41:333-344.

45. Moyle RG, Fuchs J, Pasquet E, Marks BD: Feeding behavior, toe count and the phylogenetic relationships among alcedine kingfishers. J Avian Biol 2007, 28:317-326.

46. Pasquet E, Pons JM, Fuchs J, Cruaud C, Bretagnolle V: Evolutionary history and biogeography of the drongos (Dicruridae), a tropical Old World clade of corvoid passerines. Mol Phylogenet Evol 2007, 45:158-167.

47. Voelker G: Systematics and historical biogeography of wagtails: dispersal versus vicariance revisited. Condor 2002, 104:725-739.

48. Warren BH, Bermingham E, Bowie RCK, Prys-Jones RP, Thébaud C: Molecular phylogeography reveals island colonization history and diversification of western Indian Ocean sunbirds (Nectarinia: Nectariniidae). Mol Phylogenet Evol 2003, 29:67-85.

49. Warren BH, Bermingham E, Prys-Jones RP, Thébaud C: Immigration, species radiation and extinction in a highly diverse songbird lineage: white-eyes on Indian Ocean islands. Mol Ecol 2006, 15:3769-3786.

50. Johnson KP, Sorenson MD: Phylogeny and biogeography of the dabbling ducks (Genus: Anas): A comparison of molecular and morphological evidence. Auk 1999, 116:792-805.

51. Emerick CM, Duncan RA: Age progressive volcanism in the Comores Archipelago, western Indian Ocean and implications for Somali plate tectonics. Earth Planet Sci Lett 1982, 60:415-428.

52. Nougier J, Cantagrel JM, Karche JP: The Comores archipelago in the western Indian Ocean: Volcanology, geochronology, and geodynamic setting. J Afr Earth Sci 1986, 5:135-145.

53. Bowie RCK, Voelker G, Fjeldså J, Lens L, Hackett SJ, Crowe TM: Systematics of the Olive Thrush Turdus olivaceus species complex with reference to the taxonomic status of the endangered Taita Thrush T. helleri. J Avian Biol 2005, 36:391-404.

54. Voelker G, Rohwer S, Bowie RCK, Outlaw DC: Molecular systematics of a speciose, cosmopolitan songbird genus: Defining the limits of, and relationships among, the Turdus thrushes. Mol Phylogenet Evol 2007, 42(2):422-434.

55. Louette M, Meirte D, Jocqué R: La faune terrestre de l'archipel des Comores. In Studies in Afrotropical Zoology 293 Tervuren, Belgium: Musée Royale de l'Afrique Centrale; 2004.

56. Haq BU, Hardenbol J, Vail PR: Chronology of fluctuating sea levels since the Triassic. Science 1987, 235:1156-1167.

57. Colonna M, Casanova J, Dullo WC, Camoin G: Sea-level changes and 18O record for the past 34 000 yr from Mayotte reef, Indian Ocean. Quaternary Res 1996, 46:335-339.

58. Rohling EJ, Fenton M, Jorissen FJ, Bertrand P, Ganssen G, Caulet JP: Magnitudes of sealevel lowstands of the past 500 000 years. Nature 1998, 394:162-165.

59. Mourer-Chauviré C, Bour R, Moutou F, Ribes S: Mascarenotus nov. gen. (Aves, Strigiformes) genre endémique éteint des Mascareignes et M.grucheti sp. nov., espèce éteinte de la Réunion. C R Acad Sci Paris 1994, 318:1699-1706.

60. Gaston KJ, Blackburn TM: Birds, body size and the threat of extinction. Phil Trans R Soc Lond B 1995, 347:205-212.

61. Warakagoda DH, Rasmussen PC: A new species of scops-owl from Sri Lanka. Bull Br Ornithol Club 2004, 124:85-105.

62. Ericson PGP, Anderson CL, Britton T, Elzanowski A, Johansson US, Källersjö M, Ohlson JI, Parsons TJ, Zuccon D, Mayr G: Diversification of Neoaves: integration of molecular sequence data and fossils. Biol Lett 2006, 2:543-547.

63. Livezey BC, Zusi RL: Higher-order phylogeny of modern birds (Theropoda, Aves: Neornithes) based on comparative anatomy. II. Analysis and discussion. Zool J Linn Soc 2007, 149:1-95.

64. Winnepenninckx B, Backeljau T, De Wachter R: Extraction of high molecular weight DNA from mollusks. Trends Genet 1993, 9:407.

65. Sorenson MD, Ast JC, Dimcheff DE, Yuri T, Mindell DP: Primers for a PCR-based approach to mitochondrial genome sequencing

in birds and other vertebrates. Mol Phylogenet Evol 1999, 12:105-114.

66. Eberhard JR, Bermingham E: Phylogeny and biogeography of the *Amazona ochrocephala* (Aves: Psittacidae) complex. *Auk* 2004, 121:318-332.

67. Chesser RT: Molecular systematics of the rhinocryptid genus *Pterotochos*. Condor 1999, 101:439-446.

68. Slade RW, Moritz C, Heideman A, Hale PT: Rapid assessment of single-copy nuclear DNA variation in diverse species. Mol Ecol 1993, 2:359-373.

69. Heslewood MM, Elphinstone MS, Tidemann SC, Baverstock PR: *Myoglobin* intron variation in the *Gouldian Finch Erythrina gouldiae* assessed by temperature gradient gel electrophoresis. Electrophoresis 1998, 19:142-151.

70. Bures S, Dvornik P, Saetre G-P: Hybridization and apparent hybridization between meadow pipit (*Anthus pratensis*) and water pipit (*A. spinoletta*). Hereditas 2002, 136:254-256.

71. Guindon S, Gascuel O: A simple, fast and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 2003, 52:696-704.

72. Huelsenbeck JP, Ronquist F: MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics 2001, 17:754-755.

73. Huelsenbeck JP, Ronquist F: MrBayes: A program for the Bayesian inference of phylogeny. Version 3.0b4 [http://mrbayes.csit.fsu.edu/].

74. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19:1572-1574.

75. Sofford DP: PAUP*. Phylogenetic analysis using parsimony (* and other methods), version 4.0b10. Sinauer Associates; 2003.

76. Nylander JAA: MrModeltest2.2. [http://www.ebc.uu.se/systzoo/staff/nylander.html].

77. Akaike H: Information theory as an extension of the maximum likelihood principle. In Second international symposium on information theory Edited by: Petrov BN, Csaki F. Budapest, Akademiai Kiado; 1973.

78. Felsenstein J: Confidence limits on phylogenies: An approach using bootstrap. Evolution 1985, 39:783-791.

79. Kass RE, Raferty AE: Bayes factors. J Am Stat Assoc 1995, 90:773-795.

80. Brown JM, Lemmon AR: The importance of data partitioning and the utility of Bayes Factors in Bayesian phylogenetics. Syst Biol 2007, 56:643-665.

81. Hillis DM, Bull JJ: An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst Biol 1993, 42:182-192.

82. Mlíkovský J: Cenozoic birds of the World. Part 1: Europe. Ninox Press, Prague; 2002.

83. Ford NL: A systematic study of the owls based on comparative osteology. In PhD thesis University of Michigan; 1967.

84. Mlíkovský J: Two new owls (Aves: Strigidae) from the Early Miocene of the Czech Republic, with comments on the fossil history of the Subfamily Striginae. Buteo 1998, 10:5-22.

85. Drummond AJ, Nicholls GK, Rodrigo AG, Solomon W: Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data. Genetics 2002, 161:1307-1320.

86. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biology 2006, 4:e88.

87. Drummond AJ, Rambaut A: BEAST v1.4.6. [http://beast.bio.ed.ac.uk/]. Last accessed 4 October 2007.

88. Ho SYW: Calibrating molecular estimates of substitution rates and divergence times in birds. J Avian Biol 2007, 38:409-414.