Little Adaptive Potential in a Threatened Passerine Bird

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

A high-density linkage map reveals broad- and fine-scale sex differences in recombination in the hihi (stitchbird; Notiomystis cincta)

Recombination, the process of DNA exchange between homologous chromosomes during meiosis, plays a major role in genomic diversity and evolutionary change. Variation in recombination rate is widespread despite recombination often being essential for progression of meiosis. One such variation is heterochiasmy, where recombination rates differ between sexes. Heterochiasmy has been observed across broad taxonomic groups, yet it remains an evolutionary enigma. We used Lep-MAP3, a pedigree-based software that is efficient in handling large datasets, to generate linkage maps for the hihi or stitchbird (Notiomystis cincta), utilising information from >36 K SNPs and 36 families. We constructed 29 linkage maps, including for the previously unscaffolded Z chromosome. The hihi is an endangered passerine endemic to Aotearoa New Zealand that is sexually dimorphic and exhibits high levels of sexual conflict, including sperm competition. Patterns in recombination in the hihi are consistent with those in other birds, including higher recombination rates in micro-chromosomes. Heterochiasmy in the hihi is male-biased, in line with predictions of the Haldane-Huxley rule, with the male linkage map being 15% longer. Micro-chromosomes exhibit heterochiasmy to a greater extent, contrary to that reported in other birds. At the intra-chromosomal level, heterochiasmy is higher nearer to chromosome ends and in gene-rich regions. Regions of extreme heterochiasmy are enriched for genes implicated in cell structure. This study adds an important contribution in assessing evolutionary theories of heterochiasmy and provides a framework for future studies investigating fine-scale heterochiasmy.

Centrality to the metapopulation is more important for population genetic diversity than habitat area or fragmentation

Drift and gene flow affect genetic diversity. Given that the strength of genetic drift increases as population size decreases, management activities have focused on increasing population size through preserving habitats to preserve genetic diversity. Few studies have empirically evaluated the impacts of drift and gene flow on genetic diversity. Kryptolebias marmoratus, henceforth 'rivulus', is a small killifish restricted to fragmented New World mangrove forests with gene flow primarily associated with ocean currents. Rivulus form distinct populations across patches, making them a well-suited system to test the extent to which habitat area, fragmentation and connectivity are associated with genetic diversity. Using over 1000 individuals genotyped at 32 microsatellite loci, high-resolution landcover data and oceanographic simulations with graph theory, we demonstrate that centrality (connectivity) to the metapopulation is more strongly associated with genetic diversity than habitat area or fragmentation. By comparing models with and without centrality standardized by the source population's genetic diversity, our results suggest that metapopulation centrality is critical to genetic diversity regardless of the diversity of adjacent populations. While we find evidence that habitat area and fragmentation are related to genetic diversity, centrality is always a significant predictor with a larger effect than any measure of habitat configuration.

Genomic signatures of inbreeding depression for a threatened Aotearoa New Zealand passerine

For small and isolated populations, the increased chance of mating between related individuals can result in a substantial reduction in individual and population fitness. Despite the increasing availability of genomic data to measure inbreeding accurately across the genome, inbreeding depression studies for threatened species are still scarce due to the difficulty of measuring fitness in the wild. Here, we investigate inbreeding and inbreeding depression for the extensively monitored Tiritiri Mātangi island population of a threatened Aotearoa New Zealand passerine, the hihi (Notiomystis cincta). First, using a custom 45 k single nucleotide polymorphism (SNP) array, we explore genomic inbreeding patterns by inferring homozygous segments across the genome. Although all individuals have similar levels of ancient inbreeding, highly inbred individuals are affected by recent inbreeding, which can probably be explained by bottleneck effects such as habitat loss after European arrival and their translocation to the island in the 1990s. Second, we investigate genomic inbreeding effects on fitness, measured as lifetime reproductive success, and its three components, juvenile survival, adult annual survival and annual reproductive success, in 363 hihi. We find that global inbreeding significantly affects juvenile survival but none of the remaining fitness traits. Finally, we employ a genome-wide association approach to test the locus-specific effects of inbreeding on fitness, and identify 13 SNPs significantly associated with lifetime reproductive success. Our findings suggest that inbreeding depression does impact hihi, but at different genomic scales for different traits, and that purging has therefore failed to remove all variants with deleterious effects from this population of conservation concern.

Madagascar’s extraordinary biodiversity: Threats and opportunities

Madagascar’s unique biota is heavily affected by human activity and is under intense threat. Here, we review the current state of knowledge on the conservation status of Madagascar’s terrestrial and freshwater biodiversity by presenting data and analyses on documented and predicted species-level conservation statuses, the most prevalent and relevant threats, ex situ collections and programs, and the coverage and comprehensiveness of protected areas. The existing terrestrial protected area network in Madagascar covers 10.4% of its land area and includes at least part of the range of the majority of described native species of vertebrates with known distributions (97.1% of freshwater fishes, amphibians, reptiles, birds, and mammals combined) and plants (67.7%). The overall figures are higher for threatened species (97.7% of threatened vertebrates and 79.6% of threatened plants occurring within at least one protected area). International Union for Conservation of Nature (IUCN) Red List assessments and Bayesian neural network analyses for plants identify overexploitation of biological resources and unsustainable agriculture as the most prominent threats to biodiversity. We highlight five opportunities for action at multiple levels to ensure that conservation and ecological restoration objectives, programs, and activities take account of complex underlying and interacting factors and produce tangible benefits for the biodiversity and people of Madagascar.

Telomere length is heritable and genetically correlated with lifespan in a wild bird

Telomeres are protective caps at the end of eukaryotic chromosomes that shorten with age and in response to stressful or resource-demanding conditions. Their length predicts individual health and lifespan across a wide range of animals, but whether the observed positive association between telomere length and lifespan is environmentally induced, or set at conception due to a shared genetic basis, has not been tested in wild animals. We applied quantitative genetic "animal models" to longitudinal telomere measurements collected over a 10-year period from individuals of a wild seabird (common tern; Sterna hirundo) with known pedigree. We found no variation in telomere shortening with age among individuals at the phenotypic and genetic level, and only a small permanent environmental effect on adult telomere length. Instead, we found telomere length to be highly heritable and strongly positively genetically correlated with lifespan. Such heritable differences between individuals that are set at conception may present a hitherto underappreciated component of variation in somatic state.

Variation in shape and consistency of selection between populations of the threatened Hihi (Notiomystis cincta)

The shape and intensity of natural selection can vary between years, potentially resulting in a chronic reduction of fitness as individuals need to track a continually changing optimum of fitness (i.e., a "lag load"). In endangered species, often characterized by small population size, the lack of genetic diversity is expected to limit the response to this constant need to adjust to fluctuating selection, increasing the fitness burden and thus the risk of extinction. Here, we use long-term monitoring data to assess whether the type of selection for a key fitness trait (i.e., lay date) differs between two reintroduced populations of a threatened passerine bird, the hihi (Notiomystis cincta). We apply recent statistical developments to test for the presence or absence of fluctuation in selection in both the Tiritiri Mātangi Island and the Kārori sanctuary populations. Our results support the presence of stabilizing selection in Tiritiri Mātangi with a potential moving optimum for lay date. In Kārori our results favour a regime of directional selection. Although the shape of selection may differ, for both populations an earlier lay date generally increases fitness in both environments. Further, the moving optimum models of lay date on Tiritiri Mātangi, suggesting that selection varies between years, imply a substantial lag load in addition to the fitness burden caused by the population laying too late. Our results highlight the importance of characterizing the form and temporal variation of selection for each population to predict the effects of environmental change and to inform management.

The quantitative genetics of fitness in a wild seabird

Additive genetic variance in fitness is a prerequisite for adaptive evolution, as a trait must be genetically correlated with fitness to evolve. Despite its relevance, additive genetic variance in fitness has not often been estimated in nature. Here, we investigate additive genetic variance in lifetime and annual fitness components in common terns (Sterna hirundo). Using 28 years of data comprising approximately 6000 pedigreed individuals, we find that additive genetic variances in the zero-inflated and Poisson components of lifetime fitness were effectively zero but estimated with high uncertainty. Similarly, additive genetic variances in adult annual reproductive success and survival did not differ from zero but were again associated with high uncertainty. Simulations suggested that we would be able to detect additive genetic variances as low as 0.05 for the zero-inflated component of fitness but not for the Poisson component, for which adequate statistical power would require approximately two more decades (four tern generations) of data collection. As such, our study suggests heritable variance in common tern fitness to be rather low if not zero, shows how studying the quantitative genetics of fitness in natural populations remains challenging, and highlights the importance of maintaining long-term individual-based studies of natural populations.

Genetic variance in fitness indicates rapid contemporary adaptive evolution in wild animals

The rate of adaptive evolution, the contribution of selection to genetic changes that increase mean fitness, is determined by the additive genetic variance in individual relative fitness. To date, there are few robust estimates of this parameter for natural populations, and it is therefore unclear whether adaptive evolution can play a meaningful role in short-term population dynamics. We developed and applied quantitative genetic methods to long-term datasets from 19 wild bird and mammal populations and found that, while estimates vary between populations, additive genetic variance in relative fitness is often substantial and, on average, twice that of previous estimates. We show that these rates of contemporary adaptive evolution can affect population dynamics and hence that natural selection has the potential to partly mitigate effects of current environmental change.

Who are you? A framework to identify and report genetic sample mix-ups

Sample mix-ups occur when samples have accidentally been duplicated, mislabelled or swapped. When samples are subsequently genotyped or sequenced, this can lead to individual IDs being incorrectly linked to genetic data, resulting in incorrect or biased research results, or reduced power to detect true biological patterns. We surveyed the community and found that almost 80% of responding researchers have encountered sample mix-ups. However, many recent studies in the field of molecular ecology do not appear to systematically report individual assignment checks as part of their publications. Although checks may be done, lack of consistent reporting means that it is difficult to assess whether sample mix-ups have occurred or been detected. Here, we present an easy-to-follow sample verification framework that can utilise existing metadata, including species, population structure, sex and pedigree information. We demonstrate its application to a dataset representing individuals of a threatened Aotearoa New Zealand bird species, the hihi, genotyped on a 50K SNP array. We detected numerous incorrect genotype-ID associations when comparing observed and genetic sex or comparing to relationships in a verified microsatellite pedigree. The framework proposed here helped to confirm 488 individuals (39%), correct another 20 bird-genotype links, and detect hundreds of incorrect sample IDs, emphasizing the value of routinely checking genetic and genomic datasets for their accuracy. We therefore promote the implementation and reporting of this simple yet effective sample verification framework as a standardized quality control step for studies in the field of molecular ecology.

Population structure, genomic diversity and demographic history of Komodo dragons inferred from whole-genome sequencing

Population and conservation genetics studies have greatly benefited from the development of new techniques and bioinformatic tools associated with next-generation sequencing. Analysis of extensive data sets from whole-genome sequencing of even a few individuals allows the detection of patterns of fine-scale population structure and detailed reconstruction of demographic dynamics through time. In this study, we investigated the population structure, genomic diversity and demographic history of the Komodo dragon (Varanus komodoensis), the world's largest lizard, by sequencing the whole genomes of 24 individuals from the five main Indonesian islands comprising the entire range of the species. Three main genomic groups were observed. The populations of the Island of Komodo and the northern coast of Flores, in particular, were identified as two distinct conservation units. Degrees of genomic divergence among island populations were interpreted as a result of changes in sea level affecting connectivity across islands. Demographic inference suggested that Komodo dragons probably experienced a relatively steep population decline over the last million years, reaching a relatively stable N_e during the Saalian glacial cycle (400-150 thousand years ago) followed by a rapid N_e decrease. Genomic diversity of Komodo dragons was similar to that found in endangered or already extinct reptile species. Overall, this study provides an example of how whole-genome analysis of a few individuals per population can help define population structure and intraspecific demographic dynamics. This is particularly important when applying population genomics data to conservation of rare or elusive endangered species.

The crucial role of genome-wide genetic variation in conservation

The unprecedented rate of extinction calls for efficient use of genetics to help conserve biodiversity. Several recent genomic and simulation-based studies have argued that the field of conservation biology has placed too much focus on conserving genome-wide genetic variation, and that the field should instead focus on managing the subset of functional genetic variation that is thought to affect fitness. Here, we critically evaluate the feasibility and likely benefits of this approach in conservation. We find that population genetics theory and empirical results show that conserving genome-wide genetic variation is generally the best approach to prevent inbreeding depression and loss of adaptive potential from driving populations toward extinction. Focusing conservation efforts on presumably functional genetic variation will only be feasible occasionally, often misleading, and counterproductive when prioritized over genome-wide genetic variation. Given the increasing rate of habitat loss and other environmental changes, failure to recognize the detrimental effects of lost genome-wide genetic variation on long-term population viability will only worsen the biodiversity crisis.

New developments in the field of genomic technologies and their relevance to conservation management

Recent technological advances in the field of genomics offer conservation managers and practitioners new tools to explore for conservation applications. Many of these tools are well developed and used by other life science fields, while others are still in development. Considering these technological possibilities, choosing the right tool(s) from the toolbox is crucial and can pose a challenging task. With this in mind, we strive to inspire, inform and illuminate managers and practitioners on how conservation efforts can benefit from the current genomic and biotechnological revolution. With inspirational case studies we show how new technologies can help resolve some of the main conservation challenges, while also informing how implementable the different technologies are. We here focus specifically on small population management, highlight the potential for genetic rescue, and discuss the opportunities in the field of gene editing to help with adaptation to changing environments. In addition, we delineate potential applications of gene drives for controlling invasive species. We illuminate that the genomic toolbox offers added benefit to conservation efforts, but also comes with limitations for the use of these novel emerging techniques.

Phenotypic plasticity drives phenological changes in a Mediterranean blue tit population

Earlier phenology induced by climate change, such as the passerines' breeding time, is observed in many natural populations. Understanding the nature of such changes is key to predict the responses of wild populations to climate change. Genetic changes have been rarely investigated for laying date, though it has been shown to be heritable and under directional selection, suggesting that the trait could evolve. In a Corsican blue tit population, the birds' laying date has significantly advanced over 40 years, and we here determine whether this response is of plastic or evolutionary origin, by comparing the predictions of the breeder's and the Robertson-Price (STS) equations, to the observed genetic changes. We compare the results obtained for two fitness proxies (fledgling and recruitment success), using models accounting for their zero inflation. Because the trait appears heritable and under directional selection, the breeder's equation predicts that genetic changes could drive a significant part of the phenological change observed. We, however, found that fitness proxies and laying date are not genetically correlated. The STS, therefore, predicts no evolution of the breeding time, predicting correctly the absence of trend in breeding values. Our results also emphasize that when investigating selection on a plastic trait under fluctuating selection, part of the fitness-trait phenotypic covariance can be due to within individual covariance. In the case of repeated measurements, splitting within and between individual covariance can shift our perspective on the actual intensity of selection over multiple selection episodes, shedding light on the potential for the trait to evolve.

Genomic signatures of inbreeding in a critically endangered parrot, the kākāpō

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

The relevance of pedigrees in the conservation genomics era

Over the past 50 years conservation genetics has developed a substantive toolbox to inform species management. One of the most long‐standing tools available to manage genetics—the pedigree—has been widely used to characterize diversity and maximize evolutionary potential in threatened populations. Now, with the ability to use high throughput sequencing to estimate relatedness, inbreeding, and genome‐wide functional diversity, some have asked whether it is warranted for conservation biologists to continue collecting and collating pedigrees for species management. In this perspective, we argue that pedigrees remain a relevant tool, and when combined with genomic data, create an invaluable resource for conservation genomic management. Genomic data can address pedigree pitfalls (e.g., founder relatedness, missing data, uncertainty), and in return robust pedigrees allow for more nuanced research design, including well‐informed sampling strategies and quantitative analyses (e.g., heritability, linkage) to better inform genomic inquiry. We further contend that building and maintaining pedigrees provides an opportunity to strengthen trusted relationships among conservation researchers, practitioners, Indigenous Peoples, and Local Communities.

The design and application of a 50 K SNP chip for a threatened Aotearoa New Zealand passerine, the hihi

Next-generation sequencing has transformed the fields of ecological and evolutionary genetics by allowing for cost-effective identification of genome-wide variation. Single nucleotide polymorphism (SNP) arrays, or "SNP chips", enable very large numbers of individuals to be consistently genotyped at a selected set of these identified markers, and also offer the advantage of being able to analyse samples of variable DNA quality. We used reduced representation restriction-aided digest sequencing (RAD-seq) of 31 birds of the threatened hihi (Notiomystis cincta; stitchbird) and low-coverage whole genome sequencing (WGS) of 10 of these birds to develop an Affymetrix 50 K SNP chip. We overcame the limitations of having no hihi reference genome and a low quantity of sequence data by separate and pooled de novo assembly of each of the 10 WGS birds. Reads from all individuals were mapped back to these de novo assemblies to identify SNPs. A subset of RAD-seq and WGS SNPs were selected for inclusion on the chip, prioritising SNPs with the highest quality scores whose flanking sequence uniquely aligned to the zebra finch (Taeniopygia guttata) genome. Of the 58,466 SNPs manufactured on the chip, 72% passed filtering metrics and were polymorphic. By genotyping 1,536 hihi on the array, we found that SNPs detected in multiple assemblies were more likely to successfully genotype, representing a cost-effective approach to identify SNPs for genotyping. Here, we demonstrate the utility of the SNP chip by describing the high rates of linkage disequilibrium in the hihi genome, reflecting the history of population bottlenecks in the species.

Repeatability and heritability of social reaction norms in a wild agamid lizard

In the evolutionary transition from solitary to group living, it should be adaptive for animals to respond to the environment and choose when to socialize to reduce conflict and maximize access to resources. Due to the associated proximate mechanisms (e.g. neural network, endocrine system), it is likely that this behavior varies between individuals according to genetic and non-genetic factors. We used long-term behavioral and genetic data from a population of eastern water dragons (Intellagama lesueurii) to explore variation in plasticity of social behavior, in response to sex ratio and density. To do so, we modeled individual variation in social reaction norms, which describe individuals' mean behavior and behavioral responses to changes in their environment, and partitioned variance into genetic and non-genetic components. We found that reaction norms were repeatable over multiple years, suggesting that individuals consistently differed in their behavioral responses to changes in the social environment. Despite high repeatability of reaction norm components, trait heritability was below our limit of detection based on power analyses (h² < 0.12), leading to very little power to detect heritability of plasticity. This was in contrast to a relatively greater amount of variance associated with environmental effects. This could suggest that mechanisms such as social learning and frequency-dependence may shape variance in reaction norms, which will be testable as the dataset grows.

Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi

Inbreeding can lead to a loss of heterozygosity in a population and when combined with genetic drift may reduce the adaptive potential of a species. However, there is uncertainty about whether resequencing data can provide accurate and consistent inbreeding estimates. Here, we performed an in-depth inbreeding analysis for hihi (Notiomystis cincta), an endemic and nationally vulnerable passerine bird of Aotearoa New Zealand. We first focused on subsampling variants from a reference genome male, and found that low-density data sets tend to miss runs of homozygosity (ROH) in some places and overestimate ROH length in others, resulting in contrasting homozygosity landscapes. Low-coverage resequencing and 50 K SNP array densities can yield comparable inbreeding results to high-coverage resequencing approaches, but the results for all data sets are highly dependent on the software settings employed. Second, we extended our analysis to 10 hihi where low-coverage whole genome resequencing, RAD-seq and SNP array genotypes are available. We inferred ROH and individual inbreeding to evaluate the relative effects of sequencing depth versus SNP density on estimating inbreeding coefficients and found that high rates of missingness downwardly bias both the number and length of ROH. In summary, when using genomic data to evaluate inbreeding, studies must consider that ROH estimates are heavily dependent on analysis parameters, data set density and individual sequencing depth.

Genetic load has potential in large populations but is realized in small inbred populations

Populations with higher genetic diversity and larger effective sizes have greater evolutionary capacity (i.e., adaptive potential) to respond to ecological stressors. We are interested in how the variation captured in protein-coding genes fluctuates relative to overall genomic diversity and whether smaller populations suffer greater costs due to their genetic load of deleterious mutations compared with larger populations. We analyzed individual whole-genome sequences (N = 74) from three different populations of Montezuma quail (Cyrtonyx montezumae), a small ground-dwelling bird that is sustainably harvested in some portions of its range but is of conservation concern elsewhere. Our historical demographic results indicate that Montezuma quail populations in the United States exhibit low levels of genomic diversity due in large part to long-term declines in effective population sizes over nearly a million years. The smaller and more isolated Texas population is significantly more inbred than the large Arizona and the intermediate-sized New Mexico populations we surveyed. The Texas gene pool has a significantly smaller proportion of strongly deleterious variants segregating in the population compared with the larger Arizona gene pool. Our results demonstrate that even in small populations, highly deleterious mutations are effectively purged and/or lost due to drift. However, we find that in small populations the realized genetic load is elevated because of inbreeding coupled with a higher frequency of slightly deleterious mutations that are manifested in homozygotes. Overall, our study illustrates how population genomics can be used to proactively assess both neutral and functional aspects of contemporary genetic diversity in a conservation framework while simultaneously considering deeper demographic histories.

Adaptive Radiation Genomics of Two Ecologically Divergent Hawai'ian Honeycreepers: The 'akiapōlā'au and the Hawai'i 'amakihi

The Hawai'ian honeycreepers (drepanids) are a classic example of adaptive radiation: they adapted to a variety of novel dietary niches, evolving a wide range of bill morphologies. Here we investigated genomic diversity, demographic history, and genes involved in bill morphology phenotypes in 2 honeycreepers: the 'akiapōlā'au (Hemignathus wilsoni) and the Hawai'i 'amakihi (Chlorodrepanis virens). The 'akiapōlā'au is an endangered island endemic, filling the "woodpecker" niche by using a unique bill morphology, while the Hawai'i 'amakihi is a dietary generalist common on the islands of Hawai'i and Maui. We de novo sequenced the 'akiapōlā'au genome and compared it to the previously sequenced 'amakihi genome. The 'akiapōlā'au is far less heterozygous and has a smaller effective population size than the 'amakihi, which matches expectations due to its smaller census population and restricted ecological niche. Our investigation revealed genomic islands of divergence, which may be involved in the honeycreeper radiation. Within these islands of divergence, we identified candidate genes (including DLK1, FOXB1, KIF6, MAML3, PHF20, RBP1, and TIMM17A) that may play a role in honeycreeper adaptations. The gene DLK1, previously shown to influence Darwin's finch bill size, may be related to honeycreeper bill morphology evolution, while the functions of the other candidates remain unknown.

Genome insights give cause for optimism in the ongoing battle to save the vaquita

Whilst the global threat to biodiversity is accelerating, recent advances in sequencing and assembly technologies, along with decreasing costs, are ushering in a golden age for biodiversity genomics. In a From the Cover article in this issue of Molecular Ecology Resources, Morin et al. report an exceptionally high-quality draft genome assembly of the vaquita (Phocoena sinus), a porpoise first described in 1958 and currently numbering fewer than 20 individuals in the wild. On the brink of extinction and described as critically endangered by the International Union for the Conservation of Nature (IUCN), a precipitous population decline since the 1990 s has been due to bycatch in gillnets in the Upper Gulf of California.

Population genomics for wildlife conservation and management

Biodiversity is under threat worldwide. Over the past decade, the field of population genomics has developed across nonmodel organisms, and the results of this research have begun to be applied in conservation and management of wildlife species. Genomics tools can provide precise estimates of basic features of wildlife populations, such as effective population size, inbreeding, demographic history and population structure, that are critical for conservation efforts. Moreover, population genomics studies can identify particular genetic loci and variants responsible for inbreeding depression or adaptation to changing environments, allowing for conservation efforts to estimate the capacity of populations to evolve and adapt in response to environmental change and to manage for adaptive variation. While connections from basic research to applied wildlife conservation have been slow to develop, these connections are increasingly strengthening. Here we review the primary areas in which population genomics approaches can be applied to wildlife conservation and management, highlight examples of how they have been used, and provide recommendations for building on the progress that has been made in this field.

Polygenic basis for adaptive morphological variation in a threatened Aotearoa | New Zealand bird, the hihi (Notiomystis cincta)

To predict if a threatened species can adapt to changing selective pressures, it is crucial to understand the genetic basis of adaptive traits, especially in species historically affected by severe bottlenecks. We estimated the heritability of three hihi (Notiomystis cincta) morphological traits known to be under selection (nestling tarsus length, body mass and head-bill length) using 523 individuals and 39 699 single nucleotide polymorphisms (SNPs) from a 50 K Affymetrix SNP chip. We then examined the genetic architecture of the traits via chromosome partitioning analyses and genome-wide association scans (GWAS). Heritabilities estimated using pedigree relatedness or genomic relatedness were low. For tarsus length, the proportion of genetic variance explained by each chromosome was positively correlated with its size, and more than one chromosome explained significant variation for body mass and head-bill length. Finally, GWAS analyses suggested many loci of small effect contributing to trait variation for all three traits, although one locus (an SNP within an intron of the transcription factor HEY2) was tentatively associated with tarsus length. Our findings suggest a polygenic nature for the morphological traits, with many small effect size loci contributing to the majority of the variation, similar to results from many other wild populations. However, the small effective population size, polygenic architecture and already low heritabilities suggest that both the total response and rate of response to selection are likely to be limited in hihi.

High MHC gene copy number maintains diversity despite homozygosity in a Critically Endangered single-island endemic bird, but no evidence of MHC-based mate choice

Small population sizes can, over time, put species at risk due to the loss of genetic variation and the deleterious effects of inbreeding. Losing diversity in the major histocompatibility complex (MHC) could be particularly harmful, given its key role in the immune system. Here, we assess MHC class I (MHC-I) diversity and its effects on mate choice and survival in the Critically Endangered Raso lark Alauda razae, a species restricted to the 7 km² islet of Raso, Cape Verde, since ~1460, whose population size has dropped as low as 20 pairs. Exhaustively genotyping 122 individuals, we find no effect of MHC-I genotype/diversity on mate choice or survival. However, we demonstrate that MHC-I diversity has been maintained through extreme bottlenecks by retention of a high number of gene copies (at least 14), aided by cosegregation of multiple haplotypes comprising 2-8 linked MHC-I loci. Within-locus homozygosity is high, contributing to low population-wide diversity. Conversely, each individual had comparably many alleles, 6-16 (average 11), and the large and divergent haplotypes occur at high frequency in the population, resulting in high within-individual MHC-I diversity. This functional immune gene diversity will be of critical importance for this highly threatened species' adaptive potential.

A comparison of pedigree, genetic and genomic estimates of relatedness for informing pairing decisions in two critically endangered birds: Implications for conservation breeding programmes worldwide

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimize relatedness between individuals in an effort to avoid inbreeding, maximize diversity and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome-wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites and SNPs from whole-genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange-fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent-offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree- and SNP-based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

Genetic diversity, demographic history and neo-sex chromosomes in the Critically Endangered Raso lark

Small effective population sizes could expose island species to inbreeding and loss of genetic variation. Here, we investigate factors shaping genetic diversity in the Raso lark, which has been restricted to a single islet for approximately 500 years, with a population size of a few hundred. We assembled a reference genome for the related Eurasian skylark and then assessed diversity and demographic history using RAD-seq data (75 samples from Raso larks and two related mainland species). We first identify broad tracts of suppressed recombination in females, indicating enlarged neo-sex chromosomes. We then show that genetic diversity across autosomes in the Raso lark is lower than in its mainland relatives, but inconsistent with long-term persistence at its current population size. Finally, we find that genetic signatures of the recent population contraction are overshadowed by an ancient expansion and persistence of a very large population until the human settlement of Cape Verde. Our findings show how genome-wide approaches to study endangered species can help avoid confounding effects of genome architecture on diversity estimates, and how present-day diversity can be shaped by ancient demographic events.

Evolutionary genomics can improve prediction of species' responses to climate change

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco-evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large-scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco-evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions.

Conservation: Bye-Bye to the Hihi?

Theoretically, small populations should be less capable of adapting to environmental change, yet empirical evidence is mixed. A new study on a rare New Zealand bird, the Hihi, uses genomics and long-term demographic data to reveal low adaptive potential.

Additive genetic variance for lifetime fitness and the capacity for adaptation in an annual plant

The immediate capacity for adaptation under current environmental conditions is directly proportional to the additive genetic variance for fitness, V_A (W). Mean absolute fitness, W ¯ , is predicted to change at the rate V A ( W ) W ¯ , according to Fisher's Fundamental Theorem of Natural Selection. Despite ample research evaluating degree of local adaptation, direct assessment of V_A (W) and the capacity for ongoing adaptation is exceedingly rare. We estimated V_A (W) and W ¯ in three pedigreed populations of annual Chamaecrista fasciculata, over three years in the wild. Contrasting with common expectations, we found significant V_A (W) in all populations and years, predicting increased mean fitness in subsequent generations (0.83 to 6.12 seeds per individual). Further, we detected two cases predicting "evolutionary rescue," where selection on standing V_A (W) was expected to increase fitness of declining populations ( W ¯ < 1.0) to levels consistent with population sustainability and growth. Within populations, inter-annual differences in genetic expression of fitness were striking. Significant genotype-by-year interactions reflected modest correlations between breeding values across years, indicating temporally variable selection at the genotypic level that could contribute to maintaining V_A (W). By directly estimating V_A (W) and total lifetime W ¯ , our study presents an experimental approach for studies of adaptive capacity in the wild.