Population genetics of Capercaillie (Tetrao urogallus) in the Jura and the Pyrenees: a non-invasive approach to avian conservation genetics

RÉSUMÉ Le Grand tétras est un galliforme de montagne apparenté au faisan et au tétras lyre. Il est distribué de manière continue à travers la toundra et les montagnes de moyenne altitude en Europe de l'ouest. Toutefois, les populations d'Europe de l'ouest ont subi un déclin constant au cours des derniers siècles. Les causes de ce déclin sont probablement liées à l'activité humaine, telle .que l'élevage ou le tourisme, qui ont engendré une modification et une fragmentation de l'habitat de l'espèce. Malheureusement, les populations soumises à de forts déclins démographiques peuvent subir des effets génétiques (augmentation de la consanguinité et perte de diversité génétique) pouvant diminuer leur potentiel de reproduction et conduire irrémédiablement à l'extinction. Cette thèse présente les analyses conduites dans le but d'estimer l'impact du déclin démographique des populations de Grand tétras sur l'étendue et la distribution de leur variabilité génétique dans le Jura et dans les Pyrénées. Du fait de la législation locale protégeant les tétraonidés en général, mais également en raison de la biologie très cryptique du Grand tétras, l'ensemble des analyses de cette étude a été réalisé à partir de matériel génétique extrait des fientes (ou échantillonnage génétique non invasif). Dans la première partie de l'étude, je détaille les protocoles d'extraction. d'ADN et d'amplification par PCR modifiés à partir des protocoles classiques utilisant des échantillons conventionnels, riches en ADN. L'utilisation d'ADN fécal impose des contraintes dues à la mauvaise qualité et à la faible quantité du matériel génétique à disposition dans les fientes. Ces contraintes ont pu être partiellement contournées en réalisant des répétitions multiples du génotypage afin d'obtenir un degré de fiabilité suffisante. J'ai également analysé les causes de la dégradation de l'ADN dans les excréments. Parmi les causes les plus communes, telles que l'activité bactérienne, l'hydrolyse spontanée et la dégradation enzymatique par les DNases libres, c'est ce dernier facteur qui apparaît comme étant la cause majeure et la plus rapide responsable de la dégradation de la qualité des échantillons. La rapidité de l'action enzymatique suggère que les plans d'échantillonnages de excréments sur le terrain pourraient être optimisés en les réalisant dans des conditions climatiques froides et sèches, favorisant ainsi l'inhibition des DNases. La seconde partie de la thèse est une étude par simulation visant à déterminer la capacité du logiciel Structure à identifier les structures génétiques complexes et hiérarchiques fréquemment rencontrées dans les populations naturelles, et ce en utilisant différents types de marqueurs génétiques. Les troisième et quatrième parties de cette thèse décrivent le statut génétique des populations résiduelles du Jura et des Pyrénées à partir de l'analyse de 11 loci microsatellites. Nous n'avons pas pu mettre en évidence dans les deux populations des effets liés à la consanguinité ou à la réduction de la diversité génétique. De plus, la différenciation génétique entre les patches d'habitats favorables reste modérée et corrélée à la distance géographique, ce qui suggère que la dispersion d'individus entre les patches a été importante au moins pendant ces dernières générations. La comparaison des paramètres de la diversité génétique avec ceux d'autres populations de Grand tétras, ou d'autres espèces proches, indique que la population du Jura a retenu une proportion importante de sa diversité originelle. Ces résultats suggèrent que le déclin récent des populations a jusqu'ici eu un impact modéré sur les facteurs génétiques et que ces populations semblent avoir conservé le potentiel génétique nécessaire à leur survie à long terme. Finalement, en cinquième partie, l'analyse de l'apparentement entre les mâles qui participent à la parade sur les places de chant (leks) indique que ces derniers sont distribués en agrégats de manière non aléatoire, préférentiellement entre individus apparentés. De plus, la corrélation entre les distances génétique et géographique entre les leks est en accord avec les motifs d'isolement par la distance mis en évidence à d'autres niveaux hiérarchiques (entre patches d'habitat et populations), ainsi qu'avec les études menées sur d'autres espèces ayant choisi ce même système de reproduction. En conclusion, cette première étude basée uniquement sur de l'ADN nucléaire aviaire extrait à partir de fèces a fourni des informations nouvelles qui n'auraient pas pu être obtenues par une méthode d'observation sur le terrain ou d'échantillonnage génétique classique. Aucun oiseau n'a été dérangé ou capturé, et les résultats sont comparables à d'autres études concernant des espèces proches. Néanmoins, la taille de ces populations approche des niveaux au-dessous desquels la survie à long terme est fortement incertaine. La persistance de la diversité génétique pour les prochaines générations reste en conséquence liée à la survie des adultes et à une reprise du succès de la reproduction. ABSTRACT Capercaillie (Tetrao urogallus) is a large grouse that is continuously distributed across the tundra and the mid-high mountains of Western Europe. However, the populations in Western Europe have been showing a constant decline during the last decades. The causes for this decline are possibly related to human activities, such as cattle breeding and tourism that have both led to habitat modification and fragmentation. Unfortunately, populations that have undergone drastic demographic bottlenecks often go through genetic processes of inbreeding and loss of diversity that decrease their fitness and eventually lead to extinction. This thesis presents the investigations conducted to estimate the impact of the demographic decline of capercaillie populations on the extent and distribution of their genetic variability in the Jura and in the Pyrenees mountains. Because grouse are protected by wildlife legislation, and also because of the cryptic behaviour of capercaillie, all DNA material used in this study was extracted from faeces (non-invasive genetic sampling). In the first part of my thesis, I detail the protocols of DNA extraction and PCR amplification adapted from classical methods using conventional DNA-rich samples. The use of faecal DNA imposes specific constraints due to the low quantity and the highly degraded genetic material available. These constraints are partially overcome by performing multiple genotyping repetitions to obtain sufficient reliability. I also investigate the causes of DNA degradation in faeces. Among the main degraders, namely bacterial activity, spontaneous hydrolysis, and free-¬DNase activities, the latter was pointed out as the most important according to our experiments. These enzymes degrade DNA very rapidly, and, as a consequence, faeces sampling schemes must be planned preferably in cold and dry weather conditions, allowing for enzyme activity inhibition. The second part of the thesis is a simulation study aiming to assess the capacity of the software Structure to detect population structure in hierarchical models relevant to situations encountered in wild populations, using several genetic markers. The methods implemented in Structure appear efficient in detecting the highest hierarchical structure. The third and fourth parts of the thesis describe the population genetics status of the remaining Jura and Pyrenees populations using 11 microsatellite loci. In either of these populations, no inbreeding nor reduced genetic diversity was detected. Furthermore, the genetic differentiation between patches defined by habitat suitability remains moderate and correlated with geographical distance, suggesting that significant dispersion between patches was at work at least until the last generations. The comparison of diversity indicators with other species or other populations of capercaillie indicate that population in the Jura has retained a large part of its original genetic diversity. These results suggest that the recent decline has had so forth a moderate impact on• genetic factors and that these populations might have retained the potential for long term survival, if the decline is stopped. Finally, in the fifth part, the analysis of relatedness between males participating in the reproduction parade, or lek, indicate that capercaillie males, like has been shown for some other grouse species, gather on leks• among individuals that are more related than the average of the population. This pattern appears to be due to both population structure and kin-association. As a conclusion, this first study relying exclusively on nuclear DNA extracted from faeces has provided novel information that was not available through field observation or classical genetic sampling. No bird has been captured or disturbed, and the results are consistent with other studies of closely related species. However, the size of these populations is approaching thresholds below which long-term survival is unlikely. The persistence of genetic diversity for the forthcoming generations remains therefore bond to adult survival and to the increase of reproduction success.

Phylogéographie des principales vipères européennes (Vipera ammodytes, V. aspis et V. berus), structuration génétique et multipaternité chez Vipera berus

Résumé: Les vipères du genre Vipera sont des serpents venimeux distribués dans la totalité du Paléarctique. Malgré cette répartition considérable, elles sont extrêmement menacées, leur déclin étant principalement dû à la destruction et à la fragmentation de leur habitat ainsi qu'à la persécution humaine. Afin d'apporter de nouveaux éléments dans le contexte de la protection de ce groupe de reptiles, nous avons utilisé durant ce travail de thèse différents marqueurs moléculaires pour étudier la structuration génétique à petite et à large échelle chez trois espèces appartenant au genre Vipera. La première étude, une phylogéographie moléculaire de la vipère ammodytes (Vipera ammodytes), a montré dans l'ensemble de l'aire de répartition une forte structuration génétique provenant d'isolements antérieures au Pléistocène. La présence d'un nombre important de clades dans le centre des Balkans suggère que cette région a fourni de nombreux refuges isolés durant les glaciations. Ces dernières ont également eu un impact considérable sur la diversité génétique au sein de la majorité des clades, suite à d'importants goulots d'étranglement durant le Pléistocène. L'étude de la phylogéographie de la vipère aspic (Vipera aspis) a montré une différenciation génétique entre les populations présentes de chaque côté des Alpes, mais également une forte structuration interne avec la mise en évidence d'un refuge en France. Cette étude est la première à établir clairement l'utilisation d'un refuge français pour un vertébré terrestre. La troisième partie de cette thèse a étudié la phylogéographie de la vipère péliade (Vipera berus), espèce cible de ce travail. En plus de la mise en évidence d'un groupe génétique inattendu (localisé dans le nord de l'Italie, le sud de l'Autriche, le nord de la Slovénie et l'extrême sud-est de la Suisse), la variabilité génétique au sein du groupe nordique (comprenant les animaux de l'entier de l'aire de répartition de l'espèce à l'exception des individus du groupe italien et les animaux provenant des Balkans) est suffisamment importante pour conclure à l'utilisation de refuges glaciaires nordiques durant les dernières glaciations, en complément des refuges habituellement décrits pour la majorité des espèces animales (soit les péninsules ibérique, italienne et balakanique). Ces résultats nous ont conduit à effectuer une étude morphologique (quatrième partie) comparant les vipères péliades du "clade italien" et du "clade nordique" décrits ci-dessus. Seules de petites différences morphologiques ont pu être mises en évidence, malgré une séparation de ces groupes estimée à plus d'un million d'années. Une étude à plus petite échelle, centrée sur le Massif jurassien et certaines populations alpines et françaises, a été entreprise afin d'estimer leur diversité génétique et d'évaluer la structuration génétique entre les populations à l'aide de marqueurs microsatellites (cinquième partie). Une importante structuration a été observée entre les populations distantes de plus de 3 kilomètres, la structuration entre les populations plus proches étant plus limitée. De plus, une diversité génétique plus faible dans les populations jurassiennes et alpines comparativement aux populations du massif central et de la côte atlantique a été constatée, probablement due à une perte de diversité génétique lors de la recolonisation post-glaciaire. La sixième étude s'est intéressée au succès reproducteur des mâles de vipères péliades en conditions naturelles. Une corrélation entre la taille des mâles et leur succès reproducteur a été relevée, les individus de plus grande taille ayant un succès reproducteur plus élevé. Le taux de multipaternité a aussi été investigué, démontrant que la proportion de pontes issues de plusieurs pères est élevée (69%) malgré la faible densité de vipères observée sur le site étudié. Finalement, aucun lien entre le nombre de pères au sein d'une ponte et la mortalité des jeunes à la naissance n'a pu être mis en évidence, contrastant avec des travaux précédents. En conclusion, l'observation de la structuration très marquée chez les vipères péliades devrait permettre d'affiner les méthodes de protection de l'espèce dans le massif jurassien. A plus large échelle, l'importante structuration génétique observée chez les vipères ammodytes, aspic et péliade résultant de l'utilisation de nombreux refuges glaciaires, complémentaires aux refuges habituellement utilisés par les espèces animales, démontre l'intérêt de l'analyse phylogéographique des reptiles pour la compréhension des phénomènes de colonisation et d' extinction des populations durant la fin du Tertiaire et le Quaternaire. La mise en évidence chez les différentes espèces de vipères étudiées de nombreux groupes génétiques distincts (ESUs) devrait conduire à des modifications de la taxonomie ainsi qu'au statut de protection de ces espèces. Abstract: The vipers of the genus Vipera are venomous snakes widespread throughout the Palaearctic regions. Despite a large distribution area, several species are extremely threatened, especially due to the destruction and fragmentation of their habitats, as well as by human persecution. In order to increase the knowledge on these species and to improve their protection, several molecular markers have been used to investigate the genetic structure on small and large scales, within three species of the genus Vipera. The first study, a molecular phylogeography of the nose-horned viper (Vipera ammodytes), showed a considerable structuring throughout the distribution area, due to isolation into refugia before the Pleistocene. A high number of clades in the centre of the Balkans suggests that this region harboured numerous isolated glacial refugia during the last glaciation. Moreover, low genetic diversity within several clades implies that most populations of nose-horned vipers have suffered bottlenecks during the Pleistocene. The study of the phylogeography of the asp viper (Vipera aspis) showed genetic differentiation between populations on each side of the Alps, as well as considerable internal genetic structure, suggesting the use of a glacial refugium in France. This study is the first to establish firmly the occurrence of a French refugia for a terrestrial vertebrate. The third part of this work involved a phylogeographic study of the adder (Vipera berus), the target species of this thesis. Three clades were revealed: a Balkan clade (corresponding to the subspecies V. b. sachalinensis), an unexpected Italian clade (limited to northern Italy, southern Austria, northern Slovenia and southeasternmost corner of Switzerland) and a Northern clade clade (including adders of the whole distribution area excepted animals from the Balkan and the Italian clades). The genetic variability within the Northern clade is sufficiently high to conclude that a northern glacial refugia during the last glaciation, in addition to those refugia already described for the main species (Iberian, Italian and Balkan peninsula). These results motivated a morphological study (part four) comparing the adders from the Italian and the Northern clades describe above. Only small morphological differences have been found, despite the split between these two clades have taken place more than 1 million years ago. A study on a local scale, focused on the Jura Mountains, on a few populations in the Alps and France was, performed to estimate the genetic diversity and the genetic structure between populations using microsatellite markers (part five). Considerable structure was observed between populations separated by more than 3 kilometres, whereas the structure between closer populations is less marked. Moreover, lower genetic diversity in the populations from Jura Mountains and Alps was noticed compared to populations from Massif Central of Atlantic coast. Such loss of genetic variation probably followed post-glacial recolonisation. The sixth study focused on the reproductive success of male adders in the wild. A positive correlation between body length and reproductive success was observed. Multiple paternity was also observed in most of clutches (69%) despite the low density of adders in the study area. Finally, no relationship was found between the number of fathers in a clutch and the survival of offspring at birth, contradicting previous studies. To conclude, the observation of a significant genetic structure in Vipera berus will enable recommendations to be made to improve protection of this species in the Jura Mountain. On a larger scale, the considerable genetic structure found within Vipera ammdoytes, V. aspis and V. berus, resulting from isolation in additional glacial refugia to those already described for other species, demonstrates the relevance of phylogeographic studies of reptiles to better understand the colonisation and disappearance during the last Tertiary and the Quaternary. The observation of several groups of evolutionary significant units (ESUs) within the three studied species might lead to a revision of the taxonomy, as well as their conservation status.

Live trapping design for the harvest mouse (Micromys minutus) in its summer habitat

documented accurately since 1960. Most records are based on nest findings and there have been few direct observations or captures, mainly because live trapping of this species is not simple. Therefore, an efficient trapping technique is needed for population studies and to facilitate the management of its habitat. By combining the methods used to capture very small (Suncus etruscus) and climbing (Muscardinus avellanarius) mammals, we developed a design using Longworth traps with mouse excluders set on suspended platforms. This allowed us to trap more harvest mice in four field sessions of 60 trap-nights than have ever been caught previously since its discovery in Switzerland.

EndoG links Bnip3-induced mitochondrial damage and caspase-independent DNA fragmentation in ischemic cardiomyocytes

Mitochondrial dysfunction, caspase activation and caspase-dependent DNA fragmentation are involved in cell damage in many tissues. However, differentiated cardiomyocytes repress the expression of the canonical apoptotic pathway and their death during ischemia is caspase-independent. The atypical BH3-only protein Bnip3 is involved in the process leading to caspase-independent DNA fragmentation in cardiomyocytes. However, the pathway by which DNA degradation ensues following Bnip3 activation is not resolved. To identify the mechanism involved, we analyzed the interdependence of Bnip3, Nix and EndoG in mitochondrial damage and DNA fragmentation during experimental ischemia in neonatal rat ventricular cardiomyocytes. Our results show that the expression of EndoG and Bnip3 increases in the heart throughout development, while the caspase-dependent machinery is silenced. TUNEL-positive DNA damage, which depends on caspase activity in other cells, is caspase-independent in ischemic cardiomyocytes and ischemia-induced DNA high and low molecular weight fragmentation is blocked by repressing EndoG expression. Ischemia-induced EndoG translocation and DNA degradation are prevented by silencing the expression of Bnip3, but not Nix, or by overexpressing Bcl-xL. These data establish a link between Bnip3 and EndoG-dependent, TUNEL-positive, DNA fragmentation in ischemic cardiomyocytes in the absence of caspases, defining an alternative cell death pathway in postmitotic cells.

Sustained sleep fragmentation induces sleep homeostasis in mice.

STUDY OBJECTIVES: Sleep fragmentation (SF) is an integral feature of sleep apnea and other prevalent sleep disorders. Although the effect of repetitive arousals on cognitive performance is well documented, the effects of long-term SF on electroencephalography (EEG) and molecular markers of sleep homeostasis remain poorly investigated. To address this question, we developed a mouse model of chronic SF and characterized its effect on EEG spectral frequencies and the expression of genes previously linked to sleep homeostasis including clock genes, heat shock proteins, and plasticity-related genes. DESIGN: N/A. SETTING: Animal sleep research laboratory. PARTICIPANTS: Sixty-six C57BL6/J adult mice. INTERVENTIONS: Instrumental sleep disruption at a rate of 60/h during 14 days. MEASUREMENTS AND RESULTS: Locomotor activity and EEG were recorded during 14 days of SF followed by recovery for 2 days. Despite a dramatic number of arousals and decreased sleep bout duration, SF minimally reduced total quantity of sleep and did not significantly alter its circadian distribution. Spectral analysis during SF revealed a homeostatic drive for slow wave activity (SWA; 1-4 Hz) and other frequencies as well (4-40 Hz). Recordings during recovery revealed slow wave sleep consolidation and a transient rebound in SWA, and paradoxical sleep duration. The expression of selected genes was not induced following chronic SF. CONCLUSIONS: Chronic SF increased sleep pressure confirming that altered quality with preserved quantity triggers core sleep homeostasis mechanisms. However, it did not induce the expression of genes induced by sleep loss, suggesting that these molecular pathways are not sustainably activated in chronic diseases involving SF.

The dynamics of coexistence: habitat sharing versus segregation patterns among three sympatric montane vipers

Contact zones of closely related and ecologically similar species constitute rare opportunities to study the evolutionary consequences of past speciation processes. They represent natural laboratories in which strong competition could lead to the exclusion of one species, or the various species may switch into distinct ecological niches. Alternatively, if reproductive isolation has not yet been achieved, they may hybridize. We elucidate the degree of taxon integrity by comparing genetics and habitat use of three similar-sized congeneric viper species, Vipera ammodytes, Viperaaspis, and Viperaberus, of Nadiza Valley in western Slovenia. No hybridization was detected for either mitochondrial or nuclear genomes. Similarly, external intermediacy by a single prestudy viper (probably V.ammodytesxV. aspis) indicates that hybridization occasionally occurs, but should be very rare. Populations of the three related viperids are partially allopatric in Nadiza Valley, but they also coexist in a narrow contact zone in the montane grassland along the south-exposed slope of Mount Stol (1673m a.s.l.). Here, the three species that occupy areas in or near patches of rocky microhabitats (e.g. stone piles, slides, and walls) live in syntopy. However, fine-scale measurements of structural components show partial habitat segregation, in which V.berus becomes more dominant at elevations above 1400m and occupies mostly the mountain ridge and north-exposed slopes of Mount Stol, V.aspis occurs below 1300m and is the only species to inhabit stoneless patches of grass and bushes around 1000m and lower, and V.ammodytes occurs at all elevations up to 1500m, but is restricted to a rocky microhabitat. We suggest that a high degree of microstructure divergence, slightly different environmental niches, and a generally favourable habitat for all three viper species, keep the pressure for mis-mating and hybridization low, although mechanisms such as reduced hybrid inferiority and temporal mating segregation cannot yet be excluded.

Stable isotopes unveil habitat partitioning among the marine mammals off NW Africa and reveal unique trophic niches for two globally threatened species.

Stable isotope abundances of carbon (δ13C) and nitrogen (δ15N) in the bone of 13 species of marine mammals from the northwest coast of Africa were investigated to assess their positions in the local trophic web and their preferred habitats. Also, samples of primary producers and potential prey species from the study area were collected to characterise the local isotopic landscape. This characterisation indicated that δ13C values increased from offshore to nearshore and that δ15N was a good proxy for trophic level. Therefore, the most coastal species were Monachus monachus and Sousa teuszii, whereas the most pelagic were Physeter macrocephalus and Balaenoptera acutorostrata. δ15N values indicated that marine mammals located at the lowest trophic level were B. acutorostrata, Stenella coeruleoalba and Delphinus sp., and those occupying the highest trophic level were M. monachus and P. macrocephalus. The trophic level of Orcinus orca was similar to that of M. monachus, suggesting that O. orca preys on fish. Conservation of coastal and threatened species (M. monachus and S. teuszii) off NW Africa should be a priority because these species, as the main apex predators, cannot be replaced by other marine mammals.

The Armc10/SVH gene: Genome context, regulation of mitochondrial dynamics and protection against Aβ-induced mitochondrial fragmentation

Mitochondrial function and dynamics are essential for neurotransmission, neural function and neuronal viability. Recently, we showed that the eutherian-specific Armcx gene cluster (Armcx1-6 genes), located in the X chromosome, encodes for a new family of proteins that localise to mitochondria, regulating mitochondrial trafficking. The Armcx gene cluster evolved by retrotransposition of the Armc10 gene mRNA, which is present in all vertebrates and is considered to be the ancestor gene. Here we investigate the genomic organisation, mitochondrial functions and putative neuroprotective role of the Armc10 ancestor gene. The genomic context of the Armc10 locus shows considerable syntenic conservation among vertebrates, and sequence comparisons and CHIP-data suggest the presence of at least three conserved enhancers. We also show that the Armc10 protein localises to mitochondria and that it is highly expressed in the brain. Furthermore, we show that Armc10 levels regulate mitochondrial trafficking in neurons, but not mitochondrial aggregation, by controlling the number of moving mitochondria. We further demonstrate that the Armc10 protein interacts with the KIF5/Miro1-2/Trak2 trafficking complex. Finally, we show that overexpression of Armc10 in neurons prevents A beta-induced mitochondrial fission and neuronal death. Our data suggest both conserved and differential roles of the Armc10/Armcx gene family in regulating mitochondrial dynamics in neurons, and underscore a protective effect of the Armc10 gene against A beta-induced toxicity. Overall, our findings support a further degree of regulation of mitochondrial dynamics in the brain of more evolved mammals.

The Armc10/SVH gene: Genome context, regulation of mitochondrial dynamics and protection against Aβ-induced mitochondrial fragmentation

Mitochondrial function and dynamics are essential for neurotransmission, neural function and neuronal viability. Recently, we showed that the eutherian-specific Armcx gene cluster (Armcx1-6 genes), located in the X chromosome, encodes for a new family of proteins that localise to mitochondria, regulating mitochondrial trafficking. The Armcx gene cluster evolved by retrotransposition of the Armc10 gene mRNA, which is present in all vertebrates and is considered to be the ancestor gene. Here we investigate the genomic organisation, mitochondrial functions and putative neuroprotective role of the Armc10 ancestor gene. The genomic context of the Armc10 locus shows considerable syntenic conservation among vertebrates, and sequence comparisons and CHIP-data suggest the presence of at least three conserved enhancers. We also show that the Armc10 protein localises to mitochondria and that it is highly expressed in the brain. Furthermore, we show that Armc10 levels regulate mitochondrial trafficking in neurons, but not mitochondrial aggregation, by controlling the number of moving mitochondria. We further demonstrate that the Armc10 protein interacts with the KIF5/Miro1-2/Trak2 trafficking complex. Finally, we show that overexpression of Armc10 in neurons prevents A beta-induced mitochondrial fission and neuronal death. Our data suggest both conserved and differential roles of the Armc10/Armcx gene family in regulating mitochondrial dynamics in neurons, and underscore a protective effect of the Armc10 gene against A beta-induced toxicity. Overall, our findings support a further degree of regulation of mitochondrial dynamics in the brain of more evolved mammals.

Fragmentation of contaminant and endogenous DNA in ancient samples determined by shotgun sequencing; prospects for human palaeogenomics

Despite the successful retrieval of genomes from past remains, the prospects for human palaeogenomics remain unclear because of the difficulty of distinguishing contaminant from endogenous DNA sequences. Previous sequence data generated on high-throughput sequencing platforms indicate that fragmentation of ancient DNA sequences is a characteristic trait primarily arising due to depurination processes that create abasic sites leading to DNA breaks.