Evolution of Chromosomal Inversions across an Avian Radiation.

Chromosomal inversions are structural mutations that can play a prominent role in adaptation and speciation. Inversions segregating across species boundaries (trans-species inversions) are often taken as evidence for ancient balancing selection or adaptive introgression, but can also be due to incomplete lineage sorting. Using whole-genome resequencing data from 18 populations of 11 recognized munia species in the genus Lonchura (N = 176 individuals), we identify four large para- and pericentric inversions ranging in size from 4 to 20 Mb. All four inversions cosegregate across multiple species and predate the numerous speciation events associated with the rapid radiation of this clade across the prehistoric Sahul (Australia, New Guinea) and Bismarck Archipelago. Using coalescent theory, we infer that trans-specificity is improbable for neutrally segregating variation despite substantial incomplete lineage sorting characterizing this young radiation. Instead, the maintenance of all three autosomal inversions (chr1, chr5, and chr6) is best explained by selection acting along ecogeographic clines not observed for the collinear parts of the genome. In addition, the sex chromosome inversion largely aligns with species boundaries and shows signatures of repeated positive selection for both alleles. This study provides evidence for trans-species inversion polymorphisms involved in both adaptation and speciation. It further highlights the importance of informing selection inference using a null model of neutral evolution derived from the collinear part of the genome.

Genetic architecture facilitates then constrains adaptation in a host-parasite coevolutionary arms race.

In coevolutionary arms races, interacting species impose selection on each other, generating reciprocal adaptations and counter adaptations. This process is typically enhanced by genetic recombination and heterozygosity, but these sources of evolutionary novelty may be secondarily lost when uniparental inheritance evolves to ensure the integrity of sex-linked adaptations. We demonstrate that host-specific egg mimicry in the African cuckoo finch Anomalospiza imberbis is maternally inherited, confirming the validity of an almost century-old hypothesis. We further show that maternal inheritance not only underpins the mimicry of different host species but also additional mimetic diversification that approximates the range of polymorphic egg "signatures" that have evolved within host species as an escalated defense against parasitism. Thus, maternal inheritance has enabled the evolution and maintenance of nested levels of mimetic specialization in a single parasitic species. However, maternal inheritance and the lack of sexual recombination likely disadvantage cuckoo finches by stifling further adaptation in the ongoing arms races with their individual hosts, which we show have retained biparental inheritance of egg phenotypes. The inability to generate novel genetic combinations likely prevents cuckoo finches from mimicking certain host phenotypes that are currently favored by selection (e.g., the olive-green colored eggs laid by some tawny-flanked prinia, Prinia subflava, females). This illustrates an important cost of coding coevolved adaptations on the nonrecombining sex chromosome, which may impede further coevolutionary change by effectively reversing the advantages of sexual reproduction in antagonistic coevolution proposed by the Red Queen hypothesis.

Chronically implantable LED arrays for behavioral optogenetics in primates.

Optogenetic methods have been widely used in rodent brains, but remain relatively under-developed for nonhuman primates such as rhesus macaques, an animal model with a large brain expressing sophisticated sensory, motor and cognitive behaviors. To address challenges in behavioral optogenetics in large brains, we developed Opto-Array, a chronically implantable array of light-emitting diodes for high-throughput optogenetic perturbation. We demonstrated that optogenetic silencing in the macaque primary visual cortex with the help of the Opto-Array results in reliable retinotopic visual deficits in a luminance discrimination task. We separately confirmed that Opto-Array illumination results in local neural silencing, and that behavioral effects are not due to tissue heating. These results demonstrate the effectiveness of the Opto-Array for behavioral optogenetic applications in large brains.

Inter- and intra-archipelago dynamics of population structure and gene flow in a Polynesian bird.

Islands are separated by natural barriers that prevent gene flow between terrestrial populations and promote allopatric diversification. Birds in the South Pacific are an excellent model to explore the interplay between isolation and gene flow due to the region's numerous archipelagos and well-characterized avian communities. The wattled honeyeater complex (Foulehaio spp.) comprises three allopatric species that are widespread and common across Fiji, Tonga, Samoa, and Wallis and Futuna. Here, we explored patterns of diversification within and among these lineages using genomic and morphometric data. We found support for three clades of Foulehaio corresponding to three recognized species. Within F. carunculatus, population genetic analyses identified nine major lineages, most of which were composed of sub-lineages that aligned nearly perfectly to individual island populations. Despite genetic structure and great geographic distance between populations, we found low levels of gene flow between populations in adjacent archipelagos. Additionally, body size of F. carunculatus varied randomly with respect to evolutionary history (as Ernst Mayr predicted), but correlated negatively with island size, consistent with the island rule. Our findings support a hypothesis that widespread taxa can show population structure between immediately adjacent islands, and likely represent many independent lineages loosely connected by gene flow.

ddRAD-seq data reveal significant genome-wide population structure and divergent genomic regions that distinguish the mallard and close relatives in North America.

Recently evolved species typically share genetic variation across their genomes due to incomplete lineage sorting and/or ongoing gene flow. Given only subtle allele frequency differences at most loci and the expectation that divergent selection may affect only a tiny fraction of the genome, distinguishing closely related species based on multi-locus data requires substantial genomic coverage. In this study, we used ddRAD-seq to sample the genomes of five recently diverged, New World "mallards" (Anas spp.), a group of dabbling duck species characterized by diagnosable phenotypic differences but minimal genetic differentiation. With increased genomic sampling, we aimed to characterize population structure within this group and identify genomic regions that may have experienced divergent selection during speciation. We analyzed 3,017 autosomal ddRAD-seq loci and 177 loci from the Z-chromosome. In contrast to previous studies, the ddRAD-seq data were sufficient to assign individuals to their respective species or subspecies and to generate estimates of gene flow in a phylogenetic framework. We find limited evidence of contemporary gene flow between the dichromatic mallard and several monochromatic taxa, but find evidence for historical gene flow between some monochromatic species pairs. We conclude that the overall genetic similarity of these taxa likely reflects retained ancestral polymorphism rather than recent and extensive gene flow. Thus, despite recurring cases of hybridization in this group, our results challenge the current dogma predicting the genetic extinction of the New World monochromatic dabbling ducks via introgressive hybridization with mallards. Moreover, ddRAD-seq data were sufficient to identify previously unknown outlier regions across the Z-chromosome and several autosomal chromosomes that may have been involved in the diversification of species in this recent radiation.

Persistence of an endangered native duck, feral mallards, and multiple hybrid swarms across the main Hawaiian Islands.

Interspecific hybridization is recognized as an important process in the evolutionary dynamics of both speciation and the reversal of speciation. However, our understanding of the spatial and temporal patterns of hybridization that erode versus promote species boundaries is incomplete. The endangered, endemic koloa maoli (or Hawaiian duck, Anas wyvilliana) is thought to be threatened with genetic extinction through ongoing hybridization with an introduced congener, the feral mallard (A. platyrhynchos). We investigated spatial and temporal variation in hybrid prevalence in populations throughout the main Hawaiian Islands, using genomic data to characterize population structure of koloa, quantify the extent of hybridization, and compare hybrid proportions over time. To accomplish this, we genotyped 3,308 double-digest restriction-site-associated DNA (ddRAD) loci in 425 putative koloa, mallards, and hybrids from populations across the main Hawaiian Islands. We found that despite a population decline in the last century, koloa genetic diversity is high. There were few hybrids on the island of Kaua'i, home to the largest population of koloa. By contrast, we report that sampled populations outside of Kaua'i can now be characterized as hybrid swarms, in that all individuals sampled were of mixed koloa × mallard ancestry. Further, there is some evidence that these swarms are stable over time. These findings demonstrate spatial variation in the extent and consequences of interspecific hybridization, and highlight how islands or island-like systems with small population sizes may be especially prone to genetic extinction when met with a congener that is not reproductively isolated.

Phylogenomics clarifies biogeographic and evolutionary history, and conservation status of West Indian tremblers and thrashers (Aves: Mimidae).

The West Indian avifauna has provided fundamental insights into island biogeography, taxon cycles, and the evolution of avian behavior. Our interpretations, however, should rely on robust hypotheses of evolutionary relationships and consistent conclusions about taxonomic status in groups with many endemic island populations. Here we present a phylogenetic study of the West Indian thrashers, tremblers, and allies, an assemblage of at least 5 species found on 29 islands, including what is considered the Lesser Antilles' only avian radiation. We improve on previous phylogenetic studies of this group by using double-digest restriction site-associated DNA sequencing (ddRAD-seq) to broadly sample loci scattered across the nuclear genome. A variety of analyses, based on either nucleotide variation in 2223 loci recovered in all samples or at 13,282 loci confidently scored as present or absent in all samples, converged on a single well-supported phylogenetic hypothesis. Results indicate that the resident West Indian taxa form a monophyletic group, exclusive of the Neotropical-Nearctic migratory Gray Catbird Dumetella carolinensis, which breeds in North America; this outcome differs from earlier studies suggesting that Gray Catbird was nested within a clade of island resident species. Thus, our findings imply a single colonization of the West Indies without the need to invoke a subsequent 'reverse colonization' of the mainland by West Indian taxa. Additionally, our study is the first to sample both endemic subspecies of the endangered White-breasted Thrasher Ramphocinclus brachyurus. We find that these subspecies have a long history of evolutionary independence with no evidence of gene flow, and are as genetically divergent from each other as other genera in the group. These findings support recognition of R. brachyurus (restricted to Martinique) and the Saint Lucia Thrasher R. sanctaeluciae as two distinct, single-island endemic species, and indicate the need to re-evaluate conservation plans for these taxa. Our results demonstrate the utility of phylogenomic datasets for generating robust systematic hypotheses.

ddRAD-seq phylogenetics based on nucleotide, indel, and presence-absence polymorphisms: Analyses of two avian genera with contrasting histories.

Genotype-by-sequencing (GBS) methods have revolutionized the field of molecular ecology, but their application in molecular phylogenetics remains somewhat limited. In addition, most phylogenetic studies based on large GBS data sets have relied on analyses of concatenated data rather than species tree methods that explicitly account for genealogical stochasticity among loci. We explored the utility of "double-digest" restriction site-associated DNA sequencing (ddRAD-seq) for phylogenetic analyses of the Lagonosticta firefinches (family Estrildidae) and the Vidua brood parasitic finches (family Viduidae). As expected, the number of homologous loci shared among samples was negatively correlated with genetic distance due to the accumulation of restriction site polymorphisms. Nonetheless, for each genus, we obtained data sets of ∼3000 loci shared in common among all samples, including a more distantly related outgroup taxon. For all samples combined, we obtained >1000 homologous loci despite ∼20my divergence between estrildid and parasitic finches. In addition to nucleotide polymorphisms, the ddRAD-seq data yielded large sets of indel and locus presence-absence polymorphisms, all of which had higher consistency indices than mtDNA sequence data in the context of concatenated parsimony analyses. Species tree methods, using individual gene trees or single nucleotide polymorphisms as input, generated results broadly consistent with analyses of concatenated data, particularly for Lagonosticta, which appears to have a well resolved, bifurcating history. Results for Vidua were also generally consistent across methods and data sets, although nodal support and results from different species tree methods were more variable. Lower gene tree congruence in Vidua is likely the result of its unique evolutionary history, which includes rapid speciation by host shift and occasional hybridization and introgression due to incomplete reproductive isolation. We conclude that ddRAD-seq is a cost-effective method for generating robust phylogenetic data sets, particularly for analyses of closely related species and genera.

De novo assembly of the dual transcriptomes of a polymorphic raptor species and its malarial parasite.

BACKGROUND: Studies of non-model species are important for understanding the molecular processes underpinning phenotypic variation under natural ecological conditions. The common buzzard (Buteo buteo; Aves: Accipitriformes) is a widespread and common Eurasian raptor with three distinct plumage morphs that differ in several fitness-related traits, including parasite infestation. To provide a genomic resource for plumage polymorphic birds in general and to search for candidate genes relating to fitness, we generated a transcriptome from a single dead buzzard specimen plus easily accessible, minimally invasive samples from live chicks. RESULTS: We not only de novo assembled a near-complete buzzard transcriptome, but also obtained a significant fraction of the transcriptome of its malaria-like parasite, Leucocytozoon buteonis. By identifying melanogenesis-related transcripts that are differentially expressed in light ventral and dark dorsal feathers, but which are also expressed in other regions of the body, we also identified a suite of candidate genes that could be associated with fitness differences among the morphs. These include several immune-related genes, providing a plausible link between melanisation and parasite load. qPCR analysis of a subset of these genes revealed significant differences between ventral and dorsal feathers and an additional effect of morph. CONCLUSION: This new resource provides preliminary insights into genes that could be involved in fitness differences between the buzzard colour morphs, and should facilitate future studies of raptors and their malaria-like parasites.

Population genetic structure of a common host predicts the spread of white-nose syndrome, an emerging infectious disease in bats.

Landscape complexity influences patterns of animal dispersal, which in turn may affect both gene flow and the spread of pathogens. White-nose syndrome (WNS) is an introduced fungal disease that has spread rapidly throughout eastern North America, causing massive mortality in bat populations. We tested for a relationship between the population genetic structure of the most common host, the little brown myotis (Myotis lucifugus), and the geographic spread of WNS to date by evaluating logistic regression models of WNS risk among hibernating colonies in eastern North America. We hypothesized that risk of WNS to susceptible host colonies should increase with both geographic proximity and genetic similarity, reflecting historical connectivity, to infected colonies. Consistent with this hypothesis, inclusion of genetic distance between infected and susceptible colonies significantly improved models of disease spread, capturing heterogeneity in the spatial expansion of WNS despite low levels of genetic differentiation among eastern populations. Expanding our genetic analysis to the continental range of little brown myotis reveals strongly contrasting patterns of population structure between eastern and western North America. Genetic structure increases markedly moving westward into the northern Great Plains, beyond the current distribution of WNS. In western North America, genetic differentiation of geographically proximate populations often exceeds levels observed across the entire eastern region, suggesting infrequent and/or locally restricted dispersal, and thus relatively limited opportunities for pathogen introduction in western North America. Taken together, our analyses suggest a possibly slower future rate of spread of the WNS pathogen, at least as mediated by little brown myotis.

Speciation genomics and a role for the Z chromosome in the early stages of divergence between Mexican ducks and mallards.

Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and Mexican duck (A. [p.] diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z-linked sex chromosome loci. Between the two taxa, Z-linked loci (ΦST  = 0.088) were 5.2 times more differentiated than autosomal DNA (ΦST  = 0.017) but comparable to mtDNA (ΦST  = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z-linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (ΦST ) and absolute (dXY ) estimates of divergence. In contrast, the ratio of Z-linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (ΦST  = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.

Rapid diversification associated with ecological specialization in Neotropical Adelpha butterflies.

Rapid diversification is often associated with morphological or ecological adaptations that allow organisms to radiate into novel niches. Neotropical Adelpha butterflies, which comprise over 200 species and subspecies, are characterized by extraordinary breadth in host plant use and wing colour patterns compared to their closest relatives. To examine the relationship between phenotypic and species diversification, we reconstructed the phylogenetic history of Adelpha and its temperate sister genus Limenitis using genomewide restriction-site-associated DNA (RAD) sequencing. Despite a declining fraction of shared markers with increasing evolutionary distance, the RAD-Seq data consistently generated well-supported trees using a variety of phylogenetic methods. These well-resolved phylogenies allow the identification of an ecologically important relationship with a toxic host plant family, as well as the confirmation of widespread, convergent wing pattern mimicry throughout the genus. Taken together, our results support the hypothesis that evolutionary innovations in both larvae and adults have permitted the colonization of novel host plants and fuelled adaptive diversification within this large butterfly radiation.

Ancient host specificity within a single species of brood parasitic bird.

Parasites that exploit multiple hosts often experience diversifying selection for host-specific adaptations. This can result in multiple strains of host specialists coexisting within a single parasitic species. A long-standing conundrum is how such sympatric host races can be maintained within a single parasitic species in the face of interbreeding among conspecifics specializing on different hosts. Striking examples are seen in certain avian brood parasites such as cuckoos, many of which show host-specific differentiation in traits such as host egg mimicry. Exploiting a Zambian egg collection amassed over several decades and supplemented by recent fieldwork, we show that the brood parasitic Greater Honeyguide Indicator indicator exhibits host-specific differentiation in both egg size and egg shape. Genetic analysis of honeyguide eggs and chicks show that two highly divergent mitochondrial DNA lineages are associated with ground- and tree-nesting hosts, respectively, indicating perfect fidelity to two mutually exclusive sets of host species for millions of years. Despite their age and apparent adaptive diversification, however, these ancient lineages are not cryptic species; a complete lack of differentiation in nuclear genes shows that mating between individuals reared by different hosts is sufficiently frequent to prevent speciation. These results indicate that host specificity is maternally inherited, that host-specific adaptation among conspecifics can be maintained without reproductive isolation, and that host specificity can be remarkably ancient in evolutionary terms.

A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat.

One of the most widely distributed bats in the New World, the big brown bat (Eptesicus fuscus) exhibits well-documented geographic variation in morphology and life history traits, suggesting the potential for significant phylogeographic structure as well as adaptive differentiation among populations. In a pattern broadly consistent with morphologically defined subspecies, we found deeply divergent mitochondrial lineages restricted to different geographic regions. In contrast, sequence data from two nuclear loci suggest a general lack of regional genetic structure except for peripheral populations in the Caribbean and Mexico/South America. Coalescent analyses suggest that the striking difference in population structure between genomes cannot be attributed solely to different rates of lineage sorting, but is likely due to male-mediated gene flow homogenizing nuclear genetic diversity across most of the continental range. Despite this ongoing gene flow, selection has apparently been effective in producing and maintaining adaptive differentiation among populations, while strong female site fidelity, maintained over the course of millions of years, has produced remarkably deep divergence among geographically isolated matrilines. Our results highlight the importance of evaluating multiple genetic markers for a more complete understanding of population structure and history.

Leveraging genomics to understand threats to migratory birds.

Understanding how risk factors affect populations across their annual cycle is a major challenge for conserving migratory birds. For example, disease outbreaks may happen on the breeding grounds, the wintering grounds, or during migration and are expected to accelerate under climate change. The ability to identify the geographic origins of impacted individuals, especially outside of breeding areas, might make it possible to predict demographic trends and inform conservation decision-making. However, such an effort is made more challenging by the degraded state of carcasses and resulting low quality of DNA available. Here, we describe a rapid and low-cost approach for identifying the origins of birds sampled across their annual cycle that is robust even when DNA quality is poor. We illustrate the approach in the common loon (Gavia immer), an iconic migratory aquatic bird that is under increasing threat on both its breeding and wintering areas. Using 300 samples collected from across the breeding range, we develop a panel of 158 single-nucleotide polymorphisms (SNP) loci with divergent allele frequencies across six genetic subpopulations. We use this SNP panel to identify the breeding grounds for 142 live nonbreeding individuals and carcasses. For example, genetic assignment of loons sampled during botulism outbreaks in parts of the Great Lakes provides evidence for the significant role the lakes play as migratory stopover areas for loons that breed across wide swaths of Canada, and highlights the vulnerability of a large segment of the breeding population to botulism outbreaks that are occurring in the Great Lakes with increasing frequency. Our results illustrate that the use of SNP panels to identify breeding origins of carcasses collected during the nonbreeding season can improve our understanding of the population-specific impacts of mortality from disease and anthropogenic stressors, ultimately allowing more effective management.

Sexual imprinting misguides species recognition in a facultative interspecific brood parasite.

Sexual reproduction relies on the recognition of conspecifics for breeding. Most experiments in birds have implicated a critical role for early social learning in directing subsequent courtship behaviours and mating decisions. This classical view of avian sexual imprinting is challenged, however, by studies of megapodes and obligate brood parasites, species in which reliable recognition is achieved despite the lack of early experience with conspecifics. By rearing males with either conspecific or heterospecific brood mates, we experimentally tested the effect of early social experience on the association preferences and courtship behaviours of two sympatrically breeding ducks. We predicted that redheads (Aythya americana), which are facultative interspecific brood parasites, would show a diminished effect of early social environment on subsequent courtship preferences when compared with their host and congener, the canvasback (Aythya valisineria). Contrary to expectations, cross-fostered males of both species courted heterospecific females and preferred them in spatial association tests, whereas control males courted and associated with conspecific females. These results imply that ontogenetic constraints on species recognition may be a general impediment to the initial evolution of interspecific brood parasitism in birds. Under more natural conditions, a variety of mechanisms may mitigate or counteract the effects of early imprinting for redheads reared in canvasback broods.

Phylogenetic and structural analysis of the HbA (alphaA/betaA) and HbD (alphaD/betaA) hemoglobin genes in two high-altitude waterfowl from the Himalayas and the Andes: Bar-headed goose (Anser indicus) and Andean goose (Chloephaga melanoptera).

Two species of waterfowl living at high altitude provide a prominent example of parallel adaptation at the molecular level. The bar-headed goose (Anser indicus) breeds at high elevations in central Asia and migrates across the Himalayas, where the partial pressure of oxygen (O(2)) is one-third of sea level. In South America, the distantly related Andean goose (Chloephaga melanoptera) is endemic to the high Andes. Both species exhibit increased blood-O(2) affinity, which has been attributed to the effects of single amino acid substitutions in the major hemoglobin. Here we present phylogenetic analyses of the swans and geese (Anserinae) and South American sheldgeese (Anatinae) using the three genes that encode the major (HbA) and minor (HbD) hemoglobin isoforms. We sought to determine whether two amino acid substitutions that have been the focus of extensive biochemical analysis (Ala-alpha(A)119 and Ser-beta(A)55) are uniquely derived in bar-headed goose and Andean goose, respectively, and to examine evidence of molecular adaptation at other positions in hemoglobin genes by comparing these two high-altitude taxa to their closest relatives. Bayesian analysis of the alpha(A)-, alpha(D)-, and beta(A)-subunit genes produced well-resolved phylogenies, with high posterior probabilities and bootstrap values for most genera. The bar-headed goose is likely sister to all other Anser species. Andean goose, the sole highland representative of the South American sheldgeese is either sister to the other Chloephaga species or sister to Neochen. In the bar-headed goose, four derived substitutions were observed in HbA (alpha(A)12, 18, 63, 119) and two in HbD (alpha(D)2, 47). Four derived substitutions in Andean goose include three in HbA (alpha(A)8, 77; beta(A)86) and two in HbD (alpha(D)9; beta(A)86). Considering both highland species, four substitutions (Ala-alpha(A)8, Ala-alpha(A)12, Ser-alpha(A)18, Leu-alpha(D)9) were located at adjacent positions on the A helix (or AB corner) of the alpha-chains, three others (Thr-alpha(A)77, Ser-beta(A)86, Ser-alpha(D)2) were in close proximity to inositolpentaphosphate (IP(5)) binding sites, and Ala-alpha(A)119 occurred at an alphabeta intersubunit contact. Ser-beta(A)55, which is involved in the same alphabeta intersubunit contact and was previously shown to increase Hb-O(2) affinity, is not unique to Andean goose, but is a synapomorphy of the South American sheldgeese, a clade of predominantly lowland waterfowl. Our findings illustrate the importance of understanding phylogenetic relationships and polarity of character-state changes when making inferences about adaptive evolution.

Speciation by host switch in brood parasitic indigobirds.

A growing body of empirical and theoretical work supports the plausibility of sympatric speciation, but there remain few examples in which all the essential components of the process are well understood. The African indigobirds Vidua spp. are host-specific brood parasites. Indigobird nestlings are reared along with host young, and mimic the mouth markings of their respective hosts. As adults, male indigobirds mimic host song, whereas females use these songs to choose both their mates and the nests they parasitize. These behavioural mechanisms promote the cohesion of indigobird populations associated with a given host species, and provide a mechanism for reproductive isolation after a new host is colonized. Here we show that all indigobird species are similar genetically, but are significantly differentiated in both mitochondrial haplotype and nuclear allele frequencies. These data support a model of recent sympatric speciation. In contrast to the cuckoo Cuculus canorus, in which only female lineages are faithful to specific hosts, host switches have led to speciation in indigobirds because both males and females imprint on their hosts.

Does coevolution promote species richness in parasitic cuckoos?

Why some lineages have diversified into larger numbers of species than others is a fundamental but still relatively poorly understood aspect of the evolutionary process. Coevolution has been recognized as a potentially important engine of speciation, but has rarely been tested in a comparative framework. We use a comparative approach based on a complete phylogeny of all living cuckoos to test whether parasite-host coevolution is associated with patterns of cuckoo species richness. There are no clear differences between parental and parasitic cuckoos in the number of species per genus. However, a cladogenesis test shows that brood parasitism is associated with both significantly higher speciation and extinction rates. Furthermore, subspecies diversification rate estimates were over twice as high in parasitic cuckoos as in parental cuckoos. Among parasitic cuckoos, there is marked variation in the severity of the detrimental effects on host fitness; chicks of some cuckoo species are raised alongside the young of the host and others are more virulent, with the cuckoo chick ejecting or killing the eggs/young of the host. We show that cuckoos with a more virulent parasitic strategy have more recognized subspecies. In addition, cuckoo species with more recognized subspecies have more hosts. These results hold after controlling for confounding geographical effects such as range size and isolation in archipelagos. Although the power of our analyses is limited by the fact that brood parasitism evolved independently only three times in cuckoos, our results suggest that coevolutionary arms races with hosts have contributed to higher speciation and extinction rates in parasitic cuckoos.

Incomplete reproductive isolation following host shift in brood parasitic indigobirds.

Behavioural and molecular studies suggest that brood parasitic indigobirds (Vidua spp.) rapidly diversified through a process of speciation by host shift. However, behavioural imprinting on host song, the key mechanism promoting speciation in this system, may also lead to hybridization and gene flow among established indigobird species when and if female indigobirds parasitize hosts already associated with other indigobird species. It is therefore not clear to what extent the low level of genetic differentiation among indigobird species is due to recent common ancestry versus ongoing gene flow. We tested for reproductive isolation among three indigobird species in Cameroon, one of which comprises two morphologically indistinguishable host races. Mimicry of host songs corresponded with plumage colour in 184 male indigobirds, suggesting that females rarely parasitize the host of another indigobird species. Paternity analyses, however, suggest that imperfect specificity in host and/or mate choice allows for continuing gene flow between recently formed host races of the Cameroon Indigobird Vidua camerunensis; while 63 pairs of close relatives were associated with the same host, two strongly supported father-son pairs included males mimicking the songs of the two different hosts of V. camerunensis. Thus, complete reproductive isolation is not necessarily an automatic consequence of host shifts, a result that suggests an important role for natural and/or sexual selection in indigobird speciation.