Towards complete and error-free genome assemblies of all vertebrate species.

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

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.

Ancient mitochondrial genomes unveil the origins and evolutionary history of New Zealand's enigmatic takahe and moho.

Many avian species endemic to Aotearoa New Zealand were driven to extinction or reduced to relict populations following successive waves of human arrival, due to hunting, habitat destruction and the introduction of mammalian predators. Among the affected species were the large flightless South Island takahe (Porphyrio hochstetteri) and the moho (North Island takahe; P. mantelli), with the latter rendered extinct and the former reduced to a single relictual population. Little is known about the evolutionary history of these species prior to their decline and/or extinction. Here we sequenced mitochondrial genomes from takahe and moho subfossils (12 takahe and 4 moho) and retrieved comparable sequence data from takahe museum skins (n = 5) and contemporary individuals (n = 17) to examine the phylogeny and recent evolutionary history of these species. Our analyses suggest that prehistoric takahe populations lacked deep phylogeographic structure, in contrast to moho, which exhibited significant spatial genetic structure, albeit based on limited sample sizes (n = 4). Temporal genetic comparisons show that takahe have lost much of their mitochondrial genetic diversity, likely due to a sudden demographic decline soon after human arrival (~750 years ago). Time-calibrated phylogenetic analyses strongly support a sister species relationship between takahe and moho, suggesting these flightless taxa diverged around 1.5 million years ago, following a single colonisation of New Zealand by a flighted Porphyrio ancestor approximately 4 million years ago. This study highlights the utility of palaeogenetic approaches for informing the conservation and systematic understanding of endangered species whose ranges have been severely restricted by anthropogenic impacts.

Candidate gene polymorphisms are linked to dispersive and migratory behaviour: Searching for a mechanism behind the "paradox of the great speciators".

The "paradox of the great speciators" has puzzled evolutionary biologists for over half a century. A great speciator requires excellent dispersal propensity to explain its occurrence on multiple islands, but reduced dispersal ability to explain its high number of subspecies. A rapid reduction in dispersal ability is often invoked to solve this apparent paradox, but a proximate mechanism has not been identified yet. Here, we explored the role of six genes linked to migration and animal personality differences (CREB1, CLOCK, ADCYAP1, NPAS2, DRD4, and SERT) in 20 South Pacific populations of silvereye (Zosterops lateralis) that range from highly sedentary to partially migratory, to determine if genetic variation is associated with dispersal propensity and migration. We detected genetic associations in three of the six genes: (i) in a partial migrant population, migrant individuals had longer microsatellite alleles at the CLOCK gene compared to resident individuals from the same population; (ii) CREB1 displayed longer average microsatellite allele lengths in recently colonized island populations (<200 years), compared to evolutionarily older populations. Bayesian broken stick regression models supported a reduction in CREB1 length with time since colonization; and (iii) like CREB1, DRD4 showed differences in polymorphisms between recent and old colonizations but a larger sample is needed to confirm. ADCYAP1, SERT, and NPAS2 were variable but that variation was not associated with dispersal propensity. The association of genetic variants at three genes with migration and dispersal ability in silvereyes provides the impetus for further exploration of genetic mechanisms underlying dispersal shifts, and the prospect of resolving a long-running evolutionary paradox through a genetic lens.

Genomic trajectories of a near-extinction event in the Chatham Island black robin.

BACKGROUND: Understanding the micro--evolutionary response of populations to demographic declines is a major goal in evolutionary and conservation biology. In small populations, genetic drift can lead to an accumulation of deleterious mutations, which will increase the risk of extinction. However, demographic recovery can still occur after extreme declines, suggesting that natural selection may purge deleterious mutations, even in extremely small populations. The Chatham Island black robin (Petroica traversi) is arguably the most inbred bird species in the world. It avoided imminent extinction in the early 1980s and after a remarkable recovery from a single pair, a second population was established and the two extant populations have evolved in complete isolation since then. Here, we analysed 52 modern and historical genomes to examine the genomic consequences of this extreme bottleneck and the subsequent translocation. RESULTS: We found evidence for two-fold decline in heterozygosity and three- to four-fold increase in inbreeding in modern genomes. Moreover, there was partial support for temporal reduction in total load for detrimental variation. In contrast, compared to historical genomes, modern genomes showed a significantly higher realised load, reflecting the temporal increase in inbreeding. Furthermore, the translocation induced only small changes in the frequency of deleterious alleles, with the majority of detrimental variation being shared between the two populations. CONCLUSION: Our results highlight the dynamics of mutational load in a species that recovered from the brink of extinction, and show rather limited temporal changes in mutational load. We hypothesise that ancestral purging may have been facilitated by population fragmentation and isolation on several islands for thousands of generations and may have already reduced much of the highly deleterious load well before human arrival and introduction of pests to the archipelago. The majority of fixed deleterious variation was shared between the modern populations, but translocation of individuals with low mutational load could possibly mitigate further fixation of high-frequency deleterious variation.

Phylogenomic Discordance in the Eared Seals is best explained by Incomplete Lineage Sorting following Explosive Radiation in the Southern Hemisphere.

The phylogeny and systematics of fur seals and sea lions (Otariidae) have long been studied with diverse data types, including an increasing amount of molecular data. However, only a few phylogenetic relationships have reached acceptance because of strong gene-tree species tree discordance. Divergence times estimates in the group also vary largely between studies. These uncertainties impeded the understanding of the biogeographical history of the group, such as when and how trans-equatorial dispersal and subsequent speciation events occurred. Here, we used high-coverage genome-wide sequencing for 14 of the 15 species of Otariidae to elucidate the phylogeny of the family and its bearing on the taxonomy and biogeographical history. Despite extreme topological discordance among gene trees, we found a fully supported species tree that agrees with the few well-accepted relationships and establishes monophyly of the genus Arctocephalus. Our data support a relatively recent trans-hemispheric dispersal at the base of a southern clade, which rapidly diversified into six major lineages between 3 and 2.5 Ma. Otaria diverged first, followed by Phocarctos and then four major lineages within Arctocephalus. However, we found Zalophus to be nonmonophyletic, with California (Zalophus californianus) and Steller sea lions (Eumetopias jubatus) grouping closer than the Galapagos sea lion (Zalophus wollebaeki) with evidence for introgression between the two genera. Overall, the high degree of genealogical discordance was best explained by incomplete lineage sorting resulting from quasi-simultaneous speciation within the southern clade with introgresssion playing a subordinate role in explaining the incongruence among and within prior phylogenetic studies of the family. [Hybridization; ILS; phylogenomics; Pleistocene; Pliocene; monophyly.].

Historical reconstruction unveils the risk of mass mortality and ecosystem collapse during pancontinental megadrought.

An important new hypothesis in landscape ecology is that extreme, decade-scale megadroughts can be potent drivers of rapid, macroscale ecosystem degradation and collapse. If true, an increase in such events under climate change could have devastating consequences for global biodiversity. However, because few megadroughts have occurred in the modern ecological era, the taxonomic breadth, trophic depth, and geographic pattern of these impacts remain unknown. Here we use ecohistorical techniques to quantify the impact of a record, pancontinental megadrought period (1891 to 1903 CE) on the Australian biota. We show that during this event mortality and severe stress was recorded in >45 bird, mammal, fish, reptile, and plant families in arid, semiarid, dry temperate, and Mediterranean ecosystems over at least 2.8 million km2 (36%) of the Australian continent. Trophic analysis reveals a bottom-up pattern of mortality concentrated in primary producer, herbivore, and omnivore guilds. Spatial and temporal reconstruction of premortality rainfall shows that mass mortality and synchronous ecosystem-wide collapse emerged in multiple geographic hotspots after 2 to 4 y of severe (>40%) and intensifying rainfall deficits. However, the presence of hyperabundant herbivores significantly increased the sensitivity of ecosystems to overgrazing-induced meltdown and permanent ecosystem change. The unprecedented taxonomic breadth and spatial scale of these impacts demonstrate that continental-scale megadroughts pose a major future threat to global biodiversity, especially in ecosystems affected by intensive agricultural use, trophic simplification, and invasive species.

Evolution of the "world's only alpine parrot": Genomic adaptation or phenotypic plasticity, behaviour and ecology?

Climate warming, in particular in island environments, where opportunities for species to disperse are limited, may become a serious threat to cold adapted alpine species. In order to understand how alpine species may respond to a warming world, we need to understand the drivers that have shaped their habitat specialisation and the evolutionary adaptations that allow them to utilize alpine habitats. The endemic, endangered New Zealand kea (Nestor notabilis) is considered the only alpine parrot in the world. As a species commonly found in the alpine zone it may be highly susceptible to climate warming. But is it a true alpine specialist? Is its evolution driven by adaptation to the alpine zone, or is the kea an open habitat generalist that simply uses the alpine zone to, for example, avoid lower lying anthropogenic landscapes? We use whole genome data of the kea and its close, forest adapted sister species, the kaka (Nestor meridionalis) to reconstruct the evolutionary history of both species and identify the functional genomic differences that underlie their habitat specialisations. Our analyses do not identify major functional genomic differences between kea and kaka in pathways associated with high-altitude. Rather, we found evidence that selective pressures on adaptations commonly found in alpine species are present in both Nestor species, suggesting that selection for alpine adaptations has not driven their divergence. Strongly divergent demographic responses to past climate warming between the species nevertheless highlight potential future threats to kea survival in a warming world.

An island-hopping bird reveals how founder events shape genome-wide divergence.

When populations colonize new areas, both strong selection and strong drift can be experienced due to novel environments and small founding populations, respectively. Empirical studies have predominantly focused on the phenotype when assessing the role of selection, and limited neutral-loci when assessing founder-induced loss of diversity. Consequently, the extent to which processes interact to influence evolutionary trajectories is difficult to assess. Genomic-level approaches provide the opportunity to simultaneously consider these processes. Here, we examine the roles of selection and drift in shaping genomic diversity and divergence in historically documented sequential island colonizations by the silvereye (Zosterops lateralis). We provide the first empirical demonstration of the rapid appearance of highly diverged genomic regions following population founding, the position of which are highly idiosyncratic. As these regions rarely contained loci putatively under selection, it is most likely that these differences arise via the stochastic nature of the founding process. However, selection is required to explain rapid evolution of larger body size in insular silvereyes. Reconciling our genomic data with these phenotypic patterns suggests there may be many genomic routes to the island phenotype, which vary across populations. Finally, we show that accelerated divergence associated with multiple founding steps is the product of genome-wide rather than localized differences, and that diversity erodes due to loss of rare alleles. However, even multiple founder events do not result in divergence and diversity levels seen in evolutionary older subspecies, and therefore do not provide a shortcut to speciation as proposed by founder-effect speciation models.

Marine mammal population decline linked to obscured by-catch.

Declines of marine megafauna due to fisheries by-catch are thought to be mitigated by exclusion devices that release nontarget species. However, exclusion devices may instead conceal negative effects associated with by-catch caused by fisheries (i.e., unobserved or discarded by-catch with low postrelease survival or reproduction). We show that the decline of the endangered New Zealand (NZ) sea lion (Phocarctos hookeri) is linked to latent levels of by-catch occurring in sub-Antarctic trawl fisheries. Exclusion devices have been used since 2001 but have not slowed or reversed population decline. However, 35% of the variability in NZ sea lion pup production is explained by latent by-catch, and the population would increase without this factor. Our results indicate that exclusion devices can obscure rather than alleviate fishery impacts on marine megafauna.

Understanding Maladaptation by Uniting Ecological and Evolutionary Perspectives.

Evolutionary biologists have long trained their sights on adaptation, focusing on the power of natural selection to produce relative fitness advantages while often ignoring changes in absolute fitness. Ecologists generally have taken a different tack, focusing on changes in abundance and ranges that reflect absolute fitness while often ignoring relative fitness. Uniting these perspectives, we articulate various causes of relative and absolute maladaptation and review numerous examples of their occurrence. This review indicates that maladaptation is reasonably common from both perspectives, yet often in contrasting ways. That is, maladaptation can appear strong from a relative fitness perspective, yet populations can be growing in abundance. Conversely, resident individuals can appear locally adapted (relative to nonresident individuals) yet be declining in abundance. Understanding and interpreting these disconnects between relative and absolute maladaptation, as well as the cases of agreement, is increasingly critical in the face of accelerating human-mediated environmental change. We therefore present a framework for studying maladaptation, focusing in particular on the relationship between absolute and relative fitness, thereby drawing together evolutionary and ecological perspectives. The unification of these ecological and evolutionary perspectives has the potential to bring together previously disjunct research areas while addressing key conceptual issues and specific practical problems.

Genomics detects population structure within and between ocean basins in a circumpolar seabird: The white-chinned petrel.

The Southern Ocean represents a continuous stretch of circumpolar marine habitat, but the potential physical and ecological drivers of evolutionary genetic differentiation across this vast ecosystem remain unclear. We tested for genetic structure across the full circumpolar range of the white-chinned petrel (Procellaria aequinoctialis) to unravel the potential drivers of population differentiation and test alternative population differentiation hypotheses. Following range-wide comprehensive sampling, we applied genomic (genotyping-by-sequencing or GBS; 60,709 loci) and standard mitochondrial-marker approaches (cytochrome b and first domain of control region) to quantify genetic diversity within and among island populations, test for isolation by distance, and quantify the number of genetic clusters using neutral and outlier (non-neutral) loci. Our results supported the multi-region hypothesis, with a range of analyses showing clear three-region genetic population structure, split by ocean basin, within two evolutionary units. The most significant differentiation between these regions confirmed previous work distinguishing New Zealand and nominate subspecies. Although there was little evidence of structure within the island groups of the Indian or Atlantic oceans, a small set of highly-discriminatory outlier loci could assign petrels to ocean basin and potentially to island group, though the latter needs further verification. Genomic data hold the key to revealing substantial regional genetic structure within wide-ranging circumpolar species previously assumed to be panmictic.

Complete genomes of two extinct New Zealand passerines show responses to climate fluctuations but no evidence for genomic erosion prior to extinction.

Human intervention, pre-human climate change (or a combination of both), as well as genetic effects, contribute to species extinctions. While many species from oceanic islands have gone extinct due to direct human impacts, the effects of pre-human climate change and human settlement on the genomic diversity of insular species and the role that loss of genomic diversity played in their extinctions remains largely unexplored. To address this question, we sequenced whole genomes of two extinct New Zealand passerines, the huia (Heteralocha acutirostris) and South Island kokako (Callaeas cinereus). Both species showed similar demographic trajectories throughout the Pleistocene. However, the South Island kokako continued to decline after the last glaciation, while the huia experienced some recovery. Moreover, there was no indication of inbreeding resulting from recent mating among closely related individuals in either species. This latter result indicates that population fragmentation associated with forest clearing by Maori may not have been strong enough to lead to an increase in inbreeding and exposure to genomic erosion. While genomic erosion may not have directly contributed to their extinctions, further habitat fragmentation and the introduction of mammalian predators by Europeans may have been an important driver of extinction in huia and South Island kokako.

Reciprocal translocation of small numbers of inbred individuals rescues immunogenetic diversity.

Genetic rescue can reduce inbreeding depression and increase fitness of small populations, even when the donor populations are highly inbred. In a recent experiment involving two inbred island populations of the New Zealand South Island robin, Petroica australis, reciprocal translocations improved microsatellite diversity and individual fitness. While microsatellite loci may reflect patterns of genome-wide diversity, they generally do not indicate the specific genetic regions responsible for increased fitness. We tested the effectiveness of this reciprocal translocation for rescuing diversity of two immunogenetic regions: Toll-like receptor (TLR) and major histocompatibility complex (MHC) genes. We found that the relatively small number of migrants (seven and ten per island) effectively brought the characteristic TLR gene diversity of each source population into the recipient population. However, when migrants transmitted TLR alleles that were already present at high frequency in the recipient population, it was possible for offspring of mixed heritage to have decreased gene diversity compared to recipient population diversity prior to translocation. In contrast to TLRs, we did not observe substantial changes in MHC allelic diversity following translocation, with limited evidence of a decrease in differentiation, perhaps because most MHC alleles were observed at both sites prior to the translocation. Overall, we conclude that small numbers of migrants may successfully restore the diversity of immunogenetic loci with few alleles, but that translocating larger numbers of animals would provide additional opportunity for the genetic rescue of highly polymorphic immunity regions, such as the MHC, even when the source population is inbred.

Perennial grasslands enhance biodiversity and multiple ecosystem services in bioenergy landscapes.

Agriculture is being challenged to provide food, and increasingly fuel, for an expanding global population. Producing bioenergy crops on marginal lands--farmland suboptimal for food crops--could help meet energy goals while minimizing competition with food production. However, the ecological costs and benefits of growing bioenergy feedstocks--primarily annual grain crops--on marginal lands have been questioned. Here we show that perennial bioenergy crops provide an alternative to annual grains that increases biodiversity of multiple taxa and sustain a variety of ecosystem functions, promoting the creation of multifunctional agricultural landscapes. We found that switchgrass and prairie plantings harbored significantly greater plant, methanotrophic bacteria, arthropod, and bird diversity than maize. Although biomass production was greater in maize, all other ecosystem services, including methane consumption, pest suppression, pollination, and conservation of grassland birds, were higher in perennial grasslands. Moreover, we found that the linkage between biodiversity and ecosystem services is dependent not only on the choice of bioenergy crop but also on its location relative to other habitats, with local landscape context as important as crop choice in determining provision of some services. Our study suggests that bioenergy policy that supports coordinated land use can diversify agricultural landscapes and sustain multiple critical ecosystem services.

Human-mediated extirpation of the unique Chatham Islands sea lion and implications for the conservation management of remaining New Zealand sea lion populations.

While terrestrial megafaunal extinctions have been well characterized worldwide, our understanding of declines in marine megafauna remains limited. Here, we use ancient DNA analyses of prehistoric (<1450-1650 AD) sea lion specimens from New Zealand's isolated Chatham Islands to assess the demographic impacts of human settlement. These data suggest there was a large population of sea lions, unique to the Chatham Islands, at the time of Polynesian settlement. This distinct mitochondrial lineage became rapidly extinct within 200 years due to overhunting, paralleling the extirpation of a similarly large endemic mainland population. Whole mitogenomic analyses confirm substantial intraspecific diversity among prehistoric lineages. Demographic models suggest that even low harvest rates would likely have driven rapid extinction of these lineages. This study indicates that surviving Phocarctos populations are remnants of a once diverse and widespread sea lion assemblage, highlighting dramatic human impacts on endemic marine biodiversity. Our findings also suggest that Phocarctos bycatch in commercial fisheries may contribute to the ongoing population decline.

Myth or relict: Does ancient DNA detect the enigmatic Upland seal?

The biological status of the so-called 'Upland seal' has remained contentious ever since historical records described a distinct seal from the uplands of New Zealand's (NZ) remote sub-Antarctic islands. Subsequent genetic surveys of the NZ fur seal (Arctocephalus forsteri) detected two highly-divergent mtDNA clades, hypothesized to represent a post-sealing hybrid swarm between 'mainland' (Australia-NZ; A. forsteri) and sub-Antarctic (putative 'Upland'; A. snaresensis) lineages. We present ancient-DNA analyses of prehistoric mainland NZ and sub-Antarctic fur seals, revealing that both of these genetic lineages were already widely distributed across the region at the time of human arrival. These findings indicate that anthropogenic factors did not contribute to the admixture of these lineages, and cast doubt on the validity of the Upland seal. Human-mediated impacts on Arctocephalus genetic diversity are instead highlighted by a dramatic temporal haplotype frequency-shift due to genetic drift in heavily bottlenecked populations following the cessation of industrial-scale harvesting. These extinction-recolonisation dynamics add to a growing picture of human-mediated change in NZ's coastal and marine ecosystems.

Phylogenetic affinities of the Fregetta storm-petrels are not black and white.

The Fregetta storm-petrels generally are regarded to comprise two species: black-bellied storm-petrels F. tropica (monotypic) breed at Antarctic and sub-Antarctic islands (46-63°S), and white-bellied storm-petrels F. grallaria breed at south temperate islands (28-37°S), with four recognized subspecies. Confusion surrounds the status of birds at Gough Island (40°S), central South Atlantic, which have been attributed usually to a white-bellied form of black-bellied storm-petrel F. t. melanoleuca. We use cytochrome b and nuclear ß-fibrinogen gene sequences to show that F. t. melanoleuca are present during the breeding season at Gough and islands in the nearby Tristan da Cunha archipelago (37°S), exhibiting limited divergence from F. t. tropica. We also show that there is greater diversity among F. grallaria populations, with eastern South Pacific F. g. segethi and F. g. titan differing by c. 0.011, and both differing from western South Pacific nominate F. g. grallaria by c. 0.059. The Tristan archipelago supports a population of F. grallaria closely allied to the nominate form, as well as a distinct form identified as F. g. leucogaster. Further research is needed to assess how F. grallaria and F. tropica segregate in sympatry at Tristan and Gough, and why this is the only location where both species have white-bellies.