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.

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.

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.

MHC variation reflects the bottleneck histories of New Zealand passerines.

Most empirical evidence suggests that balancing selection does not counter the effects of genetic drift in shaping postbottleneck major histocompatibility complex (MHC) genetic diversity when population declines are severe or prolonged. However, few studies have been able to include data from historical specimens, or to compare populations/species with different bottleneck histories. In this study, we examined MHC class II B and microsatellite diversity in four New Zealand passerine (songbird) species that experienced moderate to very severe declines. We compared diversity from historical samples (collected c. 1884-1938) to present-day populations. Using a Bayesian framework, we found that the change in genetic diversity from historical to contemporary samples was affected by three main factors: (i) whether the data were based on MHC or microsatellite markers, (ii) species (as a surrogate for bottleneck severity) and (iii) whether the comparison between historical and contemporary samples was made using historical samples originating from the mainland, or using historical samples originating from islands. The greatest losses in genetic diversity occurred for the most severely bottlenecked species, particularly between historical mainland and contemporary samples. Additionally, where loss of diversity occurred, the change was greater for MHC genes compared to microsatellite loci.

Evidence for Bergmann's Rule and Not Allopatric Subspeciation in the Threatened Kaka (Nestor meridionalis).

Species of conservation concern characterized by small and declining populations greatly benefit from proactive management approaches such as population translocations. Because they often show intra-specific genetic and phenotypic variation, which can result from drift or differential selective pressures between habitats, understanding the distribution of such variation and its underlying processes is a prerequisite to develop effective management guidelines. Indeed, translocations among genetically differentiated populations potentially locally adapted are discouraged in order to avoid outbreeding depression, while translocations among populations characterized by high gene flow with no evidence for local adaptation are encouraged. Here, we first test whether 2 recognized subspecies, the North Island kaka (Nestor meridionalis septentrionalis) and South Island kaka (Nestor meridionalis meridionalis) of New Zealand fit a scenario of allopatric subspeciation following the separation of the North and South Islands at the end of the Pleistocene using 1 mtDNA (n = 96) and 9 microsatellite markers (n = 126). We then test whether morphological differences among the 2 subspecies support a pattern of local adaptation, comparing phenotypic divergence (P ST) and the level of divergence by drift alone (F ST) among populations. We find little population structure between islands, ruling out allopatric subspeciation in kaka. Further, P ST exceeds F ST, supporting an adaptive latitudinal size cline consistent with Bergmann's rule. These results therefore suggest that using neutral genetic diversity alone can be misleading when identifying management units and that the nature of phenotypic variation should be considered in translocations efforts. We finally discuss North and South Island management units but suggest that cross-island translocation be allowed.

Evidence for multiple MHC class II ß loci in New Zealand's critically endangered kakapo, Strigops habroptilus.

Immunologically important genes of the major histocompatibility complex (MHC) have been characterized in a number of avian species with the general finding of considerable variation in size and structural organization among organisms. A range of nonpasserines which represent early-diverging Neoave lineages have been described as having only one MHC class II ß locus potentially leading to the conclusion that this is the ancestral condition. Here, we examine the monotypic, early-diverging, critically endangered kakapo, Strigops habroptilus, for allelic variation at MHC class II ß exon 2, as part of species' recovery efforts. We found two to four confirmed sequence variants per individual indicating the presence of more than one MHC class II ß locus. Given the kakapo's basal evolutionary status, evidence for multiple MHC class II ß loci seems to counter the proposed mono-locus history of modern birds. However, MHC gene duplication, maintenance, and loss among and within bird species may confound avian relationships making it difficult to elucidate the ancestral state. This study adds essential data for disentangling the course of MHC structural evolution in birds.

Characterization of MHC class II B polymorphism in bottlenecked New Zealand saddlebacks reveals low levels of genetic diversity.

The major histocompatibility complex (MHC) is integral to the vertebrate adaptive immune system. Characterizing diversity at functional MHC genes is invaluable for elucidating patterns of adaptive variation in wild populations, and is particularly interesting in species of conservation concern, which may suffer from reduced genetic diversity and compromised disease resilience. Here, we use next generation sequencing to investigate MHC class II B (MHCIIB) diversity in two sister taxa of New Zealand birds: South Island saddleback (SIS), Philesturnus carunculatus, and North Island saddleback (NIS), Philesturnus rufusater. These two species represent a passerine family outside the more extensively studied Passerida infraorder, and both have experienced historic bottlenecks. We examined exon 2 sequence data from populations that represent the majority of genetic diversity remaining in each species. A high level of locus co-amplification was detected, with from 1 to 4 and 3 to 12 putative alleles per individual for South and North Island birds, respectively. We found strong evidence for historic balancing selection in peptide-binding regions of putative alleles, and we identified a cluster combining non-classical loci and pseudogene sequences from both species, although no sequences were shared between the species. Fewer total alleles and fewer alleles per bird in SIS may be a consequence of their more severe bottleneck history; however, overall nucleotide diversity was similar between the species. Our characterization of MHCIIB diversity in two closely related species of New Zealand saddlebacks provides an important step in understanding the mechanisms shaping MHC diversity in wild, bottlenecked populations.

Genetic drift outweighs natural selection at toll-like receptor (TLR) immunity loci in a re-introduced population of a threatened species.

During population establishment, genetic drift can be the key driver of changes in genetic diversity, particularly while the population is small. However, natural selection can also play a role in shaping diversity at functionally important loci. We used a well-studied, re-introduced population of the threatened Stewart Island robin (N = 722 pedigreed individuals) to determine whether selection shaped genetic diversity at innate immunity toll-like receptor (TLR) genes, over a 9-year period of population growth following establishment with 12 genetic founders. We found no evidence for selection operating with respect to TLR diversity on first-year overwinter survival for the majority of loci, genotypes and alleles studied. However, survival of individuals with TLR4BE genotype was significantly improved: these birds were less than half as likely to die prior to maturity compared with all other TLR4 genotypes. Furthermore, the population frequency of this genotype, at a two-fold excess over Hardy-Weinberg expectation, was increased by nonrandom mating. Near-complete sampling and full pedigree and reproductive data enabled us to eliminate other potential causes of these patterns including inbreeding, year effects, density dependence, selection on animals at earlier life history stages or genome-level association of the TLR4E allele with 'good genes'. However, comparison of observed levels of gene diversity to predictions under simulated genetic drift revealed results consistent with neutral expectations for all loci, including TLR4. Although selection favoured TLR4BE heterozygotes in this population, these effects were insufficient to outweigh genetic drift. This is the first empirical study to show that genetic drift can overwhelm natural selection in a wild population immediately following establishment.

Simulating retention of rare alleles in small populations to assess management options for species with different life histories.

Preserving allelic diversity is important because it provides the capacity for adaptation and thus enables long-term population viability. Allele retention is difficult to predict in animals with overlapping generations, so we used a new computer model to simulate retention of rare alleles in small populations of 3 species with contrasting life-history traits: North Island Brown Kiwi (Apteryx mantelli; monogamous, long-lived), North Island Robins (Petroica longipes; monogamous, short-lived), and red deer (Cervus elaphus; polygynous, moderate lifespan). We simulated closed populations under various demographic scenarios and assessed the amounts of artificial immigration needed to achieve a goal of retaining 90% of selectively neutral rare alleles (frequency in the source population = 0.05) after 10 generations. The number of immigrants per generation required to meet the genetic goal ranged from 11 to 30, and there were key similarities and differences among species. None of the species met the genetic goal without immigration, and red deer lost the most allelic diversity due to reproductive skew among polygynous males. However, red deer required only a moderate rate of immigration relative to the other species to meet the genetic goal because nonterritorial breeders had a high turnover. Conversely, North Island Brown Kiwi needed the most immigration because the long lifespan of locally produced territorial breeders prevented a large proportion of immigrants from recruiting. In all species, the amount of immigration needed generally decreased with an increase in carrying capacity, survival, or reproductive output and increased as individual variation in reproductive success increased, indicating the importance of accurately quantifying these parameters to predict the effects of management. Overall, retaining rare alleles in a small, isolated population requires substantial investment of management effort. Use of simulations to explore strategies optimized for the populations in question will help maximize the value of this effort..

Disentangling the roles of natural selection and genetic drift in shaping variation at MHC immunity genes.

The major histocompatibility complex (MHC) forms an integral component of the vertebrate immune response and, due to strong selection pressures, is one of the most polymorphic regions of the entire genome. Despite over 15 years of research, empirical studies offer highly contradictory explanations of the relative roles of different evolutionary forces, selection and genetic drift, acting on MHC genes during population bottlenecks. Here, we take a meta-analytical approach to quantify the results of studies into the effects of bottlenecks on MHC polymorphism. We show that the consequences of selection acting on MHC loci prior to a bottleneck event, combined with drift during the bottleneck, will result in overall loss of MHC polymorphism that is ∼15% greater than loss of neutral genetic diversity. These results are counter to general expectations that selection should maintain MHC polymorphism, but do agree with the results of recent simulation models and at least two empirical studies. Notably, our results suggest that negative frequency-dependent selection could be more important than overdominance for maintaining high MHC polymorphism in pre-bottlenecked populations.

The imprecision of heterozygosity-fitness correlations hinders the detection of inbreeding and inbreeding depression in a threatened species.

In nonpedigreed wild populations, inbreeding depression is often quantified through the use of heterozygosity-fitness correlations (HFCs), based on molecular estimates of relatedness. Although such correlations are typically interpreted as evidence of inbreeding depression, by assuming that the marker heterozygosity is a proxy for genome-wide heterozygosity, theory predicts that these relationships should be difficult to detect. Until now, the vast majority of empirical research in this area has been performed on generally outbred, nonbottlenecked populations, but differences in population genetic processes may limit extrapolation of results to threatened populations. Here, we present an analysis of HFCs, and their implications for the interpretation of inbreeding, in a free-ranging pedigreed population of a bottlenecked species: the endangered takahe (Porphyrio hochstetteri). Pedigree-based inbreeding depression has already been detected in this species. Using 23 microsatellite loci, we observed only weak evidence of the expected relationship between multilocus heterozygosity and fitness at individual life-history stages (such as survival to hatching and fledging), and parameter estimates were imprecise (had high error). Furthermore, our molecular data set could not accurately predict the inbreeding status of individuals (as 'inbred' or 'outbred', determined from pedigrees), nor could we show that the observed HFCs were the result of genome-wide identity disequilibrium. These results may be attributed to high variance in heterozygosity within inbreeding classes. This study is an empirical example from a free-ranging endangered species, suggesting that even relatively large numbers (>20) of microsatellites may give poor precision for estimating individual genome-wide heterozygosity. We argue that pedigree methods remain the most effective method of quantifying inbreeding in wild populations, particularly those that have gone through severe bottlenecks.

Founder effects, inbreeding, and loss of genetic diversity in four avian reintroduction programs.

The number of individuals translocated and released as part of a reintroduction is often small, as is the final established population, because the reintroduction site is typically small. Small founder and small resulting populations can result in population bottlenecks, which are associated with increased rates of inbreeding and loss of genetic diversity, both of which can affect the long-term viability of reintroduced populations. I used information derived from pedigrees of four monogamous bird species reintroduced onto two different islands (220 and 259 ha) in New Zealand to compare the pattern of inbreeding and loss of genetic diversity among the reintroduced populations. Although reintroduced populations founded with few individuals had higher levels of inbreeding, as predicted, other factors, including biased sex ratio and skewed breeding success, contributed to high levels of inbreeding and loss of genetic diversity. Of the 10-58 individuals released, 4-25 genetic founders contributed at least one living descendent and yielded approximately 3-11 founder-genome equivalents (number of genetic founders assuming an equal contribution of offspring and no random loss of alleles across generations) after seven breeding seasons. This range is much lower than the 20 founder-genome equivalents recommended for captive-bred populations. Although the level of inbreeding in one reintroduced population initially reached three times that of a closely related species, the long-term estimated rate of inbreeding of this one population was approximately one-third that of the other species due to differences in carrying capacities of the respective reintroduction sites. The increasing number of reintroductions to suitable areas that are smaller than those I examined here suggests that it might be useful to develop long-term strategies and guidelines for reintroduction programs, which would minimize inbreeding and maintain genetic diversity.

Inbreeding depression accumulation across life-history stages of the endangered Takahe.

Studies evaluating the impact of inbreeding depression on population viability of threatened species tend to focus on the effects of inbreeding at a single life-history stage (e.g., juvenile survival). We examined the effects of inbreeding across the full life-history continuum, from survival up to adulthood, to subsequent reproductive success, and to the recruitment of second-generation offspring, in wild Takahe ( Porphyrio hochstetteri ) by analyzing pedigree and fitness data collected over 21 breeding seasons. Although the effect size of inbreeding at individual life-history stages was small, inbreeding depression accumulated across multiple life-history stages and ultimately reduced long-term fitness (i.e., successful recruitment of second-generation offspring). The estimated total lethal equivalents (2B) summed across all life-history stages were substantial (16.05, 95% CI 0.08-90.8) and equivalent to an 88% reduction in recruitment of second-generation offspring for closely related pairs (e.g., sib-sib pairings) relative to unrelated pairs (according to the pedigree). A history of small population size in the Takahe could have contributed to partial purging of the genetic load and the low level of inbreeding depression detected at each single life-history stage. Nevertheless, our results indicate that such "purged" populations can still exhibit substantial inbreeding depression, especially when small but negative fitness effects accumulate across the species' life history. Because inbreeding depression can ultimately affect population viability of small, isolated populations, our results illustrate the importance of measuring the effects of inbreeding across the full life-history continuum.

Unexpected positive and negative effects of continuing inbreeding in one of the world's most inbred wild animals.

Inbreeding depression, the reduced fitness of offspring of related individuals, is a central theme in evolutionary biology. Inbreeding effects are influenced by the genetic makeup of a population, which is driven by any history of genetic bottlenecks and genetic drift. The Chatham Island black robin represents a case of extreme inbreeding following two severe population bottlenecks. We tested whether inbreeding measured by a 20-year pedigree predicted variation in fitness among individuals, despite the high mean level of inbreeding and low genetic diversity in this species. We found that paternal and maternal inbreeding reduced fledgling survival and individual inbreeding reduced juvenile survival, indicating that inbreeding depression affects even this highly inbred population. Close inbreeding also reduced survival for fledglings with less-inbred mothers, but unexpectedly improved survival for fledglings with highly inbred mothers. This counterintuitive interaction could not be explained by various potentially confounding variables. We propose a genetic mechanism, whereby a highly inbred chick with a highly inbred parent inherits a "proven" genotype and thus experiences a fitness advantage, which could explain the interaction. The positive and negative effects we found emphasize that continuing inbreeding can have important effects on individual fitness, even in populations that are already highly inbred.

Episodic positive selection in the evolution of avian toll-like receptor innate immunity genes.

Toll-like receptors (TLRs) are a family of conserved pattern-recognition molecules responsible for initiating innate and acquired immune responses. Because they play a key role in host defence, these genes have received increasing interest in the evolutionary and population genetics literature, as their variation represents a potential target of adaptive evolution. However, the role of pathogen-mediated selection (i.e. episodic positive selection) in the evolution of these genes remains poorly known and has not been examined outside of mammals. A recent increase in the number of bird species for which TLR sequences are available has enabled us to examine the selective processes that have influenced evolution of the 10 known avian TLR genes. Specifically, we tested for episodic positive selection to identify codons that experience purifying selection for the majority of their evolution, interspersed with bursts of positive selection that may occur only in restricted lineages. We included up to 23 species per gene (mean = 16.0) and observed that, although purifying selection was evident, an average of 4.5% of codons experienced episodic positive selection across all loci. For four genes in which sequence coverage traversed both the extracellular leucine-rich repeat region (LRR) and transmembrane/intracellular domains of the proteins, increased positive selection was observed at the extracellular domain, consistent with theoretical predictions. Our results provide evidence that episodic positive selection has played an important role in the evolution of most avian TLRs, consistent with the role of these loci in pathogen recognition and a mechanism of host-pathogen coevolution.

Severe inbreeding depression and no evidence of purging in an extremely inbred wild species--the Chatham Island black robin.

Although evidence of inbreeding depression in wild populations is well established, the impact of genetic purging in the wild remains controversial. The contrasting effects of inbreeding depression, fixation of deleterious alleles by genetic drift, and the purging of deleterious alleles via natural selection mean that predicting fitness outcomes in populations subjected to prolonged bottlenecks is not straightforward. We report results from a long-term pedigree study of arguably the world's most inbred wild species of bird: the Chatham Island black robin Petroica traversi, in which conditions were ideal for purging to occur. Contrary to expectations, black robins showed a strong, negative relationship between inbreeding and juvenile survival, yielding lethal equivalents (2B) of 6.85. We also determined that the negative relationship between inbreeding and survival did not appear to be mediated by levels of ancestral inbreeding and may be attributed in part to unpurged lethal recessives. Although the black robin demographic history provided ideal conditions for genetic purging, our results show no clear evidence of purging in the major life-history trait of juvenile survival. Our results also show no evidence of fixation of deleterious alleles in juvenile survival, but do confirm that continued high levels of contemporary inbreeding in a historically inbred population could lead to additional severe inbreeding depression.

Inbreeding influences within-brood heterozygosity-fitness correlations (HFCS) in an isolated passerine population.

Molecular estimates of inbreeding may be made using genetic markers such as microsatellites, however the interpretation of resulting heterozygosity-fitness correlations (HFCs) with respect to inbreeding depression is not straightforward. We investigated the relationship between pedigree-determined inbreeding coefficients (f) and HFCs in a closely monitored, reintroduced population of Stewart Island robins (Petroica australis rakiura) on Ulva Island, New Zealand. Using a full sibling design, we focused on differences in juvenile survival associated specifically with individual sibling variation in standardized multilocus heterozygosity (SH) when expected f was identical. We found that within broods, siblings with higher SH at microsatellite loci experienced a higher probability of juvenile survival. This effect, however, was detected primarily within broods that experienced inbreeding or when inbreeding had occurred in their pedigree histories (i.e., at the parents' level). Thus we show, for the first time in a wild population, that the strength of an HFC is partially dependent on the presence of inbreeding events in the recent pedigree history. Our results illustrate the importance of realized effects of inbreeding on genetic variation and fitness and the value of full-sibling designs for the study of HFCs in the context of small, inbred populations.

How does the 50/500 rule apply to MVPs?

The 50/500 rule has been used as a guiding principle in conservation for assessing minimum viable effective population size (N(e)). There is much confusion in the recent literature about how the 500 value should be applied to assess extinction risk and set priorities in conservation biology. Here, we argue that the confusion arises when the genetic basis for a short-term N(e) of 50 to avoid inbreeding depression is used to justify a long-term N(e) of 500 to maintain evolutionary potential. This confusion can result in misleading conclusions about how genetic arguments alone are sufficient to set minimum viable population (MVP) thresholds for assessing the extinction risk of threatened species, especially those that emphasize that MVPs should be in the thousands to maintain evolutionary potential.