Search

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 Endangered 3 ), 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; Vulnerable 3 ) and the turkey vulture (Cathartes aura; Least Concern 3 ). 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., Competing Interests: Declaration of interests J.D.W. receives research funding from Sierra Pacific Industries. E.L.A. is Scientific Advisory Board co-chair and consultant for HolyHaid Lab Corporation (Shenzhen, China), whose parent company is Hollyhigh International Capital (Beijing and Shanghai, China)., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

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.

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., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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., (Copyright © 2016 Elsevier Inc. All rights reserved.)

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., (Copyright © 2015, American Association for the Advancement of Science.)

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., (Copyright © 2015 Hanna et al.)

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.

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., (Published by Elsevier Inc.)

To better determine the history of modern birds, we performed a genome-scale phylogenetic analysis of 48 species representing all orders of Neoaves using phylogenomic methods created to handle genome-scale data. We recovered a highly resolved tree that confirms previously controversial sister or close relationships. We identified the first divergence in Neoaves, two groups we named Passerea and Columbea, representing independent lineages of diverse and convergently evolved land and water bird species. Among Passerea, we infer the common ancestor of core landbirds to have been an apex predator and confirm independent gains of vocal learning. Among Columbea, we identify pigeons and flamingoes as belonging to sister clades. Even with whole genomes, some of the earliest branches in Neoaves proved challenging to resolve, which was best explained by massive protein-coding sequence convergence and high levels of incomplete lineage sorting that occurred during a rapid radiation after the Cretaceous-Paleogene mass extinction event about 66 million years ago., (Copyright © 2014, American Association for the Advancement of Science.)

Birds are the most species-rich class of tetrapod vertebrates and have wide relevance across many research fields. We explored bird macroevolution using full genomes from 48 avian species representing all major extant clades. The avian genome is principally characterized by its constrained size, which predominantly arose because of lineage-specific erosion of repetitive elements, large segmental deletions, and gene loss. Avian genomes furthermore show a remarkably high degree of evolutionary stasis at the levels of nucleotide sequence, gene synteny, and chromosomal structure. Despite this pattern of conservation, we detected many non-neutral evolutionary changes in protein-coding genes and noncoding regions. These analyses reveal that pan-avian genomic diversity covaries with adaptations to different lifestyles and convergent evolution of traits., (Copyright © 2014, American Association for the Advancement of Science.)

Competing Interests: The authors have declared that no competing interests exist.

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.

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., (Copyright © 2011 Elsevier Inc. All rights reserved.)

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.

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.

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.

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.

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.

Little is known about the possible influence that past glacial events have had on the phylogeography and population structure of avian predators in the Arctic and sub-Arctic. In this study, we use microsatellite and mitochondrial control region DNA variation to investigate the population genetic structure of gyrfalcons (Falco rusticolus) throughout a large portion of their circumpolar distribution. In most locations sampled, the mtDNA data revealed little geographic structure; however, five out of eight mtDNA haplotypes were unique to a particular geographic area (Greenland, Iceland, or Alaska) and the Iceland population differed from others based on haplotype frequency differences (F(ST)). With the microsatellite results, significant population structure (F(ST), principal components analysis, and cluster analysis) was observed identifying Greenland and Iceland as separate populations, while Norway, Alaska and Canada were identified as a single population consistent with contemporary gene flow across Russia. Within Greenland, differing levels of gene flow between western and eastern sampling locations was indicated with apparent asymmetric dispersal in western Greenland from north to south. This dispersal bias is in agreement with the distribution of plumage colour variants with white gyrfalcons in much higher proportion in northern Greenland. Lastly, because the mtDNA control region sequence differed by only one to four nucleotides from a common haplotype among all gyrfalcons, we infer that the observed microsatellite population genetic structure has developed since the last glacial maximum. This conclusion is further supported by our finding that a closely related species, the saker falcon (Falco cherrug), has greater genetic heterogeneity, including mtDNA haplotypes differing by 1-16 nucleotide substitutions from a common gyrfalcon haplotype. This is consistent with gyrfalcons having expanded rapidly from a single glacial-age refugium to their current circumpolar distribution. Additional sampling of gyrfalcons from Fennoscandia and Russia throughout Siberia is necessary to test putative gene flow between Norway and Alaska and Canada as suggested by this study.

The evolution of penguins has been investigated extensively, although inconclusively, by morphologists, biogeographers and molecular phylogeneticists. We investigated this issue using retroposon analysis of insertions of CR1, which is a member of the LINE (long interspersed element) family, in the genomes of penguins and penguin relatives. The retroposon method is a powerful tool for identifying monophyletic groups. Because retroposons often show different relative frequencies of retroposition during evolution, it is first necessary to identify a certain subgroup that was specifically active during the period when the species in question diverged. Hence, we systematically analyzed many CR1 members isolated from penguin and penguin-related genomes. These CR1s are divided into at least three distinct subgroups that share diagnostic nucleotide insertions and/or deletions, namely, penguin CR1 Sph I, Sph II type A and Sph II type B. The analysis of the inserted retroposons by PCR revealed that different CR1 subfamilies or types had amplified at different rates among different periods during penguin evolution. Namely, the penguin CR1 Sph I subfamily had higher rates of retroposition in a common ancestor of all orders examined in this study or at least in a common ancestor of all extant penguins, and the subfamily Sph II type A also had the same tendency. Therefore, these CR1 members can be used to elucidate the phylogenetic relationships of Sphenisciformes (penguins) among different avian orders. In contrast, the penguin CR1 Sph II type B subfamily had higher rates of retroposition just before and after the emergence of the extant genera in Spheniscidae, suggesting that they are useful for elucidating the intra-relationships among extant penguins. This is the first report for the characterization among the members of CR1 family in avian genomes excluding those of chickens. Hence, this work will be a cornerstone for elucidating the phylogenetic relationships in penguin evolution using the retroposon method.

We assessed phylogenetic relationships for birds of prey in the family Accipitridae using molecular sequence from two mitochondrial genes (1047 bases ND2 and 1041 bases cyt-b) and one nuclear intron (1074 bases beta-fibrinogen intron 7). We sampled representatives of all 14 Accipitridae subfamilies, focusing on four subfamilies of eagles (booted eagles, sea eagles, harpy eagles, and snake eagles) and two subfamilies of Old World vultures (Gypaetinae and Aegypiinae) with nearly all known species represented. Multiple well-supported relationships among accipitrids identified with DNA differ from those traditionally recognized based on morphology or life history traits. Monophyly of sea eagles (Haliaeetinae) and booted eagles (Aquilinae) was supported; however, harpy eagles (Harpiinae), snake eagles (Circaetinae), and Old World vultures were found to be non-monophyletic. The Gymnogene (Polyboroides typus) and the Crane Hawk (Geranospiza caerulescens) were not found to be close relatives, presenting an example of convergent evolution for specialized limb morphology enabling predation on cavity nesting species. Investigation of named subspecies within Hieraaetus fasciatus and H. morphnoides revealed significant genetic differentiation or non-monophyly supporting recognition of H. spilogaster and H. weiskei as distinctive species.

The Cape Verde kite (Milvus milvus fasciicauda) is considered to be one of the rarest birds of prey in the world and at significant risk of extinction. For this reason there is great interest in both the taxonomic and the population status of this group. To help resolve its taxonomic status, we provide phylogenetic analyses based on three mitochondrial genes for a sampling of kites in the genus Milvus, including a broad geographical sampling of black kites (Milvus migrans), red kites (Milvus milvus), Cape Verde kite museum specimens collected between 1897 and 1924, and five kites trapped on the Cape Verde Islands during August 2002. We found that the historical Cape Verde kites, including the type specimen, were non-monophyletic and scattered within a larger red kite clade. The recently trapped kites from the Cape Verde Islands were all phylogenetically diagnosed as black kites. Our findings suggest that the traditional Cape Verde kite is not a distinctive evolutionary unit, and the case for species status, as recently suggested by others, is not supported. We do find support for recognition of at least one clade of yellow-billed kites, traditionally considered as a black kite subspecies, as a distinctive phylogenetic species.

We provide phylogenetic analyses for primary Reptilia lineages including, for the first time, Sphenodon punctatus (tuatara) using data from whole mitochondrial genomes. Our analyses firmly support a sister relationship between Sphenodon and Squamata, which includes lizards and snakes. Using Sphenodon as an outgroup for select squamates, we found evidence indicating a sister relationship, among our study taxa, between Serpentes (represented by Dinodon) and Varanidae. Our analyses support monophyly of Archosauria, and a sister relationship between turtles and archosaurs. This latter relationship is congruent with a growing set of morphological and molecular analyses placing turtles within crown Diapsida and recognizing them as secondarily anapsid (lacking a skull fenestration). Inclusion of Sphenodon, as the only surviving member of Sphenodontia (with fossils from the mid-Triassic), helps to fill a sampling gap within previous analyses of reptilian phylogeny. We also report a unique configuration for the mitochondrial genome of Sphenodon, including two tRNA(Lys) copies and an absence of ND5, tRNA(His), and tRNA(Thr) genes.

We investigated the phylogenetic relationships of species and subspecies of the cosmopolitan genus Pica using 813 bp of the mitochondrial genome (including portions of 16s rDNA, tRNA-Leu, and ND1). The phylogenetic relationships within the genus Pica revealed in our molecular analyses can be summarized as follows: (1). the Korean magpie (Pica pica sericea) appears basal within the genus Pica; (2). the European magpie (Pica pica pica) shows a close relationship to the Kamchatkan magpie (Pica pica camtschatica); (3). two North American species (Pica hudsonia and Pica nuttalli) shows a sister-group relationship; (4). most importantly, the European+Kamchatkan clade appears more closely related to the North American clade than to Korean magpies. Based on these results and genetic distance data, it is possible that members of an ancestral magpie lineage in East Asia initially moved north to form Kamchatkan magpies and then crossed the Bering land bridge to found North American taxa. At a later date, a group might have split off from Kamchatkan magpies and migrated west to form the Eurasian subspecies. The divergence between the two North American taxa appears to have happened no later than the divergence of Eurasian subspecies and both processes appear to have been relatively rapid. Rather than the formation of P. hudsonia by re-colonization from an Asian magpie ancestor, as suggested by, our data suggest a shared ancestry between P. hudsonia and P. nuttalli. Based on the above findings, including phylogenetic placement of P. hudsonia and P. nuttalli as nested within the larger Pica pica clade, and the lack of evidence suggesting reproductive isolation within the genus Pica, we believe that the current classification may be inaccurate. A more conservative classification would recognize one monophyletic species (i.e., P. pica) and treat P. nuttalli and P. hudsonia as subspecies (i.e., P. p. nuttalli and P. p. hudsonia). More extensive studies on the population genetics and biogeography of magpies should be conducted to better inform any taxonomic decisions.

The apparently rapid and ancient diversification of many avian orders complicates the resolution of their relationships using molecular data. Recent studies based on complete mitochondrial DNA (mtDNA) sequences or shorter lengths of nuclear sequence have helped corroborate the basic structure of the avian tree (e.g., a basal split between Paleognathae and Neognathae) but have made relatively little progress in resolving relationships among the many orders within Neoaves. We explored the potential of a moderately sized mtDNA data set ( approximately 5000 bp for each of 41 taxa), supplemented with data from a nuclear intron ( approximately 700 bp per taxon), to resolve relationships among avian orders. Our sampling of taxa addresses two issues: (1). the sister relationship and monophyly, respectively, of Anseriformes and Galliformes and (2). relationships of the enigmatic hoatzin Opisthocomus hoazin. Our analyses support a basal split between Galloanserae and Neoaves within Neognathae and monophyly of both Galliformes and Anseriformes. Within Galliformes, megapodes and then cracids branch basally. Within Anseriformes, mitochondrial data support a screamer (Anhimidae) plus magpie goose (Anseranatidae) clade. This result, however, may be an artifact of divergent base composition in one of the two anatids we sampled. With deletion of the latter taxon, Anseranas is sister to anatids as in traditional arrangements and recent morphological studies. Although our data provide limited resolution of relationships within Neoaves, we find no support for a sister relationship between either cuckoos (Cuculiformes) or turacos (Musophagiformes) and hoatzin. Both mitochondrial and nuclear data are consistent with a relationship between hoatzin and doves (Columbiformes), although this result is weakly supported. We also show that mtDNA sequences reported in another recent study included pervasive errors that biased the analysis towards finding a sister relationship between hoatzin and turacos.

The sudden appearance and potential lethality of severe acute respiratory syndrome associated coronavirus (SARS-CoV) in humans has focused attention on understanding its origins. Here, we assess phylogenetic relationships for the SARS-CoV lineage as well as the history of host-species shifts for SARS-CoV and other coronaviruses. We used a Bayesian phylogenetic inference approach with sliding window analyses of three SARS-CoV proteins: RNA dependent RNA polymerase (RDRP), nucleocapsid (N) and spike (S). Conservation of RDRP allowed us to use a set of Arteriviridae taxa to root the Coronaviridae phylogeny. We found strong evidence for a recombination breakpoint within SARS-CoV RDRP, based on different, well supported trees for a 5' fragment (supporting SARS-CoV as sister to a clade including all other coronaviruses) and a 3' fragment (supporting SARS-CoV as sister to group three avian coronaviruses). These different topologies are statistically significant: the optimal 5' tree could be rejected for the 3' region, and the optimal 3' tree could be rejected for the 5' region. We did not find statistical evidence for recombination in analyses of N and S, as there is little signal to differentiate among alternative trees. Comparison of phylogenetic trees for 11 known host-species and 36 coronaviruses, representing coronavirus groups 1-3 and SARS-CoV, based on N showed statistical incongruence indicating multiple host-species shifts for coronaviruses. Inference of host-species associations is highly sensitive to sampling and must be considered cautiously. However, current sampling suggests host-species shifts between mouse and rat, chicken and turkey, mammals and manx shearwater, and humans and other mammals. The sister relationship between avian coronaviruses and the 3' RDRP fragment of SARS-CoV suggests an additional host-species shift. Demonstration of recombination in the SARS-CoV lineage indicates its potential for rapid unpredictable change, a potentially important challenge for public health management and for drug and vaccine development.

Until recently, none of the diverse elements bearing reverse transcriptase (retroids) have been known from Archaea. However, in the recently published genomes of the acetate-utilizing archaeal methanogens, Methanosarcina acetivorans and M. mazei, several open reading frames (ORFs) are annotated as reverse transcriptase (RT). These annotations led us to the characterization of a retron and 13 retrointrons, including three twintrons, clustered at seven loci of the M. acetivorans genome, and four retrointrons at two loci of the M. mazei genome. Based on a phylogeny of the RT ORFs, we infer four lateral gene transfers (LGT) of these retroids from Bacteria to Archaea and of retrointron mobility within the Archaea genomes. Our phylogenetic analysis also identifies several novel retrons from GenBank in the bacterial groups Firmicutes, Fusobacteria, Cyanobacteria and beta-Proteobacteria, as well as in M. acetivorans. The discovery of retrointrons in Archaea as a consequence of LGT from Bacteria suggests that they did not originate in the progenote and parallels the "mitochondrial seed" theory of the origin of spliceosomes. Extending the known phylogenetic distribution of retroids to Archaea is consistent with the view that they have played a significant role in evolution of genomes throughout the tree of life.

Recent molecular studies addressing the phylogenetic relationships of avian orders have had conflicting results. While studies using nuclear DNA sequences tend to support traditional taxonomic views, also supported by morphological data [(paleognaths (galloanseres (all other birds)))], with songbirds forming a clade within Neoaves (all other birds), analyses with complete mtDNA genomes have resulted in topologies that place songbirds as one of the earliest-diverging avian lineages. Considering that over half of the extant bird species are songbirds, these different results have very different implications for our understanding of avian evolution. We analyzed data sets comprising nearly 4 kb of mitochondrial DNA (mtDNA) (complete 12S, ND1, ND2, and cytochrome b) plus 600 bp of the nuclear gene c-mos for 15 birds that were chosen to represent all major avian clades and to minimize potential long-branch attraction problems; we used a partition-specific maximum likelihood approach. Our results show congruence with respect to the ingroup among phylogenies obtained with mtDNA and the nuclear gene c-mos, separately or combined. The data sets support a traditional avian taxonomy, with paleognaths (ratites and tinamous) occupying a basal position and with songbirds more derived and forming a monophyletic group. We also show that, for mtDNA studies, turtles may be a better outgroup for birds than crocodilians because of their slower rate of sequence evolution.

A gastroenterologist was convicted of attempted second-degree murder by injecting his former girlfriend with blood or blood-products obtained from an HIV type 1 (HIV-1)-infected patient under his care. Phylogenetic analyses of HIV-1 sequences were admitted and used as evidence in this case, representing the first use of phylogenetic analyses in a criminal court case in the United States. Phylogenetic analyses of HIV-1 reverse transcriptase and env DNA sequences isolated from the victim, the patient, and a local population sample of HIV-1-positive individuals showed the victim's HIV-1 sequences to be most closely related to and nested within a lineage comprised of the patient's HIV-1 sequences. This finding of paraphyly for the patient's sequences was consistent with the direction of transmission from the patient to the victim. Analysis of the victim's viral reverse transcriptase sequences revealed genotypes consistent with known mutations that confer resistance to AZT, similar to those genotypes found in the patient. A priori establishment of the patient and victim as a suspected transmission pair provided a clear hypothesis for phylogenetic testing. All phylogenetic models and both genes examined strongly supported the close relationship between the HIV-1 sequences of the patient and the victim. Resampling of blood from the suspected transmission pair and independent sequencing by different laboratories provided precaution against laboratory error.

The avian subfamily Tetraoninae (grouse and ptarmigan) is a Holarctic group in the order Galliformes distinguished by morphological adaptations to cold environments and behavioral traits associated with elaborate courtship. Here we investigate the relationships of 17 tetraonines and 12 other galliform species using mitochondrial 12S and ND2 sequence data. We found support for the recent phylogenetic classification that separates the genus Dendragapus into two genera, Falcipennis and Dendragapus. In addition, we found support for a tetraonine clade in which the first divergence is between Bonasa umbellus and all others, followed by divergence between a Bonasa bonasia/Bonasa sewerzowi clade and the remaining tetraonines. Falcipennis canadensis is sister to a clade with four Tetrao species, and the genus Centrocercus is sister to a Dendragapus obscurus/Tympanuchus clade. Our data indicate a basal position for Cracidae and Megapodiidae among the five recognized galliform families. We also found strong support for the monophyly of Phasianidae, although the relative positions of Numididae and Odontiphoridae remains unresolved. We use a maximum likelihood approach to infer ages of 37mya for divergence of Numididae and Phasianidae and 28mya for the divergence of Tetraoninae and Meleagris gallopavo. These estimates must be viewed as tentative as they depend on tests of rates of molecular evolution and accurate fossil dates.

Retroviral and other reverse transcriptase (RT)-containing sequences may be subject to unique evolutionary pressures, and models of molecular sequence evolution developed using other kinds of sequences may not be optimal. Here we develop and present a new substitution matrix for maximum likelihood (ML) phylogenetic analysis which has been optimized on a dataset of 33 amino acid sequences from the retroviral Pol proteins. When compared to other matrices, this model (rtREV) yields higher log-likelihood values on a range of datasets including lentiviruses, spumaviruses, betaretroviruses, gammaretroviruses, and other elements containing reverse transcriptase. We provide evidence that rtREV is a more realistic evolutionary model for analyses of the pol gene, although it is inapplicable to analyses involving the gag gene.

Systematic studies of Fringillidae have long been problematic due to their apparent recent and explosive diversification. We present phylogenetic hypotheses of 44 fringillids that represent the overall diversity of the family, based on 3.2 kb of mitochondrial DNA sequences, and phylogenetic analyses for a subset of fringillids based on new and published mitochondrial cytochrome b sequences. Monophyly of Fringillidae and its two constituent subfamilies, Fringillinae and Emberizinae, was consistently supported with the exceptions of Peucedramus being placed outside of Fringillinae and Euphonia being placed within Fringillinae instead of within Emberizinae. Within Emberizinae, Thraupini (tanagers), Cardinalini (cardinals and grosbeaks), and Emberizini (New World sparrows) did not form separate monophyletic groups. Our results indicate that Emberizinae consists of three clades, each with a different overall geographical distribution. Several taxa traditionally considered members of Thraupini fall outside of the thraupine clade, including the only North American genus, Piranga. Consequently, the thraupine clade includes only Neotropical species. Increasing evidence suggests that Fringillidae, often called "New World nine-primaried oscines," does not in fact have a New World origin.

In a previous study, we found avian sarcoma and leukosis virus (ASLV) gag genes in 19 species of birds in the order Galliformes including all grouse and ptarmigan (Tetraoninae) surveyed. Our data suggested that retroviruses had been transmitted horizontally among some host species. To further investigate these elements, we sequenced a replication-defective retrovirus, here named tetraonine endogenous retrovirus (TERV), from Bonasa umbellus (ruffed grouse). This is the first report of a complete, replication-defective ASLV provirus sequence from any bird other than the domestic chicken. We found a replication-defective proviral sequence consisting of putative Gag and Env proteins flanked by long terminal repeats. Reverse transcription-PCR analysis showed that retroviral gag sequences closely related to TERV are transcribed, supporting the hypothesis that TERV is an active endogenous retrovirus. Phylogenetic analyses suggest that TERV may have arisen via recombination between different retroviral lineages infecting birds. Southern blotting using gag probes showed that TERV occurs in tetraonines but not in chickens or ducks, suggesting that integration occurred after the earliest phasianid divergences but prior to the radiation of tetraonine birds.

Maximum likelihood analysis, accounting for site-heterogeneity in evolutionary rate with the Gamma-distribution model, was carried out with amino acid sequences of 12 mitochondrial proteins and nucleotide sequences of mitochondrial 12S and 16S rRNAs from three turtles, one squamate, one crocodile, and eight birds. The analysis strongly suggests that turtles are closely related to archosaurs (birds+crocodilians), and it supports both Tree-2: (((birds, crocodilians), turtles), squamates) and Tree-3: ((birds, (crocodilians, turtles)), squamates). A more traditional Tree-1: (((birds, crocodilians), squamates), turtles) and a tree in which turtles are basal to other amniotes were rejected with high statistical significance. Tree-3 has recently been proposed by Hedges and Poling [Science 283 (1999) 998-1001] based mainly on nuclear genes. Therefore, we re-analyzed their data using the maximum likelihood method, and evaluated the total evidence of the analyses of mitochondrial and nuclear data sets. Tree-1 was again rejected strongly. The most likely hypothesis was Tree-3, though Tree-2 remained a plausible candidate.

We describe a previously unrecognized form of gene homology using the term "gametology," which we define as homology arising through lack of recombination and subsequent differentiation of sex chromosomes. We demonstrate use of gametologous genes to root each other in phylogenetic analyses of sex-specific avian Chromo-helicase-DNA binding gene (CHD) sequences. Phylogenetic analyses of a set of neognath bird sequences yield monophyletic groups for CHD-W and CHD-Z gametologs, as well as congruent relationships between these two clades and between them and current views of avian taxonomy. Phylogenetic analyses including paleognath bird CHD sequences and rooting with crocodilian CHD sequences, suggest an early divergence for paleognath CHD within the avian CHD clade. Based on our CHD analyses calibrated with avian fossil dates, we estimate the divergence between CHD-W and CHD-Z at 123 MYA, suggesting an early differentiation of sex chromosomes that predates most extant avian orders. In agreement with the notion of male-driven evolution, we find a faster rate of change in male-linked CHD-Z sequences.

In a study of the evolution and distribution of avian retroviruses, we found avian sarcoma and leukosis virus (ASLV) gag genes in 26 species of galliform birds from North America, Central America, eastern Europe, Asia, and Africa. Nineteen of the 26 host species from whom ASLVs were sequenced were not previously known to contain ASLVs. We assessed congruence between ASLV phylogenies based on a total of 110 gag gene sequences and ASLV-host phylogenies based on mitochondrial 12S ribosomal DNA and ND2 sequences to infer coevolutionary history for ASLVs and their hosts. Widespread distribution of ASLVs among diverse, endemic galliform host species suggests an ancient association. Congruent ASLV and host phylogenies for two species of Perdix, two species of Gallus, and Lagopus lagopus and L. mutus also indicate an old association with vertical transmission and cospeciation for these ASLVs and hosts. An inference of horizontal transmission of ASLVs among some members of the Tetraoninae subfamily (grouse and ptarmigan) is supported by ASLV monophyletic groups reflecting geographic distribution and proximity of hosts rather than host species phylogeny. We provide a preliminary phylogenetic taxonomy for the new ASLVs, in which named taxa denote monophyletic groups.

An adjustable fitness model for amino acid site substitutions is investigated. This model, a generalization of previously developed evolutionary models, has several distinguishing characteristics: it separately accounts for the processes of mutation and substitution, allows for heterogeneity among substitution rates and among evolutionary constraints, and does not make any prior assumptions about which sites or characteristics of proteins are important to molecular evolution. While the model has fewer adjustable parameters than the general reversible mtREV model, when optimized it outperforms mtREV in likelihood analysis on protein-coding mitochondrial genes. In addition, the optimized fitness parameters of the model show correspondence to some biophysical characteristics of amino acids.

A PCR-based approach to sequencing complete mitochondrial genomes is described along with a set of 86 primers designed primarily for avian mitochondrial DNA (mtDNA). This PCR-based approach allows an accurate determination of complete mtDNA sequences that is faster than sequencing cloned mtDNA. The primers are spaced at about 500-base intervals along both DNA strands. Many of the primers incorporate degenerate positions to accommodate variation in mtDNA sequence among avian taxa and to reduce the potential for preferential amplification of nuclear pseudogenes. Comparison with published vertebrate mtDNA sequences suggests that many of the primers will have broad taxonomic utility. In addition, these primers should make available a wider variety of mitochondrial genes for studies based on smaller data sets., (Copyright 1999 Academic Press.)

Using the set of all vertebrate mtDNA protein sequences published as of May 1998, plus unpublished examples for elephant and birds, we examined divergence times in Placentalia and Aves. Using a parsimony-based test, we identified a subset of slower evolutionary rate placental sequences that do not appear to violate the clock assumption. Analyzing just these sequences decreases support for Marsupionta and the carnivore + perissodactyl group but increases support for armadillo diverging earlier than rabbit (which may represent the whole Glires group). A major theme of the paper is to use more comprehensive estimates of divergence time standard error (SE). From the well-studied horse/rhino split, estimated to be 55 million years before present (mybp), the splitting time within carnivores is confidently shown to be older than 50 million years. Some of our estimates of divergence times within placentals are relatively old, at up to 169 million years, but are within 2 SE of other published estimates. The whale/cow split at 65 mybp may be older than commonly assumed. All the sampled splits between the main groups of fereuungulates (the clade of carnivores, cetartiodactyls, perissodactyls, and pholidotes) seem to be distinctly before the Cretaceous/Tertiary boundary. Analyses suggest a close relationship between elephants (representing Afrotheria) and armadillos (Xenarthra), and our timing of this splitting is coincident with the opening of the South Atlantic, a major vicariant event. Recalibrating with this event (at 100 mybp), we obtain younger estimates for the earliest splits among placentals. Divergence times within birds are also assessed by using previously unpublished sequences. We fail to reject a clock for all bird taxa available. Unfortunately, available deep calibration points for birds are questionable, so a new calibration based on the age of the Anseriform stem lineage is estimated. The divergence time of rhea and ostrich may be much more recent than commonly assumed, while that of passerines may be older. Our major concern is the rooting point of the bird subtree, as the nearest outgroup (alligator) is very distant.

Several different groups of birds have been proposed as being the oldest or earliest diverging extant lineage within the avian phylogenetic tree, particularly ratites (Struthioniformes), waterfowl (Anseriformes), and shorebirds (Charadriiformes). Difficulty in resolving this issue stems from several factors, including the relatively rapid radiation of primary (ordinal) bird lineages and the lack of characters from an extant outgroup for birds that is closely related to them by measure of time. To help resolve this question, we have sequenced entire mitochondrial genomes for five birds (a rhea, a duck, a falcon, and two perching birds), one crocodilian, and one turtle. Maximum parsimony and maximum likelihood analyses of these new sequences together with published sequences (18 taxa total) yield the same optimal tree topology, in which a perching bird (Passeriformes) is sister to all the other bird taxa. A basal position for waterfowl among the bird study taxa is rejected by maximum likelihood analyses. However, neither the conventional view, in which ratites (including rhea) are basal to other birds, nor tree topologies with falcon or chicken basal among birds could be rejected in the same manner. In likelihood analyses of a subset of seven birds, alligator, and turtle (9 taxa total), we find that increasing the number of parameters in the model shifts the optimal topology from one with a perching bird basal among birds to the conventional view with ratites diverging basally; moreover, likelihood scores for the two trees are not significantly different. Thus, although our largest set of taxa and characters supports a tree with perching birds diverging basally among birds, the position of this earliest divergence among birds appears unstable. Our analyses indicate a sister relationship between a waterfowl/chicken clade and ratites, relative to perching birds and falcon. We find support for a sister relationship between turtles and a bird/crocodilian clade, and for rejecting both the Haemothermia hypothesis (birds and mammals as sister taxa) and the placement of turtles as basal within the phylogenetic tree for amniote animals.

HIV-1 subtype phylogeny is investigated using a previously developed computational model of natural amino acid site substitutions. This model, based on Boltzmann statistics and Metropolis kinetics, involves an order of magnitude fewer adjustable parameters than traditional substitution matrices and deals more effectively with the issue of protein site heterogeneity. When optimized for sequences of HIV-1 envelope (env) proteins from a few specific subtypes, our model is more likely to describe the evolutionary record for other subtypes than are methods using a single substitution matrix, even a matrix optimized over the same data. Pairwise distances are calculated between various probabilistic ancestral subtype sequences, and a distance matrix approach is used to find the optimal phylogenetic tree. Our results indicate that the relationships between subtypes B, C, and D and those between subtypes A and H may be closer than previously thought.

Mitochondrial genomes of all vertebrate animals analyzed to date have the same 37 genes, whose arrangement in the circular DNA molecule varies only in the relative position of a few genes. This relative conservation suggests that mitochondrial gene order characters have potential utility as phylogenetic markers for higher-level vertebrate taxa. We report discovery of a mitochondrial gene order that has had multiple independent originations within birds, based on sampling of 137 species representing 13 traditionally recognized orders. This provides evidence of parallel evolution in mitochondrial gene order for animals. Our results indicate operation of physical constraints on mitochondrial gene order changes and support models for gene order change based on replication error. Bird mitochondria have a displaced OL (origin of light-strand replication site) as do various other Reptilia taxa prone to gene order changes. Our findings point to the need for broad taxonomic sampling in using mitochondrial gene order for phylogenetic analyses. We found, however, that the alternative mitochondrial gene orders distinguish the two primary groups of songbirds (order Passeriformes), oscines and suboscines, in agreement with other molecular as well as morphological data sets. Thus, although mitochondrial gene order characters appear susceptible to some parallel evolution because of mechanistic constraints, they do hold promise for phylogenetic studies.