Search

Large structural variants in the genome, such as inversions, may play an important role in producing population structure and local adaptation to the environment through suppression of recombination. However, relatively few studies have linked inversions to phenotypic traits that are sexually selected and may play a role in reproductive isolation. Here, we found that geographic differences in the sexually selected plumage of a warbler, the common yellowthroat (Geothlypis trichas), are largely due to differences in the Z (sex) chromosome (males are ZZ), which contains at least one putative inversion spanning 40% (31/77 Mb) of its length. The inversions on the Z chromosome vary dramatically east and west of the Appalachian Mountains, which provides evidence of cryptic population structure within the range of the most widespread eastern subspecies (G. t. trichas). In an eastern (New York) and western (Wisconsin) population of this subspecies, female prefer different male ornaments; larger black facial masks are preferred in Wisconsin and larger yellow breasts are preferred in New York. The putative inversion also contains genes related to vision, which could influence mating preferences. Thus, structural variants on the Z chromosome are associated with geographic differences in male ornaments and female choice, which may provide a mechanism for maintaining different patterns of sexual selection in spite of gene flow between populations of the same subspecies., (© 2024 John Wiley & Sons Ltd.)

In humans, gut microbiome (GM) differences are often correlated with, and sometimes causally implicated in, ageing. However, it is unclear how these findings translate in wild animal populations. Studies that investigate how GM dynamics change within individuals, and with declines in physiological condition, are needed to fully understand links between chronological age, senescence and the GM, but have rarely been done. Here, we use longitudinal data collected from a closed population of Seychelles warblers (Acrocephalus sechellensis) to investigate how bacterial GM alpha diversity, composition and stability are associated with host senescence. We hypothesised that GM diversity and composition will differ, and become more variable, in older adults, particularly in the terminal year prior to death, as the GM becomes increasingly dysregulated due to senescence. However, GM alpha diversity and composition remained largely invariable with respect to adult age and did not differ in an individual's terminal year. Furthermore, there was no evidence that the GM became more heterogenous in senescent age groups (individuals older than 6 years), or in the terminal year. Instead, environmental variables such as season, territory quality and time of day, were the strongest predictors of GM variation in adult Seychelles warblers. These results contrast with studies on humans, captive animal populations and some (but not all) studies on non-human primates, suggesting that GM deterioration may not be a universal hallmark of senescence in wild animal species. Further work is needed to disentangle the factors driving variation in GM-senescence relationships across different host taxa., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

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., (© 2023 John Wiley & Sons Ltd.)

Exposure to rising sublethal temperatures can affect development and somatic condition, and thereby Darwinian fitness. In the context of climate warming, these changes could have implications for population viability, but they can be subtle and consequently difficult to quantify. Using telomere length (TL) as a known biomarker of somatic condition in early life, we investigated the impact of pre-hatching and nestling climate on six cohorts of wild nestling superb fairy wrens (Malurus cyaneus) in temperate south-eastern Australia. Models incorporating only climate information from the nestling phase were best supported compared to those including the (pre-)laying to incubation phase (previously shown to affect mass) or both phases combined. This implies that nestling TL is most sensitive to ambient climate in the nestling phase. The top model showed a negative relationship between early-life TL and nestling mean daily minimum temperature when rainfall was low which gradually became positive with increasing rainfall. In addition, there was a positive relationship between TL and the frequency of hot days (daily maximum temperature ≥35°C), although these temperatures were rare and short-term. Including other pre-hatching and nestling period, climate variables (e.g., mean daily maximum temperature and mean diurnal temperature variability) did not improve the prediction of nestling TL. Overall, our results suggest that cooler nights when conditions are dry and short-term temperature spikes above 35°C during development are conducive for somatic maintenance. While these findings indicate a potential pathway for climate warming to impact wildlife fitness, they emphasize the need to elucidate the mechanisms underlying these complex associations., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Understanding the geographic linkages among populations across the annual cycle is an essential component for understanding the ecology and evolution of migratory species and for facilitating their effective conservation. While genetic markers have been widely applied to describe migratory connections, the rapid development of new sequencing methods, such as low-coverage whole genome sequencing (lcWGS), provides new opportunities for improved estimates of migratory connectivity. Here, we use lcWGS to identify fine-scale population structure in a widespread songbird, the American Redstart (Setophaga ruticilla), and accurately assign individuals to genetically distinct breeding populations. Assignment of individuals from the nonbreeding range reveals population-specific patterns of varying migratory connectivity. By combining migratory connectivity results with demographic analysis of population abundance and trends, we consider full annual cycle conservation strategies for preserving numbers of individuals and genetic diversity. Notably, we highlight the importance of the Northern Temperate-Greater Antilles migratory population as containing the largest proportion of individuals in the species. Finally, we highlight valuable considerations for other population assignment studies aimed at using lcWGS. Our results have broad implications for improving our understanding of the ecology and evolution of migratory species through conservation genomics approaches., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Across the animal kingdom, newly independent juveniles form social associations that influence later fitness, mate choice and gene flow, but little is known about the ontogeny of social environments, particularly in wild populations. Here we test whether associations among young animals form randomly or are influenced by environmental or genetic conditions established by parents. Parents' decisions determine natal birth sites, which could affect who independent young initially encounter; secondly, mate choice determines genetic condition (e.g. inbreeding) of young and the parental care they receive, which can affect sociability. However, genetic and environmental factors are confounded unless related offspring experience different natal environments. Therefore, we used a long-term genetic pedigree, breeding records and social network data from three cohorts of a songbird with high extra-pair paternity (hihi, Notiomystis cincta) to disentangle (1) how nest location and relatedness contribute to association structure once juveniles disperse away from birth sites, and (2) if juvenile and/or parental inbreeding predicts individual sociability. We detected positive spatial autocorrelation: hihi that fledged closer by were more likely to associate even after dispersing, irrespective of genetic relatedness. Juvenile inbreeding did not predict sociability, but those raised by more inbred fathers formed more, stronger, associations, which did not depend on whether that male was the genetic parent or not. These results suggest that the natal environment created by parents, rather than focal genetic condition, establishes the foundation for social associations. Overall, we highlight how social inheritance may play an important role in population dynamics and evolutionary potential in wild animals., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Telomeres are well known for their associations with lifespan and ageing across diverse taxa. Early-life telomere length can be influenced by developmental conditions and has been shown positively affect lifetime reproductive success in a limited number of studies. Whether these effects are caused by a change in lifespan, reproductive rate or perhaps most importantly reproductive senescence is unclear. Using long-term data on female breeding success from a threatened songbird (the hihi, Notiomystis cincta), we show that the early-life telomere length of individuals predicts the presence and rate of future senescence of key reproductive traits: clutch size and hatching success. In contrast, senescence of fledging success is not associated with early-life telomere length, which may be due to the added influence of biparental care at this stage. Early-life telomere length does not predict lifespan or lifetime reproductive success in this species. Females may therefore change their reproductive allocation strategy depending on their early developmental conditions, which we hypothesise are reflected in their early-life telomere length. Our results offer new insights on the role that telomeres play in reproductive senescence and individual fitness and suggest telomere length can be used as a predictor for future life history in threatened species., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Suboptimal conditions during development can shorten telomeres, the protective DNA caps on the end of chromosomes. Shorter early-life telomere length (TL) can indicate reduced somatic maintenance, leading to lower survival and shorter lifespan. However, despite some clear evidence, not all studies show a relationship between early-life TL and survival or lifespan, which may be due to differences in biology or study design (e.g., survival period measured). In superb fairy-wrens (Malurus cyaneus), we assessed whether early-life TL predicts mortality across different life-history stages (fledgling, juvenile, adult). However, in contrast to a similar study on a congener, early-life TL did not predict mortality across any life stage in this species. We then performed a meta-analysis including 32 effect sizes from 23 studies (15 birds and three mammals) to quantify the effect of early-life TL on mortality whilst taking into consideration potential sources of biological and methodological variation. Overall, the effect of early-life TL on mortality was significant, corresponding to a 15% reduction in mortality risk with each standard deviation increase in TL. However, the effect became weaker when correcting for publication bias. Contrary to our predictions, there was no evidence that effects of early-life TL on mortality varied with species lifespan or the period over which survival was measured. However, negative effects of early-life TL on mortality risk were pervasive throughout life. These results imply that effects of early-life TL on mortality are more likely to be context-dependent than age-dependent, although substantial power and publication bias issues highlight the need for more research., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Organisms often acquire physiological and morphological modifications to conquer ecological challenges when colonizing new environments which lead to their adaptive evolution. However, deciphering the genomic mechanism of ecological adaptation is difficult because ecological environments are often too complex for straightforward interpretation. Thus, we examined the adaptation of a widespread songbird-the rufous-capped babbler (Cyanoderma ruficeps)-to a relatively simple system: distinct environments across elevational gradients on the mountainous island of Taiwan. We focused on the genomic sequences of 43 birds from five populations to show that the Taiwan group split from its sister group in mainland China around 1-2 million years ago (Ma) and colonized the montane habitats of Taiwan at least twice around 0.03-0.22 Ma. The montane and lowland Taiwan populations diverged with gene flow between them, suggesting strong selection associated with different elevations. We found that the montane babblers had smaller beaks than the lowland ones, consistent with Allen's rule, and identified candidate genes-COL9A1 and SOX11-underlying the beak size changes. We also found that altitudinally divergent mutations were mostly located in noncoding regions and tended to accumulate in chromosomal inversions and autosomes. The altitudinally divergent mutations might regulate genes related to haematopoietic, metabolic, immune, auditory and vision functions, as well as cerebrum morphology and plumage development. The results reveal the genomic bases of morphological and physiological adaptation in this species to the low temperature, hypoxia and high UV light environment at high elevation. These findings improve our understanding of how ecological adaptation drives population divergence from the perspective of genomic architecture., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

The seasonal migration of birds is a fascinating natural wonder. Avian migratory behaviour changes are common and are probably a polygenic process as avian migration is governed by multiple correlated components with a variable genetic basis. However, the genetic and phenotypic changes involving migration changes are poorly studied. Using one annotated near-chromosomal level de novo genome assembly, 50 resequenced genomes, hundreds of morphometric data and species distribution information, we investigated population structure and genomic and phenotypic differences associated with differences in migratory behaviour in a songbird species, Yellow-throated Bunting Emberiza elegans (Aves: Emberizidae). Population genomic analyses reveal extensive gene flow between the southern resident and the northern migratory populations of this species. The hand-wing index is significantly lower in the resident populations than in the migratory populations, indicating reduced flight efficiency of the resident populations. Here, we discuss the possibility that nonmigratory populations may have originated from migratory populations though migration loss. We further infer that the alterations of genes related to energy metabolism, nervous system and circadian rhythm may have played major roles in regulating migration change. Our study sheds light on phenotypic and polygenic changes involving migration change., (© 2022 John Wiley & Sons Ltd.)

Telomere length and shortening rate are increasingly being used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between-population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the nonbreeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north-south gradient (Finland, Estonia, England and Spain) predicting a negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from midlatitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Reproductive investment often comes at a cost to longevity, but the mechanisms that underlie these long-term effects are not well understood. In male vertebrates, elevated testosterone has been shown to increase reproductive success, but simultaneously to decrease survival. One factor that may contribute to or serve as a biomarker of these long-term effects of testosterone on longevity is telomeres, which are often positively related to lifespan and have been shown to shorten in response to reproduction. In this longitudinal study, we measured the effects of experimentally elevated testosterone on telomere shortening in free-living, male dark-eyed juncos (Junco hyemalis carolinensis), a system in which the experimental elevation of testosterone has previously been shown to increase reproductive success and reduce survival. We found a small, significant effect of testosterone treatment on telomeres, with testosterone-treated males exhibiting significantly greater telomere shortening with age than controls. These results are consistent with the hypothesis that increased telomere shortening may be a long-term cost of elevated testosterone exposure. As both testosterone and telomeres are conserved physiological mechanisms, our results suggest that their interaction may apply broadly to the long-term costs of reproduction in male vertebrates., (© 2021 John Wiley & Sons Ltd.)

Telomeres are protective, nucleoprotein structures at the end of chromosomes that have been associated with lifespan across taxa. However, the extent to which these associations can be attributed to absolute length vs. the rate of telomere shortening prior to sampling remains unresolved. In a longitudinal study, we examined the relationship between lifespan, telomere length and the rate of telomere shortening in wild, purple-crowned fairy-wrens (Malurus coronatus coronatus). To this end, we measured telomere length using quantitative polymerase chain reaction in the blood of 59 individuals sampled as nestlings and 4-14 months thereafter, and in 141 known-age individuals sampled on average three times across adulthood. We applied within-subject centring analyses to simultaneously test for associations between lifespan and average telomere length and telomere shortening. We reveal that the rate of telomere shortening and to a lesser extent telomere length in the first year of life independently predicted lifespan, with individuals with faster shortening rates and/or shorter telomeres living less long. In contrast, in adulthood neither telomere shortening nor telomere length predicted lifespan, despite a considerably larger data set. Our results suggest that telomere length measured very early in life (during development) and longitudinal assessments of telomere shortening during the first year of life constitute more useful biomarkers of total life expectancy than either telomere length measured after development, or telomere shortening later in adulthood., (© 2021 John Wiley & Sons Ltd.)

While an increasing number of studies are adopting molecular and chemical methods for dietary characterization, these studies often employ only one of these laboratory-based techniques; this approach may yield an incomplete, or even biased, understanding of diet due to each method's inherent limitations. To explore the utility of coupling molecular and chemical techniques for dietary characterizations, we applied DNA metabarcoding alongside stable isotope analysis to characterize the dietary niche of breeding Louisiana waterthrush (Parkesia motacilla), a migratory songbird hypothesized to preferentially provision its offspring with pollution-intolerant, aquatic arthropod prey. While DNA metabarcoding was unable to determine if waterthrush provision aquatic and terrestrial prey in different abundances, we found that specific aquatic taxa were more likely to be detected in successive seasons than their terrestrial counterparts, thus supporting the aquatic specialization hypothesis. Our isotopic analysis added greater context to this hypothesis by concluding that breeding waterthrush provisioned Ephemeroptera and Plecoptera, two pollution-intolerant, aquatic orders, in higher quantities than other prey groups, and expanded their functional trophic niche when such prey were not abundantly provisioned. Finally, we found that the dietary characterizations from each approach were often uncorrelated, indicating that the results gleaned from a diet study can be particularly sensitive to the applied methodologies. Our findings contribute to a growing body of work indicating the importance of high-quality, aquatic habitats for both consumers and their pollution-intolerant prey, while also demonstrating how the application of multiple, laboratory-based techniques can provide insights not offered by either technique alone., (© 2022 John Wiley & Sons Ltd.)

In a rapidly warming world, exposure to high temperatures may impact fitness, but the gene regulatory mechanisms that link sublethal heat to sexually selected traits are not well understood, particularly in endothermic animals. Our experiment used zebra finches (Taeniopygia guttata), songbirds that experience extreme temperature fluctuations in their native Australia. We exposed captive males to an acute thermal challenge (43°C) compared with thermoneutral (35°C) and lower (27°C) temperatures. We found significantly more heat dissipation behaviours at 43°C, a temperature previously shown to reduce song production and fertility, and more heat retention behaviours at 27°C. Next, we characterized transcriptomic responses in tissues important for mating effort-the posterior telencephalon, for its role in song production, and the testis, for its role in fertility and hormone production. Differential expression of hundreds of genes in the testes, but few in the brain, suggests the brain is less responsive to extreme temperatures. Nevertheless, gene network analyses revealed that expression related to dopaminergic signalling in the brain covaried with heat dissipation behaviours, providing a mechanism by which temporary thermal challenges may alter motivational circuits for song production. In both brain and testis, we observed correlations between thermally sensitive gene networks and individual differences in thermoregulatory behaviour. Although we cannot directly relate these gene regulatory changes to mating success, our results suggest that individual variation in response to thermal challenges could impact sexually selected traits in a warming world., (© 2022 John Wiley & Sons Ltd.)

Colonization of a novel environment by a few individuals can lead to rapid evolutionary change, yet there is scarce evidence of the relative contributions of neutral and selective factors in promoting divergence during the early stages of colonization. Here we explore the role of neutral and selective forces in the divergence of a unique urban population of the dark-eyed junco (Junco hyemalis), which became established on the campus of the University of California at San Diego (UCSD) in the early 1980s. Previous studies based on microsatellite loci documented significant genetic differentiation of the urban population as well as divergence in phenotypic traits relative to nearby montane populations, yet the geographical origin of the colonization and the contributing factors remained uncertain. Our genome-wide single nucleotide polymorphism data set confirmed the marked genetic differentiation of the UCSD population, and we identified the coastal subspecies pinosus from central California as its sister group instead of the neighbouring mountain population. Demographic inference recovered a separation from pinosus as recent as 20-32 generations ago after a strong bottleneck, suggesting a role for drift in genetic differentiation. However, we also found significant associations between habitat variables and genome-wide variants linked to functional genes, some of which have been reported as potentially adaptive in birds inhabiting modified environments. These results suggest that the interplay between founder events and selection may result in rapid shifts in neutral and adaptive loci across the genome, and reveal the UCSD junco population as a case of contemporary evolutionary divergence in an anthropogenic environment., (© 2022 John Wiley & Sons Ltd.)

Understanding the genomic landscape of adaptation is central to understanding microevolution in wild populations. Genomic targets of selection and the underlying genomic mechanisms of adaptation can be elucidated by genome-wide scans for past selective sweeps or by scans for direct fitness associations. We sequenced and assembled 150 haplotypes of 75 blue tits (Cyanistes caeruleus) of a single Central European population by a linked-read technology. We used these genome data in combination with coalescent simulations (i) to estimate an historical effective population size of ~250,000, which recently declined to ~10,000, and (ii) to identify genome-wide distributed selective sweeps of beneficial variants probably originating from standing genetic variation (soft sweeps). The genes linked to these soft sweeps, but also those linked to hard sweeps based on new beneficial mutants, showed a significant enrichment for functions associated with gene expression and transcription regulation. This emphasizes the importance of regulatory evolution in the population's adaptive history. Soft sweeps were further enriched for genes related to axon and synapse development, indicating the significance of neuronal connectivity changes in the brain potentially linked to behavioural adaptations. A previous scan of heterozygosity-fitness correlations revealed a consistent negative effect on arrival date at the breeding site for a single microsatellite in the MDGA2 gene. Here, we used the haplotype structure around this microsatellite to explain the effect as a local and direct outbreeding effect of a gene involved in synapse development., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Genetic structure and phenotypic variation among populations are affected by both geographic distance and environmental variation across species' distributions. Understanding the relative contributions of isolation by distance (IBD) and isolation by environment (IBE) is important for elucidating population dynamics across habitats and ecological gradients. In this study, we compared phenotypic and genetic variation among Horned Lark (Eremophila alpestris) populations from 10 sites encompassing an elevational gradient from low-elevation desert scrub in Death Valley (285 a.s.l.) to high-elevation meadows in the White Mountains of the Sierra Nevada of California (greater than 3000 m a.s.l.). Using a ddRAD data set of 28,474 SNPs aligned to a high-quality reference genome, we compared genetic structure with elevational, environmental, and spatial distance to quantify how different aspects of the landscape drive genomic and phenotypic differentiation in Horned Larks. We found larger-bodied birds were associated with sites that had less seasonality and higher annual precipitation, and longer spurs occurred in soils with more clay and silt content, less sand, and finer fragments. Larks have large neo-sex chromosomes, and we found that associations with elevation and environmental variation were much stronger among neo-sex chromosomes compared to autosomes. Furthermore, we found that putative chromosomal translocations, fusions, and inversions were associated with elevation and may underlie local adaptation across an elevational gradient in Horned Larks. Our results suggest that genetic variation in Horned Larks is affected more by IBD than IBE, but specific phenotypes and genomic regions-particually on neo-sex chromosomes-bear stronger associations with the environment., (© 2022 John Wiley & Sons Ltd.)

The genetic basis of bird migration has been the focus of several studies. Two willow warbler subspecies (Phylloscopus trochilus trochilus and Phylloscopus trochilus acredula) follow different migratory routes to wintering grounds in Africa. Their breeding populations overlap in contact areas or "migratory divides" located in central Scandinavia and in eastern Poland. Earlier analyses demonstrated that the genetic differences between these two migratory phenotypes are few and cluster on chromosomes 1 and 5. In addition, an amplified fragment length polymorphism-derived biallelic marker (known as WW2) presents steep clines across both migratory divides but failed to be mapped in the genome. Here, we characterize the WW2 marker and describe its two variants (WW2 ancestral and WW2 derived) as portions of long terminal repeat retrotransposons originating from an ancient infection by an endogenous retrovirus. We used quantitative polymerase chain reaction techniques to quantify copy numbers of the WW2 derived variant in the two subspecies and their hybrids. This, together with genome analyses revealed that WW2 derived variants are much more abundant in P. t. acredula and appear embedded in a large repeat-rich region (>12 Mbp), not associated with the divergent regions of chromosomes 1 or 5. However, it might interact with genetic elements controlling migration direction. Testing this hypothesis further will require knowing the exact location of this region, such as by obtaining more complete genome assemblies preferably in combination with techniques like fluorescence in situ hybridization applied to a willow warbler karyotype, and finally to investigate the copy number of this marker in hybrids with known migratory tracks., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Researchers seeking to generate genomic data for non-model organisms are faced with a number of trade-offs when deciding which method to use. The selection of reduced representation approaches versus whole genome resequencing will ultimately affect the marker density, sequencing depth, and the number of individuals that can multiplexed. These factors can affect researchers' ability to accurately characterize certain genomic features, such as landscapes of divergence-how F ST varies across the genomes. To provide insight into the effect of sequencing method on the estimation of divergence landscapes, we applied an identical bioinformatic pipeline to three generations of sequencing data (GBS, ddRAD, and WGS) produced for the same system, the yellow-rumped warbler species complex. We compare divergence landscapes generated using each method for the myrtle warbler (Setophaga coronata coronata) and the Audubon's warbler (S. c. auduboni), and for Audubon's warblers with deeply divergent mtDNA resulting from mitochondrial introgression. We found that most high-F ST peaks were not detected in the ddRAD data set, and that while both GBS and WGS were able to identify the presence of large peaks, WGS was superior at a finer scale. Comparing Audubon's warblers with divergent mitochondrial haplotypes, only WGS allowed us to identify small (10-20 kb) regions of elevated differentiation, one of which contained the nuclear-encoded mitochondrial gene NDUFAF3. We calculated the cost per base pair for each method and found it was comparable between GBS and WGS, but significantly higher for ddRAD. These comparisons highlight the advantages of WGS over reduced representation methods when characterizing landscapes of divergence., (© 2021 John Wiley & Sons Ltd.)

Geographic contact between sister lineages often occurs near the final stages of speciation, but its role in speciation's completion remains debated. Reproductive isolation may be essentially complete prior to secondary contact. Alternatively, costly interactions between partially reproductively isolated species - such as maladaptive hybridization or competition for resources - may select for divergence, increasing reproductive isolation and driving speciation toward completion. Here, we use coalescent demographic modelling and whole-genome data sets to show that a period of contact and elevated hybridization between sympatric eastern North American populations of two cryptic bird species preceded a major increase in reproductive isolation between these populations within the last 10,000 years. In contrast, substantial introgression continues to the present in a western contact zone where geographic overlap is much narrower and probably of more recent origin. In the sympatric eastern region where reproductive isolation has increased, it is not accompanied by character displacement in key morphometric traits, plumage coloration, or ecological traits. While the precise trait and underlying mechanism driving increased reproductive isolation remains unknown, we discuss several possibilities and outline avenues for future research. Overall, our results highlight how demographic models can reveal the geographic context in which reproductive isolation was completed, and demonstrate how contact can accelerate the final stages of speciation., (© 2021 John Wiley & Sons Ltd.)

Contact zones between recently diverged taxa provide opportunities to examine the causes of reproductive isolation and the processes that determine whether two species can coexist over a broad region. The Pacific wren (Troglodytes pacificus) and winter wren (Troglodytes hiemalis) are two morphologically similar songbirds that started diverging about 4 million years ago, older than most sister species pairs of temperate songbirds. The ranges of these species come into narrow contact in western Canada, where the two species remain distinct. To assess evidence for differentiation, hybridization and introgression in this system, we examined variation in over 250,000 single nucleotide polymorphism markers distributed across the genome. The two species formed highly divergent genetic clusters, consistent with long-term differentiation. In a set of 75 individuals, two first-generation hybrids (i.e., F 1 's) were detected, indicating only moderate levels of assortative mating between these taxa. We found no recent backcrosses or other evidence of recent breeding success of F 1 's, indicating very low or zero fitness of F 1  hybrids. Examination of genomic variation shows evidence for only a single backcrossing event many generations ago. The moderate rate of hybridization combined with very low F 1  hybrid fitness is expected to result in a population sink in the contact zone, largely explaining the narrow overlap of the two species. If such dynamics are common in nature, they could explain the narrow range overlap often observed between pairs of closely related species., (© 2021 John Wiley & Sons Ltd.)

Species with a circannual life cycle need to match the timing of their life history events to the environment to maximize fitness. However, our understanding of how circannual traits such as timing of reproduction are regulated on a molecular level remains limited. Recent studies have implicated that epigenetic mechanisms can be an important part in the processes that regulate circannual traits. Here, we explore the role of DNA methylation in mediating reproductive timing in a seasonally breeding bird species, the great tit (Parus major), using genome-wide DNA methylation data from individual females that were blood sampled repeatedly throughout the breeding season. We demonstrate rapid and directional changes in DNA methylation within the promoter region of several genes, including a key transcription factor (NR5A1) known from earlier studies to be involved in the initiation of timing of reproduction. Interestingly, the observed changes in DNA methylation at NR5A1 identified here are in line with earlier gene expression studies of reproduction in chicken, indicating that the observed shifts in DNA methylation at this gene can have a regulatory role. Our findings provide an important step towards elucidating the genomic mechanism that mediates seasonal timing of a key life history traits and provide support for the idea that epigenetic mechanisms may play an important role in circannual traits., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Due to its central importance to many aspects of evolutionary biology and population genetics, the long-term effective population size (N e ) has been estimated for numerous species and populations. However, estimating contemporary N e is difficult and in practice this parameter is often unknown. In principle, contemporary N e can be estimated using either analyses of temporal changes in allele frequencies, or the extent of linkage disequilibrium (LD) between unlinked markers. We applied these approaches to estimate contemporary N e of a relatively recently founded island population of collared flycatchers (Ficedula albicollis). We sequenced the genomes of 85 birds sampled in 1993 and 2015, and applied several temporal methods to estimate N e at a few thousand (4000-7000). The approach based on LD provided higher estimates of N e (20,000-32,000) and was associated with high variance, often resulting in infinite N e . We conclude that whole-genome sequencing data offers new possibilities to estimate high (>1000) contemporary N e , but also note that such estimates remain challenging, in particular for LD-based methods for contemporary N e estimation., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Morphological traits have served generations of biologists as a taxonomic indicator, and have been the main basis for defining and classifying species diversity for centuries. A quantitative integration of behavioural characters, such as vocalizations, in studies on biotic differentiation has arisen more recently, and the relative importance of these different traits in the diversification process remains poorly understood. To provide a framework within which to interpret the evolutionary interplay between morphological and behavioural traits, we generated a draft genome of a cryptic Southeast Asian songbird, the limestone wren-babbler Napothera crispifrons. We resequenced whole genomes of multiple individuals of all three traditional subspecies and of a distinct leucistic population. We demonstrate strong genomic and mitochondrial divergence among all three taxa, pointing to the existence of three species-level lineages. Despite its great phenotypic distinctness, the leucistic population was characterized by shallow genomic differentiation from its neighbour, with only a few localized regions emerging as highly diverged. Quantitative bioacoustic analysis across multiple traits revealed deep differences especially between the two taxa characterized by limited plumage differentiation. Our study demonstrates that differentiation in these furtive songbirds has resulted in a complex mosaic of colour-based and bioacoustic differences among populations. Extreme colour differences can be anchored in few genomic loci and may therefore arise and subside rapidly., (© 2020 John Wiley & Sons Ltd.)

Elucidating forces capable of driving species diversification in the face of gene flow remains a key goal in evolutionary biology. Song sparrows, Melospiza melodia, occur as 25 subspecies in diverse habitats across North America, are among the continent's most widespread vertebrate species, and are exemplary of many highly variable species for which the conservation of locally adapted populations may be critical to their range-wide persistence. We focus here on six morphologically distinct subspecies resident in the San Francisco Bay region, including three salt-marsh endemics and three residents in upland and riparian habitats adjacent to the Bay. We used reduced-representation sequencing to generate 2,773 SNPs to explore genetic differentiation, spatial population structure, and demographic history. Clustering separated individuals from each of the six subspecies, indicating subtle differentiation at microgeographic scales. Evidence of limited gene flow and low nucleotide diversity across all six subspecies further supports a hypothesis of isolation among locally adapted populations. We suggest that natural selection for genotypes adapted to salt marsh environments and changes in demography over the past century have acted in concert to drive the patterns of diversification reported here. Our results offer evidence of microgeographic specialization in a highly polytypic bird species long discussed as a model of sympatric speciation and rapid adaptation, and they support the hypothesis that conserving locally adapted populations may be critical to the range-wide persistence of similarly highly variable species., (© 2020 John Wiley & Sons Ltd.)

Multi-cellular organisms are under constant attack from parasites, making immune defence a critical aspect of fitness. In vertebrate animals, genes of the major histocompatibility complex (MHC) determine the breadth of pathogens to which individuals can respond. Having many MHC alleles can confer better protection against infectious disease, and balancing selection at MHC is widespread. Indeed, MHC loci are famously variable, with some populations harbouring thousands of alleles (Biedrzycka et al., 2018; Robinson, Soormally, Hayhurst, & Marsh, 2016). MHC has also long fascinated behavioural ecologists because mate choice-for example, preferring MHC-dissimilar partners-may amplify the effects of natural selection (Penn & Potts, 1999). But despite keen interest in the evolutionary ecology of MHC, extensive duplication (Minias, Pikus, Whittingham, & Dunn, 2019) has made these genes challenging to study. In a From the Cover article in this issue of Molecular Ecology, Stervander, Dierickx, Thorley, Brooke, and Westerdahl (2020) characterizes class I MHC in a Critically Endangered songbird, relating genotype to mate choice and survivorship. By inferring copy number and patterns of allelic co-segregation, the authors pave the way to elucidating the genomic architecture of MHC in this bottlenecked population. These insights help reconcile apparently counterintuitive findings: no effect of MHC genotype on mate choice or survival, and high MHC diversity within individuals despite low diversity at the population level. The latter finding is cause for optimism regarding conservation prospects. Moreover, these results suggest that ancient duplication events can have longstanding effects on the adaptive landscapes of natural and sexual selection., (© 2020 John Wiley & Sons Ltd.)

Extra-pair paternity (EPP) is often linked to male age in socially monogamous vertebrates; that is, older males are more likely to gain EPP and less likely to be cuckolded. However, whether this occurs because males improve at gaining paternity as they grow older, or because "higher quality" males that live longer are preferred by females, has rarely been tested, despite being central to our understanding of the evolutionary drivers of female infidelity. Moreover, how extra-pair reproduction changes with age within females has received even less attention. Using 18 years of longitudinal data from an individually marked population of Seychelles warblers (Acrocephalus sechellensis), we found considerable within-individual changes in extra-pair reproduction in both sexes: an early-life increase and a late-life decline. Furthermore, males were cuckolded less as they aged. Our results indicate that in this species age-related patterns of extra-pair reproduction are determined by within-individual changes with age, rather than differences among individuals in longevity. These results challenge the hypothesis-based on longevity reflecting intrinsic quality-that the association between male age and EPP is due to females seeking high-quality paternal genes for offspring. Importantly, EPP accounted for up to half of male reproductive success, emphasizing the male fitness benefits of this reproductive strategy. Finally, the occurrence of post-peak declines in extra-pair reproduction provides explicit evidence of senescence in infidelity in both males and females., (© 2020 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

The evolutionary consequences of individual genetic diversity are frequently studied by assessing heterozygosity-fitness correlations (HFCs). The prevalence of positive and negative HFCs and the predominance of general versus local effects in wild populations are far from understood, partly because comprehensive studies testing for both inbreeding and outbreeding depression are lacking. We studied a genetically diverse population of blue tits in southern Germany using a genome-wide set of 87 microsatellites to investigate the relationship between proxies of reproductive success and measures of multilocus and single-locus individual heterozygosity (MLH and SLH). We used complimentary measures of MLH and partitioned markers into functional categories according to their position in the blue tit genome. HFCs based on MLH were consistently negative for functional loci, whereas correlations were rather inconsistent for loci found in nonfunctional areas of the genome. Clutch size was the only reproductive variable showing a general effect. We found evidence for local effects for three measures of reproductive success: arrival date at the breeding site, the probability of breeding at the study site and male reproductive success. For these, we observed consistent, and relatively strong, negative effects at one functional locus. Remarkably, this marker had a similar effect in another blue tit population from Austria (~400 km to the east). We suggest that a genetic local effect on timing of arrival might be responsible for most negative HFCs detected, with carry-over effects on other reproductive traits. This effect could reflect individual differences in the distance between overwintering areas and breeding sites., (© 2020 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Several factors have been proposed as drivers of species diversification in the Neotropics, including environmental heterogeneity, the development of drainage systems and historical changes in forest distribution due to climatic oscillations. Here, we investigate which drivers contributed to the evolutionary history and current patterns of diversity of a polymorphic songbird (Arremon taciturnus) that is widely distributed in Amazonian and Atlantic forests as well as in Cerrado gallery and seasonally-dry forests. We use genomic, phenotypic and habitat heterogeneity data coupled with climatic niche modelling. Results suggest the evolutionary history of the species is mainly related to paleoclimatic changes, although changes in the strength of the Amazon river as a barrier to dispersal, current habitat heterogeneity and geographic distance were also relevant. We propose an ancestral distribution in the Guyana Shield, and recent colonization of areas south of the Amazon river at ~380 to 166 kya, and expansion of the distribution to southern Amazonia, Cerrado and the Atlantic Forest. Since then, populations south of the Amazon River have been subjected to cycles of isolation and possibly secondary contact due to climatic changes that affected habitat heterogeneity and population connectivity. Most Amazonian rivers are not associated with long lasting isolation of populations, but some might act as secondary barriers, susceptible to crossing under specific climatic conditions. Morphological variation, while stable in some parts of the distribution, is not a reliable indicator of genetic structure or phylogenetic relationships., (© 2020 John Wiley & Sons Ltd.)

Quaternary climate oscillations are a well-known driver of animal diversification, but their effects are most well studied in areas where glaciations lead to habitat fragmentation. In large areas of the planet, however, glaciations have had the opposite effect, but here their impacts are much less well understood. This is especially true in Southeast Asia, where cyclical changes in land distribution have generated enormous land expansions during glacial periods. In this study, we selected a panel of five songbird species complexes covering a range of ecological specificities to investigate the effects Quaternary land bridges have had on the connectivity of Southeast Asian forest biota. Specifically, we combined morphological and bioacoustic analysis with an arsenal of population genomic and modelling approaches applied to thousands of genome-wide DNA markers across a total of more than 100 individuals. Our analyses show that species dependent on forest understorey exhibit deep differentiation between Borneo and western Sundaland, with no evidence of gene flow during the land bridges accompanying the last 1-2 ice ages. In contrast, dispersive canopy species and habitat generalists have experienced more recent gene flow. Our results argue that there remains much cryptic species-level diversity to be discovered in Southeast Asia even in well-known animal groups such as birds, especially in nondispersive forest understorey inhabitants. We also demonstrate that Quaternary land bridges have not been equally suitable conduits of gene flow for all species complexes and that life history is a major factor in predicting relative population divergence time across Quaternary climate fluctuations., (© 2020 John Wiley & Sons Ltd.)

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome-wide markers to investigate the diversification of the Reunion grey white-eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping-by-sequencing and pooled RAD-seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z-linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process., (© 2020 John Wiley & Sons Ltd.)

Detecting positive selection using genomic data is critical to understanding the role of adaptive evolution. Of particular interest in this context is sex chromosomes since they are thought to play a special role in local adaptation and speciation. We sought to circumvent the challenges associated with statistical phasing when using haplotype-based statistics in sweep scans by benefitting from that whole chromosome haplotypes of the sex chromosomes can be obtained by resequencing of individuals of the hemizygous sex. We analyzed whole Z chromosome haplotypes from 100 females from several populations of four black and white flycatcher species (in birds, females are ZW and males ZZ). Based on integrated haplotype score (iHS) and number of segregating sites by length (nSL) statistics, we found strong and frequent haplotype structure in several regions of the Z chromosome in each species. Most of these sweep signals were population-specific, with essentially no evidence for regions under selection shared among species. Some completed sweeps were revealed by the cross-population extended haplotype homozygosity (XP-EHH) statistic. Importantly, by using statistically phased Z chromosome data from resequencing of males, we failed to recover the signals of selection detected in analyses based on whole chromosome haplotypes from females; instead, what likely represent false signals of selection were frequently seen. This highlights the power issues in statistical phasing and cautions against conclusions from selection scans using such data. The detection of frequent selective sweeps on the avian Z chromosome supports a large role of sex chromosomes in adaptive evolution., (© 2019 The Authors. Molecular Ecology Published by John Wiley & Sons Ltd.)

Understanding migratory connectivity is essential for determining the drivers behind population dynamics and for implementing effective conservation strategies for migratory species. Genetic markers provide a means to describe migratory connectivity; however, they can be uninformative for species with weak population genetic structure, which has limited their application. Here, we demonstrated a genomic approach to describing migratory connectivity in the prothonotary warbler, Protonotaria citrea, a Neotropical songbird of conservation concern. Using 26,189 single nucleotide polymorphisms (SNPs), we revealed regional genetic structure between the Mississippi River Valley and the Atlantic Seaboard with overall weak genetic differentiation among populations (F ST  = 0.0055; 95% CI: 0.0051-0.0059). Genetic variation had a stronger association with geographic rather than environmental factors, with each explaining 14.5% and 8.2% of genetic variation, respectively. By varying the numbers of genomic markers used in population assignment models with individuals of known provenance, we identified a maximum assignment accuracy (89.7% to site, 94.3% to region) using a subset of 600 highly differentiated SNPs. We then assigned samples from nonbreeding sites to breeding region and found low migratory connectivity. Our results highlight the importance of filtering markers for informative loci in models of population assignment. Quantifying migratory connectivity for weakly structured species will be useful for expanding studies to a wider range of migratory species across taxonomic groups and may contribute to a deeper understanding of the evolution of migratory strategies., (© 2019 John Wiley & Sons Ltd.)

Annual reproductive success is often highest in individuals that initiate breeding early, yet relatively few individuals start breeding during this apparently optimal time. This suggests that individuals, particularly females who ultimately dictate when offspring are born, incur costs by initiating reproduction early in the season. We hypothesized that increases in the ageing rate of somatic cells may be one such cost. Telomeres, the repetitive DNA sequences on the ends of chromosomes, may be good proxies of biological wear and tear as they shorten with age and in response to stress. Using historical data from a long-term study population of dark-eyed juncos (Junco hyemalis), we found that telomere loss between years was greater in earlier breeding females, regardless of chronological age. There was no relationship between telomere loss and the annual number of eggs laid or chicks that reached independence. However, telomere loss was greater when temperatures were cooler, and cooler temperatures generally occur early in the season. This suggests that environmental conditions could be the primary cause of accelerated telomere loss in early breeders., (© 2018 John Wiley & Sons Ltd.)

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., (© 2018 John Wiley & Sons Ltd.)

Detailed evaluations of genomic variation between sister species often reveal distinct chromosomal regions of high relative differentiation (i.e., "islands of differentiation" in F ST ), but there is much debate regarding the causes of this pattern. We briefly review the prominent models of genomic islands of differentiation and compare patterns of genomic differentiation in three closely related pairs of New World warblers with the goal of evaluating support for the four models. Each pair (MacGillivray's/mourning warblers; Townsend's/black-throated green warblers; and Audubon's/myrtle warblers) consists of forms that were likely separated in western and eastern North American refugia during cycles of Pleistocene glaciations and have now come into contact in western Canada, where each forms a narrow hybrid zone. We show strong differences between pairs in their patterns of genomic heterogeneity in F ST , suggesting differing selective forces and/or differing genomic responses to similar selective forces among the three pairs. Across most of the genome, levels of within-group nucleotide diversity (π Within ) are almost as large as levels of between-group nucleotide distance (π Between ) within each pair, suggesting recent common ancestry and/or gene flow. In two pairs, a pattern of the F ST peaks having low π Between suggests that selective sweeps spread between geographically differentiated groups, followed by local differentiation. This "sweep-before-differentiation" model is consistent with signatures of gene flow within the yellow-rumped warbler species complex. These findings add to our growing understanding of speciation as a complex process that can involve phases of adaptive introgression among partially differentiated populations., (© 2018 John Wiley & Sons Ltd.)

Colour plays a prominent role in species recognition; therefore, understanding the proximate basis of pigmentation can provide insight into reproductive isolation and speciation. Colour differences between taxa may be the result of regulatory differences or be caused by mutations in coding regions of the expressed genes. To investigate these two alternatives, we studied the pigment composition and the genetic basis of coloration in two divergent dark-eyed junco (Junco hyemalis) subspecies, the slate-coloured and Oregon juncos, which have evolved marked differences in plumage coloration since the Last Glacial Maximum. We used HPLC and light microscopy to investigate pigment composition and deposition in feathers from four body areas. We then used RNA-seq to compare the relative roles of differential gene expression in developing feathers and sequence divergence in transcribed loci under common-garden conditions. Junco feathers differed in eumelanin and pheomelanin content and distribution. Within subspecies, in lighter feathers melanin synthesis genes were downregulated (including PMEL, TYR, TYRP1, OCA2 and MLANA), and ASIP was upregulated. Feathers from different body regions also showed differential expression of HOX and WNT genes. Feathers from the same body regions that differed in colour between the two subspecies showed differential expression of ASIP and three other genes (MFSD12, KCNJ13 and HAND2) associated with pigmentation in other taxa. Sequence variation in the expressed genes was not related to colour differences. Our findings support the hypothesis that differential regulation of a few genes can account for marked differences in coloration, a mechanism that may facilitate the rapid phenotypic diversification of juncos., (© 2018 John Wiley & Sons Ltd.)

Theoretical work suggests that sexual conflict should promote the maintenance of genetic diversity by the opposing directions of selection on males and females. If such conflict is pervasive, it could potentially lead to genomic heterogeneity in levels of genetic diversity an idea that so far has not been empirically tested on a genomewide scale. We used large-scale population genomic and transcriptomic data from the collared flycatcher (Ficedula albicollis) to analyse how sexual conflict, for which we use sex-biased gene expression as a proxy, relates to genetic variability. Here, we demonstrate that the extent of sex-biased gene expression of both male-biased and female-biased genes is significantly correlated with levels of nucleotide diversity in gene sequences and that this correlation extends to diversity levels also in intergenic DNA and introns. We find signatures of balancing selection in sex-biased genes but also note that relaxed purifying selection could potentially explain part of the observed patterns. The finding of significant genetic differentiation between males and females for male-biased (and gonad-specific) genes indicates ongoing sexual conflict and sex-specific viability selection, potentially driven by sexual selection. Our results thus indicate that sexual antagonism could potentially be considered as one viable explanation to the long-standing question in evolutionary biology of how genomes can remain so genetically variable in face of strong natural and sexual selection., (© 2018 John Wiley & Sons Ltd.)

Several hypotheses predict that the major histocompatibility complex (MHC) drives mating preference in females. Olfactory, colour or morphological traits are often found as reliable signals of the MHC profile, but the role of avian song mediating MHC-based female choice remains largely unexplored. We investigated the relationship between several MHC and acoustic features in the collared flycatcher (Ficedula albicollis), a European passerine with complex songs. We screened a fragment of the class IIB second exon of the MHC molecule, of which individuals harbour 4-15 alleles, while considerable sequence diversity is maintained at the population level. To make statistical inferences from a large number of comparisons, we adopted both null-hypothesis testing and effect size framework in combination with randomization procedures. After controlling for potential confounding factors, neither MHC allelic diversity nor the presence of particular alleles was associated remarkably with the investigated qualitative and quantitative song traits. Furthermore, genetic similarity among males based on MHC sequences was not reflected by the similarity in their song based on syllable content. Overall, these results suggest that the relationship between features of song and the allelic composition and diversity of MHC is not strong in the studied species. However, a biologically motivated analysis revealed that individuals that harbour an MHC allele that impairs survival perform songs with broader frequency range. This finding suggests that certain aspects of the song may bear reliable information concerning the MHC profile of the individuals, which can be used by females to optimize mate choice., (© 2018 John Wiley & Sons Ltd.)

Diversifying selection between populations that inhabit different environments can promote lineage divergence within species and ultimately drive speciation. The mitochondrial genome (mitogenome) encodes essential proteins of the oxidative phosphorylation (OXPHOS) system and can be a strong target for climate-driven selection (i.e., associated with inhabiting different climates). We investigated whether Pleistocene climate changes drove mitochondrial selection and evolution within Australian birds. First, using phylogeographic analyses of the mitochondrial ND2 gene for 17 songbird species, we identified mitochondrial clades (mitolineages). Second, using distance-based redundancy analyses, we tested whether climate predicts variation in intraspecific genetic divergence beyond that explained by geographic distances and geographic position. Third, we analysed 41 complete mitogenome sequences representing each mitolineage of 17 species using codon models in a phylogenetic framework and a biochemical approach to identify signals of selection on OXPHOS protein-coding genes and test for parallel selection in mitolineages of different species existing in similar climates. Of 17 species examined, 13 had multiple mitolineages (range: 2-6). Climate was a significant predictor of mitochondrial variation in eight species. At least two amino acid replacements in OXPHOS complex I could have evolved under positive selection in specific mitolineages of two species. Protein homology modelling showed one of these to be in the loop region of the ND6 protein channel and the other in the functionally critical helix HL region of ND5. These findings call for direct tests of the functional and evolutionary significance of mitochondrial protein candidates for climate-associated selection., (© 2018 John Wiley & Sons Ltd.)

Hybrid sterility is a common first step in the evolution of postzygotic reproductive isolation. According to Haldane's Rule, it affects predominantly the heterogametic sex. While the genetic basis of hybrid male sterility in organisms with heterogametic males has been studied for decades, the genetic basis of hybrid female sterility in organisms with heterogametic females has received much less attention. We investigated the genetic basis of reproductive isolation in two closely related avian species, the common nightingale (Luscinia megarhynchos) and the thrush nightingale (L. luscinia), that hybridize in a secondary contact zone and produce viable hybrid progeny. In accordance with Haldane's Rule, hybrid females are sterile, while hybrid males are fertile, allowing gene flow to occur between the species. Using transcriptomic data from multiple individuals of both nightingale species, we identified genomic islands of high differentiation (F ST ) and of high divergence (D xy ), and we analysed gene content and patterns of molecular evolution within these islands. Interestingly, we found that these islands were enriched for genes related to female meiosis and metabolism. The islands of high differentiation and divergence were also characterized by higher levels of linkage disequilibrium than the rest of the genome in both species indicating that they might be situated in genomic regions of low recombination. This study provides one of the first insights into genetic basis of hybrid female sterility in organisms with heterogametic females., (© 2018 John Wiley & Sons Ltd.)

It is becoming increasingly clear that a significant proportion of the functional sequence within eukaryotic genomes is noncoding. However, since the identification of conserved elements (CEs) has been restricted to a limited number of model organisms, the dynamics and evolutionary character of the genomic landscape of conserved, and hence likely functional, sequence is poorly understood in most species. Moreover, identification and analysis of the full suite of functional sequence are particularly important for the understanding of the genetic basis of trait loci identified in genome scans or quantitative trait locus mapping efforts. We report that ~6.6% of the collared flycatcher genome (74.0 Mb) is spanned by ~1.28 million CEs, a higher proportion of the genome but a lower total amount of conserved sequence than has been reported in mammals. We identified >200,000 CEs specific to either the archosaur, avian, neoavian or passeridan lineages, constituting candidates for lineage-specific adaptations. Importantly, no less than ~71% of CE sites were nonexonic (52.6 Mb), and conserved nonexonic sequence density was negatively correlated with functional exonic density at local genomic scales. Additionally, nucleotide diversity was strongly reduced at nonexonic conserved sites (0.00153) relative to intergenic nonconserved sites (0.00427). By integrating deep transcriptome sequencing and additional genome annotation, we identified novel protein-coding genes, long noncoding RNA genes and transposon-derived (exapted) CEs. The approach taken here based on the use of a progressive cactus whole-genome alignment to identify CEs should be readily applicable to nonmodel organisms in general and help to reveal the rich repertoire of putatively functional noncoding sequence as targets for selection., (© 2017 The Authors. Molecular Ecology Published by John Wiley & Sons Ltd.)

Transposable elements (TEs) are genomic parasites capable of inserting virtually anywhere in the host genome, with manifold consequences for gene expression, DNA methylation and genomic stability. Notably, they can contribute to phenotypic variation and hence be associated with, for example, local adaptation and speciation. However, some organisms such as birds have been widely noted for the low densities of TEs in their genomes and this has been attributed to a potential dearth in transposition during their evolution. Here, we show that avian evolution witnessed diverse and abundant transposition on very recent timescales. First, we made an in-depth repeat annotation of the collared flycatcher genome, including identification of 23 new, retrovirus-like LTR retrotransposon families. Then, using whole-genome resequencing data from 200 Ficedula flycatchers, we detected 11,888 polymorphic TE insertions (TE presence/absence variations, TEVs) that segregated within and among species. The density of TEVs was one every 1.5-2.5 Mb per individual, with heterozygosities of 0.12-0.16. The majority of TEVs belonged to some 10 different LTR families, most of which are specific to the flycatcher lineage. TEVs were validated by tracing the segregation of hundreds of TEVs across a three-generation pedigree of collared flycatchers and also by their utility as markers recapitulating the phylogenetic relationships among flycatcher species. Our results suggest frequent germline invasions of songbird genomes by novel retroviruses as a rich source of structural variation, which may have had underappreciated phenotypic consequences for the diversification of this species-rich group of birds., (© 2017 The Authors. Molecular Ecology Published by John Wiley & Sons Ltd.)

Red carotenoid colours in birds are widely assumed to be sexually selected quality indicators, but this rests on a very incomplete understanding of genetic mechanisms and honesty-mediating costs. Recent progress was made by the implication of the gene CYP2J19 as an avian carotenoid ketolase, catalysing the synthesis of red C4-ketocarotenoids from yellow dietary precursors, and potentially a major mechanism behind red coloration in birds. Here, we investigate the role of CYP2J19 in the spectacular colour diversification of African weaverbirds (Ploceidae), represented by five genera and 16 species: eight red, seven yellow and one without carotenoid coloration. All species had a single copy of CYP2J19, unlike the duplication found in the zebra finch, with high expression in the retina, confirming its function in colouring red oil droplets. Expression was weak or undetected in skin and follicles of pigment-depositing feather buds, as well as in beaks and tarsi, including those of the red-billed quelea. In contrast, the hepatic (liver) expression of CYP2J19 was consistently higher (>14-fold) in seven species with C4-ketocarotenoid coloration than in species without (including one red species), an association strongly supported by a phylogenetic comparative analysis. The results suggest a critical role of the candidate ketolase, CYP2J19, in the evolution of red C4-ketocarotenoid colour variation in ploceids. As ancestral state reconstruction suggests that ketocarotenoid coloration has evolved twice in this group (once in Euplectes and once in the Quelea/Foudia clade), we argue that while CYP2J19 has retained its ancestral role in the retina, it has likely been co-opted for red coloration independently in the two lineages, via increased hepatic expression., (© 2017 John Wiley & Sons Ltd.)

Extra-pair paternity (EPP), where offspring are sired by a male other than the social male, varies enormously both within and among species. Trying to explain this variation has proved difficult because the majority of the interspecific variation is phylogenetically based. Ideally, variation in EPP should be investigated in closely related species, but clades with sufficient variation are rare. We present a comprehensive multifactorial test to explain variation in EPP among individuals in 20 populations of nine species over 89 years from a single bird family (Maluridae). Females had higher EPP in the presence of more helpers, more neighbours or if paired incestuously. Furthermore, higher EPP occurred in years with many incestuous pairs, populations with many helpers and species with high male density or in which males provide less care. Altogether, these variables accounted for 48% of the total and 89% of the interspecific and interpopulation variation in EPP. These findings indicate why consistent patterns in EPP have been so challenging to detect and suggest that a single predictor is unlikely to account for the enormous variation in EPP across levels of analysis. Nevertheless, it also shows that existing hypotheses can explain the variation in EPP well and that the density of males in particular is a good predictor to explain variation in EPP among species when a large part of the confounding effect of phylogeny is excluded., (© 2017 John Wiley & Sons Ltd.)

The little greenbul, a common rainforest passerine from sub-Saharan Africa, has been the subject of long-term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest-savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise F ST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation., (© 2017 John Wiley & Sons Ltd.)

Recombination rate is heterogeneous across the genome of various species and so are genetic diversity and differentiation as a consequence of linked selection. However, we still lack a clear picture of the underlying mechanisms for regulating recombination. Here we estimated fine-scale population recombination rate based on the patterns of linkage disequilibrium across the genomes of multiple populations of two closely related flycatcher species (Ficedula albicollis and F. hypoleuca). This revealed an overall conservation of the recombination landscape between these species at the scale of 200 kb, but we also identified differences in the local rate of recombination despite their recent divergence (<1 million years). Genetic diversity and differentiation were associated with recombination rate in a lineage-specific manner, indicating differences in the extent of linked selection between species. We detected 400-3,085 recombination hotspots per population. Location of hotspots was conserved between species, but the intensity of hotspot activity varied between species. Recombination hotspots were primarily associated with CpG islands (CGIs), regardless of whether CGIs were at promoter regions or away from genes. Recombination hotspots were also associated with specific transposable elements (TEs), but this association appears indirect due to shared preferences of the transposition machinery and the recombination machinery for accessible open chromatin regions. Our results suggest that CGIs are a major determinant of the localization of recombination hotspots, and we propose that both the distribution of TEs and fine-scale variation in recombination rate may be associated with the evolution of the epigenetic landscape., (© 2017 The Authors. Molecular Ecology Published by John Wiley & Sons Ltd.)

Differential introgression of mitochondrial vs. nuclear DNA generates discordant patterns of geographic variation and can promote population divergence and speciation. We examined a potential case of mitochondrial introgression leading to two perpendicular axes of differentiation. The Eastern Yellow Robin Eopsaltria australis, a widespread Australian bird, shows a deep mitochondrial split that is perpendicular to north-south nuclear DNA and plumage colour differentiation. We propose a scenario to explain this pattern: (i) first, both nuclear and mitochondrial genomes differentiated in concert during north-south population divergence; (ii) later, their histories disconnected after two mitochondrial introgression events resulting in a deep mitochondrial split perpendicular to the nuclear DNA structure. We explored this scenario by coalescent modelling of ten mitochondrial genes and 400 nuclear DNA loci. Initial mitochondrial and nuclear genome divergences were estimated to have occurred in the early Pleistocene, consistent with the proposed scenario. Subsequent climatic transitions may have driven later mitochondrial introgression. We consider neutral introgression unlikely and instead propose that the evidence is more consistent with adaptive mitochondrial introgression and selection against incompatible mitochondrial-nuclear combinations. This likely generated an axis of coastal-inland mitochondrial differentiation in the face of nuclear gene flow, perpendicular to the initial north-south axis of differentiation (reflected in genomewide nuclear DNA and colour variation)., (© 2017 John Wiley & Sons Ltd.)

Seasonal migration is a widespread phenomenon, which is found in many different lineages of animals. This spectacular behaviour allows animals to avoid seasonally adverse environmental conditions to exploit more favourable habitats. Migration has been intensively studied in birds, which display astonishing variation in migration strategies, thus providing a powerful system for studying the ecological and evolutionary processes that shape migratory behaviour. Despite intensive research, the genetic basis of migration remains largely unknown. Here, we used state-of-the-art radio-tracking technology to characterize the migratory behaviour of a partially migratory population of European blackbirds (Turdus merula) in southern Germany. We compared gene expression of resident and migrant individuals using high-throughput transcriptomics in blood samples. Analyses of sequence variation revealed a nonsignificant genetic structure between blackbirds differing by their migratory phenotype. We detected only four differentially expressed genes between migrants and residents, which might be associated with hyperphagia, moulting and enhanced DNA replication and transcription. The most pronounced changes in gene expression occurred between migratory birds depending on when, in relation to their date of departure, blood was collected. Overall, the differentially expressed genes detected in this analysis may play crucial roles in determining the decision to migrate, or in controlling the physiological processes required for the onset of migration. These results provide new insights into, and testable hypotheses for, the molecular mechanisms controlling the migratory phenotype and its underlying physiological mechanisms in blackbirds and other migratory bird species., (© 2017 John Wiley & Sons Ltd.)