Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Species-specific traits mediate avian demographic responses under past climate change.

Anticipating species' responses to environmental change is a pressing mission in biodiversity conservation. Despite decades of research investigating how climate change may affect population sizes, historical context is lacking, and the traits that mediate demographic sensitivity to changing climate remain elusive. We use whole-genome sequence data to reconstruct the demographic histories of 263 bird species over the past million years and identify networks of interacting morphological and life history traits associated with changes in effective population size (Ne) in response to climate warming and cooling. Our results identify direct and indirect effects of key traits representing dispersal, reproduction and survival on long-term demographic responses to climate change, thereby highlighting traits most likely to influence population responses to ongoing climate warming.

Changes in the functional diversity of modern bird species over the last million years.

Despite evidence of declining biosphere integrity, we currently lack understanding of how the functional diversity associated with changes in abundance among ecological communities has varied over time and before widespread human disturbances. We combine morphological, ecological, and life-history trait data for >260 extant bird species with genomic-based estimates of changing effective population size (Ne) to quantify demographic-based shifts in avian functional diversity over the past million years and under pre-anthropogenic climate warming. We show that functional diversity was relatively stable over this period, but underwent significant changes in some key areas of trait space due to changing species abundances. Our results suggest that patterns of population decline over the Pleistocene have been concentrated in particular regions of trait space associated with extreme reproductive strategies and low dispersal ability, consistent with an overall erosion of functional diversity. Further, species most sensitive to climate warming occupied a relatively narrow region of functional space, indicating that the largest potential population increases and decreases under climate change will occur among species with relatively similar trait sets. Overall, our results identify fluctuations in functional space of extant species over evolutionary timescales and represent the demographic-based vulnerability of different regions of functional space among these taxa. The integration of paleodemographic dynamics with functional trait data enhances our ability to quantify losses of biosphere integrity before anthropogenic disturbances and attribute contemporary biodiversity loss to different drivers over time.

Draft genome assemblies of four manakins.

Manakins are a family of small suboscine passerine birds characterized by their elaborate courtship displays, non-monogamous mating system, and sexual dimorphism. This family has served as a good model for the study of sexual selection. Here we present genome assemblies of four manakin species, including Cryptopipo holochlora, Dixiphia pipra (also known as Pseudopipra pipra), Machaeropterus deliciosus and Masius chrysopterus, generated by Single-tube Long Fragment Read (stLFR) technology. The assembled genome sizes ranged from 1.10 Gb to 1.19 Gb, with average scaffold N50 of 29 Mb and contig N50 of 169 Kb. On average, 12,055 protein-coding genes were annotated in the genomes, and 9.79% of the genomes were annotated as repetitive elements. We further identified 75 Mb of Z-linked sequences in manakins, containing 585 to 751 genes and an ~600 Kb pseudoautosomal region (PAR). One notable finding from these Z-linked sequences is that a possible Z-to-autosome/PAR reversal could have occurred in M. chrysopterus. These de novo genomes will contribute to a deeper understanding of evolutionary history and sexual selection in manakins.

Historical specimens and the limits of subspecies phylogenomics in the New World quails (Odontophoridae).

As phylogenomics focuses on comprehensive taxon sampling at the species and population/subspecies levels, incorporating genomic data from historical specimens has become increasingly common. While historical samples can fill critical gaps in our understanding of the evolutionary history of diverse groups, they also introduce additional sources of phylogenomic uncertainty, making it difficult to discern novel evolutionary relationships from artifacts caused by sample quality issues. These problems highlight the need for improved strategies to disentangle artifactual patterns from true biological signal as historical specimens become more prevalent in phylogenomic datasets. Here, we tested the limits of historical specimen-driven phylogenomics to resolve subspecies-level relationships within a highly polytypic family, the New World quails (Odontophoridae), using thousands of ultraconserved elements (UCEs). We found that relationships at and above the species-level were well-resolved and highly supported across all analyses, with the exception of discordant relationships within the two most polytypic genera which included many historical specimens. We examined the causes of discordance and found that inferring phylogenies from subsets of taxa resolved the disagreements, suggesting that analyzing subclades can help remove artifactual causes of discordance in datasets that include historical samples. At the subspecies-level, we found well-resolved geographic structure within the two most polytypic genera, including the most polytypic species in this family, Northern Bobwhites (Colinus virginianus), demonstrating that variable sites within UCEs are capable of resolving phylogenetic structure below the species level. Our results highlight the importance of complete taxonomic sampling for resolving relationships among polytypic species, often through the inclusion of historical specimens, and we propose an integrative strategy for understanding and addressing the uncertainty that historical samples sometimes introduce to phylogenetic analyses.

Population connectivity across a highly fragmented distribution: Phylogeography of the Chalcophaps doves.

Chalcophaps is a morphologically conserved genus of ground-walking doves distributed from India to mainland China, south to Australia, and across the western Pacific to Vanuatu. Here, we reconstruct the evolutionary history of this genus using DNA sequence data from two nuclear genes and one mitochondrial gene, sampled from throughout the geographic range of Chalcophaps. We find support for three major evolutionary lineages in our phylogenetic reconstruction, each corresponding to the three currently recognized Chalcophaps species. Despite this general concordance, we identify discordant mitochondrial and nuclear ancestries in the subspecies C. longirostris timorensis, raising further questions about the evolutionary history of this Timor endemic population. Within each of the three species, we find evidence for isolation by distance or hierarchical population structure, indicating an important role for geography in the diversification of this genus. Despite being distributed broadly across a highly fragmented geographic region known as a hotspot for avian diversification, the Chalcophaps doves show modest levels of phenotypic and genetic diversity, a pattern potentially explained by strong population connectivity owing to high overwater dispersal capability.

Divergence time estimation of Galliformes based on the best gene shopping scheme of ultraconserved elements.

BACKGROUND: Divergence time estimation is fundamental to understanding many aspects of the evolution of organisms, such as character evolution, diversification, and biogeography. With the development of sequence technology, improved analytical methods, and knowledge of fossils for calibration, it is possible to obtain robust molecular dating results. However, while phylogenomic datasets show great promise in phylogenetic estimation, the best ways to leverage the large amounts of data for divergence time estimation has not been well explored. A potential solution is to focus on a subset of data for divergence time estimation, which can significantly reduce the computational burdens and avoid problems with data heterogeneity that may bias results. RESULTS: In this study, we obtained thousands of ultraconserved elements (UCEs) from 130 extant galliform taxa, including representatives of all genera, to determine the divergence times throughout galliform history. We tested the effects of different "gene shopping" schemes on divergence time estimation using a carefully, and previously validated, set of fossils. Our results found commonly used clock-like schemes may not be suitable for UCE dating (or other data types) where some loci have little information. We suggest use of partitioning (e.g., PartitionFinder) and selection of tree-like partitions may be good strategies to select a subset of data for divergence time estimation from UCEs. Our galliform time tree is largely consistent with other molecular clock studies of mitochondrial and nuclear loci. With our increased taxon sampling, a well-resolved topology, carefully vetted fossil calibrations, and suitable molecular dating methods, we obtained a high quality galliform time tree. CONCLUSIONS: We provide a robust galliform backbone time tree that can be combined with more fossil records to further facilitate our understanding of the evolution of Galliformes and can be used as a resource for comparative and biogeographic studies in this group.

When good mitochondria go bad: Cyto-nuclear discordance in landfowl (Aves: Galliformes).

Mitochondrial sequences were among the first molecular data collected for phylogenetic studies and they are plentiful in DNA sequence archives. However, the future value of mitogenomic data in phylogenetics is uncertain, because its phylogenetic signal sometimes conflicts with that of the nuclear genome. A thorough understanding of the causes and prevalence of cyto-nuclear discordance would aid in reconciling different results owing to sequence data type, and provide a framework for interpreting megaphylogenies when taxa which lack substantial nuclear data are placed using mitochondrial data. Here, we examine the prevalence and possible causes of cyto-nuclear discordance in the landfowl (Aves: Galliformes), leveraging 47 new mitogenomes assembled from off-target reads recovered as part of a target-capture study. We evaluated two hypotheses, that cyto-nuclear discordance is "genuine" and a result of biological processes such as incomplete lineage sorting or introgression, and that cyto-nuclear discordance is an artifact of inaccurate mitochondrial tree estimation (the "inaccurate estimation" hypothesis). We identified seven well-supported topological differences between the mitogenomic tree and trees based on nuclear data. These well-supported topological differences were robust to model selection. An examination of sites suggests these differences were driven by small number of sites, particularly from third-codon positions, suggesting that they were not confounded by convergent directional selection. Hence, the hypothesis of genuine discordance was supported.

A phylogenomic supermatrix of Galliformes (Landfowl) reveals biased branch lengths.

Building taxon-rich phylogenies is foundational for macroevolutionary studies. One approach to improve taxon sampling beyond individual studies is to build supermatricies of publicly available data, incorporating taxa sampled across different studies and utilizing different loci. Most existing supermatrix studies have focused on loci commonly sequenced with Sanger technology ("legacy" markers, such as mitochondrial data and small numbers of nuclear loci). However, incorporating phylogenomic studies into supermatrices allows problem nodes to be targeted and resolved with considerable amounts of data, while improving taxon sampling with legacy data. Here we estimate phylogeny from a galliform supermatrix which includes well-known model and agricultural species such as the chicken and turkey. We assembled a supermatrix comprising 4500 ultra-conserved elements (UCEs) collected as part of recent phylogenomic studies in this group and legacy mitochondrial and nuclear (intron and exon) sequences. Our resulting phylogeny included 88% of extant species and recovered well-accepted relationships with strong support. However, branch lengths, which are particularly important in down-stream macroevolutionary studies, appeared vastly skewed. Taxa represented only by rapidly evolving mitochondrial data had high proportions of missing data and exhibited long terminal branches. Conversely, taxa sampled for slowly evolving UCEs with low proportions of missing data exhibited substantially shorter terminal branches. We explored several branch length re-estimation methods with particular attention to terminal branches and conclude that re-estimation using well-sampled mitochondrial sequences may be a pragmatic approach to obtain trees suitable for macroevolutionary analysis.

Phylogenomics of manakins (Aves: Pipridae) using alternative locus filtering strategies based on informativeness.

Target capture sequencing effectively generates molecular marker arrays useful for molecular systematics. These extensive data sets are advantageous where previous studies using a few loci have failed to resolve relationships confidently. Moreover, target capture is well-suited to fragmented source DNA, allowing data collection from species that lack fresh tissues. Herein we use target capture to generate data for a phylogeny of the avian family Pipridae (manakins), a group that has been the subject of many behavioral and ecological studies. Most manakin species feature lek mating systems, where males exhibit complex behavioral displays including mechanical and vocal sounds, coordinated movements of multiple males, and high speed movements. We analyzed thousands of ultraconserved element (UCE) loci along with a smaller number of coding exons and their flanking regions from all but one species of Pipridae. We examined three different methods of phylogenetic estimation (concatenation and two multispecies coalescent methods). Phylogenetic inferences using UCE data yielded strongly supported estimates of phylogeny regardless of analytical method. Exon probes had limited capability to capture sequence data and resulted in phylogeny estimates with reduced support and modest topological differences relative to the UCE trees, although these conflicts had limited support. Two genera were paraphyletic among all analyses and data sets, with Antilophia nested within Chiroxiphia and Tyranneutes nested within Neopelma. The Chiroxiphia-Antilophia clade was an exception to the generally high support we observed; the topology of this clade differed among analyses, even those based on UCE data. To further explore relationships within this group, we employed two filtering strategies to remove low-information loci. Those analyses resulted in distinct topologies, suggesting that the relationships we identified within Chiroxiphia-Antilophia should be interpreted with caution. Despite the existence of a few continuing uncertainties, our analyses resulted in a robust phylogenetic hypothesis of the family Pipridae that provides a comparative framework for future ecomorphological and behavioral studies.

Dense sampling of bird diversity increases power of comparative genomics.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

Taxonomic evaluation of the Grallaria rufula (Rufous Antpitta) complex (Aves: Passeriformes: Grallariidae) distinguishes sixteen species.

Populations in the Rufous Antpitta (Grallaria rufula) complex occupy humid montane forests of the Andes from northern Colombia and adjacent Venezuela to central Bolivia. Their tawny to cinnamon-colored plumages are generally uniform, featuring subtle variation in hue and saturation across this range. In contrast to their conservative plumage, substantial vocal differences occur among geographically isolated or parapatric populations. Working within the framework of a comprehensive molecular phylogeny, we reexamined species limits in the G. rufula complex, basing taxonomic recommendations on diagnostic differences in vocalizations and considering identifiable differences in plumage where pertinent. We identified 16 populations for species designation, including seven populations previously described as subspecies and, remarkably, six new species described herein. Within one of these species, we identified less robust vocal differences between populations that we designate as subspecies. Geographic variation exists within another species, but its critical evaluation requires additional material. Taxonomic revisions of groups consisting of cryptic species, like the Grallaria rufula complex, are imperative for their conservation. Rather than widespread species as currently defined, these complexes can comprise many range-restricted taxa at higher risk of extinction given the continuing human pressures on their habitats.

Acknowledging uncertainty in evolutionary reconstructions of ecological niches.

Reconstructing ecological niche evolution can provide insight into the biogeography and diversification of evolving lineages. However, comparative phylogenetic methods may infer the history of ecological niche evolution inaccurately because (a) species' niches are often poorly characterized; and (b) phylogenetic comparative methods rely on niche summary statistics rather than full estimates of species' environmental tolerances. Here, we propose a new framework for coding ecological niches and reconstructing their evolution that explicitly acknowledges and incorporates the uncertainty introduced by incomplete niche characterization. Then, we modify existing ancestral state inference methods to leverage full estimates of environmental tolerances. We provide a worked empirical example of our method, investigating ecological niche evolution in the New World orioles (Aves: Passeriformes: Icterus spp.). Temperature and precipitation tolerances were generally broad and conserved among orioles, with niche reduction and specialization limited to a few terminal branches. Tools for performing these reconstructions are available in a new R package called nichevol.

Genomic differentiation in an endemic Philippine genus (Aves: Sarcophanops) owing to geographical isolation on recently disassociated islands.

Phylogeographical studies of Philippine vertebrates have demonstrated that genetic variation is broadly partitioned by Pleistocene island aggregation. Contemporary island discontinuity is expected to influence genetic differentiation but remains relatively undocumented, perhaps because the current episode of island isolation started in relatively recent times. We investigated inter- and intra-island population structure in a Philippine endemic bird genus (Sarcophanops) to determine whether genetic differentiation has evolved during the recent period of isolation. We sequenced thousands of genome-wide restriction site associated DNA (RAD) markers from throughout the Mindanao group to assess fine-scale genetic structure across islands. Specifically, we investigated patterns of gene flow and connectivity within and between taxonomic and geographical bounds. A previous assessment of mitochondrial DNA detected deep structure between Sarcophanops samarensis and a sister species, Sarcophanops steerii, but was insufficient to detect differentiation within either species. Analysis of RAD markers, however, revealed structure within S. samarensis between the islands of Samar/Leyte and Bohol. This genetic differentiation probably demonstrates an effect of recent geographical isolation (after the Last Glacial Maximum) on the genetic structure of Philippine avifauna. We suggest that the general lack of evidence for differentiation between recently isolated populations is a failure to detect subtle population structure owing to past genetic sampling constraints, rather than the absence of such structure.

Correction to: Progress on bringing together raptor collections in Europe for contaminant research and monitoring in relation to chemicals regulation.

The correct affiliation of Sabrina Lo Brutto is shown in this paper.

Earth history and the passerine superradiation.

Avian diversification has been influenced by global climate change, plate tectonic movements, and mass extinction events. However, the impact of these factors on the diversification of the hyperdiverse perching birds (passerines) is unclear because family level relationships are unresolved and the timing of splitting events among lineages is uncertain. We analyzed DNA data from 4,060 nuclear loci and 137 passerine families using concatenation and coalescent approaches to infer a comprehensive phylogenetic hypothesis that clarifies relationships among all passerine families. Then, we calibrated this phylogeny using 13 fossils to examine the effects of different events in Earth history on the timing and rate of passerine diversification. Our analyses reconcile passerine diversification with the fossil and geological records; suggest that passerines originated on the Australian landmass ∼47 Ma; and show that subsequent dispersal and diversification of passerines was affected by a number of climatological and geological events, such as Oligocene glaciation and inundation of the New Zealand landmass. Although passerine diversification rates fluctuated throughout the Cenozoic, we find no link between the rate of passerine diversification and Cenozoic global temperature, and our analyses show that the increases in passerine diversification rate we observe are disconnected from the colonization of new continents. Taken together, these results suggest more complex mechanisms than temperature change or ecological opportunity have controlled macroscale patterns of passerine speciation.

Molecular systematics of swifts of the genus Chaetura (Aves: Apodiformes: Apodidae).

Phylogenetic relationships among swifts of the morphologically conservative genus Chaetura were studied using mitochondrial and nuclear DNA sequences. Taxon sampling included all species and 21 of 30 taxa (species and subspecies) within Chaetura. Our results indicate that Chaetura is monophyletic and support the division of the genus into the two subgenera previously identified using plumage characters. However, our genetic data, when considered in combination with phenotypic data, appear to be at odds with the current classification of some species of Chaetura. We recommend that C. viridipennis, currently generally treated as specifically distinct from C. chapmani, be returned to its former status as C. chapmani viridipennis, and that C. andrei, now generally regarded as synonymous with C. vauxi aphanes, again be recognized as a valid species. Widespread Neotropical species C. spinicaudus is paraphyletic with respect to more range-restricted species C. fumosa, C. egregia, and C. martinica. Geographically structured genetic variation within some other species of Chaetura, especially notable in C. cinereiventris, suggests that future study may lead to recognition of additional species in this genus. Biogeographic analysis indicated that Chaetura originated in South America and identified several dispersal events to Middle and North America following the formation of the Isthmus of Panama.

Colonization and diversification of the white-browed shortwing (Aves: Muscicapidae: Brachypteryx montana) in the Philippines.

Molecular phylogenetic approaches have greatly improved our knowledge of the pattern and process of biological diversification across the globe; however, many regions remain poorly documented, even for well-studied vertebrate taxa. The Philippine archipelago, one of the least-studied 'biodiversity hotspots', is an ideal natural laboratory for investigating the factors driving diversification in an insular and geologically dynamic setting. We investigated the history and geography of diversification of the Philippine populations of a widespread montane bird, the White-browed Shortwing (Brachypteryx montana). Leveraging dense archipelago-wide sampling, we generated a multi-locus genetic dataset (one nuclear and two mtDNA markers), which we analyzed using phylogenetic, population genetic, and coalescent-based methods. Our results demonstrate that Philippine shortwings (1) likely colonized the Philippines from the Sunda Shelf to Mindanao in the late Miocene or Pliocene, (2) diversified across inter-island barriers into three divergent lineages during the Pliocene and early Pleistocene, (3) have not diversified within the largest island, Luzon, contrary to patterns observed in other montane taxa, and (4) colonized Palawan from the oceanic Philippines rather than from Borneo, challenging the assumption of Palawan functioning exclusively as a biogeographic extension of the Sunda Shelf. Additionally, our finding that divergent (c. 4.0 mya) lineages are coexisting in secondary sympatry on Mindanao without apparent gene flow suggests that the speciation process is likely complete for these shortwing lineages. Overall, these investigations provide insight into how topography and island boundaries influence diversification within remote oceanic archipelagos and echo the results of many other studies in demonstrating that taxonomic diversity continues to be underestimated in the Philippines.