Postglacial expansion and not human influence best explains the population structure in the endangered kea (Nestor notabilis)

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

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

Exploring the endangerment mechanisms of Hipposideros pomona based on molecular phylogeographic methods

The endangerment mechanisms of various species are a focus of studies on biodiversity and conservation biology. Hipposideros pomona is an endangered species, but the reasons behind its endangerment remain unclear. We investigated the endangerment mechanisms of H. pomona using mitochondrial DNA, nuclear DNA, and microsatellite loci markers. The results showed that the nucleotide diversity of mitochondria DNA and heterozygosity of microsatellite markers were high (π = 0.04615, H O = 0.7115), whereas the nucleotide diversity of the nuclear genes was low (THY: π = 0.00508, SORBS2: π = 0.00677, ACOX2: π = 0.00462, COPS7A: π = 0.00679). The phylogenetic tree and median-joining network based on mitochondrial DNA sequences clustered the species into three clades, namely North Vietnam-Fujian, Myanmar-West Yunnan, and Laos-Hainan clades. However, joint analysis of nuclear genes did not exhibit clustering. Analysis of molecular variance revealed a strong population genetic structure; IMa2 analysis did not reveal significant gene flow between all groups (p > .05), and isolation-by-distance analysis revealed a significant positive correlation between genetic and geographic distances (p < .05). The mismatch distribution analysis, neutral test, and Bayesian skyline plots revealed that the H. pomona population were relatively stable and exhibited a contraction trend. The results implied that H. pomona exhibits female philopatry and male-biased dispersal. The Hengduan Mountains could have acted as a geographical barrier for gene flow between the North Vietnam-Fujian clade and the Myanmar-West Yunnan clade, whereas the Qiongzhou Strait may have limited interaction between the Hainan populations and other clades. The warm climate during the second interglacial Quaternary period (c. 0.33 Mya) could have been responsible for species differentiation, whereas the cold climate during the late Quaternary last glacial maximum (c. 10 ka BP) might have caused the overall contraction of species. The lack of significant gene flow in nuclear microsatellite loci markers among the different populations investigated reflects recent habitat fragmentation due to anthropogenic activities; thus, on-site conservation of the species and restoration of gene flow corridors among populations need immediate implementation.

Interpreting past trophic ecology of a threatened alpine parrot, kea Nestor notabilis, from museum specimens

When ecosystems are under severe pressure or environments change, trophic position and intraspecific niche width may decrease or narrow, signalling that conservation action is required. In New Zealand, alpine and subalpine ecosystems have been extensively modified through farming since 19th-century European settlement, with consequences for indigenous species such as the kea Nestor notabilis. We investigated feather stable isotope values in the kea and predicted a lower trophic position in modern kea populations, to reflect reduced lowland habitat and a mixed diet with more plant material. We predicted that size and sex would influence trophic values in this sexually dimorphic species, with larger birds more likely to have a high protein diet. We examined potential dietary changes in 68 museum collected kea from 1880s to 2000s, first recording accession details including provenance and sex and measuring culmen length. We used bulk carbon and nitrogen stable isotope analyses (BSIAs) of feathers and a further feather subset using compound-specific stable isotope analyses of amino acids (CSIA-AA) to obtain isotopic values and estimate trophic position. BSIA showed δ¹⁵ N values in kea feathers declined through time and could indicate that early century kea were highly omnivorous, with δ¹⁵ N values on average higher than in modern kea. Variance in δ¹⁵ N values was greater after 1950, driven by a few individuals. Few differences between males and females were evident, although females in the south region had lower δ¹⁵ N values. There was a tendency for large male birds to have higher trophic values, perhaps reflecting dominant male bird behaviour noted in historical records. Nonetheless, CSIA-AA performed on a subset of the data suggested that variation in BSIA is likely due to baseline changes rather than relative trophic position which may be more homogenous than these data indicate. Although there was more variability in modern kea, we suggest caution in interpretation. Stable isotope data, particularly CSIA-AA, from museum specimens can reveal potential change in ecological networks as well as sexually dimorphic feeding patterns within species. The data can reveal temporal and regional variation in species trophic position and changes in ecosystem integrity to inform conservation decision-making.

Thirty years of ancient DNA and the faunal biogeography of Aotearoa New Zealand: lessons and future directions

Thirty years ago, DNA sequences were obtained from an extinct Aotearoa New Zealand animal for the first time. Since then, ancient DNA research has provided many - often unexpected - insights into the origins of New Zealand's terrestrial and marine vertebrate fauna. Because recent human activities in New Zealand have caused the decline or extinction of many endemic plant, bird, reptile, and marine mammal species, ancient DNA has been instrumental in reconstructing their identities and origins. However, most ancient DNA studies focusing on New Zealand species have been restricted to vertebrates, with small sample sizes, and/or relatively few genetic markers. This has limited their power to infer fine-scale biogeographic patterns, including (pre)historic distributions and range-shifts driven by past climate and environmental change. Recently, 'next-generation' methodological and technological advances have broadened the range of hypotheses that can feasibly be tested with ancient DNA. These advances represent an exciting opportunity for further exploring New Zealand biogeography using ancient DNA, but their promise has not yet been fully realised. In this review, we summarise the last 30 years of ancient DNA research into New Zealand faunal biogeography and highlight key objectives, challenges, and possibilities for the next 30 years and beyond.

Identifying loci under selection via explicit demographic models

Adaptive genetic variation is a function of both selective and neutral forces. To accurately identify adaptive loci, it is thus critical to account for demographic history. Theory suggests that signatures of selection can be inferred using the coalescent, following the premise that genealogies of selected loci deviate from neutral expectations. Here, we build on this theory to develop an analytical framework to identify loci under selection via explicit demographic models (LSD). Under this framework, signatures of selection are inferred through deviations in demographic parameters, rather than through summary statistics directly, and demographic history is accounted for explicitly. Leveraging the property of demographic models to incorporate directionality, we show that LSD can provide information on the environment in which selection acts on a population. This can prove useful in elucidating the selective processes underlying local adaptation, by characterizing genetic trade-offs and extending the concepts of antagonistic pleiotropy and conditional neutrality from ecological theory to practical application in genomic data. We implement LSD via approximate Bayesian computation and demonstrate, via simulations, that LSD (a) has high power to identify selected loci across a large range of demographic-selection regimes, (b) outperforms commonly applied genome-scan methods under complex demographies and (c) accurately infers the directionality of selection for identified candidates. Using the same simulations, we further characterize the behaviour of isolation-with-migration models conducive to the study of local adaptation under regimes of selection. Finally, we demonstrate an application of LSD by detecting loci and characterizing genetic trade-offs underlying flower colour in Antirrhinum majus.

MUCOPOLYSACCHARIDOSIS II (MPS II) IN A FREE-LIVING KAKA (NESTOR MERIDIONALIS) IN NEW ZEALAND

A lysosomal storage disease, identified as a mucopolysaccharidosis (MPS), was diagnosed in a free-living Kaka (Nestor meridionalis), an endemic New Zealand parrot, which exhibited weakness, incoordination, and seizures. Histopathology showed typical colloid-like cytoplasmic inclusions in Purkinje cells and many other neurons throughout the brain. Electron microscopy revealed that storage bodies contained a variety of linear, curved, or circular membranous profiles and electron-dense bodies. Because the bird came from a small isolated population of Kaka in the northern South Island, a genetic cause was deemed likely. Tandem mass spectrometry revealed increased levels of heparan sulfate-derived disaccharides in the brain and liver compared with tissues from controls. Enzymatic assays documented low levels of iduronate-2-sulfatase activity, which causes a lysosomal storage disorder called MPS type II or Hunter syndrome. A captive breeding program is currently in progress, and the possibility of detecting carriers of this disorder warrants further investigation.

Genetic divergence between isolated populations of the North Island New Zealand Rifleman (Acanthisitta chloris granti) implicates ancient biogeographic impacts rather than recent habitat fragmentation

This research investigates the extent and causal mechanisms of genetic population divergence in a poorly flighted passerine, the North Island Rifleman or Titipounamu (Acanthisitta chloris granti). While this species has a historically widespread distribution, anthropogenic forest clearance has resulted in a highly fragmented current distribution. We conducted analyses of mitochondrial DNA (COI and Control Region) and 12 nuclear DNA microsatellites to test for population divergence and estimate times of divergence. diyabc and biogeobears were then used to assess likely past dispersal scenarios based on both mtDNA and nDNA. The results reveal several significantly divergent lineages across the North Island of New Zealand and indicate that some populations have been isolated for extensive periods of time (0.7-4.9 mya). Modeling indicated a dynamic history of population connectivity, with a drastic restriction in gene flow between three geographic regions, followed by a more recent re-establishment of connectivity. Our analyses indicate the dynamic influence of key geological and climatological events on the distribution of genetic diversity in this species, including support for the genetic impact of old biogeographic boundaries such as the Taupo Line and Cockayne's Line, rather than recent anthropogenic habitat fragmentation. These findings present a rare example of an avian species with a genetic history more like that of flightless taxa and so provide new general insights into vicariant processes affecting populations of passerines with limited dispersal.

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

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

Reduced representation sequencing detects only subtle regional structure in a heavily exploited and rapidly recolonizing marine mammal species

Next-generation reduced representation sequencing (RRS) approaches show great potential for resolving the structure of wild populations. However, the population structure of species that have shown rapid demographic recovery following severe population bottlenecks may still prove difficult to resolve due to high gene flow between subpopulations. Here, we tested the effectiveness of the RRS method Genotyping-By-Sequencing (GBS) for describing the population structure of the New Zealand fur seal (NZFS, Arctocephalus forsteri), a species that was heavily exploited by the 19th century commercial sealing industry and has since rapidly recolonized most of its former range from a few isolated colonies. Using 26,026 neutral single nucleotide polymorphisms (SNPs), we assessed genetic variation within and between NZFS colonies. We identified low levels of population differentiation across the species range (<1% of variation explained by regional differences) suggesting a state of near panmixia. Nonetheless, we observed subtle population substructure between West Coast and Southern East Coast colonies and a weak, but significant (p = 0.01), isolation-by-distance pattern among the eight colonies studied. Furthermore, our demographic reconstructions supported severe bottlenecks with potential 10-fold and 250-fold declines in response to Polynesian and European hunting, respectively. Finally, we were able to assign individuals treated as unknowns to their regions of origin with high confidence (96%) using our SNP data. Our results indicate that while it may be difficult to detect population structure in species that have experienced rapid recovery, next-generation markers and methods are powerful tools for resolving fine-scale structure and informing conservation and management efforts.

Full Mitogenomes in the Critically Endangered Kākāpō Reveal Major Post-Glacial and Anthropogenic Effects on Neutral Genetic Diversity

Understanding how species respond to population declines is a central question in conservation and evolutionary biology. Population declines are often associated with loss of genetic diversity, inbreeding and accumulation of deleterious mutations, which can lead to a reduction in fitness and subsequently contribute to extinction. Using temporal approaches can help us understand the effects of population declines on genetic diversity in real time. Sequencing pre-decline as well as post-decline mitogenomes representing all the remaining mitochondrial diversity, we estimated the loss of genetic diversity in the critically endangered kākāpō (Strigops habroptilus). We detected a signal of population expansion coinciding with the end of the Pleistocene last glacial maximum (LGM). Also, we found some evidence for northern and southern lineages, supporting the hypothesis that the species may have been restricted to isolated northern and southern refugia during the LGM. We observed an important loss of neutral genetic diversity associated with European settlement in New Zealand but we could not exclude a population decline associated with Polynesian settlement in New Zealand. However, we did not find evidence for fixation of deleterious mutations. We argue that despite high pre-decline genetic diversity, a rapid and range-wide decline combined with the lek mating system, and life-history traits of kākāpō contributed to a rapid loss of genetic diversity following severe population declines.

When new human-modified habitats favour the expansion of an amphibian pioneer species: Evolutionary history of the natterjack toad (Bufo calamita) in a coal basin

Human activities affect microevolutionary dynamics by inducing environmental changes. In particular, land cover conversion and loss of native habitats decrease genetic diversity and jeopardize the adaptive ability of populations. Nonetheless, new anthropogenic habitats can also promote the successful establishment of emblematic pioneer species. We investigated this issue by examining the population genetic features and evolutionary history of the natterjack toad (Bufo [Epidalea] calamita) in northern France, where populations can be found in native coastal habitats and coalfield habitats shaped by European industrial history, along with an additional set of European populations located outside this focal area. We predicted contrasting patterns of genetic structure, with newly settled coalfield populations departing from migration-drift equilibrium. As expected, coalfield populations showed a mosaic of genetically divergent populations with short-range patterns of gene flow, and native coastal populations indicated an equilibrium state with an isolation-by-distance pattern suggestive of postglacial range expansion. However, coalfield populations exhibited (i) high levels of genetic diversity, (ii) no evidence of local inbreeding or reduced effective population size and (iii) multiple maternal mitochondrial lineages, a genetic footprint depicting independent colonization events. Furthermore, approximate Bayesian computations suggested several evolutionary trajectories from ancient isolation in glacial refugia during the Pleistocene, with biogeographical signatures of recent expansion probably confounded by human-mediated mixing of different lineages. From an evolutionary and conservation perspective, this study highlights the ecological value of industrial areas, provided that ongoing regional gene flow is ensured within the existing lineage boundaries.

Semen Collection and Spermatozoa Characteristics in the Kea Parrot (Nestor notabilis)

We describe the seminal characteristics of the kea parrot (Nestor notabilis), an endangered species endemic to the South Island of New Zealand. Semen was collected in the full breeding season from 6 birds in the collection of an amateur aviculturist. The manual massage technique was used. A total of 25 ejaculates was collected and evaluated for volume, degree of contamination, and spermatozoa concentration; motility and kinetic parameters were assessed on diluted samples (modified Tyrode's Albumin Lactate Pyruvate, pH 8.2, temperature 37.5°C) with a computer-aided sperm analyzer. Four ejaculates were not analyzed because of an excessively high degree of contamination. Semen color ranged from transparent or turbid yellow to whitish. The geometric mean of spermatozoa number/ejaculate was 765.9 ± 2084.7 × 10⁶. Total and progressive motility were 71.7% ± 20.0% and 59.8% ± 22.0%, respectively. Great variability was observed both among birds and among different ejaculates of the same subject. The seminal characteristics of kea are worth further investigation, with the aim of relating semen quality to fertility and defining a minimum inseminating dose for breeding purposes. A deeper knowledge of male reproductive biology also can increase the success of breeding programs and enable the use of the kea as a model species for other more threatened species, such as the kaka ( Nestor meridionalis ) and the kakapo (Strigops habroptila).

European Colonization, Not Polynesian Arrival, Impacted Population Size and Genetic Diversity in the Critically Endangered New Zealand Kākāpō

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 kākāpō (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 kākāpō decline, using microsatellites and mitochondrial DNA (mtDNA) to characterize recent kākāpō 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 kākāpō, 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 kākāpō. Our findings inform contemporary kākāpō 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.

Likelihood-Free Inference in High-Dimensional Models

Methods that bypass analytical evaluations of the likelihood function have become an indispensable tool for statistical inference in many fields of science. These so-called likelihood-free methods rely on accepting and rejecting simulations based on summary statistics, which limits them to low-dimensional models for which the value of the likelihood is large enough to result in manageable acceptance rates. To get around these issues, we introduce a novel, likelihood-free Markov chain Monte Carlo (MCMC) method combining two key innovations: updating only one parameter per iteration and accepting or rejecting this update based on subsets of statistics approximately sufficient for this parameter. This increases acceptance rates dramatically, rendering this approach suitable even for models of very high dimensionality. We further derive that for linear models, a one-dimensional combination of statistics per parameter is sufficient and can be found empirically with simulations. Finally, we demonstrate that our method readily scales to models of very high dimensionality, using toy models as well as by jointly inferring the effective population size, the distribution of fitness effects (DFE) of segregating mutations, and selection coefficients for each locus from data of a recent experiment on the evolution of drug resistance in influenza.

Ancient and contemporary DNA reveal a pre-human decline but no population bottleneck associated with recent human persecution in the kea (Nestor notabilis)

The impact of population bottlenecks is an important factor to consider when assessing species survival. Population declines can considerably limit the evolutionary potential of species and make them more susceptible to stochastic events. New Zealand has a well documented history of decline of endemic avifauna related to human colonization. Here, we investigate the genetic effects of a recent population decline in the endangered kea (Nestor notabilis). Kea have undergone a long-lasting persecution between the late 1800s to 1970s where an estimated 150,000 kea were culled under a governmental bounty scheme. Kea now number 1,000–5,000 individuals in the wild and it is likely that the recent population decline may have reduced the genetic diversity of the species. Comparison of contemporary (n = 410), historical (n = 15) and fossil samples (n = 4) showed a loss of mitochondrial diversity since the end of the last glaciation (Otiran Glacial) but no loss of overall genetic diversity associated with the cull. Microsatellite data indicated a recent bottleneck for only one population and a range-wide decline in N_e dating back some 300 – 6,000 years ago, a period predating European arrival in NZ. These results suggest that despite a recent human persecution, kea might have experienced a large population decline before stabilizing in numbers prior to human settlement of New Zealand in response to Holocene changes in habitat distribution. Our study therefore highlights the need to understand the respective effects of climate change and human activities on endangered species dynamics when proposing conservation guidelines.