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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

Variation at innate immunity Toll-like receptor genes in a bottlenecked population of a New Zealand robin.

Toll-like receptors (TLRs) are an ancient family of genes encoding transmembrane proteins that bind pathogen-specific molecules and initiate both innate and adaptive aspects of the immune response. Our goal was to determine whether these genes show sufficient genetic diversity in a bottlenecked population to be a useful addition or alternative to the more commonly employed major histocompatibility complex (MHC) genotyping in a conservation genetics context. We amplified all known avian TLR genes in a severely bottlenecked population of New Zealand's Stewart Island robin (Petroica australis rakiura), for which reduced microsatellite diversity was previously observed. We genotyped 17-24 birds from a reintroduced island population (including the 12 founders) for nine genes, seven of which were polymorphic. We observed a total of 24 single-nucleotide polymorphisms overall, 15 of which were non-synonymous, representing up to five amino-acid variants at a locus. One locus (TLR1LB) showed evidence of past directional selection. Results also confirmed a passerine duplication of TLR7. The levels of TLR diversity that we observe are sufficient to justify their further use in addressing conservation genetic questions, even in bottlenecked 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.

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.

Dye shift: a neglected source of genotyping error in molecular ecology.

Molecular ecologists must be vigilant in detecting and accounting for genotyping error, yet potential errors stemming from dye-induced mobility shift (dye shift) may be frequently neglected and largely unknown to researchers who employ 3-primer systems with automated genotyping. When left uncorrected, dye shift can lead to mis-scoring alleles and even to falsely calling new alleles if different dyes are used to genotype the same locus in subsequent reactions. When we used four different fluorophore labels from a standard dye set to genotype the same set of loci, differences in the resulting size estimates for a single allele ranged from 2.07 bp to 3.68 bp. The strongest effects were associated with the fluorophore PET, and relative degree of dye shift was inversely related to locus size. We found little evidence in the literature that dye shift is regularly accounted for in 3-primer studies, despite knowledge of this phenomenon existing for over a decade. However, we did find some references to erroneous standard correction factors for the same set of dyes that we tested. We thus reiterate the need for strict quality control when attempting to reduce possible sources of genotyping error, and in cases where different dyes are applied to a single locus, perhaps mistakenly, we strongly discourage researchers from assuming generic correction patterns.