Rapid postglacial diversification and long-term stasis within the songbird genus Junco: phylogeographic and phylogenomic evidence

Genetic evidence for widespread population size expansion in North American boreal birds prior to the Last Glacial Maximum

Pleistocene climate cycles are well documented to have shaped contemporary species distributions and genetic diversity. Northward range expansions in response to deglaciation following the Last Glacial Maximum (LGM; approximately 21 000 years ago) are surmised to have led to population size expansions in terrestrial taxa and changes in seasonal migratory behaviour. Recent findings, however, suggest that some northern temperate populations may have been more stable than expected through the LGM. We modelled the demographic history of 19 co-distributed boreal-breeding North American bird species from full mitochondrial gene sets and species-specific molecular rates. We used these demographic reconstructions to test how species with different migratory strategies were affected by glacial cycles. Our results suggest that effective population sizes increased in response to Pleistocene deglaciation earlier than the LGM, whereas genetic diversity was maintained throughout the LGM despite shifts in geographical range. We conclude that glacial cycles prior to the LGM have most strongly shaped contemporary genetic diversity in these species. We did not find a relationship between historic population dynamics and migratory strategy, contributing to growing evidence that major switches in migratory strategy during the LGM are unnecessary to explain contemporary migratory patterns.

Punctuated evolution in the learned songs of African sunbirds

Learned traits are thought to be subject to different evolutionary dynamics than other phenotypes, but their evolutionary tempo and mode has received little attention. Learned bird song has been thought to be subject to rapid and constant evolution. However, we know little about the evolutionary modes of learned song divergence over long timescales. Here, we provide evidence that aspects of the territorial songs of Eastern Afromontane sky island sunbirds Cinnyris evolve in a punctuated fashion, with periods of stasis of the order of hundreds of thousands of years or more, broken up by evolutionary pulses. Stasis in learned songs is inconsistent with learned traits being subject to constant or frequent change, as would be expected if selection does not constrain song phenotypes over evolutionary timescales. Learned song may instead follow a process resembling peak shifts on adaptive landscapes. While much research has focused on the potential for rapid evolution in bird song, our results suggest that selection can tightly constrain the evolution of learned songs over long timescales. More broadly, these results demonstrate that some aspects of highly variable, plastic traits can exhibit punctuated evolution, with stasis over long time periods.

Sequential colonization of oceanic archipelagos led to a species-level radiation in the common chaffinch complex (Aves: Fringilla coelebs)

Oceanic archipelagos are excellent systems for studying speciation, yet inference of evolutionary process requires that the colonization history of island organisms be known with accuracy. Here, we used phylogenomics and patterns of genetic diversity to infer the sequence and timing of colonization of Macaronesia by mainland common chaffinches (Fringilla coelebs), and assessed whether colonization of the different archipelagos has resulted in a species-level radiation. To reconstruct the evolutionary history of the complex we generated a molecular phylogeny based on genome-wide SNP loci obtained from genotyping-by-sequencing, we ran ancestral range biogeographic analyses, and assessed fine-scale genetic structure between and within archipelagos using admixture analysis. To test for a species-level radiation, we applied a probabilistic tree-based species delimitation method (mPTP) and an integrative taxonomy approach including phenotypic differences. Results revealed a circuitous colonization pathway in Macaronesia, from the mainland to the Azores, followed by Madeira, and finally the Canary Islands. The Azores showed surprisingly high genetic diversity, similar to that found on the mainland, and the other archipelagos showed the expected sequential loss of genetic diversity. Species delimitation methods supported the existence of several species within the complex. We conclude that the common chaffinch underwent a rapid radiation across Macaronesia that was driven by the sequential colonization of the different archipelagos, resulting in phenotypically and genetically distinct, independent evolutionary lineages. We recommend a taxonomic revision of the complex that takes into account its genetic and phenotypic diversity.

Temperature heterogeneity correlates with intraspecific variation in physiological flexibility in a small endotherm

Phenotypic flexibility allows individuals to reversibly modify trait values and theory predicts an individual's relative degree of flexibility positively correlates with the environmental heterogeneity it experiences. We test this prediction by integrating surveys of population genetic and physiological variation with thermal acclimation experiments and indices of environmental heterogeneity in the Dark-eyed Junco (Junco hyemalis) and its congeners. We combine field measures of thermogenic capacity for 335 individuals, 22,006 single nucleotide polymorphisms genotyped in 181 individuals, and laboratory acclimations replicated on five populations. We show that Junco populations: (1) differ in their thermogenic responses to temperature variation in the field; (2) harbor allelic variation that also correlates with temperature heterogeneity; and (3) exhibit intra-specific variation in thermogenic flexibility in the laboratory that correlates with the heterogeneity of their native thermal environment. These results provide comprehensive support that phenotypic flexibility corresponds with environmental heterogeneity and highlight its importance for coping with environmental change.

Females as the Gatekeepers to Seasonal Breeding: What We Can Learn by Studying Reproductive Mechanisms in Both Sexes

Seasonal reproduction is a widespread adaptation in vertebrates, such that individuals time their reproductive efforts to match peak resource abundance. Individuals rely on environmental cues to regulate hormonal mechanisms governing timing of breeding. Historically, studies on physiological mechanisms of seasonal reproduction, specifically in birds, have disproportionately focused on males compared to females. For this review, I conducted a literature search of the last decade of avian research and found a persistent sex bias in the field of physiological mechanisms of seasonal reproduction. Using work conducted with the dark-eyed junco (Junco hyemalis) as a case study, I present a possible solution to combat the sex bias: natural comparisons of populations that differ in reproductive timing to investigate mechanisms of reproduction in both sexes. Populations of dark-eyed juncos that differ in migratory behavior (i.e., migrant and resident) exhibit overlapping ranges during winter and early spring; residents begin breeding in early spring prior to the departure of migrants. This system, and others like it, provides an opportunity to compare mechanisms of reproduction in populations that differ in reproductive timing despite experiencing the same environmental conditions in early spring. In juncos, migrant and resident females and males exhibit similar patterns of hypothalamic regulation of reproduction in early spring, but sex differences in gonadal sensitivity between the populations could be an important distinction that partially explains sex differences in reproductive development. Comparing mechanisms of reproduction in free-living populations and in captivity can reveal important mechanisms that determine the onset of reproductive development, as well as potential sex differences in these mechanisms. Understanding the mechanisms of reproductive phenology has important implications for understanding how species will survive and reproduce in a changing climate.

Molecular signatures of long-distance oceanic dispersal and the colonization of Pacific islands in Lycium carolinianum

PREMISE

Long-distance dispersal has been important in explaining the present distributions of many plant species. Despite being infrequent, such dispersal events have considerable evolutionary consequences, because bottlenecks during colonization can result in reduced genetic diversity. We examined the phylogeographic history of Lycium carolinianum, a widespread taxon that ranges from southeastern North America to several Pacific islands, with intraspecific diversity in sexual and mating systems.

METHODS

We used Bayesian, likelihood, and coalescent approaches with nuclear and plastid sequence data and genome-wide single nucleotide polymorphisms to reconstruct the dispersal history of this species. We also compared patterns of genetic variation in mainland and island populations using single nucleotide polymorphisms and allelic diversity at the S-RNase mating system gene.

RESULTS

Lycium carolinianum is monophyletic and dispersed once from the North American mainland, colonizing the Pacific islands ca. 40,100 years ago. This dispersal was accompanied by a loss of genetic diversity in SNPs and the S-RNase locus due to a colonization bottleneck and the loss of self-incompatibility. Additionally, we documented at least two independent transitions to gynodioecy: once following the colonization of the Hawaiian Islands and loss of self-incompatibility, and a second time associated with polyploidy in the Yucatán region of Mexico.

CONCLUSIONS

Long-distance dispersal via fleshy, bird dispersed fruits best explains the unusually widespread distribution of L. carolinianum. The collapse of diversity at the S-RNase locus in island populations suggests that self-fertilization may have facilitated the subsequent colonization of Pacific islands following a single dispersal from mainland North America.

Change in sexual signalling traits outruns morphological divergence across an ecological gradient in the post-glacial radiation of the songbird genus Junco

The relative roles of natural and sexual selection in promoting evolutionary lineage divergence remains controversial and difficult to assess in natural systems. Local adaptation through natural selection is known to play a central role in promoting evolutionary divergence, yet secondary sexual traits can vary widely among species in recent radiations, suggesting that sexual selection may also be important in the early stages of speciation. Here, we compare rates of divergence in ecologically relevant traits (morphology) and sexually selected signalling traits (coloration) relative to neutral structure in genome-wide molecular markers and examine patterns of variation in sexual dichromatism to explore the roles of natural and sexual selection in the diversification of the songbird genus Junco (Aves: Passerellidae). Juncos include divergent lineages in Central America and several dark-eyed junco (J. hyemalis) lineages that diversified recently as the group recolonized North America following the last glacial maximum (ca. 18,000 years ago). We found an accelerated rate of divergence in sexually selected characters relative to ecologically relevant traits. Moreover, sexual dichromatism measurements suggested a positive relationship between the degree of colour divergence and the strength of sexual selection when controlling for neutral genetic distance. We also found a positive correlation between dichromatism and latitude, which coincides with the geographic axis of decreasing lineage age in juncos but also with a steep ecological gradient. Finally, we found significant associations between genome-wide variants linked to functional genes and proxies of both sexual and natural selection. These results suggest that the joint effects of sexual and ecological selection have played a prominent role in the junco radiation.

What Do Ecology, Evolution, and Behavior Have in Common? The Organism in the Middle

Biologists who publish in The American Naturalist are drawn to its unifying mission of covering research in the fields of ecology, evolution, behavior, and integrative biology. Presented here is one scientist's attempt to straddle these fields by focusing on a single organism. It is also an account of how time spent in the field stimulates a naturalist to wonder "why did that animal just do that?" and how research is guided by chance and intention interacting with the scientific literature and the people one meets along the way. With respect to the science, the examples come from bird migration, hormones and their connection to phenotypic integration, sexual and natural selection, and urban ecology. They also come from research on the impact of environmental change on the timing of reproduction and the potential for allochrony in migratory species to influence population divergence.

Adaptive genetic diversity and evidence of population genetic structure in the endangered Sierra Madre Sparrow (Xenospiza baileyi)

The magnitude and distribution of genetic diversity through space and time can provide useful information relating to evolutionary potential and conservation status in threatened species. In assessing genetic diversity in species that are of conservation concern, several studies have focused on the use of Toll-like receptors (TLRs). TLRs are innate immune genes related to pathogen resistance, and polymorphisms may reflect not only levels of functional diversity, but may also be used to assess genetic diversity within and among populations. Here, we combined four potentially adaptive markers (TLRs) with one mitochondrial (COI) marker to evaluate genetic variation in the endangered Sierra Madre Sparrow (Xenospiza baileyi). This species offers an ideal model to investigate population and evolutionary genetic processes that may be occurring in a habitat restricted endangered species with disjunct populations (Mexico City and Durango), the census sizes of which differ by an order of magnitude. TLRs diversity in the Sierra Madre Sparrow was relatively high, which was not expected given its two small, geographically isolated populations. Genetic diversity was different (but not significantly so) between the two populations, with less diversity seen in the smaller Durango population. Population genetic structure between populations was due to isolation and different selective forces acting on different TLRs; population structure was also evident in COI. Reduction of genetic diversity in COI was observed over 20 years in the Durango population, a result likely caused by habitat loss, a factor which may be the main cause of diversity decline generally. Our results provide information related to the ways in which adaptive variation can be altered by demographic changes due to human-mediated habitat alterations. Furthermore, our findings may help to guide conservation schemes for both populations and their restricted habitat.

A simple dynamic model explains the diversity of island birds worldwide

Colonization, speciation and extinction are dynamic processes that influence global patterns of species richness^1-6. Island biogeography theory predicts that the contribution of these processes to the accumulation of species diversity depends on the area and isolation of the island^7,8. Notably, there has been no robust global test of this prediction for islands where speciation cannot be ignored⁹, because neither the appropriate data nor the analytical tools have been available. Here we address both deficiencies to reveal, for island birds, the empirical shape of the general relationships that determine how colonization, extinction and speciation rates co-vary with the area and isolation of islands. We compiled a global molecular phylogenetic dataset of birds on islands, based on the terrestrial avifaunas of 41 oceanic archipelagos worldwide (including 596 avian taxa), and applied a new analysis method to estimate the sensitivity of island-specific rates of colonization, speciation and extinction to island features (area and isolation). Our model predicts-with high explanatory power-several global relationships. We found a decline in colonization with isolation, a decline in extinction with area and an increase in speciation with area and isolation. Combining the theoretical foundations of island biogeography^7,8 with the temporal information contained in molecular phylogenies¹⁰ proves a powerful approach to reveal the fundamental relationships that govern variation in biodiversity across the planet.

Evaluating testosterone as a phenotypic integrator: From tissues to individuals to species

Hormones have the potential to bring about rapid phenotypic change; however, they are highly conserved over millions of years of evolution. Here, we examine the evolution of hormone-mediated phenotypes, and the extent to which regulation is achieved via independence or integration of the many components of endocrine systems. We focus on the sex steroid testosterone (T), its cognate receptor (androgen receptor) and related endocrine components. We pose predictions about the mechanisms underlying phenotypic integration, including coordinated sensitivity to T within and among tissues and along the HPG axis. We then assess these predictions with case studies from wild birds, asking whether gene expression related to androgenic signaling naturally co-varies among individuals in ways that would promote phenotypic integration. Finally, we review how mechanisms of integration and independence vary over developmental or evolutionary time, and we find limited support for integration.

Genomic and mitochondrial evidence of ancient isolations and extreme introgression in the four-lined snake

Comparing mitochondrial and genomic phylogenies is an essential tool for investigating speciation processes, because each genome carries different inheritance properties and evolutionary characteristics. Furthermore, mitonuclear discordance may arise from ecological adaptation, historic isolation, population size changes, and sex-biased dispersal. Closely related taxa are expected to experience gene flow; however, this may not be true for insular populations or populations isolated in refugia. The four-lined snake Elaphe quatuorlineata has a fragmented distribution, separating populations of the Italian and Balkan Peninsulas, whereas several insular Aegean populations of significantly smaller body size (Cyclades island group and Skyros Island, Greece) are currently considered distinct subspecies. We constructed the species-tree phylogeny of this species utilizing genome-wide single nucleotide polymorphisms and a gene-tree based on complete cytochrome b sequences, aiming to detect convergence and discrepancies between biparentally and maternally inherited genomes. Population structuring, phylogenetic patterns and migration events among geographically defined lineages supported our hypothesis of isolation in multiple sub-refugia. Where biogeographical barriers did not restrict subsequent dispersal, extensive genetic exchange occurred between mainland Balkan populations. This process has led to the mitochondrial sweep of an ancestral mitolineage that survived only in peripheral (East Greece) and insular populations (North Cyclades and Skyros). The Central Cyclades represent an ancient lineage for both molecular markers that emerged almost 3.3 Mya. Considering their distinct morphology, insular E. quatuorlineata populations should be the future focus of an extensive sampling, especially since the mitonuclear discordance observed in this species could be related to ecological adaptations, such as the island-dwarfism phenomenon.

Genome-wide signals of drift and local adaptation during rapid lineage divergence in a songbird

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (a) drift in geographic isolation, (b) differentiation along continuous selective gradients and (c) isolation-by-adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci.

Male courtship preference during seasonal sympatry may maintain population divergence

Animal migration can lead to a population distribution known as seasonal sympatry, in which closely-related migrant and resident populations of the same species co-occur in sympatry during part of the year, but are otherwise allopatric. During seasonal sympatry in early spring, residents may initiate reproduction before migrants depart, presenting an opportunity for gene flow. Differences in reproductive timing between migrant and resident populations may favor residents that exhibit preferences for potential mates of similar migratory behavior and reproductive timing, thus maintaining population divergence. We studied dark-eyed juncos (Junco hyemalis), a songbird that exhibits seasonal sympatry. We conducted simulated courtship interactions in which we presented free-living resident males with either a caged migrant or resident female and quantified courtship behavior prior to the departure of the migrants. We found that resident males preferred to court resident females: they sang more short-range songs and exhibited more visual displays associated with courtship when presented with resident females. We conclude that males distinguish between migrant and resident females during seasonal sympatry when the risk of interacting with non-reproductive, migrant females is high. Male mate choice in seasonal sympatry is likely adaptive for male reproductive success. As a secondary effect, male mating preference could act to maintain or promote divergence between populations that differ in migratory strategy.

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.

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection

Background

Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome.

Results

We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima’s D. However, we found few genomic regions that showed both elevated polymorphism and reduced F_ST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested.

Conclusions

Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4806-7) contains supplementary material, which is available to authorized users.