The rate and nature of mitochondrial DNA mutations in human pedigrees.

We examined the rate and nature of mitochondrial DNA (mtDNA) mutations in humans using sequence data from 64,806 contemporary Icelanders from 2,548 matrilines. Based on 116,663 mother-child transmissions, 8,199 mutations were detected, providing robust rate estimates by nucleotide type, functional impact, position, and different alleles at the same position. We thoroughly document the true extent of hypermutability in mtDNA, mainly affecting the control region but also some coding-region variants. The results reveal the impact of negative selection on viable deleterious mutations, including rapidly mutating disease-associated 3243A>G and 1555A>G and pre-natal selection that most likely occurs during the development of oocytes. Finally, we show that the fate of new mutations is determined by a drastic germline bottleneck, amounting to an average of 3 mtDNA units effectively transmitted from mother to child.

Temporal dynamics of woolly mammoth genome erosion prior to extinction.

A number of species have recently recovered from near-extinction. Although these species have avoided the immediate extinction threat, their long-term viability remains precarious due to the potential genetic consequences of population declines, which are poorly understood on a timescale beyond a few generations. Woolly mammoths (Mammuthus primigenius) became isolated on Wrangel Island around 10,000 years ago and persisted for over 200 generations before becoming extinct around 4,000 years ago. To study the evolutionary processes leading up to the mammoths' extinction, we analyzed 21 Siberian woolly mammoth genomes. Our results show that the population recovered quickly from a severe bottleneck and remained demographically stable during the ensuing six millennia. We find that mildly deleterious mutations gradually accumulated, whereas highly deleterious mutations were purged, suggesting ongoing inbreeding depression that lasted for hundreds of generations. The time-lag between demographic and genetic recovery has wide-ranging implications for conservation management of recently bottlenecked populations.

Hominin population bottleneck coincided with migration from Africa during the Early Pleistocene ice age transition.

Two recently published analyses make cases for severe bottlenecking of human populations occurring in the late Early Pleistocene, one case at about 0.9 Mya based on a genomic analysis of modern human populations and the low number of hominin sites of this age in Africa and the other at about 1.1 Mya based on an age inventory of sites of hominin presence in Eurasia. Both models point to climate change as the bottleneck trigger, albeit manifested at very different times, and have implications for human migrations as a mechanism to elude extinction at bottlenecking. Here, we assess the climatic and chronologic components of these models and suggest that the several hundred-thousand-year difference is largely an artifact of biases in the chronostratigraphic record of Eurasian hominin sites. We suggest that the best available data are consistent with the Galerian hypothesis expanded from Europe to Eurasia as a major migration pulse of fauna including hominins in the late Early Pleistocene as a consequence of the opening of land routes from Africa facilitated by a large sea level drop associated with the first major ice age of the Pleistocene and concurrent with widespread aridity across Africa that occurred during marine isotope stage 22 at ~0.9 Mya. This timing agrees with the independently dated bottleneck from genomic analysis of modern human populations and allows speculations about the relative roles of climate forcing on the survival of hominins.

Increased homozygosity due to endogamy results in fitness consequences in a human population.

Recessive alleles have been shown to directly affect both human Mendelian disease phenotypes and complex traits. Pedigree studies also suggest that consanguinity results in increased childhood mortality and adverse health phenotypes, presumably through penetrance of recessive mutations. Here, we test whether the accumulation of homozygous, recessive alleles decreases reproductive success in a human population. We address this question among the Namibian Himba, an endogamous agro-pastoralist population, who until very recently practiced natural fertility. Using a sample of 681 individuals, we show that Himba exhibit elevated levels of "inbreeding," calculated as the fraction of the genome in runs of homozygosity (FROH). Many individuals contain multiple long segments of ROH in their genomes, indicating that their parents had high kinship coefficients. However, we do not find evidence that this is explained by first-cousin consanguinity, despite a reported social preference for cross-cousin marriages. Rather, we show that elevated haplotype sharing in the Himba is due to a bottleneck, likely in the past 60 generations. We test whether increased recessive mutation load results in observed fitness consequences by assessing the effect of FROH on completed fertility in a cohort of postreproductive women (n = 69). We find that higher FROH is significantly associated with lower fertility. Our data suggest a multilocus genetic effect on fitness driven by the expression of deleterious recessive alleles, especially those in long ROH. However, these effects are not the result of consanguinity but rather elevated background identity by descent.

Contrasting demographic history and mutational load in three threatened whitebark pines (Pinus subsect. Gerardianae): implications for conservation.

Demographic history and mutational load are of paramount importance for the adaptation of the endangered species. However, the effects of population evolutionary history and genetic load on the adaptive potential in endangered conifers remain unclear. Here, using population transcriptome sequencing, whole chloroplast genomes and mitochondrial DNA markers, combined with niche analysis, we determined the demographic history and mutational load for three threatened whitebark pines having different endangered statuses, Pinus bungeana, P. gerardiana and P. squamata. Demographic inference indicated that severe bottlenecks occurred in all three pines at different times, coinciding with periods of major climate and geological changes; in contrast, while P. bungeana experienced a recent population expansion, P. gerardiana and P. squamata maintained small population sizes after bottlenecking. Abundant homozygous-derived variants accumulated in the three pines, particularly in P. squamata, while the species with most heterozygous variants was P. gerardiana. Abundant moderately and few highly deleterious variants accumulated in the pine species that have experienced the most severe demographic bottlenecks (P. gerardiana and P. squamata), most likely because of purging effects. Finally, niche modeling showed that the distribution of P. bungeana might experience a significant expansion in the future, and the species' identified genetic clusters are also supported by differences in the ecological niche. The integration of genomic, demographic and niche data has allowed us to prove that the three threatened pines have contrasting patterns of demographic history and mutational load, which may have important implications in their adaptive potential and thus are also key for informing conservation planning.

Severe Bottleneck Impacted the Genomic Structure of Egg-Eating Cichlids in Lake Victoria.

Within 15,000 years, the explosive adaptive radiation of haplochromine cichlids in Lake Victoria, East Africa, generated 500 endemic species. In the 1980s, the upsurge of Nile perch, a carnivorous fish artificially introduced to the lake, drove the extinction of more than 200 endemic cichlids. The Nile perch predation particularly harmed piscivorous cichlids, including paedophages, cichlids eat eggs and fries, which is an example of the unique trophic adaptation seen in African cichlids. Here, aiming to investigate past demographic events possibly triggered by the invasion of Nile perch and the subsequent impacts on the genetic structure of cichlids, we conducted large-scale comparative genomics. We discovered evidence of recent bottleneck events in 4 species, including 2 paedophages, which began during the 1970s to 1980s, and population size rebounded during the 1990s to 2000s. The timing of the bottleneck corresponded to the historical records of endemic haplochromines" disappearance and later resurgence, which is likely associated with the introduction of Nile perch by commercial demand to Lake Victoria in the 1950s. Interestingly, among the 4 species that likely experienced bottleneck, Haplochromis sp. "matumbi hunter," a paedophagous cichlid, showed the most severe bottleneck signatures. The components of shared ancestry inferred by ADMIXTURE suggested a high genetic differentiation between matumbi hunter and other species. In contrast, our phylogenetic analyses highly supported the monophyly of the 5 paedophages, consistent with the results of previous studies. We conclude that high genetic differentiation of matumbi hunter occurred due to the loss of shared genetic components among haplochromines in Lake Victoria caused by the recent severe bottleneck.

Transmission Bottleneck Size Estimation from De Novo Viral Genetic Variation.

Sequencing of viral infections has become increasingly common over the last decade. Deep sequencing data in particular have proven useful in characterizing the roles that genetic drift and natural selection play in shaping within-host viral populations. They have also been used to estimate transmission bottleneck sizes from identified donor-recipient pairs. These bottleneck sizes quantify the number of viral particles that establish genetic lineages in the recipient host and are important to estimate due to their impact on viral evolution. Current approaches for estimating bottleneck sizes exclusively consider the subset of viral sites that are observed as polymorphic in the donor individual. However, these approaches have the potential to substantially underestimate true transmission bottleneck sizes. Here, we present a new statistical approach for instead estimating bottleneck sizes using patterns of viral genetic variation that arise de novo within a recipient individual. Specifically, our approach makes use of the number of clonal viral variants observed in a transmission pair, defined as the number of viral sites that are monomorphic in both the donor and the recipient but carry different alleles. We first test our approach on a simulated dataset and then apply it to both influenza A virus sequence data and SARS-CoV-2 sequence data from identified transmission pairs. Our results confirm the existence of extremely tight transmission bottlenecks for these 2 respiratory viruses.

Genome-Wide Allele Frequency Changes Reveal That Dynamic Metapopulations Evolve Differently.

Two important characteristics of metapopulations are extinction-(re)colonization dynamics and gene flow between subpopulations. These processes can cause strong shifts in genome-wide allele frequencies that are generally not observed in "classical" (large, stable, and panmictic) populations. Subpopulations founded by one or a few individuals, the so-called propagule model, are initially expected to show intermediate allele frequencies at polymorphic sites until natural selection and genetic drift drive allele frequencies toward a mutation-selection-drift equilibrium characterized by a negative exponential-like distribution of the site frequency spectrum. We followed changes in site frequency spectrum distribution in a natural metapopulation of the cyclically parthenogenetic pond-dwelling microcrustacean Daphnia magna using biannual pool-seq samples collected over a 5-yr period from 118 ponds occupied by subpopulations of known age. As expected under the propagule model, site frequency spectra in newly founded subpopulations trended toward intermediate allele frequencies and shifted toward right-skewed distributions as the populations aged. Immigration and subsequent hybrid vigor altered this dynamic. We show that the analysis of site frequency spectrum dynamics is a powerful approach to understand evolution in metapopulations. It allowed us to disentangle evolutionary processes occurring in a natural metapopulation, where many subpopulations evolve in parallel. Thereby, stochastic processes like founder and immigration events lead to a pattern of subpopulation divergence, while genetic drift leads to converging site frequency spectrum distributions in the persisting subpopulations. The observed processes are well explained by the propagule model and highlight that metapopulations evolve differently from classical populations.

Sorting of mitochondrial and plastid heteroplasmy in Arabidopsis is extremely rapid and depends on MSH1 activity.

The fate of new mitochondrial and plastid mutations depends on their ability to persist and spread among the numerous organellar genome copies within a cell (heteroplasmy). The extent to which heteroplasmies are transmitted across generations or eliminated through genetic bottlenecks is not well understood in plants, in part because their low mutation rates make these variants so infrequent. Disruption of MutS Homolog 1 (MSH1), a gene involved in plant organellar DNA repair, results in numerous de novo point mutations, which we used to quantitatively track the inheritance of single nucleotide variants in mitochondrial and plastid genomes in Arabidopsis. We found that heteroplasmic sorting (the fixation or loss of a variant) was rapid for both organelles, greatly exceeding rates observed in animals. In msh1 mutants, plastid variants sorted faster than those in mitochondria and were typically fixed or lost within a single generation. Effective transmission bottleneck sizes (N) for plastids and mitochondria were N ∼ 1 and 4, respectively. Restoring MSH1 function further increased the rate of heteroplasmic sorting in mitochondria (N ∼ 1.3), potentially because of its hypothesized role in promoting gene conversion as a mechanism of DNA repair, which is expected to homogenize genome copies within a cell. Heteroplasmic sorting also favored GC base pairs. Therefore, recombinational repair and gene conversion in plant organellar genomes can potentially accelerate the elimination of heteroplasmies and bias the outcome of this sorting process.

Hole Relaxation Bottlenecks in CdSe/CdTe/CdSe Lateral Heterostructures Lead to Bicolor Emission.

Concentric lateral CdSe/CdTe/CdSe heterostructures show bicolor photoluminescence from both a red charge transfer band of the CdSe/CdTe interface and a green fluorescence from CdSe. This work uses visible and near-infrared transient spectroscopy measurements to demonstrate that the deviation from Kasha's rule arises from a hole relaxation bottleneck from CdSe to CdTe. Hole transfer can take up to 1 ns, which permits radiative relaxation of excitons remaining in CdSe. Simulations indicate that the hole relaxation bottleneck arises due to the sparse density of states and poor spatial overlap of hole states at energies near the CdSe band edge. The divergent kinetics of transfer for band edge and hot holes is exploited to vary the ratio of green and red photoluminescence with excitation wavelength, providing another knob to control emission color. These findings support the use of lateral heterojunctions as a method for slowing carrier relaxation in two-dimensional materials.

Excitonic Thermalization Bottleneck in Twisted TMD Heterostructures.

Twisted van der Waals heterostructures show intriguing interface exciton physics, including hybridization effects and emergence of moiré potentials. Recent experiments have revealed that moiré-trapped excitons exhibit remarkable dynamics, where excited states show lifetimes that are several orders of magnitude longer than in monolayers. The origin of this behavior is still under debate. Based on a microscopic many-particle approach, we investigate the phonon-driven relaxation cascade of nonequilibrium moiré excitons in the exemplary MoSe2-WSe2 heterostructure. We track exciton relaxation pathways across different moiré mini-bands and identify the phonon-scattering channels assisting the spatial redistribution of excitons into low-energy pockets of the moiré potential. We unravel a phonon bottleneck in the flat band structure at low twist angles preventing excitons from fully thermalizing into the lowest state, explaining the measured enhanced emission intensity and lifetime of excited moiré excitons. Overall, our work provides important insights into exciton relaxation dynamics in flat-band exciton materials.

Unveiling Spatiotemporal Diffusion of Hot Carriers Influenced by Spatial Nonuniform Hot Phonon Bottleneck Effect in Monolayer MoS2.

The exceptional semiconducting properties of two-dimensional (2D) transition metal dichalcogenides (TMDs) have made them highly promising for the development of future electronic and optoelectronic devices. Extensive studies of TMDs are partly associated with their ability to generate 2D-confined hot carriers above the conduction band edges, enabling potential applications that rely on such transient excited states. In this work, room-temperature spatiotemporal hot carrier dynamics in monolayer MoS2 is studied by transient absorption microscopy (TAM), featuring an initial ultrafast expansion followed by a rapid negative diffusion, and ultimately a slow long-term expansion of the band edge C-excitons. We provide direct experimental evidence to identify the abnormal negative diffusion process as a spatial contraction of the hot carriers resulting from spatial variation in the hot phonon bottleneck effect due to the Gaussian intensity distribution of the pump laser beam.

Genomics Reveals Complex Population History and Unexpected Diversity of Eurasian Otters (Lutra lutra) in Britain Relative to Genetic Methods.

Conservation genetic analyses of many endangered species have been based on genotyping of microsatellite loci and sequencing of short fragments of mtDNA. The increase in power and resolution afforded by whole genome approaches may challenge conclusions made on limited numbers of loci and maternally inherited haploid markers. Here, we provide a matched comparison of whole genome sequencing versus microsatellite and control region (CR) genotyping for Eurasian otters (Lutra lutra). Previous work identified four genetically differentiated "stronghold" populations of otter in Britain, derived from regional populations that survived the population crash of the 1950s-1980s. Using whole genome resequencing data from 45 samples from across the British stronghold populations, we confirmed some aspects of population structure derived from previous marker-driven studies. Importantly, we showed that genomic signals of the population crash bottlenecks matched evidence from otter population surveys. Unexpectedly, two strongly divergent mitochondrial lineages were identified that were undetectable using CR fragments, and otters in the east of England were genetically distinct and surprisingly variable. We hypothesize that this previously unsuspected variability may derive from past releases of Eurasian otters from other, non-British source populations in England around the time of the population bottleneck. Our work highlights that even reasonably well-studied species may harbor genetic surprises, if studied using modern high-throughput sequencing methods.

Genus-Wide Transcriptional Landscapes Reveal Correlated Gene Networks Underlying Microevolutionary Divergence in Diatoms.

Marine microbes like diatoms make up the base of marine food webs and drive global nutrient cycles. Despite their key roles in ecology, biogeochemistry, and biotechnology, we have limited empirical data on how forces other than adaptation may drive diatom diversification, especially in the absence of environmental change. One key feature of diatom populations is frequent extreme reductions in population size, which can occur both in situ and ex situ as part of bloom-and-bust growth dynamics. This can drive divergence between closely related lineages, even in the absence of environmental differences. Here, we combine experimental evolution and transcriptome landscapes (t-scapes) to reveal repeated evolutionary divergence within several species of diatoms in a constant environment. We show that most of the transcriptional divergence can be captured on a reduced set of axes, and that repeatable evolution can occur along a single major axis of variation defined by core ortholog expression comprising common metabolic pathways. Previous work has associated specific transcriptional changes in gene networks with environmental factors. Here, we find that these same gene networks diverge in the absence of environmental change, suggesting these pathways may be central in generating phenotypic diversity as a result of both selective and random evolutionary forces. If this is the case, these genes and the functions they encode may represent universal axes of variation. Such axes that capture suites of interacting transcriptional changes during diversification improve our understanding of both global patterns in local adaptation and microdiversity, as well as evolutionary forces shaping algal cultivation.

A Selective Bottleneck During Host Entry Drives the Evolution of New Legume Symbionts.

During the emergence of new host-microbe symbioses, microbial fitness results from the ability to complete the different steps of symbiotic life cycles, where each step imposes specific selective pressures. However, the relative contribution of these different selective pressures to the adaptive trajectories of microbial symbionts is still poorly known. Here, we characterized the dynamics of phenotypic adaptation to a simplified symbiotic life cycle during the experimental evolution of a plant pathogenic bacterium into a legume symbiont. We observed that fast adaptation was predominantly explained by improved competitiveness for host entry, which outweighed adaptation to within-host proliferation. Whole-population sequencing of bacteria at regular time intervals along this evolution experiment revealed the continuous accumulation of new mutations (fuelled by a transient hypermutagenesis phase occurring at each cycle before host entry, a phenomenon described in previous work) and sequential sweeps of cohorts of mutations with similar temporal trajectories. The identification of adaptive mutations within the fixed mutational cohorts showed that several adaptive mutations can co-occur in the same cohort. Moreover, all adaptive mutations improved competitiveness for host entry, while only a subset of those also improved within-host proliferation. Computer simulations predict that this effect emerges from the presence of a strong selective bottleneck at host entry occurring before within-host proliferation and just after the hypermutagenesis phase in the rhizosphere. Together, these results show how selective bottlenecks can alter the relative influence of selective pressures acting during bacterial adaptation to multistep infection processes.

Patterns of genetic connectedness between modern and medieval Estonian genomes reveal the origins of a major ancestry component of the Finnish population.

The Finnish population is a unique example of a genetic isolate affected by a recent founder event. Previous studies have suggested that the ancestors of Finnic-speaking Finns and Estonians reached the circum-Baltic region by the 1st millennium BC. However, high linguistic similarity points to a more recent split of their languages. To study genetic connectedness between Finns and Estonians directly, we first assessed the efficacy of imputation of low-coverage ancient genomes by sequencing a medieval Estonian genome to high depth (23×) and evaluated the performance of its down-sampled replicas. We find that ancient genomes imputed from >0.1× coverage can be reliably used in principal-component analyses without projection. By searching for long shared allele intervals (LSAIs; similar to identity-by-descent segments) in unphased data for >143,000 present-day Estonians, 99 Finns, and 14 imputed ancient genomes from Estonia, we find unexpectedly high levels of individual connectedness between Estonians and Finns for the last eight centuries in contrast to their clear differentiation by allele frequencies. High levels of sharing of these segments between Estonians and Finns predate the demographic expansion and late settlement process of Finland. One plausible source of this extensive sharing is the 8th-10th centuries AD migration event from North Estonia to Finland that has been proposed to explain uniquely shared linguistic features between the Finnish language and the northern dialect of Estonian and shared Christianity-related loanwords from Slavic. These results suggest that LSAI detection provides a computationally tractable way to detect fine-scale structure in large cohorts.

Overexpression of BAPT and DBTNBT genes in Taxus baccata in vitro cultures to enhance the biotechnological production of paclitaxel.

Paclitaxel is one of the most effective anticancer drugs ever developed. Although the most sustainable approach to its production is provided by plant cell cultures, the yield is limited by bottleneck enzymes in the taxane biosynthetic pathway: baccatin-aminophenylpropanoyl-13-O-transferase (BAPT) and 3'-N-debenzoyltaxol N-benzoyltransferase (DBTNBT). With the aim of enhancing paclitaxel production by overcoming this bottleneck, we obtained distinct lines of Taxus baccata in vitro roots, each independently overexpressing either of the two flux-limiting genes, BAPT or DBTNBT, through a Rhizobium rhizogenes A4-mediated transformation. Due to the slow growth rate of the transgenic Taxus roots, they were dedifferentiated to obtain callus lines and establish cell suspensions. The transgenic cells were cultured in a two-stage system and stimulated for taxane production by a dual elicitation treatment with 1 µm coronatine plus 50 mm of randomly methylated-ß-cyclodextrins. A high overexpression of BAPT (59.72-fold higher at 48 h) and DBTNBT (61.93-fold higher at 72 h) genes was observed in the transgenic cell cultures, as well as an improved taxane production. Compared to the wild type line (71.01 mg/L), the DBTNBT line produced more than four times higher amounts of paclitaxel (310 mg/L), while the content of this taxane was almost doubled in the BAPT line (135 mg/L). A transcriptional profiling of taxane biosynthetic genes revealed that GGPPS, TXS and DBAT genes were the most reactive to DBTNBT overexpression and the dual elicitation, their expression increasing gradually and constantly. The same genes exhibited a pattern of isolated peaks of expression in the elicited BAPT-overexpressing line.

Automated MRI-based segmentation of intracranial arterial calcification by restricting feature complexity.

PURPOSE: To develop an automated deep learning model for MRI-based segmentation and detection of intracranial arterial calcification. METHODS: A novel deep learning model under the variational autoencoder framework was developed. A theoretically grounded dissimilarity loss was proposed to refine network features extracted from MRI and restrict their complexity, enabling the model to learn more generalizable MR features that enhance segmentation accuracy and robustness for detecting calcification on MRI. RESULTS: The proposed method was compared with nine baseline methods on a dataset of 113 subjects and showed superior performance (for segmentation, Dice similarity coefficient: 0.620, area under precision-recall curve [PR-AUC]: 0.660, 95% Hausdorff Distance: 0.848 mm, Average Symmetric Surface Distance: 0.692 mm; for slice-wise detection, F1 score: 0.823, recall: 0.764, precision: 0.892, PR-AUC: 0.853). For clinical needs, statistical tests confirmed agreement between the true calcification volumes and predicted values using the proposed approach. Various MR sequences, namely T1, time-of-flight, and SNAP, were assessed as inputs to the model, and SNAP provided unique and essential information pertaining to calcification structures. CONCLUSION: The proposed deep learning model with a dissimilarity loss to reduce feature complexity effectively improves MRI-based identification of intracranial arterial calcification. It could help establish a more comprehensive and powerful pipeline for vascular image analysis on MRI.

Stochastic organelle genome segregation through Arabidopsis development and reproduction.

Organelle DNA (oDNA) in mitochondria and plastids is vital for plant (and eukaryotic) life. Selection against damaged oDNA is mediated in part by segregation - sorting different oDNA types into different cells in the germline. Plants segregate oDNA very rapidly, with oDNA recombination protein MSH1 a key driver of this segregation, but we have limited knowledge of the dynamics of this segregation within plants and between generations. Here, we reveal how oDNA evolves through Arabidopsis thaliana development and reproduction. We combine stochastic modelling, Bayesian inference, and model selection with new and existing tissue-specific oDNA measurements from heteroplasmic Arabidopsis plant lines through development and between generations. Segregation proceeds gradually but continually during plant development, with a more rapid increase between inflorescence formation and the next generation. When MSH1 is compromised, the majority of observed segregation can be achieved through partitioning at cell divisions. When MSH1 is functional, mtDNA segregation is far more rapid; we show that increased oDNA gene conversion is a plausible mechanism quantitatively explaining this acceleration. These findings reveal the quantitative, time-dependent details of oDNA segregation in Arabidopsis. We also discuss the support for different models of the plant germline provided by these observations.

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

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