Louisiana State University LSU Digital Commons Faculty Publications Department of Biological Sciences 11-12-2020 Dense sampling of bird diversity increases power of comparative genomics Shaohong Feng China National Genebank Josefin Stiller Københavns Universitet Yuan Deng China National Genebank Joel Armstrong University of California, Santa Cruz Qi Fang China National Genebank See next page for additional authors Follow this and additional works at: https://digitalcommons.lsu.edu/biosci_pubs Recommended Citation Feng, S., Stiller, J., Deng, Y., Armstrong, J., Fang, Q., Reeve, A., Xie, D., Chen, G., Guo, C., Faircloth, B., Petersen, B., Wang, Z., Zhou, Q., Diekhans, M., Chen, W., Andreu-Sánchez, S., Margaryan, A., Howard, J., Parent, C., Pacheco, G., Sinding, M., Puetz, L., Cavill, E., Ribeiro, Â., Eckhart, L., Fjeldså, J., Hosner, P., Brumfield, R., Christidis, L., Bertelsen, M., Sicheritz-Ponten, T., Tietze, D., Robertson, B., & Song, G. (2020). Dense sampling of bird diversity increases power of comparative genomics. Nature, 587 (7833), 252-257. https://doi.org/10.1038/s41586-020-2873-9 This Article is brought to you for free and open access by the Department of Biological Sciences at LSU Digital Commons. It has been accepted for inclusion in Faculty Publications by an authorized administrator of LSU Digital Commons. For more information, please contact [email protected]. Authors Shaohong Feng, Josefin Stiller, Yuan Deng, Joel Armstrong, Qi Fang, Andrew Hart Reeve, Duo Xie, Guangji Chen, Chunxue Guo, Brant C. Faircloth, Bent Petersen, Zongji Wang, Qi Zhou, Mark Diekhans, Wanjun Chen, Sergio Andreu-Sánchez, Ashot Margaryan, Jason Travis Howard, Carole Parent, George Pacheco, Mikkel Holger S. Sinding, Lara Puetz, Emily Cavill, Ângela M. Ribeiro, Leopold Eckhart, Jon Fjeldså, Peter A. Hosner, Robb T. Brumfield, Les Christidis, Mads .F Bertelsen, Thomas Sicheritz-Ponten, Dieter Thomas Tietze, Bruce C. Robertson, and Gang Song This article is available at LSU Digital Commons: https://digitalcommons.lsu.edu/biosci_pubs/487 Article Dense sampling of bird diversity increases power of comparative genomics https://doi.org/10.1038/s41586-020-2873-9 A list of authors and affiliations appears at the end of the paper. Received: 9 August 2019 Whole-genome sequencing projects are increasingly populating the tree of life and Accepted: 27 July 2020 characterizing biodiversity1–4. Sparse taxon sampling has previously been proposed Published online: 11 November 2020 to confound phylogenetic inference5, and captures only a fraction of the genomic Open access 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 Check for updates 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 confdently 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-specifc variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will ofer new perspectives on evolutionary processes in cross-species comparative analyses and assist in eforts to conserve species. Comparative genomics is rapidly growing, fuelled by the advancement (Pedionomus torquatus) and the Bali myna (Leucopsar rothschildi, which of sequencing technologies. Many large-scale initiatives have been has fewer than 50 adults remaining in the wild13). proposed with a core mission of producing genomes for hundreds of Two hundred and sixty-seven of the 363 species represented in our species, representing the phylogenetic diversity of particular taxa6–8. genome data are newly released, comprising 18.4 trillion base pairs Although the generation of genomes is now more routine, an immediate (bp) of raw data and 284 billion bp of assemblies. The assemblies are challenge is how to efficiently compare large numbers of genomes in comparable in quality to previously published bird genomes9,10, but an evolutionary context. A critical first step is the accurate detection of vary in contiguity (average scaffold N50 = 1.42 megabases (Mb), contig orthologous sequences. In this study, we release a large-scale dataset of N50 = 42.57 kilobases (kb); see interactive supplementary figure 1, bird genomes, which we use to establish a framework for comparative hosted at https://genome-b10k.herokuapp.com/main). The sequencing analysis. We provide insight on how scaling-up genome sampling assists coverage ranged from 35× (blue-throated roller (Eurystomus gularis) in our understanding of avian genomic diversity and in the detection and yellowhead (Mohoua ochrocephala)) to 368× (song sparrow (Melo- of signals of natural selection down to individual bases. spiza melodia)) and genomic completeness was high (average 95.8%). The B10K Project began with the Avian Phylogenomics Consortium We annotated all 363 genomes using a homology-based method with a (phase I), which analysed 48 genomes from representatives of most uniform gene set that included gene models from chicken, zebra finch bird orders9,10. Here we report the genome sequencing outcomes from and human, and published transcriptomes (Supplementary Tables 6– phase II of the project: these outcomes include a total of 363 species in 8), to predict an average of 15,464 protein-coding genes for each species 92.4% (218 out of 23611,12) of avian families (Supplementary Tables 1–5). (Supplementary Table 1). We also assembled mitochondrial genomes Species were selected to span the overall diversity and to subdivide long for 336 species, with 216 species fully circularized and 228 species branches, when possible (Fig. 1, Supplementary Data). Our sampling with a complete mitochondrial annotation (Supplementary Table 1). covers bird species from every continent (Extended Data Fig. 1) and Bird genomes at the ordinal level were previously found to contain a more than triples the previous taxonomic coverage of avian genome low proportion of transposable elements, except for the downy wood- sequencing; to our knowledge, 155 bird families are represented here pecker (Picoides pubescens) in Piciformes10. Consistent with these for the first time. We chose short-read sequencing as our main strategy findings, 96.1% of birds at the family level had a transposable element for generating data, which enabled us to use older samples (the oldest of content lower than 15%—but we found additional outliers (Extended which was collected in 1982) and access rare museum specimens—such Data Fig. 2a, Supplementary Table 1). In particular, long interspersed as one of the few vouchered tissues of the Henderson crake (Zapornia nuclear elements were prevalent in all nine sequenced species in Pici- atra), which occurs on a single island. We incorporated 68 species of formes, which suggests an ancestral expansion in this lineage (24% on concern on the International Union for Conservation of Nature (IUCN) average, Welch two-sample t-test, P = 9.98 × 10−5) (Extended Data Fig. 2b, Red List of Threatened Species (Supplementary Table 1); these include d). The common scimitarbill (Rhinopomastus cyanomelas) and com- 12 endangered and 2 critically endangered species—the plains-wanderer mon hoopoe (Upupa epops) in Bucerotiformes also had exceptionally 252 | Nature | Vol 587 | 12 November 2020 Uria lomvia Arenaria interpres Fregetta grallaria Aquila chrysaetos Ramphastos sulfuratus Alopecoenas beccarii Agapornis roseicollis Piaya cayana Calypte anna * Odontophorus gujanensis Onychorhynchus coronatus Apteryx owenii Orthonyx spaldingii Dyaphorophyia castanea ** Agelaius phoeniceus Chaetops frenatus Tichodroma muraria Leptocoma aspasia Hemignathus wilsoni Fig. 1 | Newly sequenced genomes densely cover the bird tree of life. The per sequenced family (purple branches). The grey arc marks the diverse 10,135 bird species11,12 are shown on a draft phylogeny that synthesizes Passeriformes radiation, with 6,063 species, of which 173 species have genome taxonomic and phylogenetic information36 (Supplementary Data). In total, assemblies now. Chicken (*) and zebra finch (**) are marked for orientation. 363 species, covering 92.4% of all families, now have at least 1 genome assembly Paintings illustrate examples of sequenced species. high transposable element content (23% and 18%, respectively) owing Project phase I. Nonetheless, 498 genes remained absent across all to recent expansions of long interspersed nuclear elements, whereas 363 bird species, which adds to evidence that these genes were truly two hornbill species in the same order exhibited the typical low propor- lost in the common ancestor. tions (Extended Data Fig. 2e). We also investigated a number of genes that were previously associ- Previous studies have suggested that hundreds of genes were lost ated with phenotypes and physiological pathways. For example, we in the ancestor of birds10,14. Gene-loss inference is complicated by found that rhodopsin (encoded by RH1) and the