Genomics of natural bird populations: a gene-based set of reference markers evenly spread across the avian genome

Fine-Scale Map Reveals Highly Variable Recombination Rates Associated with Genomic Features in the Eurasian Blackcap

Recombination is responsible for breaking up haplotypes, influencing genetic variability, and the efficacy of selection. Bird genomes lack the protein PR domain-containing protein 9, a key determinant of recombination dynamics in most metazoans. Historical recombination maps in birds show an apparent stasis in positioning recombination events. This highly conserved recombination pattern over long timescales may constrain the evolution of recombination in birds. At the same time, extensive variation in recombination rate is observed across the genome and between different species of birds. Here, we characterize the fine-scale historical recombination map of an iconic migratory songbird, the Eurasian blackcap (Sylvia atricapilla), using a linkage disequilibrium-based approach that accounts for population demography. Our results reveal variable recombination rates among and within chromosomes, which associate positively with nucleotide diversity and GC content and negatively with chromosome size. Recombination rates increased significantly at regulatory regions but not necessarily at gene bodies. CpG islands are associated strongly with recombination rates, though their specific position and local DNA methylation patterns likely influence this relationship. The association with retrotransposons varied according to specific family and location. Our results also provide evidence of heterogeneous intrachromosomal conservation of recombination maps between the blackcap and its closest sister taxon, the garden warbler. These findings highlight the considerable variability of recombination rates at different scales and the role of specific genomic features in shaping this variation. This study opens the possibility of further investigating the impact of recombination on specific population-genomic features.

X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

UCEasy: A software package for automating and simplifying the analysis of ultraconserved elements (UCEs)

Background

The use of Ultraconserved Elements (UCEs) as genetic markers in phylogenomics has become popular and has provided promising results. Although UCE data can be easily obtained from targeted enriched sequencing, the protocol for in silico analysis of UCEs consist of the execution of heterogeneous and complex tools, a challenge for scientists without training in bioinformatics. Developing tools with the adoption of best practices in research software can lessen this problem by improving the execution of computational experiments, thus promoting better reproducibility.

New information

We present UCEasy, an easy-to-install and easy-to-use software package with a simple command line interface that facilitates the computational analysis of UCEs from sequencing samples, following the best practices of research software. UCEasy is a wrapper that standardises, automates and simplifies the quality control of raw reads, assembly and extraction and alignment of UCEs, generating at the end a data matrix with different levels of completeness that can be used to infer phylogenetic trees. We demonstrate the functionalities of UCEasy by reproducing the published results of phylogenomic studies of the bird genus Turdus (Aves) and of Adephaga families (Coleoptera) containing genomic datasets to efficiently extract UCEs.

A global phylogeny of turtles reveals a burst of climate-associated diversification on continental margins

Living turtles are characterized by extraordinarily low species diversity given their age. The clade's extensive fossil record indicates that climate and biogeography may have played important roles in determining their diversity. We investigated this hypothesis by collecting a molecular dataset for 591 individual turtles that, together, represent 80% of all turtle species, including representatives of all families and 98% of genera, and used it to jointly estimate phylogeny and divergence times. We found that the turtle tree is characterized by relatively constant diversification (speciation minus extinction) punctuated by a single threefold increase. We also found that this shift is temporally and geographically associated with newly emerged continental margins that appeared during the Eocene-Oligocene transition about 30 million years before present. In apparent contrast, the fossil record from this time period contains evidence for a major, but regional, extinction event. These seemingly discordant findings appear to be driven by a common global process: global cooling and drying at the time of the Eocene-Oligocene transition. This climatic shift led to aridification that drove extinctions in important fossil-bearing areas, while simultaneously exposing new continental margin habitat that subsequently allowed for a burst of speciation associated with these newly exploitable ecological opportunities.

Complete mitochondrial genome of Chroicocephalus brunnicephalus from India: phylogeny with other Larids

The complete mitogenome sequence of the brown-headed gull, Chroicocephalus brunnicephalus was determined in this study. The 16,771 bp genome consists of 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes, and a control region (CR). The decoded mitogenome was AT-rich (54.77%) with nine overlapping and 17 intergenic spacer regions. Most of the PCGs were started by a typical ATG initiation codon except for cox1 and nad3. Further, the usual termination codons (AGG, TAG, TAA, and AGA) were used by 11 PCGs except for cox3 and nad4. The concatenated PCGs based Bayesian phylogeny clearly discriminates all the Laridae species and reflects the sister relationship of C. brunnicephalus with C. ridibundus. The present mitogenome-based phylogeny was congruent with the earlier hypothesis and confirmed the evolutionary position of the brown-headed gull as masked species. The generated mitogenome of C. brunnicephalus is almost identical to the previously generated mitogenome from China except for two base pairs in CR. To visualize the population structure of this migratory species, we propose more sampling from different geographical locations and the generation of additional molecular data to clarify the reality.

Phylogeography, Population Structure, and Species Delimitation in Rockhopper Penguins (Eudyptes chrysocome and Eudyptes moseleyi)

Rockhopper penguins are delimited as 2 species, the northern rockhopper (Eudyptes moseleyi) and the southern rockhopper (Eudyptes chrysocome), with the latter comprising 2 subspecies, the western rockhopper (Eudyptes chrysocome chrysocome) and the eastern rockhopper (Eudyptes chrysocome filholi). We conducted a phylogeographic study using multilocus data from 114 individuals sampled across 12 colonies from the entire range of the northern/southern rockhopper complex to assess potential population structure, gene flow, and species limits. Bayesian and likelihood methods with nuclear and mitochondrial DNA, including model testing and heuristic approaches, support E. moseleyi and E. chrysocome as distinct species lineages with a divergence time of 0.97 Ma. However, these analyses also indicated the presence of gene flow between these species. Among southern rockhopper subspecies, we found evidence of significant gene flow and heuristic approaches to species delimitation based on the genealogical diversity index failed to delimit them as species. The best-supported population models for the southern rockhoppers were those where E. c. chrysocome and E. c. filholi were combined into a single lineage or 2 lineages with bidirectional gene flow. Additionally, we found that E. c. filholi has the highest effective population size while E. c. chrysocome showed similar effective population size to that of the endangered E. moseleyi. We suggest that the current taxonomic definitions within rockhopper penguins be upheld and that E. chrysocome populations, all found south of the subtropical front, should be treated as a single taxon with distinct management units for E. c. chrysocome and E. c. filholi.

Phylogenomics and biogeography of the world's thrushes (Aves, Turdus): new evidence for a more parsimonious evolutionary history

To elucidate the relationships and spatial range evolution across the world of the bird genus Turdus (Aves), we produced a large genomic dataset comprising ca 2 million nucleotides for ca 100 samples representing 53 species, including over 2000 loci. We estimated time-calibrated maximum-likelihood and multispecies coalescent phylogenies and carried out biogeographic analyses. Our results indicate that there have been considerably fewer trans-oceanic dispersals within the genus Turdus than previously suggested, such that the Palaearctic clade did not originate in America and the African clade was not involved in the colonization of the Americas. Instead, our findings suggest that dispersal from the Western Palaearctic via the Antilles to the Neotropics might have occurred in a single event, giving rise to the rich Neotropical diversity of Turdus observed today, with no reverse dispersals to the Palaearctic or Africa. Our large multilocus dataset, combined with dense species-level sampling and analysed under probabilistic methods, brings important insights into historical biogeography and systematics, even in a scenario of fast and spatially complex diversification.

High genetic diversity and low differentiation retained in the European fragmented and declining Greater Spotted Eagle (Clanga clanga) population

Characterising genetic diversity and structure of populations is essential for effective conservation of threatened species. The Greater Spotted Eagle (Clanga clanga), a large and globally vulnerable raptor, is extinct or in severe decline in most of its previous range in Europe. We assessed whether the remnants of European population are genetically impoverished, and isolated from each other. We evaluated levels of genetic diversity and population structuring by sequencing mitochondrial pseudo-control region and 10 introns from various nuclear genes, and estimated length diversity in 23 microsatellite markers. The European population has expanded since the late Pleistocene, and does not exhibit signs of a recent population bottleneck. The global genetic diversity in Europe was rather similar to that detected in other similar species. Microsatellites suggested shallow but significant differentiation between the four extant populations in Estonia, Poland, Belarus and Russia (Upper Volga region) populations, but introns and mtDNA showed that only the Estonian population differed from the others. Mitochondrial diversity was highest in the northernmost Estonian population, introns suggested lower diversity in Upper Volga, microsatellites indicated equal diversity among populations. A recent bottleneck was detected in Poland, which is consistent with the observed repopulation of the region. We conclude that significant gene flow and high genetic diversity are retained in the fragmented Greater Spotted Eagle populations; there is currently no need for genetic augmentation in Europe.

The niches of nuthatches affect their lineage evolution differently across latitude

Ecological niche evolution can promote or hinder the differentiation of taxa and determine their distribution. Niche-mediated evolution may differ among climatic regimes, and thus, species that occur across a wide latitudinal range offer a chance to test these heterogeneous evolutionary processes. In this study, we examine (a) how many lineages have evolved across the continent-wide range of the Eurasian nuthatch (Sitta europaea), (b) whether the lineages' niches are significantly divergent or conserved and (c) how their niche evolution explains their geographic distribution. Phylogenetic reconstruction and ecological niche models (ENMs) showed that the Eurasian nuthatch contained six parapatric lineages that diverged within 2 Myr and did not share identical climatic niches. However, the niche discrepancy between these distinct lineages was relatively conserved compared with the environmental differences between their ranges and thus was unlikely to drive lineage divergence. The ENMs of southern lineages tended to cross-predict with their neighbouring lineages whereas those of northern lineages generally matched with their abutting ranges. The coalescence-based analyses revealed more stable populations for the southern lineages than the northern ones during the last glaciation cycle. In contrast to the overlapping ENMs, the smaller parapatric distribution suggests that the southern lineages might have experienced competitive exclusion to prevent them from becoming sympatric. On the other hand, the northern lineages have expanded their ranges and their current abutting distribution might have resulted from lineages adapting to different climatic conditions in allopatry. This study suggests that niche evolution may affect lineage distribution in different ways across latitude.

Postcopulatory sexual selection reduces Z-linked genetic variation and might contribute to the large Z effect in passerine birds

The X and Z sex chromosomes play a disproportionately large role in intrinsic postzygotic isolation. The underlying mechanisms of this large X/Z effect are, however, still poorly understood. Here we tested whether faster rates of molecular evolution caused by more intense positive selection or genetic drift on the Z chromosome could contribute to the large Z effect in two closely related passerine birds, the Common Nightingale (Luscinia megarhynchos) and the Thrush Nightingale (L. luscinia). We found that the two species differ in patterns of molecular evolution on the Z chromosome. The Z chromosome of L. megarhynchos showed lower levels of within-species polymorphism and an excess of non-synonymous polymorphisms relative to non-synonymous substitutions. This is consistent with increased levels of genetic drift on this chromosome and may be attributed to more intense postcopulatory sexual selection acting on L. megarhynchos males as was indicated by significantly longer sperm and higher between-male variation in sperm length in L. megarhynchos compared to L. luscinia. Interestingly, analysis of interspecific gene flow on the Z chromosome revealed relatively lower levels of introgression from L. megarhynchos to L. luscinia than vice versa, indicating that the Z chromosome of L. megarhynchos accumulated more hybrid incompatibilities. Our results are consistent with the view that postcopulatory sexual selection may reduce the effective population size of the Z chromosome and thus lead to stronger genetic drift on this chromosome in birds. This can result in relatively faster accumulation of hybrid incompatibilities on the Z and thus contribute to the large Z effect.

The first set of universal nuclear protein-coding loci markers for avian phylogenetic and population genetic studies

Multiple nuclear markers provide genetic polymorphism data for molecular systematics and population genetic studies. They are especially required for the coalescent-based analyses that can be used to accurately estimate species trees and infer population demographic histories. However, in avian evolutionary studies, these powerful coalescent-based methods are hindered by the lack of a sufficient number of markers. In this study, we designed PCR primers to amplify 136 nuclear protein-coding loci (NPCLs) by scanning the published Red Junglefowl (Gallus gallus) and Zebra Finch (Taeniopygia guttata) genomes. To test their utility, we amplified these loci in 41 bird species representing 23 Aves orders. The sixty-three best-performing NPCLs, based on high PCR success rates, were selected which had various mutation rates and were evenly distributed across 17 avian autosomal chromosomes and the Z chromosome. To test phylogenetic resolving power of these markers, we conducted a Neoavian phylogenies analysis using 63 concatenated NPCL markers derived from 48 whole genomes of birds. The resulting phylogenetic topology, to a large extent, is congruence with results resolved by previous whole genome data. To test the level of intraspecific polymorphism in these makers, we examined the genetic diversity in four populations of the Kentish Plover (Charadrius alexandrinus) at 17 of NPCL markers chosen at random. Our results showed that these NPCL markers exhibited a level of polymorphism comparable with mitochondrial loci. Therefore, this set of pan-avian nuclear protein-coding loci has great potential to facilitate studies in avian phylogenetics and population genetics.

Developing genome-wide SNPs and constructing an ultrahigh-density linkage map in oil palm

Oil palm (Elaeis guineensis Jacq.) is the leading oil-producing crops and the most important edible oil resource worldwide. DNA markers and genetic linkage maps are essential resources for marker-assisted selection to accelerate genetic improvement. We conducted RAD-seq on an Illumina NextSeq500 to discover genome-wide SNPs, and used the SNPs to construct a linkage map for an oil palm (Tenera) population derived from a cross between a Deli Dura and an AVROS Pisifera. The RAD-seq produced 1,076 million single-end reads across the breeding population containing 155 trees. Mining this dataset detected 510,251 loci. After filtering out loci with low accuracy and more than 20% missing data, 11,394 SNPs were retained. Using these SNPs, in combination with 188 anchor SNPs and 123 microsatellites, we constructed a linkage map containing 10,023 markers covering 16 chromosomes. The map length is 2,938.2 cM with an average marker space of 0.29 cM. The large number of SNPs will supply ample choices of DNA markers in analysing the genetic diversity, population structure and evolution of oil palm. This high-density linkage map will contribute to mapping quantitative trait loci (QTL) for important traits, thus accelerating oil palm genetic improvement.

Molecular phylogeny and divergence of the map turtles (Emydidae: Graptemys)

The map turtles (genus Graptemys) comprise a morphologically diverse clade that forms a major component of the southeastern US hotspot of chelonian diversity. Map turtles have experienced both recent and rapid diversification resulting in long-standing uncertainty regarding species boundaries and phylogenetic relationships within the genus as well as timing of their divergence. We present a phylogeny for the group that includes geographically representative sampling for all described species and subspecies. We make use of an empirical prior on rates of molecular evolution to estimate divergence times with a molecular clock under a coalescent framework. Together, the phylogeny and divergence time estimates suggest that diversification has been both more recent and more rapid than has so far been suspected. We provide a well-supported evolutionary framework for Graptemys that is necessary for understanding map turtle diversity, biogeography, and for conservation of this threatened clade of turtles.

Ice age unfrozen: severe effect of the last interglacial, not glacial, climate change on East Asian avifauna

BACKGROUND

The glacial-interglacial cycles in the Pleistocene caused repeated range expansion and contraction of species in several regions in the world. However, it remains uncertain whether such climate oscillations had similar impact on East Asian biota, despite its widely recognized importance in global biodiversity. Here we use both molecular and ecological niche profiles on 11 East Asian avian species with various elevational ranges to reveal their response to the late Pleistocene climate changes.

RESULTS

The ecological niche models (ENM) consistently showed that these avian species might substantially contract their ranges to the south during the Last Interglacial period (LIG) and expanded their northern range margins through the Last Glacial Maximum (LGM), leading to the LGM ranges observed for all 11 species. Consistently, coalescent simulations based on 25-30 nuclear genes retrieved signatures of significant population growth through the last glacial period across all species studied. Climate statistics suggested that high climatic variability during the LIG and a relatively mild climate at the LGM potentially explained the historical population dynamics of these birds.

CONCLUSIONS

This is the first study based on multiple species and both lines of ecological niche profiles and genetic data to characterize the unique response of East Asian biota to late Pleistocene climate. The present study highlights regional differences in the evolutionary consequence of climate change during the last glacial cycle and implies that global warming might pose a great risk to species in this region given potentially higher climatic variation in the future analogous to that during the LIG.

A workflow of massive identification and application of intron markers using snakes as a model

Relative to the commonly used mitochondrial and nuclear protein-coding genes, the noncoding intron sequences are a promising source of informative markers that have the potential to resolve difficult phylogenetic nodes such as rapid radiations and recent divergences. Yet many issues exist in the use of intron markers, which prevent their extensive application as conventional markers. We used the diverse group of snakes as an example to try paving the way for massive identification and application of intron markers. We performed a series of bioinformatics screenings which identified appropriate introns between single-copy and conserved exons from two snake genomes, adding particular constraints on sequence length variability and sequence variability. A total of 1,273 candidate intron loci were retrieved. Primers for nested polymerase chain reaction (PCR) were designed for over a hundred candidates and tested in 16 snake representatives. 96 intron markers were developed that could be amplified across a broad range of snake taxa with high PCR successful rates. The markers were then applied to 49 snake samples. The large number of amplicons was subjected to next-generation sequencing (NGS). An analytic strategy was developed to accurately recover the amplicon sequences, and approximately, 76% of the marker sequences were recovered. The average p-distances of the intron markers at interfamily, intergenus, interspecies, and intraspecies levels were .168, .052, .015, and .004, respectively, suggesting that they were useful to study snake relationships of different evolutionary depths. A snake phylogeny was constructed with the intron markers, which produced concordant results with robust support at both interfamily and intragenus levels. The intron markers provide a convenient way to explore the signals in the noncoding regions to address the controversies on the snake tree. Our improved strategy of genome screening is effective and can be applied to other animal groups. NGS coupled with appropriate sequence processing can greatly facilitate the extensive application of molecular markers.

Broad-scale recombination pattern in the primitive bird Rhea americana (Ratites, Palaeognathae)

Birds have genomic and chromosomal features that make them an attractive group to analyze the evolution of recombination rate and the distribution of crossing over. Yet, analyses are biased towards certain species, especially domestic poultry and passerines. Here we analyze for the first time the recombination rate and crossover distribution in the primitive ratite bird, Rhea americana (Rheiformes, Palaeognathae). Using a cytogenetic approach for in situ mapping of crossovers we found that the total genetic map is 3050 cM with a global recombination rate of 2.1 cM/Mb for female rheas. In the five largest macrobivalents there were 3 or more crossovers in most bivalents. Recombination rates for macrobivalents ranges between 1.8-2.1 cM/Mb and the physical length of their synaptonemal complexes is highly predictive of their genetic lengths. The crossover rate at the pseudoautosomal region is 2.1 cM/Mb, similar to those of autosomal pairs 5 and 6 and only slightly higher compared to other macroautosomes. It is suggested that the presence of multiple crossovers on the largest macrobivalents is a feature common to many avian groups, irrespective of their position throughout phylogeny. These data provide new insights to analyze the heterogeneous recombination landscape of birds.

Gene flow during glacial habitat shifts facilitates character displacement in a Neotropical flycatcher radiation

BACKGROUND

Pleistocene climatic fluctuations are known to be an engine of biotic diversification at higher latitudes, but their impact on highly diverse tropical areas such as the Andes remains less well-documented. Specifically, while periods of global cooling may have led to fragmentation and differentiation at colder latitudes, they may - at the same time - have led to connectivity among insular patches of montane tropical habitat with unknown consequences on diversification. In the present study we utilized ~5.5 kb of DNA sequence data from eight nuclear loci and one mitochondrial gene alongside diagnostic morphological and bioacoustic markers to test the effects of Pleistocene climatic fluctuations on diversification in a complex of Andean tyrant-flycatchers of the genus Elaenia.

RESULTS

Population genetic and phylogenetic approaches coupled with coalescent simulations demonstrated disparate levels of gene flow between the taxon chilensis and two parapatric Elaenia taxa predominantly during the last glacial period but not thereafter, possibly on account of downward shifts of montane forest habitat linking the populations of adjacent ridges. Additionally, morphological and bioacoustic analyses revealed a distinct pattern of character displacement in coloration and vocal traits between the two sympatric taxa albiceps and pallatangae, which were characterized by a lack of gene flow.

CONCLUSION

Our study demonstrates that global periods of cooling are likely to have facilitated gene flow among Andean montane Elaenia flycatchers that are more isolated from one another during warm interglacial periods such as the present era. We also identify a hitherto overlooked case of plumage and vocal character displacement, underpinning the complexities of gene flow patterns caused by Pleistocene climate change across the Andes.

Matching loci surveyed to questions asked in phylogeography

Although mitochondrial DNA (mtDNA) has long been used for assessing genetic variation within and between populations, its workhorse role in phylogeography has been criticized owing to its single-locus nature. The only choice for testing mtDNA results is to survey nuclear loci, which brings into contrast the difference in locus effective size and coalescence times. Thus, it remains unclear how erroneous mtDNA-based estimates of species history might be, especially for evolutionary events in the recent past. To test the robustness of mtDNA and nuclear sequences in phylogeography, we provide one of the largest paired comparisons of summary statistics and demographic parameters estimated from mitochondrial, five Z-linked and 10 autosomal genes of 30 avian species co-distributed in the Caucasus and Europe. The results suggest that mtDNA is robust in estimating inter-population divergence but not in intra-population diversity, which is sensitive to population size change. Here, we provide empirical evidence showing that mtDNA was more likely to detect population divergence than any other single locus owing to its smaller Ne and thus faster coalescent time. Therefore, at least in birds, numerous studies that have based their inferences of phylogeographic patterns solely on mtDNA should not be readily dismissed.

Multi-locus sequence data illuminate demographic drivers of Pleistocene speciation in semi-arid southern Australian birds (Cinclosoma spp.)

BACKGROUND

During the Pleistocene, shifts of species distributions and their isolation in disjunct refugia led to varied outcomes in how taxa diversified. Some species diverged, others did not. Here, we begin to address another facet of the role of the Pleistocene in generating today's diversity. We ask which processes contributed to divergence in semi-arid southern Australian birds. We isolated 11 autosomal nuclear loci and one mitochondrial locus from a total of 29 specimens of the sister species pair, Chestnut Quail-thrush Cinclosoma castanotum and Copperback Quail-thrush C. clarum.

RESULTS

A population clustering analysis confirmed the location of the current species boundary as a well-known biogeographical barrier in southern Australia, the Eyrean Barrier. Coalescent-based analyses placed the time of species divergence to the Middle Pleistocene. Gene flow between the species since divergence has been low. The analyses suggest the effective population size of the ancestor was 54 to 178 times smaller than populations since divergence. This contrasts with recent multi-locus studies in some other Australian birds (butcherbirds, ducks) where a lack of phenotypic divergence was accompanied by larger historical population sizes. Post-divergence population size histories of C. clarum and C. castanotum were inferred using the extended Bayesian skyline model. The population size of C. clarum increased substantially during the late Pleistocene and continued to increase through the Last Glacial Maximum and Holocene. The timing of this expansion across its vast range is broadly concordant with that documented in several other Australian birds. In contrast, effective population size of C. castanotum was much more constrained and may reflect its smaller range and more restricted habitat east of the Eyrean Barrier compared with that available to C. clarum to the west.

CONCLUSIONS

Our results contribute to awareness of increased population sizes, following significant contractions, as having been important in shaping diversity in Australian arid and semi-arid zones. Further, we improve knowledge of the role of Pleistocene climatic shifts in areas of the planet that were not glaciated at that time but which still experienced that period's cyclical climatic fluctuations.

Social selection parapatry in Afrotropical sunbirds

The extent of range overlap of incipient and recent species depends on the type and magnitude of phenotypic divergence that separates them, and the consequences of phenotypic divergence on their interactions. Signal divergence by social selection likely initiates many speciation events, but may yield niche-conserved lineages predisposed to limit each others' ranges via ecological competition. Here, we examine this neglected aspect of social selection speciation theory in relation to the discovery of a nonecotonal species border between sunbirds. We find that Nectarinia moreaui and Nectarinia fuelleborni meet in a ∼6 km wide contact zone, as estimated by molecular cline analysis. These species exploit similar bioclimatic niches, but sing highly divergent learned songs, consistent with divergence by social selection. Cline analyses suggest that within-species stabilizing social selection on song-learning predispositions maintains species differences in song despite both hybridization and cultural transmission. We conclude that ecological competition between moreaui and fuelleborni contributes to the stabilization of the species border, but that ecological competition acts in conjunction with reproductive interference. The evolutionary maintenance of learned song differences in a hybrid zone recommend this study system for future studies on the mechanisms of learned song divergence and its role in speciation.

Genetic consequences of population expansions and contractions in the common hippopotamus (Hippopotamus amphibius) since the Late Pleistocene

Over the past two decades, an increasing amount of phylogeographic work has substantially improved our understanding of African biogeography, in particular the role played by Pleistocene pluvial-drought cycles on terrestrial vertebrates. However, still little is known on the evolutionary history of semi-aquatic animals, which faced tremendous challenges imposed by unpredictable availability of water resources. In this study, we investigate the Late Pleistocene history of the common hippopotamus (Hippopotamus amphibius), using mitochondrial and nuclear DNA sequence variation and range-wide sampling. We documented a global demographic and spatial expansion approximately 0.1-0.3 Myr ago, most likely associated with an episode of massive drainage overflow. These events presumably enabled a historical continent-wide gene flow among hippopotamus populations, and hence, no clear continental-scale genetic structuring remains. Nevertheless, present-day hippopotamus populations are genetically disconnected, probably as a result of the mid-Holocene aridification and contemporary anthropogenic pressures. This unique pattern contrasts with the biogeographic paradigms established for savannah-adapted ungulate mammals and should be further investigated in other water-associated taxa. Our study has important consequences for the conservation of the hippo, an emblematic but threatened species that requires specific protection to curtail its long-term decline.

Developing Exon-Primed Intron-Crossing (EPIC) markers for population genetic studies in three Aedes disease vectors

Aedes aegypti, Aedes notoscriptus, and Aedes albopictus are important vectors of many arboviruses implicated in human disease such as dengue fever. Genetic markers applied across vector species can provide important information on population structure, gene flow, insecticide resistance, and taxonomy, however, robust microsatellite markers have proven difficult to develop in these species and mosquitoes generally. Here we consider the utility and transferability of 15 Ribosome protein (Rp) Exon-Primed Intron-Crossing (EPIC) markers for population genetic studies in these 3 Aedes species. Rp EPIC markers designed for Ae. aegypti also successfully amplified populations of the sister species, Ae. albopictus, as well as the distantly related species, Ae. notoscriptus. High SNP and good indel diversity in sequenced alleles plus support for amplification of the same regions across populations and species were additional benefits of these markers. These findings point to the general value of EPIC markers in mosquito population studies.

Rapid diversification and secondary sympatry in Australo-Pacific kingfishers (Aves: Alcedinidae: Todiramphus)

Todiramphus chloris is the most widely distributed of the Pacific's 'great speciators'. Its 50 subspecies constitute a species complex that is distributed over 16 000 km from the Red Sea to Polynesia. We present, to our knowledge, the first comprehensive molecular phylogeny of this enigmatic radiation of kingfishers. Ten Pacific Todiramphus species are embedded within the T. chloris complex, rendering it paraphyletic. Among these is a radiation of five species from the remote islands of Eastern Polynesian, as well as the widespread migratory taxon, Todiramphus sanctus. Our results offer strong support that Pacific Todiramphus, including T. chloris, underwent an extensive range expansion and diversification less than 1 Ma. Multiple instances of secondary sympatry have accumulated in this group, despite its recent origin, including on Australia and oceanic islands in Palau, Vanuatu and the Solomon Islands. Significant ecomorphological and behavioural differences exist between secondarily sympatric lineages, which suggest that pre-mating isolating mechanisms were achieved rapidly during diversification. We found evidence for complex biogeographic patterns, including a novel phylogeographic break in the eastern Solomon Islands that separates a Northern Melanesian clade from Polynesian taxa. In light of our results, we discuss systematic relationships of Todiramphus and propose an updated taxonomy. This paper contributes to our understanding of avian diversification and assembly on islands, and to the systematics of a classically polytypic species complex.

Differential introgression and effective size of marker type influence phylogenetic inference of a recently divergent avian group (Phasianidae: Tympanuchus)

Life history strategies can influence the effective population size (Ne) of loci differently based on their mode of inheritance. Recognizing how this may affect the rate of lineage sorting among marker types is important for studies focused on resolving phylogenetic relationships among recently divergent taxa. In this study, we use gene tree, coalescent-based species tree, and isolation-with-migration analyses to explore the differences between marker types (autosomal, Z-linked, and mitochondrial) in resolving phylogenetic relationships among North American prairie grouse (Tympanuchus). We found that Z-linked loci were more likely to identify monophyletic relationships among prairie grouse species compared to autosomal and mtDNA loci in both species and gene tree analyses, with species tree analyses outperforming gene trees. These results were further supported with isolation-with-migration analyses, where Z-linked loci largely followed a strict isolation model while autosomal loci were more likely to fit a model with gene flow between species following population divergence. While accounting for differences in inheritance pattern (or Ne) for marker type, results suggest that additional factors, such as strong sexual selection and sex-biased introgression (i.e., male-biased postzygotic hybrid behavioral isolation or "unsexy son"), may further explain the decreased diversity levels and increased rate of lineage sorting observed with the Z-linked loci relative to autosomal and mtDNA loci. In fact, to our knowledge no hybrid male prairie grouse have been observed breeding in the wild, yet hybrid females along with backcross females are known to produce viable offspring. Overall, this study highlights that more work is needed to determine how complex models of gene flow (i.e., sex biased introgression) and differences in the effective size among marker types based on differing life history strategies influence divergence date estimation and species delimitation.

The future of ancient DNA: Technical advances and conceptual shifts

Technological innovations such as next generation sequencing and DNA hybridisation enrichment have resulted in multi-fold increases in both the quantity of ancient DNA sequence data and the time depth for DNA retrieval. To date, over 30 ancient genomes have been sequenced, moving from 0.7× coverage (mammoth) in 2008 to more than 50× coverage (Neanderthal) in 2014. Studies of rapid evolutionary changes, such as the evolution and spread of pathogens and the genetic responses of hosts, or the genetics of domestication and climatic adaptation, are developing swiftly and the importance of palaeogenomics for investigating evolutionary processes during the last million years is likely to increase considerably. However, these new datasets require new methods of data processing and analysis, as well as conceptual changes in interpreting the results. In this review we highlight important areas of future technical and conceptual progress and discuss research topics in the rapidly growing field of palaeogenomics.

Genetic and phenotypic characterization of a hybrid zone between polyandrous Northern and Wattled Jacanas in Western Panama

Background

Hybridization provides a unique perspective into the ecological, genetic and behavioral context of speciation. Hybridization is common in birds, but has not yet been reported among bird species with a simultaneously polyandrous mating system; a mating system where a single female defends a harem of males who provide nearly all parental care. Unlike simple polyandry, polyandrous mating is extremely rare in birds, with only 1% of bird species employing this mating system. Although it is classically held that females are “choosy” in avian hybrid systems, nearly-exclusive male parental care raises the possibility that female selection against heterospecific matings might be reduced compared to birds with other mating systems.

Results

We describe a narrow hybrid zone in southwestern Panama between two polyandrous freshwater waders: Northern Jacana, Jacana spinosa and Wattled Jacana, J. jacana. We document coincident cline centers for three phenotypic traits, mtDNA, and one of two autosomal introns. Cline widths for these six markers varied from seven to 142 km, with mtDNA being the narrowest, and five of the six markers having widths less than 100 km. Cline tails were asymmetrical, with greater introgression of J. jacana traits extending westward into the range of J. spinosa. Likewise, within the hybrid zone, the average hybrid index of phenotypic hybrids was significantly biased towards J. spinosa. Species distribution models indicate that the hybrid zone is located at the edge of a roughly 100 km wide overlap where habitat is predicted to be suitable for both species, with more westerly areas suitable only for spinosa and eastward habitats suitable only for J. jacana.

Conclusion

The two species of New World jacanas maintain a narrow, and persistent hybrid zone in western Panama. The hybrid zone may be maintained by the behavioral dominance of J. spinosa counterbalanced by unsuitable habitat for J. spinosa east of the contact zone. Although the two parental species are relatively young, mitochondrial cline width was extremely narrow. This result suggests strong selection against maternally-inherited markers, which may indicate either mitonuclear incompatibilities and/or female choice against heterospecific matings typical of avian hybrid systems, despite jacana sex role reversal.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0227-7) contains supplementary material, which is available to authorized users.

Exonic versus intronic SNPs: contrasting roles in revealing the population genetic differentiation of a widespread bird species

Recent years have seen considerable progress in applying single nucleotide polymorphisms (SNPs) to population genetics studies. However, relatively few have attempted to use them to study the genetic differentiation of wild bird populations and none have examined possible differences of exonic and intronic SNPs in these studies. Here, using 144 SNPs, we examined population genetic differentiation in the saker falcon (Falco cherrug) across Eurasia. The position of each SNP was verified using the recently sequenced saker genome with 108 SNPs positioned within the introns of 10 fragments and 36 SNPs in the exons of six genes, comprising MHC, MC1R and four others. In contrast to intronic SNPs, both Bayesian clustering and principal component analyses using exonic SNPs consistently revealed two genetic clusters, within which the least admixed individuals were found in Europe/central Asia and Qinghai (China), respectively. Pairwise D analysis for exonic SNPs showed that the two populations were significantly differentiated and between the two clusters the frequencies of five SNP markers were inferred to be influenced by selection. Central Eurasian populations clustered in as intermediate between the two main groups, consistent with their geographic position. But the westernmost populations of central Europe showed evidence of demographic isolation. Our work highlights the importance of functional exonic SNPs for studying population genetic pattern in a widespread avian species.

Genetic, ecological and morphological divergence between populations of the endangered Mexican Sheartail hummingbird (Doricha eliza)

The Mexican Sheartail (Doricha eliza), an endangered hummingbird, is endemic to Mexico where two populations have a disjunct distribution. One population is distributed along the northern tip of the Yucatan Peninsula whereas the other is mostly restricted to central Veracruz. Despite their disjunct distribution, previous work has failed to detect morphological or behavioral differences between these populations. Here we use variation in morphology, mtDNA and nuDNA sequences to determine the degree of morphological and molecular divergence between populations, their divergence time, and historical demography. We use species distribution modeling and niche divergence tests to infer the relative roles of vicariance and dispersal in driving divergence in the genus. Our Bayesian and maximum likelihood phylogenetic analyses revealed that Doricha eliza populations form a monophyletic clade and support their sister relationship with D. enicura. We found marked genetic differentiation, with reciprocal monophyly of haplotypes and highly restricted gene flow, supporting a history of isolation over the last 120,000 years. Genetic divergence between populations is consistent with the lack of overlap in environmental space and slight morphological differences between males. Our findings indicate that the divergence of the Veracruz and Yucatan populations is best explained by a combination of a short period of isolation exacerbated by subsequent divergence in climate conditions, and that rather than vicariance, the two isolated ranges of D. eliza are the product of recent colonization and divergence in isolation.

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

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

Development of rapidly evolving intron markers to estimate multilocus species trees of rodents

One of the major challenges in the analysis of closely related species, speciation and phylogeography is the identification of variable sequence markers that allow the determination of genealogical relationships in multiple genomic regions using coalescent and species tree approaches. Rodent species represent nearly half of the mammalian diversity, but so far no systematic study has been carried out to detect suitable informative markers for this group. Here, we used a bioinformatic pipeline to extract intron sequences from rodent genomes available in databases and applied a series of filters that allowed the identification of 208 introns that adequately fulfilled several criteria for these studies. The main required characteristics of the introns were that they had the maximum possible mutation rates, that they were part of single-copy genes, that they had an appropriate sequence length for amplification, and that they were flanked by exons with suitable regions for primer design. In addition, in order to determine the validity of this approach, we chose ten of these introns for primer design and tested them in a panel of eleven rodent species belonging to different representative families. We show that all these introns can be amplified in the majority of species and that, overall, 79% of the amplifications worked with minimum optimization of the annealing temperature. In addition, we confirmed for a pair of sister species the relatively high level of sequence divergence of these introns. Therefore, we provide here a set of adequate intron markers that can be applied to different species of Rodentia for their use in studies that require significant sequence variability.

The advantages of going large: genome-wide SNPs clarify the complex population history and systematics of the threatened western pond turtle

As the field of phylogeography has matured, it has become clear that analyses of one or a few genes may reveal more about the history of those genes than the populations and species that are the targets of study. To alleviate these concerns, the discipline has moved towards larger analyses of more individuals and more genes, although little attention has been paid to the qualitative or quantitative gains that such increases in scale and scope may yield. Here, we increase the number of individuals and markers by an order of magnitude over previously published work to comprehensively assess the phylogeographical history of a well-studied declining species, the western pond turtle (Emys marmorata). We present a new analysis of 89 independent nuclear SNP markers and one mitochondrial gene sequence scored for rangewide sampling of >900 individuals, and compare these to smaller-scale, rangewide genetic and morphological analyses. Our enlarged SNP data fundamentally revise our understanding of evolutionary history for this lineage. Our results indicate that the gains from greatly increasing both the number of markers and individuals are substantial and worth the effort, particularly for species of high conservation concern such as the pond turtle, where accurate assessments of population history are a prerequisite for effective management.

Distinguishing the effects of selection from demographic history in the genetic variation of two sister passerines based on mitochondrial-nuclear comparison

Determining the mechanisms responsible for the distribution of genetic diversity in natural populations has occupied a central role in molecular evolution. Our study was motivated by the unprecedented observation that a widespread Eurasian flycatcher, Ficedula albicilla, exhibited no variation at the mitochondrial DNA (mtDNA) ND2 gene in 75 individuals sampled over a 5000-km distance. In contrast, its sister species, F. parva, had low but considerably higher levels of mtDNA variation. We assessed whether natural selection or demographic factors could explain the absence of mtDNA variation in F. albicilla. Eighteen nuclear genes were sequenced to estimate the two species' phylogeographic histories, and for comparison to the mtDNA data. Multilocus coalescence analyses suggested that F. albicilla experienced a population expansion perhaps following a population bottleneck. Simulations based on this demographic history, however, did not replicate the extremely low level of mtDNA variation. Historical range changes based on ecological niche models also failed to explain the observed mtDNA patterns. Neutrality tests (DHEW and ML-HKA) suggested a non-neutral pattern in the mtDNA of F. albicilla. We found a transmembrane-skewed distribution of nonsynonymous substitutions between the two species, three of which caused functional change; the results implied that positive selection could have targeted mtDNA. Several lines of evidence support selection rather than demographic history as the main force influencing the patterns of mtDNA variation. Despite the influence of natural selection, many of the phylogeographic inferences derived from mtDNA were robust, including species limits and a high level of gene flow among populations within species.

An empirical comparison of character-based and coalescent-based approaches to species delimitation in a young avian complex

The process of discovering species is a fundamental responsibility of systematics. Recently, there has been a growing interest in coalescent-based methods of species delimitation aimed at objectively identifying species early in the divergence process. However, few empirical studies have compared these new methods with character-based approaches for discovering species. In this study, we applied both a character-based and a coalescent-based approaches to delimit species in a closely related avian complex, the light-vented/Taiwan bulbul (Pycnonotus sinensis/Pycnonotus taivanus). Population aggregation analyses of plumage, mitochondrial and 13 nuclear intron character data sets produced conflicting species hypotheses with plumage data suggesting three species, mitochondrial data suggesting two species, and nuclear intron data suggesting one species. Such conflict is expected among recently diverged species, and by integrating all sources of data, we delimited three species verified with independently congruent character evidence as well as a more weakly supported fourth species identified by a single character. Attempts to validate species hypothesis using Bayesian Phylogenetics and Phylogeography (BPP), a coalescent-based method of species delimitation, revealed several issues that can seemingly affect statistical support for species recognition. We found that θ priors had a dramatic impact on speciation probabilities, with lower values consistently favouring splitting and higher values consistently favouring lumping. More resolved guide trees also resulted in overall higher speciation probabilities. Finally, we found suggestive evidence that BPP is sensitive to the divergent effects of nonrandom mating caused by intraspecific processes such as isolation-with-distance, and therefore, BPP may not be a conservative method for delimiting independently evolving population lineages. Based on these concerns, we questioned the reliability of BPP results and based our conclusions about species limits exclusively on character data.

A multi-locus inference of the evolutionary diversification of extant flamingos (Phoenicopteridae)

Background

Modern flamingos (Phoenicopteridae) occupy a highly specialized ecology unique among birds and represent a potentially powerful model system for informing the mechanisms by which a lineage of birds adapts and radiates. However, despite a rich fossil record and well-studied feeding morphology, molecular investigations of the evolutionary progression among modern flamingos have been limited. Here, using three mitochondrial (mtDNA) markers, we present the first DNA sequence-based study of population genetic variation in the widely distributed Chilean Flamingo and, using two mtDNA and 10 nuclear (nDNA) markers, recover the species tree and divergence time estimates for the six extant species of flamingos. Phylogenetic analyses include likelihood and Bayesian frameworks and account for potential gene tree discordance. Analyses of divergence times are fossil calibrated at the divergence of Mirandornithes (flamingos + grebes) and the divergence of crown grebes.

Results

mtDNA sequences confirmed the presence of a single metapopulation represented by two minimally varying mtDNA barcodes in Chilean flamingos. Likelihood and Bayesian methods recovered identical phylogenies with flamingos falling into shallow-keeled (comprising the Greater, American and Chilean Flamingos) and deep-keeled (comprising the Lesser, Andean and James’s Flamingos) sub-clades. The initial divergence among flamingos occurred at or shortly after the Mio-Pliocene boundary (6–3 Ma) followed by quick consecutive divergences throughout the Plio-Pleistocene. There is significant incongruence between the ages recovered by the mtDNA and nDNA datasets, likely due to mutational saturation occurring in the mtDNA loci.

Conclusion

The finding of a single metapopulation in the widespread Chilean Flamingo confirms similar findings in other widespread flamingo species. The robust species phylogeny is congruent with previous classifications of flamingos based on feeding morphology. Modern phoenicopterids likely originated in the New World with each sub-clade dispersing across the Atlantic at least once. Our divergence time estimates place flamingos among the youngest families of birds, counter to the classical notion of flamingos as among the oldest based on biogeography and the fossil record. Finally, we designate ‘Phoeniconaias’ as a junior synonym of ‘Phoenicoparrus’ and redefine the latter genus as containing all flamingos more closely related to Phoenicoparrus andinus than Phoenicopterus roseus.

The ecological and geographic context of morphological and genetic divergence in an understorey-dwelling bird

Advances in understanding the process of species formation require an integrated perspective that includes the evaluation of spatial, ecological and genetic components. One approach is to focus on multiple stages of divergence within the same species. Species that comprise phenotypically different populations segregated in apparently distinct habitats, in which range is presently continuous but was putatively geographically isolated provide an interesting system to study the mechanisms of population divergence. Here, we attempt to elucidate the role of ecology and geography in explaining observed morphological and genetic variation in an understorey-dwelling bird endemic to southeastern Africa, where two subspecies are recognized according to phenotype and habitat affinity. We carried out a range-wide analysis of climatic requirements, morphological and genetic variation across southeast Africa to test the hypothesis that the extent of gene flow among populations of the brown scrub-robin are influenced by their distinct climatic niches. We recovered two distinct trends depending on whether our analyses were hierarchically structured at the subspecies or at the within subspecies level. Between subspecies we found pronounced morphological differentiation associated with strong reproductive isolation (no gene flow) between populations occupying divergent climatic niches characterized by changes in the temperature of the warmest and wettest month. In contrast, within subspecies, we recovered continuous morphological variation with extensive gene flow among populations inhabiting the temperate and sub-tropical forests of southern Africa, despite divergence along the climate axis that is mainly determined by minimum temperature and precipitation of the coldest months. Our results highlight the role of niche divergence as a diversifying force that can promote reproductive isolation in vertebrates.

Coalescent analyses show isolation without migration in two closely related tropical orioles: the case of Icterus graduacauda and Icterus chrysater

The Isthmus of Tehuantepec has played an important role in shaping the avian diversity of Mexico, as well as the rest of the Western Hemisphere. It has been both a barrier and a land connector between North and South America for many groups of birds. Furthermore, climatic change over the Pleistocene has resulted in ecological fluctuations that led to periods of connection and isolation of the highlands in this area. Here we studied the divergence of two species of orioles whose distribution in the highlands is separated by the lowlands of the Isthmus of Tehuantepec: Icterus graduacauda (west of the Isthmus) and Icterus chrysater (east of the Isthmus). We sequenced multiple loci (one mitochondrial gene and six nuclear introns) and performed coalescent analyses (Isolation with Migration) to test whether their divergence resulted from prior occupancy of the ancestral area followed by a vicariant event or recent dispersal from one side or the other of this Isthmus. Results strongly indicate a vicariant event roughly 300,000 years ago in the Pleistocene followed by little or no gene flow. Both mitochondrial and nuclear genes show that the Isthmus of Tehuantepec is a strong barrier to gene flow. Thus, these two species appear to not exchange genes despite their recent divergence and the close geographic proximity of their ranges.