Demographic History of the Genus Pan Inferred from Whole Mitochondrial Genome Reconstructions

Population dynamics and genetic connectivity in recent chimpanzee history

Knowledge on the population history of endangered species is critical for conservation, but whole-genome data on chimpanzees (Pan troglodytes) is geographically sparse. Here, we produced the first non-invasive geolocalized catalog of genomic diversity by capturing chromosome 21 from 828 non-invasive samples collected at 48 sampling sites across Africa. The four recognized subspecies show clear genetic differentiation correlating with known barriers, while previously undescribed genetic exchange suggests that these have been permeable on a local scale. We obtained a detailed reconstruction of population stratification and fine-scale patterns of isolation, migration, and connectivity, including a comprehensive picture of admixture with bonobos (Pan paniscus). Unlike humans, chimpanzees did not experience extended episodes of long-distance migrations, which might have limited cultural transmission. Finally, based on local rare variation, we implement a fine-grained geolocalization approach demonstrating improved precision in determining the origin of confiscated chimpanzees.

Evolutionary history of the extinct Sardinian dhole

The Sardinian dhole (Cynotherium sardous)¹ was an iconic and unique canid species that was endemic to Sardinia and Corsica until it became extinct at the end of the Late Pleistocene.^2-5 Given its peculiar dental morphology, small body size, and high level of endemism, several extant canids have been proposed as possible relatives of the Sardinian dhole, including the Asian dhole and African hunting dog ancestor.^3,6-9 Morphometric analyses^3,6,8-12 have failed to clarify the evolutionary relationship with other canids.We sequenced the genome of a ca-21,100-year-old Sardinian dhole in order to understand its genomic history and clarify its phylogenetic position. We found that it represents a separate taxon from all other living canids from Eurasia, Africa, and North America, and that the Sardinian dhole lineage diverged from the Asian dhole ca 885 ka. We additionally detected historical gene flow between the Sardinian and Asian dhole lineages, which ended approximately 500-300 ka, when the land bridge between Sardinia and mainland Italy was already broken, severing their population connectivity. Our sample showed low genome-wide diversity compared to other extant canids-probably a result of the long-term isolation-that could have contributed to the subsequent extinction of the Sardinian dhole.

Characterization of nuclear mitochondrial insertions in the whole genomes of primates

The transfer and integration of whole and partial mitochondrial genomes into the nuclear genomes of eukaryotes is an ongoing process that has facilitated the transfer of genes and contributed to the evolution of various cellular pathways. Many previous studies have explored the impact of these insertions, referred to as NumtS, but have focused primarily on older events that have become fixed and are therefore present in all individual genomes for a given species. We previously developed an approach to identify novel Numt polymorphisms from next-generation sequence data and applied it to thousands of human genomes. Here, we extend this analysis to 79 individuals of other great ape species including chimpanzee, bonobo, gorilla, orang-utan and also an old world monkey, macaque. We show that recent Numt insertions are prevalent in each species though at different apparent rates, with chimpanzees exhibiting a significant increase in both polymorphic and fixed Numt sequences as compared to other great apes. We further assessed positional effects in each species in terms of evolutionary time and rate of insertion and identified putative hotspots on chromosome 5 for Numt integration, providing insight into both recent polymorphic and older fixed reference NumtS in great apes in comparison to human events.

Newly discovered cichlid fish biodiversity threatened by hybridization with non-native species

Invasive freshwater fishes are known to readily hybridize with indigenous congeneric species, driving loss of unique and irreplaceable genetic resources. Here we reveal that newly discovered (2013-2016) evolutionarily significant populations of Korogwe tilapia (Oreochromis korogwe) from southern Tanzania are threatened by hybridization with the larger invasive Nile tilapia (Oreochromis niloticus). We use a combination of morphology, microsatellite allele frequencies and whole genome sequences to show that O. korogwe from southern lakes (Nambawala, Rutamba and Mitupa) are distinct from geographically disjunct populations in northern Tanzania (Zigi River and Mlingano Dam). We also provide genetic evidence of O. korogwe × niloticus hybrids in three southern lakes and demonstrate heterogeneity in the extent of admixture across the genome. Finally, using the least admixed genomic regions we estimate that the northern and southern O. korogwe populations most plausibly diverged ~140,000 years ago, suggesting that the geographical separation of the northern and southern groups is not a result of a recent translocation, and instead these populations represent independent evolutionarily significant units. We conclude that these newly discovered and phenotypically unique cichlid populations are already threatened by hybridization with an invasive species, and propose that these irreplaceable genetic resources would benefit from conservation interventions.

Similar patterns of genetic diversity and linkage disequilibrium in Western chimpanzees (Pan troglodytes verus) and humans indicate highly conserved mechanisms of MHC molecular evolution

BACKGROUND

Many species are threatened with extinction as their population sizes decrease with changing environments or face novel pathogenic threats. A reduction of genetic diversity at major histocompatibility complex (MHC) genes may have dramatic effects on populations' survival, as these genes play a key role in adaptive immunity. This might be the case for chimpanzees, the MHC genes of which reveal signatures of an ancient selective sweep likely due to a viral epidemic that reduced their population size a few million years ago. To better assess how this past event affected MHC variation in chimpanzees compared to humans, we analysed several indexes of genetic diversity and linkage disequilibrium across seven MHC genes on four cohorts of chimpanzees and we compared them to those estimated at orthologous HLA genes in a large set of human populations.

RESULTS

Interestingly, the analyses uncovered similar patterns of both molecular diversity and linkage disequilibrium across the seven MHC genes in chimpanzees and humans. Indeed, in both species the greatest allelic richness and heterozygosity were found at loci A, B, C and DRB1, the greatest nucleotide diversity at loci DRB1, DQA1 and DQB1, and both significant global linkage disequilibrium and the greatest proportions of haplotypes in linkage disequilibrium were observed at pairs DQA1 ~ DQB1, DQA1 ~ DRB1, DQB1 ~ DRB1 and B ~ C. Our results also showed that, despite some differences among loci, the levels of genetic diversity and linkage disequilibrium observed in contemporary chimpanzees were globally similar to those estimated in small isolated human populations, in contrast to significant differences compared to large populations.

CONCLUSIONS

We conclude, first, that highly conserved mechanisms shaped the diversity of orthologous MHC genes in chimpanzees and humans. Furthermore, our findings support the hypothesis that an ancient demographic decline affecting the chimpanzee populations - like that ascribed to a viral epidemic - exerted a substantial effect on the molecular diversity of their MHC genes, albeit not more pronounced than that experienced by HLA genes in human populations that underwent rapid genetic drift during humans' peopling history. We thus propose a model where chimpanzees' MHC genes regenerated molecular variation through recombination/gene conversion and/or balancing selection after the selective sweep.

Sex-biased patterns shaped the genetic history of Roma

The Roma population is a European ethnic minority characterized by recent and multiple dispersals and founder effects. After their origin in South Asia around 1,500 years ago, they migrated West. In Europe, they diverged into ethnolinguistically distinct migrant groups that spread across the continent. Previous genetic studies based on genome-wide data and uniparental markers detected Roma founder events and West-Eurasian gene flow. However, to the best of our knowledge, it has not been assessed whether these demographic processes have equally affected both sexes in the population. The present study uses the largest and most comprehensive dataset of complete mitochondrial and Y chromosome Roma sequences to unravel the sex-biased patterns that have shaped their genetic history. The results show that the Roma maternal genetic pool carries a higher lineage diversity from South Asia, as opposed to a single paternal South Asian lineage. Nonetheless, the European gene flow events mainly occurred through the maternal lineages; however, a signal of this gene flow is also traceable in the paternal lineages. We also detect a higher female migration rate among European Roma groups. Altogether, these results suggest that sociocultural factors influenced the emergence of sex-biased genetic patterns at global and local scales in the Roma population through time.

The Relationship Between Tool Use and Prey Availability in Chimpanzees (Pan troglodytes schweinfurthii) of Northern Democratic Republic of Congo

A key feature of human behavioral diversity is that it can be constrained by cultural preference (“cultural override”); that is, population-specific preferences can override resource availability. Here we investigate whether a similar phenomenon can be found in one of our closest relatives, as well as the potential impacts of ecological differences on feeding behavior. Our study subjects were different subpopulations of Eastern chimpanzees (Pan troglodytes schweinfurthii) occupying two very different habitats, moist tropical lowland forests vs. moist tropical forest–savanna mosaic on opposite sides of a major river. Given differences in encounter rates of different kinds of tool sites on both sides of the Uele River, we predicted that these subpopulations would differ in their likelihood of using tools to prey on two insect species despite similar availability. In surveys conducted over a 9-year period at 19 different survey regions in northern Democratic Republic of Congo (10 in lowland forest and 9 in mosaic), we collected and analyzed data on chimpanzee tool-assisted exploitation of insects. To determine the availability of insect species eaten by the chimpanzees, we counted insects and their mounds on transects and recces at 12 of these sites. For stick tools used to harvest epigaeic Dorylus and ponerine ants, we evaluated seasonal, geographical, and prey-availability factors that might influence their occurrence, using nest encounter rate as a proxy to control for chimpanzee abundance. Across the 19 survey regions spanning both sides of the Uele, we found little difference in the availability of epigaeic Dorylus and ponerine ants. Despite this, tool encounter rates for epigaeic Dorylus, but not ponerine, ants were significantly higher in the mosaic to the north of the Uele. Furthermore, we found no evidence for termite fishing anywhere, despite the availability of Macrotermes mounds throughout the region and the fact that chimpanzees at a number of other study sites use tools to harvest these termites. Instead, the chimpanzees of this region used a novel percussive technique to harvest two other types of termites, Cubitermes sp. and Thoracotermes macrothorax. This mismatch between prey availability and predation is consistent with cultural override, but given the different habitats on the two sides of the Uele River, we cannot fully rule out the influence of ecological factors. Comparing our findings with those of similar studies of other chimpanzee populations promises to contribute to our understanding of the evolution of behavioral diversity in humans and our closest cousins.

Taxonomic variation in the supraorbital region of catarrhine primates

OBJECTIVES

This study aimed to test the taxonomic utility of the catarrhine supraorbital region using 3D geometric morphometrics, with the aim of establishing its potential use in elucidating the position of more debated hominin groups.

MATERIALS AND METHODS

230 3D coordinates were used to record the supraorbital morphology of two datasets: one containing 460 non-hominin catarrhine primates from species and subspecies of Gorilla, Pan, Papio, and Macaca; and the other containing 55 Pleistocene hominins from Homo, Australopithecus, and Paranthropus. Principal component analyses in tangent, form, and allometry-free shape space were used to assess differentiation of taxa, with biological distinctiveness of taxa being established using step-wise discriminant analysis with subsampling.

RESULTS

Results indicated that the recorded supraorbital morphology could be used to separate non-hominin catarrhine primate genera, species, and subspecies, although accuracy was found to decrease with decreasing Linnaean rank. In addition, analyses in tangent space were found to produce the highest accuracy when classifying primates of known taxonomy. Biological distinctiveness of the middle and later Homo species was comparable to or higher than that of the non-hominin primates, and relatively lower for the earlier groups of Homo.

DISCUSSION

This study indicates that the supraorbital region preserves taxonomic information that can be used to delineate between closely related groups, both within hominins and wider catarrhine primates. Therefore, this region may be used to provide insight when assessing the taxonomic affiliation of disputed hominin specimens.

Humans and Chimpanzees Display Opposite Patterns of Diversity in Arylamine N-Acetyltransferase Genes

Among the many genes involved in the metabolism of therapeutic drugs, human arylamine N-acetyltransferases (NATs) genes have been extensively studied, due to their medical importance both in pharmacogenetics and disease epidemiology. One member of this small gene family, NAT2, is established as the locus of the classic human acetylation polymorphism in drug metabolism. Current hypotheses hold that selective processes favoring haplotypes conferring lower NAT2 activity have been operating in modern humans’ recent history as an adaptation to local chemical and dietary environments. To shed new light on such hypotheses, we investigated the genetic diversity of the three members of the NAT gene family in seven hominid species, including modern humans, Neanderthals and Denisovans. Little polymorphism sharing was found among hominids, yet all species displayed high NAT diversity, but distributed in an opposite fashion in chimpanzees and bonobos (Pan genus) compared to modern humans, with higher diversity in Pan species at NAT1 and lower at NAT2, while the reverse is observed in humans. This pattern was also reflected in the results returned by selective neutrality tests, which suggest, in agreement with the predicted functional impact of mutations detected in non-human primates, stronger directional selection, presumably purifying selection, at NAT1 in modern humans, and at NAT2 in chimpanzees. Overall, the results point to the evolution of divergent functions of these highly homologous genes in the different primate species, possibly related to their specific chemical/dietary environment (exposome) and we hypothesize that this is likely linked to the emergence of controlled fire use in the human lineage.

Interspecific Gene Flow Shaped the Evolution of the Genus Canis

The evolutionary history of the wolf-like canids of the genus Canis has been heavily debated, especially regarding the number of distinct species and their relationships at the population and species level [1-6]. We assembled a dataset of 48 resequenced genomes spanning all members of the genus Canis except the black-backed and side-striped jackals, encompassing the global diversity of seven extant canid lineages. This includes eight new genomes, including the first resequenced Ethiopian wolf (Canis simensis), one dhole (Cuon alpinus), two East African hunting dogs (Lycaon pictus), two Eurasian golden jackals (Canis aureus), and two Middle Eastern gray wolves (Canis lupus). The relationships between the Ethiopian wolf, African golden wolf, and golden jackal were resolved. We highlight the role of interspecific hybridization in the evolution of this charismatic group. Specifically, we find gene flow between the ancestors of the dhole and African hunting dog and admixture between the gray wolf, coyote (Canis latrans), golden jackal, and African golden wolf. Additionally, we report gene flow from gray and Ethiopian wolves to the African golden wolf, suggesting that the African golden wolf originated through hybridization between these species. Finally, we hypothesize that coyotes and gray wolves carry genetic material derived from a "ghost" basal canid lineage.

Significant loss of mitochondrial diversity within the last century due to extinction of peripheral populations in eastern gorillas

Species and populations are disappearing at an alarming rate as a direct result of human activities. Loss of genetic diversity associated with population decline directly impacts species' long-term survival. Therefore, preserving genetic diversity is of considerable conservation importance. However, to assist in conservation efforts, it is important to understand how genetic diversity is spatially distributed and how it changes due to anthropogenic pressures. In this study, we use historical museum and modern faecal samples of two critically endangered eastern gorilla taxa, Grauer's (Gorilla beringei graueri) and mountain gorillas (Gorilla beringei beringei), to directly infer temporal changes in genetic diversity within the last century. Using over 100 complete mitochondrial genomes, we observe a significant decline in haplotype and nucleotide diversity in Grauer's gorillas. By including historical samples from now extinct populations we show that this decline can be attributed to the loss of peripheral populations rather than a decrease in genetic diversity within the core range of the species. By directly quantifying genetic changes in the recent past, our study shows that human activities have severely impacted eastern gorilla genetic diversity within only four to five generations. This rapid loss calls for dedicated conservation actions, which should include preservation of the remaining peripheral populations.

The impact of endogenous content, replicates and pooling on genome capture from faecal samples

Target-capture approach has improved over the past years, proving to be very efficient tool for selectively sequencing genetic regions of interest. These methods have also allowed the use of noninvasive samples such as faeces (characterized by their low quantity and quality of endogenous DNA) to be used in conservation genomic, evolution and population genetic studies. Here we aim to test different protocols and strategies for exome capture using the Roche SeqCap EZ Developer kit (57.5 Mb). First, we captured a complex pool of DNA libraries. Second, we assessed the influence of using more than one faecal sample, extract and/or library from the same individual, to evaluate its effect on the molecular complexity of the experiment. We validated our experiments with 18 chimpanzee faecal samples collected from two field sites as a part of the Pan African Programme: The Cultured Chimpanzee. Those two field sites are in Kibale National Park, Uganda (N = 9) and Loango National Park, Gabon (N = 9). We demonstrate that at least 16 libraries can be pooled, target enriched through hybridization, and sequenced allowing for the genotyping of 951,949 exome markers for population genetic analyses. Further, we observe that molecule richness, and thus, data acquisition, increase when using multiple libraries from the same extract or multiple extracts from the same sample. Finally, repeated captures significantly decrease the proportion of off-target reads from 34.15% after one capture round to 7.83% after two capture rounds, supporting our conclusion that two rounds of target enrichment are advisable when using complex faecal samples.

Cellular and population level processes influence the rate, accumulation and observed frequency of inherited and somatic mtDNA mutations

Mitochondria are found in all animals and have the unique feature of containing multiple copies of their own small, circular DNA genome (mtDNA). The rate and pattern of mutation accumulation in the mtDNA are influenced by molecular, cellular and population level processes. We distinguish between inherited and somatic mtDNA mutations and review evidence for the often-made assumption that mutations accumulate at a higher rate in mtDNA than in nuclear DNA (nDNA). We conclude that the whole genome mutation accumulation rate is higher for mtDNA than for nDNA but include the caveat that rates overlap considerably between the individual mtDNA- and nDNA-encoded genes. Next, we discuss the postulated causal mechanisms for the high rate of mtDNA mutation accumulation in both inheritance and in somatic cells. Perhaps unexpectedly, mtDNA is resilient to many mutagens of nDNA but is prone to errors of replication. We then consider the influence of maternal inheritance, recombination and selection on the observed accumulation pattern of inherited mtDNA mutations. Finally, we discuss environmental influences of temperature and diet on the observed frequency of inherited and somatic mtDNA mutations. We conclude that it is necessary to understand the cellular processes to fully interpret the pattern of mutations and how they influence our interpretations of evolution and disease.

Limited MHC class I intron 2 repertoire variation in bonobos

Common chimpanzees (Pan troglodytes) experienced a selective sweep, probably caused by a SIV-like virus, which targeted their MHC class I repertoire. Based on MHC class I intron 2 data analyses, this selective sweep took place about 2-3 million years ago. As a consequence, common chimpanzees have a skewed MHC class I repertoire that is enriched for allotypes that are able to recognise conserved regions of the SIV proteome. The bonobo (Pan paniscus) shared an ancestor with common chimpanzees approximately 1.5 to 2 million years ago. To investigate whether the signature of this selective sweep is also detectable in bonobos, the MHC class I gene repertoire of two bonobo panels comprising in total 29 animals was investigated by Sanger sequencing. We identified 14 Papa-A, 20 Papa-B and 11 Papa-C alleles, of which eight, five and eight alleles, respectively, have not been reported previously. Within this pool of MHC class I variation, we recovered only 2 Papa-A, 3 Papa-B and 6 Papa-C intron 2 sequences. As compared to humans, bonobos appear to have an even more diminished MHC class I intron 2 lineage repertoire than common chimpanzees. This supports the notion that the selective sweep may have predated the speciation of common chimpanzees and bonobos. The further reduction of the MHC class I intron 2 lineage repertoire observed in bonobos as compared to the common chimpanzee may be explained by a founding effect or other subsequent selective processes.

The Y chromosomes of the great apes

The great apes (orangutans, gorillas, chimpanzees, bonobos and humans) descended from a common ancestor around 13 million years ago, and since then their sex chromosomes have followed very different evolutionary paths. While great-ape X chromosomes are highly conserved, their Y chromosomes, reflecting the general lability and degeneration of this male-specific part of the genome since its early mammalian origin, have evolved rapidly both between and within species. Understanding great-ape Y chromosome structure, gene content and diversity would provide a valuable evolutionary context for the human Y, and would also illuminate sex-biased behaviours, and the effects of the evolutionary pressures exerted by different mating strategies on this male-specific part of the genome. High-quality Y-chromosome sequences are available for human and chimpanzee (and low-quality for gorilla). The chromosomes differ in size, sequence organisation and content, and while retaining a relatively stable set of ancestral single-copy genes, show considerable variation in content and copy number of ampliconic multi-copy genes. Studies of Y-chromosome diversity in other great apes are relatively undeveloped compared to those in humans, but have nevertheless provided insights into speciation, dispersal, and mating patterns. Future studies, including data from larger sample sizes of wild-born and geographically well-defined individuals, and full Y-chromosome sequences from bonobos, gorillas and orangutans, promise to further our understanding of population histories, male-biased behaviours, mutation processes, and the functions of Y-chromosomal genes.

Chimpanzee genomic diversity reveals ancient admixture with bonobos

Our closest living relatives, chimpanzees and bonobos, have a complex demographic history. We analyzed the high-coverage whole genomes of 75 wild-born chimpanzees and bonobos from 10 countries in Africa. We found that chimpanzee population substructure makes genetic information a good predictor of geographic origin at country and regional scales. Multiple lines of evidence suggest that gene flow occurred from bonobos into the ancestors of central and eastern chimpanzees between 200,000 and 550,000 years ago, probably with subsequent spread into Nigeria-Cameroon chimpanzees. Together with another, possibly more recent contact (after 200,000 years ago), bonobos contributed less than 1% to the central chimpanzee genomes. Admixture thus appears to have been widespread during hominid evolution.