Explosive ice age diversification of kiwi

Morphological and molecular description of a novel species of Eimeria (Apicomplexa) that infects extraintestinal tissues of kiwi (Aves: Apteryx spp.)

Coccidia (Apiconmplexa) are naturally occurring and occasionally detrimental parasites of kiwi (Apteryx spp.), a unique, flightless bird species dependent upon conservation efforts for survival. Using morphological and molecular data, a new coccidia species, Eimeria koka n. sp., isolated from two closely related but geographically isolated kiwi host species, Apteryx rowi Tennyson et al. (rowi) and Apteryx mantelli Bartlett (North Island brown kiwi), is described. Oocysts are oval (20.8 × 15.9 μm) with a mean L/W ratio of 1.3, and a distinctive rough, crenellated brown oocyst wall (mean 1.2 μm), an oocyst residuum, 1-2 polar granules, and no micropyle. Sporocysts are ellipsoidal (11.6 × 6.3 μm) with a Stieda body and sporocyst residuum. Phylogenetic analysis of the cytochrome C oxidase subunit 1 (CO1) placed E. koka n. sp. in a separate clade to other Eimeria species previously identified from kiwi (Coker et al., Syst Parasitol 100(3):269-281, 2023). Comparison of DNA from oocysts with infected tissues from a single juvenile North Island brown kiwi confirmed parasitism of the kidney and lung tissues. This is the first Eimeria species identified from extraintestinal tissues in kiwi. Further molecular studies are recommended to determine the tissue distribution of E. koka n. sp. and other Eimeria species in kiwi.

Plio-Pleistocene Environmental Changes Drove the Settlement of Aotearoa New Zealand by Australian Open-Habitat Bird Lineages

In a changing environment, vacant niches can be filled either by adaptation of local taxa or range-expanding invading species. The relative tempo of these patterns is of key interest in the modern age of climate change. Aotearoa New Zealand has been a hotspot of biogeographic research for decades due to its long-term isolation and dramatic geological history. An island with high levels of faunal endemicity, it is a system well suited to studying the relative effects of in situ evolution versus dispersal in determining faunal assemblages, while its turbulent climate and geological history provide valuable insights into the evolutionary impacts of environmental changes. Such investigations are of urgent importance given predicted climate change and human impacts rapidly affecting environments globally. Here, we analyse the divergence dates of nearly all endemic Aotearoa New Zealand bird species from their overseas relatives to assess the role of environmental changes in driving speciation and colonisation, with special regard to cooling climate during the Pliocene and Pleistocene. We uncover a wave of colonisation events by Australian open-habitat adapted species since the Pliocene that peaked at the beginning of the Pleistocene. Furthermore, we highlight an even distribution of divergence dates in forest-adapted taxa through time, consistent with millions of years of extensive forest cover. Finally, we note parallels to the modern-day establishment of new bird populations from Australia and suggest this is largely influenced by anthropogenic land-use patterns. This research contributes to the growing body of work recognising the long-lasting impacts of Pleistocene climate change on Aotearoa New Zealand's avifauna, and reinforces biological invasions as a key evolutionary response to changing environmental conditions.

Admixture and environmental fluctuations shape the evolutionary history of a predator radiation in East Africa's Lake Tanganyika

Top predators have oversized impacts on food webs and ecosystem dynamics, and introducing a novel predator to a naive environment can have dramatic consequences for endemic biodiversity. Lake Tanganyika is unique among African lakes in the diversity of the pelagic top predators in the genus Lates, where four species are endemic to the lake. Using a combination of reduced-representation and whole genome resequencing data, and pairing these with phylogenetic and demographic modeling approaches, we find that Lates colonization of Lake Tanganyika was much more recent (~1-2 Mya) than other major and diverse clades within the lake. Demographic modeling suggests that diversification among Lates species within the lake occurred during a time period of dramatic changes in lake levels driven by glacial-interglacial cycles, supporting a role of these fluctuations as a "species pump" for lacustrine taxa. We further find that these lake level fluctuations likely contributed to multiple bouts of admixture among Lates species during the mid- to late-Pleistocene (~90-500 Kya). Together, our findings suggest a dynamic and environmentally linked evolutionary history of the Lates radiation with the potential for dramatic ecosystem consequences for the taxa already present in Lake Tanganyika prior to Lates colonization and diversification.

Significance Statement

When introduced to novel ecosystems, top predators can cause major alterations to biodiversity and food webs. Species introductions to novel habitats can also provide invading taxa with ecological opportunities that facilitate evolutionary diversification. Here, we find evidence that the radiation of endemic top predators in East Africa's Lake Tanganyika originated surprisingly recently, and that these species have experienced periods of hybridization with a widespread riverine relative throughout their history. These findings have major implications for the history of the lake and suggest that the introduction of Nile perch into Lake Victoria, which caused dramatic ecosystem and food web changes, may be a contemporary analog to the historical events in Lake Tanganyika.

The ghost of ice ages past: Impact of Last Glacial Maximum landscapes on modern biodiversity

Modeled modern and Last Glacial Maximum (LGM) climate ranges for 47 genetically confirmed small Holarctic land snails documented profound landscape dynamism over the last 21,000 years. Following deglaciation, range areas tended to increase by 50% while isolating barrier widths were cut in half. At the same time, the nature of isolating barriers underwent profound change, with the North American continental ice sheet becoming as important in the LGM as the Atlantic Ocean is today in separating Nearctic and Palearctic faunas. Because appropriate modern climate occurs for these species throughout the Holarctic, with no clear barriers being present-especially for such efficient passive dispersers-the current >90% turnover observed between Eurasian and North American species pools appears at least in part related to the LGM landscape. Understanding current and predicting potential future biodiversity patterns thus requires consideration of the landscape template across at least 15,000 years time scales.

Linking the spatial and genomic structure of adaptive potential for conservation management: a review

We unified the recent literature with the goal to contribute to the discussion on how genetic diversity might best be conserved. We argue that this decision will be guided by how genomic variation is distributed among manageable populations (i.e., its spatial structure), the degree to which adaptive potential is best predicted by variation across the entire genome or the subset of that variation that is identified as putatively adaptive (i.e., its genomic structure), and whether we are managing species as single entities or as collections of diversifying lineages. The distribution of genetic variation and our ultimate goal will have practical implications for on-the-ground management. If adaptive variation is largely polygenic or responsive to change, its spatial structure might be broadly governed by the forces determining genome-wide variation (linked selection, drift, and gene flow), making measurement and prioritization straightforward. If we are managing species as single entities, then population-level prioritization schemes are possible so as to maximize future pooled genetic variation. We outline one such scheme based on the popular Shapley value from cooperative game theory that considers the relative genetic contribution of a population to an unknown future collection of populations.

Introgression Underlies Phylogenetic Uncertainty But Not Parallel Plumage Evolution in a Recent Songbird Radiation

Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails," a group of 5 species in the songbird genus Motacilla: 1 species, Motacilla alba, shows wide intra-specific plumage variation, while the 4r others form 2 pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the 2 species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome to quantify the extent of discordance across genomic regions and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investing genome-wide patterns of gene tree heterogeneity to help understand the mechanisms underlying phenotypic evolution. [Gene tree heterogeneity; incomplete lineage sorting; introgression; parallel evolution; phylogenomics; plumage evolution; wagtails.].

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.

Pleistocene diversification of unifoliolate-leaved Lupinus (Leguminosae: Papilionoideae) in Florida

The importance and prevalence of recent ice-age and post-glacial speciation and species diversification during the Pleistocene across many organismal groups and physiographic settings are well established. However, the extent to which Pleistocene diversification can be attributed to climatic oscillations and their effects on distribution ranges and population structure remains debatable. In this study, we use morphologic, geographic and genetic (RADseq) data to document Pleistocene speciation and intra-specific diversification of the unifoliolate-leaved clade of Florida Lupinus, a small group of species largely restricted to inland and coastal sand ridges across the Florida peninsula and panhandle. Phylogenetic and demographic analyses alongside morphological and geographic evidence suggest that recent speciation and intra-specific divergence within this clade were driven by a combination of non-adaptive allopatric divergence caused by edaphic niche conservatism and opportunities presented by the emergence of new post-glacial sand ridge habitats. These results highlight the central importance of even modest geographic isolation and short periods of allopatric divergence following range expansion in the emergence of new taxa and add to the growing evidence that Pleistocene climatic oscillations may contribute to rapid diversification in a myriad of physiographic settings. Furthermore, our results shed new light on long-standing taxonomic debate surrounding the number of species in the Florida unifoliate Lupinus clade providing support for recognition of five species and a set of intra-specific variants. The important conservation implications for the narrowly restricted, highly endangered species Lupinus aridorum, which we show to be genetically distinct from its sister species Lupinus westianus, are discussed.

Beringia as a high-latitude engine of avian speciation

Beringia is a biogeographically dynamic region that extends from northeastern Asia into northwestern North America. This region has affected avian divergence and speciation in three important ways: (i) by serving as a route for intercontinental colonisation between Asia and the Americas; (ii) by cyclically splitting (and often reuniting) populations, subspecies, and species between these continents; and (iii) by providing isolated refugia through glacial cycles. The effects of these processes can be seen in taxonomic splits of shallow to increasing depths and in the presence of regional endemics. We review the taxa involved in the latter two processes (splitting-reuniting and isolation), with a focus on three research topics: avian diversity, time estimates of the generation of that diversity, and the regions within Beringia that might have been especially important. We find that these processes have generated substantial amounts of avian diversity, including 49 pairs of avian subspecies or species whose breeding distributions largely replace one another across the divide between the Old World and the New World in Beringia, and 103 avian species and subspecies endemic to this region. Among endemics, about one in three is recognised as a full biological species. Endemic taxa in the orders Charadriiformes (shorebirds, alcids, gulls, and terns) and Passeriformes (perching birds) are particularly well represented, although they show very different levels of diversity through evolutionary time. Endemic Beringian Charadriiformes have a 1.31:1 ratio of species to subspecies. In Passeriformes, endemic taxa have a 0.09:1 species-to-subspecies ratio, suggesting that passerine (and thus terrestrial) endemism might be more prone to long-term extinction in this region, although such 'losses' could occur through their being reconnected with wider continental populations during favourable climatic cycles (e.g. subspecies reintegration with other populations). Genetic evidence suggests that most Beringian avian taxa originated over the past 3 million years, confirming the importance of Quaternary processes. There seems to be no obvious clustering in their formation through time, although there might be temporal gaps with lower rates of diversity generation. For at least 62 species, taxonomically undifferentiated populations occupy this region, providing ample potential for future evolutionary diversification.

Microstructural and crystallographic evolution of palaeognath (Aves) eggshells

The avian palaeognath phylogeny has been recently revised significantly due to the advancement of genome-wide comparative analyses and provides the opportunity to trace the evolution of the microstructure and crystallography of modern dinosaur eggshells. Here, eggshells of all major clades of Palaeognathae (including extinct taxa) and selected eggshells of Neognathae and non-avian dinosaurs are analysed with electron backscatter diffraction. Our results show the detailed microstructures and crystallographies of (previously) loosely categorized ostrich-, rhea-, and tinamou-style morphotypes of palaeognath eggshells. All rhea-style eggshell appears homologous, while respective ostrich-style and tinamou-style morphotypes are best interpreted as homoplastic morphologies (independently acquired). Ancestral state reconstruction and parsimony analysis additionally show that rhea-style eggshell represents the ancestral state of palaeognath eggshells both in microstructure and crystallography. The ornithological and palaeontological implications of the current study are not only helpful for the understanding of evolution of modern and extinct dinosaur eggshells, but also aid other disciplines where palaeognath eggshells provide useful archive for comparative contrasts (e.g. palaeoenvironmental reconstructions, geochronology, and zooarchaeology).

Genomic insights into the evolutionary relationships and demographic history of kiwi

Kiwi are a unique and emblematic group of birds endemic to New Zealand. Deep-time evolutionary relationships among the five extant kiwi species have been difficult to resolve, in part due to the absence of pre-Quaternary fossils to inform speciation events. Here, we utilise single representative nuclear genomes of all five extant kiwi species (great spotted kiwi, little spotted kiwi, Okarito brown kiwi, North Island brown kiwi, and southern brown kiwi) and investigate their evolutionary histories with phylogenomic, genetic diversity, and deep-time (past million years) demographic analyses. We uncover relatively low levels of gene-tree phylogenetic discordance across the genomes, suggesting clear distinction between species. However, we also find indications of post-divergence gene flow, concordant with recent reports of interspecific hybrids. The four species for which unbiased levels of genetic diversity could be calculated, due to the availability of reference assemblies (all species except the southern brown kiwi), show relatively low levels of genetic diversity, which we suggest reflects a combination of older environmental as well as more recent anthropogenic influence. In addition, we suggest hypotheses regarding the impact of known past environmental events, such as volcanic eruptions and glacial periods, on the similarities and differences observed in the demographic histories of the five kiwi species over the past million years.

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.

Phylogeography and Cryptic Species Structure of a Locally Adapted Parasite in New Zealand

Phylogeographic patterns of many taxa on New Zealand's South Island are characterised by disjunct distributions that have been attributed to Pleistocene climatic cycles and the formation of the Southern Alps. Pleistocene glaciation has also been implicated in shaping the contemporary genetic differentiation between populations of the aquatic snail Potamopyrgus antipodarum. We investigated whether similar phylogeographic patterns exist for the snail's locally adapted trematode parasite, Atriophallophorus winterbourni. We found evidence for a barrier to gene-flow in sympatry between cryptic, but ecologically divergent species. When focusing on the most common of these species, disjunct geographic distributions are found for mitochondrial lineages that diverged during the Pleistocene. The boundary between these distributions is found in the central part of the South Island and reinforced by low cross-alpine migration. Further support for a vicariant origin of the phylogeographic pattern was found when assessing nuclear multilocus SNP data. Nuclear and mitochondrial population differentiation was concordant in pattern, except for populations in a potential secondary contact zone. Additionally, we found larger than expected differentiation between nuclear- and mitochondrial-based empirical Bayes F_ST estimates (global F_ST : 0.02 versus 0.39 for nuclear and mitochondrial data, respectively). Population subdivision is theoretically expected to be stronger for mitochondrial genomes due to a smaller effective population size, but the strong difference here, together with mito-nuclear discordance in a putative contact zone, is potentially indicative of divergent gene flow of nuclear and mitochondrial genomes.

Legacy of supervolcanic eruptions on population genetic structure of brown kiwi

Supervolcanoes are volcanoes capable of mega-colossal eruptions that emit more than 1,000 km³ of ash and other particles.¹ The earth's most recent mega-colossal eruption was the Oruanui eruption of the Taupo supervolcano 25,580 years before present (YBP) on the central North Island of New Zealand.² This eruption blanketed major swaths of the North Island in thick layers of ash and igneous rock,^2,3 devastating habitats and likely causing widespread population extinctions.^4-7 An additional devastating super-colossal eruption (>100 km³) of the Taupo supervolcano occurred approximately 1,690 YBP.⁸ The impacts of such massive but ephemeral natural disasters on contemporary population genetic structure remain underexplored. Here, we combined data for 4,951 SNPs with spatially explicit demographic and coalescent models within an approximate Bayesian computation framework to test the drivers of genetic structure in brown kiwi (Apteryx mantelli). Our results strongly support the importance of eruptions of the Taupo supervolcano in restructuring pre-existing geographic patterns of population differentiation and genetic diversity. Range shifts due to climatic oscillations-a frequent explanation for genetic structure⁹-are insufficient to fully explain the empirical data. Meanwhile, recent range contraction and fragmentation due to historically documented anthropogenic habitat alteration adds no explanatory power to our models. Our results support a major role for cycles of destruction and post-volcanic recolonization in restructuring the population genomic landscape of brown kiwi and highlight how ancient and ephemeral mega-disasters may leave a lasting legacy on patterns of intraspecific genetic variation.

Glacial vicariance and oceanic circulation shape population structure of the coastal legume Derris trifoliata in the Indo-West Pacific

PREMISE

The phylogeography of coastal plant species is shaped by contemporary and historical biogeographic processes. In this study, we aim to decipher the phylogeography of Derris trifoliata, a woody legume of relatively recent origin and wide distribution, in coastal areas in the Indo-West Pacific (IWP) region.

METHODS

Genetic diversity and population structure were assessed by analyzing six nuclear and three chloroplast DNA sequences from 30 populations across the species' range. Phylogeography was inferred by estimating gene flow, divergence time, historical population size changes, and historical habitat suitability using paleoclimatic niche modeling.

RESULTS

High genetic diversity was observed at the species level. The populations of three oceanic regions included in this study (i.e., Indian Ocean, South China Sea, and Pacific Ocean) formed distinct clades and likely diverged during the late Pleistocene. Potential barriers to gene flow were identified, including the Sunda and Sahul shelves, geographic distance, and current patterns of oceanic circulation. Analysis of changes in population size supported the bottleneck model, which was strengthened by estimates of habitat suitability across paleoclimatic conditions.

CONCLUSIONS

The once widespread distribution of D. trifoliata was fragmented by changes in climatic suitability and biogeographic barriers that arose following sea-level changes during the Pleistocene. In addition, contemporary patterns of oceanic circulation and geographic distance between populations appear to maintain genetic differentiation across its distribution in the IWP.

Phylogeography reveals the complex impact of the Last Glacial Maximum on New Zealand's terrestrial biota

We review the major phylogeographic patterns in Aotearoa New Zealand's terrestrial flora and fauna that have been associated with the Ōtira Glaciation of the Pleistocene, the end of which coincides with the global Last Glacial Maximum (LGM). We focus on (1) the complexity of biogeographic histories of New Zealand species, with LGM-driven phenomena overlaying the impacts of mountain-building and other drivers of phylogeographic structure; (2) the locations of glacial refugia and sets of taxa which may have shared refugia; and (3) the role of glaciation in driving diversification. We end with a brief focus on the next directions, including what can we learn about New Zealand's glacial history by expanding our phylogeographic toolbox to include genomic methods and hypothesis-driven inference methods. We provide follow-up questions which take advantage of the wealth of phylogeographic data for New Zealand.

Genetic evidence for post-glacial expansion from a southern refugium in the eastern moa (Emeus crassus)

Cycles of glacial expansion and contraction throughout the Pleistocene drove increases and decreases, respectively, in the geographical range and population size of many animal species. Genetic data have revealed that during glacial maxima the distribution of many Eurasian animals was restricted to small refugial areas, from which species expanded to reoccupy parts of their former range as the climate warmed. It has been suggested that the extinct eastern moa (Emeus crassus)-a large, flightless bird from New Zealand-behaved analogously during glacial maxima, possibly surviving only in a restricted area of lowland habitat in the southern South Island of New Zealand during the Last Glacial Maximum (LGM). However, previous studies have lacked the power and geographical sampling to explicitly test this hypothesis using genetic data. Here we analyse 46 ancient mitochondrial genomes from Late Pleistocene and Holocene bones of the eastern moa from across their post-LGM distribution. Our results are consistent with a post-LGM increase in the population size and genetic diversity of eastern moa. We also demonstrate that genetic diversity was higher in eastern moa from the southern extent of their range, supporting the hypothesis that they expanded from a single glacial refugium following the LGM.

Insights into Aotearoa New Zealand's biogeographic history provided by the study of natural hybrid zones

Hybridisation is commonly observed in geographical zones of contact between distinct lineages. These contact zones have long been of interest for biogeographers because they provide insight into the evolutionary and ecological processes that influence the distribution of species as well as the process of speciation. Here we review research on hybrid zones and zones of past introgression, both terrestrial and marine, in Aotearoa New Zealand. Many of New Zealand's hybrid zones occur between lineages or species that diverged prior to the Last Glacial Maximum (LGM), with numerous divergences dating to the early Pleistocene or Pliocene. Few secondary contact zones have been detected in terrestrial plants and in marine taxa. This may reflect a lack of the intensive sampling required to detect hybrid zones in these groups but for plants may also indicate widespread Pleistocene survival across the country. Lastly, we suggest avenues for research into New Zealand hybrid zones that are likely to be fruitful.

Demographic decline and lineage-specific adaptations characterize New Zealand kiwi

Small and fragmented populations may become rapidly differentiated due to genetic drift, making it difficult to distinguish whether neutral genetic structure is a signature of recent demographic events, or of long-term evolutionary processes that could have allowed populations to adaptively diverge. We sequenced 52 whole genomes to examine Holocene demographic history and patterns of adaptation in kiwi (Apteryx), and recovered 11 strongly differentiated genetic clusters corresponding to previously recognized lineages. Demographic models suggest that all 11 lineages experienced dramatic population crashes relative to early- or mid-Holocene levels. Small population size is associated with low genetic diversity and elevated genetic differentiation (FST), suggesting that population declines have strengthened genetic structure and led to the loss of genetic diversity. However, population size is not correlated with inbreeding rates. Eight lineages show signatures of lineage-specific selective sweeps (284 sweeps total) that are unlikely to have been caused by demographic stochasticity. Overall, these results suggest that despite strong genetic drift associated with recent bottlenecks, most kiwi lineages possess unique adaptations and should be recognized as separate adaptive units in conservation contexts. Our work highlights how whole-genome datasets can address longstanding uncertainty about the evolutionary and conservation significance of small and fragmented populations of threatened species.

Genome Stability Is in the Eye of the Beholder: CR1 Retrotransposon Activity Varies Significantly across Avian Diversity

Since the sequencing of the zebra finch genome it has become clear that avian genomes, while largely stable in terms of chromosome number and gene synteny, are more dynamic at an intrachromosomal level. A multitude of intrachromosomal rearrangements and significant variation in transposable element (TE) content have been noted across the avian tree. TEs are a source of genome plasticity, because their high similarity enables chromosomal rearrangements through nonallelic homologous recombination, and they have potential for exaptation as regulatory and coding sequences. Previous studies have investigated the activity of the dominant TE in birds, chicken repeat 1 (CR1) retrotransposons, either focusing on their expansion within single orders, or comparing passerines with nonpasserines. Here, we comprehensively investigate and compare the activity of CR1 expansion across orders of birds, finding levels of CR1 activity vary significantly both between and within orders. We describe high levels of TE expansion in genera which have speciated in the last 10 Myr including kiwis, geese, and Amazon parrots; low levels of TE expansion in songbirds across their diversification, and near inactivity of TEs in the cassowary and emu for millions of years. CR1s have remained active over long periods of time across most orders of neognaths, with activity at any one time dominated by one or two families of CR1s. Our findings of higher TE activity in species-rich clades and dominant families of TEs within lineages mirror past findings in mammals and indicate that genome evolution in amniotes relies on universal TE-driven processes.

Major population splits coincide with episodes of rapid climate change in a forest-dependent bird

Climate change influences population demography by altering patterns of gene flow and reproductive isolation. Direct mutation rates offer the possibility for accurate dating on the within-species level but are currently only available for a handful of vertebrate species. Here, we use the first directly estimated mutation rate in birds to study the evolutionary history of pied flycatchers (Ficedula hypoleuca). Using a combination of demographic inference and species distribution modelling, we show that all major population splits in this forest-dependent system occurred during periods of increased climate instability and rapid global temperature change. We show that the divergent Spanish subspecies originated during the Eemian-Weichselian transition 115-104 thousand years ago (kya), and not during the last glacial maximum (26.5-19 kya), as previously suggested. The magnitude and rates of climate change during the glacial-interglacial transitions that preceded population splits in pied flycatchers were similar to, or exceeded, those predicted to occur in the course of the current, human-induced climate crisis. As such, our results provide a timely reminder of the strong impact that episodes of climate instability and rapid temperature changes can have on species' evolutionary trajectories, with important implications for the natural world in the Anthropocene.

Intraskeletal bone growth patterns in the North Island Brown Kiwi (Apteryx mantelli): Growth mark discrepancy and implications for extinct taxa

Osteohistology, the study of bone microstructure, provides an important avenue for assessing extinct and extant vertebrate growth and life history. Cortical vascularity and collagen fibre organization are direct reflections of growth rate, while bone growth marks are indicative of absolute age. However, each skeletal element has its own ontogenetic trajectory and microstructure of certain bones may not be a true representation of whole body growth. Extensive comparative study of modern taxa is required to resolve intraskeletal discrepancies among age, vascularity and tissue organization in extinct vertebrates. Despite their comparative utility, studies of bone microstructure in modern taxa are severely lacking. Here, we add to a growing comparative osteohistological database by describing (1) bone tissue organization, (2) growth mark count, (3) sexually dimorphic bone (e.g. medullary bone) and (4) secondary cortical reconstruction in the bone microstructure of a 14-year-old male and 5-year-old female North Island Brown Kiwi (Apteryx mantelli). Transverse and longitudinal histological ground sections were processed and described for femora, tibiotarsi, tarsometatarsi, humeri, ulnae and radii in both kiwis. Cortical bone can generally be described as parallel-fibered tissue, interrupted by cyclical growth marks, with vascular canals oriented longitudinally within primary and secondary osteons. Tissue morphologically resembling medullary bone is present in the hindlimbs of the female, and coarse compacted cancellous bone (CCCB) is found sporadically in the male and female hindlimbs. Lines of arrested growth (LAGs) are present in all hindlimb bones of both kiwi, but remodelling has obliterated all LAGs in the male ulnae and radii. LAG count varies intraskeletally, but large weight bearing elements such as femora and tibiotarsi have less remodelling and, thus, higher number of LAGs. LAG count did not match absolute age in any skeletal element; a maximum of seven LAGs are present in the male kiwi and a maximum of seven LAGs in the female kiwi. The tissue organization within the forelimbs and hindlimbs is reflective of the protracted growth strategy of the North Island Brown Kiwi and congruent with previous studies of the kiwi. LAGs were highly variable throughout the skeleton of the kiwi and a decoupling of age and LAG deposition is apparent from the male kiwi samples. Excess LAGs in the 5-year-old female kiwi may be a product of hatching, egg laying or captivity. Regardless, LAG count variation in the kiwi stresses the importance of intraskeletal sampling when assessing growth patterns of extinct taxa. An extensive ontogenetic sampling of kiwi is necessary for future investigations of bone growth patterns, CCCB formation, medullary bone and LAG deposition and obliteration in these elusive birds.

Seasonal migration patterns and the maintenance of evolutionary diversity in a cryptic bird radiation

Morphological differentiation associated with evolutionary diversification is often explained with adaptive benefits but the processes and mechanisms maintaining cryptic diversity are still poorly understood. Using genome‐wide data, we show here that the pale sand martin Riparia diluta in Central and East Asia consists of three genetically deeply differentiated lineages which vary only gradually in morphology but broadly reflect traditional taxonomy. We detected no signs of gene flow along the eastern edge of the Qinghai‐Tibetan plateau between lowland south‐eastern Chinese R. d. fohkienensis and high‐altitude R. d. tibetana. Largely different breeding and migration timing between these low and high altitude populations as indicated by phenology data suggests that allochrony might act as prezygotic isolation mechanism in the area where their ranges abut. Mongolian populations of R. d. tibetana, however, displayed signs of limited mixed ancestries with Central Asian R. d. diluta. Their ranges meet in the area of a well‐known avian migratory divide, where western lineages take a western migration route around the Qinghai‐Tibetan plateau to winter quarters in South Asia, and eastern lineages take an eastern route to Southeast Asia. This might also be the case between western R. d. diluta and eastern R. d. tibetana as indicated by differing wintering grounds. We hypothesize that hybrids might have nonoptimal intermediate migration routes and selection against them might restrict gene flow. Although further potential isolation mechanisms might exist in the pale sand martin, our study points towards contrasting migration behaviour as an important factor in maintaining evolutionary diversity under morphological stasis.

Mixed Mating in a Multi-Origin Population Suggests High Potential for Genetic Rescue in North Island Brown Kiwi, Apteryx mantelli

Reinforcement translocations are increasingly utilised in conservation with the goal of achieving genetic rescue. However, concerns regarding undesirable results, such as genetic homogenisation or replacement, are widespread. One factor influencing translocation outcomes is the rate at which the resident and the introduced individuals interbreed. Consequently, post-release mate choice is a key behaviour to consider in conservation planning. Here we studied mating, and its consequences for genomic admixture, in the North Island brown kiwi Apteryx mantelli population on Ponui Island which was founded by two translocation events over 50 years ago. The two source populations used are now recognised as belonging to two separate management units between which birds differ in size and are genetically differentiated. We examined the correlation between male and female morphometrics for 17 known pairs and quantified the relatedness of 20 pairs from this admixed population. In addition, we compared the genetic similarity and makeup of 106 Ponui Island birds, including 23 known pairs, to birds representing the source populations for the original translocations. We found no evidence for size-assortative mating. On the contrary, genomic SNP data suggested that kiwi of one feather did not flock together, meaning that mate choice resulted in pairing between individuals that were less related than expected by random chance. Furthermore, the birds in the current Ponui Island population were found to fall along a gradient of genomic composition consistent with non-clustered representation of the two parental genomes. These findings indicate potential for successful genetic rescue in future Apteryx reinforcement translocations, a potential that is currently under utilised due to restrictive translocation policies. In light of our findings, we suggest that reconsideration of these policies could render great benefits for the future diversity of this iconic genus in New Zealand.

Evaluating the role of reference-genome phylogenetic distance on evolutionary inference

When a high-quality genome assembly of a target species is unavailable, an option to avoid the costly de novo assembly process is a mapping-based assembly. However, mapping shotgun data to a distant relative may lead to biased or erroneous evolutionary inference. Here, we used short-read data from a mammal (beluga whale) and a bird species (rowi kiwi) to evaluate whether reference genome phylogenetic distance can impact downstream demographic (Pairwise Sequentially Markovian Coalescent) and genetic diversity (heterozygosity, runs of homozygosity) analyses. We mapped to assemblies of species of varying phylogenetic distance (from conspecific to genome-wide divergence of >7%), and de novo assemblies created using cross-species scaffolding. We show that while reference genome phylogenetic distance has an impact on demographic analyses, it is not pronounced until using a reference genome with >3% divergence from the target species. When mapping to cross-species scaffolded assemblies, we are unable to replicate the original beluga demographic results, but are able with the rowi kiwi, presumably reflecting the more fragmented nature of the beluga assemblies. We find that increased phylogenetic distance has a pronounced impact on genetic diversity estimates; heterozygosity estimates deviate incrementally with increasing phylogenetic distance. Moreover, runs of homozygosity are largely undetectable when mapping to any nonconspecific assembly. However, these biases can be reduced when mapping to a cross-species scaffolded assembly. Taken together, our results show that caution should be exercised when selecting reference genomes. Cross-species scaffolding may offer a way to avoid a costly, traditional de novo assembly, while still producing robust, evolutionary inference.

The effects of drift and selection on latitudinal genetic variation in Scandinavian common toads (Bufo bufo) following postglacial recolonisation

Clinal variation is paramount for understanding the factors shaping genetic diversity in space and time. During the last glacial maximum, northern Europe was covered by glacial ice that rendered the region uninhabitable for most taxa. Different evolutionary processes during and after the recolonisation of this area from different glacial refugia have affected the genetic landscape of the present day European flora and fauna. In this study, we focus on the common toad (Bufo bufo) in Sweden and present evidence suggesting that these processes have resulted in two separate lineages of common toad, which colonised Sweden from two directions. Using ddRAD sequencing data for demographic modelling, structure analyses, and analysis of molecular variance (AMOVA), we provide evidence of a contact zone located between Uppland and Västerbotten in central Sweden. Genetic diversity was significantly higher in southern Sweden compared to the north, in accordance with a pattern of decreased genetic diversity with increasing distance from glacial refugia. Candidate genes under putative selection are identified through outlier detection and gene-environment association methods. We provide evidence of divergent selection related to stress response and developmental processes in these candidate genes. The colonisation of Sweden by two separate lineages may have implications for how future conservation efforts should be directed by identifying management units and putative local adaptations.

The eggshell structure in apteryx; form, function, and adaptation

Apteryx is a genus of flightless birds endemic to New Zealand known to lay very large eggs in proportion to body weight. The eggshell of Apteryx is unusually thin and less porous than allometrically expected possibly as a compensation for a very long incubation period. Past studies have been carried out on Apteryx australis, a species which once comprised all kiwi with brown plumage, now separated into three distinct species. These species use different habitats and live at different latitudes and altitudes, therefore generating a need to revise our knowledge of the attributes of their eggshells. In this study, we measured the physical characteristics and water conductance on eggshell fragments of these three species and Great-spotted Kiwi and relate them to the environmental conditions of their respective environments; we also measured the water vapor conductance of Brown Kiwi eggs of late stages of incubation. We found that several trade-offs exist between incubation behavior, environmental conditions, and eggshell structure. We found differences between species in eggshell water vapor conductance seemingly related to altitude; Brown Kiwi and Rowi generally inhabiting lower altitudes had the highest conductance and Tokoeka, generally living in montane environments, the lowest. This is achieved by an increased eggshell thickness rather than a pore area reduction. Finally, the water vapor conductance late in incubation was 58% higher than infertile unincubated eggs, suggesting a drastic increase in conductance throughout the long incubation period. Using the values previously reported, we calculated the embryonic eggshell thinning to be 32.5% at the equatorial region of the eggshell. We describe several new features, such as triangular mineral particles in the cuticle, reported for the extinct Trigonoolithus amoei, and confirmed the existence of plugged pores. We suggest that these structures provide microbial protection needed by a burrow nesting species with a long incubation period.

Detecting Lineage-Specific Shifts in Diversification: A Proper Likelihood Approach

The branching patterns of molecular phylogenies are generally assumed to contain information on rates of the underlying speciation and extinction processes. Simple birth-death models with constant, time-varying, or diversity-dependent rates have been invoked to explain these patterns. They have one assumption in common: all lineages have the same set of diversification rates at a given point in time. It seems likely, however, that there is variability in diversification rates across subclades in a phylogenetic tree. This has inspired the construction of models that allow multiple rate regimes across the phylogeny, with instantaneous shifts between these regimes. Several methods exist for calculating the likelihood of a phylogeny under a specified mapping of diversification regimes and for performing inference on the most likely diversification history that gave rise to a particular phylogenetic tree. Here, we show that the likelihood computation of these methods is not correct. We provide a new framework to compute the likelihood correctly and show, with simulations of a single shift, that the correct likelihood indeed leads to parameter estimates that are on average in much better agreement with the generating parameters than the incorrect likelihood. Moreover, we show that our corrected likelihood can be extended to multiple rate shifts in time-dependent and diversity-dependent models. We argue that identifying shifts in diversification rates is a nontrivial model selection exercise where one has to choose whether shifts in now-extinct lineages are taken into account or not. Hence, our framework also resolves the recent debate on such unobserved shifts. [Diversification; macroevolution; phylogeny; speciation].

Genetic Rescue and the Plight of Ponui Hybrids

Long-term sustainable and resilient populations is a key goal of conservation. How to best achieve this is controversial. There are, for instance, polarized views concerning the fitness and conservation value of hybrid populations founded through multi-origin translocations. A classic example concerns Apteryx (kiwi) in New Zealand. The A. mantelli of Ponui Island constitute a hybrid population where the birds are highly successful in their island habitat. A key dilemma for managers is understanding the reason for this success. Are the hybrid birds of Ponui Island of “no future conservation value” as recently asserted, or do they represent an outstanding example of genetic rescue and an important resource for future translocations? There has been a paradigm shift in scientific thinking concerning hybrids, but the ecological significance of admixed genomes remains difficult to assess. This limits what we can currently predict in conservation science. New understanding from genome science challenges the sufficiency of population genetic models to inform decision making and suggests instead that the contrasting outcomes of hybridization, “outbreeding depression” and “heterosis,” require understanding additional factors that modulate gene and protein expression and how these factors are influenced by the environment. We discuss these findings and the investigations that might help us to better understand the birds of Ponui, inform conservation management of kiwi and provide insight relevant for the future survival of Apteryx.

Dispersal barriers and opportunities drive multiple levels of phylogeographic concordance in the Southern Alps of New Zealand

Phylogeographic concordance, or the sharing of phylogeographic patterns among codistributed species, suggests similar responses to topography or climatic history. While the orientation and timing of breaks between lineages are routinely compared, spatial dynamics within regions occupied by individual lineages provide a second opportunity for comparing responses to past events. In environments with complex topography and glacial history, such as New Zealand's South Island, geographically nested comparisons can identify the processes leading to phylogeographic concordance between and within regional genomic clusters. Here, we used single nucleotide polymorphisms (obtained via ddRADseq) for two codistributed forest beetle species, Agyrtodes labralis (Leiodidae) and Brachynopus scutellaris (Staphylinidae), to evaluate the role of climate change and topography in shaping phylogeographic concordance at two, nested spatial scales: do species diverge over the same geographic barriers, with similar divergence times? And within regions delimited by these breaks, do species share similar spatial dynamics of directional expansion or isolation-by-distance? We found greater congruence of phylogeographic breaks between regions divided by the strongest dispersal barriers (i.e., the Southern Alps). However, these shared breaks were not indicative of shared spatial dynamics within the regions they delimit, and the most similar spatial dynamics between species occurred within regions with the strongest gradients in historical climatic stability. Our results indicate that lack of concordance as traditionally detected by lineage turnover does not rule out the possibility of shared histories, and variation in the presence and type of concordance may provide insights into the different processes shaping phylogeographic patterns across geologically dynamic regions.

Deep Macroevolutionary Impact of Humans on New Zealand's Unique Avifauna

Islands are at the frontline of the anthropogenic extinction crisis [1]. A vast number of island birds have gone extinct since human colonization [2], and an important proportion is currently threatened with extinction [3]. While the number of lost or threatened avian species has often been quantified [4], the macroevolutionary consequences of human impact on island biodiversity have rarely been measured [5]. Here, we estimate the amount of evolutionary time that has been lost or is under threat due to anthropogenic activity in a classic example, New Zealand. Half of its bird taxa have gone extinct since humans arrived [6, 7] and many are threatened [8], including lineages forming highly distinct branches in the avian tree of life [9-11]. Using paleontological and ancient DNA information, we compiled a dated phylogenetic dataset for New Zealand's terrestrial avifauna. We extend the method DAISIE developed for island biogeography [12] to allow for the fact that many of New Zealand's birds are evolutionarily isolated and use it to estimate natural rates of speciation, extinction, and colonization. Simulating under a range of human-induced extinction scenarios, we find that it would take approximately 50 million years (Ma) to recover the number of species lost since human colonization of New Zealand and up to 10 Ma to return to today's species numbers if currently threatened species go extinct. This study puts into macroevolutionary perspective the impact of humans in an isolated fauna and reveals how conservation decisions we take today will have repercussions for millions of years.

Glacial cycles drive rapid divergence of cryptic field vole species

Understanding the factors that contribute to the generation of reproductively isolated forms is a fundamental goal of evolutionary biology. Cryptic species are an especially interesting challenge to study in this context since they lack obvious morphological differentiation that provides clues to adaptive divergence that may drive reproductive isolation. Geographical isolation in refugial areas during glacial cycling is known to be important for generating genetically divergent populations, but its role in the origination of new species is still not fully understood and likely to be situation dependent. We combine analysis of 35,434 single-nucleotide polymorphisms (SNPs) with environmental niche modeling (ENM) to investigate genomic and ecological divergence in three cryptic species formerly classified as the field vole (Microtus agrestis). The SNPs demonstrate high genomic divergence (pairwise F ST values of 0.45-0.72) and little evidence of gene flow among the three field vole cryptic species, and we argue that genetic drift may have been a particularly important mechanism for divergence in the group. The ENM reveals three areas as potential glacial refugia for the cryptic species and differing climatic niches, although with spatial overlap between species pairs. This evidence underscores the role that glacial cycling has in promoting genetic differentiation and reproductive isolation by subdivision into disjunct distributions at glacial maxima in areas relatively close to ice sheets. Future investigation of the intrinsic barriers to gene flow between the field vole cryptic species is required to fully assess the mechanisms that contribute to reproductive isolation. In addition, the Portuguese field vole (M. rozianus) shows a high inbreeding coefficient and a restricted climatic niche, and warrants investigation into its conservation status.

Factors affecting abundance of different stages of the endophilic tick Ixodes anatis in brown kiwi (Apteryx mantelli) shelters

Ixodes anatis is a species of endophilic (nidicolous) tick species parasitizing brown kiwi (Apteryx mantelli). Even though they are endemic to New Zealand like their host, very little is known about these ticks or their population dynamics and relationships with their hosts. We conducted a study from May 2013 to June 2014 to evaluate the effect of shelter location (one of three gullies), habitat (forest, scrub and pasture) and type (tree, soil and surface) on the abundance of the different life stages of I. anatis. In total, 12,172 ticks were collected from 63 shelters, which were sampled monthly for 11 months over the 14 month period. Un-engorged larvae predominated over other stages accounting for 87.2% of the samples collected. We found that location, habitat in which the shelters were located, and the type of shelter were significant predictors of I. anatis abundance. Tree shelters in forests had significantly higher tick abundance than those in scrub and pasture. Tree and soil shelters in general had significantly more ticks than surface shelters. Shelters located in Kauri Bush a drier site, had higher abundances than those in wetter sites. While some of these changes can be explained with the movement of the host, we believe more research needs to be done on the effect of shelters' microclimate on I. anatis' life cycle.

Pleistocene glacial cycles drive isolation, gene flow and speciation in the high-elevation Andes

Mountain ranges are amongst the most species-rich habitats, with many large and rapid evolutionary radiations. The tempo and mode of diversification in these systems are key unanswered questions in evolutionary biology. Here we study the Andean Lupinus radiation to understand the processes driving very rapid diversification in montane systems. We use genomic and transcriptomic data of multiple species and populations, and apply phylogenomic and demographic analyses to test whether diversification proceeded without interspecific gene flow - as expected if Andean orogeny and geographic isolation were the main drivers of diversification - or if diversification was accompanied by gene flow, in which case other processes were probably involved. We uncover several episodes of gene flow between species, including very recent events likely to have been prompted by changes in habitat connectivity during Pleistocene glacial cycles. Furthermore, we find that gene flow between species was heterogeneously distributed across the genome. We argue that exceptionally fast diversification of Andean Lupinus was partly a result of Late Pleistocene glacial cycles, with associated cycles of expansion and contraction driving geographic isolation or secondary contact of species. Furthermore, heterogeneous gene flow across the genome suggests a role for selection and ecological speciation in rapid diversification in this system.

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.

Impact of Model Violations on the Inference of Species Boundaries Under the Multispecies Coalescent

The use of genetic data for identifying species-level lineages across the tree of life has received increasing attention in the field of systematics over the past decade. The multispecies coalescent model provides a framework for understanding the process of lineage divergence and has become widely adopted for delimiting species. However, because these studies lack an explicit assessment of model fit, in many cases, the accuracy of the inferred species boundaries are unknown. This is concerning given the large amount of empirical data and theory that highlight the complexity of the speciation process. Here, we seek to fill this gap by using simulation to characterize the sensitivity of inference under the multispecies coalescent (MSC) to several violations of model assumptions thought to be common in empirical data. We also assess the fit of the MSC model to empirical data in the context of species delimitation. Our results show substantial variation in model fit across data sets. Posterior predictive tests find the poorest model performance in data sets that were hypothesized to be impacted by model violations. We also show that while the inferences assuming the MSC are robust to minor model violations, such inferences can be biased under some biologically plausible scenarios. Taken together, these results suggest that researchers can identify individual data sets in which species delimitation under the MSC is likely to be problematic, thereby highlighting the cases where additional lines of evidence to identify species boundaries are particularly important to collect. Our study supports a growing body of work highlighting the importance of model checking in phylogenetics, and the usefulness of tailoring tests of model fit to assess the reliability of particular inferences. [Populations structure, gene flow, demographic changes, posterior prediction, simulation, genetics.].

Extraction of DNA from captive-sourced feces and molted feathers provides a novel method for conservation management of New Zealand kiwi (Apteryx spp.)

Although some taxa are increasing in number due to active management and predator control, the overall number of kiwi (Apteryx spp.) is declining. Kiwi are cryptic and rare, meaning current monitoring tools, such as call counts, radio telemetry, and surveys using detection dogs are labor-intensive, yield small datasets, and require substantial resources or provide inaccurate estimates of population sizes. A noninvasive genetic approach could help the conservation effort. We optimized a panel of 23 genetic markers (22 autosomal microsatellite loci and an allosomal marker) to discriminate between all species of kiwi and major lineages within species, while simultaneously determining sex. Markers successfully amplified from both fecal and shed feather DNA samples collected in captivity. We found that DNA extraction was more efficient from shed feathers, but DNA quality was greater with feces, although this was sampling dependent. Our microsatellite panel was able to distinguish between contemporary kiwi populations and lineages and provided PI values in the range of 4.3 × 10-5 to 2.0 × 10-19, which in some cases were sufficient for individualization and mark-recapture studies. As such, we have tested a wide-reaching, noninvasive molecular approach that will improve conservation management by providing better parameter estimates associated with population ecology and demographics such as abundance, growth rates, and genetic diversity.

Impact of the HIV-1 genetic background and HIV-1 population size on the evolution of raltegravir resistance

BACKGROUND

Emergence of resistance against integrase inhibitor raltegravir in human immunodeficiency virus type 1 (HIV-1) patients is generally associated with selection of one of three signature mutations: Y143C/R, Q148K/H/R or N155H, representing three distinct resistance pathways. The mechanisms that drive selection of a specific pathway are still poorly understood. We investigated the impact of the HIV-1 genetic background and population dynamics on the emergence of raltegravir resistance. Using deep sequencing we analyzed the integrase coding sequence (CDS) in longitudinal samples from five patients who initiated raltegravir plus optimized background therapy at viral loads > 5000 copies/ml. To investigate the role of the HIV-1 genetic background we created recombinant viruses containing the viral integrase coding region from pre-raltegravir samples from two patients in whom raltegravir resistance developed through different pathways. The in vitro selections performed with these recombinant viruses were designed to mimic natural population bottlenecks.

RESULTS

Deep sequencing analysis of the viral integrase CDS revealed that the virological response to raltegravir containing therapy inversely correlated with the relative amount of unique sequence variants that emerged suggesting diversifying selection during drug pressure. In 4/5 patients multiple signature mutations representing different resistance pathways were observed. Interestingly, the resistant population can consist of a single resistant variant that completely dominates the population but also of multiple variants from different resistance pathways that coexist in the viral population. We also found evidence for increased diversification after stronger bottlenecks. In vitro selections with low viral titers, mimicking population bottlenecks, revealed that both recombinant viruses and HXB2 reference virus were able to select mutations from different resistance pathways, although typically only one resistance pathway emerged in each individual culture.

CONCLUSIONS

The generation of a specific raltegravir resistant variant is not predisposed in the genetic background of the viral integrase CDS. Typically, in the early phases of therapy failure the sequence space is explored and multiple resistance pathways emerge and then compete for dominance which frequently results in a switch of the dominant population over time towards the fittest variant or even multiple variants of similar fitness that can coexist in the viral population.

Strong population structure in a species manipulated by humans since the Neolithic: the European fallow deer (Dama dama dama)

Species that have been translocated and otherwise manipulated by humans may show patterns of population structure that reflect those interactions. At the same time, natural processes shape populations, including behavioural characteristics like dispersal potential and breeding system. In Europe, a key factor is the geography and history of climate change through the Pleistocene. During glacial maxima throughout that period, species in Europe with temperate distributions were forced south, becoming distributed among the isolated peninsulas represented by Anatolia, Italy and Iberia. Understanding modern patterns of diversity depends on understanding these historical population dynamics. Traditionally, European fallow deer (Dama dama dama) are thought to have been restricted to refugia in Anatolia and possibly Sicily and the Balkans. However, the distribution of this species was also greatly influenced by human-mediated translocations. We focus on fallow deer to better understand the relative influence of these natural and anthropogenic processes. We compared modern fallow deer putative populations across a broad geographic range using microsatellite and mitochondrial DNA loci. The results revealed highly insular populations, depauperate of genetic variation and significantly differentiated from each other. This is consistent with the expectations of drift acting on populations founded by small numbers of individuals, and reflects known founder populations in the north. However, there was also evidence for differentiation among (but not within) physically isolated regions in the south, including Iberia. In those regions we find evidence for a stronger influence from natural processes than may be expected for a species with such strong, known anthropogenic influence.