Postglacial colonisation patterns and the role of isolation and expansion in driving diversification in a passerine bird

Evidence of Site-Specific and Male-Biased Germline Mutation Rate in a Wild Songbird

Germline mutations are the ultimate source of genetic variation and the raw material for organismal evolution. Despite their significance, the frequency and genomic locations of mutations, as well as potential sex bias, are yet to be widely investigated in most species. To address these gaps, we conducted whole-genome sequencing of 12 great reed warblers (Acrocephalus arundinaceus) in a pedigree spanning 3 generations to identify single-nucleotide de novo mutations (DNMs) and estimate the germline mutation rate. We detected 82 DNMs within the pedigree, primarily enriched at CpG sites but otherwise randomly located along the chromosomes. Furthermore, we observed a pronounced sex bias in DNM occurrence, with male warblers exhibiting three times more mutations than females. After correction for false negatives and adjusting for callable sites, we obtained a mutation rate of 7.16 × 10-9 mutations per site per generation (m/s/g) for the autosomes and 5.10 × 10-9 m/s/g for the Z chromosome. To demonstrate the utility of species-specific mutation rates, we applied our autosomal mutation rate in models reconstructing the demographic history of the great reed warbler. We uncovered signs of drastic population size reductions predating the last glacial period (LGP) and reduced gene flow between western and eastern populations during the LGP. In conclusion, our results provide one of the few direct estimates of the mutation rate in wild songbirds and evidence for male-driven mutations in accordance with theoretical expectations.

Contrasting mtDNA and microsatellite data of great reed warbler Acrocephalus arundinaceus breeding populations on a small geographic scale

The great reed warbler has two genetically distinguishable haplogroups: "Clade A" occurs in higher proportions in Western Europe and Kazakhstan, and colonised Europe and Asia from a refugium in South-West Europe; and "Clade B", which is more common in Eastern Europe, and colonised parts of Europe from a refugium in the Middle East. Our aims were (i) to analyse the rate of differentiation in Hungarian breeding populations in order to see whether European-scale pattern is visible or not on as a small scale as the territory of Hungary and (ii) to compare the results obtained with mtDNA sequencing and microsatellite markers. To analyse the genetic differentiation, the mtDNA control region II was sequenced in 68 adult breeding birds, and 51 were fingerprinted at 11 microsatellite loci, while both analyses were performed on 36 birds (a total of 83 birds). The microsatellite data gave a better resolution and represented the fine-scale pattern of the suspected recolonisation. The lack of genetic differentiation among the breeding populations based on mitochondrial data seems to support this finding, because the admixture of the clades in this particular geographic region obliterates differentiation. Accordingly, the F_st values from different branches are significantly based on microsatellite data only. The mtDNA methods only give reliable results when a geographic and ecological factor plays a role in the population subdivision, but in the case of an intermixing population larger-scale studies are needed.

The Mitogenome Relationships and Phylogeography of Barn Swallows (Hirundo rustica)

The barn swallow (Hirundo rustica) poses a number of fascinating scientific questions, including the taxonomic status of postulated subspecies. Here we obtained and assessed the sequence variation of 411 complete mitogenomes, mainly from the European H. r. rustica, but other subspecies as well. In almost every case, we observed subspecies-specific haplogroups, which we employed together with estimated radiation times to postulate a model for the geographical and temporal worldwide spread of the species. The female barn swallow carrying the Hirundo rustica ancestral mitogenome left Africa (or its vicinity) around 280 thousand years ago (kya), and her descendants expanded first into Eurasia and then, at least 51 kya, into the Americas, from where a relatively recent (< 20 kya) back migration to Asia took place. The exception to the haplogroup subspecies specificity is represented by the sedentary Levantine H. r. transitiva that extensively shares haplogroup A with the migratory European H. r. rustica and, to a lesser extent, haplogroup B with the Egyptian H. r. savignii. Our data indicate that rustica and transitiva most likely derive from a sedentary Levantine population source that split at the end of the Younger Dryas (11.7 kya). Since then, however, transitiva received genetic inputs from and admixed with both the closely related rustica and the adjacent savignii. Demographic analyses confirm this species' strong link with climate fluctuations and human activities making it an excellent indicator for monitoring and assessing the impact of current global changes on wildlife.

Molecular biogeography of the fungus-dwelling saproxylic beetle Bolitophagus reticulatus indicates rapid expansion from glacial refugia

The geographical distributions of species associated with European temperate broadleaf forests have been significantly influenced by glacial–interglacial cycles. During glacial periods, these species persisted in Mediterranean and extra-Mediterranean refugia and later, during interglacial periods, expanded northwards. The widespread saproxylic beetle Bolitophagus reticulatus depends closely on European temperate broadleaf forests. It usually develops in the tinder fungus Fomes fomentarius, a major decomposer of broadleaf-wood. We sampled B. reticulatus in sporocarps from European beech (Fagus sylvatica) and Oriental beech (Fagus orientalis) across Europe and the Caucasus region. We analysed mitochondrial gene sequences (cox1, cox2, cob) and 17 microsatellites to reconstruct the geographical distribution of glacial refugia and postglacial recolonization pathways. We found only marginal genetic differentiation of B. reticulatus, except for a significant split between populations of the Caucasus region and Europe. This indicates the existence of past refugia south of the Great Caucasus, and a contact zone with European populations in the Crimean region. Further potential refugia might have been located at the foothills of the Pyrenees and in the Balkan region. Our genetic data suggest a phalanx-wise recolonization of Europe, a reflection of the high mobility of B. reticulatus.

Common patterns in the molecular phylogeography of western palearctic birds: a comprehensive review

A plethora of studies have offered crucial insights in the phylogeographic status of Western Palearctic bird species. However, an overview integrating all this information and analyzing the combined results is still missing. In this study, we compiled all published peer-reviewed and grey literature available on the phylogeography of Western Palearctic bird species. Our literature review indicates a total number of 198 studies, with the overwhelming majority published as journal articles (n = 186). In total, these literature items offer information on 145 bird species. 85 of these species are characterized by low genetic differentiation, 46 species indicate genetic variation but no geographic structuring i.e. panmixia, while 14 species show geographically distinct lineages and haplotypes. Majority of bird species inhabiting the Western Palearctic display genetic admixture. The glaciation cycles in the past few million years were pivotal factors in shaping this situation: during warm periods many species expanded their distribution range to the north over wide areas of Eurasia; whereas, during ice ages most areas were no longer suitable and species retreated to refugia, where lineages mixed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10336-021-01893-x.

Unravelling population processes over the Late Pleistocene driving contemporary genetic divergence in Palearctic buzzards

Population range expansions and contractions as a response to climate and habitat change throughout the Quaternary are known to have contributed to complex phylogenetic and population genetic events. Speciation patterns and processes in Palearctic buzzards (genus Buteo) are a long-standing example of morphological and genetic data incongruence, attributed to panmixia, habitat range shifts, contact zones, and climate change. Here we assess the systematics, phylogeography and population genetic structure of three nominal species of Palearctic buzzards, Buteo buteo (including B. b. vulpinus), B. rufinus (including B. r. cirtensis) and B. hemilasius. Phylogenetic analyses inferred from mitochondrial data recover B. hemilasius as sister to the sister clades B. r. rufinus and B. buteo complex (B. b. buteo, B. b. vulpinus, but also including B. r. cirtensis). In contrast, we find an unresolved genetic delimitation inferred from four nuclear loci, suggesting an ancestral genetic pool for all species. Time-trees suggest population contractions and expansions throughout the Pleistocene, which likely reflect habitat change and contrasting ecological niche requirements between species. Microsatellite-based extended Bayesian skyline plots reveal relatively constant population sizes for B. hemilasius, B. r. rufinus, and B. b. vulpinus, in contrast to a dramatic population expansion in B. r. cirtensis within the last 3 kya. Overall, our study illustrates how complex population processes over the Late Pleistocene have shaped the patterns of genetic divergence in Palearctic buzzards, due to the joint effects of shared ancestral polymorphisms, population expansions and contractions, with hybridization at contact zones leading to admixture and introgression.

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.

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

Natural systems composed of closely related taxa that vary in the degree of phenotypic divergence and geographic isolation provide an opportunity to investigate the rate of phenotypic diversification and the relative roles of selection and drift in driving lineage formation. The genus Junco (Aves: Emberizidae) of North America includes parapatric northern forms that are markedly divergent in plumage pattern and colour, in contrast to geographically isolated southern populations in remote areas that show moderate phenotypic divergence. Here, we quantify patterns of phenotypic divergence in morphology and plumage colour and use mitochondrial DNA genes, a nuclear intron, and genomewide SNPs to reconstruct the demographic and evolutionary history of the genus to infer relative rates of evolutionary divergence among lineages. We found that geographically isolated populations have evolved independently for hundreds of thousands of years despite little differentiation in phenotype, in sharp contrast to phenotypically diverse northern forms, which have diversified within the last few thousand years as a result of the rapid postglacial recolonization of North America. SNP data resolved young northern lineages into reciprocally monophyletic lineages, indicating low rates of gene flow even among closely related parapatric forms, and suggesting a role for strong genetic drift or multifarious selection acting on multiple loci in driving lineage divergence. Juncos represent a compelling example of speciation in action, where the combined effects of historical and selective factors have produced one of the fastest cases of speciation known in vertebrates.

Intralocus sexual conflict over wing length in a wild migratory bird

Intralocus sexual conflict (ISC) occurs when males and females have different adaptive peaks but are constrained from evolving sexual dimorphism because of shared genes. Implications of this conflict on evolutionary dynamics in wild populations have not been investigated in detail. In comprehensive analyses of selection, heritability, and genetic correlations, we found evidence for an ISC over wing length, a key trait for flight performance and migration, in a long-term study of wild great reed warblers (Acrocephalus arundinaceus). We found moderate sexual dimorphism, high heritability, moderate sexually antagonistic selection, and strong positive cross-sex genetic correlation in wing length, together supporting the presence of ISC. A negative genetic correlation between male wing length and female fitness indicated that females inheriting alleles for longer wings from their male relatives also inherited lower fitness. Moreover, cross-sex genetic correlations imposed constraint on the predicted microevolutionary trajectory of wing length (based on selection gradients), especially in females where the predicted response was reversed. The degree of sexual dimorphism in wing length did not change over time, suggesting no sign of conflict resolution. Our study provides novel insight into how an ISC over a fitness trait can affect microevolution in a wild population under natural selection.

Genetic patterns of habitat fragmentation and past climate-change effects in the Mediterranean high-mountain plant Armeria caespitosa (Plumbaginaceae)

PREMISE OF THE STUDY

Mountain plants are among the species most vulnerable to global warming, because of their isolation, narrow geographic distribution, and limited geographic range shifts. Stochastic and selective processes can act on the genome, modulating genetic structure and diversity. Fragmentation and historical processes also have a great influence on current genetic patterns, but the spatial and temporal contexts of these processes are poorly known. We aimed to evaluate the microevolutionary processes that may have taken place in Mediterranean high-mountain plants in response to changing historical environmental conditions.

METHODS

Genetic structure, diversity, and loci under selection were analyzed using AFLP markers in 17 populations distributed over the whole geographic range of Armeria caespitosa, an endemic plant that inhabits isolated mountains (Sierra de Guadarrama, Spain). Differences in altitude, geographic location, and climate conditions were considered in the analyses, because they may play an important role in selective and stochastic processes.

KEY RESULTS

Bayesian clustering approaches identified nine genetic groups, although some discrepancies in assignment were found between alternative analyses. Spatially explicit analyses showed a weak relationship between genetic parameters and spatial or environmental distances. However, a large proportion of outlier loci were detected, and some outliers were related to environmental variables.

CONCLUSIONS

A. caespitosa populations exhibit spatial patterns of genetic structure that cannot be explained by the isolation-by-distance model. Shifts along the altitude gradient in response to Pleistocene climatic oscillations and environmentally mediated selective forces might explain the resulting structure and genetic diversity values found.

Phylogeography of a habitat specialist with high dispersal capability: the Savi's Warbler Locustella luscinioides

In order to describe the influence of Pleistocene glaciations on the genetic structure and demography of a highly mobile, but specialized, passerine, the Savi’s Warbler (Locustella luscinioides), mitochondrial DNA sequences (ND2) and microsatellites were analysed in c.330 individuals of 17 breeding and two wintering populations. Phylogenetic, population genetics and coalescent methods were used to describe the genetic structure, determine the timing of the major splits and model the demography of populations. Savi’s Warblers split from its sister species c.8 million years ago and have two major haplotype groups that diverged in the early/middle Pleistocene. One of these clades originated in the Balkans and is currently widespread, showing strong evidence for population expansion; whereas the other is restricted to Iberia and remained stable. Microsatellites agreed with a genetic break around the Pyrenees, but showed considerable introgression and a weaker genetic structure. Both genetic markers showed an isolation-by-distance pattern associated with the population expansion of the eastern clade. Breeding populations seem to be segregated at the wintering sites, but results on migratory connectivity are preliminary. Savi’s Warbler is the only known migratory bird species in which Iberian birds did not expand beyond the Pyrenees after the last glaciation. Despite the long period of independent evolution of western and eastern populations, complete introgression occurred when these groups met in Iberia. Mitochondrial sequences indicated the existence of refugia-within-refugia in the Iberian Peninsula during the last glacial period, which is surprising given the high dispersal capacity of this species. Plumage differences of eastern subspecies seemed to have evolved recently through natural selection, in agreement with the glacial expansion hypothesis. This study supports the great importance of the Iberian Peninsula and its role for the conservation of genetic variation.

Multilocus analysis of the evolutionary dynamics of rainforest bird populations in Southeast Asia

Sundaland has a dynamic geographic history because its landmasses were periodically interconnected when sea levels fell during glacial periods. Superimposed on this geographic dynamism were environmental changes related to climatic oscillations. To investigate how tropical taxa responded to such changes, we studied the divergence and demographic history of two co-distributed rainforest passerine species, Arachnothera longirostra and Malacocincla malaccensis. We sampled birds primarily from Borneo and the Malay Peninsula, which straddle the now-submerged Sunda shelf, and analysed multilocus DNA data with a variety of coalescent and gene genealogy methods. Cross-shelf divergence in both species occurred well before the last glacial maximum, i.e., before the most recent land connection. However, post-divergence gene flow occurred, and it was more pronounced in A. longirostra (a highly vagile nectarivore/insectivore) than in M. malaccensis (an understory insectivore). Despite current habitat continuity on Borneo, the population of M. malaccensis in northeastern Borneo is substantially divergent from that on the rest of the island. The NE population experienced dramatic demographic fluctuations, probably because of competition with the other population, which expanded from western Borneo after the mid-Pleistocene. In contrast, the Borneo population of A. longirostra has little structure and appears to have experienced demographic expansion 16 kya, long after it had diverged from the Malay Peninsula population (630-690 kya). Malay Peninsula populations of both species have remained relatively stable. Overall, the most recent glacial period was not the chief determinant of the evolutionary dynamics of our study species, and in this respect, they are different from temperate species.

Environmental and climatic determinants of molecular diversity and genetic population structure in a coenagrionid damselfly

Identifying environmental factors that structure intraspecific genetic diversity is of interest for both habitat preservation and biodiversity conservation. Recent advances in statistical and geographical genetics make it possible to investigate how environmental factors affect geographic organisation and population structure of molecular genetic diversity within species. Here we present a study on a common and wide ranging insect, the blue tailed damselfly Ischnuraelegans, which has been the target of many ecological and evolutionary studies. We addressed the following questions: (i) Is the population structure affected by longitudinal or latitudinal gradients?; (ii) Do geographic boundaries limit gene flow?; (iii) Does geographic distance affect connectivity and is there a signature of past bottlenecks?; (iv) Is there evidence of a recent range expansion and (vi) what is the effect of geography and climatic factors on population structure? We found low to moderate genetic sub-structuring between populations (mean F(ST) = 0.06, D(est) = 0.12), and an effect of longitude, but not latitude, on genetic diversity. No significant effects of geographic boundaries (e.g. water bodies) were found. F(ST)-and D(est)-values increased with geographic distance; however, there was no evidence for recent bottlenecks. Finally, we did not detect any molecular signatures of range expansions or an effect of geographic suitability, although local precipitation had a strong effect on genetic differentiation. The population structure of this small insect has probably been shaped by ecological factors that are correlated with longitudinal gradients, geographic distances, and local precipitation. The relatively weak global population structure and high degree of genetic variation within populations suggest that I. elegans has high dispersal ability, which is consistent with this species being an effective and early coloniser of new habitats.

Mutation rate is linked to diversification in birds

How does genome evolution affect the rate of diversification of biological lineages? Recent studies have suggested that the overall rate of genome evolution is correlated with the rate of diversification. If true, this claim has important consequences for understanding the process of diversification, and implications for the use of DNA sequence data to reconstruct evolutionary history. However, the generality and cause of this relationship have not been established. Here, we test the relationship between the rate of molecular evolution and net diversification with a 19-gene, 17-kb DNA sequence dataset from 64 families of birds. We show that rates of molecular evolution are positively correlated to net diversification in birds. Using a 7.6-kb dataset of protein-coding DNA, we show that the synonymous substitution rate, and therefore the mutation rate, is correlated to net diversification. Further analysis shows that the link between mutation rates and net diversification is unlikely to be the indirect result of correlations with life-history variables that may influence both quantities, suggesting that there might be a causal link between mutation rates and net diversification.

Influence of quaternary sea-level variations on a land bird endemic to Pacific atolls

Little is known about the effect of quaternary climate variations on organisms that inhabited carbonate islands of the Pacific Ocean, although it has been suggested that one or several uplifted islands provided shelter for terrestrial birds when sea-level reached its highest. To test this hypothesis, we investigated the history of colonization of the Tuamotu reed-warbler (Acrocephalus atyphus) in southeastern Polynesia, and found high genetic structure between the populations of three elevated carbonate islands. Estimates of time since divergence support the hypothesis that these islands acted as refugia during the last interglacial maximum. These findings are particularly important for defining conservation priorities on atolls that endure the current trend of sea-level rise owing to global warming.

Population dynamic of the extinct European aurochs: genetic evidence of a north-south differentiation pattern and no evidence of post-glacial expansion

Background

The aurochs (Bos primigenius) was a large bovine that ranged over almost the entirety of the Eurasian continent and North Africa. It is the wild ancestor of the modern cattle (Bos taurus), and went extinct in 1627 probably as a consequence of human hunting and the progressive reduction of its habitat. To investigate in detail the genetic history of this species and to compare the population dynamics in different European areas, we analysed Bos primigenius remains from various sites across Italy.

Results

Fourteen samples provided ancient DNA fragments from the mitochondrial hypervariable region. Our data, jointly analysed with previously published sequences, support the view that Italian aurochsen were genetically similar to modern bovine breeds, but very different from northern/central European aurochsen. Bayesian analyses and coalescent simulations indicate that the genetic variation pattern in both Italian and northern/central European aurochsen is compatible with demographic stability after the last glaciation. We provide evidence that signatures of population expansion can erroneously arise in stable aurochsen populations when the different ages of the samples are not taken into account.

Conclusions

Distinct groups of aurochsen probably inhabited Italy and northern/central Europe after the last glaciation, respectively. On the contrary, Italian and Fertile Crescent aurochsen likely shared several mtDNA sequences, now common in modern breeds. We argue that a certain level of genetic homogeneity characterized aurochs populations in Southern Europe and the Middle East, and also that post-glacial recolonization of northern and central Europe advanced, without major demographic expansions, from eastern, and not southern, refugia.

A strong quantitative trait locus for wing length on chromosome 2 in a wild population of great reed warblers

Wing length is a key character for essential behaviours related to bird flight such as migration and foraging. In the present study, we initiate the search for the genes underlying wing length in birds by studying a long-distance migrant, the great reed warbler (Acrocephalus arundinaceus). In this species wing length is an evolutionary interesting trait with pronounced latitudinal gradient and sex-specific selection regimes in local populations. We performed a quantitative trait locus (QTL) scan for wing length in great reed warblers using phenotypic, genotypic, pedigree and linkage map data from our long-term study population in Sweden. We applied the linkage analysis mapping method implemented in GridQTL (a new web-based software) and detected a genome-wide significant QTL for wing length on chromosome 2, to our knowledge, the first detected QTL in wild birds. The QTL extended over 25 cM and accounted for a substantial part (37%) of the phenotypic variance of the trait. A genome scan for tarsus length (a body-size-related trait) did not show any signal, implying that the wing-length QTL on chromosome 2 was not associated with body size. Our results provide a first important step into understanding the genetic architecture of avian wing length, and give opportunities to study the evolutionary dynamics of wing length at the locus level.

Egg eviction imposes a recoverable cost of virulence in chicks of a brood parasite

Background

Chicks of virulent brood parasitic birds eliminate their nestmates and avoid costly competition for foster parental care. Yet, efforts to evict nest contents by the blind and naked common cuckoo Cuculus canorus hatchling are counterintuitive as both adult parasites and large older cuckoo chicks appear to be better suited to tossing the eggs and young of the foster parents.

Methodology/Principal Findings

Here we show experimentally that egg tossing imposed a recoverable growth cost of mass gain in common cuckoo chicks during the nestling period in nests of great reed warbler Acrocephalus arundinaceus hosts. Growth rates of skeletal traits and morphological variables involved in the solicitation of foster parental care remained similar between evictor and non-evictor chicks throughout development. We also detected no increase in predation rates for evicting nests, suggesting that egg tossing behavior by common cuckoo hatchlings does not increase the conspicuousness of nests.

Conclusion

The temporary growth cost of egg eviction by common cuckoo hatchlings is the result of constraints imposed by rejecter host adults and competitive nestmates on the timing and mechanism of parasite virulence.