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

Evidence for the smallest fossil Pongo in southern China

The rarity of Pongo fossils with precise absolute dating from the Middle Pleistocene hampers our understanding of the taxonomy and spatiotemporal distribution of Quaternary orangutans in southern China. Here, we report a newly discovered sample of 113 isolated teeth of fossil Pongo from Zhongshan Cave in the Bubing Basin, Guangxi, southern China. We describe the Pongo specimens from Zhongshan Cave and compare them metrically to other samples of fossil Pongo species (i.e., Pongo weidenreichi, Pongo devosi, Pongo duboisi, Pongo palaeosumatrensis, Pongo javensis, and Pongo sp.) and to extant orangutans (i.e., Pongo pygmaeus and Pongo abelii). The Zhongshan Pongo assemblage is dated using U-series and coupled electron spin resonance/U-series methods. Our results reasonably constrain the Zhongshan Pongo assemblage to 184 ± 16 ka, which is consistent with the biostratigraphic evidence. The Zhongshan Pongo teeth are only 6.5% larger on average than those of extant Pongo. The Zhongshan teeth are smaller overall than those of Pongo from all other cave sites in southern China, and they currently represent the smallest fossil orangutans in southern China. Based on their dental size, and the presence of a well-developed lingual pillar and lingual cingulum on the upper and lower incisors, an intermediate frequency of lingual cingulum remnants on the upper molars, and a higher frequency of moderate to heavy wrinkling on the upper and lower molars, we provisionally assign the Zhongshan fossils to P. devosi. Our results confirm earlier claims that P. weidenreichi is replaced by a smaller species in southern China, P. devosi, by the late Middle Pleistocene. The occurrence of P. devosi in Zhongshan Cave further extends its spatial and temporal distribution. The Pongo specimens from Zhongshan provide important new evidence to demonstrate that the dental morphological features of Pongo in southern China changed substantially during the late Middle Pleistocene.

Phylogenomic analyses reveal incongruences between divergence times and fossil records of freshwater snails in East Asia

Fossils provide important insight into our understanding of phylogenetic history by serving as calibration points for divergence time estimation. However, uncertainties in the fossil record due to parallel evolution and convergent evolution can critically affect estimates of node ages. Here, we compare and contrast estimates of phylogenetic divergence with geologic and fossil history for two freshwater snail genera of the family Viviparidae in East Asia (Cipangopaludina and Margarya). Cipangopaludina species are commonly widely distributed species in East Asia, but extant Margarya species are endemic to the ancient lakes in Yunnan, China. According to some previous studies, parallel evolution or convergent evolution of shell morphology has occurred in the family several times which may affect divergence time estimation using fossil records. In this study, we used SNP data derived from ddRAD-seq loci to investigate population demographic history of both genera. Our results show a common pattern of lake endemic lineages diversifying from widely distributed lineages in the Miocene, and multiple colonization to a single ancient lake occurred in the Pleistocene. Our results indicate substantial incongruence among estimated phylogenomic divergence times, some fossil records, and formation ages of ancient lakes. These findings suggest some fossil records may be misidentified in these groups and highlight the need to carefully evaluate geological evidence and fossil records when using these for divergence time estimation.

How ancient forest fragmentation and riparian connectivity generate high levels of genetic diversity in a microendemic Malagasy tree

Understanding landscape changes is central to predicting evolutionary trajectories and defining conservation practices. While human-driven deforestation is intense throughout Madagascar, exceptions in areas such as the Loky-Manambato region (north) raise questions regarding the causes and age of forest fragmentation. The Loky-Manambato region also harbours a rich and endemic flora, whose evolutionary origin remains poorly understood. We assessed the genetic diversity of an endangered microendemic Malagasy olive species (Noronhia spinifolia Hong-Wa) to better understand the vegetation dynamics in the Loky-Manambato region and its influence on past evolutionary processes. We characterized 72 individuals sampled across eight forests through nuclear and mitochondrial restriction-associated DNA sequencing data and chloroplast microsatellites. Combined population and landscape genetics analyses indicate that N. spinifolia diversity is largely explained by the current forest cover, highlighting a long-standing habitat mosaic in the region. This sustains a major and long-term role of riparian corridors in maintaining connectivity across these antique mosaic habitats, calling for the study of organismal interactions that promote gene flow.

The anthropogenic effect of land use on population genetics of Malcus inconspicuus

Since the beginning of the Holocene era, human activities have seriously impacted animal habitats and vegetative environments. Species that are dependent on natural habitats or with narrow niches might be more severely affected by habitat changes. Malcus inconspicuus is distributed in subtropical China and highly dependent on the mountain environment. Our study investigated the role of the mountainous landscape in the historical evolution of M. inconspicuus and the impact of Holocene human activities on it. A phylogeographical approach was implemented with integrative datasets including double-digest restriction site-associated DNA (ddRAD), mitochondrial data, and distribution data. Three obvious clades and an east-west phylogeographical pattern were found in subtropical China. Mountainous landscape has "multifaceted" effects on the evolutionary history of M. inconspicuus, it has contributed to population differentiation, provided glacial refuges, and provided population expansion corridors during the postglacial period. The effective population size (Ne) of M. inconspicuus showed a sharp decline during the Holocene era, which revealed a significantly negative correlation with the development of cropland in a hilly area at the same time and space. It supported that the species which are highly dependent on natural habitats might undergo greater impact when the habitat was damaged by agricultural activities and we should pay more attention to them, especially in the land development of their distribution areas.

Evidence for the latest fossil Pongo in southern China

Pongo fossils with precise absolute age brackets are rare, limiting our understanding of their taxonomy and spatiotemporal distribution in southern China during the Late Pleistocene. Twenty-four isolated teeth of fossil orangutans were recently discovered during excavations at Yicun Cave in Guangxi Zhuang Autonomous Region, southern China. Here, we dated the fossil-bearing layer using Uranium-series dating of the associated flowstone and soda straw stalactites. Our results date the Yicun orangutan fossils to between 66 ± 0.32 ka and 57 ± 0.26 ka; thus, these fossils currently represent the last appearance datum of Pongo in southern China. We further conducted a detailed morphological comparison of the Yicun fossil teeth with large samples of fossil (n = 2454) and extant (n = 441) orangutans from mainland and island Southeast Asia to determine their taxonomic position. Compared to other fossil and extant orangutan samples, the Yicun Pongo assemblage has larger teeth and displays greater variation in occlusal structure. Based on the high frequency of cingular remnants and light to moderate enamel wrinkling of the molars, we assigned the Yicun fossils to Pongo weidenreichi, a species that was widespread in southern China throughout the Pleistocene. Lastly, we used published stable carbon isotope data of Early to Late Pleistocene mammalian fossil teeth from mainland Southeast Asia to reconstruct changes in the paleoenvironment and to interpret dental size variation of Pongo assemblages in a broader temporal and environmental context. The carbon isotope data show that dental size reduction in Pongo is associated with environmental changes. These morphological changes in Pongo appear to coincide with the expansion of savannah biomes and the contraction of forest habitats from the Middle Pleistocene onward. The variation in dental size of forest-dwelling Pongo in mainland Southeast Asia may have resulted from habitat differentiation during the Pleistocene.

Global genetic diversity status and trends: towards a suite of Essential Biodiversity Variables (EBVs) for genetic composition

Biodiversity underlies ecosystem resilience, ecosystem function, sustainable economies, and human well‐being. Understanding how biodiversity sustains ecosystems under anthropogenic stressors and global environmental change will require new ways of deriving and applying biodiversity data. A major challenge is that biodiversity data and knowledge are scattered, biased, collected with numerous methods, and stored in inconsistent ways. The Group on Earth Observations Biodiversity Observation Network (GEO BON) has developed the Essential Biodiversity Variables (EBVs) as fundamental metrics to help aggregate, harmonize, and interpret biodiversity observation data from diverse sources. Mapping and analyzing EBVs can help to evaluate how aspects of biodiversity are distributed geographically and how they change over time. EBVs are also intended to serve as inputs and validation to forecast the status and trends of biodiversity, and to support policy and decision making. Here, we assess the feasibility of implementing Genetic Composition EBVs (Genetic EBVs), which are metrics of within‐species genetic variation. We review and bring together numerous areas of the field of genetics and evaluate how each contributes to global and regional genetic biodiversity monitoring with respect to theory, sampling logistics, metadata, archiving, data aggregation, modeling, and technological advances. We propose four Genetic EBVs: (i) Genetic Diversity; (ii) Genetic Differentiation; (iii) Inbreeding; and (iv) Effective Population Size (N_e). We rank Genetic EBVs according to their relevance, sensitivity to change, generalizability, scalability, feasibility and data availability. We outline the workflow for generating genetic data underlying the Genetic EBVs, and review advances and needs in archiving genetic composition data and metadata. We discuss how Genetic EBVs can be operationalized by visualizing EBVs in space and time across species and by forecasting Genetic EBVs beyond current observations using various modeling approaches. Our review then explores challenges of aggregation, standardization, and costs of operationalizing the Genetic EBVs, as well as future directions and opportunities to maximize their uptake globally in research and policy. The collection, annotation, and availability of genetic data has made major advances in the past decade, each of which contributes to the practical and standardized framework for large‐scale genetic observation reporting. Rapid advances in DNA sequencing technology present new opportunities, but also challenges for operationalizing Genetic EBVs for biodiversity monitoring regionally and globally. With these advances, genetic composition monitoring is starting to be integrated into global conservation policy, which can help support the foundation of all biodiversity and species' long‐term persistence in the face of environmental change. We conclude with a summary of concrete steps for researchers and policy makers for advancing operationalization of Genetic EBVs. The technical and analytical foundations of Genetic EBVs are well developed, and conservation practitioners should anticipate their increasing application as efforts emerge to scale up genetic biodiversity monitoring regionally and globally.

Impact of model assumptions on demographic inferences: the case study of two sympatric mouse lemurs in northwestern Madagascar

BACKGROUND

Quaternary climate fluctuations have been acknowledged as major drivers of the geographical distribution of the extraordinary biodiversity observed in tropical biomes, including Madagascar. The main existing framework for Pleistocene Malagasy diversification assumes that forest cover was strongly shaped by warmer Interglacials (leading to forest expansion) and by cooler and arid glacials (leading to forest contraction), but predictions derived from this scenario for forest-dwelling animals have rarely been tested with genomic datasets.

RESULTS

We generated genomic data and applied three complementary demographic approaches (Stairway Plot, PSMC and IICR-simulations) to infer population size and connectivity changes for two forest-dependent primate species (Microcebus murinus and M. ravelobensis) in northwestern Madagascar. The analyses suggested major demographic changes in both species that could be interpreted in two ways, depending on underlying model assumptions (i.e., panmixia or population structure). Under panmixia, the two species exhibited larger population sizes across the Last Glacial Maximum (LGM) and towards the African Humid Period (AHP). This peak was followed by a population decline in M. ravelobensis until the present, while M. murinus may have experienced a second population expansion that was followed by a sharp decline starting 3000 years ago. In contrast, simulations under population structure suggested decreasing population connectivity between the Last Interglacial and the LGM for both species, but increased connectivity during the AHP exclusively for M. murinus.

CONCLUSION

Our study shows that closely related species may differ in their responses to climatic events. Assuming that Pleistocene climatic conditions in the lowlands were similar to those in the Malagasy highlands, some demographic dynamics would be better explained by changes in population connectivity than in population size. However, changes in connectivity alone cannot be easily reconciled with a founder effect that was shown for M. murinus during its colonization of the northwestern Madagascar in the late Pleistocene. To decide between the two alternative models, more knowledge about historic forest dynamics in lowland habitats is necessary. Altogether, our study stresses that demographic inferences strongly depend on the underlying model assumptions. Final conclusions should therefore be based on a comparative evaluation of multiple approaches.

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.

Evidence for the Pleistocene Arc Hypothesis from genome-wide SNPs in a Neotropical dry forest specialist, the Rufous-fronted Thornbird (Furnariidae: Phacellodomus rufifrons)

South American dry forests have a complex and poorly understood biogeographic history. Based on the fragmented distribution of many Neotropical dry forest species, it has been suggested that this biome was more widely distributed and contiguous under drier climate conditions in the Pleistocene. To test this scenario, known as the Pleistocene Arc Hypothesis, we studied the phylogeography of the Rufous-fronted Thornbird (Phacellodomus rufifrons), a widespread dry forest bird with a disjunct distribution closely matching that of the biome itself. We sequenced mtDNA and used ddRADseq to sample 7,167 genome-wide single-nucleotide polymorphisms from 74 P. rufifrons individuals across its range. We found low genetic differentiation over two prominent geographic breaks - particularly across a 1,000 km gap between populations in Bolivia and Northern Peru. Using demographic analyses of the joint site frequency spectrum, we found evidence of recent divergence without subsequent gene flow across those breaks. By contrast, parapatric morphologically distinct populations in northeastern Brazil show high genetic divergence with evidence of recent gene flow. These results, in combination with our paleoclimate species distribution modelling, support the idea that currently disjunct patches of dry forest were more connected in the recent past, probably during the Middle and Late Pleistocene. This notion fits the major predictions of the Pleistocene Arc Hypothesis and illustrates the importance of comprehensive genomic and geographic sampling for examining biogeographic and evolutionary questions in complex ecosystems like Neotropical dry forests.

Population genetics of wild Macaca fascicularis with low-coverage shotgun sequencing of museum specimens

OBJECTIVES

Long-tailed macaques (Macaca fascicularis) are widely distributed throughout the mainland and islands of Southeast Asia, making them a useful model for understanding the complex biogeographical history resulting from drastic changes in sea levels throughout the Pleistocene. Past studies based on mitochondrial genomes (mitogenomes) of long-tailed macaque museum specimens have traced their colonization patterns throughout the archipelago, but mitogenomes trace only the maternal history. Here, our objectives were to trace phylogeographic patterns of long-tailed macaques using low-coverage nuclear DNA (nDNA) data from museum specimens.

METHODS

We performed population genetic analyses and phylogenetic reconstruction on nuclear single nucleotide polymorphisms (SNPs) from shotgun sequencing of 75 long-tailed macaque museum specimens from localities throughout Southeast Asia.

RESULTS

We show that shotgun sequencing of museum specimens yields sufficient genome coverage (average ~1.7%) for reconstructing population relationships using SNP data. Contrary to expectations of divergent results between nuclear and mitochondrial genomes for a female philopatric species, phylogeographical patterns based on nuclear SNPs proved to be closely similar to those found using mitogenomes. In particular, population genetic analyses and phylogenetic reconstruction from the nDNA identify two major clades within M. fascicularis: Clade A includes all individuals from the mainland along with individuals from northern Sumatra, while Clade B consists of the remaining island-living individuals, including those from southern Sumatra.

CONCLUSIONS

Overall, we demonstrate that low-coverage sequencing of nDNA from museum specimens provides enough data for examining broad phylogeographic patterns, although greater genome coverage and sequencing depth would be needed to distinguish between very closely related populations, such as those throughout the Philippines.

Genetic analyses reveal population structure and recent decline in leopards (Panthera pardus fusca) across the Indian subcontinent

Background

Large carnivores maintain the stability and functioning of ecosystems. Currently, many carnivore species face declining population sizes due to natural and anthropogenic pressures. The leopard, Panthera pardus, is probably the most widely distributed and highly adaptable large felid globally, still persisting in most of its historic range. However, we lack subspecies-level data on country or regional scale on population trends, as ecological monitoring approaches are difficult to apply on such wide-ranging species. We used genetic data from leopards sampled across the Indian subcontinent to investigate population structure and patterns of demographic decline.

Methods

We collected faecal samples from the Terai-Arc landscape of northern India and identified 56 unique individuals using a panel of 13 microsatellite markers. We merged this data with already available 143 leopard individuals and assessed genetic structure at country scale. Subsequently, we investigated the demographic history of each identified subpopulations and compared genetic decline analyses with countrywide local extinction probabilities.

Results

Our genetic analyses revealed four distinct subpopulations corresponding to Western Ghats, Deccan Plateau-Semi Arid, Shivalik and Terai region of the north Indian landscape, each with high genetic variation. Coalescent simulations with microsatellite loci revealed a possibly human-induced 75–90% population decline between ∼120–200 years ago across India. Population-specific estimates of genetic decline are in concordance with ecological estimates of local extinction probabilities in these subpopulations obtained from occupancy modeling of the historic and current distribution of leopards in India.

Conclusions

Our results confirm the population decline of a widely distributed, adaptable large carnivore. We re-iterate the relevance of indirect genetic methods for such species in conjunction with occupancy assessment and recommend that detailed, landscape-level ecological studies on leopard populations are critical to future conservation efforts. Our approaches and inference are relevant to other widely distributed, seemingly unaffected carnivores such as the leopard.

ABLE: blockwise site frequency spectra for inferring complex population histories and recombination

We introduce ABLE (Approximate Blockwise Likelihood Estimation), a novel simulation-based composite likelihood method that uses the blockwise site frequency spectrum to jointly infer past demography and recombination. ABLE is explicitly designed for a wide variety of data from unphased diploid genomes to genome-wide multi-locus data (for example, RADSeq) and can also accommodate arbitrarily large samples. We use simulations to demonstrate the accuracy of this method to infer complex histories of divergence and gene flow and reanalyze whole genome data from two species of orangutan. ABLE is available for download at https://github.com/champost/ABLE.

Orangutans venture out of the rainforest and into the Anthropocene

Conservation benefits from understanding how adaptability and threat interact to determine a taxon's vulnerability. Recognizing how interactions with humans have shaped taxa such as the critically endangered orangutan (Pongo spp.) offers insights into this relationship. Orangutans are viewed as icons of wild nature, and most efforts to prevent their extinction have focused on protecting minimally disturbed habitat, with limited success. We synthesize fossil, archeological, genetic, and behavioral evidence to demonstrate that at least 70,000 years of human influence have shaped orangutan distribution, abundance, and ecology and will likely continue to do so in the future. Our findings indicate that orangutans are vulnerable to hunting but appear flexible in response to some other human activities. This highlights the need for a multifaceted, landscape-level approach to orangutan conservation that leverages sound policy and cooperation among government, private sector, and community stakeholders to prevent hunting, mitigate human-orangutan conflict, and preserve and reconnect remaining natural forests. Broad cooperation can be encouraged through incentives and strategies that focus on the common interests and concerns of different stakeholders. Orangutans provide an illustrative example of how acknowledging the long and pervasive influence of humans can improve strategies to preserve biodiversity in the Anthropocene.

FOXP2 variation in great ape populations offers insight into the evolution of communication skills

The gene coding for the forkhead box protein P2 (FOXP2) is associated with human language disorders. Evolutionary changes in this gene are hypothesized to have contributed to the emergence of speech and language in the human lineage. Although FOXP2 is highly conserved across most mammals, humans differ at two functional amino acid substitutions from chimpanzees, bonobos and gorillas, with an additional fixed substitution found in orangutans. However, FOXP2 has been characterized in only a small number of apes and no publication to date has examined the degree of natural variation in large samples of unrelated great apes. Here, we analyzed the genetic variation in the FOXP2 coding sequence in 63 chimpanzees, 11 bonobos, 48 gorillas, 37 orangutans and 2 gibbons and observed undescribed variation in great apes. We identified two variable polyglutamine microsatellites in chimpanzees and orangutans and found three nonsynonymous single nucleotide polymorphisms, one in chimpanzees, one in gorillas and one in orangutans with derived allele frequencies of 0.01, 0.26 and 0.29, respectively. Structural and functional protein modeling indicate a biochemical effect of the substitution in orangutans, and because of its presence solely in the Sumatran orangutan species, the mutation may be associated with reported population differences in vocalizations.

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.

Tracing the phylogeographic history of Southeast Asian long-tailed macaques through mitogenomes of museum specimens

The biogeographical history of Southeast Asia is complicated due to the continuous emergences and disappearances of land bridges throughout the Pleistocene. Here, we use long-tailed macaques (Macaca fascicularis), which are widely distributed throughout the mainland and islands of Southeast Asia, asa model for better understanding the biogeographical patterns of diversification in this geographically complex region. A reliable intraspecific phylogeny including individuals from localities on oceanic islands, continental islands, and the mainland is needed to trace relatedness along with the pattern and timing of colonization in this region. We used high-throughput sequencing techniques to sequence mitochondrial genomes (mitogenomes) from 95 Southeast Asian M. fascicularis specimens housed at natural history museums around the world. To achieve a comprehensive picture, we more than tripled the mitogenome sample size for M. fascicularis from previous studies, and for the first time included documented samples from the Philippines and several small Indonesian islands. Confirming the result from a previous, recent intraspecific phylogeny for M. fascicularis, the newly reconstructed phylogeny of 135 specimens divides the samples into two major clades: Clade A includes haplotypes from the mainland and some from northern Sumatra, while Clade B includes all insular haplotypes along with lineages from southern Sumatra. This study resolves a previous disparity by revealing a disjunction in the origin of Sumatran macaques, with separate lineages originating within the two major clades, suggesting that at least two major migrations to Sumatra occurred. However, our dated phylogeny reveals that the two major clades split ∼1.88Ma, which is earlier than in previously published phylogenies. Our new data reveal that most Philippine macaque lineages diverged from the Borneo stock within the last ∼0.06-0.43Ma. Finally, our study provides insight into successful sequencing of DNA across museums and shotgun sequencing of DNA specimens asa method to sequence the mitogenome.

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.

Evidence that the rate of strong selective sweeps increases with population size in the great apes

Quantifying the number of selective sweeps and their combined effects on genomic diversity in humans and other great apes is notoriously difficult. Here we address the question using a comparative approach to contrast diversity patterns according to the distance from genes in all great ape taxa. The extent of diversity reduction near genes compared with the rest of intergenic sequences is greater in a species with larger effective population size. Also, the maximum distance from genes at which the diversity reduction is observed is larger in species with large effective population size. In Sumatran orangutans, the overall genomic diversity is ∼30% smaller than diversity levels far from genes, whereas this reduction is only 9% in humans. We show by simulation that selection against deleterious mutations in the form of background selection is not expected to cause these differences in diversity among species. Instead, selective sweeps caused by positive selection can reduce diversity level more severely in a large population if there is a higher number of selective sweeps per unit time. We discuss what can cause such a correlation, including the possibility that more frequent sweeps in larger populations are due to a shorter waiting time for the right mutations to arise.

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.

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short-range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19-29 individuals from each population. The sample includes two pairs of low- and high-altitude populations from two different valleys. High-altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low-altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low- and high-altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low- and high-altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low- and high-altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.

Reconstructing the demographic history of divergence between European river and brook lampreys using approximate Bayesian computations

Inferring the history of isolation and gene flow during species divergence is a central question in evolutionary biology. The European river lamprey (Lampetra fluviatilis) and brook lamprey (L. planeri) show a low reproductive isolation but have highly distinct life histories, the former being parasitic-anadromous and the latter non-parasitic and freshwater resident. Here we used microsatellite data from six replicated population pairs to reconstruct their history of divergence using an approximate Bayesian computation framework combined with a random forest model. In most population pairs, scenarios of divergence with recent isolation were outcompeted by scenarios proposing ongoing gene flow, namely the Secondary Contact (SC) and Isolation with Migration (IM) models. The estimation of demographic parameters under the SC model indicated a time of secondary contact close to the time of speciation, explaining why SC and IM models could not be discriminated. In case of an ancient secondary contact, the historical signal of divergence is lost and neutral markers converge to the same equilibrium as under the less parameterized model allowing ongoing gene flow. Our results imply that models of secondary contacts should be systematically compared to models of divergence with gene flow; given the difficulty to discriminate among these models, we suggest that genome-wide data are needed to adequately reconstruct divergence history.

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.

Derivation of induced pluripotent stem cells from orangutan skin fibroblasts

Background

Orangutans are an endangered species whose natural habitats are restricted to the Southeast Asian islands of Borneo and Sumatra. Along with the African great apes, orangutans are among the closest living relatives to humans. For potential species conservation and functional genomics studies, we derived induced pluripotent stem cells (iPSCs) from cryopreserved somatic cells obtained from captive orangutans.

Results

Primary skin fibroblasts from two Sumatran orangutans were transduced with retroviral vectors expressing the human OCT4, SOX2, KLF4, and c-MYC factors. Candidate orangutan iPSCs were characterized by global gene expression and DNA copy number analysis. All were consistent with pluripotency and provided no evidence of large genomic insertions or deletions. In addition, orangutan iPSCs were capable of producing cells derived from all three germ layers in vitro through embryoid body differentiation assays and in vivo through teratoma formation in immune-compromised mice.

Conclusions

We demonstrate that orangutan skin fibroblasts are capable of being reprogrammed into iPSCs with hallmark molecular signatures and differentiation potential. We suggest that reprogramming orangutan somatic cells in genome resource banks could provide new opportunities for advancing assisted reproductive technologies relevant for species conservation efforts. Furthermore, orangutan iPSCs could have applications for investigating the phenotypic relevance of genomic changes that occurred in the human, African great ape, and/or orangutan lineages. This provides opportunities for orangutan cell culture models that would otherwise be impossible to develop from living donors due to the invasive nature of the procedures required for obtaining primary cells.

Electronic supplementary material

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