Evidence from Cameroon reveals differences in the genetic structure and histories of chimpanzee populations

Genetic tracing of the illegal trade of the white-bellied pangolin (Phataginus tricuspis) in western Central Africa

The white-bellied pangolin is subject to intense trafficking, feeding both local and international trade networks. In order to assess its population genetics and trace its domestic trade, we genotyped 562 pangolins from local to large bushmeat markets in western central Africa. We show that the two lineages described from the study region (WCA and Gab) were overlapping in ranges, with limited introgression in southern Cameroon. There was a lack of genetic differentiation across WCA and a significant signature of isolation-by-distance possibly due to unsuspected dispersal capacities involving a Wahlund effect. We detected a c. 74.1-82.5% decline in the effective population size of WCA during the Middle Holocene. Private allele frequency tracing approach indicated up to 600 km sourcing distance by large urban markets from Cameroon, including Equatorial Guinea. The 20 species-specific microsatellite loci provided individual-level genotyping resolution and should be considered as valuable resources for future forensic applications. Because admixture was detected between lineages, we recommend a multi-locus approach for tracing the pangolin trade. The Yaoundé market was the main hub of the trade in the region, and thus should receive specific monitoring to mitigate pangolins' domestic trafficking. Our study also highlighted the weak implementation of CITES regulations at European borders.

Morphological variation of the Pan talus relative to that of Gorilla

OBJECTIVES

Differences in talar articular morphology relative to locomotion have recently been found within Pan and Gorilla. Whole-bone talar morphology within, and shared variation among, Pan and Gorilla (sub)species, however, has yet to be investigated. Here we separately analyze talar external shape within Pan (P. t. troglodytes, P. t. schweinfurthii, P. t. verus, P. paniscus) and Gorilla (G. g. gorilla, G. b. beringei, G. b. graueri) relative to degree of arboreality and body size. Pan and Gorilla are additionally analyzed together to determine if consistent shape differences exist within the genera.

MATERIALS AND METHODS

Talar external shape was quantified using a weighted spherical harmonic analysis. Shape variation both within and among Pan and Gorilla was described using principal component analyses. Root mean square distances were calculated between taxon averages, and resampling statistics conducted to test for pairwise differences.

RESULTS

P. t. verus (most arboreal Pan) talar shape significantly differs from other Pan taxa (p < 0.05 for pairwise comparisons) driven by more asymmetrical trochlear rims and a medially-set talar head. P. t. troglodytes, P. t. schweinfurthii, and P. paniscus do not significantly differ (p > 0.05 for pairwise comparisons). All gorilla taxa exhibit significantly different talar morphologies (p < 0.007 for pairwise comparisons). The more terrestrial subspecies of G. beringei and P. troglodytes exhibit a superoinferiorly taller talar head/neck complex.

DISCUSSION

P. t. verus exhibits talar morphologies that have been previously related to more frequent arboreality. The adaptations in the more terrestrial G. beringei and P. troglodytes subspecies may serve to facilitate load transmission.

Evidence for ecotone speciation across an African rainforest-savanna gradient

Accelerating climate change and habitat loss make it imperative that plans to conserve biodiversity consider species' ability to adapt to changing environments. However, in biomes where biodiversity is highest, the evolutionary mechanisms responsible for generating adaptative variation and, ultimately, new species are frequently poorly understood. African rainforests represent one such biome, as decadal debates continue concerning the mechanisms generating African rainforest biodiversity. These debates hinge on the relative importance of geographic isolation versus divergent natural selection across environmental gradients. Hindering progress is a lack of robust tests of these competing hypotheses. Because African rainforests are severely at-risk due to climate change and other anthropogenic activities, addressing this long-standing debate is critical for making informed conservation decisions. We use demographic inference and allele frequency-environment relationships to investigate mechanisms of diversification in an African rainforest skink, Trachylepis affinis, a species inhabiting the gradient between rainforest and rainforest-savanna mosaic (ecotone). We provide compelling evidence of ecotone speciation, in which gene flow has all but ceased between rainforest and ecotone populations, at a level consistent with infrequent hybridization between sister species. Parallel patterns of genomic, morphological, and physiological divergence across this environmental gradient and pronounced allele frequency-environment correlation indicate speciation is mostly probably driven by ecological divergence, supporting a central role for divergent natural selection. Our results provide strong evidence for the importance of ecological gradients in African rainforest speciation and inform conservation strategies that preserve the processes that produce and maintain biodiversity.

Morphological variation of the maxilla in modern humans and African apes

Differences in morphology among modern humans and African apes are frequently used when assessing whether hominin fossils should be attributed to a single species or represent evidence for taxic diversity. A good understanding of the degree and structure of the intergeneric, interspecific, and intraspecific variation, including aspects such as sexual dimorphism and age, are key in this context. Here we explore the variation and differences shown by the maxilla of extant hominines, as maxillary morphology is central in the diagnosis of several hominin taxa. Our sample includes adults of all currently recognized hominine species and subspecies, with a balanced species sex ratio. In addition, we compared the adults with a small sample of late juveniles. The morphology of the maxillae was captured using three-dimensional landmarks, and the size and shape were analyzed using geometric morphometric methods. Key observations are that 1) the maxillae of all extant hominine species and subspecies show statistically significant differences, but complete separation in shape is only seen at the genus level; 2) the degree of variation is not consistent between genera, with subspecies of Gorilla being more different from each other than are species of Pan; 3) the pattern of sexual shape dimorphism is different in Pan, Gorilla, and Homo, often showing opposite trends; and 4) differentiation between maxillary shapes is increased after adjustment for static intraspecific allometry. These results provide a taxonomically up-to-date comparative morphological framework to help interpret the hominin fossil record, and we discuss the practical implications in that context.

Evolutionary Genetic Signatures of Selection on Bone-Related Variation within Human and Chimpanzee Populations

Bone strength and the incidence and severity of skeletal disorders vary significantly among human populations, due in part to underlying genetic differentiation. While clinical models predict that this variation is largely deleterious, natural population variation unrelated to disease can go unnoticed, altering our perception of how natural selection has shaped bone morphologies over deep and recent time periods. Here, we conduct the first comparative population-based genetic analysis of the main bone structural protein gene, collagen type I α 1 (COL1A1), in clinical and 1000 Genomes Project datasets in humans, and in natural populations of chimpanzees. Contrary to predictions from clinical studies, we reveal abundant COL1A1 amino acid variation, predicted to have little association with disease in the natural population. We also find signatures of positive selection associated with intron haplotype structure, linkage disequilibrium, and population differentiation in regions of known gene expression regulation in humans and chimpanzees. These results recall how recent and deep evolutionary regimes can be linked, in that bone morphology differences that developed among vertebrates over 450 million years of evolution are the result of positive selection on subtle type I collagen functional variation segregating within populations over time.

Quantitative estimates of glacial refugia for chimpanzees (Pan troglodytes) since the Last Interglacial (120,000 BP)

Paleoclimate reconstructions have enhanced our understanding of how past climates have shaped present-day biodiversity. We hypothesize that the geographic extent of Pleistocene forest refugia and suitable habitat fluctuated significantly in time during the late Quaternary for chimpanzees (Pan troglodytes). Using bioclimatic variables representing monthly temperature and precipitation estimates, past human population density data, and an extensive database of georeferenced presence points, we built a model of changing habitat suitability for chimpanzees at fine spatio-temporal scales dating back to the Last Interglacial (120,000 BP). Our models cover a spatial resolution of 0.0467° (approximately 5.19 km² grid cells) and a temporal resolution of between 1000 and 4000 years. Using our model, we mapped habitat stability over time using three approaches, comparing our modeled stability estimates to existing knowledge of Afrotropical refugia, as well as contemporary patterns of major keystone tropical food resources used by chimpanzees, figs (Moraceae), and palms (Arecacae). Results show habitat stability congruent with known glacial refugia across Africa, suggesting their extents may have been underestimated for chimpanzees, with potentially up to approximately 60,000 km² of previously unrecognized glacial refugia. The refugia we highlight coincide with higher species richness for figs and palms. Our results provide spatio-temporally explicit insights into the role of refugia across the chimpanzee range, forming the empirical foundation for developing and testing hypotheses about behavioral, ecological, and genetic diversity with additional data. This methodology can be applied to other species and geographic areas when sufficient data are available.

Structure of Chimpanzee Gut Microbiomes across Tropical Africa

Understanding variation in host-associated microbial communities is important given the relevance of microbiomes to host physiology and health. Using 560 fecal samples collected from wild chimpanzees (Pan troglodytes) across their range, we assessed how geography, genetics, climate, vegetation, and diet relate to gut microbial community structure (prokaryotes, eukaryotic parasites) at multiple spatial scales. We observed a high degree of regional specificity in the microbiome composition, which was associated with host genetics, available plant foods, and potentially with cultural differences in tool use, which affect diet. Genetic differences drove community composition at large scales, while vegetation and potentially tool use drove within-region differences, likely due to their influence on diet. Unlike industrialized human populations in the United States, where regional differences in the gut microbiome are undetectable, chimpanzee gut microbiomes are far more variable across space, suggesting that technological developments have decoupled humans from their local environments, obscuring regional differences that could have been important during human evolution.

IMPORTANCE Gut microbial communities are drivers of primate physiology and health, but the factors that influence the gut microbiome in wild primate populations remain largely undetermined. We report data from a continent-wide survey of wild chimpanzee gut microbiota and highlight the effects of genetics, vegetation, and potentially even tool use at different spatial scales on the chimpanzee gut microbiome, including bacteria, archaea, and eukaryotic parasites. Microbial community dissimilarity was strongly correlated with chimpanzee population genetic dissimilarity, and vegetation composition and consumption of algae, honey, nuts, and termites were potentially associated with additional divergence in microbial communities between sampling sites. Our results suggest that host genetics, geography, and climate play a far stronger role in structuring the gut microbiome in chimpanzees than in humans.

WP, Sterkhova VV, Kusamba C, Rödel M, Penner J, Peterson AT, Brown RM. Rivers, not refugia, drove diversification in arboreal, sub-Saharan African snakes

The relative roles of rivers versus refugia in shaping the high levels of species diversity in tropical rainforests have been widely debated for decades. Only recently has it become possible to take an integrative approach to test predictions derived from these hypotheses using genomic sequencing and paleo-species distribution modeling. Herein, we tested the predictions of the classic river, refuge, and river-refuge hypotheses on diversification in the arboreal sub-Saharan African snake genus Toxicodryas. We used dated phylogeographic inferences, population clustering analyses, demographic model selection, and paleo-distribution modeling to conduct a phylogenomic and historical demographic analysis of this genus. Our results revealed significant population genetic structure within both Toxicodryas species, corresponding geographically to river barriers and divergence times from the mid-Miocene to Pliocene. Our demographic analyses supported the interpretation that rivers are indications of strong barriers to gene flow among populations since their divergence. Additionally, we found no support for a major contraction of suitable habitat during the last glacial maximum, allowing us to reject both the refuge and river-refuge hypotheses in favor of the river-barrier hypothesis. Based on conservative interpretations of our species delimitation analyses with the Sanger and ddRAD data sets, two new cryptic species are identified from east-central Africa. This study highlights the complexity of diversification dynamics in the African tropics and the advantages of integrative approaches to studying speciation in tropical regions.

Hybrid apes in the Anthropocene: Burden or asset for conservation?

Conservationists often view hybrid animals as problematic, at least if anthropogenic influence caused the intermixing to occur. However, critics propose that humans should respect non-human autonomy, reject and accept the creatures they have helped to create.Based on two case studies of our own ethological, genetic and ethnographic research about chimpanzee and orangutan subspecies hybrids, we assess what, if anything, should be done about such animals. We consider problems posed by cross-bred apes relating to: (a) Breeding-Do hybrids really experience reduced reproductive success? How are population-level concerns and welfare of individual animals balanced in conservation breeding? (b) Essentialism-Are anti-hybrid arguments based on essentialist or purist thinking? Does essentialism vary by conservation context? (c) Pragmatism-How do socio-economic circumstances influence whether hybrids are embraced or ignored? Does the erosion of 'untouched nature' render hybrids more important?We show that answers to these questions are complex and context-specific, and that therefore decisions should be made on a case-by-case basis. For example, we find that anti-hybrid arguments are essentialist in some cases (e.g. ape management in zoos) but not in others (e.g. ape reintroduction). Thus, rather than present recommendations, we conclude by posing nine questions that conservationists should ask themselves when making decisions about taxonomic hybrids. ​ A free Plain Language Summary can be found within the Supporting Information of this article.

Tectonics, climate and the diversification of the tropical African terrestrial flora and fauna

Tropical Africa is home to an astonishing biodiversity occurring in a variety of ecosystems. Past climatic change and geological events have impacted the evolution and diversification of this biodiversity. During the last two decades, around 90 dated molecular phylogenies of different clades across animals and plants have been published leading to an increased understanding of the diversification and speciation processes generating tropical African biodiversity. In parallel, extended geological and palaeoclimatic records together with detailed numerical simulations have refined our understanding of past geological and climatic changes in Africa. To date, these important advances have not been reviewed within a common framework. Here, we critically review and synthesize African climate, tectonics and terrestrial biodiversity evolution throughout the Cenozoic to the mid-Pleistocene, drawing on recent advances in Earth and life sciences. We first review six major geo-climatic periods defining tropical African biodiversity diversification by synthesizing 89 dated molecular phylogeny studies. Two major geo-climatic factors impacting the diversification of the sub-Saharan biota are highlighted. First, Africa underwent numerous climatic fluctuations at ancient and more recent timescales, with tectonic, greenhouse gas, and orbital forcing stimulating diversification. Second, increased aridification since the Late Eocene led to important extinction events, but also provided unique diversification opportunities shaping the current tropical African biodiversity landscape. We then review diversification studies of tropical terrestrial animal and plant clades and discuss three major models of speciation: (i) geographic speciation via vicariance (allopatry); (ii) ecological speciation impacted by climate and geological changes, and (iii) genomic speciation via genome duplication. Geographic speciation has been the most widely documented to date and is a common speciation model across tropical Africa. We conclude with four important challenges faced by tropical African biodiversity research: (i) to increase knowledge by gathering basic and fundamental biodiversity information; (ii) to improve modelling of African geophysical evolution throughout the Cenozoic via better constraints and downscaling approaches; (iii) to increase the precision of phylogenetic reconstruction and molecular dating of tropical African clades by using next generation sequencing approaches together with better fossil calibrations; (iv) finally, as done here, to integrate data better from Earth and life sciences by focusing on the interdisciplinary study of the evolution of tropical African biodiversity in a wider geodiversity context.

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.

Environmental variability supports chimpanzee behavioural diversity

Large brains and behavioural innovation are positively correlated, species-specific traits, associated with the behavioural flexibility animals need for adapting to seasonal and unpredictable habitats. Similar ecological challenges would have been important drivers throughout human evolution. However, studies examining the influence of environmental variability on within-species behavioural diversity are lacking despite the critical assumption that population diversification precedes genetic divergence and speciation. Here, using a dataset of 144 wild chimpanzee (Pan troglodytes) communities, we show that chimpanzees exhibit greater behavioural diversity in environments with more variability — in both recent and historical timescales. Notably, distance from Pleistocene forest refugia is associated with the presence of a larger number of behavioural traits, including both tool and non-tool use behaviours. Since more than half of the behaviours investigated are also likely to be cultural, we suggest that environmental variability was a critical evolutionary force promoting the behavioural, as well as cultural diversification of great apes.

Environmental variability is one potential driver of behavioural and cultural diversity in humans and other animals. Here, the authors show that chimpanzee behavioural diversity is higher in habitats that are more seasonal and historically unstable, and in savannah woodland relative to forested sites.

A review of research in primate sanctuaries

While non-human primate studies have long been conducted in laboratories, and more recently at zoological parks, sanctuaries are increasingly considered a viable setting for research. Accredited sanctuaries in non-range countries house thousands of primates formerly used as subjects of medical research, trained performers or personal pets. In range countries, however, sanctuaries typically house orphaned primates confiscated from illegal poaching and the bushmeat and pet trafficking trades. Although the primary mission of these sanctuaries is to rescue and rehabilitate residents, many of these organizations are increasingly willing to participate in non-invasive research. Notably, from a scientific standpoint, most sanctuaries provide potential advantages over traditional settings, such as large, naturalistic physical and social environments which may result in more relevant models of primates' free-ranging wild counterparts than other captive settings. As a result, an impressive scope of research in the fields of primate behaviour, cognition, veterinary science, genetics and physiology have been studied in sanctuaries. In this review, we examine the range and form of research that has been conducted at accredited sanctuaries around the world. We also describe the potential challenges of sanctuary-based work and the considerations that external researchers may face when deciding to collaborate with primate sanctuaries on their research projects.

Phylogenetic analyses reveal that simian foamy virus isolated from Japanese Yakushima macaques (Macaca fuscata yakui) is distinct from most of Japanese Hondo macaques (Macaca fuscata fuscata)

Japanese macaque (Macaca fuscata) is an indigenous Old World monkey (OWM) species that inhabits the Japanese archipelago. There are two subspecies of Japanese macaque: Yakushima macaque (M. f. yakui) which inhabits Yakushima Island exclusively, and Hondo macaque (M. f. fuscata) which inhabits the mainland of Japan. Yakushima macaque is considered to be branched off from a certain parental macaque group that had inhabited the mainland of Japan. However, the process of sub-speciation of the Yakushima macaque is still unclear at present. In this study, to gain new insight into the process of sub-speciation of Japanese macaque, we utilized the simian foamy virus (SFV) as a marker. SFVs are found in virtually all primates except humans and undergo species-specific cospeciation with the hosts. The phylogenetic analysis of conserved regions of the env gene in SFVs remarkably resembled that of the OWMs with high statistical confidence. The phylogenetic analyses also indicated that there are four (1-4) genotypes among Asian OWMs investigated. SFVs derived from Asian OWMs except Yakushima macaque were classified as genotypes 1-3, whereas SFVs isolated from all Yakushima macaques and one Hondo macaque were classified as genotype 4. Interestingly, genotype 4 was firstly branched off from the rest of the genotypes, which might indicate that the macaques infected with genotype 4 SFV were derived from the "older" population of Japanese macaques. The high prevalence of genotype 4 SFVs among Yakushima macaque might reflect the possibility that they are a descendant of the population settled earlier, which has been geographically isolated in Yakushima Island.

Subspecific identity and a comparison of genetic diversity between wild and ex situ wildebeest

The original North American ex situ wildebeest population was believed to originate from the white-bearded wildebeest (Connochaetes taurinus albojubatus), which is both morphologically distinct and geographically separated from the brindled wildebeest (C. t. taurinus). However, after an import of wildebeest into North America in 2001, managers have suspected that white-bearded and brindled wildebeest were mixed in herds at multiple institutions. We sequenced the mitochondrial control region (d-loop) from a portion of the managed North American population and compared our sequences with previously published sequences from wild individuals to determine the subspecific identity and genetic diversity of our ex situ population. We were able to confidently identify C. t. albojubatus as the subspecies identity of the sampled portion of our population. Within our population, haplotype and nucleotide diversity were low (0.169 and 0.001, respectively) with a single common haplotype (H1) containing 41 of the 45 individuals sequenced, while two rare haplotypes (H2 and H3) were derived from three individuals and a single individual, respectively. Nucleotide and haplotype diversity were greater overall in the wild populations compared with our managed population. However, C. t. albojubatus was found to exhibit lower nucleotide diversity in both wild and ex situ populations when compared to other wild subspecies. Though the overall goal of the North American wildebeest population is for public education and not reintroduction, maintaining genetic diversity is vital for the long-term viability of this managed population, which may benefit from periodic supplementation of wild animals.

Phylogeographic evidence for two species of muriqui (genus Brachyteles)

The taxonomy of muriquis, the largest extant primates in the New World, is controversial. While some specialists argue for a monotypic genus (Brachyteles arachnoides), others favor a two-species classification, splitting northern muriquis (Brachyteles hypoxanthus) from southern muriquis (B. arachnoides). This uncertainty affects how we study the differences between these highly endangered and charismatic primates, as well as the design of more effective conservation programs. To address this issue, between 2003 and 2017 we collected over 230 muriqui fecal samples across the genus' distribution in the Brazilian Atlantic Forest, extracted DNA from these samples, and sequenced 423 base pairs of the mitochondrial DNA (mtDNA) control region. Phylogenetic and species delimitation analyses of our sequence dataset robustly support two reciprocally monophyletic groups corresponding to northern and southern muriquis separated by an average 12.7% genetic distance. The phylogeographic break between these lineages seems to be associated with the Paraíba do Sul River and coincides with the transition between the north and south Atlantic Forest biogeographic zones. Published divergence estimates from whole mitochondrial genomes and nuclear loci date the split between northern and southern muriquis to the Early Pleistocene (ca. 2.0 mya), and our new mtDNA dataset places the coalescence time for each of these two clades near the last interglacial (ca. 120-80 kya). Our results, together with both phenotypic and ecological differences, support recognizing northern and southern muriquis as sister species that should be managed as distinct evolutionarily significant units. Given that only a few thousand muriquis remain in nature, it is imperative that conservation strategies are tailored to protect both species from extinction.

Reconstructing the ancestral phenotypes of great apes and humans (Homininae) using subspecies-level phylogenies

Owing to their close affinity, the African great apes are of interest in the study of human evolution. Although numerous researchers have described the ancestors we share with these species with reference to extant great apes, few have done so with phylogenetic comparative methods. One obstacle to the application of these techniques is the within-species phenotypic variation found in this group. Here, we leverage this variation, modelling common ancestors using ancestral state reconstructions (ASRs) with reference to subspecies-level trait data. A subspecies-level phylogeny of the African great apes and humans was estimated from full-genome mitochondrial DNA sequences and used to implement ASRs for 14 continuous traits known to vary between great ape subspecies. Although the inclusion of within-species phenotypic variation increased the phylogenetic signal for our traits and improved the performance of our ASRs, whether this was done through the inclusion of subspecies phylogeny or through the use of existing methods made little difference. Our ASRs corroborate previous findings that the last common ancestor of humans, chimpanzees and bonobos was a chimp-like animal, but also suggest that the last common ancestor of humans, chimpanzees, bonobos and gorillas was an animal unlike any extant African great ape.

A mathematical landmark-based method for measuring worn molars in hominoid systematics

Worn teeth pose a major limitation to researchers in the fields of extinct and extant hominoid systematics because they lack clearly identifiable anatomical landmarks needed to take measurements on the crown enamel surface and are typically discarded from a study. This is particularly detrimental when sample sizes for some groups are already characteristically low, if there is an imbalance between samples representing populations, sexes or dietary strategies, or if the worn teeth in question are type specimens of fossil species or other key specimens. This study proposes a methodology based predominantly on mathematically-derived landmarks for measuring size and shape features of molars, irrespective of wear. With 110 specimens of lower second molars from five species of extant hominoids (Pan troglodytes, P. paniscus, Gorilla gorilla, G. beringei, Homo sapiens), n ≥ 20 per species, n ≥ 10 per subspecies, good species separation in morphospace is achieved in a principal components analysis. Classification accuracy in a discriminant function analysis is 96.4% at the species level and 88.2% at the subspecies level (92.7% and 79.1%, respectively, on cross-validation). The classification accuracy compares favorably to that achieved by anatomically-derived measurements based on published research (94% and 84% at the species and subspecies level respectively; 91% and 76% on cross-validation). The mathematical landmarking methodology is rapid and uncomplicated. The results support the use of mathematical landmarks to enable the inclusion of worn molar teeth in dental studies so as to maximize sample sizes and restore balance between populations and/or sexes in hominoid systematic studies.

CD4 receptor diversity in chimpanzees protects against SIV infection

Human and simian immunodeficiency viruses (HIV/SIVs) use CD4 as the primary receptor to enter target cells. Here, we show that the chimpanzee CD4 is highly polymorphic, with nine coding variants present in wild populations, and that this diversity interferes with SIV envelope (Env)-CD4 interactions. Testing the replication fitness of SIVcpz strains in CD4+ T cells from captive chimpanzees, we found that certain viruses were unable to infect cells from certain hosts. These differences were recapitulated in CD4 transfection assays, which revealed a strong association between CD4 genotypes and SIVcpz infection phenotypes. The most striking differences were observed for three substitutions (Q25R, Q40R, and P68T), with P68T generating a second N-linked glycosylation site (N66) in addition to an invariant N32 encoded by all chimpanzee CD4 alleles. In silico modeling and site-directed mutagenesis identified charged residues at the CD4-Env interface and clashes between CD4- and Env-encoded glycans as mechanisms of inhibition. CD4 polymorphisms also reduced Env-mediated cell entry of monkey SIVs, which was dependent on at least one D1 domain glycan. CD4 allele frequencies varied among wild chimpanzees, with high diversity in all but the western subspecies, which appeared to have undergone a selective sweep. One allele was associated with lower SIVcpz prevalence rates in the wild. These results indicate that substitutions in the D1 domain of the chimpanzee CD4 can prevent SIV cell entry. Although some SIVcpz strains have adapted to utilize these variants, CD4 diversity is maintained, protecting chimpanzees against infection with SIVcpz and other SIVs to which they are exposed.

Habitat differentiation among three Nigeria-Cameroon chimpanzee (Pan troglodytes ellioti) populations

Ecological niche models (ENMs) are often used to predict species distribution patterns from datasets that describe abiotic and biotic factors at coarse spatial scales. Ground-truthing ENMs provide important information about how these factors relate to species-specific requirements at a scale that is biologically relevant for the species. Chimpanzees are territorial and have a predominantly frugivorous diet. The spatial and temporal variation in fruit availability for different chimpanzee populations is thus crucial, but rarely depicted in ENMs. The genetic and geographic distinction within Nigeria-Cameroon chimpanzee (Pan troglodytes ellioti) populations represents a unique opportunity to understand fine scale species-relevant ecological variation in relation to ENMs. In Cameroon, P. t. ellioti is composed of two genetically distinct populations that occupy different niches: rainforests in western Cameroon and forest-woodland-savanna mosaic (ecotone) in central Cameroon. We investigated habitat variation at three representative sites using chimpanzee-relevant environmental variables, including fruit availability, to assess how these variables distinguish these niches from one another. Contrary to the assumption of most ENM studies that intact forest is essential for the survival of chimpanzees, we hypothesized that the ecotone and human-modified habitats in Cameroon have sufficient resources to sustain large chimpanzee populations. Rainfall, and the diversity, density, and size of trees were higher at the rainforest. The ecotone had a higher density of terrestrial herbs and lianas. Fruit availability was higher at Ganga (ecotone) than at Bekob and Njuma. Seasonal variation in fruit availability was highest at Ganga, and periods of fruit scarcity were longer than at the rainforest sites. Introduced and secondary forest species linked with anthropogenic modification were common at Bekob, which reduced seasonality in fruit availability. Our findings highlight the value of incorporating fine scale species-relevant ecological data to create more realistic models, which have implications for local conservation planning efforts.

Rediscovery of Red Wolf Ghost Alleles in a Canid Population Along the American Gulf Coast

Rediscovering species once thought to be extinct or on the edge of extinction is rare. Red wolves have been extinct along the American Gulf Coast since 1980, with their last populations found in coastal Louisiana and Texas. We report the rediscovery of red wolf ghost alleles in a canid population on Galveston Island, Texas. We analyzed over 7000 single nucleotide polymorphisms (SNPs) in 60 canid representatives from all legally recognized North American Canis species and two phenotypically ambiguous canids from Galveston Island. We found notably high Bayesian cluster assignments of the Galveston canids to captive red wolves with extensive sharing of red wolf private alleles. Today, the only known extant wild red wolves persist in a reintroduced population in North Carolina, which is dwindling amongst political and taxonomic controversy. Our rediscovery of red wolf ancestry after almost 40 years introduces both positive opportunities for additional conservation action and difficult policy challenges.

Phylogeny Estimation by Integration over Isolation with Migration Models

Phylogeny estimation is difficult for closely related populations and species, especially if they have been exchanging genes. We present a hierarchical Bayesian, Markov-chain Monte Carlo method with a state space that includes all possible phylogenies in a full Isolation-with-Migration model framework. The method is based on a new type of genealogy augmentation called a "hidden genealogy" that enables efficient updating of the phylogeny. This is the first likelihood-based method to fully incorporate directional gene flow and genetic drift for estimation of a species or population phylogeny. Application to human hunter-gatherer populations from Africa revealed a clear phylogenetic history, with strong support for gene exchange with an unsampled ghost population, and relatively ancient divergence between a ghost population and modern human populations, consistent with human/archaic divergence. In contrast, a study of five chimpanzee populations reveals a clear phylogeny with several pairs of populations having exchanged DNA, but does not support a history with an unsampled ghost population.

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.

When Viruses Don't Go Viral: The Importance of Host Phylogeographic Structure in the Spatial Spread of Arenaviruses

Many emerging infections are RNA virus spillovers from animal reservoirs. Reservoir identification is necessary for predicting the geographic extent of infection risk, but rarely are taxonomic levels below the animal species considered as reservoir, and only key circumstances in nature and methodology allow intrinsic virus-host associations to be distinguished from simple geographic (co-)isolation. We sampled and genetically characterized in detail a contact zone of two subtaxa of the rodent Mastomys natalensis in Tanzania. We find two distinct arenaviruses, Gairo and Morogoro virus, each spatially confined to a single M. natalensis subtaxon, only co-occurring at the contact zone’s centre. Inter-subtaxon hybridization at this centre and a continuum of quality habitat for M. natalensis show that both viruses have the ecological opportunity to spread into the other substaxon’s range, but do not, strongly suggesting host-intrinsic barriers. Such barriers could explain why human cases of another M. natalensis-borne arenavirus, Lassa virus, are limited to West Africa.

Reservoirs of zoonotic viruses are usually equated with a particular wildlife species. It is rarely assessed whether genetic groups below the species level may instead represent the actual reservoir, though this would have major implications on estimations of the zoonosis’ spatial distribution. Here we investigate whether geographically and genetically distinct subtaxa of the widespread African rodent Mastomys natalensis carry distinct arenaviruses, by sampling in detail across a contact zone of two of these subtaxa. Ongoing hybridization shows that individuals of the subtaxa are in direct physical contact, in principle allowing viral exchange, yet neither of the two arenaviruses -Gairo and Morogoro virus- were found to have crossed the zone. Such intraspecific genetic barriers to arenavirus spatial spread have important implications for our understanding of the related Lassa arenavirus, a pathogen potentially lethal to humans of which Mastomys natalensis is also the main reservoir. Although Lassa virus appears to infect several secondary hosts, its distribution is restricted to West Africa and matches that of another M. natalensis subtaxon. Our data thus indicates that it is because of M. natalensis intraspecific distinctions that the human Lassa fever endemic area has not expanded to the rest of sub-Saharan Africa.

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.

Evolution and demography of the great apes

The great apes are the closest living relatives of humans. Chimpanzees and bonobos group together with humans, while gorillas and orangutans are more divergent from humans. Here, we review insights into their evolution pertaining to the topology of species and subspecies and the reconstruction of their demography based on genome-wide variation. These advances have only become possible recently through next-generation sequencing technologies. Given the close relationship to humans, they provide an important evolutionary context for human genetics.

Why Are Nigeria-Cameroon Chimpanzees (Pan troglodytes ellioti) Free of SIVcpz Infection?

Simian immunodeficiency virus (SIV) naturally infects two subspecies of chimpanzee: Pan troglodytes troglodytes from Central Africa (SIVcpzPtt) and P. t. schweinfurtii from East Africa (SIVcpzPts), but is absent in P. t. verus from West Africa and appears to be absent in P. t. ellioti inhabiting Nigeria and western Cameroon. One explanation for this pattern is that P. t. troglodytes and P. t schweinfurthii may have acquired SIVcpz after their divergence from P. t. verus and P. t. ellioti. However, all of the subspecies, except P. t. verus, still occasionally exchange migrants making the absence of SIVcpz in P. t. ellioti puzzling. Sampling of P. t. ellioti has been minimal to date, particularly along the banks of the Sanaga River, where its range abuts that of P. t. troglodytes. This study had three objectives. First, we extended the sampling of SIVcpz across the range of chimpanzees north of the Sanaga River to address whether under-sampling might account for the absence of evidence for SIVcpz infection in P. t. ellioti. Second, we investigated how environmental variation is associated with the spread and prevalence of SIVcpz in the two chimpanzee subspecies inhabiting Cameroon since environmental variation has been shown to contribute to their divergence from one another. Finally, we compared the prevalence and distribution of SIVcpz with that of Simian Foamy Virus (SFV) to examine the role of ecology and behavior in shaping the distribution of diseases in wild host populations. The dataset includes previously published results on SIVcpz infection and SFVcpz as well as newly collected data, and represents over 1000 chimpanzee fecal samples from 41 locations across Cameroon. Results revealed that none of the 181 P. t. ellioti fecal samples collected across the range of P. t. ellioti tested positive for SIVcpz. In addition, species distribution models suggest that environmental variation contributes to differences in the distribution and prevalence of SIVcpz and SFVcpz. The ecological niches of these two viruses are largely non-overlapping, although stronger statistical support for this conclusion will require more sampling. Overall this study demonstrates that SIVcpz infection is absent or very rare in P. t. ellioti, despite multiple opportunities for transmission. The reasons for its absence remain unclear, but might be explained by one or more factors, including environmental variation, viral competition, and/or local adaptation—all of which should be explored in greater detail through continued surveillance of this region.

Y-Chromosome Structural Diversity in the Bonobo and Chimpanzee Lineages

The male-specific regions of primate Y-chromosomes (MSY) are enriched for multi-copy genes highly expressed in the testis. These genes are located in large repetitive sequences arranged as palindromes, inverted-, and tandem repeats termed amplicons. In humans, these genes have critical roles in male fertility and are essential for the production of sperm. The structure of human and chimpanzee amplicon sequences show remarkable difference relative to the remainder of the genome, a difference that may be the result of intense selective pressure on male fertility. Four subspecies of common chimpanzees have undergone extended periods of isolation and appear to be in the early process of subspeciation. A recent study found amplicons enriched for testis-expressed genes on the primate X-chromosome the target of hard selective sweeps, and male-fertility genes on the Y-chromosome may also be the targets of selection. However, little is understood about Y-chromosome amplicon diversity within and across chimpanzee populations. Here, we analyze nine common chimpanzee (representing three subspecies: Pan troglodytes schweinfurthii, Pan troglodytes ellioti, and Pan troglodytes verus) and two bonobo (Pan paniscus) male whole-genome sequences to assess Y ampliconic copy-number diversity across the Pan genus. We observe that the copy number of Y chromosome amplicons is variable among chimpanzees and bonobos, and identify several lineage-specific patterns, including variable copy number of azoospermia candidates RBMY and DAZ. We detect recurrent switchpoints of copy-number change along the ampliconic tracts across chimpanzee populations, which may be the result of localized genome instability or selective forces.

Multigenomic Delineation of Plasmodium Species of the Laverania Subgenus Infecting Wild-Living Chimpanzees and Gorillas

Plasmodium falciparum, the major cause of malaria morbidity and mortality worldwide, is only distantly related to other human malaria parasites and has thus been placed in a separate subgenus, termed Laverania Parasites morphologically similar to P. falciparum have been identified in African apes, but only one other Laverania species, Plasmodium reichenowi from chimpanzees, has been formally described. Although recent studies have pointed to the existence of additional Laverania species, their precise number and host associations remain uncertain, primarily because of limited sampling and a paucity of parasite sequences other than from mitochondrial DNA. To address this, we used limiting dilution polymerase chain reaction to amplify additional parasite sequences from a large number of chimpanzee and gorilla blood and fecal samples collected at two sanctuaries and 30 field sites across equatorial Africa. Phylogenetic analyses of more than 2,000 new sequences derived from the mitochondrial, nuclear, and apicoplast genomes revealed six divergent and well-supported clades within the Laverania parasite group. Although two of these clades exhibited deep subdivisions in phylogenies estimated from organelle gene sequences, these sublineages were geographically defined and not present in trees from four unlinked nuclear loci. This greatly expanded sequence data set thus confirms six, and not seven or more, ape Laverania species, of which P. reichenowi, Plasmodium gaboni, and Plasmodium billcollinsi only infect chimpanzees, whereas Plasmodium praefalciparum, Plasmodium adleri, and Pladmodium blacklocki only infect gorillas. The new sequence data also confirm the P. praefalciparum origin of human P. falciparum.

How did bonobos come to range south of the congo river? Reconsideration of the divergence of Pan paniscus from other Pan populations

While investigating the genetic structure in wild bonobos,(1) we realized that the widely accepted scenario positing that the Pleistocene appearance of the Congo River separated the common ancestor of chimpanzees (Pan troglodytes) and bonobos (P. paniscus) into two species is not supported by recent geographical knowledge about the formation of the Congo River. We explored the origin of bonobos using a broader biogeographical perspective by examining local faunas in the central African region. The submarine Congo River sediments and paleotopography of central Africa show that the Congo River has functioned as a geographical barrier for the last 34 million years. This evidence allows us to hypothesize that when the river was first formed, the ancestor of bonobos did not inhabit the current range of the species on the left bank of the Congo River but that, during rare times when the Congo River discharge decreased during the Pleistocene, one or more founder populations of ancestral Pan paniscus crossed the river to its left bank. The proposed scenario for formation of the Congo River and the corridor hypothesis for an ancestral bonobo population is key to understanding the distribution of great apes and their evolution.

Reintroduction of confiscated and displaced mammals risks outbreeding and introgression in natural populations, as evidenced by orang-utans of divergent subspecies

Confiscated and displaced mammals are often taken to sanctuaries, where the explicit goal may be reintroduction to the wild. By inadvertently collecting animals from different source populations, however, such efforts risk reintroducing individuals that have not been in genetic contact for significant periods of time. Using genetic analyses and 44 years of data from Camp Leakey, an orang-utan rehabilitation site on Borneo, we determined the minimum extent to which orang-utans representing non-native, geographically and reproductively isolated taxa were reintroduced into the surrounding wild population. We found two reintroduced females were from a non-native subspecies, and have since produced at least 22 hybridized and introgressed descendants to date, of which at least 15 are living. Given that Bornean orang-utan subspecies are thought to have diverged from a common ancestor around 176,000 years ago, with marked differentiation over the last 80,000 years, we highlight the need for further evaluation of the effects of hybridizing orang-utans of different taxa--particularly in light of the ~1500 displaced orang-utans awaiting urgent reintroduction. As endangered mammals are increasing in number in sanctuaries worldwide, we stress the need for re-examination of historical reintroductions, to assess the extent and effects of de facto translocations in the past.

Visualizing spatial population structure with estimated effective migration surfaces

Genetic data often exhibit patterns broadly consistent with 'isolation by distance'-a phenomenon where genetic similarity decays with geographic distance. In a heterogeneous habitat, this may occur more quickly in some regions than in others: for example, barriers to gene flow can accelerate differentiation between neighboring groups. We use the concept of 'effective migration' to model the relationship between genetics and geography. In this paradigm, effective migration is low in regions where genetic similarity decays quickly. We present a method to visualize variation in effective migration across a habitat from geographically indexed genetic data. Our approach uses a population genetic model to relate effective migration rates to expected genetic dissimilarities. We illustrate its potential and limitations using simulations and data from elephant, human and Arabidopsis thaliana populations. The resulting visualizations highlight important spatial features of population structure that are difficult to discern using existing methods for summarizing genetic variation.

Inference of purifying and positive selection in three subspecies of chimpanzees (Pan troglodytes) from exome sequencing

We study genome-wide nucleotide diversity in three subspecies of extant chimpanzees using exome capture. After strict filtering, Single Nucleotide Polymorphisms and indels were called and genotyped for greater than 50% of exons at a mean coverage of 35× per individual. Central chimpanzees (Pan troglodytes troglodytes) are the most polymorphic (nucleotide diversity, θ_w = 0.0023 per site) followed by Eastern (P. t. schweinfurthii) chimpanzees (θ_w = 0.0016) and Western (P. t. verus) chimpanzees (θ_w = 0.0008). A demographic scenario of divergence without gene flow fits the patterns of autosomal synonymous nucleotide diversity well except for a signal of recent gene flow from Western into Eastern chimpanzees. The striking contrast in X-linked versus autosomal polymorphism and divergence previously reported in Central chimpanzees is also found in Eastern and Western chimpanzees. We show that the direction of selection statistic exhibits a strong nonmonotonic relationship with the strength of purifying selection S, making it inappropriate for estimating S. We instead use counts in synonymous versus nonsynonymous frequency classes to infer the distribution of S coefficients acting on nonsynonymous mutations in each subspecies. The strength of purifying selection we infer is congruent with the differences in effective sizes of each subspecies: Central chimpanzees are undergoing the strongest purifying selection followed by Eastern and Western chimpanzees. Coding indels show stronger selection against indels changing the reading frame than observed in human populations.

Environmental variation and rivers govern the structure of chimpanzee genetic diversity in a biodiversity hotspot

Background

The mechanisms that underlie the diversification of tropical animals remain poorly understood, but new approaches that combine geo-spatial modeling with spatially explicit genetic data are providing fresh insights on this topic. Data about the diversification of tropical mammals remain particularly sparse, and vanishingly few opportunities exist to study endangered large mammals that increasingly exist only in isolated pockets. The chimpanzees of Cameroon represent a unique opportunity to examine the mechanisms that promote genetic differentiation in tropical mammals because the region is home to two chimpanzee subspecies: Pan troglodytes ellioti and P. t. trogolodytes. Their ranges converge in central Cameroon, which is a geographically, climatically and environmentally complex region that presents an unparalleled opportunity to examine the roles of rivers and/or environmental variation in influencing the evolution of chimpanzee populations.

Results

We analyzed microsatellite genotypes and mtDNA HVRI sequencing data from wild chimpanzees sampled at a fine geographic scale across Cameroon and eastern Nigeria using a spatially explicit approach based upon Generalized Dissimilarity Modeling. Both the Sanaga River and environmental variation were found to contribute to driving separation of the subspecies. The importance of environmental variation differed among subspecies. Gene-environment associations were weak in P. t. troglodytes, whereas environmental variation was found to play a much larger role in shaping patterns of genetic differentiation in P. t. ellioti.

Conclusions

We found that both the Sanaga River and environmental variation likely play a role in shaping patterns of chimpanzee genetic diversity. Future studies using single nucleotide polymorphism (SNP) data are necessary to further understand how rivers and environmental variation contribute to shaping patterns of genetic variation in chimpanzees.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0274-0) contains supplementary material, which is available to authorized users.

Chimpanzee population structure in Cameroon and Nigeria is associated with habitat variation that may be lost under climate change

Background

The Nigeria-Cameroon chimpanzee (Pan troglodytes ellioti) is found in the Gulf of Guinea biodiversity hotspot located in western equatorial Africa. This subspecies is threatened by habitat fragmentation due to logging and agricultural development, hunting for the bushmeat trade, and possibly climate change. Although P. t. ellioti appears to be geographically separated from the neighboring central chimpanzee (P. t. troglodytes) by the Sanaga River, recent population genetics studies of chimpanzees from across this region suggest that additional factors may also be important in their separation. The main aims of this study were: 1) to model the distribution of suitable habitat for P. t. ellioti across Cameroon and Nigeria, and P. t. troglodytes in southern Cameroon, 2) to determine which environmental factors best predict their optimal habitats, and 3) to compare modeled niches and test for their levels of divergence from one another. A final aim of this study was to examine the ways that climate change might impact suitable chimpanzee habitat across the region under various scenarios.

Results

Ecological niche models (ENMs) were created using the software package Maxent for the three populations of chimpanzees that have been inferred to exist in Cameroon and eastern Nigeria: (i) P. t. troglodytes in southern Cameroon, (ii) P. t. ellioti in northwestern Cameroon, and (iii) P. t. ellioti in central Cameroon. ENMs for each population were compared using the niche comparison test in ENMtools, which revealed complete niche divergence with very little geographic overlap of suitable habitat between populations.

Conclusions

These findings suggest that a positive relationship may exist between environmental variation and the partitioning of genetic variation found in chimpanzees across this region. ENMs for each population were also projected under three different climate change scenarios for years 2020, 2050, and 2080. Suitable habitat of P. t. ellioti in northwest Cameroon / eastern Nigeria is expected to remain largely unchanged through 2080 in all considered scenarios. In contrast, P. t. ellioti in central Cameroon, which represents half of the population of this subspecies, is expected to experience drastic reductions in its ecotone habitat over the coming century.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0275-z) contains supplementary material, which is available to authorized users.