Waves of Colonization and Gene Flow in a Great Speciator.

Secondary contact between previously allopatric lineages offers a test of reproductive isolating mechanisms that may have accrued in isolation. Such instances of contact can produce stable hybrid zones-where reproductive isolation can further develop via reinforcement or phenotypic displacement-or result in the lineages merging. Ongoing secondary contact is most visible in continental systems, where steady input from parental taxa can occur readily. In oceanic island systems, however, secondary contact between closely related species of birds is relatively rare. When observed on sufficiently small islands, relative to population size, secondary contact likely represents a recent phenomenon. Here, we examine the dynamics of a group of birds whose apparent widespread hybridization influenced Ernst Mayr's foundational work on allopatric speciation: the whistlers of Fiji (Aves: Pachycephala). We demonstrate two clear instances of secondary contact within the Fijian archipelago, one resulting in a hybrid zone on a larger island, and the other resulting in a wholly admixed population on a smaller, adjacent island. We leveraged low genome-wide divergence in the hybrid zone to pinpoint a single genomic region associated with observed phenotypic differences. We use genomic data to present a new hypothesis that emphasizes rapid plumage evolution and post-divergence gene flow.

Systematics and biogeography of the whistlers (Aves: Pachycephalidae) inferred from ultraconserved elements and ancestral area reconstruction.

The utility of islands as natural laboratories of evolution is exemplified in the patterns of differentiation in widespread, phenotypically variable lineages. The whistlers (Aves: Pachycephalidae) are one of the most complex avian radiations, with a combination of widespread and locally endemic taxa spanning the vast archipelagos of the Indo-Pacific, making them an ideal group to study patterns and processes of diversification on islands. Here, we present a robust, species-level phylogeny of all five genera and 85% of species within Pachycephalidae, based on thousands of ultraconserved elements (UCEs) generated with a target-capture approach and high-throughput sequencing. We clarify phylogenetic relationships within Pachycephala and report on divergence timing and ancestral range estimation. We explored multiple biogeographic coding schemes that incorporated geological uncertainty in this complex region. The biogeographic origin of this group was difficult to discern, likely owing to aspects of dynamic Earth history in the Indo-Pacific. The Australo-Papuan region was the likely origin of crown-group whistlers, but the specific ancestral area could not be identified more precisely than Australia or New Guinea, and Wallacea may have played a larger role than previously realized in the evolutionary history of whistlers. Multiple independent colonizations of island archipelagos across Melanesia, Wallacea, and the Philippines contributed to the relatively high species richness of extant whistlers. This work refines our understanding of one of the regions' most celebrated bird lineages and adds to our growing knowledge about the patterns and processes of diversification in the Indo-Pacific.

Detecting turnover among complex communities using null models: a case study with sky-island haemosporidian parasites.

Turnover in species composition between sites, or beta diversity, is a critical component of species diversity that is typically influenced by geography, environment, and biotic interactions. Quantifying turnover is particularly challenging, however, in multi-host, multi-parasite assemblages where undersampling is unavoidable, resulting in inflated estimates of turnover and uncertainty about its spatial scale. We developed and implemented a framework using null models to test for community turnover in avian haemosporidian communities of three sky islands in the southwestern United States. We screened 776 birds for haemosporidian parasites from three genera (Parahaemoproteus, Plasmodium, and Leucocytozoon) by amplifying and sequencing a mitochondrial DNA barcode. We detected infections in 280 birds (36.1%), sequenced 357 infections, and found a total of 99 parasite haplotypes. When compared to communities simulated from a regional pool, we observed more unique, single-mountain haplotypes and fewer haplotypes shared among three mountain ranges than expected, indicating that haemosporidian communities differ to some degree among adjacent mountain ranges. These results were robust even after pruning datasets to include only identical sets of host species, and they were consistent for two of the three haemosporidian genera. The two more distant mountain ranges were more similar to each other than the one located centrally, suggesting that the differences we detected were due to stochastic colonization-extirpation dynamics. These results demonstrate that avian haemosporidian communities of temperate-zone forests differ on relatively fine spatial scales between adjacent sky islands. Null models are essential tools for testing the spatial scale of turnover in complex, undersampled, and poorly known systems.