Translator

 

Calibrando o relógio genômico de aves modernas usando fósseis

Santiago Claramunt https://orcid.org/0000-0002-8926-5974 s.claramunt@utoronto.ca , Edward L. Braun https://orcid.org/0000-0003-1643-5212 , Joel Cracraft https:/ /orcid.org/0000-0001-7587-8342 , +4 e Josefin Stiller https://orcid.org /0000-0001-6009-9581 Informações e afiliações dos autores

16 de setembro de 2024

121 ( 39 ) e2405887121

https://doi.org/10.1073/pnas.2405887121

2.368

Métricas

Total de visualizações 2.368

• Últimos 12 meses 2.368

• 
  

• Vol. 121 | não 39

  • Agradecimentos
  • Referências
• • Informações e Autores
  • Métricas e citações
  • Ver opções
  • Referências
  • Mídia
  • Compartilhar

Wu et al. ( 1 ) conduziram uma análise filogenômica resultando em um novo período de tempo para a diversificação das aves modernas, concluindo que a rápida radiação de Neoaves ocorreu bem antes do evento de extinção do Cretáceo-Paleógeno (K-Pg) e que este evento catastrófico não afetou sua dinâmica de diversificação. Aqui mostramos que os tempos de divergência obtidos por Wu et al. foram comprometidos por problemas com suas escolhas de fósseis e estratégia de calibração. 

 

Das 20 restrições de idade baseadas em fósseis utilizadas por Wu et al., 11 são problemáticas. A única restrição máxima no clado das aves, aplicada à idade da coroa de Neornithes, baseou-se no fóssil mais antigo conhecido de Ichthyornis dispar . Um fóssil pode fornecer um limite máximo para a idade de um nó se for um ancestral direto do nó, mas o fóssil mais antigo de um grupo irmão restringe a idade mínima do caule, não a idade máxima do clado da coroa ( Fig. 1 A ). ( 2 , 3 ). Como Ichthyornis não é um ancestral direto nem um grupo irmão de Neornithes ( 4 ), seu fóssil mais antigo não pode informar diretamente os limites de idade para as aves modernas. 

 

 

 

 

 

Figura 1.

Árvores do tempo de pássaros modernos a partir de calibrações fósseis. ( A ) Ichthyornis restringe a idade mínima possível de Ornithurae porque o clado não pode ser mais jovem que seu fóssil mais antigo ( Topo ), mas, não sendo um ancestral direto, Ichthyornis não restringe a idade máxima possível de Neornithes, que pode ser mais antigo que Ichthyornis ( Fundo ). ( B ) Topo : Árvore temporal obtida por análise de máxima verossimilhança com base nos limites de calibração (colchetes vermelhos) de Wu et al. ( 1 ) (uma calibração alternativa usando uma restrição de idade máxima de 130 Ma para Neornithes com base em seus resultados é mostrada em azul claro). Parte inferior : Árvore de tempo baseada na mesma árvore de máxima verossimilhança, mas informações de calibração da ref. 5 , utilizando as medianas das densidades de calibração como idades fixas. A árvore de máxima verossimilhança inicial foi gerada usando IQ-TREE ( 6 ) para analisar os 1000 loci de codificação semelhantes a relógio de Wu et al. ( 1 ), utilizando o modelo GTR + gama + sítios invariáveis ​​particionados por posição de códon. As árvores de tempo foram obtidas a partir da árvore de máxima verossimilhança por reescalonamento da árvore de tempo de máxima verossimilhança usando um relógio molecular relaxado com cinco categorias de taxa discreta na função cronos em R ( 7 ). A linha amarela representa a fronteira Cretáceo-Paleógeno (K – Pg). Ilustração do crânio de Ichthyornis modificado de OC Marsh 1886, Domínio Público, via Wikimedia Commons. Alinhamentos de sequência, código, análises alternativas e árvores resultantes estão disponíveis em https://doi.org/10.5281/zenodo.11074217 .

 

 Dez restrições de idade mínima também foram problemáticas. Em oito, Wu et al. ignoraram fósseis bem conhecidos, mais antigos e de maior qualidade, de acordo com as melhores práticas atuais ( 8 ) ( Tabela 1 ). Para três deles, eles ignoraram os fósseis do clado irmão, que restringem a idade do caule tanto quanto os fósseis do clado focal ( 3 ). Como resultado, a subestimação das idades mínimas variou de -1,7 a 26,3 milhões de anos entre as calibrações ( Tabela 1 ). No caso dos Aequornithes, os tempos de divergência estimados foram inferiores às idades mínimas indicadas por dois fósseis de alta qualidade e datados com precisão ( Tabela 1 ). 

 

Tabela 1.

Restrições problemáticas de calibração de idade mínima em Wu et al. ( 1 )

	Mínimo (Ma)
Nó de calibração	Wu et al.	correto	Problema
Idade mínima incorreta:
Caule Corvidae	7.2	17.2	Baseado em Corvus larteri [sic] (= Miocorvus larteti Milne-Edwards, 1871) do Mioceno médio de Sansan, França, idade MN 6, portanto, mínimo de 13,7 Ma. Além disso, como nenhum outro membro da infraordem Corvides foi incluído, o nó “tronco Corvidae” na árvore é o ancestral comum mais recente das infraordens Corvides e Passerides, para as quais existem numerosos fósseis mais antigos, incluindo Kurrartapu johnnguyeni Nguyen, 2013 (Corvides : Artamidae, 23,0 a 14,8 Ma) e Certhiops rummeli Manegold, 2008 (Passerides: Certhioidea, 19,5 a 17,2 Ma) ( 9 ).
Suboscines divididos - oscines	13.6	30.0	Baseado em Miocitta galbreathi Brodkorb, 1972, fragmento distal do úmero referido a Corvidae. Mas, para o mesmo nó, existem numerosos fósseis de alta qualidade (esqueletos quase completos) que têm mais do dobro da idade, como o suboscino Wieslochia weissi Mayr e Manegold, 2006 de Frauenweiler (32 a 30 Ma), Alemanha ( 4 , 5 , 9 ).
Caule Psittaciformes	53.5	51.8	Baseado em Pulchrapollia gracilis Dyke e Cooper, 2000, que é mais provavelmente um membro do caule de Psittacopasseria, portanto não restringindo o caule de Psittaciformes, mas em vez disso pode ser restringido pelo caule-Passeriformes Eozygodactylus americanus Weidig, 2010, da Formação Green River (51,8 Mãe), EUA ( 9 ).
Caule Coraciidae + Brachypteraciidae	51.6	54.6	Baseado em Primobucco mcgrewi Brodkorb, 1970, mas Septencoracias morsensis Bourdon et al., 2016 da Formação Fur (54,6 Ma), Dinamarca, é mais antigo ( 4 , 9 ).
Coliiformes-tronco	56.2	62.2	Baseado em Sandcoleus copiosus Houde e Olson, 1992, mas Tsidiiyazhi abini Ksepka et al., 2017 da Formação Nacimiento (62,2 Ma), EUA, é mais antigo ( 4 , 9 ).
Caule Gruoidea	28.3	54.6	Baseado em Parvigrus pohli Mayr, 2005, mas a idade do caule de Gruoidea é a idade da coroa dos Gruiformes, para os quais existem fósseis muito mais antigos, como Pellornis mikkelseni Bertelli et al., 2011 da Formação Fur (54,6 Ma), Dinamarca ( 4 , 9 ).
Caule Apodiformes	51.6	54.6	Baseado em Eocypselus rowei Ksepka et al., 2013, mas Eocypselus vincenti Harrison, 1984 da Formação Fur (54,6 Ma), Dinamarca, é mais antigo ( 4 , 9 ).
Caule Galliformes	51.6	65.7	Baseado em Gallinuloides wyomingensis Eastman, 1900, mas a idade do caule de Galliformes é a idade da coroa de Galloanseres, para a qual existem vários fósseis mais antigos, como o anseriforme Conflicto antarcticus Tambussi et al., 2019 da Formação López de Bertodano (66,0 a 65,7 Ma), Antártica, estendendo seu limite mínimo pelo menos até o limite K – Pg ( 9 ).
Idade estimada mais jovem que o fóssil mais antigo:
Fragatas-tronco	40.6	51.6	A divergência estimada entre Fregatidae e Phalacrocoracidae era mais jovem do que um fóssil de alta qualidade neste clado: Limnofregata azygosternon Olson, 1977, da Formação Green River, EUA ( 4 , 9 ).
Caule Esfenisciformes	45.7	60.5	A divergência estimada entre Procellariiformes e Sphenisciformes era mais jovem do que o fóssil mais antigo deste clado: Waimanu Manningi Slack et al., 2006, da Formação Waipara Greensand, Nova Zelândia ( 4 , 9 ).

 

 

 

Apesar das idades mínimas serem muito jovens, a idade dos Neornithes estimada por Wu et al. (130 Ma) era muito mais antigo do que o limite máximo assumido de 94,3 Ma. Nossa análise de máxima verossimilhança de sequências de codificação semelhantes a relógio de Wu et al. resultou em datas igualmente antigas ao usar suas restrições de calibração mínimo-máximo ( Fig. 1 B ). Em contraste, quando usamos informações de densidades de calibração derivadas do registro fóssil ( 5 , 10 ), encontramos idades mais jovens apoiando uma rápida radiação de Neoaves perto do limite K – Pg ( Fig. 1 C ). 

 

Obtivemos resultados semelhantes com diferentes tipos de dados genômicos e esquemas de particionamento (análises adicionais disponíveis em https://doi.org/10.5281/zenodo.11074217 ) e usando uma análise Bayesiana mais completa de um conjunto de dados filogenômicos maior ( 9 ). Portanto, além dos problemas com os fósseis, as idades antigas encontradas por Wu et al. pode ser o resultado da utilização de limites mínimos e máximos mal definidos que não representam adequadamente a informação de calibração que pode ser derivada dos fósseis.

Portanto, a conclusão de Wu et al. que a rápida diversificação de Neoaves ocorreu no Cretáceo Superior médio, com o evento de extinção K – Pg tendo pouca influência, não tem fundamento. Quando as informações do registro fóssil são usadas de forma mais completa, uma rápida radiação das aves modernas é evidente em torno da fronteira K – Pg ( 5 , 9 ).

Agradecimentos

Contribuições do autor

SC, ELB, JC, JF, SYWH, PH, JMTN e JS projetaram pesquisas; SC, ELB, JC, JF, SYWH, PH, JMTN e JS realizaram pesquisas; SC, ELB, JC, JF, SYWH, PH, JMTN e JS analisaram dados; e SC e SH escreveram o artigo.

Interesses conflitantes

Os autores declaram não haver interesses conflitantes.

Referências

1

S. Wu et al., Genomas, fósseis e o surgimento simultâneo de pássaros modernos e plantas com flores no Cretáceo Superior. Processo. Nacional. Acad. Ciência. EUA 121 , e2319696121 (2024).

Referência cruzada

PubMed

Google Acadêmico

2

SB Hedges, S. Kumar, M. van Tuinen, Restringindo calibrações fósseis para relógios moleculares. Bioensaios 28 , 770–771 (2006).

Ir para referência

Referência cruzada

PubMed

Google Acadêmico

3

SYW Ho, MJ Phillips, Contabilização da incerteza de calibração na estimativa filogenética de tempos de divergência evolutiva. Sist. Biol. 58 , 367–380 (2009).

Referência cruzada

PubMed

Google Acadêmico

4

G. Mayr, Evolução Aviária: O Registro Fóssil de Aves e Seu Significado Paleobiológico (Wiley, 2017).

Google Acadêmico

5

S. Claramunt, J. Cracraft, Uma nova árvore do tempo revela a marca da história da Terra na evolução dos pássaros modernos. Ciência. Av. 1 , e1501005 (2015).

Referência cruzada

PubMed

Google Acadêmico

6

MQ Bui et al., IQ-TREE 2: Novos modelos e métodos eficientes para inferência filogenética na era genômica. Mol. Biol. Evol. 37 , 1530–1534 (2020).

Ir para referência

Referência cruzada

PubMed

Google Acadêmico

7

E. Paradis, Datação molecular de filogenias por métodos de verossimilhança: Uma comparação de modelos e um novo critério de informação. Mol. Filogeneta. Evol. 67 , 436–444 (2013).

Ir para referência

Referência cruzada

PubMed

Google Acadêmico

8

JF Parham et al., Melhores práticas para justificar calibrações fósseis. Sist. Biol. 61 , 346–359 (2012).

Ir para referência

Referência cruzada

PubMed

Google Acadêmico

9

J. Stiller et al., Complexidade da evolução aviária revelada por genomas de nível familiar. Natureza 629 , 851–860 (2024).

Referência cruzada

PubMed

Google Acadêmico

10

S. Claramunt, CladeDate: Gerador de informações de calibração para estimativa de tempo de divergência. Métodos Eco. Evol. 13 , 2331–2338 (2022).

Ir para referência

Referência cruzada

Google Acadêmico

 

 

 

 

Phylogenomics of white-eyes, a ‘great speciator’, reveals Indonesian archipelago as the center of lineage diversity 

Chyi Yin Gwee,

1. Kritika M Garg,
2. Balaji Chattopadhyay,
3. Keren R Sadanandan,
4. Dewi M Prawiradilaga,
5. Martin Irestedt,
6. Fumin Lei,
7. Luke M Bloch,
8. Jessica GH Lee

Research Article Dec 22, 2020

• Cited 0
• Views 220
• Annotations

• 
  

Cite

as: eLife 2020;9:e62765 doi: 10.7554/eLife.62765

• Article
• Figures and data
• Side by side
• 
  
• Abstract
• Introduction
• Results
• Discussion
• Materials and methods
• References
• Decision letter
• Author response
• Article and author information
• Metrics

Abstract

Os arquipélagos servem como importantes "laboratórios naturais" que facilitam o estudo das radiações das ilhas e contribuem para a compreensão dos processos evolutivos. O gênero Zosterops de olhos brancos é um exemplo clássico de um "grande especiador", compreendendo c. 100 espécies de todo o Velho Mundo, a maioria delas insulares. Obtivemos uma extensa amostragem geográfica de DNA de Zosterops usando espécimes históricos e amostras coletadas recentemente. Usando mais de 700 loci de todo o genoma em conjunto com métodos de árvore de espécies coalescentes e abordagens de detecção de fluxo gênico, desvendamos a história evolutiva reticulada de Zosterops, que compreende três clados principais centrados na Indo-África, Ásia e Australásia, respectivamente. A introgressão genética entre as espécies permeia a filogenia de Zosterops, independentemente de quão distantemente relacionadas as espécies são. Crucialmente, identificamos o arquipélago indonésio, e especificamente Bornéu, como o principal centro de diversidade e a única área onde todos os três clados principais se sobrepõem, atestando a importância evolucionária desta região.

Introduction

Os arquipélagos são cenários para desvendar padrões evolutivos complexos, pois constituem laboratórios naturais para o estudo dos fatores que contribuem para a especiação, permitindo um exame da evolução das linhagens isoladas (MacArthur e Wilson, 2001; Whittaker e Fernández-Palacios, 2007; Lohman et al., 2011). Entre os grupos de vertebrados que ocorrem em arquipélagos, as radiações insulares de pássaros são as mais bem estudadas (Lerner et al., 2011; Lamichhaney et al., 2015). Esses modelos aviários exibem uma grande variabilidade em suas taxas de diversificação entre as ilhas, que estão fundamentalmente ligadas à capacidade das espécies de se dispersarem na água (Diamond et al., 1976). Em particular, os chamados 'grandes especiadores', inicialmente caracterizados por Diamond et al., 1976, destacam-se de todas as outras aves pela capacidade paradoxal de se dispersar amplamente e colonizar arquipélagos inteiros enquanto, ao mesmo tempo, diversificam em múltiplos. espécies filhas, apesar de um potencial contínuo para o fluxo gênico sobre a água (Cai et al., 2020).

Um dos poucos exemplos clássicos de 'grandes especiadores' identificados por Diamond et al., 1976 é o gênero de ave canora Zosterops, ou olhos brancos, que são dispersores capazes de se diferenciar rapidamente das populações de origem (Clegg et al., 2002; Moyle et al., 2009). O gênero Zosterops compreende c. 100 espécies que se irradiaram pelo Velho Mundo e Oceania nos últimos 1-3,5 milhões de anos, refletindo uma das taxas de diversificação mais rápidas entre os vertebrados (Warren et al., 2006; Moyle et al., 2009; Jetz et al., 2012 ; Leroy et al., 2019; Cai et al., 2020). Uma proporção esmagadora (mais de 70%) dessas espécies ocorre exclusivamente em arquipélagos distribuídos nos oceanos Atlântico, Índico e Pacífico (Figura 1). Como tal, a radiação dos olhos brancos serve como um sistema modelo para explorar a teoria da biogeografia de ilhas (Diamond et al., 1976; Moyle et al., 2009).

Figure 1

In order to achieve an understanding of the underlying processes driving the white-eye radiation, its phylogeny first needs to be resolved to provide a reliable backbone for hypothesis testing. In Zosterops, however, traditional methods that rely on morphological tools to infer how species are related to one another have proven to be unreliable, as plumage features of ecologically distinct and geographically disjunct Zosterops species are often indistinguishable (Mees, 1957; Mayr, 1965). Although a more recent application of genetic methods has helped disentangle the white-eye radiation to some extent, most studies have concentrated on Afrotropical, Melanesian, and Indian Ocean members of the genus (Slikas et al., 2000; Warren et al., 2006; Moyle et al., 2009; Cox et al., 2014; Linck et al., 2016; Wickramasinghe et al., 2017; Manthey et al., 2020; Martins et al., 2020). There continues to be a dearth of knowledge on this radiation across the core of its Asian distribution due to limited sampling and lack of genetic data. In particular, it is crucial to unravel the phylogenetic affinities of white-eyes distributed across the Indonesian archipelago, which – comprising more than 17,000 islands – is the largest archipelago in the world and harbors about 20 endemic Zosterops species (Figure 1), including two species that were discovered in the last two decades and remain undescribed (Eaton et al., 2016; O’Connell et al., 2019). The high density of Zosterops species across the Indonesian archipelago hints at the possible importance of this region in white-eye evolution.

Apart from incomplete geographic sampling, the lack of resolution of the white-eye radiation has largely been a consequence of sparse genomic sampling: most phylogenetic studies of white-eyes have been restricted to one or a few genetic markers, resulting in trees that are plagued by unresolved polytomies, hampering useful evolutionary inference (Slikas et al., 2000; Warren et al., 2006; Moyle et al., 2009; Oatley et al., 2012; Á.S and Joseph, 2013; Cox et al., 2014; Husemann et al., 2016; Linck et al., 2016; Round et al., 2017; Wickramasinghe et al., 2017; Shakya et al., 2018; Cai et al., 2019; Lim et al., 2019; O’Connell et al., 2019; Martins et al., 2020). Disentangling relationships within rapid and recent radiations, such as white-eyes, requires overcoming the challenges of heterogeneous gene trees due to biological factors such as incomplete lineage sorting (Edwards et al., 2005; Song et al., 2012). The multispecies coalescent (MSC) model offers a promising avenue to overcoming gene tree discordance by allowing the evolutionary histories of each locus to be inferred independently (Song et al., 2012; Liu et al., 2015).

An additional source of bias in reconstructing the phylogeny of rapid and recent radiations such as white-eyes is the potential for member species to engage in secondary gene flow, defined as post-speciation gene flow, or genetic introgression long after a speciation event has occurred (Rheindt and Edwards, 2011; Edwards et al., 2016). Such introgression will be reflected in the phylogenetic signal of a varying proportion of loci, thereby leading astray efforts to search for the true species tree. Multiple analytical approaches have been devised to account for secondary gene flow, such as through tree-based analysis as implemented in Phylogeographic Inference using Approximate Likelihoods (PHRAPL) (Jackson et al., 2017a; Jackson et al., 2017b), or through SNP-based analysis like the ABBA-BABA test, which detects an excess of shared derived alleles between populations (Green et al., 2010; Patterson et al., 2012).

In this study, we used historical specimens and recently collected samples to represent 33 white-eye species across the Southern hemisphere, especially from the understudied Indonesian archipelago (Supplementary file 1). We designed RNA probes using Z. lateralis (Cornetti et al., 2015) as a reference genome, targeting 832 loci at high coverage to overcome the limitations of missing data expected from degraded DNA of historical samples, thereby recovering a comparable set of loci across both historical and fresh samples (Templeton et al., 2013). Our target capture methodology addresses the recalcitrant persistence of unresolved polytomies in the phylogeny of this rapid radiation by making use of a large set of loci (Cai et al., 2019). We employed three different coalescent species tree methods to assess topological incongruence across tree-building approaches (Liu et al., 2009; Liu et al., 2010; Vachaspati and Warnow, 2015). Recognizing that gene flow is commonly observed in recent radiations, we conducted PHRAPL (Jackson et al., 2017b) analysis and performed ABBA-BABA tests (Patterson et al., 2012) to assess introgression between closely related species with incongruent topologies and ultimately elucidate the likely evolutionary history of this complex radiation.

Results

Congruent phylogenetic trees reveal three distinct lineages

To shed light on the phylogenetic relationships of Zosterops species, we employed both concatenation methods, in which sequence data from individual loci are combined into one larger sequence, as well as MSC approaches (MP-EST [Liu et al., 2010], STAR [Liu et al., 2009], and ASTRID [Vachaspati and Warnow, 2015]), which account for individual gene tree stochasticity in a coalescent framework (Edwards et al., 2007). All four tree inference methods produced a similar phylogeny with a congruent tree topology for highly supported nodes, except for the placement of a Sundaic group consisting of Z. atricapilla and Z. auriventer (Figure 2). These two Sundaic taxa are embedded within the Australasian clade in the concatenated tree, but emerged within the Asiatic clade in the species trees constructed with MP-EST (henceforth our baseline species tree) and ASTRID (Figure 2). The other species tree method, STAR, shows a weak bootstrap support for an unresolved placement of these Sundaic taxa.

Figure 2 with 3 supplements

All methods reveal a white-eye radiation divided into three main clades consisting of an Indo-African, Asiatic, and Australasian group (Figure 2). Our taxon sampling covered 33 out of 108 white-eye species across the global radiation (Supplementary file 2). In order to expand clade assignment to well-studied white-eye species outside of our sampling regime, we examined an additional 30 species shown to have high bootstrap support for placement within any one of the three main clades based on previously published papers up until 2019 (Warren et al., 2006; Moyle et al., 2009; Cox et al., 2014; Cornetti et al., 2015; Shakya et al., 2018; Cai et al., 2019; O’Connell et al., 2019; Figure 2—figure supplement 2). These additional clade assignments were not used in the construction of our phylogenetic trees, but directly examined from the trees constructed by various studies. For example, Cai et al., 2019 provide high bootstrap support (>90%) for the position of Z. mouroniensis as a descendant of the most recent common ancestor of two unequivocal members of the Indo-African clade (i.e., Z. palpebrosus and Z. borbonicus; Figure 2); thus the species is allocated accordingly, and the breeding distribution of Z. mouroniensis (Mt Karthala on Grande Comore Island) is shaded yellow on the global range map (Figure 2a).

A majority of the added taxa are distributed across Africa, where all Zosterops species fall within the Indo-African clade, while five of these newly added taxa are distributed across Melanesia, where our species coverage allowed us to detect the presence of only the Australasian clade. Therefore, our mapping suggests that the Afrotropical and the Australo-Papuan regions are depauperate in deeper Zosterops lineage diversity (Figure 2a; Supplementary file 2). Similarly, our results reveal that most areas in continental Asia generally harbor only one of the three main Zosterops clades, except East Asia where two clades co-occur in a narrow zone of overlap between the Indo-African Z. palpebrosus and the Asiatic Z. simplex (orange in Figure 2a). In contrast, the Indonesian archipelago emerged as a center of modern-day diversity for Zosterops, with all three main clades represented on Java and Borneo, and two main clades on many other islands (Figure 2a).

Presence of secondary gene flow

The genetic signal of recent and rapid radiations is often convoluted by the presence of secondary gene flow, leading to heterogeneous gene trees which deviate from the true phylogeny. We assessed the presence of secondary gene flow specifically between members of a Sundaic species pair (Z. auriventer and Z. atricapilla) characterized by a shifting and incongruent placement across trees with representatives from the Asiatic (i.e. Z. simplex) and Australasian (i.e. Z. emiliae and Z. melanurus) clades (Figure 3a). The top two demographic models inferred by PHRAPL simulations show that Sundaic Z. auriventer is more closely related to Asiatic Z. simplex than to Australasian Z. emiliae, but inconclusive in relation to Z. melanurus due to ancestral and/or secondary gene flow between all three taxa (Figure 3b). Additionally, ancestral and/or secondary gene flow was detected between all taxa in both comparisons involving Z. atricapilla, Z. simplex, and either Z. emiliae or Z. melanurus (Figure 3—figure supplement 1). The genealogical divergence index (gdi) of the top demographic model inferred for each combination is relatively high, ranging between 0.321 and 0.751, suggesting high divergence between sister species despite secondary gene flow.

Figure 3 with 1 supplement

The ABBA-BABA approach corroborates that the conflicting placement of the Sundaic species pair consisting of Z. atricapilla and Z. auriventer may be attributed to secondary gene flow (Figure 3). In both concatenated and species trees, Z. emiliae emerges as basal to all other members of the Australasian clade (Figures 2 and 3a). Therefore, in the absence of introgression, the Sundaic Z. atricapilla and Z. auriventer should be equally closely related to both Z. emiliae and other members of the Australasian clade (tree topology in Figure 3a). The ABBA-BABA statistics revealed a significant excess of allele sharing between the Sundaic taxa (Z. atricapilla and Z. auriventer) and two Australasian taxa (Z. melanurus and Z. citrinella) (Figure 3c, Table 1). This result reflects that ancestral introgression occurred between the Sundaic taxa and the Australasian lineage after Z. emiliae had diverged (Figure 3c).

Table 1

D-statistics of a selection of species combinations to test if two species (H2, H3) exhibit an excess of derived allele sharing.

The ABBA-BABA test was restricted to species combinations with conflicting tree topologies observed in this study (set 1), conflicting tree topologies between this study and Cai et al., 2019 (set …

Set	H1	H2	H3	D-stat	Z	No. of ABBA	No. of BABA
1	emiliae	melanurus	atricapilla	0.129	3.807	192.84	148.63
	emiliae	melanurus	auriventer	0.163	5.168	190.17	136.97
	emiliae	citrinella	atricapilla	0.110	3.086	178.55	143.09
	emiliae	citrinella	auriventer	0.125	3.788	172.21	133.99
2	citrinella	melanurus	simplex	0.093	3.074	142.55	118.22
3	simplex	japonicus	erythropleurus	0.215	5.272	171.87	111.15
4	erythropleurus	simplex	auriventer	0.108	3.169	166.20	133.74
5	melanurus	emiliae	simplex	−0.053	−1.456	152.33	169.34
	melanurus	emiliae	japonicus	−0.066	−1.729	149.05	170.02
	citrinella	emiliae	simplex	0.027	0.706	154.77	146.74
	citrinella	emiliae	japonicus	−0.027	−0.686	146.72	154.82

In addition, secondary gene flow was also detected in several pairs of species which overlap in present-day distribution (sets 2 to 4 in Table 1). For instance, we found that Z. melanurus shares significantly more alleles with Z. simplex than with Z. erythropleurus, suggesting occasional hybridization between Z. melanurus buxtoni and Z. simplex erwini on Sumatra, where both occur. In a similar vein, Z. erythropleurus and Z. japonicus displayed excess allele sharing, indicating potential introgression in areas of overlap in Korea, and Z. auriventer and Z. simplex exhibited excess allele sharing that hints at occasional hybridization in parts of peninsular Malaysia and Borneo where they overlap.

Discussion

Phylogeny of Zosterops and presence of secondary gene flow

The evolutionary history of Zosterops has received a fair amount of scientific attention, but mostly by means of single mitochondrial or few nuclear loci, therefore resulting in trees plagued by unresolved polytomies (e.g. Figure 2—figure supplement 3; Degnan and Moritz, 1992; Degnan, 1993; Slikas et al., 2000; Warren et al., 2006; Moyle et al., 2009; Oatley et al., 2012; Á.S and Joseph, 2013; Cox et al., 2014; Husemann et al., 2016; Linck et al., 2016; Round et al., 2017; Wickramasinghe et al., 2017; Shakya et al., 2018; Cai et al., 2019; Lim et al., 2019; O’Connell et al., 2019). Using more than 700 genome-wide loci with a dense species sampling, our study produced an improved phylogeny of Zosterops and reveals the existence of three discrete main clades characterized by an Indo-African, Asiatic, and Australasian core of distribution, respectively (Figure 2). However, there was no strong support to unravel the sequence of diversification events among those three main clades. This lack of basal resolution could be attributed to a quick succession of divergence events at the time, and/or reticulated evolutionary history unresolved by species tree and concatenated tree approaches.

Both concatenated and species tree analyses resulted in congruent topologies for well-supported nodes, except for the placement of a Sundaic species pair comprising Z. atricapilla and Z. auriventer (Figure 2). Uncertainties in the affinity of these two Sundaic species appear to be due to secondary gene flow between the incipient Asiatic and Australasian clades (Figure 3). Such genetic introgression destabilizes the robustness of the MSC model, which does not account for secondary gene flow, thereby confounding species tree estimation (Figure 3, Table 1). At the same time, traditional phylogenetic approaches such as concatenation are equally negatively impacted by secondary gene flow and are additionally subject to the biases of incomplete lineage sorting (Liu et al., 2019).

To assess the magnitude of genetic introgression that has resulted in the controversial placement of the two Sundaic taxa, we used analytical approaches that specifically account for post-speciation gene flow, such as PHRAPL. Most of the top demographic models inferred by PHRAPL produced topologies in which the Sundaic species of controversial placement emerge as sister to the Asiatic clade (Figure 3b, Figure 3—figure supplement 1). However, conflicting demographic models were observed between the top two results of some combinations due to ancestral and/or secondary gene flow in all tested directions (Figure 3, Figure 3—figure supplement 1). For instance, the erroneous inference of Z. melanurus and Z. simplex emerging as sister species to each other – basal to Z. auriventer – may be attributed to secondary gene flow also present between Z. melanurus and Z. simplex (Figure 3b, Table 1). Convergence onto a single demographic model may require further PHRAPL simulations allowing for more parameters such as asymmetrical rates of gene flow (Morales et al., 2017). Additional demographic analyses using other programs, such as DaDi (Gutenkunst et al., 2010) and fastsimcoal (Excoffier and Foll, 2011), may also assist in fully disentangling the complex relationships within this rapid radiation.

SNP-based analysis using ABBA-BABA statistics conclusively identified introgression as an underlying cause of the conflicting placement of Z. atricapilla and Z. auriventer. Specifically, excess allele sharing between these two species and Z. melanurus, but not between them and Z. emiliae, suggests introgression between the incipient stages of the Asiatic and Australasian clades – after Z. emiliae had split off from other Australasian species (Figure 3c). Such ancient introgression generates patterns of allele sharing that would lead to the two controversial Sundaic species partly being reflected as members of the one clade or of the other (Figures 2 and 3), depending on tree inference methods.

More generally, the ABBA-BABA test detected rampant secondary gene flow between species that geographically overlap, regardless of their phylogenetic proximity, underscoring the pervasive nature of genetic introgression in rapidly evolving lineages such as white-eyes (black arrows in Figure 2a; sets 2–4 in Table 1). For example, Z. auriventer and Z. simplex display excess allele sharing indicating recent gene flow in areas of Sundaland where they overlap, even though they belong to different main clades of Zosterops (Figure 2, Table 1). The same is true for Z. melanurus and Z. simplex, which co-occur on Sumatra, and for Z. japonicus and Z. erythropleurus, which overlap in Korea (Table 1). Such introgression between species is likely recent, limited, and of the kind that the MSC model remains robust to Liu et al., 2009. The detection of rampant secondary gene flow in multiple pairs of sympatric white-eye species is in agreement with the recent discovery of introgression between various non-sister white-eye species across the Solomon Islands (Manthey et al., 2020).

Indonesian archipelago harbors all three main clades

While the highest rates of Zosterops diversification are known to have occurred on archipelagoes in general (Diamond et al., 1976; Moyle et al., 2009), the geographic distribution of deeper-level lineage diversity in this genus remains unexplored. Our phylogenetic results demonstrate that Africa, most parts of continental Asia and probably also all of the Australo-Papuan region each harbor representatives from only one of the three main Zosterops clades, respectively (Figure 2a), regardless of Zosterops species diversity. In contrast, virtually all parts of the Indonesian archipelago harbor white-eye species from two to three of the main Zosterops clades (Figure 2a). Positioned between the Sunda and Sahul shelf, the Indonesian islands are a center of syntopy of lineages from either side of Wallace’s Line (Moss and Wilson, 1998; de Bruyn et al., 2014). This is consistent with the rapid rate of tectonic change reconstructed for the Indonesian archipelago over the last 30 million years (Hall, 2002; Hall, 2012; de Bruyn et al., 2014; Nugraha and Hall, 2018), which has led to a narrowing of open sea between Asia and Australia and thereby facilitated overwater dispersal of Sundaic and/or Australo-Papuan lineages for many organismic groups (e.g. Heads, 2001; Irestedt et al., 2013; Gwee et al., 2017; Ng et al., 2017; Garg et al., 2018; Reilly et al., 2019; Oliver et al., 2020, including white-eyes). At the same time, while the narrowing of the sea gap between Australo-Papua and Asia has created numerous stepping-stone islands to facilitate overwater dispersal, most of the Wallacean region remains dominated by deep sea, and there are as yet no land bridges (Voris, 2000; Rheindt et al., 2020). This complicated archipelagic setting has likely acted as a diversification driver in white-eyes, which have the capability of colonizing and populating these deep-sea islands.

To the west of Wallace’s Line, the Sundaic islands of Borneo and Java, which constitute a large part of the Greater Sunda archipelago, each harbor a number of representatives of the Asiatic and Australasian clades but are additionally inhabited by the coastal endemic species Z. flavus. This species is phylogenetically more closely related with the Indo-African clade rather than with the geographically more proximate Asiatic and Australasian clades, demonstrating an impressive potential for dispersal capability that may in part underlie the rapid diversification rate of the genus.

Low-lying Indonesian islands to the east of Wallace’s Line, such as Sumba and Kai, generally harbor only Zosterops species from the Australasian clade, whereas Z. japonicus of the Asiatic clade additionally occurs on islands that reach montane elevations of over 1200 m. For example, Z. japonicus is present on mountainous Buru (c. 9500 km² in size) but absent on Sumba (c. 11,000 km² in size), which largely comprises savannah with a small hilly region not exceeding 1200 m in elevation, despite Sumba being a larger island (Figure 1). Although the elevation of an island contributes substantially to the number of main Zosterops clades present, it seems to have less influence on the total number of Zosterops species (Figure 1). For instance, the Kai islands, with a combined area of only approximately 1400 km² and an elevation of less than 700 m, harbor three distinct Zosterops species, including two island endemics Z. uropygialis and Z. grayi (not sampled) (Figures 1 and 2).

Our phylogenetic results reveal that the widespread Lemon-bellied White-eye Z. chloris is non-monophyletic (Figure 2; Figure 2—figure supplement 1) and several small, low-lying islands situated between the Banda Sea and Arafura Sea, including Kai and Aru, may harbor a cryptic species morphologically identical to the Lemon-bellied White-eye. White-eyes are renowned for their conservative morphology, which contributes to the confusion that has surrounded their taxonomic treatment (Mees, 1957; Mayr, 1965; Lim et al., 2019; Manthey et al., 2020). Further research is required to ascertain the evolutionary status of these and other overlooked island populations.

Borneo is a hotspot for evolution and harbors deep phylogenetic isolates

Borneo has been identified as a major source of diversification across the Southeast Asian region for a variety of organismic groups, including birds, mammals, amphibians, and plants (de Bruyn et al., 2014). Our study reveals Borneo’s unique status as the only place in the distribution of the genus where members of all three main Zosterops clades occur, and where as many as four Zosterops species co-exist within a few square kilometers of one another, rendering it the center of faunal mixing for white-eyes (Figures 1 and 2). Borneo forms the eastern part of the Sundaic region, which – at present – is splintered into multiple bigger and many smaller landmasses comprising the Greater Sunda Islands and Malay Peninsula. For the longest time over the past 400,000 years, however, these landmasses have been merged into a larger landmass, Sundaland, that has facilitated the evolution of much of Southeast Asia’s equatorial rainforest fauna (Sarr et al., 2019). The east of Sundaland (i.e. Borneo) has constituted a particularly stable part of this subcontinent, remaining above water for the longest uninterrupted time and offering a wide variety of habitats such as mangroves (Z. flavus), submontane forest (Z. atricapilla and Z. auriventer), montane forest (Z. emiliae and a still undescribed white-eye from the Meratus range Eaton et al., 2016), and coastal woodland (Z. simplex). Our results from ancestral range estimation suggest an Asian origin for the entire Zosterops radiation, and a Sundaic origin for the Australasian clade (Figure 2—figure supplement 2). Borneo constitutes the largest landmass within the Sundaic region, and its elevated count in discrete Zosterops lineages suggests that it has played a key role in the diversification of this important radiation.

While Borneo’s exceptional biodiversity has been appreciated as early as during Alfred R. Wallace’s times (Wallace, 1962), most of this diversity has traditionally been interpreted as being of a Sundaic element largely shared with Sumatra and peninsular Malaysia, and is only slowly being recognized as having attained species-level depths of differentiation (Cros et al., 2020). On the other hand, Borneo is known for hosting a number of deep phylogenetic isolates, such as the enigmatic Bristlehead Pityriasis gymnocephala (Oliveros et al., 2019). In the context of Zosterops diversification, we add the Mountain Black-eye Z. emiliae as an overlooked phylogenetic isolate (Figure 2). Although it has been shown to be embedded within Zosterops for over a decade (Moyle et al., 2009; Cai et al., 2019), it continues to be treated as a monospecific genus (Chlorocharis) by some modern sources (del Hoyo et al., 2016). Using our much-improved taxon sampling, Z. emiliae emerged as a basal sister to the Australasian clade with moderate support under various analytical regimes (Figure 2). In contrast, Cai et al., 2019 placed Z. emiliae with members of the Asiatic clade. Our use of the ABBA-BABA test to verify whether this conflicting position may be due to ancient introgression did not identify an excess of allele sharing between Z. emiliae and the Asiatic members (set five in Table 1), suggesting that the incongruent placement of Z. emiliae by Cai et al., 2019 is unlikely to be a result of secondary gene flow. Instead, incomplete lineage sorting may have generated such phylogenetic conflict as the divergence of Z. emiliae likely fell within a time of rapid diversification within the genus, leading to its recalcitrance to phylogenetic resolution when only few loci are applied.

Conclusions

Our study presents the application of species tree methods on a large set of genome-wide markers across a comprehensive sampling of members of a rapid radiation of a classic ‘great speciator’. Using approaches to account for secondary gene flow, our study demonstrates the pervasive presence of genetic introgression across this explosive radiation. The resultant phylogeny of Zosterops white-eyes reveals that the Indonesian archipelago, and Borneo in particular, are an evolutionary hotspot for the diversification of the genus. This archipelagic region harbors members of clades centred in the neighboring Asian and Australo-Papuan landmasses, and even from a geographically distant Indo-African clade. The western Indonesian archipelago is the sunken remnant of a subcontinent – Sundaland – that only started to be periodically submerged starting from ~400,000 years ago, and offers potential for differentiation. The identification of areas in western Indonesia as a major center of modern phylogenetic diversity not only contributes to their recognition as an arena of important evolutionary processes, but also elevates their status as a region of global conservation relevance.

Materials and methods

Taxon sampling

Request a detailed protocol

A total of 48 historical toepad samples and 52 fresh samples were acquired from various museums and through fieldwork conducted across peninsular Malaysia and the Indonesian archipelago (Ashari et al., 2019; Supplementary file 1). In total, 33 white-eye species were represented [following the taxonomy by del Hoyo et al., 2016 with more recent revisions by Lim et al., 2019 and Martins et al., 2020 (see Supplementary file 2)].

Probe design for target capture

Request a detailed protocol

Target enrichment protocols have been shown to be highly effective at capturing historical DNA for phylogenomic studies (Bryson et al., 2016; van der Valk et al., 2017; Chattopadhyay et al., 2019; Baveja et al., 2020). We designed loci specifically targeting both conserved exons and variable intronic regions of the Zosterops genome (Chattopadhyay et al., 2019). We first used EvolMarkers (Li et al., 2012) to identify conserved single copy coding sequences in the genomes of Z. lateralis (accession no. GCA_001281735) (Cornetti et al., 2015), Ficedula albicollis (accession no. GCA_000247815.1) (Ellegren et al., 2012), and Taeniopygia guttata (accession no. GCF_003957565.1; released by the Vertebrate Genomes Project). To identify conserved exons, EvolMarkers performs a BLAST search (Altschul et al., 1990), for which we set a minimum of 55% identity and e-value of less than 10E-15. Only single-hit exons longer than 500 bp were used for further downstream analysis. Then we isolated 500 bp upstream and downstream of these conserved exons from the Z. lateralis genome to include variable intronic regions using bedtools 2.28.0 (Quinlan and Hall, 2010). We further checked for overlapping targets and merged all overlapping loci in bedtools, removing any loci with GC content less than 40% or more than 60%. Loci comprising repeat elements were identified using RepeatMasker 4.0.7 (Smit et al., 2015) and removed. After filtering, our design retained 832 loci, which were used by Arbor Biosciences (USA) to design a total of 63,244 RNA baits. Each locus was targeted with 4X tiling density of overlapping baits, each bait of 100 bp, for in-solution target enrichment.

Laboratory procedures

Request a detailed protocol

Fresh DNA was extracted following the manufacturer’s protocol using the DNeasy Blood and Tissue Kit (Qiagen, Germany). The DNA of historical toepad samples was extracted under sterile conditions inside a dedicated ancient DNA facility, and extractions were performed inside a biosafety cabinet with laminar air-flow. The ancient DNA facility room was subject to at least 12 hr of UV light and thoroughly cleaned with bleach in between each session of historical DNA extractions. We used the same kit for extraction of historical DNA with slight modifications (Chattopadhyay et al., 2019). Extraction negatives were included to ensure absence of contamination. Double-stranded DNA concentrations were ascertained using a Qubit 2.0 high sensitivity DNA Assay kit (Invitrogen, USA), and fragment sizes were assessed using an AATI Fragment Analyzer (Agilent, USA). The negatives were also quantified using a Qubit 2.0 assay and AATI to ensure absence of DNA.

Fresh DNA was sheared into a targeted size of 250 bp using a Bioruptor Pico sonication device (Diagenode, Belgium) with 13 cycles of sonication prior to library preparation. Each cycle consisted of 30 s of sonication followed by 30 s of rest. We used NEBNext Ultra II DNA Library Prep Kits for Illumina (New England BioLabs, USA) and NEBNext 8 bp dual indexes (New England BioLabs, USA) for both fresh and historical library preparation. The libraries using fresh tissue were size selected for an insert size of 250 bp with AMPure XP (Beckman Coulter, USA) beads, giving an expected final library size of ~370 bp with adapters and primers included. Size selection was omitted during library preparation of historical samples to reduce DNA loss. A total of 12 cycles of polymerase chain reaction (PCR) was applied and negative controls were carried out for both fresh and historical libraries.

The library preparation of historical DNA was conducted inside a dedicated PCR cabinet with laminar air-flow. The PCR cabinet was subject to 1 hr of UV light and thoroughly cleaned with bleach in between each batch of library preparations. We added NEBNext FFPE DNA repair mix (New England BioLabs, USA) to the historical DNA prior to library preparation to reduce deamination of cytosine to uracil, repair nicks, and fill in 5’ overhangs of the damaged DNA. DNA quantification and assessments of libraries’ fragment sizes were conducted as above. The peak fragment sizes of each library prepared with historical samples ranged between 200 bp and 300 bp, whereas the peak fragment sizes of each library prepared with fresh samples ranged between 330 bp and 420 bp. The negatives were also quantified by Qubit 2.0 and assessed by AATI to ensure that only adapters and primer-dimer DNA (single peak at ~55 bp and ~150 bp respectively) were present.

Target capture was performed on all historical and fresh samples using a MYbaits kit version 3 (Arbor Biosciences, USA), with a modified protocol following Chattopadhyay et al., 2019. In brief, we diluted the volume of baits to 1.85 µL per historical DNA sample (~3X dilution) and 1.1 µL per fresh DNA sample (~5X dilution). Biotinylated RNA baits and target sequences were hybridized at 60°C for 40 hr for historical samples and at 65°C for 20 hr for fresh samples. Following hybridization, the samples were cleaned according to the myBaits manual, and PCR was conducted using IS5 and IS6 primers with 20 cycles for historical samples and 15 cycles for fresh samples (Fortes and Paijmans, 2015). DNA quantification and assessments of libraries’ fragment sizes were conducted as above (see extraction). The peak fragment sizes of each historical target capture library ranged between 250 bp and 400 bp, whereas the peak fragment sizes of each fresh target capture library ranged between 370 bp and 410 bp. The target capture libraries were sequenced using the Illumina HiSeq 2500 and HiSeq 4000 platforms with 150 bp paired-end runs for all samples. Fresh and historical samples were run on separate lanes.

We additionally sequenced the whole genomes of nine Zosterops individuals. Their DNA was extracted using the DNeasy Blood and Tissue Kit (Qiagen, Germany) according to the manufacturer’s protocol. The samples were then prepared using a Nextera Library Prep Kit (Illumina, USA) with dual indexes. The whole genome libraries were sequenced on the Illumina X10 or NovaSeq platforms at Medgenome (Foster City, California) with 150 bp paired-end runs.

Sequence assembly

Request a detailed protocol

We removed adapter sequences using Trimmomatic 0.38 (Bolger et al., 2014) and duplicates with FastUniq 1.1 (Xu et al., 2012). Paired trimmed reads of historical samples were examined with mapDamage 2.0.9 (Jónsson et al., 2013) to assess DNA deamination in read ends. We further trimmed 5 bp from the 3’ ends of both forward and reverse reads as the mapDamage results show a high amount (>0.1 probability) of G to A misincorporation in the read ends, and reran mapDamage to ensure the average probability of misincorporation across samples remained below 0.1. We ran HybPiper 1.3.1 (Johnson et al., 2016) to extract target sequences. Following the HybPiper pipeline, the trimmed reads were aligned to each target gene using BWA 0.7.17 (Li, 2013). We conducted a de novo assembly of target sequences using SPAdes 3.13 (Bankevich et al., 2012), applying a sequencing depth cut-off of at least 16X coverage per contig. The contigs generated by SPAdes were re-aligned against the target sequences using Exonerate 2.4.0 (Slater and Birney, 2005) to assemble coding sequence regions (including intronic regions) and the resulting DNA sequence of each locus was extracted for downstream analyses. The length of each locus assembled for each sample was examined using the following python scripts in the HybPiper package: get_seq_lengths.py and hybpiper_stats.py. We removed four historical samples due to high amount of missing data and kept 770 loci after visual quality checks across all samples, ensuring each locus is present in at least 85% of individuals and contains less than 30% missing nucleotides.

An additional 12 Zosterops samples were included by extracting target sequences from their whole genomes using blastn in BLAST+ 2.6.0 (Camacho et al., 2009). The whole genomes of Z. lateralis (accession no. GCA_001281735.1) (Cornetti et al., 2015), Z. pallidus (accession no. GCA_007556475.1) (Leroy et al., 2019), and Z. borbonicus (accession no. GCA_007252995.1) (Leroy et al., 2019) were obtained from NCBI. Whole genomes of nine individuals were resequenced for this study (Supplementary file 1). Target sequences were also extracted from the Mixornis gularis whole genome (Tan et al., 2018) to be used as an outgroup for some downstream analyses. For the nine whole genomes that were generated for the present study, we first cleaned the raw reads using a modified Perl script to remove exact PCR duplicates and low complexity reads (Bi et al., 2012; Singhal, 2013), and used Trimmomatic 0.30 (Bolger et al., 2014) for adapter removal. Contaminants were removed by aligning the raw reads to the reference genomes of potential contaminant sources (such as ribosomal RNA, human, and bacterial DNA) with Bowtie 2.0.1 (Langmead and Salzberg, 2012), followed by another cleanup of reads using Cutadapt 1.16 (Martin, 2011). We merged overlapping paired reads using FLASH 1.2.11 (Magoč and Salzberg, 2011), and aligned the cleaned paired reads against the Zosterops lateralis genome using BWA-MEM in BWA 0.7.8 (Li, 2013). The data were then converted into bam file format and sorted using SAMtools 1.5 (Li et al., 2009). We used mpileup in BCFtools (Li et al., 2009) to calculate the genotype likelihoods of each site of the reads, and generated a consensus sequence in fasta format using BCFtools. For all 12 whole genome resequenced samples, we used blastn in BLAST+ 2.6.0 (Camacho et al., 2009) to extract the same set of loci as that used in target capture.

Phylogenetic analyses

Request a detailed protocol

MAFFT 7.0 (Katoh and Standley, 2013) was run to ensure each locus direction was consistent throughout all samples. All 770 loci were concatenated, giving a final alignment of 1,635,155 bp with 7.60% gaps. RAxML 8.2.12 (Stamatakis, 2014) was run to construct a Maximum Likelihood (ML) tree using a GTR + GAMMA + Invariant Sites model with 100 rapid bootstraps to search for the best-scoring ML tree, and the tree was rooted with M. gularis.

Best-scoring ML gene trees were inferred for each locus with RAxML 8.2.12 (Stamatakis, 2014) with 20 independent searches from a random starting tree and using a GTR + GAMMA + Invariant Sites substitution model. Each gene tree was run with 100 bootstraps for node support. A total of 768 gene trees were rooted with M. gularis, while no outgroup sequence was present for two gene trees which had to be discarded for species tree inference.

We adopted the MSC model using three different algorithm methods: MP-EST 2.0 (Liu et al., 2010), STAR (Liu et al., 2009), and ASTRID 1.4 (Vachaspati and Warnow, 2015). All three species tree methods were run with the 768 best-scoring ML gene trees to infer the species tree topology, and with 100 different sets of input gene trees to infer bootstrap support. Nodes with less than 68% bootstrap support were collapsed.

We assessed the level of congruence in the phylogenetic placement of Zosterops species between the present and previously published datasets following the species-level classification by del Hoyo et al., 2016 and the more recent taxonomic revisions within the Asiatic and Afrotropical white-eye complexes by Lim et al., 2019 and Martins et al., 2020, respectively (see Supplementary file 2). Specifically, we assigned species to one of three main clades that emerged within the genus in our analyses. We expanded our clade assignment to 30 Zosterops taxa not sampled in our study but found to be embedded within one of the three main clades with high bootstrap support (>90%) by at least one study up until 2019 (Supplementary file 2). These additional clade assignments were not used in the construction of our phylogenetic trees, but directly examined from the trees constructed by the respective studies indicated in Supplementary file 2.

Mitochondrial tree

Request a detailed protocol

We observed non-specific binding during hybridization of probes and sample DNA, allowing us to assemble mitochondrial DNA from the raw reads of each individual prepared by target capture. The raw reads were first mapped to the reference mitogenome of Z. lateralis (accession no. NC029146) using bwa 0.7.17 (Li, 2013), and converted to bam files using SAMtools 1.9 (Li et al., 2009). The bam files were then imported into CLC Genomics Workbench 7.0.4, remapped to the same reference mitogenome, and locally re-aligned. A consensus mitogenome of each individual was extracted with a minimum coverage of five, otherwise an ambiguous base ‘N’ was inserted. Finally, we extracted 1041 bp of ND2 sequence from each individual by aligning each assembled mitogenome to the ND2 sequence of Z. lateralis. Some samples were removed due to extensive missing nucleotides, and the ND2 sequences of 68 individuals were retained. As a means to assess the presence of artifacts from DNA damage, especially in toepad DNA, we compared the ND2 sequences generated in this study with the ND2 sequences of the same taxa deposited on GenBank by previously published studies. We also added the Genbank sequences of 16 Zosterops species not represented in our sampling regime. A maximum likelihood tree was generated using RAxML with 10,000 bootstrap replicates under the GTR + GAMMA model.

Testing gene flow with ABBA-BABA statistics

Request a detailed protocol

We computed SNP-based ABBA-BABA statistics to test for gene flow among a subset of taxa with topological incongruence among trees, as well as populations which overlap in distribution. The bam files generated during locus assembly were used as input for SNP calling for the fresh samples, while the bam files of the historical samples were first subjected to mapDamage 2.0.9 (Jónsson et al., 2013) to rescale the quality scores of possibly deaminated sites. We used ANGSD 0.923 (Korneliussen et al., 2014) to call SNPs with the following filters applied: minimum depth of 20, block size of 50,000 bp, remove transitions, minimum mapping quality of 30, and minimum base quality of 20. We computed D statistics and used the jackknife.R script in ANGSD 0.923 (Korneliussen et al., 2014) to compute critical values (Z) and test for significance. A positive critical value Z > 3, corresponding to a p-value below 0.0013, suggests a significant excess of ABBA-like alleles as compared to BABA-like alleles (Green et al., 2010; Patterson et al., 2012). This critical value is widely applied as a threshold for detecting introgression and reduces the likelihood of false positives (Zheng and Janke, 2018).

Demographic analysis

Request a detailed protocol

We additionally assessed the presence of secondary gene flow using PHRAPL 0.62 (Jackson et al., 2017b) by simulating the probability of observing a set of gene trees across various demographic models. We conducted simulations on a smaller subset of populations exhibiting topological incongruence: Z. auriventer (A), Z. atricapilla (B), Z. simplex (C), Z. emiliae (D) and Z. melanurus (E). A total of four combinations (ACD, ACE, BCD and BCE) were tested with 770 gene trees, each dataset consisting of three populations and an outgroup, Z. senegalensis. Each population was randomly subsampled to at most four individuals with ten replicates per locus: four out of five Z. auriventer individuals, four out of ten Z. simplex individuals, four out of seven Z. melanurus individuals, one out of two Z. atricapilla samples, and one Z. emiliae individual. We generated 48 possible demographic models under the following settings: an overall maximum of three free parameters (K = 3), a maximum of two coalescent events (K = 2), either complete isolation or migration event(s) with a single rate (K = 1), no variation in population size and growth (K = 1), only fully resolved trees were assumed, and only symmetrical migration between populations was assumed. Each dataset was simulated under these 48 different models with 10,000 trees using grid search. The probability of observing the set of gene trees under each model was assessed using Akaike Information Criterion (AIC) scores and the top two models with the lowest AIC values for each dataset were plotted in R 3.6.1 (R Development Core Team, 2019). We included the second-best model for each combination in our assessment because a previous study has shown that the second-best model may emerge as the correct demographic model with further testing (Morales et al., 2017). Additionally, we computed the genealogical divergence index (gdi) of each combination to assess the overall level of divergence between the sister taxa inferred from each model with the combined effects of genetic drift and gene flow (Jackson et al., 2017a). A gdi index close to 0 suggests panmixia, while an index close to one suggests strong divergence.

Ancestral range estimation

Request a detailed protocol

We ran BioGeoBEARS 1.1.2 (Matzke, 2014) in R 3.6.1 (R Development Core Team, 2019) to estimate the ancestral range of each internal node of the unrooted ML concatenated tree. We used ETE Toolkit 3.1.1 (Huerta-Cepas et al., 2016) python script to set the root node between the Indo-African clade and the other clades. We assigned each species to one or more out of ten geographical areas: the Afrotropical region, West Indian Ocean islands, South Asia, mainland Southeast Asia, East Asia, Philippines, the Sundaic region, Wallacea, Melanesia, and Australia. We applied all six models (DEC, DEC+j, DIVA, DIVA+j, BayArea and BayArea+j) and selected the best-fitting model (DEC+j) based on AIC values.

References

1. 1
   Comparative phylogeography of Australo-Papuan mangrove-restricted and mangrove-associated avifaunas

   1. N Á.S
   2. L Joseph

   (2013)

   Biological Journal of the Linnean Society 109:574–598.

   https://doi.org/10.1111/bij.12082

   • Google Scholar
2. 2
   Basic local alignment search tool

   1. SF Altschul
   2. W Gish
   3. W Miller
   4. EW Myers
   5. DJ Lipman

   (1990)

   Journal of Molecular Biology 215:403–410.

   https://doi.org/10.1016/S0022-2836(05)80360-2

   • PubMed
   • Google Scholar
3. 3
   New records and range extensions of birds from Timor, Alor and Rote

   1. H Ashari
   2. DM Prawiradilaga
   3. JA Eaton
   4. S Suparno
   5. FE Rheindt

   (2019)

   TREUBIA 45:47–64.

   https://doi.org/10.14203/treubia.v45i0.3628

   • Google Scholar
4. 4
   SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing

   1. A Bankevich
   2. S Nurk
   3. D Antipov
   4. AA Gurevich
   5. M Dvorkin
   6. AS Kulikov
   7. VM Lesin
   8. SI Nikolenko
   9. S Pham
   10. AD Prjibelski
   11. AV Pyshkin
   12. AV Sirotkin
   13. N Vyahhi
   14. G Tesler
   15. MA Alekseyev
   16. PA Pevzner

   (2012)

   Journal of Computational Biology 19:455–477.

   https://doi.org/10.1089/cmb.2012.0021

   • PubMed
   • Google Scholar
5. 5
   Using historical genome‐wide DNA to unravel the confused taxonomy in a songbird lineage that is extinct in the wild

   1. P Baveja
   2. KM Garg
   3. B Chattopadhyay
   4. KR Sadanandan
   5. DM Prawiradilaga
   6. P Yuda
   7. JGH Lee
   8. FE Rheindt

   (2020)

   Evolutionary Applications 4:1–12.

   https://doi.org/10.1111/eva.13149

   • Google Scholar
6. 6
   Transcriptome-based exon capture enables highly cost-effective comparative genomic data collection at moderate evolutionary scales

   1. K Bi
   2. D Vanderpool
   3. S Singhal
   4. T Linderoth
   5. C Moritz
   6. JM Good

   (2012)

   BMC Genomics 13:403.

   https://doi.org/10.1186/1471-2164-13-403

   • Google Scholar
7. 7
   Trimmomatic: a flexible trimmer for Illumina sequence data

   1. AM Bolger
   2. M Lohse
   3. B Usadel

   (2014)

   Bioinformatics 30:2114–2120.

   https://doi.org/10.1093/bioinformatics/btu170

   • Google Scholar
8. 8
   Target enrichment of thousands of ultraconserved elements sheds new light on early relationships within New World sparrows (Aves: Passerellidae)

   1. RW Bryson
   2. BC Faircloth
   3. WLE Tsai
   4. JE McCormack
   5. J Klicka

   (2016)

   The Auk 133:451–458.

   https://doi.org/10.1642/AUK-16-26.1

   • Google Scholar
9. 9
   Near-complete phylogeny and taxonomic revision of the world’s babblers (Aves: Passeriformes)

   1. T Cai
   2. A Cibois
   3. P Alström
   4. RG Moyle
   5. JD Kennedy
   6. S Shao
   7. R Zhang
   8. M Irestedt
   9. PGP Ericson
   10. M Gelang
   11. Y Qu
   12. F Lei
   13. J Fjeldså

   (2019)

   Molecular Phylogenetics and Evolution 130:346–356.

   https://doi.org/10.1016/j.ympev.2018.10.010

   • Google Scholar
10. 10
    The role of evolutionary time, diversification rates and dispersal in determining the global diversity of a large radiation of passerine birds

    1. T Cai
    2. S Shao
    3. JD Kennedy
    4. P Alström
    5. RG Moyle
    6. Y Qu
    7. F Lei
    8. J Fjeldså

    (2020)

    Journal of Biogeography 47:1612–1625.

    https://doi.org/10.1111/jbi.13823

    • Google Scholar
11. 11
    BLAST+: architecture and applications

    1. C Camacho
    2. G Coulouris
    3. V Avagyan
    4. N Ma
    5. J Papadopoulos
    6. K Bealer
    7. TL Madden

    (2009)

    BMC Bioinformatics 10:421.

    https://doi.org/10.1186/1471-2105-10-421

    • Google Scholar
12. 12
    Historic DNA reveals Anthropocene threat to a tropical urban fruit bat

    1. B Chattopadhyay
    2. KM Garg
    3. IH Mendenhall
    4. FE Rheindt

    (2019)

    Current Biology 29:R1299–R1300.

    https://doi.org/10.1016/j.cub.2019.11.013

    • Google Scholar
13. 13
    Microevolution in island forms: the roles of drift and directional selection in morphological divergence of a passerine bird

    1. SM Clegg
    2. SM Degnan
    3. C Moritz
    4. A Estoup
    5. J Kikkawa
    6. IPF Owens

    (2002)

    Evolution 56:2090–2099.

    https://doi.org/10.1111/j.0014-3820.2002.tb00134.x

    • Google Scholar
14. 14
    The genome of the "Great Speciator" Provides Insights into Bird Diversification

    1. L Cornetti
    2. LM Valente
    3. LT Dunning
    4. X Quan
    5. RA Black
    6. O Hébert
    7. V Savolainen

    (2015)

    Genome Biology and Evolution 7:2680–2691.

    https://doi.org/10.1093/gbe/evv168

    • PubMed
    • Google Scholar
15. 15
    Niche divergence promotes rapid diversification of East African sky island white-eyes (Aves: Zosteropidae)

    1. SC Cox
    2. RP Prys-Jones
    3. JC Habel
    4. BA Amakobe
    5. JJ Day

    (2014)

    Molecular Ecology 23:4103–4118.

    https://doi.org/10.1111/mec.12840

    • Google Scholar
16. 16
    Quaternary land bridges have not been universal conduits of gene flow

    1. E Cros
    2. B Chattopadhyay
    3. KM Garg
    4. NSR Ng
    5. S Tomassi
    6. S Benedick
    7. DP Edwards
    8. FE Rheindt

    (2020)

    Molecular Ecology 29:2692–2706.

    https://doi.org/10.1111/mec.15509

    • Google Scholar
17. 17
    Borneo and Indochina are major evolutionary hotspots for southeast asian biodiversity

    1. M de Bruyn
    2. B Stelbrink
    3. RJ Morley
    4. R Hall
    5. GR Carvalho
    6. CH Cannon
    7. G van den Bergh
    8. E Meijaard
    9. I Metcalfe
    10. L Boitani
    11. L Maiorano
    12. R Shoup
    13. T von Rintelen

    (2014)

    Systematic Biology 63:879–901.

    https://doi.org/10.1093/sysbio/syu047

    • Google Scholar
18. 18
    The perils of single gene trees — mitochondrial versus single-copy nuclear DNA variation in white-eyes (Aves: Zosteropidae)

    1. SM Degnan

    (1993)

    Molecular Ecology 2:219–225.

    https://doi.org/10.1111/j.1365-294X.1993.tb00011.x

    • Google Scholar
19. 19
    Phylogeography of mitochondrial DNA in two species of White-Eyes in Australia

    1. SM Degnan
    2. C Moritz

    (1992)

    The Auk 109:800–811.

    https://doi.org/10.2307/4088154

    • Google Scholar
20. 20
    Illustrated Checklist of the Birds of the World, Volume 2: Passerines

    1. J del Hoyo
    2. NJ Collar
    3. DA Christie
    4. A Elliott
    5. LDC Fishpool
    6. P Boesman
    7. GM Kirwan

    (2016)

    Cambridge: Lynx Edicions, Barcelona & BirdLife International.

    • Google Scholar
21. 21
    Species-distance relation for birds of the solomon archipelago, and the paradox of the great speciators

    1. JM Diamond
    2. ME Gilpin
    3. E Mayr

    (1976)

    PNAS 73:2160–2164.

    https://doi.org/10.1073/pnas.73.6.2160

    • PubMed
    • Google Scholar
22. 22
    A short survey of the meratus mountains, South Kalimantan Province, Indonesia: two undescribed avian species discovered

    1. JA Eaton
    2. SL Mitchell
    3. CNG Bocos
    4. FE Rheindt

    (2016)

    BirdingASIA 26:107–113.

    • Google Scholar
23. 23
    Speciation in birds: genes, geography, and sexual selection

    1. SV Edwards
    2. SB Kingan
    3. JD Calkins
    4. CN Balakrishnan
    5. WB Jennings
    6. WJ Swanson
    7. MD Sorenson

    (2005)

    PNAS 102:6550–6557.

    https://doi.org/10.1073/pnas.0501846102

    • PubMed
    • Google Scholar
24. 24
    High-resolution species trees without concatenation

    1. SV Edwards
    2. L Liu
    3. DK Pearl

    (2007)

    PNAS 104:5936–5941.

    https://doi.org/10.1073/pnas.0607004104

    • PubMed
    • Google Scholar
25. 25
    Reticulation, divergence, and the phylogeography-phylogenetics continuum

    1. SV Edwards
    2. S Potter
    3. CJ Schmitt
    4. JG Bragg
    5. C Moritz

    (2016)

    PNAS 113:8025–8032.

    https://doi.org/10.1073/pnas.1601066113

    • PubMed
    • Google Scholar
26. 26
    The genomic landscape of species divergence in Ficedula flycatchers

    1. H Ellegren
    2. L Smeds
    3. R Burri
    4. PI Olason
    5. N Backström
    6. T Kawakami
    7. A Künstner
    8. H Mäkinen
    9. K Nadachowska-Brzyska
    10. A Qvarnström
    11. S Uebbing
    12. JBW Wolf

    (2012)

    Nature 491:756–760.

    https://doi.org/10.1038/nature11584

    • Google Scholar
27. 27
    fastsimcoal: a continuous-time coalescent simulator of genomic diversity under arbitrarily complex evolutionary scenarios

    1. L Excoffier
    2. M Foll

    (2011)

    Bioinformatics 27:1332–1334.

    https://doi.org/10.1093/bioinformatics/btr124

    • Google Scholar
28. 28
    Analysis of whole mitogenomes from ancient samples

    1. GG Fortes
    2. JLA Paijmans

    (2015)

    In: T Kroneis, editors. Whole Genome Amplification. Springer. pp. 179–195.

    https://doi.org/10.1007/978-1-4939-2990-0_13

    • Google Scholar
29. 29
    Pleistocene land bridges act as semipermeable agents of avian gene flow in wallacea

    1. KM Garg
    2. B Chattopadhyay
    3. PR Wilton
    4. D Malia Prawiradilaga
    5. FE Rheindt

    (2018)

    Molecular Phylogenetics and Evolution 125:196–203.

    https://doi.org/10.1016/j.ympev.2018.03.032

    • PubMed
    • Google Scholar
30. 30
    A draft sequence of the neandertal genome

    1. RE Green
    2. J Krause
    3. AW Briggs
    4. T Maricic
    5. U Stenzel
    6. M Kircher
    7. N Patterson
    8. H Li
    9. W Zhai
    10. MH Fritz
    11. NF Hansen
    12. EY Durand
    13. AS Malaspinas
    14. JD Jensen
    15. T Marques-Bonet
    16. C Alkan
    17. K Prüfer
    18. M Meyer
    19. HA Burbano
    20. JM Good
    21. R Schultz
    22. A Aximu-Petri
    23. A Butthof
    24. B Höber
    25. B Höffner
    26. M Siegemund
    27. A Weihmann
    28. C Nusbaum
    29. ES Lander
    30. C Russ
    31. N Novod
    32. J Affourtit
    33. M Egholm
    34. C Verna
    35. P Rudan
    36. D Brajkovic
    37. Ž Kucan
    38. I Gušic
    39. VB Doronichev
    40. LV Golovanova
    41. C Lalueza-Fox
    42. M de la Rasilla
    43. J Fortea
    44. A Rosas
    45. RW Schmitz
    46. PLF Johnson
    47. EE Eichler
    48. D Falush
    49. E Birney
    50. JC Mullikin
    51. M Slatkin
    52. R Nielsen
    53. J Kelso
    54. M Lachmann
    55. D Reich
    56. S Pääbo

    (2010)

    Science 328:710–722.

    https://doi.org/10.1126/science.1188021

    • PubMed
    • Google Scholar
31. 31
    Diffusion approximations for demographic inference: dadi

    1. R Gutenkunst
    2. R Hernandez
    3. S Williamson
    4. C Bustamante

    (2010)

    Nature Precedings 1:1.

    https://doi.org/10.1038/npre.2010.4594.1

    • Google Scholar
32. 32
    Bioacoustic and multi-locus DNA data of Ninox owls support high incidence of extinction and recolonisation on small, low-lying islands across Wallacea

    1. CY Gwee
    2. L Christidis
    3. JA Eaton
    4. JA Norman
    5. CR Trainor
    6. P Verbelen
    7. FE Rheindt

    (2017)

    Molecular Phylogenetics and Evolution 109:246–258.

    https://doi.org/10.1016/j.ympev.2016.12.024

    • Google Scholar
33. 33
    Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, model and animations

    1. R Hall

    (2002)

    Journal of Asian Earth Sciences 20:353–431.

    https://doi.org/10.1016/S1367-9120(01)00069-4

    • Google Scholar
34. 34
    Late Jurassic–Cenozoic reconstructions of the Indonesian region and the Indian Ocean

    1. R Hall

    (2012)

    Tectonophysics 570-571:1–41.

    https://doi.org/10.1016/j.tecto.2012.04.021

    • Google Scholar
35. 35
    Birds of paradise, biogeography and ecology in New Guinea: a review

    1. M Heads

    (2001)

    Journal of Biogeography 28:893–925.

    https://doi.org/10.1046/j.1365-2699.2001.00600.x

    • Google Scholar
36. 36
    ETE 3: reconstruction, analysis, and visualization of phylogenomic data

    1. J Huerta-Cepas
    2. F Serra
    3. P Bork

    (2016)

    Molecular Biology and Evolution 33:1635–1638.

    https://doi.org/10.1093/molbev/msw046

    • PubMed
    • Google Scholar
37. 37
    Four new mitochondrial genomes of the genus zosterops (aves: passeriformes: zosteropidae) from East Africa with a phylogenetic evaluation of the group

    1. M Husemann
    2. S Sturm
    3. M Curto
    4. H Meimberg
    5. JC Habel

    (2016)

    Mitochondrial DNA Part B 1:544–548.

    https://doi.org/10.1080/23802359.2016.1198937

    • Google Scholar
38. 38
    The spatio-temporal colonization and diversification across the Indo-Pacific by a ‘great speciator’ (Aves, Erythropitta erythrogaster)

    1. M Irestedt
    2. P-H Fabre
    3. H Batalha-Filho
    4. KA Jønsson
    5. CS Roselaar
    6. G Sangster
    7. PGP Ericson

    (2013)

    Proceedings of the Royal Society B: Biological Sciences 280:20130309.

    https://doi.org/10.1098/rspb.2013.0309

    • Google Scholar
39. 39
    Species delimitation with gene flow

    1. ND Jackson
    2. BC Carstens
    3. AE Morales
    4. BC O'Meara

    (2017a)

    Systematic Biology 66:799–812.

    https://doi.org/10.1093/sysbio/syw117

    • PubMed
    • Google Scholar
40. 40
    PHRAPL: phylogeographic inference using approximate likelihoods

    1. ND Jackson
    2. AE Morales
    3. BC Carstens
    4. BC O'Meara

    (2017b)

    Systematic Biology 66:1045–1053.

    https://doi.org/10.1093/sysbio/syx001

    • PubMed
    • Google Scholar
41. 41
    The global diversity of birds in space and time

    1. W Jetz
    2. GH Thomas
    3. JB Joy
    4. K Hartmann
    5. AO Mooers

    (2012)

    Nature 491:444–448.

    https://doi.org/10.1038/nature11631

    • Google Scholar
42. 42
    HybPiper: extracting coding sequence and introns for phylogenetics from High-Throughput sequencing reads using target enrichment

    1. MG Johnson
    2. EM Gardner
    3. Y Liu
    4. R Medina
    5. B Goffinet
    6. AJ Shaw
    7. NJC Zerega
    8. NJ Wickett

    (2016)

    Applications in Plant Sciences 4:1600016.

    https://doi.org/10.3732/apps.1600016

    • Google Scholar
43. 43
    mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters

    1. H Jónsson
    2. A Ginolhac
    3. M Schubert
    4. PLF Johnson
    5. L Orlando

    (2013)

    Bioinformatics 29:1682–1684.

    https://doi.org/10.1093/bioinformatics/btt193

    • Google Scholar
44. 44
    MAFFT multiple sequence alignment software version 7: improvements in performance and usability

    1. K Katoh
    2. DM Standley

    (2013)

    Molecular Biology and Evolution 30:772–780.

    https://doi.org/10.1093/molbev/mst010

    • PubMed
    • Google Scholar
45. 45
    ANGSD: analysis of next generation sequencing data

    1. TS Korneliussen
    2. A Albrechtsen
    3. R Nielsen

    (2014)

    BMC Bioinformatics 15:356.

    https://doi.org/10.1186/s12859-014-0356-4

    • PubMed
    • Google Scholar
46. 46
    Evolution of Darwin’s finches and their beaks revealed by genome sequencing

    1. S Lamichhaney
    2. J Berglund
    3. MS Almén
    4. K Maqbool
    5. M Grabherr
    6. A Martinez-Barrio
    7. M Promerová
    8. C-J Rubin
    9. C Wang
    10. N Zamani
    11. BR Grant
    12. PR Grant
    13. MT Webster
    14. L Andersson

    (2015)

    Nature 518:371–375.

    https://doi.org/10.1038/nature14181

    • Google Scholar
47. 47
    Fast gapped-read alignment with bowtie 2

    1. B Langmead
    2. SL Salzberg

    (2012)

    Nature Methods 9:357–359.

    https://doi.org/10.1038/nmeth.1923

    • PubMed
    • Google Scholar
48. 48
    Multilocus resolution of phylogeny and timescale in the extant adaptive radiation of hawaiian honeycreepers

    1. HR Lerner
    2. M Meyer
    3. HF James
    4. M Hofreiter
    5. RC Fleischer

    (2011)

    Current Biology 21:1838–1844.

    https://doi.org/10.1016/j.cub.2011.09.039

    • PubMed
    • Google Scholar
49. 49
    A bird’s white-eye view on neosex chromosome evolution

    1. T Leroy
    2. Y Anselmetti
    3. M-K Tilak
    4. S Bérard
    5. L Csukonyi
    6. M Gabrielli
    7. C Scornavacca
    8. B Milá
    9. C Thébaud
    10. B Nabholz

    (2019)

    bioRxiv.

    https://doi.org/10.1101/505610

    • Google Scholar
50. 50
    The sequence alignment/Map format and SAMtools

    1. H Li
    2. B Handsaker
    3. A Wysoker
    4. T Fennell
    5. J Ruan
    6. N Homer
    7. G Marth
    8. G Abecasis
    9. R Durbin
    10. 1000 Genome Project Data Processing Subgroup

    (2009)

    Bioinformatics 25:2078–2079.

    https://doi.org/10.1093/bioinformatics/btp352

    • PubMed
    • Google Scholar
51. 51
    EvolMarkers: a database for mining exon and intron markers for evolution, ecology and conservation studies

    1. C Li
    2. JJ Riethoven
    3. GJ Naylor

    (2012)

    Molecular Ecology Resources 12:967–971.

    https://doi.org/10.1111/j.1755-0998.2012.03167.x

    • PubMed
    • Google Scholar
52. 52
    Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM

    1. H Li

    (2013)

    arXiv.

    • Google Scholar
53. 53
    Molecular evidence suggests radical revision of species limits in the great speciator white-eye genus Zosterops

    1. BTM Lim
    2. KR Sadanandan
    3. C Dingle
    4. YY Leung
    5. DM Prawiradilaga
    6. M Irham
    7. H Ashari
    8. JGH Lee
    9. FE Rheindt

    (2019)

    Journal of Ornithology 160:1–16.

    https://doi.org/10.1007/s10336-018-1583-7

    • Google Scholar
54. 54
    Genetic differentiation within a widespread "supertramp" taxon: Molecular phylogenetics of the Louisiade White-eye (Zosterops griseotinctus)

    1. E Linck
    2. S Schaack
    3. JP Dumbacher

    (2016)

    Molecular Phylogenetics and Evolution 94:113–121.

    https://doi.org/10.1016/j.ympev.2015.08.018

    • PubMed
    • Google Scholar
55. 55
    Estimating species phylogenies using coalescence times among sequences

    1. L Liu
    2. L Yu
    3. DK Pearl
    4. SV Edwards

    (2009)

    Systematic Biology 58:468–477.

    https://doi.org/10.1093/sysbio/syp031

    • PubMed
    • Google Scholar
56. 56
    A maximum pseudo-likelihood approach for estimating species trees under the coalescent model

    1. L Liu
    2. L Yu
    3. SV Edwards

    (2010)

    BMC Evolutionary Biology 10:302.

    https://doi.org/10.1186/1471-2148-10-302

    • PubMed
    • Google Scholar
57. 57
    Coalescent methods are robust to the simultaneous effects of long branches and incomplete lineage sorting

    1. L Liu
    2. Z Xi
    3. CC Davis

    (2015)

    Molecular Biology and Evolution 32:791–805.

    https://doi.org/10.1093/molbev/msu331

    • PubMed
    • Google Scholar
58. 58
    Modern Phylogenomics: Building Phylogenetic Trees Using the Multispecies Coalescent Model

    1. L Liu
    2. C Anderson
    3. D Pearl
    4. SV Edwards

    (2019)

    In: M Anisimova, editors. Evolutionary Genomics: Statistical and Computational Methods. New York: Springer. pp. 211–239.

    https://doi.org/10.1007/978-1-4939-9074-0_7

    • Google Scholar
59. 59
    Biogeography of the Indo-Australian archipelago

    1. DJ Lohman
    2. M de Bruyn
    3. T Page
    4. K von Rintelen
    5. R Hall
    6. PKL Ng
    7. H-T Shih
    8. GR Carvalho
    9. T von Rintelen

    (2011)

    Annual Review of Ecology, Evolution, and Systematics 42:205–226.

    https://doi.org/10.1146/annurev-ecolsys-102710-145001

    • Google Scholar
60. 60
    The Theory of Island Biogeography

    1. RH MacArthur
    2. EO Wilson

    (2001)

    New Jersey: Princeton University Press.

    • Google Scholar
61. 61
    FLASH: fast length adjustment of short reads to improve genome assemblies

    1. T Magoč
    2. SL Salzberg

    (2011)

    Bioinformatics 27:2957–2963.

    https://doi.org/10.1093/bioinformatics/btr507

    • PubMed
    • Google Scholar
62. 62
    Gene flow and rapid differentiation characterize a rapid insular radiation in the southwest Pacific (Aves: Zosterops)

    1. JD Manthey
    2. CH Oliveros
    3. MJ Andersen
    4. CE Filardi
    5. RG Moyle

    (2020)

    Evolution 74:1788–1803.

    https://doi.org/10.1111/evo.14043

    • Google Scholar
63. 63
    Cutadapt removes adapter sequences from high-throughput sequencing reads

    1. M Martin

    (2011)

    EMBnet.journal 17:10–12.

    https://doi.org/10.14806/ej.17.1.200

    • Google Scholar
64. 64
    A comprehensive molecular phylogeny of Afrotropical white-eyes (Aves: Zosteropidae) highlights prior underestimation of mainland diversity and complex colonisation history

    1. FC Martins
    2. SC Cox
    3. M Irestedt
    4. RP Prŷs-Jones
    5. JJ Day

    (2020)

    Molecular Phylogenetics and Evolution 149:106843.

    https://doi.org/10.1016/j.ympev.2020.106843

    • Google Scholar
65. 65
    Model selection in historical biogeography reveals that founder-event speciation is a crucial process in island clades

    1. NJ Matzke

    (2014)

    Systematic Biology 63:951–970.

    https://doi.org/10.1093/sysbio/syu056

    • PubMed
    • Google Scholar
66. 66
    Relationships among Indo-Australian Zosteropidae (Aves)

    1. E Mayr

    (1965)

    In: N Wright, editors. Breviora. Cambridge, Massachusetts, USA: Museum of Comparative Zoology, Harvard University. pp. 1–6.

    • Google Scholar
67. 67
    A systematic review of the Indo-Australian Zosteropidae (Part I)

    1. GF Mees

    (1957)

    Zoologische Verhandelingen 35:1–204.

    • Google Scholar
68. 68
    Speciation with gene flow in north american Myotis bats

    1. AE Morales
    2. ND Jackson
    3. TA Dewey
    4. BC O'Meara
    5. BC Carstens

    (2017)

    Systematic Biology 66:440–452.

    https://doi.org/10.1093/sysbio/syw100

    • PubMed
    • Google Scholar
69. 69
    Biogeographic implications of the Tertiary palaeogeographic evolution of Sulawesi and Borneo

    1. SJ Moss
    2. MEJ Wilson

    (1998)

    Biogeography and Geological Evolution of SE Asia 133:163.

    • Google Scholar
70. 70
    Explosive pleistocene diversification and hemispheric expansion of a "great speciator"

    1. RG Moyle
    2. CE Filardi
    3. CE Smith
    4. J Diamond

    (2009)

    PNAS 106:1863–1868.

    https://doi.org/10.1073/pnas.0809861105

    • PubMed
    • Google Scholar
71. 71
    The effects of pleistocene climate change on biotic differentiation in a montane songbird clade from wallacea

    1. NSR Ng
    2. PR Wilton
    3. DM Prawiradilaga
    4. YC Tay
    5. M Indrawan
    6. KM Garg
    7. FE Rheindt

    (2017)

    Molecular Phylogenetics and Evolution 114:353–366.

    https://doi.org/10.1016/j.ympev.2017.05.007

    • PubMed
    • Google Scholar
72. 72
    Late cenozoic palaeogeography of Sulawesi, Indonesia

    1. AMS Nugraha
    2. R Hall

    (2018)

    Palaeogeography, Palaeoclimatology, Palaeoecology 490:191–209.

    https://doi.org/10.1016/j.palaeo.2017.10.033

    • Google Scholar
73. 73
    A multi-locus phylogeny reveals a complex pattern of diversification related to climate and habitat heterogeneity in southern African white-eyes

    1. G Oatley
    2. G Voelker
    3. TM Crowe
    4. RCK Bowie

    (2012)

    Molecular Phylogenetics and Evolution 64:633–644.

    https://doi.org/10.1016/j.ympev.2012.05.022

    • Google Scholar
74. 74
    Oligocene divergence of frogmouth birds (Podargidae) across Wallace's Line

    1. PM Oliver
    2. H Heiniger
    3. AF Hugall
    4. L Joseph
    5. KJ Mitchell

    (2020)

    Biology Letters 16:20200040.

    https://doi.org/10.1098/rsbl.2020.0040

    • Google Scholar
75. 75
    Earth history and the passerine superradiation

    1. CH Oliveros
    2. DJ Field
    3. DT Ksepka
    4. FK Barker
    5. A Aleixo
    6. MJ Andersen
    7. P Alström
    8. BW Benz
    9. EL Braun
    10. MJ Braun
    11. GA Bravo
    12. RT Brumfield
    13. RT Chesser
    14. S Claramunt
    15. J Cracraft
    16. AM Cuervo
    17. EP Derryberry
    18. TC Glenn
    19. MG Harvey
    20. PA Hosner
    21. L Joseph
    22. RT Kimball
    23. AL Mack
    24. CM Miskelly
    25. AT Peterson
    26. MB Robbins
    27. FH Sheldon
    28. LF Silveira
    29. BT Smith
    30. ND White
    31. RG Moyle
    32. BC Faircloth

    (2019)

    PNAS 116:7916–7925.

    https://doi.org/10.1073/pnas.1813206116

    • PubMed
    • Google Scholar
76. 76
    A sympatric pair of undescribed white-eye species (Aves: zosteropidae: Zosterops) with different origins

    1. DP O’Connell
    2. DJ Kelly
    3. N Lawless
    4. K O’Brien
    5. FÓ Marcaigh
    6. A Karya
    7. K Analuddin
    8. NM Marples

    (2019)

    Zoological Journal of the Linnean Society 186:701–724.

    https://doi.org/10.1093/zoolinnean/zlz022

    • Google Scholar
77. 77
    Ancient admixture in human history

    1. N Patterson
    2. P Moorjani
    3. Y Luo
    4. S Mallick
    5. N Rohland
    6. Y Zhan
    7. T Genschoreck
    8. T Webster
    9. D Reich

    (2012)

    Genetics 192:1065–1093.

    https://doi.org/10.1534/genetics.112.145037

    • PubMed
    • Google Scholar
78. 78
    BEDTools: a flexible suite of utilities for comparing genomic features

    1. AR Quinlan
    2. IM Hall

    (2010)

    Bioinformatics 26:841–842.

    https://doi.org/10.1093/bioinformatics/btq033

    • Google Scholar
79. 79
    R: A Language and Environment for Statistical Computing Computer Program, version 2.6.2

    1. R Development Core Team

    (2019)

    R Foundation for Statistical Computing, Vienna, Austria.
80. 80
    Leap‐frog dispersal and mitochondrial introgression: Phylogenomics and biogeography of Limnonectes fanged frogs in the Lesser Sundas Archipelago of Wallacea

    1. SB Reilly
    2. AL Stubbs
    3. BR Karin
    4. K Bi
    5. E Arida
    6. DT Iskandar
    7. JA McGuire

    (2019)

    Journal of Biogeography 46:757–769.

    https://doi.org/10.1111/jbi.13526

    • Google Scholar
81. 81
    A lost world in Wallacea: description of a montane archipelagic avifauna

    1. FE Rheindt
    2. DM Prawiradilaga
    3. H Ashari
    4. Suparno
    5. CY Gwee
    6. GWX Lee
    7. MY Wu
    8. NSR Ng

    (2020)

    Science 367:167–170.

    https://doi.org/10.1126/science.aax2146

    • PubMed
    • Google Scholar
82. 82
    Genetic Introgression: An Integral but neglected component of speciation in birds

    1. FE Rheindt
    2. SV Edwards

    (2011)

    The Auk 128:620–632.

    https://doi.org/10.1525/auk.2011.128.4.620

    • Google Scholar
83. 83
    Disentangling avian diversity: South-East Asian mainland Oriental white-eye Zosterops palpebrosus constitutes two distinct lineages

    1. PD Round
    2. S Manawattana
    3. J Khudamrongsawat
    4. S Thunhikorn
    5. M Safoowong
    6. T Bhummakasikara

    (2017)

    Forktail 33:103–115.

    • Google Scholar
84. 84
    Subsiding sundaland

    1. A-C Sarr
    2. L Husson
    3. P Sepulchre
    4. A-M Pastier
    5. K Pedoja
    6. M Elliot
    7. C Arias-Ruiz
    8. T Solihuddin
    9. S Aribowo
    10. Susilohadi

    (2019)

    Geology 47:119–122.

    https://doi.org/10.1130/G45629.1

    • Google Scholar
85. 85
    Preliminary assessment of community composition and phylogeographic relationships of the birds of the Meratus Mountains, south-east borneo, Indonesia

    1. SB Shakya
    2. T Haryoko
    3. RC Burner
    4. DM Prawiradilaga
    5. FH Sheldon

    (2018)

    Bulletin of the British Ornithologists' Club 138:45–66.

    https://doi.org/10.25226/bboc.v138i1.2018.a6

    • Google Scholar
86. 86
    De novo transcriptomic analyses for non-model organisms: an evaluation of methods across a multi-species data set

    1. S Singhal

    (2013)

    Molecular Ecology Resources 13:403–416.

    https://doi.org/10.1111/1755-0998.12077

    • Google Scholar
87. 87
    Automated generation of heuristics for biological sequence comparison

    1. G Slater
    2. E Birney

    (2005)

    BMC Bioinformatics 6:31.

    https://doi.org/10.1186/1471-2105-6-31

    • Google Scholar
88. 88
    Phylogenetic relationships of micronesian White-Eyes based on mitochondrial sequence data

    1. B Slikas
    2. IB Jones
    3. SR Derrickson
    4. RC Fleischer

    (2000)

    The Auk 117:355–365.

    https://doi.org/10.1093/auk/117.2.355

    • Google Scholar
89. 89
    RepeatMasker 4.0

    1. A Smit
    2. R Hubley
    3. P Green

    (2015)

    Seattle, WA: Institute for Systems Biology.

    • Google Scholar
90. 90
    Resolving conflict in eutherian mammal phylogeny using phylogenomics and the multispecies coalescent model

    1. S Song
    2. L Liu
    3. SV Edwards
    4. S Wu

    (2012)

    PNAS 109:14942–14947.

    https://doi.org/10.1073/pnas.1211733109

    • PubMed
    • Google Scholar
91. 91
    RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies

    1. A Stamatakis

    (2014)

    Bioinformatics 30:1312–1313.

    https://doi.org/10.1093/bioinformatics/btu033

    • Google Scholar
92. 92
    Novel genome and genome-wide SNPs reveal early fragmentation effects in an edge-tolerant songbird population across an urbanized tropical metropolis

    1. DJX Tan
    2. B Chattopadhyay
    3. KM Garg
    4. E Cros
    5. PGP Ericson
    6. M Irestedt
    7. FE Rheindt

    (2018)

    Scientific Reports 8:1–12.

    https://doi.org/10.1038/s41598-018-31074-5

    • Google Scholar
93. 93
    DNA capture and next-generation sequencing can recover whole mitochondrial genomes from highly degraded samples for human identification

    1. JEL Templeton
    2. PM Brotherton
    3. B Llamas
    4. J Soubrier
    5. W Haak
    6. A Cooper
    7. JJ Austin

    (2013)

    Investigative Genetics 4:26.

    https://doi.org/10.1186/2041-2223-4-26

    • Google Scholar
94. 94
    ASTRID: accurate species TRees from internode distances

    1. P Vachaspati
    2. T Warnow

    (2015)

    BMC Genomics 16:S3.

    https://doi.org/10.1186/1471-2164-16-S10-S3

    • PubMed
    • Google Scholar
95. 95
    Whole mitochondrial genome capture from faecal samples and museum-preserved specimens

    1. T van der Valk
    2. F Lona Durazo
    3. L Dalén
    4. K Guschanski

    (2017)

    Molecular Ecology Resources 17:e111–e121.

    https://doi.org/10.1111/1755-0998.12699

    • Google Scholar
96. 96
    Maps of Pleistocene sea levels in Southeast Asia: shorelines, river systems and time durations

    1. HK Voris

    (2000)

    Journal of Biogeography 27:1153–1167.

    https://doi.org/10.1046/j.1365-2699.2000.00489.x

    • Google Scholar
97. 97
    The Malay Archipelago: The Land of the Orang-Utan and the Bird of Paradise; a Narrative of Travel, with Studies of Man and Nature

    1. AR Wallace

    (1962)

    New York: Courier Corporation.

    • Google Scholar
98. 98
    Immigration, species radiation and extinction in a highly diverse songbird lineage: white-eyes on indian ocean islands

    1. BH Warren
    2. E Bermingham
    3. RP Prys-Jones
    4. C Thébaud

    (2006)

    Molecular Ecology 15:3769–3786.

    https://doi.org/10.1111/j.1365-294X.2006.03058.x

    • PubMed
    • Google Scholar
99. 99
    Island Biogeography: Ecology, Evolution, and Conservation

    1. RJ Whittaker
    2. JM Fernández-Palacios

    (2007)

    Oxford: Oxford University Press.

    https://doi.org/10.1038/sj.hdy.6885841

    • Google Scholar
100. 100
     Non-sister sri lankan white-eyes (genus Zosterops) are a result of independent colonizations

     1. N Wickramasinghe
     2. VV Robin
     3. U Ramakrishnan
     4. S Reddy
     5. SS Seneviratne

     (2017)

     PLOS ONE 12:e0181441.

     https://doi.org/10.1371/journal.pone.0181441

     • PubMed
     • Google Scholar
101. 101
     FastUniq: a fast de novo duplicates removal tool for paired short reads

     1. H Xu
     2. X Luo
     3. J Qian
     4. X Pang
     5. J Song
     6. G Qian
     7. J Chen
     8. S Chen

     (2012)

     PLOS ONE 7:e52249.

     https://doi.org/10.1371/journal.pone.0052249

     • PubMed
     • Google Scholar
102. 102
     Gene flow analysis method, the D-statistic, is robust in a wide parameter space

     1. Y Zheng
     2. A Janke

     (2018)

     BMC Bioinformatics 19:10.

     https://doi.org/10.1186/s12859-017-2002-4

     • Google Scholar

Decision letter

1. Rosalyn Gloag

   Reviewing Editor; University of Sidney, Australia
2. George H Perry

   Senior Editor; Pennsylvania State University, United States
3. Leo Joseph

   Reviewer; CSIRO, Australia

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

Acceptance summary:

This study makes a valuable contribution to our understanding of the evolution of white-eyes, and more broadly the processes by which animals diverge and speciate. In particular, it highlights ways in which the biogeography of the Indonesian archipelago has generated the remarkable biodiversity of that region.

Decision letter after peer review:

Thank you for submitting your article "Phylogenomics of a 'Great Speciator' Reveals Indonesian Archipelago as the Center of Lineage Diversity" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and George Perry as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Leo Joseph (Reviewer #2)

We would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). Specifically, we are asking editors to accept without delay manuscripts, like yours, that they judge can stand as eLife papers without additional data, even if they feel that they would make the manuscript stronger. Thus the revisions requested below largely address clarity and presentation. In a few places, we suggest better justifying the current analyses and conclusions (e.g. explaining why an alternative analyses was not used).

All reviewers agreed that this study provides valuable insight into the phylogeography of Zosterops. It address the phylogenetic conflicts within the genus to resolve three main biogeographic clades. It further uses genomic data to provide evidence for introgression between species, highlighting that introgression is pervasive in rapidly diverging radiations. An impressive sample collection is used, spanning field and museum sampling.

Reviewers identified one main concern that should be addressed in a revision:

1) Borneo is posited as the origin of the radiation, but this notion could be better supported. Two suggestions were made:

i) Consider the region's tectonic history and include could some kind of modern biogeographic analysis, or a defense of why it is not included. For example, what about a BioGeoBears analysis?

ii) Consider alternative scenarios that might explain Borneo's species diversity; the apparent diversity in habitats of Borneo could have facilitated continental invaders to settle in the island. Perhaps, more clarification of the distribution of each species could provide better support to the authors' claims.

https://doi.org/10.7554/eLife.62765.sa1

Author response

Reviewers identified one main concern that should be addressed in a revision:

1) Borneo is posited as the origin of the radiation, but this notion could be better supported. Two suggestions were made:

i) Consider the region's tectonic history and include could some kind of modern biogeographic analysis, or a defense of why it is not included. For example, what about a BioGeoBears analysis?

We have accommodated this concern in multiple ways: (1) We have added a BioGeoBears analysis that demonstrates an Asian origin for white-eyes and an Indonesian origin for specific white-eye clades. (2) We have added more discussion on Earth history, specifically the Quaternary history of Borneo with its changing land connections during the sea level fluctuations of the Pleistocene. (Please note that ‘tectonic’ history per se is less relevant to Zosterops as significant tectonic changes have occurred at a pace that’s much slower than the relatively recent diversification of this genus). (3) We have strengthened explanations throughout the manuscript to make clear that we are not focusing on Borneo as the origin of the radiation, but as the center of modern-day diversity of this genus.

ii) Consider alternative scenarios that might explain Borneo's species diversity; the apparent diversity in habitats of Borneo could have facilitated continental invaders to settle in the island. Perhaps, more clarification of the distribution of each species could provide better support to the authors' claims.

We have accommodated this concern in multiple ways: 1) We have added a BioGeo Bears analysis that demonstrates an Asian origin for white-eyes and an Indonesian origin for specific white-eye clades. 2) We have added more discussion on Earth history, specifically the Quaternary history of Borneo with its changing land connections during the sea level fluctuations of the Pleistocene. (Please note that ‘tectonic’ history per se is less relevant to Zosterops as significant tectonic changes have occurred at a pace that’s much slower than the relatively recent diversification of this genus). 3) We have strengthened explanations throughout the manuscript to make clear that we are not focusing on Borneo as the origin of the radiation, but as the center of modern-day diversity of this genus.

https://doi.org/10.7554/eLife.62765.sa2

Article and author information

Author details

1. Chyi Yin Gwee

   National University of Singapore, Department of Biological Sciences, Singapore, Singapore

   Present address

   Division of Evolutionary Biology, Faculty of Biology, LMU Munich, Munich, Germany

   Contribution

   Resources, Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing

   Competing interests

   No competing interests declared

   

19. 0000-0001-8946-7085

Funding

Ministry of Education - Singapore (R-154-000-A59-112)

• Frank E Rheindt

Wildlife Reserves Singapore Conservation Fund (R-154-000-A99-592)

• Frank E Rheindt

Croeni Foundation (R-154-000-A05-592)

• Frank E Rheindt

SEABIG (R-154-000-648-646)

• Balaji Chattopadhyay

SEABIG (R-154-000-648-733)

• Balaji Chattopadhyay

University of Southampton (511206105)

• Kelvin S-H Peh

National Science Foundation (DEB-1441652)

• Rauri CK Bowie

National Science Foundation (DEB-1457845)

• Rauri CK Bowie

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We are indebted to the following staff and museums for loaning toepads and tissue material to us: P Sweet (American Museum of Natural History, New York, USA); A Drew and L Joseph (Australian National Wildlife Collection – Commonwealth Scientific and Industrial Research Organisation, Canberra, Australia); SM Birks (Burke Museum of Natural History and Culture, Seattle, USA); KP Lim (Lee Kong Chian Natural History Museum, Singapore); Museum of Vertebrate Zoology (Berkeley, USA); Museums Victoria (Melbourne, Australia); Museum Zoologicum Bogoriense (Bogor, West Java, Indonesia); P Kamminga (Naturalis Biodiversity Center, Leiden, Netherlands); M Penck (South Australian Museum, Adelaide, Australia); J Stigenberg and U Johansson (Swedish Museum of Natural History, Stockholm, Sweden); Western Australian Museum (Perth, Australia); and K Zyskowski (Yale Peabody Museum of Natural History, New Haven, USA). The Ministry of Research and Technology of Indonesia issued the following research permits: 10/TKPIPA/FRP/SM/X/2013 to FER; 314/SIP/FRP/E5/Dit.KI/X/2018 to FER and CYG; and 233/SIP/FRP/SM/VI/2013 and 147/SIP/FRP/SM/V/2014 to JM McGuire and RCKB. The Economic Planning Unit of Malaysia issued the research permit UPE 40/200/19/3295 to MCKS and KS-HP. We thank BTM Lim, EYX Ng, P Baveja, YF Chung, RYC Teo, E Arida, KC Rowe and JM McGuire for field, lab and/or analytical support, and the reviewers for their valuable comments. KMG acknowledges support from the Ramanujan Fellowship of SERB.

Senior Editor

1. George H Perry, Pennsylvania State University, United States

Reviewing Editor

1. Rosalyn Gloag, University of Sidney, Australia

Reviewer

1. Leo Joseph, CSIRO, Australia

Publication history

1. Received: September 3, 2020
2. Accepted: December 21, 2020
3. Accepted Manuscript published: December 22, 2020 (version 1)
4. Version of Record published: December 31, 2020 (version 2)

Copyright

© 2020, Gwee et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

National University of Singapore, Department of Biological Sciences, Singapore;

Max Planck Institute for Ornithology, Germany;

Division of Zoology, Research Center for Biology, Indonesian Institute of Sciences (LIPI), Cibinong Science Center, Indonesia;

Department of Bioinformatics and Genetics, Swedish Museum of Natural History, Sweden;

Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, China;

Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, China;

Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, Berkeley, United States;

Wildlife Reserves Singapore, Singapore;

University of Western Australia, School of Biological Sciences, Australia

Species density map of the genus Zosterops across the Indonesian archipelago (main map) and across the entire distribution range (bottom left inset).

Islands referred to in the text are specifically labeled on the main map. The total number of Zosterops species is shown in brackets beside each labeled region on the inset. We adopted del Hoyo et …

Phylogeny of Zosterops. Trees generated by (a) MP-EST, (b) ASTRID, (c) STAR and (d) the concatenation method (refer to Figure 2—figure supplement 1 for a full concatenated tree and Figure 2—figure supplement 2 for ancestral range estimation).

All nodes are supported by a bootstrap value of 100 unless otherwise stated. Nodes with less than 68% bootstrap support were collapsed. The three main clades are color-coded blue (Australasian), red …

Detection of secondary gene flow in a Sundaic species pair of unresolved placement (gray) with members of either the Asiatic (red) or Australasian (blue) clade.

(a) Placement of the Sundaic pair (Z. atricapilla and Z. auriventer) conflicts between MP-EST species tree and concatenated tree. (b) The top two demographic models in PHRAPL simulations show that …

0000-0003-1706-0520

Kritika M Garg

National University of Singapore, Department of Biological Sciences, Singapore, Singapore

Present address

Institute of Bioinformatics and Applied Biotechnology, Electronics City, India

Contribution

Resources, Formal analysis, Investigation, Methodology, Project administration, Writing - review and editing

Contributed equally with

Balaji Chattopadhyay

Competing interests

No competing interests declared

Balaji Chattopadhyay

National University of Singapore, Department of Biological Sciences, Singapore, Singapore

Present address

Trivedi School of Biosciences, Ashoka University, Sonipat, India

Contribution

Resources, Funding acquisition, Methodology, Project administration, Writing - review and editing

Contributed equally with

Kritika M Garg

Competing interests

No competing interests declared

Keren R Sadanandan

1. National University of Singapore, Department of Biological Sciences, Singapore, Singapore
2. Max Planck Institute for Ornithology, Seewiesen, Germany

Contribution

Resources, Investigation, Project administration, Writing - review and editing

Competing interests

No competing interests declared

Dewi M Prawiradilaga

Division of Zoology, Research Center for Biology, Indonesian Institute of Sciences (LIPI), Cibinong Science Center, Cibinong, Indonesia

Contribution

Resources, Project administration, Writing - review and editing

Competing interests

No competing interests declared

Martin Irestedt

Department of Bioinformatics and Genetics, Swedish Museum of Natural History, Stockholm, Sweden

Contribution

Resources, Writing - review and editing

Competing interests

No competing interests declared

Fumin Lei

1. Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
2. Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, China

Contribution

Resources, Writing - review and editing

Competing interests

No competing interests declared

Luke M Bloch

Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, Berkeley, Berkeley, United States

Contribution

Resources, Writing - review and editing

Competing interests

No competing interests declared

Jessica GH Lee

Wildlife Reserves Singapore, Singapore, Singapore

Contribution

Resources, Funding acquisition, Writing - review and editing

Competing interests

No competing interests declared

Mohammad Irham

Division of Zoology, Research Center for Biology, Indonesian Institute of Sciences (LIPI), Cibinong Science Center, Cibinong, Indonesia

Contribution

Resources, Writing - review and editing

Competing interests

No competing interests declared

Tri Haryoko

Division of Zoology, Research Center for Biology, Indonesian Institute of Sciences (LIPI), Cibinong Science Center, Cibinong, Indonesia

Contribution

Resources, Writing - review and editing

Competing interests

No competing interests declared

Malcolm CK Soh

University of Western Australia, School of Biological Sciences, Perth, Australia

Contribution

Resources, Writing - review and editing

Competing interests

No competing interests declared

Kelvin S-H Peh

University of Southampton, School of Biological Sciences, University, Southampton, United Kingdom

Contribution

Resources, Funding acquisition, Writing - review and editing

Competing interests

No competing interests declared

0000-0002-2921-1341

Karen MC Rowe

Sciences Department, Museums Victoria, Melbourne, Australia

Contribution

Resources, Writing - review and editing

Competing interests

No competing interests declared

0000-0002-6131-6418

Teuku Reza Ferasyi

1. Faculty of Veterinary Medicine, Universitas Syiah Kuala, Darussalam, Indonesia
2. Jiangsu Key Laboratory of Phylogenomics and Comparative Genomics, School of Life Sciences, Jiangsu Normal University, Xuzhou, China

Contribution

Resources, Writing - review and editing

Competing interests

No competing interests declared

Shaoyuan Wu

1. Department of Biochemistry and Molecular Biology, 2011 Collaborative Innovation Center of Tianjin for Medical Epigenetics, Tianjin Key Laboratory of Medical Epigenetics, School of Basic Medical Sciences, Tianjin Medical University, Tianjin, China
2. Center for Tropical Veterinary Studies – One Health Collaboration Center, Universitas Syiah Kuala, Darussalam, Indonesia

Contribution

Resources, Writing - review and editing

Competing interests

No competing interests declared

Guinevere OU Wogan

Museum of Vertebrate Zoology and Department of Environmental Science, Policy, and Management, University of California, Berkeley, Berkeley, United States

Contribution

Resources, Formal analysis, Methodology, Writing - review and editing

Competing interests

No competing interests declared

Rauri CK Bowie

Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, Berkeley, Berkeley, United States

Contribution

Resources, Funding acquisition, Investigation, Writing - review and editing

Competing interests

No competing interests declared

0000-0001-8328-6021

Frank E Rheindt

National University of Singapore, Department of Biological Sciences, Singapore, Singapore

Contribution

Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing

For correspondence

dbsrfe@nus.edu.sg

Competing interests

No competing interests declared