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.].

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.

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.

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.

Comparison of beak and feather disease virus prevalence and immunity-associated genetic diversity over time in an island population of red-crowned parakeets.

Pathogen outbreaks in the wild can contribute to a population's extinction risk. Concern over the effects of pathogen outbreaks in wildlife is amplified in small, threatened populations, where degradation of genetic diversity may hinder natural selection for enhanced immunocompetence. Beak and feather disease virus (BFDV) was detected for the first time in an island population of red-crowned parakeets (Cyanoramphus novaezelandiae) in 2008 on Little Barrier Island (Hauturu-o-Toi) of New Zealand. By 2013, the prevalence of the viral infection had significantly decreased within the population. We tested whether the population of red-crowned parakeets showed a selective response to BFDV, using neutral microsatellite and two immunity-associated genetic markers, the major histocompatibility complex (MHC) and Toll-like receptors (TLRs). We found evidence for selection at viral-associated TLR3; however, the ability of TLR3 to elicit an immune response in the presence of BFDV warrants confirmation. Alternatively, because red-crowned parakeet populations are prone to fluctuations in size, the decrease in BFDV prevalence over time may be attributed to the Little Barrier Island population dropping below the density threshold for viral maintenance. Our results highlight that natural processes such as adaptation for enhanced immunocompetence and/or density fluctuations are efficient mechanisms for reducing pathogen prevalence in a threatened, isolated population.

European Colonization, Not Polynesian Arrival, Impacted Population Size and Genetic Diversity in the Critically Endangered New Zealand Kakapo.

Island endemic species are often vulnerable to decline and extinction following human settlement, and the genetic study of historical museum specimens can be useful in understanding these processes. The kakapo (Strigops habroptilus) is a critically endangered New Zealand parrot that was formerly widespread and abundant. It is well established that both Polynesian and European colonization of New Zealand impacted the native avifauna, but the timeframe and severity of impacts have differed depending on species. Here, we investigated the relative importance of the 2 waves of human settlement on kakapo decline, using microsatellites and mitochondrial DNA (mtDNA) to characterize recent kakapo genetic and demographic history. We analyzed samples from 49 contemporary individuals and 54 museum specimens dating from 1884 to 1985. Genetic diversity decreased significantly between historical and contemporary kakapo, with a decline in mean number of microsatellite alleles from 6.15 to 3.08 and in number of mtDNA haplotypes from 17 to 3. Modeling of demographic history indicated a recent population bottleneck linked to the period of European colonization (approximately 5 generations ago) but did not support a major decline linked to Polynesian settlement. Effective population size estimates were also larger for historical than contemporary kakapo. Our findings inform contemporary kakapo management by indicating the timeframe and possible cause of the bottleneck, which has implications for the management of extant genetic diversity. We demonstrate the broader utility of a historical perspective in understanding causes of decline and managing extinction risk in contemporary endangered species.

Genetic Evidence of a Population Bottleneck and Inbreeding in the Endangered New Zealand Sea Lion, Phocarctos hookeri.

The New Zealand sea lion (NZSL) is of high conservation concern due to its limited distribution and its declining population size. Historically, it occupied most of coastal New Zealand, but is now restricted to a few coastal sites in southern mainland New Zealand and the sub-Antarctic Islands. NZSLs have experienced a recent reduction in population size due to sealing in the 1900s, which is expected to have resulted in increased inbreeding and a loss of genetic variation, potentially reducing the evolutionary capacity of the species and negatively impacting on its long-term prospects for survival. We used 17 microsatellite loci, previously shown to have cross-species applications in pinnipeds, to determine locus- and population-specific statistics for 1205 NZSLs from 7 consecutive breeding seasons. We show that the NZSL population has a moderate level of genetic diversity in comparison to other pinnipeds. We provide genetic evidence for a population reduction, likely caused by historical sealing, and a measure of allele sharing/parental relatedness (internal relatedness) that is suggestive of increased inbreeding in pups that died during recent epizootic episodes. We hypothesize that population bottlenecks and nonrandom mating have impacted on the population genetic architecture of NZSLs, affecting its population recovery.

Low Spatial Genetic Differentiation Associated with Rapid Recolonization in the New Zealand Fur Seal Arctocephalus forsteri.

Population declines resulting from anthropogenic activities are of major consequence for the long-term survival of species because the resulting loss of genetic diversity can lead to extinction via the effects of inbreeding depression, fixation of deleterious mutations, and loss of adaptive potential. Otariid pinnipeds have been exploited commercially to near extinction with some species showing higher demographic resilience and recolonization potential than others. The New Zealand fur seal (NZFS) was heavily impacted by commercial sealing between the late 18th and early 19th centuries, but has recolonized its former range in southern Australia. The species has also recolonized its former range in New Zealand, yet little is known about the pattern of recolonization. Here, we first used 11 microsatellite markers (n = 383) to investigate the contemporary population structure and dispersal patterns in the NZFS (Arctocephalus forsteri). Secondly, we model postsealing recolonization with 1 additional mtDNA cytochrome b (n = 261) marker. Our data identified 3 genetic clusters: an Australian, a subantarctic, and a New Zealand one, with a weak and probably transient subdivision within the latter cluster. Demographic history scenarios supported a recolonization of the New Zealand coastline from remote west coast colonies, which is consistent with contemporary gene flow and with the species' high resilience. The present data suggest the management of distinct genetic units in the North and South of New Zealand along a genetic gradient. Assignment of individuals to their colony of origin was limited (32%) with the present data indicating the current microsatellite markers are unlikely sufficient to assign fisheries bycatch of NZFSs to colonies.