Haldane's rule in an avian system: using cline theory and divergence population genetics to test for differential introgression of mitochondrial, autosomal, and sex-linked loci across the Passerina bunting hybrid zone

Using cline fitting and divergence population genetics, we tested a prediction of Haldane's rule: autosomal alleles should introgress more than z-linked alleles or mitochondrial haplotypes across the Passerina amoena/Passerina cyanea (Aves: Cardinalidae) hybrid zone. We screened 222 individuals collected along a transect in the Great Plains of North America that spans the contact zone for mitochondrial (two genes), autosomal (four loci) and z-linked (two loci) markers. Maximum-likelihood cline widths estimated from the mitochondrial (223 km) and z-linked (309 km) datasets were significantly narrower on average than the autosomal cline widths (466 km). We also found that mean coalescent-based estimates of introgression were larger for the autosomal loci (0.63 genes/generation, scaled to the mutation rate mu) than for both the mitochondrial (0.27) and z-linked loci (0.59). These patterns are consistent with Haldane's rule, but the among-locus variation also suggests many independently segregating loci are required to investigate introgression patterns across the genome. These results provide the first comprehensive comparison of mitochondrial, sex-linked, and autosomal loci across an avian hybrid zone and add to the body of evidence suggesting that sex chromosomes play an important role in the formation and maintenance of reproductive isolation between closely related species.     

Behavioural response to song and genetic divergence in two subspecies of white‐crowned sparrows (Zonotrichia leucophrys)

Divergence in sexual signals may drive reproductive isolation between lineages, but behavioural barriers can weaken in contact zones. Here, we investigate the role of song as a behavioural and genetic barrier in a contact zone between two subspecies of white‐crowned sparrows (Zonotrichia leucophrys). We employed a reduced genomic data set to assess population structure and infer the history underlying divergence, gene flow and hybridization. We also measured divergence in song and tested behavioural responses to song using playback experiments within and outside the contact zone. We found that the subspecies form distinct genetic clusters, and demographic inference supported a model of secondary contact. Song phenotype, particularly length of the first note (a whistle), was a significant predictor of genetic subspecies identity and genetic distance along the hybrid zone, suggesting a close link between song and genetic divergence in this system. Individuals from both parental and admixed localities responded significantly more strongly to their own song than to the other subspecies song, supporting song as a behavioural barrier. Putative parental and admixed individuals were not significantly different in their strength of discrimination between own and other songs; however, individuals from admixed localities tended to discriminate less strongly, and this difference in discrimination strength was explained by song dissimilarity as well as genetic distance. Therefore, we find that song acts as a reproductive isolating mechanism that is potentially weakening in a contact zone between the subspecies. Our findings also support the hypothesis that intraspecific song variation can reduce gene flow between populations.     

Niche evolution and diversification in a Neotropical radiation of birds (Aves: Furnariidae)

Rapid diversification may be caused by ecological adaptive radiation via niche divergence. In this model, speciation is coupled with niche divergence and lineage diversification is predicted to be correlated with rates of niche evolution. Studies of the role of niche evolution in diversification have generally focused on ecomorphological diversification but climatic‐niche evolution may also be important. We tested these alternatives using a phylogeny of 298 species of ovenbirds (Aves: Furnariidae). We found that within Furnariidae, variation in species richness and diversification rates of subclades were best predicted by rate of climatic‐niche evolution than ecomorphological evolution. Although both are clearly important, univariate regression and multivariate model averaging more consistently supported the climatic‐niche as the best predictor of lineage diversification. Our study adds to the growing body of evidence, suggesting that climatic‐niche divergence may be an important driver of rapid diversification in addition to ecomorphological evolution. However, this pattern may depend on the phylogenetic scale at which rate heterogeneity is examined.     

Comparison of species tree methods for reconstructing the phylogeny of bearded manakins (Aves: Pipridae, Manacus) from multilocus sequence data

Although the power of multi-locus data in estimating species trees is apparent, it is also clear that the analytical methodologies for doing so are still maturing. For example, of the methods currently available for estimating species trees from multilocus data, the Bayesian method introduced by Liu and Pearl (2007; BEST) is the only one that provides nodal support values. Using gene sequences from five nuclear loci, we explored two analytical methods (deep coalescence and BEST) to reconstruct the species tree of the five primary Manacus OTUs: M. aurantiacus, M. candei, M. vitellinus, populations of M. manacus from west of the Andes (M. manacus (w)), and populations of M. manacus from east of the Andes (M. manacus (e)). Both BEST and deep coalescence supported a sister relationship between M. vitellinus and M. manacus (w). A lower probability tree from the BEST analysis and one of the most parsimonious deep coalescence trees also supported a sister relationship between M. candei and M. aurantiacus. Because hybrid zones connect the distributions of most Manacus species, we examined the potential influence of post-divergence gene flow on the sister relationship of parapatrically distributed M. vitellinus and M. manacus (w). An isolation-with-migration (IM) analysis found relatively high levels of gene flow between M. vitellinus and M. manacus (w). Whether the gene flow is obscuring a true sister relationship between M. manacus (w) and M. manacus (e) remained unclear, pointing to the need for more detailed models accommodating multispecies, multilocus DNA sequence data.     

Sequence variation in the coding region of the melanocortin-1 receptor gene (MC1R) is not associated with plumage variation in the blue-crowned manakin (Lepidothrix coronata)

Avian plumage traits are the targets of both natural and sexual selection. Consequently, genetic changes resulting in plumage variation among closely related taxa might represent important evolutionary events. The molecular basis of such differences, however, is unknown in most cases. Sequence variation in the melanocortin-1 receptor gene (MC1R) is associated with melanistic phenotypes in many vertebrate taxa, including several avian species. The blue-crowned manakin (Lepidothrix coronata), a widespread, sexually dichromatic passerine, exhibits striking geographic variation in male plumage colour across its range in southern Central America and western Amazonia. Northern males are black with brilliant blue crowns whereas southern males are green with lighter blue crowns. We sequenced 810 bp of the MC1R coding region in 23 individuals spanning the range of male plumage variation. The only variable sites we detected among L. coronata sequences were four synonymous substitutions, none of which were strictly associated with either plumage type. Similarly, comparative analyses showed that L. coronata sequences were monomorphic at the three amino acid sites hypothesized to be functionally important in other birds. These results demonstrate that genes other than MC1R underlie melanic plumage polymorphism in blue-crowned manakins.     

hzar: hybrid zone analysis using an R software package

We present a new software package (hzar ) that provides functions for fitting molecular genetic and morphological data from hybrid zones to classic equilibrium cline models using the Metropolis–Hastings Markov chain Monte Carlo (MCMC) algorithm. The software applies likelihood functions appropriate for different types of data, including diploid and haploid genetic markers and quantitative morphological traits. The modular design allows flexibility in fitting cline models of varying complexity. To facilitate hypothesis testing, an autofit function is included that allows automated model selection from a set of nested cline models. Cline parameter values, such as cline centre and cline width, are estimated and may be compared statistically across clines. The package is written in the R language and is available through the Comprehensive R Archive Network (CRAN; http://cran.r-project.org/ ). Here, we describe hzar and demonstrate its use with a sample data set from a well‐studied hybrid zone in western Panama between white‐collared (Manacus candei) and golden‐collared manakins (M. vitellinus). Comparisons of our results with previously published results for this hybrid zone validate the hzar software. We extend analysis of this hybrid zone by fitting additional models to molecular data where appropriate.