Recombination rate variation shapes genomic variability of phylogeographic structure in a widespread North American songbird (Aves: Certhia americana).

The nonrandom distribution of chromosomal characteristics and functional elements-genomic architecture-impacts the relative strengths and impacts of population genetic processes across the genome. Due to this relationship, genomic architecture has the potential to shape variation in population genetic structure across the genome. Population genetic structure has been shown to vary across the genome in a variety of taxa, but this body of work has largely focused on pairwise population genomic comparisons between closely related taxa. Here, we used whole genome sequencing of seven phylogeographically structured populations of a North American songbird, the Brown Creeper (Certhia americana), to determine the impacts of genomic architecture on phylogeographic structure variation across the genome. Using multiple methods to infer phylogeographic structure-ordination, clustering, and phylogenetic methods-we found that recombination rate variation explained a large proportion of phylogeographic structure variation. Genomic regions with low recombination showed phylogeographic structure consistent with the genome-wide pattern. In regions with high recombination, we found strong phylogeographic structure, but with discordant patterns relative to the genome-wide pattern. In regions with high recombination rate, we found that populations with small effective population sizes evolve relatively more rapidly than larger populations, leading to discordant signatures of phylogeographic structure. These results suggest that the interplay between recombination rate variation and effective population sizes shape the relative impacts of selection and genetic drift in different parts of the genome. Overall, the combined interactions of population genetic processes, genomic architecture, and effective population sizes shape patterns of variability in phylogeographic structure across the genome of the Brown Creeper.

Genomic divergence and introgression between cryptic species of a widespread North American songbird.

Analysis of genomic variation among related populations can sometimes reveal distinct species that were previously undescribed due to similar morphological appearances, and close examination of such cases can provide much insight regarding speciation. Genomic data can also reveal the role of reticulate evolution in differentiation and speciation. White-breasted nuthatches (Sitta carolinensis) are widely distributed North American songbirds that are currently classified as a single species but have been suspected to represent a case of cryptic speciation. Previous genetic analyses suggested four divergent groups, but it was unclear whether these represented multiple reproductively isolated species. Using extensive genomic sampling of over 350 white-breasted nuthatches from across North America and a new chromosome-level reference genome, we asked if white-breasted nuthatches are comprised of multiple species and whether introgression has occurred between divergent populations. Genomic variation of over 300,000 loci revealed four highly differentiated populations (Pacific, n = 45; Eastern, n = 23; Rocky Mountains North, n = 138; and Rocky Mountains South, n = 150) with geographic ranges that are adjacent. We observed a moderate degree of admixture between Rocky Mountain populations but only a small number of hybrids between the Rockies and the Eastern population. The rarity of hybrids together with high levels of differentiation between populations is supportive of populations having some level of reproductive isolation. Between populations, we show evidence for introgression from a divergent ghost lineage of white-breasted nuthatches into the Rocky Mountains South population, which is otherwise closely related to Rocky Mountains North. We conclude that white-breasted nuthatches are best considered at least three species and that ghost lineage introgression has contributed to differentiation between the two Rocky Mountain populations. White-breasted nuthatches provide a dramatic case of morphological similarity despite high genomic differentiation, and the varying levels of reproductive isolation among the four groups provide an example of the speciation continuum.

Estimating phylogenies from genomes: A beginners review of commonly used genomic data in vertebrate phylogenomics.

Despite the increasing feasibility of sequencing whole genomes from diverse taxa, a persistent problem in phylogenomics is the selection of appropriate genetic markers or loci for a given taxonomic group or research question. In this review, we aim to streamline the decision-making process when selecting specific markers to use in phylogenomic studies by introducing commonly used types of genomic markers, their evolutionary characteristics, and their associated uses in phylogenomics. Specifically, we review the utilities of ultraconserved elements (including flanking regions), anchored hybrid enrichment loci, conserved nonexonic elements, untranslated regions, introns, exons, mitochondrial DNA, single nucleotide polymorphisms, and anonymous regions (nonspecific regions that are evenly or randomly distributed across the genome). These various genomic elements and regions differ in their substitution rates, likelihood of neutrality or of being strongly linked to loci under selection, and mode of inheritance, each of which are important considerations in phylogenomic reconstruction. These features may give each type of marker important advantages and disadvantages depending on the biological question, number of taxa sampled, evolutionary timescale, cost effectiveness, and analytical methods used. We provide a concise outline as a resource to efficiently consider key aspects of each type of genetic marker. There are many factors to consider when designing phylogenomic studies, and this review may serve as a primer when weighing options between multiple potential phylogenomic markers.

Phylogenomic Data Reveal Widespread Introgression Across the Range of an Alpine and Arctic Specialist.

Understanding how gene flow affects population divergence and speciation remains challenging. Differentiating one evolutionary process from another can be difficult because multiple processes can produce similar patterns, and more than one process can occur simultaneously. Although simple population models produce predictable results, how these processes balance in taxa with patchy distributions and complicated natural histories is less certain. These types of populations might be highly connected through migration (gene flow), but can experience stronger effects of genetic drift and inbreeding, or localized selection. Although different signals can be difficult to separate, the application of high-throughput sequence data can provide the resolution necessary to distinguish many of these processes. We present whole-genome sequence data for an avian species group with an alpine and arctic tundra distribution to examine the role that different population genetic processes have played in their evolutionary history. Rosy-finches inhabit high elevation mountaintop sky islands and high-latitude island and continental tundra. They exhibit extensive plumage variation coupled with low levels of genetic variation. Additionally, the number of species within the complex is debated, making them excellent for studying the forces involved in the process of diversification, as well as an important species group in which to investigate species boundaries. Total genomic variation suggests a broadly continuous pattern of allele frequency changes across the mainland taxa of this group in North America. However, phylogenomic analyses recover multiple distinct, well supported, groups that coincide with previously described morphological variation and current species-level taxonomy. Tests of introgression using D-statistics and approximate Bayesian computation reveal significant levels of introgression between multiple North American taxa. These results provide insight into the balance between divergent and homogenizing population genetic processes and highlight remaining challenges in interpreting conflict between different types of analytical approaches with whole-genome sequence data. [ABBA-BABA; approximate Bayesian computation; gene flow; phylogenomics; speciation; whole-genome sequencing.].

Dense Geographic and Genomic Sampling Reveals Paraphyly and a Cryptic Lineage in a Classic Sibling Species Complex.

Incomplete or geographically biased sampling poses significant problems for research in phylogeography, population genetics, phylogenetics, and species delimitation. Despite the power of using genome-wide genetic markers in systematics and related fields, approaches such as the multispecies coalescent remain unable to easily account for unsampled lineages. The Empidonax difficilis/Empidonax occidentalis complex of small tyrannid flycatchers (Aves: Tyrannidae) is a classic example of widely distributed species with limited phenotypic geographic variation that was broken into two largely cryptic (or "sibling") lineages following extensive study. Though the group is well-characterized north of the US Mexico border, the evolutionary distinctiveness and phylogenetic relationships of southern populations remain obscure. In this article, we use dense genomic and geographic sampling across the majority of the range of the E. difficilis/E. occidentalis complex to assess whether current taxonomy and species limits reflect underlying evolutionary patterns, or whether they are an artifact of historically biased or incomplete sampling. We find that additional samples from Mexico render the widely recognized species-level lineage E. occidentalis paraphyletic, though it retains support in the best-fit species delimitation model from clustering analyses. We further identify a highly divergent unrecognized lineage in a previously unsampled portion of the group's range, which a cline analysis suggests is more reproductively isolated than the currently recognized species E. difficilis and E. occidentalis. Our phylogeny supports a southern origin of these taxa. Our results highlight the pervasive impacts of biased geographic sampling, even in well-studied vertebrate groups like birds, and illustrate what is a common problem when attempting to define species in the face of recent divergence and reticulate evolution.

Extensive gene flow characterizes the phylogeography of a North American migrant bird: Black-headed Grosbeak (Pheucticus melanocephalus).

We describe range-wide phylogeographic variation in the Black-headed Grosbeak (Pheucticus melanocephalus), a songbird that is widely distributed across North American scrublands and forests. Phylogenetic analysis of mitochondrial DNA (mtDNA, n=424) revealed three geographically structured clades. One widespread clade occurs throughout the Rocky Mountains, Great Basin, and Mexican Plateau, a second clade is found on the Pacific coast and in coastal ranges; and, a third in the Sierra Madre del Sur of Oaxaca and Guerrero. Some geographical structuring occurs in Mexican Plateau and Sierra Madre Oriental mtDNA clade, presumably because these populations have been more stable over time than northern populations. Multiple mitochondrial groups are found sympatrically in the Okanogan River Valley in Washington, the eastern Sierra Nevada, and the Transvolcanic Belt across central Mexico, indicating that there is a potential for introgression. Analyses of 12 nuclear loci did not recover the same geographically structured clades. Population analyses show high levels of gene flow in nucDNA from the Interior into the Sierra Madre del Sur and Pacific population groups, possibly indicating expansion of the Interior population at the expense of peripheral populations.

Comparative morphology of the avian bony columella.

In birds, the columella is the only bony element of the sound conducting apparatus, conveying vibrations of the cartilaginous extracolumella to the fluid of the inner ear. Although avian columellar morphology has attracted some attention over the past century, it nonetheless remains poorly described in the literature. The few existing studies mostly focus on morphological descriptions in relatively few taxa, with no taxonomically broad surveys yet published. Here we use observations of columellae from 401 extant bird species to provide a comprehensive survey of columellar morphology in a phylogenetic context. We describe the columellae of several taxa for the first time and identify derived morphologies characterizing higher-level clades based on current phylogenies. In particular, we identify a derived columellar morphology diagnosing a major subclade of Accipitridae. Within Suliformes, we find that Fregatidae, Sulidae, and Phalacrocoracidae share a derived morphology that is absent in Anhingidae, suggesting a secondary reversal. Phylogenetically informed comparisons allow recognition of instances of homoplasy, including the distinctive bulbous columellae in suboscine passerines and taxa belonging to Eucavitaves, and bulging footplates that appear to have evolved at least twice independently in Strigiformes. We consider phylogenetic and functional factors influencing avian columellar morphology, finding that aquatic birds possess small footplates relative to columellar length, possibly related to hearing function in aquatic habitats. By contrast, the functional significance of the distinctive bulbous basal ends of the columellae of certain arboreal landbird taxa remains elusive.

Multilocus phylogeny and biogeography of the New World Pheucticus grosbeaks (Aves: Cardinalidae).

Using a multilocus approach, we investigated the tempo and pattern of diversification in a widely distributed New World songbird, the cardinalid genus Pheucticus. Each of the three geographic groups recovered (North American, Middle American, and South American) was comprised of a pair of currently recognized species, and four, three, and three geographically and genetically distinct phylogeographic lineages respectively. Diversification within Pheucticus appears to have occurred at a relatively constant pace throughout the Pleistocene and evenly across a broad latitudinal distribution. The Isthmus of Panama completion and Pleistocene glacial cycles both appear to have played prominent roles in the diversification of this group.

Survey of haemosporidian parasite infections in an endangered high alpine bird.

The Brown-capped Rosy-Finch is an endangered high alpine specialist that breeds on cliffs in the Rocky Mountains of North America. We know little about the parasites that infect these birds but filling this knowledge gap is important for understanding their population decline. The aim of this study was to survey haemosporidian parasite infections in Brown-capped Rosy-Finches. We sampled 104 Brown-capped Rosy-Finches during their breeding season at six sampling sites spanning the Colorado Rocky Mountains where they are nearly endemic. We used nested PCR methods to screen birds for parasite infections, and Sanger sequencing data were used to identify parasite lineages. Four of the sampling sites had birds with infections. Females were more often infected than males (prevalence = 15.6% and 9.7%, respectively). We observed twice as many infected birds sampled in July compared to birds sampled in June (prevalence = 12.5% and 6.25%, respectively). The percent of infected birds by sampling site ranged from 0 to over 20%. In total we identified eight different genetic lineages of haemosporidian parasites infecting Brown-capped Rosy-Finches-seven were in the genus Leucocytozoon, and one was in the genus Haemoproteus. Network analysis clearly separates Haemoproteus from Leucocytozoon, with Leucocytozoon lineages comprising two major clusters. Based on reports made to the MalAvi database, all lineages in Cluster 1 have the same reported insect vector, Simulium silvestre. We report that Brown-capped Rosy-Finches experience infections with haemosporidian parasites and that the percentage of infected birds differed markedly between sampling sites. We hypothesize that vector ecology and associated variation in climate at sampling sites explain our observations.

The genetic basis of plumage coloration and elevation adaptation in a clade of recently diverged alpine and arctic songbirds.

Trait genetic architecture plays an important role in the probability that variation in that trait leads to divergence and speciation. In some cases, speciation may be driven by the generation of novel phenotypes through the recombination of genes associated with traits that are important for local adaptation or sexual selection. Here, we investigate the genetic basis of three plumage color traits, and one ecological trait, breeding elevation, in a recent avian radiation, the North American rosy-finches (Leucosticte spp.). We identify unique genomic regions associated with each trait and highlight 11 candidate genes. Among these are well-characterized melanogenesis genes, including Mitf and Tyrp1, and previously reported hypoxia-related genes including Egln1. Additionally, we use mitochondrial data to date the divergence of rosy-finch clades which appear to have diverged within the past 250 ky. Given the low levels of genome-wide differentiation among rosy-finch taxa, and evidence for extensive introgression in North America, plumage coloration and adaptation to high elevations have likely played large roles in generating the observed patterns of lineage divergence. The relative independence of these candidate regions across the genome suggests that recombination might have led to multiple phenotypes, and subsequent rosy-finch speciation, over short periods of time.

Speciation in the White-breasted Nuthatch (Sitta carolinensis): a multilocus perspective.

Inferring the evolutionary and ecological processes that have shaped contemporary species distributions using the geographic distribution of gene lineages is the principal goal of phylogeographic research. Researchers in the field have recognized that inferences made from a single gene, often mitochondrial, can be informative regarding the pattern of diversification but lack conclusive information regarding the evolutionary mechanisms that led to the observed patterns. Here, we use a multilocus (20 loci) data set to explore the evolutionary history of the White-breasted Nuthatch (Sitta carolinensis). A previous single-locus study found S. carolinensis is comprised of four reciprocally monophyletic clades geographically restricted to the pine and oak forests of: (i) eastern North America, (ii) southern Rocky Mountain and Mexican Mountain ranges, (iii) Eastern Sierra Nevada and Northern Rocky Mountains and (iv) Pacific slope of North America. The diversification of the clades was attributed to the fragmentation of North American pine and oak woodlands in the Pliocene with subsequent divergences owing to the Pleistocene glacial cycles. Principal component, clustering and species tree analyses of the multilocus data resolved the same four groups or lineages found in the single-locus study. Coalescent analyses and hypothesis testing of nested isolation and migration models indicate that isolation and not gene flow has been the major evolutionary mechanism responsible for shaping genetic variation, and all the divergence events within S. carolinensis have occurred in response to the Pleistocene glacial cycles.

Isolation-driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae).

Lineage, or true 'species', trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene-tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single-locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re-evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent-based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.

Cryptic diversity in a widespread North American songbird: phylogeography of the Brown Creeper (Certhia americana).

The identification of species via morphological characteristics has traditionally left cryptic species undescribed in taxa under selection for morphological conservation (or a lack of selection for morphological change). Treecreepers (Genus: Certhia) have a conserved morphological appearance, making it difficult to ascertain relationships in the genus based on morphology alone. Recent genetic and song structure studies of Eurasian Treecreepers identified cryptic species within Old World Certhia that were previously undescribed using morphological characteristics. Here, we use mtDNA to investigate cryptic diversity and patterns of diversification in the Brown Creeper (Certhia americana), the single described Certhia species in the Americas. Phylogenetic analyses identified six well-supported geographically-structured clades; the basal divergence separates a northern and a southern lineage in the Brown Creeper, likely cryptic species previously characterized as many subspecies. Sympatry is prevalent between clades in western North America, where possible contact zones warrant further investigation. Allopatry appears to be the primary driver of deep phylogeographic structure within the Brown Creeper; however, within clade diversity is highly correlated with the life history traits of the populations that comprise the geographically structured phylogroups.

The Genomic Signature of Allopatric Speciation in a Songbird Is Shaped by Genome Architecture (Aves: Certhia americana).

The genomic signature of speciation with gene flow is often attributed to the strength of divergent selection and recombination rate in regions harboring targets for selection. In contrast, allopatric speciation provides a different geographic context and evolutionary scenario, whereby introgression is limited by isolation rather than selection against gene flow. Lacking shared divergent selection or selection against hybridization, we would predict the genomic signature of allopatric speciation would largely be shaped by genomic architecture-the nonrandom distribution of functional elements and chromosomal characteristics-through its role in affecting the processes of selection and drift. Here, we built and annotated a chromosome-scale genome assembly for a songbird (Passeriformes: Certhia americana). We show that the genomic signature of allopatric speciation between its two primary lineages is largely shaped by genomic architecture. Regionally, gene density and recombination rate variation explain a large proportion of variance in genomic diversity, differentiation, and divergence. We identified a heterogeneous landscape of selection and neutrality, with a large portion of the genome under the effects of indirect selection. We found higher proportions of small chromosomes under the effects of indirect selection, likely because they have relatively higher gene density. At the chromosome scale, differential genomic architecture of macro- and microchromosomes shapes the genomic signatures of speciation: chromosome size has: 1) a positive relationship with genetic differentiation, genetic divergence, rate of lineage sorting in the contact zone, and proportion neutral evolution and 2) a negative relationship with genetic diversity and recombination rate.

Scaling of the avian middle ear.

This article presents a comparative study of morphology of the avian middle ear. The general morphology of the columella shows considerable variation across species, yet few studies have attempted to provide quantitative comparisons, and basic anatomical data has not been thoroughly reported. In this study, we examined the middle ear in 49 taxonomically diverse species of bird. We found significant correlations between measurements of several features (columellar length, mass, tympanic membrane area, footplate area) and interaural diameter. While scaling of columellar length with interaural diameter is consistent with isometry, masses and areas showed negative allometry, or a non-proportional scaling with interaural diameter. These observations remained true even for species with unusual middle ear morphology, such as Alcedinidae (Kingfishers) in which the basal struts of the columella form a structure almost resembling a mammalian stapes, or Tytonidae (Barn Owls) which have a highly bulbous footplate. It therefore appears that allometry cannot help explain the morphological variation in the columella.

A comparative study of avian middle ear mechanics.

Sound and hearing play an important role in the lives of many birds, and studies have been published on the acoustic habitat of various species, as well as on various aspects of their hearing system. However the function of the middle ear remains largely unexplored, with existing studies focusing on either single species, or a very narrow range of species. In this article we report measurements of the middle ear transfer function in 39 taxonomically diverse avian species. We used laser vibrometry to measure the vibrations of the columellar footplate in response to tones played in each animal's ear canal, and calculated the middle ear transfer functions. Transfer functions varied substantially across species, but appeared to follow common trends. Comparisons between the peak frequency in the transfer function and length/mass of the columella reveal a correlation between the three, however statistical analysis suggests that columellar length is a primary indicator of the peak frequency. These results provide a broad survey of avian middle ear function, and the observed trends provide a method with which to begin to predict the response of single ossicle middle ear systems.

Sedentary songbirds maintain higher prevalence of haemosporidian parasite infections than migratory conspecifics during seasonal sympatry.

Long-distance migrations influence the physiology, behavior, and fitness of migratory animals throughout their annual cycles, and fundamentally alter their interactions with parasites. Several hypotheses relating migratory behavior to the likelihood of parasitism have entered the literature, making conflicting, testable predictions. To assess how migratory behavior of hosts is associated with parasitism, we compared haemosporidian parasite infections between two closely related populations of a common North American sparrow, the dark-eyed junco, that co-occur in shared habitats during the non-breeding season. One population is sedentary and winters and breeds in the Appalachian Mountains. The other population is migratory and is found in seasonal sympatry with the sedentary population from October through April, but then flies (≥ 900 km) northwards to breed. The populations were sampled in the wild on the shared montane habitat at the beginning of winter and again after confining them in a captive common environment until the spring. We found significantly higher prevalence of haemosporidian parasite infections in the sedentary population. Among infected juncos, we found no difference in parasite densities (parasitemias) between the sedentary and migrant populations and no evidence for winter dormancy of the parasites. Our results suggest that long-distance migration may reduce the prevalence of parasite infections at the population level. Our results are inconsistent with the migratory exposure hypothesis, which posits that long-distance migration increases exposure of hosts to diverse parasites, and with the migratory susceptibility hypothesis, which posits that trade-offs between immune function and migration increase host susceptibility to parasites. However, our results are consistent with the migratory culling hypothesis, which posits that heavily infected animals are less likely to survive long-distance migration, and with the migratory escape hypothesis, which posits that long-distance migration allows host populations to seasonally escape areas of high infection risk.

Testing hypotheses of Pleistocene population history using coalescent simulations: phylogeography of the pygmy nuthatch (Sitta pygmaea).

In this paper, we use mitochondrial NADH dehydrogenase subunit 2 sequences to test Pleistocene refugial hypotheses for the pygmy nuthatch (Sitta pygmaea). Pygmy nuthatches are a common resident of long-needle pine forests in western North America and demonstrate a particular affinity with ponderosa pine (Pinus ponderosa). Palaeoecological and genetic data indicate that ponderosa pine was isolated in two Pleistocene refugia corresponding to areas in the southern Sierra Nevada in the west and southern Arizona and New Mexico in the east. We use coalescent simulations to test the hypothesis that pygmy nuthatches tracked the Pleistocene history of their preferred habitat and persisted in two refugia during the periods of glacial maxima. Coalescent simulation of population history does not support the hypothesis of two Pleistocene refugia for the pygmy nuthatch. Instead, our data are consistent with a single refuge model. Nucleotide diversity is greatest in the western populations of southern and coastal California. We suggest that the pygmy nuthatch expanded from a far western glacial refuge into its current distribution since the most recent glacial maximum.

Is gene flow promoting the reversal of pleistocene divergence in the Mountain Chickadee (Poecile gambeli)?

The Pleistocene glacial cycles left a genetic legacy on taxa throughout the world; however, the persistence of genetic lineages that diverged during these cycles is dependent upon levels of gene flow and introgression. The consequences of secondary contact among taxa may reveal new insights into the history of the Pleistocene's genetic legacy. Here, we use phylogeographic methods, using 20 nuclear loci from regional populations, to infer the consequences of secondary contact following divergence in the Mountain Chickadee (Poecile gambeli). Analysis of nuclear data identified two geographically-structured genetic groups, largely concordant with results from a previous mitochondrial DNA (mtDNA) study. Additionally, the estimated multilocus divergence times indicate a Pleistocene divergence, and are highly concordant with mtDNA. The previous mtDNA study showed a paucity of sympatry between clades, while nuclear patterns of gene flow show highly varied patterns between populations. The observed pattern of gene flow, from coalescent-based analyses, indicates southern populations in both clades exhibit little gene flow within or between clades, while northern populations are experiencing higher gene flow within and between clades. If this pattern were to persist, it is possible the historical legacy of Pleistocene divergence may be preserved in the southern populations only, and the northern populations would become a genetically diverse hybrid species.