Does Effective Population Size Govern Evolutionary Differences in Telomere Length?

Lengths of telomeres vary by an order of magnitude across mammalian species. Similarly, age- and sex-standardized telomere lengths differ by up to 1 kb (14%) across human populations. How to explain these differences? Telomeres play a central role in senescence and aging, and genes that affect telomere length are likely under weak selection (i.e. telomere length is a trait that is subject to nearly neutral evolution). Importantly, natural selection is more effective in large populations than in small populations. Here, we propose that observed differences in telomere length across species and populations are largely due to differences in effective population sizes. In this perspective, we present preliminary evolutionary genetic evidence supporting this hypothesis and highlight the need for more data.

Orangutan males make increased use of social learning opportunities, when resource availability is high

Humans' colonization of diverse habitats relied on our ancestors' abilities to innovate and share innovations with others. While ecological impacts on innovations are well studied, their effect on social learning remains poorly understood. We examined how food availability affects social learning in migrant orangutan unflanged males, who may learn from local orangutans through peering (i.e., observational social learning). We analyzed 1,384 dyadic associations, including 360 peering events, among 46 wild Sumatran orangutan and 25 Bornean orangutan males, collected over 18 years. Migrants' peering rates significantly increased with higher food availability and time spent in proximity to others. Furthermore, migrants in the more sociable Sumatran population exhibited significantly higher peering rates compared to the Borneans, suggesting intrinsic and/or developmental effects of food availability on social learning. These findings emphasize the importance of investigating ecological effects on social learning on the immediate, developmental, and intrinsic levels for our understanding of cultural evolution.

Geography vs. past climate: the drivers of population genetic structure of the Himalayan langur

Background

Contemporary species distribution, genetic diversity and evolutionary history in many taxa are shaped by both historical and current climate as well as topography. The Himalayas show a huge variation in topography and climatic conditions across its entire range, and have experienced major climatic fluctuations in the past. However, very little is known regarding how this heterogenous landscape has moulded the distribution of Himalayan fauna. A recent study examined the effect of these historical events on the genetic diversity of the Himalayan langurs in Nepal Himalaya. However, this study did not include the samples from the Indian Himalayan region (IHR). Therefore, here we revisit the questions addressed in the previous study with a near complete sampling from the IHR, along with the samples from the Nepal Himalaya. We used the mitochondrial Cytochrome-b (Cyt-b, 746 bp) region combined with multiple phylogeographic analyses and palaeodistribution modelling.

Results

Our dataset contained 144 sequences from the IHR as well as the Nepal Himalaya. Phylogenetic analysis showed a low divergent western clade nested within high divergent group of eastern lineages and in the network analysis we identified 22 haplotypes over the entire distribution range of the Himalayan langurs. Samples from the Nepal Himalaya showed geographically structured haplotypes corresponding to different river barriers, whereas samples from IHR showed star-like topology with no structure. Our statistical phylogeography analysis using diyABC supported the model of east to west colonisation of these langurs with founder event during colonisation. Analysis of demographic history showed that the effective population size of the Himalayan langurs decreased at the onset of last glacial maximum (LGM) and started increasing post LGM. The palaeodistribution modelling showed that the extent of suitable habitat shifted from low elevation central Nepal, and adjoining parts of north India, during LGM to the western Himalaya at present.

Conclusion

The current genetic diversity and distribution of Himalayan langurs in the Nepal Himalaya has been shaped by river barriers, whereas the rivers in the IHR had relatively less time to act as a strong genetic barrier after the recent colonisation event. Further, the post LGM expansion could have had confounding effect on Himalayan langur population structure in both Nepal Himalaya and IHR.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02054-1.

“Engaging the Enemy”: Orangutan (Pongo pygmaeus morio) Conservation in Human Modified Environments in the Kinabatangan floodplain of Sabah, Malaysian Borneo

Throughout the equatorial tropics, forest conversion to agriculture often fragments crucial primate habitat. In 30 years, 80% of the alluvial lowland forests along the Kinabatangan River in Sabah, Malaysian Borneo, have been supplanted by oil palm (Elaeis guineensis) plantations. Today, only about 20% of the former orangutan (Pongo pygmaeus morio) population remains in the region. Because most of the land is now under the tenure of agribusiness companies, we used a pragmatic approach of mixed biosocial methods and citizen science engagement of oil palm growers (N = 6) as active conservation partners to study orangutan use of the privately administered landscape between protected forest fragments. We found that 22 of 25 remanent forest patches (0.5 to 242 hectares) surveyed within plantations contained food or shelter resources useful for orangutans. Of these, 20 are in regular transitory use by wider-ranging adult male orangutans, and in 9 patches, females are resident and raising offspring isolated within oil palm plantations. These findings indicate that orangutans retain a measure of normal metapopulation dynamics necessary for viability at the landscape level despite drastic habitat modification. We found that barriers to in situ conservation in these agroforest matrices were due to the following misconceptions across sectors: 1) Good farming practices require exclusion of wildlife; 2) Orangutans seen in plantations must be “rescued” by people; and 3) Translocation is an appropriate conservation strategy, and nondetrimental to orangutans. Our exploratory study exemplifies the value of biosocial methods and collaboration with industrial-scale farmers to support primate resilience in forests fragmented by agriculture.

Short Tandem Repeats as a High-Resolution Marker for Capturing Recent Orangutan Population Evolution

The genus Pongo is ideal to study population genetics adaptation, given its remarkable phenotypic divergence and the highly contrasting environmental conditions it’s been exposed to. Studying its genetic variation bears the promise to reveal a motion picture of these great apes’ evolutionary and adaptive history, and also helps us expand our knowledge of the patterns of adaptation and evolution. In this work, we advance the understanding of the genetic variation among wild orangutans through a genome-wide study of short tandem repeats (STRs). Their elevated mutation rate makes STRs ideal markers for the study of recent evolution within a given population. Current technological and algorithmic advances have rendered their sequencing and discovery more accurate, therefore their potential can be finally leveraged in population genetics studies. To study patterns of population variation within the wild orangutan population, we genotyped the short tandem repeats in a population of 21 individuals spanning four Sumatran and Bornean (sub-) species and eight Southeast Asian regions. We studied the impact of sequencing depth on our ability to genotype STRs and found that the STR copy number changes function as a powerful marker, correctly capturing the demographic history of these populations, even the divergences as recent as 10 Kya. Moreover, gene ontology enrichments for genes close to STR variants are aligned with local adaptations in the two islands. Coupled with more advanced STR-compatible population models, and selection tests, genomic studies based on STRs will be able to reduce the gap caused by the missing heritability for species with recent adaptations.

Comparative population genetics of the invasive mosquito Aedes albopictus and the native mosquito Aedes flavopictus in the Korean peninsula

BACKGROUND

Aedes mosquitoes are important invasive species contributing to the spread of chikungunya, dengue fever, yellow fever, zika virus, and other dangerous vector-borne diseases. Aedes albopictus is native to southeast Asia, with rapid expansion due to human activity, showing a wide distribution in the Korean peninsula. Aedes flavopictus is considered to be native to East Asia, with a broad distribution in the region, including the Korean peninsula. A better understanding of the genetic diversity of these species is critical for establishing strategies for disease prevention and vector control.

METHODS

We obtained DNA from 148 specimens of Ae. albopictus and 166 specimens of Ae. flavopictus in Korea, and amplified two mitochondrial genes (COI and ND5) to compare the genetic diversity and structure of the two species.

RESULTS

We obtained a 658-bp sequence of COI and a 423-bp sequence of ND5 from both mosquito species. We found low diversity and a nonsignificant population genetic structure in Ae. albopictus, and high diversity and a nonsignificant structure in Ae. flavopictus for these two mitochondrial genes. Aedes albopictus had fewer haplotypes with respect to the number of individuals, and a slight mismatch distribution was confirmed. By contrast, Ae. flavopictus had a large number of haplotypes compared with the number of individuals, and a large unimodal-type mismatch distribution was confirmed. Although the genetic structure of both species was nonsignificant, Ae. flavopictus exhibited higher genetic diversity than Ae. albopictus.

CONCLUSIONS

Aedes albopictus appears to be an introduced species, whereas Ae. flavopictus is endemic to the Korean peninsula, and the difference in genetic diversity between the two species is related to their adaptability and introduction history. Further studies on the genetic structure and diversity of these mosquitos will provide useful data for vector control.

'Emptying Forests?' Conservation Implications of Past Human-Primate Interactions

Non-human primates are among the most vulnerable tropical animals to extinction and ~50% of primate species are endangered. Human hunting is considered a major cause of increasingly 'empty forests', yet archaeological data remains under-utilised in testing this assertion over the longer-term. Zooarchaeological datasets allow investigation of human exploitation of primates and the reconstruction of extinction, extirpation, and translocation processes. We evaluate the application and limitations of data from zooarchaeological studies spanning the past 45 000 years in South and Southeast Asia in guiding primate conservation efforts. We highlight that environmental change was the primary threat to many South and Southeast Asian non-human primate populations during much of the Holocene, foreshadowing human-induced land-use and environmental change as major threats of the 21st century.

The historical range and drivers of decline of the Tapanuli orangutan

The Tapanuli Orangutan (Pongo tapanuliensis) is the most threatened great ape species in the world. It is restricted to an area of about 1,000 km² of upland forest where fewer than 800 animals survive in three declining subpopulations. Through a historical ecology approach involving analysis of newspaper, journals, books and museum records from the early 1800s to 2009, we demonstrate that historically Pongo tapanuliensis inhabited a much larger area, and occurred across a much wider range of habitat types and at lower elevations than now. Its current Extent of Occurrence is 2.5% and 5.0% of the historical range in the 1890s and 1940s respectively. A combination of historical fragmentation of forest habitats, mostly for small-scale agriculture, and unsustainable hunting likely drove various populations to the south, east and west of the current population to extinction. This happened prior to the industrial-scale forest conversion that started in the 1970s. Our findings indicate how sensitive P. tapanuliensis is to the combined effects of habitat fragmentation and unsustainable take-off rates. Saving this species will require prevention of any further fragmentation and killings or other removal of animals from the remaining population. Without concerted action to achieve this, the remaining populations of P. tapanuliensis are doomed to become extinct within several orangutan generations.

Genomic targets for high-resolution inference of kinship, ancestry and disease susceptibility in orang-utans (genus: Pongo)

BACKGROUND

Orang-utans comprise three critically endangered species endemic to the islands of Borneo and Sumatra. Though whole-genome sequencing has recently accelerated our understanding of their evolutionary history, the costs of implementing routine genome screening and diagnostics remain prohibitive. Capitalizing on a tri-fold locus discovery approach, combining data from published whole-genome sequences, novel whole-exome sequencing, and microarray-derived genotype data, we aimed to develop a highly informative gene-focused panel of targets that can be used to address a broad range of research questions.

RESULTS

We identified and present genomic co-ordinates for 175,186 SNPs and 2315 Y-chromosomal targets, plus 185 genes either known or presumed to be pathogenic in cardiovascular (N = 109) or respiratory (N = 43) diseases in humans - the primary and secondary causes of captive orang-utan mortality - or a majority of other human diseases (N = 33). As proof of concept, we designed and synthesized 'SeqCap' hybrid capture probes for these targets, demonstrating cost-effective target enrichment and reduced-representation sequencing.

CONCLUSIONS

Our targets are of broad utility in studies of orang-utan ancestry, admixture and disease susceptibility and aetiology, and thus are of value in addressing questions key to the survival of these species. To facilitate comparative analyses, these targets could now be standardized for future orang-utan population genomic studies. The targets are broadly compatible with commercial target enrichment platforms and can be utilized as published here to synthesize applicable probes.

Genetic variation and phylogeographic structure of Spodoptera exigua in western China based on mitochondrial DNA and microsatellite markers

The beet armyworm, Spodoptera exigua, is a significant agricultural pest of numerous crops and has caused serious economic losses in China. To effectively control this pest, we analyzed its genetic variation, population genetic structure and demographic history. We used mitochondrial DNA (mtDNA) fragments of the cytochrome oxidase subunit I (COI) and eight nuclear microsatellite loci to investigate genetic diversity and population genetic structure of S. exigua populations at 14 sampling sites in western China. Both mtDNA and microsatellite data indicated low levels of genetic diversity among all populations. A moderate genetic differentiation among some S. exigua populations was detected. Neighbor-joining dendrograms, STRUCTURE, and principal coordinate analysis (PCoA) revealed two genetically distinct groups: the KEL group and the remaining population group. Isolation by distance (IBD) results showed a weak significant correlation between geographic distance and genetic differentiation. Haplotype networks, neutrality testing, and mismatch distribution analysis indicated that the beet armyworm experienced a recent rapid expansion without a recent genetic bottleneck in western China. Thus, the results of this population genetic study can help with the development of strategies for managing this highly migratory pest.

The slow ape: High infant survival and long interbirth intervals in wild orangutans

Orangutans (Pongo spp.) are reported to have extremely slow life histories, including the longest average interbirth intervals of all mammals. Such slow life history can be viable only when unavoidable mortality is kept low. Thus, orangutans' survivorship under natural conditions is expected to be extremely high. Previous estimates of orangutan life history were based on captive individuals living under very different circumstances or on small samples from wild populations. Here, we combine birth data from seven field sites, each with demographic data collection for at least 10 years (range 12-43 years) on wild orangutans to better document their life history. Using strict criteria for data inclusion, we calculated infant survival, interbirth intervals and female age at first reproduction, across species, subspecies and islands. We found an average closed interbirth interval of 7.6 years, as well as consistently very high pre-weaning survival for males and females. Female survival of 94% until age at first birth (at around age 15 years) was higher than reported for any other mammal species under natural conditions. Similarly, annual survival among parous females is very high, but longevity remains to be estimated. Current data suggest no major life history differences between Sumatran and Bornean orangutans. The high offspring survival is remarkable, noting that modern human populations seem to have reached the same level of survival only in the 20th century. The orangutans' slow life history illustrates what can be achieved if a hominoid bauplan is exposed to low unavoidable mortality. Their high survival is likely due to their arboreal and non-gregarious lifestyle, and has allowed them to maintain viable populations, despite living in low-productivity habitats. However, their slow life history also implies that orangutans are highly vulnerable to a catastrophic population crash in the face of drastic habitat change.

Hominoid intraspecific cranial variation mirrors neutral genetic diversity

In humans, patterns of cranial variation mirror genetic diversity globally, implicating population history as a key driver of cranial disparity. Here, we demonstrate that the magnitude of genetic diversity within 12 extant ape taxa explains a large proportion of cranial shape variation. Taxa that are more genetically diverse tend to be more cranially diverse also. Our results suggest that neutral evolutionary processes such as mutation, genetic drift, and gene flow are reflected in both genetic and cranial diversity in apes. This work provides a perspective on intraspecific cranial variation in apes which has important implications for interpreting selective and developmental pressures on the cranium and for understanding shape variation in fossil hominin crania.

Natural selection, developmental constraint, and plasticity have all been invoked as explanations for intraspecific cranial variation in humans and apes. However, global patterns of human cranial variation are congruent with patterns of genetic variation, demonstrating that population history has influenced cranial variation in humans. Here we show that this finding is not unique to Homo sapiens but is also broadly evident across extant ape species. Specifically, taxa that exhibit greater intraspecific cranial shape variation also exhibit greater genetic diversity at neutral autosomal loci. Thus, cranial shape variation within hominoid taxa reflects the population history of each species. Our results suggest that neutral evolutionary processes such as mutation, gene flow, and genetic drift have played an important role in generating cranial variation within species. These findings are consistent with previous work on human cranial morphology and improve our understanding of the evolutionary processes that generate intraspecific cranial shape diversity within hominoids. This work has implications for the analysis of selective and developmental pressures on the cranium and for interpreting shape variation in fossil hominin crania.

Genetic and genomic monitoring with minimally invasive sampling methods

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

Genetic analysis suggests dispersal among chimpanzees in a fragmented forest landscape in Uganda

Habitat fragmentation is a leading threat to global biodiversity. Dispersal plays a key role in gene flow and population viability, but the impact of fragmentation on dispersal patterns remains poorly understood. Among chimpanzees, males typically remain in their natal communities while females often disperse. However, habitat loss and fragmentation may cause severe ecological disruptions, potentially resulting in decreased fitness benefits of male philopatry and limited female dispersal ability. To investigate this issue, we genotyped nearly 900 non-invasively collected chimpanzee fecal samples across a fragmented forest habitat that may function as a corridor between two large continuous forests in Uganda, and used the spatial associations among co-sampled genotypes to attribute a total of 229 individuals to 10 distinct communities, including 9 communities in the corridor habitat and 1 in continuous forest. We then used parentage analyses to infer instances of between-community dispersal. Of the 115 parent-offspring dyads detected with confidence, members of 39% (N = 26) of mother-daughter dyads were found in different communities, while members of 10% (N = 5) of father-son dyads were found in different communities. We also found direct evidence for one dispersal event that occurred during the study period, as a female's sample found first in one community was found multiple times in another community 19 months later. These findings suggest that even in fragmented habitats, chimpanzee males remain in their natal communities while females tend to disperse. Corridor enhancement in unprotected forest fragments could help maintain gene flow in chimpanzees and other species amid anthropogenic pressures.

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.

Morphometric, Behavioral, and Genomic Evidence for a New Orangutan Species

Six extant species of non-human great apes are currently recognized: Sumatran and Bornean orangutans, eastern and western gorillas, and chimpanzees and bonobos [1]. However, large gaps remain in our knowledge of fine-scale variation in hominoid morphology, behavior, and genetics, and aspects of great ape taxonomy remain in flux. This is particularly true for orangutans (genus: Pongo), the only Asian great apes and phylogenetically our most distant relatives among extant hominids [1]. Designation of Bornean and Sumatran orangutans, P. pygmaeus (Linnaeus 1760) and P. abelii (Lesson 1827), as distinct species occurred in 2001 [1, 2]. Here, we show that an isolated population from Batang Toru, at the southernmost range limit of extant Sumatran orangutans south of Lake Toba, is distinct from other northern Sumatran and Bornean populations. By comparing cranio-mandibular and dental characters of an orangutan killed in a human-animal conflict to those of 33 adult male orangutans of a similar developmental stage, we found consistent differences between the Batang Toru individual and other extant Ponginae. Our analyses of 37 orangutan genomes provided a second line of evidence. Model-based approaches revealed that the deepest split in the evolutionary history of extant orangutans occurred ∼3.38 mya between the Batang Toru population and those to the north of Lake Toba, whereas both currently recognized species separated much later, about 674 kya. Our combined analyses support a new classification of orangutans into three extant species. The new species, Pongo tapanuliensis, encompasses the Batang Toru population, of which fewer than 800 individuals survive. VIDEO ABSTRACT.

Bayesian inference reveals ancient origin of simian foamy virus in orangutans

Simian foamy viruses (SFVs) infect most nonhuman primate species and appears to co-evolve with its hosts. This co-evolutionary signal is particularly strong among great apes, including orangutans (genus Pongo). Previous studies have identified three distinct orangutan SFV clades. The first of these three clades is composed of SFV from P. abelii from Sumatra, the second consists of SFV from P. pygmaeus from Borneo, while the third clade is mixed, comprising an SFV strain found in both species of orangutan. The existence of the mixed clade has been attributed to an expansion of P. pygmaeus into Sumatra following the Mount Toba super-volcanic eruption about 73,000years ago. Divergence dating, however, has yet to be performed to establish a temporal association with the Toba eruption. Here, we use a Bayesian framework and a relaxed molecular clock model with fossil calibrations to test the Toba hypothesis and to gain a more complete understanding of the evolutionary history of orangutan SFV. As with previous studies, our results show a similar three-clade orangutan SFV phylogeny, along with strong statistical support for SFV-host co-evolution in orangutans. Using Bayesian inference, we date the origin of orangutan SFV to >4.7 million years ago (mya), while the mixed species clade dates to approximately 1.7mya, >1.6 million years older than the Toba super-eruption. These results, combined with fossil and paleogeographic evidence, suggest that the origin of SFV in Sumatran and Bornean orangutans, including the mixed species clade, likely occurred on the mainland of Indo-China during the Late Pliocene and Calabrian stage of the Pleistocene, respectively.

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.

Effective Characterisation of the Complete Orang-Utan Mitochondrial DNA Control Region, in the Face of Persistent Focus in Many Taxa on Shorter Hypervariable Regions

The hypervariable region I (HVRI) is persistently used to discern haplotypes, to distinguish geographic subpopulations, and to infer taxonomy in a range of organisms. Numerous studies have highlighted greater heterogeneity elsewhere in the mitochondrial DNA control region, however-particularly, in some species, in other understudied hypervariable regions. To assess the abundance and utility of such potential variations in orang-utans, we characterised 36 complete control-region haplotypes, of which 13 were of Sumatran and 23 of Bornean maternal ancestry, and compared polymorphisms within these and within shorter HVRI segments predominantly analysed in prior phylogenetic studies of Sumatran (~385 bp) and Bornean (~323 bp) orang-utans. We amplified the complete control region in a single PCR that proved successful even with highly degraded, non-invasive samples. By using species-specific primers to produce a single large amplicon (~1600 bp) comprising flanking coding regions, our method also serves to better avoid amplification of nuclear mitochondrial insertions (numts). We found the number, length and position of hypervariable regions is inconsistent between orang-utan species, and that prior definitions of the HVRI were haphazard. Polymorphisms occurring outside the predominantly analysed segments were phylogeographically informative in isolation, and could be used to assign haplotypes to comparable clades concordant with geographic subpopulations. The predominantly analysed segments could discern only up to 76% of all haplotypes, highlighting the forensic utility of complete control-region sequences. In the face of declining sequencing costs and our proven application to poor-quality DNA extracts, we see no reason to ever amplify only specific 'hypervariable regions' in any taxa, particularly as their lengths and positions are inconsistent and cannot be reliably defined-yet this strategy predominates widely. Given their greater utility and consistency, we instead advocate analysis of complete control-region sequences in future studies, where any shorter segment might otherwise have proven the region of choice.

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.

Genetic diversity in humans and non-human primates and its evolutionary consequences

Genetic diversity is a key parameter in population genetics and is important for understanding the process of evolution and for the development of appropriate conservation strategies. Recent advances in sequencing technology have enabled the measurement of genetic diversity of various organisms at the nucleotide level and on a genome-wide scale, yielding more precise estimates than were previously achievable. In this review, I have compiled and summarized the estimates of genetic diversity in humans and non-human primates based on recent genome-wide studies. Although studies on population genetics demonstrated fluctuations in population sizes over time, general patterns have emerged. As shown previously, genetic diversity in humans is one of the lowest among primates; however, certain other primate species exhibit genetic diversity that is comparable to or even lower than that in humans. There exists greater than 10-fold variation in genetic diversity among primate species, and I found weak correlation with species fecundity but not with body or propagule size. I further discuss the potential evolutionary consequences of population size decline on the evolution of primate species. The level of genetic diversity negatively correlates with the ratio of non-synonymous to synonymous polymorphisms in a population, suggesting that proportionally greater numbers of slightly deleterious mutations segregate in small rather than large populations. Although population size decline is likely to promote the fixation of slightly deleterious mutations, there are molecular mechanisms, such as compensatory mutations at various molecular levels, which may prevent fitness decline at the population level. The effects of slightly deleterious mutations from theoretical and empirical studies and their relevance to conservation biology are also discussed in this review.

Prevalence of Cryptosporidium spp., Enterocytozoon bieneusi, Encephalitozoon spp. and Giardia intestinalis in Wild, Semi-Wild and Captive Orangutans (Pongo abelii and Pongo pygmaeus) on Sumatra and Borneo, Indonesia

Background

Orangutans are critically endangered primarily due to loss and fragmentation of their natural habitat. This could bring them into closer contact with humans and increase the risk of zoonotic pathogen transmission.

Aims

To describe the prevalence and diversity of Cryptosporidium spp., microsporidia and Giardia intestinalis in orangutans at seven sites on Sumatra and Kalimantan, and to evaluate the impact of orangutans’ habituation and location on the occurrence of these zoonotic protists.

Result

The overall prevalence of parasites in 298 examined animals was 11.1%. The most prevalent microsporidia was Encephalitozoon cuniculi genotype II, found in 21 animals (7.0%). Enterocytozoon bieneusi genotype D (n = 5) and novel genotype Pongo 2 were detected only in six individuals (2.0%). To the best of our knowledge, this is the first report of these parasites in orangutans. Eight animals were positive for Cryptosporidium spp. (2.7%), including C. parvum (n = 2) and C. muris (n = 6). Giardia intestinalis assemblage B, subtype MB6, was identified in a single individual. While no significant differences between the different human contact level groups (p = 0.479–0.670) or between the different islands (p = 0.992) were reported in case of E. bieneusi or E. cuniculi, Cryptosporidium spp. was significantly less frequently detected in wild individuals (p < 2×10⁻¹⁶) and was significantly more prevalent in orangutans on Kalimantan than on Sumatra (p < 2×10⁻¹⁶).

Conclusion

Our results revealed that wild orangutans are significantly less frequently infected by Cryptosporidium spp. than captive and semi-wild animals. In addition, this parasite was more frequently detected at localities on Kalimantan. In contrast, we did not detect any significant difference in the prevalence of microsporidia between the studied groups of animals. The sources and transmission modes of infections were not determined, as this would require repeated sampling of individuals, examination of water sources, and sampling of humans and animals sharing the habitat with orangutans.

Land-cover changes predict steep declines for the Sumatran orangutan (Pongo abelii)

Positive news about Sumatran orangutans is rare. The species is critically endangered because of forest loss and poaching, and therefore, determining the impact of future land-use change on this species is important. To date, the total Sumatran orangutan population has been estimated at 6600 individuals. On the basis of new transect surveys, we estimate a population of 14,613 in 2015. This higher estimate is due to three factors. First, orangutans were found at higher elevations, elevations previously considered outside of their range and, consequently, not surveyed previously. Second, orangutans were found more widely distributed in logged forests. Third, orangutans were found in areas west of the Toba Lake that were not previously surveyed. This increase in numbers is therefore due to a more wide-ranging survey effort and is not indicative of an increase in the orangutan population in Sumatra. There are evidently more Sumatran orangutans remaining in the wild than we thought, but the species remains under serious threat. Current scenarios for future forest loss predict that as many as 4500 individuals could vanish by 2030. Despite the positive finding that the population is double the size previously estimated, our results indicate that future deforestation will continue to be the cause of rapid declines in orangutan numbers. Hence, we urge that all developmental planning involving forest loss be accompanied by appropriate environmental impact assessments conforming with the current national and provincial legislations, and, through these, implement specific measures to reduce or, better, avoid negative impacts on forests where orangutans occur.

Great ape Y Chromosome and mitochondrial DNA phylogenies reflect subspecies structure and patterns of mating and dispersal

The distribution of genetic diversity in great ape species is likely to have been affected by patterns of dispersal and mating. This has previously been investigated by sequencing autosomal and mitochondrial DNA (mtDNA), but large-scale sequence analysis of the male-specific region of the Y Chromosome (MSY) has not yet been undertaken. Here, we use the human MSY reference sequence as a basis for sequence capture and read mapping in 19 great ape males, combining the data with sequences extracted from the published whole genomes of 24 additional males to yield a total sample of 19 chimpanzees, four bonobos, 14 gorillas, and six orangutans, in which interpretable MSY sequence ranges from 2.61 to 3.80 Mb. This analysis reveals thousands of novel MSY variants and defines unbiased phylogenies. We compare these with mtDNA-based trees in the same individuals, estimating time-to-most-recent common ancestor (TMRCA) for key nodes in both cases. The two loci show high topological concordance and are consistent with accepted (sub)species definitions, but time depths differ enormously between loci and (sub)species, likely reflecting different dispersal and mating patterns. Gorillas and chimpanzees/bonobos present generally low and high MSY diversity, respectively, reflecting polygyny versus multimale–multifemale mating. However, particularly marked differences exist among chimpanzee subspecies: The western chimpanzee MSY phylogeny has a TMRCA of only 13.2 (10.8–15.8) thousand years, but that for central chimpanzees exceeds 1 million years. Cross-species comparison within a single MSY phylogeny emphasizes the low human diversity, and reveals species-specific branch length variation that may reflect differences in long-term generation times.

Reconstructing the demographic history of orang-utans using Approximate Bayesian Computation

Investigating how different evolutionary forces have shaped patterns of DNA variation within and among species requires detailed knowledge of their demographic history. Orang-utans, whose distribution is currently restricted to the South-East Asian islands of Borneo (Pongo pygmaeus) and Sumatra (Pongo abelii), have likely experienced a complex demographic history, influenced by recurrent changes in climate and sea levels, volcanic activities and anthropogenic pressures. Using the most extensive sample set of wild orang-utans to date, we employed an Approximate Bayesian Computation (ABC) approach to test the fit of 12 different demographic scenarios to the observed patterns of variation in autosomal, X-chromosomal, mitochondrial and Y-chromosomal markers. In the best-fitting model, Sumatran orang-utans exhibit a deep split of populations north and south of Lake Toba, probably caused by multiple eruptions of the Toba volcano. In addition, we found signals for a strong decline in all Sumatran populations ~24 ka, probably associated with hunting by human colonizers. In contrast, Bornean orang-utans experienced a severe bottleneck ~135 ka, followed by a population expansion and substructuring starting ~82 ka, which we link to an expansion from a glacial refugium. We showed that orang-utans went through drastic changes in population size and connectedness, caused by recurrent contraction and expansion of rainforest habitat during Pleistocene glaciations and probably hunting by early humans. Our findings emphasize the fact that important aspects of the evolutionary past of species with complex demographic histories might remain obscured when applying overly simplified models.

Great ape genomics

The great ape families are the species most closely related to our own, comprising chimpanzees, bonobos, gorillas, and orangutans. They live exclusively in tropical rainforests in Central Africa and the islands of Southeast Asia. Due to their close evolutionary relationship with humans, great apes share many cognitive, physiological, and morphological similarities with humans. The members of the great ape family make obvious models to facilitate the further understanding about humans' biology and history. This review will discuss how the recent addition of genome-wide data from great apes has furthered humans' understanding of these species and humanity, especially in the realm of evolutionary genetics.

Generation of SNP datasets for orangutan population genomics using improved reduced-representation sequencing and direct comparisons of SNP calling algorithms

Background

High-throughput sequencing has opened up exciting possibilities in population and conservation genetics by enabling the assessment of genetic variation at genome-wide scales. One approach to reduce genome complexity, i.e. investigating only parts of the genome, is reduced-representation library (RRL) sequencing. Like similar approaches, RRL sequencing reduces ascertainment bias due to simultaneous discovery and genotyping of single-nucleotide polymorphisms (SNPs) and does not require reference genomes. Yet, generating such datasets remains challenging due to laboratory and bioinformatical issues. In the laboratory, current protocols require improvements with regards to sequencing homologous fragments to reduce the number of missing genotypes. From the bioinformatical perspective, the reliance of most studies on a single SNP caller disregards the possibility that different algorithms may produce disparate SNP datasets.

Results

We present an improved RRL (iRRL) protocol that maximizes the generation of homologous DNA sequences, thus achieving improved genotyping-by-sequencing efficiency. Our modifications facilitate generation of single-sample libraries, enabling individual genotype assignments instead of pooled-sample analysis. We sequenced ~1% of the orangutan genome with 41-fold median coverage in 31 wild-born individuals from two populations. SNPs and genotypes were called using three different algorithms. We obtained substantially different SNP datasets depending on the SNP caller. Genotype validations revealed that the Unified Genotyper of the Genome Analysis Toolkit and SAMtools performed significantly better than a caller from CLC Genomics Workbench (CLC). Of all conflicting genotype calls, CLC was only correct in 17% of the cases. Furthermore, conflicting genotypes between two algorithms showed a systematic bias in that one caller almost exclusively assigned heterozygotes, while the other one almost exclusively assigned homozygotes.

Conclusions

Our enhanced iRRL approach greatly facilitates genotyping-by-sequencing and thus direct estimates of allele frequencies. Our direct comparison of three commonly used SNP callers emphasizes the need to question the accuracy of SNP and genotype calling, as we obtained considerably different SNP datasets depending on caller algorithms, sequencing depths and filtering criteria. These differences affected scans for signatures of natural selection, but will also exert undue influences on demographic inferences. This study presents the first effort to generate a population genomic dataset for wild-born orangutans with known population provenance.

Genetic diversity and population structure of Sitodiplosis mosellana in Northern China

The wheat midge, Sitodiplosis mosellana, is an important pest in Northern China. We tested the hypothesis that the population structure of this species arises during a range expansion over the past 30 years. This study used microsatellite and mitochondrial loci to conduct population genetic analysis of S. mosellana across its distribution range in China. We found strong genetic structure among the 16 studied populations, including two genetically distinct groups (the eastern and western groups), broadly consistent with the geography and habitat fragmentation. These results underline the importance of natural barriers in impeding dispersal and gene flow of S. mosellana populations. Low to moderate genetic diversity among the populations and moderate genetic differentiation (F_ST = 0.117) between the two groups were also found. The populations in the western group had lower genetic diversity, higher genetic differentiation and lower gene flow (F_ST = 0.116, Nm = 1.89) than those in the eastern group (F_ST = 0.049, Nm = 4.91). Genetic distance between populations was positively and significantly correlated with geographic distance (r = 0.56, P<0.001). The population history of this species provided no evidence for population expansion or bottlenecks in any of these populations. Our data suggest that the distribution of genetic diversity, genetic differentiation and population structure of S. mosellana have resulted from a historical event, reflecting its adaptation to diverse habitats and forming two different gene pools. These results may be the outcome of a combination of restricted gene flow due to geographical and environmental factors, population history, random processes of genetic drift and individual dispersal patterns. Given the current risk status of this species in China, this study can offer useful information for forecasting outbreaks and designing effective pest management programs.

Population genomic analysis reveals a rich speciation and demographic history of orang-utans (Pongo pygmaeus and Pongo abelii)

To gain insights into evolutionary forces that have shaped the history of Bornean and Sumatran populations of orang-utans, we compare patterns of variation across more than 11 million single nucleotide polymorphisms found by previous mitochondrial and autosomal genome sequencing of 10 wild-caught orang-utans. Our analysis of the mitochondrial data yields a far more ancient split time between the two populations (~3.4 million years ago) than estimates based on autosomal data (0.4 million years ago), suggesting a complex speciation process with moderate levels of primarily male migration. We find that the distribution of selection coefficients consistent with the observed frequency spectrum of autosomal non-synonymous polymorphisms in orang-utans is similar to the distribution in humans. Our analysis indicates that 35% of genes have evolved under detectable negative selection. Overall, our findings suggest that purifying natural selection, genetic drift, and a complex demographic history are the dominant drivers of genome evolution for the two orang-utan populations.

Uncovering the cultural knowledge of sanctuary apes

Behavioral differences observed between wild communities of the same species have been called “cultures” by some researchers who aimed to underline the similarities with human cultures. However, whether these differences truly result from social learning processes is debated. Despite promising recent research, data acquired in the wild still fail to exclude genetic and ecological factors from being potential explanations for the observed behavioral differences. A potential way to address this problem is through field experiments where communities of the same subspecies are exposed to identical apparatuses. This way, genetic and ecological factors can be controlled for, although their influence cannot be fully excluded. Working with wild-born Sumatran orangutans originating from two genetically distinct populations, we recently combined field experiments with captive work to show that genetic differences could not account for differences in their knowledge of stick use. Additionally, we found evidence that our subjects arrived at the sanctuary with a knowledge that they acquired but could not express in their community of origin. These findings suggest that animal cultures must also be analyzed at the cognitive level. Only in this way can we understand the true extent of animal cultures and how they relate to human cultures.