Response to Bakker et al.

Robinson and colleagues respond to the points raised about their paper by Bakker et al.

Genome-wide diversity in the California condor tracks its prehistoric abundance and decline.

Due to their small population sizes, threatened and endangered species frequently suffer from a lack of genetic diversity, potentially leading to inbreeding depression and reduced adaptability.1 During the latter half of the twentieth century, North America's largest soaring bird,2 the California condor (Gymnogyps californianus; Critically Endangered3), briefly went extinct in the wild. Though condors once ranged throughout North America, by 1982 only 22 individuals remained. Following decades of captive breeding and release efforts, there are now >300 free-flying wild condors and ∼200 in captivity. The condor's recent near-extinction from lead poisoning, poaching, and loss of habitat is well documented,4 but much about its history remains obscure. To fill this gap and aid future management of the species, we produced a high-quality chromosome-length genome assembly for the California condor and analyzed its genome-wide diversity. For comparison, we also examined the genomes of two close relatives: the Andean condor (Vultur gryphus; Vulnerable3) and the turkey vulture (Cathartes aura; Least Concern3). The genomes of all three species show evidence of historic population declines. Interestingly, the California condor genome retains a high degree of variation, which our analyses reveal is a legacy of its historically high abundance. Correlations between genome-wide diversity and recombination rate further suggest a history of purifying selection against linked deleterious alleles, boding well for future restoration. We show how both long-term evolutionary forces and recent inbreeding have shaped the genome of the California condor, and provide crucial genomic resources to enable future research and conservation.

Raptor genomes reveal evolutionary signatures of predatory and nocturnal lifestyles.

BACKGROUND: Birds of prey (raptors) are dominant apex predators in terrestrial communities, with hawks (Accipitriformes) and falcons (Falconiformes) hunting by day and owls (Strigiformes) hunting by night. RESULTS: Here, we report new genomes and transcriptomes for 20 species of birds, including 16 species of birds of prey, and high-quality reference genomes for the Eurasian eagle-owl (Bubo bubo), oriental scops owl (Otus sunia), eastern buzzard (Buteo japonicus), and common kestrel (Falco tinnunculus). Our extensive genomic analysis and comparisons with non-raptor genomes identify common molecular signatures that underpin anatomical structure and sensory, muscle, circulatory, and respiratory systems related to a predatory lifestyle. Compared with diurnal birds, owls exhibit striking adaptations to the nocturnal environment, including functional trade-offs in the sensory systems, such as loss of color vision genes and selection for enhancement of nocturnal vision and other sensory systems that are convergent with other nocturnal avian orders. Additionally, we find that a suite of genes associated with vision and circadian rhythm are differentially expressed in blood tissue between nocturnal and diurnal raptors, possibly indicating adaptive expression change during the transition to nocturnality. CONCLUSIONS: Overall, raptor genomes show genomic signatures associated with the origin and maintenance of several specialized physiological and morphological features essential to be apex predators.

Northern Spotted Owl (Strix occidentalis caurina) Genome: Divergence with the Barred Owl (Strix varia) and Characterization of Light-Associated Genes.

We report here the assembly of a northern spotted owl (Strix occidentalis caurina) genome. We generated Illumina paired-end sequence data at 90× coverage using nine libraries with insert lengths ranging from ∼250 to 9,600 nt and read lengths from 100 to 375 nt. The genome assembly is comprised of 8,108 scaffolds totaling 1.26 × 109 nt in length with an N50 length of 3.98 × 106 nt. We calculated the genome-wide fixation index (FST) of S. o. caurina with the closely related barred owl (Strix varia) as 0.819. We examined 19 genes that encode proteins with light-dependent functions in our genome assembly as well as in that of the barn owl (Tyto alba). We present genomic evidence for loss of three of these in S. o. caurina and four in T. alba. We suggest that most light-associated gene functions have been maintained in owls and their loss has not proceeded to the same extent as in other dim-light-adapted vertebrates.

Multi-locus phylogenetic inference among New World Vultures (Aves: Cathartidae).

New World Vultures are large-bodied carrion feeding birds in the family Cathartidae, currently consisting of seven species from five genera with geographic distributions in North and South America. No study to date has included all cathartid species in a single phylogenetic analysis. In this study, we investigated the phylogenetic relationships among all cathartid species using five nuclear (nuc; 4060bp) and two mitochondrial (mt; 2165bp) DNA loci with fossil calibrated gene tree (27 outgroup taxa) and coalescent-based species tree (2 outgroup taxa) analyses. We also included an additional four nuclear loci (2578bp) for the species tree analysis to explore changes in nodal support values. Although the stem lineage is inferred to have originated ∼69 million years ago (Ma; 74.5-64.9 credible interval), a more recent basal split within Cathartidae was recovered at ∼14Ma (17.1-11.1 credible interval). Two primary clades were identified: (1) Black Vulture (Coragyps atratus) together with the three Cathartes species (Lesser C. burrovianus and Greater C. melambrotus Yellow-headed Vultures, and Turkey Vulture C. aura), and (2) King Vulture (Sarcoramphus papa), California (Gymnogyps californianus) and Andean (Vultur gryphus) Condors. Support for taxon relationships within the two basal clades were inconsistent between analyses with the exception of Black Vulture sister to a monophyletic Cathartes clade. Increased support for a yellow-headed vulture clade was recovered in the species tree analysis using the four additional nuclear loci. Overall, these results are in agreement with cathartid life history (e.g. olfaction ability and behavior) and contrasting habitat affinities among sister taxa with overlapping geographic distributions. More research is needed using additional molecular loci to further resolve the phylogenetic relationships within the two basal cathartid clades, as speciation appeared to have occurred in a relatively short period of time.

Isolation of a Complete Circular Virus Genome Sequence from an Alaskan Black-Capped Chickadee (Poecile atricapillus) Gastrointestinal Tract Sample.

We report here the genome sequence of a circular virus isolated from samples of an Alaskan black-capped chickadee (Poecile atricapillus) gastrointestinal tract. The genome is 2,152 bp in length and is most similar (30 to 44.5% amino acid identity) to the genome sequences of other single-stranded DNA (ssDNA) circular viruses belonging to the gemycircularvirus group.

Response to Comment on "Whole-genome analyses resolve early branches in the tree of life of modern birds".

Mitchell et al. argue that divergence-time estimates for our avian phylogeny were too young because of an "inappropriate" maximum age constraint for the most recent common ancestor of modern birds and that, as a result, most modern bird orders diverged before the Cretaceous-Paleogene mass extinction event 66 million years ago instead of after. However, their interpretations of the fossil record and timetrees are incorrect.

Phylogenomic analyses data of the avian phylogenomics project.

BACKGROUND: Determining the evolutionary relationships among the major lineages of extant birds has been one of the biggest challenges in systematic biology. To address this challenge, we assembled or collected the genomes of 48 avian species spanning most orders of birds, including all Neognathae and two of the five Palaeognathae orders. We used these genomes to construct a genome-scale avian phylogenetic tree and perform comparative genomic analyses. FINDINGS: Here we present the datasets associated with the phylogenomic analyses, which include sequence alignment files consisting of nucleotides, amino acids, indels, and transposable elements, as well as tree files containing gene trees and species trees. Inferring an accurate phylogeny required generating: 1) A well annotated data set across species based on genome synteny; 2) Alignments with unaligned or incorrectly overaligned sequences filtered out; and 3) Diverse data sets, including genes and their inferred trees, indels, and transposable elements. Our total evidence nucleotide tree (TENT) data set (consisting of exons, introns, and UCEs) gave what we consider our most reliable species tree when using the concatenation-based ExaML algorithm or when using statistical binning with the coalescence-based MP-EST algorithm (which we refer to as MP-EST*). Other data sets, such as the coding sequence of some exons, revealed other properties of genome evolution, namely convergence. CONCLUSIONS: The Avian Phylogenomics Project is the largest vertebrate phylogenomics project to date that we are aware of. The sequence, alignment, and tree data are expected to accelerate analyses in phylogenomics and other related areas.

Rapid diversification of falcons (Aves: Falconidae) due to expansion of open habitats in the Late Miocene.

Understanding how and why lineages diversify is central to understanding the origins of biological diversity. The avian family Falconidae (caracaras, forest-falcons, falcons) has an uneven distribution of species among multiple well-supported clades, and provides a useful system for testing hypotheses about diversification rate and correlation with environmental changes. We analyzed eight independent loci for 1-7 individuals from each of the 64 currently recognized Falconidae species, together with two fossil falconid temporal calibrations, to assess phylogeny, absolute divergence times and potential shifts in diversification rate. Our analyses supported similar diversification ages in the Early to Middle Miocene for the three traditional subfamilies, Herpetotherinae, Polyborinae and Falconinae. We estimated that divergences within the subfamily Falconinae began about 16mya and divergences within the most species-rich genus, Falco, including about 60% of all Falconidae species, began about 7.5mya. We found evidence for a significant increase in diversification rate at the basal phylogenetic node for the genus Falco, and the timing for this rate shift correlates generally with expansion of C4 grasslands beginning around the Miocene/Pliocene transition. Concomitantly, Falco lineages that are distributed primarily in grassland or savannah habitats, as opposed to woodlands, and exhibit migratory, as opposed to sedentary, behavior experienced a higher diversification rate.

Approaching a state shift in Earth's biosphere.

Localized ecological systems are known to shift abruptly and irreversibly from one state to another when they are forced across critical thresholds. Here we review evidence that the global ecosystem as a whole can react in the same way and is approaching a planetary-scale critical transition as a result of human influence. The plausibility of a planetary-scale 'tipping point' highlights the need to improve biological forecasting by detecting early warning signs of critical transitions on global as well as local scales, and by detecting feedbacks that promote such transitions. It is also necessary to address root causes of how humans are forcing biological changes.

Pliocene diversification within the South American Forest falcons (Falconidae: Micrastur).

The Neotropics are one of the most species rich regions on Earth, with over 3150 species of birds. This unrivaled biodiversity has been attributed to higher proportions of mountain ranges, tropical rain forest or rain fall in the forest than in any other major biogeographic regions. Five primary hypotheses aim to explain processes of diversification within the Neotropics; (1) the Pleistocene refuge hypothesis, (2) the riverine barrier hypothesis, (3) the Miocene marine incursions hypothesis, (4) the ecological gradient hypothesis, and (5) the impact of the last Andean uplift serving as a barrier between eastern and western population Andean populations. We assessed these hypotheses to see which best explained the species richness of the forest-falcons (Micrastur), a poorly known lineage of birds that inhabit lowland and mid-elevation humid forest. Our analyses suggest all speciation events within the genus Micrastur probably occurred in the last 2.5-3.6 myrs, at or before the Pliocene/Pleistocene boundary, with the basal split within the genus being 7 myrs old. Hence our data allow us to formerly reject the classical Pleistocene refuge for Micrastur, Our divergence time estimates are younger that dates for the Miocene marine incursions, the riverine barrier and the Andean uplift hypotheses.

The use of genetics for the management of a recovering population: temporal assessment of migratory peregrine falcons in North America.

BACKGROUND: Our ability to monitor populations or species that were once threatened or endangered and in the process of recovery is enhanced by using genetic methods to assess overall population stability and size over time. This can be accomplished most directly by obtaining genetic measures from temporally-spaced samples that reflect the overall stability of the population as given by changes in genetic diversity levels (allelic richness and heterozygosity), degree of population differentiation (F(ST) and D(EST)), and effective population size (N(e)). The primary goal of any recovery effort is to produce a long-term self-sustaining population, and these genetic measures provide a metric by which we can gauge our progress and help make important management decisions. METHODOLOGY/PRINCIPAL FINDINGS: The peregrine falcon in North America (Falco peregrinus tundrius and anatum) was delisted in 1994 and 1999, respectively, and its abundance will be monitored by the species Recovery Team every three years until 2015. Although the United States Fish and Wildlife Service makes a distinction between tundrius and anatum subspecies, our genetic results based on eleven microsatellite loci suggest limited differentiation that can be attributed to an isolation by distance relationship and warrant no delineation of these two subspecies in its northern latitudinal distribution from Alaska through Canada into Greenland. Using temporal samples collected at Padre Island, Texas during migration (seven temporal time periods between 1985-2007), no significant differences in genetic diversity or significant population differentiation in allele frequencies between time periods were observed and were indistinguishable from those obtained from tundrius/anatum breeding locations throughout their northern distribution. Estimates of harmonic mean N(e) were variable and imprecise, but always greater than 500 when employing multiple temporal genetic methods. CONCLUSIONS/SIGNIFICANCE: These results, including those from simulations to assess the power of each method to estimate N(e), suggest a stable or growing population, which is consistent with ongoing field-based monitoring surveys. Therefore, historic and continuing efforts to prevent the extinction of the peregrine falcon in North America appear successful with no indication of recent decline, at least from the northern latitude range-wide perspective. The results also further highlight the importance of archiving samples and their use for continual assessment of population recovery and long-term viability.

Population structure and plumage polymorphism: The intraspecific evolutionary relationships of a polymorphic raptor, Buteo jamaicensis harlani.

BACKGROUND: Phenotypic and molecular genetic data often provide conflicting patterns of intraspecific relationships confounding phylogenetic inference, particularly among birds where a variety of environmental factors may influence plumage characters. Among diurnal raptors, the taxonomic relationship of Buteo jamaicensis harlani to other B. jamaicensis subspecies has been long debated because of the polytypic nature of the plumage characteristics used in subspecies or species designations. RESULTS: To address the evolutionary relationships within this group, we used data from 17 nuclear microsatellite loci, 430 base pairs of the mitochondrial control region, and 829 base pairs of the melanocortin 1 receptor (Mc1r) to investigate molecular genetic differentiation among three B. jamaicensis subspecies (B. j. borealis, B. j. calurus, B. j. harlani). Bayesian clustering analyses of nuclear microsatellite loci showed no significant differences between B. j. harlani and B. j. borealis. Differences observed between B. j. harlani and B. j. borealis in mitochondrial and microsatellite data were equivalent to those found between morphologically similar subspecies, B. j. borealis and B. j. calurus, and estimates of migration rates among all three subspecies were high. No consistent differences were observed in Mc1r data between B. j. harlani and other B. jamaicensis subspecies or between light and dark color morphs within B. j. calurus, suggesting that Mc1r does not play a significant role in B. jamaicensis melanism. CONCLUSIONS: These data suggest recent interbreeding and gene flow between B. j. harlani and the other B. jamaicensis subspecies examined, providing no support for the historical designation of B. j. harlani as a distinct species.

It's not too late for the harpy eagle (Harpia harpyja): high levels of genetic diversity and differentiation can fuel conservation programs.

BACKGROUND: The harpy eagle (Harpia harpyja) is the largest Neotropical bird of prey and is threatened by human persecution and habitat loss and fragmentation. Current conservation strategies include local education, captive rearing and reintroduction, and protection or creation of trans-national habitat blocks and corridors. Baseline genetic data prior to reintroduction of captive-bred stock is essential for guiding such efforts but has not been gathered previously. METHODOLOGY/FINDINGS: We assessed levels of genetic diversity, population structure and demographic history for harpy eagles using samples collected throughout a large portion of their geographic distribution in Central America (n = 32) and South America (n = 31). Based on 417 bp of mitochondrial control region sequence data, relatively high levels of haplotype and nucleotide diversity were estimated for both Central and South America, although haplotype diversity was significantly higher for South America. Historical restriction of gene flow across the Andes (i.e. between our Central and South American subgroups) is supported by coalescent analyses, the haplotype network and significant F(ST) values, however reciprocally monophyletic lineages do not correspond to geographical locations in maximum likelihood analyses. A sudden population expansion for South America is indicated by a mismatch distribution analysis, and further supported by significant (p<0.05) negative values of Fu and Li's D(F) and F, and Fu's F(S). This expansion, estimated at approximately 60 000 years BP (99 000-36 000 years BP 95% CI), encompasses a transition from a warm and dry time period prior to 50 000 years BP to an interval of maximum precipitation (50 000-36 000 years BP). Notably, this time period precedes the climatic and habitat changes associated with the last glacial maximum. In contrast, a multimodal distribution of haplotypes was observed for Central America suggesting either population equilibrium or a recent decline. SIGNIFICANCE: High levels of mitochondrial genetic diversity in combination with genetic differentiation among subgroups within regions and between regions highlight the importance of local population conservation in order to preserve maximal levels of genetic diversity in this species. Evidence of historically restricted female-mediated gene flow is an important consideration for captive-breeding programs.

Long-term survival despite low genetic diversity in the critically endangered Madagascar fish-eagle.

The critically endangered Madagascar fish-eagle (Haliaeetus vociferoides) is considered to be one of the rarest birds of prey globally and at significant risk of extinction. In the most recent census, only 222 adult individuals were recorded with an estimated total breeding population of no more than 100-120 pairs. Here, levels of Madagascar fish-eagle population genetic diversity based on 47 microsatellite loci were compared with its sister species, the African fish-eagle (Haliaeetus vocifer), and 16 of these loci were also characterized in the white-tailed eagle (Haliaeetus albicilla) and the bald eagle (Haliaeetus leucocephalus). Overall, extremely low genetic diversity was observed in the Madagascar fish-eagle compared to other surveyed Haliaeetus species. Determining whether this low diversity is the result of a recent bottleneck or a more historic event has important implications for their conservation. Using a Bayesian coalescent-based method, we show that Madagascar fish-eagles have maintained a small effective population size for hundreds to thousands of years and that its low level of neutral genetic diversity is not the result of a recent bottleneck. Therefore, efforts made to prevent Madagascar fish-eagle extinction should place high priority on maintenance of habitat requirements and reducing direct and indirect human persecution. Given the current rate of deforestation in Madagascar, we further recommend that the population be expanded to occupy a larger geographical distribution. This will help the population persist when exposed to stochastic factors (e.g. climate and disease) that may threaten a species consisting of only 200 adult individuals while inhabiting a rapidly changing landscape.

Strong mitochondrial DNA support for a Cretaceous origin of modern avian lineages.

BACKGROUND: Determining an absolute timescale for avian evolutionary history has proven contentious. The two sources of information available, paleontological data and inference from extant molecular genetic sequences (colloquially, 'rocks' and 'clocks'), have appeared irreconcilable; the fossil record supports a Cenozoic origin for most modern lineages, whereas molecular genetic estimates suggest that these same lineages originated deep within the Cretaceous and survived the K-Pg (Cretaceous-Paleogene; formerly Cretaceous-Tertiary or K-T) mass-extinction event. These two sources of data therefore appear to support fundamentally different models of avian evolution. The paradox has been speculated to reflect deficiencies in the fossil record, unrecognized biases in the treatment of genetic data or both. Here we attempt to explore uncertainty and limit bias entering into molecular divergence time estimates through: (i) improved taxon (n = 135) and character (n = 4594 bp mtDNA) sampling; (ii) inclusion of multiple cladistically tested internal fossil calibration points (n = 18); (iii) correction for lineage-specific rate heterogeneity using a variety of methods (n = 5); (iv) accommodation of uncertainty in tree topology; and (v) testing for possible effects of episodic evolution. RESULTS: The various 'relaxed clock' methods all indicate that the major (basal) lineages of modern birds originated deep within the Cretaceous, although temporal intraordinal diversification patterns differ across methods. We find that topological uncertainty had a systematic but minor influence on date estimates for the origins of major clades, and Bayesian analyses assuming fixed topologies deliver similar results to analyses with unconstrained topologies. We also find that, contrary to expectation, rates of substitution are not autocorrelated across the tree in an ancestor-descendent fashion. Finally, we find no signature of episodic molecular evolution related to either speciation events or the K-Pg boundary that could systematically mislead inferences from genetic data. CONCLUSION: The 'rock-clock' gap has been interpreted by some to be a result of the vagaries of molecular genetic divergence time estimates. However, despite measures to explore different forms of uncertainty in several key parameters, we fail to reconcile molecular genetic divergence time estimates with dates taken from the fossil record; instead, we find strong support for an ancient origin of modern bird lineages, with many extant orders and families arising in the mid-Cretaceous, consistent with previous molecular estimates. Although there is ample room for improvement on both sides of the 'rock-clock' divide (e.g. accounting for 'ghost' lineages in the fossil record and developing more realistic models of rate evolution for molecular genetic sequences), the consistent and conspicuous disagreement between these two sources of data more likely reflects a genuine difference between estimated ages of (i) stem-group origins and (ii) crown-group morphological diversifications, respectively. Further progress on this problem will benefit from greater communication between paleontologists and molecular phylogeneticists in accounting for error in avian lineage age estimates.