The genomic diversity of the Eliurus genus in northern Madagascar with a putative new species.

Madagascar exhibits extraordinarily high level of species richness and endemism, while being severely threatened by habitat loss and fragmentation (HL&F). In front of these threats to biodiversity, conservation effort can be directed, for instance, in the documentation of species that are still unknown to science, or in investigating how species respond to HL&F. The tufted-tail rats genus (Eliurus spp.) is the most speciose genus of endemic rodents in Madagascar, with 13 described species, which occupy two major habitat types: dry or humid forests. The large species diversity and association to specific habitat types make the Eliurus genus a suitable model for investigating species adaptation to new environments, as well as response to HL&F (dry vs humid). In the present study, we investigated Eliurus spp. genomic diversity across northern Madagascar, a region covered by both dry and humid fragmented forests. From the mitochondrial DNA (mtDNA) and nuclear genomic (RAD-seq) data of 124 Eliurus individuals sampled in poorly studied forests of northern Madagascar, we identified an undescribed Eliurus taxon (Eliurus sp. nova). We tested the hypothesis of a new Eliurus species using several approaches: i) DNA barcoding; ii) phylogenetic inferences; iii) species delimitation tests based on the Multi-Species Coalescent (MSC) model, iv) genealogical divergence index (gdi); v) an ad-hoc test of isolation-by-distance within versus between sister-taxa, vi) comparisons of %GC content patterns and vii) morphological analyses. All analyses support the recognition of the undescribed lineage as a putative distinct species. In addition, we show that Eliurus myoxinus, a species known from the dry forests of western Madagascar, is, surprisingly, found mostly in humid forests in northern Madagascar. In conclusion, we discuss the implications of such findings in the context of Eliurus species evolution and diversification, and use the distribution of northern Eliurus species as a proxy for reconstructing past changes in forest cover and vegetation type in northern Madagascar.

The importance of well protected forests for the conservation genetics of West African colobine monkeys.

In tropical forests, anthropogenic activities are major drivers of the destruction and degradation of natural habitats, causing severe biodiversity loss. African colobine monkeys (Colobinae) are mainly folivore and strictly arboreal primates that require large forests to subsist, being among the most vulnerable of all nonhuman primates. The Western red colobus Piliocolobus badius and the King colobus Colobus polykomos inhabit highly fragmented West African forests, including the Cantanhez Forests National Park (CFNP) in Guinea-Bissau. Both species are also found in the largest and best-preserved West African forest-the Taï National Park (TNP) in Ivory Coast. Colobine monkeys are hunted for bushmeat in both protected areas, but these exhibit contrasting levels of forest fragmentation, thus offering an excellent opportunity to investigate the importance of well-preserved forests for the maintenance of evolutionary potential in these arboreal primates. We estimated genetic diversity, population structure, and demographic history by using microsatellite loci and mitochondrial DNA. We then compared the genetic patterns of the colobines from TNP with the ones previously obtained for CFNP and found contrasting genetic patterns. Contrary to the colobines from CFNP that showed very low genetic diversity and a strong population decline, the populations in TNP still maintain high levels of genetic diversity and we found no clear signal of population decrease in Western red colobus and a limited decrease in King colobus. These results suggest larger and historically more stable populations in TNP compared to CFNP. We cannot exclude the possibility that the demographic effects resulting from the recent increase of bushmeat hunting are not yet detectable in TNP using genetic data. Nevertheless, the fact that the TNP colobus populations are highly genetically diverse and maintain large effective population sizes suggests that well-preserved forests are crucial for the maintenance of populations, species, and probably for the evolutionary potential in colobines.

Population genomic structure in Goodman's mouse lemur reveals long-standing separation of Madagascar's Central Highlands and eastern rainforests.

Madagascar's Central Highlands are largely composed of grasslands, interspersed with patches of forest. The historical perspective was that Madagascar's grasslands had anthropogenic origins, but emerging evidence suggests that grasslands were a component of the pre-human Central Highlands vegetation. Consequently, there is now vigorous debate regarding the extent to which these grasslands have expanded due to anthropogenic pressures. Here, we shed light on the temporal dynamics of Madagascar's vegetative composition by conducting a population genomic investigation of Goodman's mouse lemur (Microcebus lehilahytsara; Cheirogaleidae). These small-bodied primates occur both in Madagascar's eastern rainforests and in the Central Highlands, making them a valuable indicator species. Population divergences among forest-dwelling mammals will reflect changes to their habitat, including fragmentation, whereas patterns of post-divergence gene flow can reveal formerly wooded migration corridors. To explore these patterns, we used RADseq data to infer population genetic structure, demographic models of post-divergence gene flow, and population size change through time. The results offer evidence that open habitats are an ancient component of the Central Highlands, and that widespread forest fragmentation occurred naturally during a period of decreased precipitation near the last glacial maximum. Models of gene flow suggest that migration across the Central Highlands has been possible from the Pleistocene through the recent Holocene via riparian corridors. Though our findings support the hypothesis that Central Highland grasslands predate human arrival, we also find evidence for human-mediated population declines. This highlights the extent to which species imminently threatened by human-mediated deforestation may already be vulnerable from paleoclimatic conditions.

Cryptic Patterns of Speciation in Cryptic Primates: Microendemic Mouse Lemurs and the Multispecies Coalescent.

Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic revision has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we investigate mouse lemur diversity in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for two pairs of sister lineages that include three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among the two pairs of lineages. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species, M. mittermeieri and M. lehilahytsara is not supported-a result that is particularly striking when using the genealogical discordance index (gdi). Nonsister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated relaxed-clock estimates. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case. [Cryptic species; effective population size; microendemism; multispecies coalescent; speciation; species delimitation.].

Inferring number of populations and changes in connectivity under the n-island model.

Inferring the demographic history of species is one of the greatest challenges in populations genetics. This history is often represented as a history of size changes, ignoring population structure. Alternatively, when structure is assumed, it is defined a priori as a population tree and not inferred. Here we propose a framework based on the IICR (Inverse Instantaneous Coalescence Rate). The IICR can be estimated for a single diploid individual using the PSMC method of Li and Durbin (2011). For an isolated panmictic population, the IICR matches the population size history, and this is how the PSMC outputs are generally interpreted. However, it is increasingly acknowledged that the IICR is a function of the demographic model and sampling scheme with limited connection to population size changes. Our method fits observed IICR curves of diploid individuals with IICR curves obtained under piecewise stationary symmetrical island models. In our models we assume a fixed number of time periods during which gene flow is constant, but gene flow is allowed to change between time periods. We infer the number of islands, their sizes, the periods at which connectivity changes and the corresponding rates of connectivity. Validation with simulated data showed that the method can accurately recover most of the scenario parameters. Our application to a set of five human PSMCs yielded demographic histories that are in agreement with previous studies using similar methods and with recent research suggesting ancient human structure. They are in contrast with the view of human evolution consisting of one ancestral population branching into three large continental and panmictic populations with varying degrees of connectivity and no population structure within each continent.

The Genetic Legacy of Zoroastrianism in Iran and India: Insights into Population Structure, Gene Flow, and Selection.

Zoroastrianism is one of the oldest extant religions in the world, originating in Persia (present-day Iran) during the second millennium BCE. Historical records indicate that migrants from Persia brought Zoroastrianism to India, but there is debate over the timing of these migrations. Here we present genome-wide autosomal, Y chromosome, and mitochondrial DNA data from Iranian and Indian Zoroastrians and neighboring modern-day Indian and Iranian populations and conduct a comprehensive genome-wide genetic analysis in these groups. Using powerful haplotype-based techniques, we find that Zoroastrians in Iran and India have increased genetic homogeneity relative to other sampled groups in their respective countries, consistent with their current practices of endogamy. Despite this, we infer that Indian Zoroastrians (Parsis) intermixed with local groups sometime after their arrival in India, dating this mixture to 690-1390 CE and providing strong evidence that Iranian Zoroastrian ancestry was maintained primarily through the male line. By making use of the rich information in DNA from ancient human remains, we also highlight admixture in the ancestors of Iranian Zoroastrians dated to 570 BCE-746 CE, older than admixture seen in any other sampled Iranian group, consistent with a long-standing isolation of Zoroastrians from outside groups. Finally, we report results, and challenges, from a genome-wide scan to identify genomic regions showing signatures of positive selection in present-day Zoroastrians that might correlate to the prevalence of particular diseases among these communities.

Genetic consequences of social structure in the golden-crowned sifaka.

Many species are structured in social groups (SGs) where individuals exhibit complex mating strategies. Yet, most population genetic studies ignore SGs either treating them as small random-mating units or focusing on a higher hierarchical level (the population). Empirical studies acknowledging SGs have found an overall excess of heterozygotes within SGs and usually invoke inbreeding avoidance strategies to explain this finding. However, there is a lack of null models against which ecological theories can be tested and inbreeding avoidance quantified. Here, we investigate inbreeding (deviation from random mating) in an endangered forest-dwelling pair-living lemur species (Propithecus tattersalli). In particular, we measure the inbreeding coefficient (FIS) in empirical data at different scales: SGs, sampling sites and forest patches. We observe high excess of heterozygotes within SGs. The magnitude of this excess is highly dependent on the sampling scheme: while offspring are characterised by a high excess of heterozygotes (FIS < 0), the reproductive pair does not show dramatic departures from Hardy-Weinberg expectations. Moreover, the heterozygosity excess disappears at larger geographic scales (sites and forests). We use a modelling framework that incorporates details of the sifaka mating system but does not include active inbreeding avoidance mechanisms. The simulated data show that, although apparent "random mating" or even inbreeding may occur at the "population" level, outbreeding is maintained within SGs. Altogether our results suggest that social structure leads to high levels of outbreeding without the need for active inbreeding avoidance mechanisms. Thus, demonstrating and measuring the existence of active inbreeding avoidance mechanisms may be more difficult than usually assumed.

Ecology and morphology of mouse lemurs (Microcebus spp.) in a hotspot of microendemism in northeastern Madagascar, with the description of a new species.

Delimitation of cryptic species is increasingly based on genetic analyses but the integration of distributional, morphological, behavioral, and ecological data offers unique complementary insights into species diversification. We surveyed communities of nocturnal mouse lemurs (Microcebus spp.) in five different sites of northeastern Madagascar, measuring a variety of morphological parameters and assessing reproductive states for 123 individuals belonging to five different lineages. We documented two different non-sister lineages occurring in sympatry in two areas. In both cases, sympatric species pairs consisted of a locally restricted (M. macarthurii or M. sp. #3) and a more widespread lineage (M. mittermeieri or M. lehilahytsara). Estimated Extents of Occurrence (EOO) of these lineages differed remarkably with 560 and 1,500 km2 versus 9,250 and 50,700 km2 , respectively. Morphometric analyses distinguished unambiguously between sympatric species and detected more subtle but significant differences among sister lineages. Tail length and body size were most informative in this regard. Reproductive schedules were highly variable among lineages, most likely impacted by phylogenetic relatedness and environmental variables. While sympatric species pairs differed in their reproductive timing (M. sp. #3/M. lehilahytsara and M. macarthurii/M. mittermeieri), warmer lowland rainforests were associated with a less seasonal reproductive schedule for M. mittermeieri and M. lehilahytsara compared with populations occurring in montane forests. Distributional, morphological, and ecological data gathered in this study support the results of genomic species delimitation analyses conducted in a companion study, which identified one lineage, M. sp. #3, as meriting formal description as a new species. Consequently, a formal species description is included. Worryingly, our data also show that geographically restricted populations of M. sp. #3 and its sister species (M. macarthurii) are at high risk of local and perhaps permanent extinction from both deforestation and habitat fragmentation.

Consequences of breed formation on patterns of genomic diversity and differentiation: the case of highly diverse peripheral Iberian cattle.

BACKGROUND: Iberian primitive breeds exhibit a remarkable phenotypic diversity over a very limited geographical space. While genomic data are accumulating for most commercial cattle, it is still lacking for these primitive breeds. Whole genome data is key to understand the consequences of historic breed formation and the putative role of earlier admixture events in the observed diversity patterns. RESULTS: We sequenced 48 genomes belonging to eight Iberian native breeds and found that the individual breeds are genetically very distinct with FST values ranging from 4 to 16% and have levels of nucleotide diversity similar or larger than those of their European counterparts, namely Jersey and Holstein. All eight breeds display significant gene flow or admixture from African taurine cattle and include mtDNA and Y-chromosome haplotypes from multiple origins. Furthermore, we detected a very low differentiation of chromosome X relative to autosomes within all analyzed taurine breeds, potentially reflecting male-biased gene flow. CONCLUSIONS: Our results show that an overall complex history of admixture resulted in unexpectedly high levels of genomic diversity for breeds with seemingly limited geographic ranges that are distantly located from the main domestication center for taurine cattle in the Near East. This is likely to result from a combination of trading traditions and breeding practices in Mediterranean countries. We also found that the levels of differentiation of autosomes vs sex chromosomes across all studied taurine and indicine breeds are likely to have been affected by widespread breeding practices associated with male-biased gene flow.

Coalescence times for three genes provide sufficient information to distinguish population structure from population size changes.

The increasing amount of genomic data currently available is expanding the horizons of population genetics inference. A wide range of methods have been published allowing to detect and date major changes in population size during the history of species. At the same time, there has been an increasing recognition that population structure can generate genetic data similar to those generated under models of population size change. Recently, Mazet et al. (Heredity 116(4):362-371, 2016) introduced the idea that, for any model of population structure, it is always possible to find a panmictic model with a particular function of population size-change having an identical distribution of [Formula: see text] (the time of the first coalescence for a sample of size two). This implies that there is an identifiability problem between a panmictic and a structured model when we base our analysis only on [Formula: see text]. In this paper, based on an analytical study of the rate matrix of the ancestral lineage process, we obtain new theoretical results about the joint distribution of the coalescence times [Formula: see text] for a sample of three haploid genes in a n-island model with constant size. Even if, for any [Formula: see text], it is always possible to find a size-change scenario for a panmictic population such that the marginal distribution of [Formula: see text] is exactly the same as in a n-island model with constant population size, we show that the joint distribution of the coalescence times [Formula: see text] for a sample of three genes contains enough information to distinguish between a panmictic population and a n-island model of constant size.

Genetic and morphological diversity of mouse lemurs (Microcebus spp.) in northern Madagascar: The discovery of a putative new species?

Tropical forests harbor extremely high levels of biological diversity and are quickly disappearing. Despite the increasingly recognized high rate of habitat loss, it is expected that new species will be discovered as more effort is put to document tropical biodiversity. Exploring under-studied regions is particularly urgent if we consider the rapid changes in habitat due to anthropogenic activities. Madagascar is known for its extraordinary biological diversity and endemicity. It is also threatened by habitat loss and fragmentation. It holds more than 100 endemic primate species (lemurs). Among these, Microcebus (mouse lemurs) is one of the more diverse genera. We sampled mouse lemurs from several sites across northern Madagascar, including forests never sampled before. We obtained morphological data from 99 Microcebus individuals; we extracted DNA from tissue samples of 42 individuals and amplified two mitochondrial loci (cytb and cox2) commonly used for species identification. Our findings update the distribution of three species (Microcebus tavaratra, Microcebus arnholdi, and Microcebus mamiratra), including a major increase in the distribution area of M. arnholdi. We also report the discovery of a new Microcebus lineage genetically related to M. arnholdi. Several complementary approaches suggest that the newly identified Microcebus lineage might correspond to a new putative species, to be confirmed or rejected with additional data. In addition, morphological analyses showed (a) clear phenotypic differences between M. tavaratra and M. arnholdi, but no clear differences between the new Microcebus lineage and the sister species M. arnholdi; and (b) a significant correlation between climatic variables and morphology, suggesting a possible relationship between species identity, morphology, and environment. By integrating morphological, climatic, genetic, and spatial data of two northern Microcebus species, we show that the spatial distribution of forest-dwelling species may be used as a proxy to reconstruct the past spatial changes in forest cover and vegetation type.

The IICR and the non-stationary structured coalescent: towards demographic inference with arbitrary changes in population structure.

In the last years, a wide range of methods allowing to reconstruct past population size changes from genome-wide data have been developed. At the same time, there has been an increasing recognition that population structure can generate genetic data similar to those produced under models of population size change. Recently, Mazet et al. (Heredity 116:362-371, 2016) showed that, for any model of population structure, it is always possible to find a panmictic model with a particular function of population size changes, having exactly the same distribution of T2 (the coalescence time for a sample of size two) as that of the structured model. They called this function IICR (Inverse Instantaneous Coalescence Rate) and showed that it does not necessarily correspond to population size changes under non-panmictic models. Besides, most of the methods used to analyse data under models of population structure tend to arbitrarily fix that structure and to minimise or neglect population size changes. Here, we extend the seminal work of Herbots (PhD thesis, University of London, 1994) on the structured coalescent and propose a new framework, the Non-Stationary Structured Coalescent (NSSC) that incorporates demographic events (changes in gene flow and/or deme sizes) to models of nearly any complexity. We show how to compute the IICR under a wide family of stationary and non-stationary models. As an example we address the question of human and Neanderthal evolution and discuss how the NSSC framework allows to interpret genomic data under this new perspective.

The IICR (inverse instantaneous coalescence rate) as a summary of genomic diversity: insights into demographic inference and model choice.

Several inferential methods using genomic data have been proposed to quantify and date population size changes in the history of species. At the same time an increasing number of studies have shown that population structure can generate spurious signals of population size change. Recently, Mazet et al. (2016) introduced, for a sample size of two, a time-dependent parameter, which they called the IICR (inverse instantaneous coalescence rate). The IICR is equivalent to a population size in panmictic models, but not necessarily in structured models. It is characterised by a temporal trajectory that suggests population size changes, as a function of the sampling scheme, even when the total population size was constant. Here, we extend the work of Mazet et al. (2016) by (i) showing how the IICR can be computed for any demographic model of interest, under the coalescent, (ii) applying this approach to models of population structure (1D and 2D stepping stone, split models, two- and three-island asymmetric gene flow, continent-island models), (iii) stressing the importance of the sampling strategy in generating different histories, (iv) arguing that IICR plots can be seen as summaries of genomic information that can thus be used for model choice or model exclusion (v) applying this approach to the question of admixture between humans and Neanderthals. Altogether these results are potentially important given that the widely used PSMC (pairwise sequentially Markovian coalescent) method of Li and Durbin (2011) estimates the IICR of the sample, not necessarily the history of the populations.

On some genetic consequences of social structure, mating systems, dispersal, and sampling.

Many species are spatially and socially organized, with complex social organizations and dispersal patterns that are increasingly documented. Social species typically consist of small age-structured units, where a limited number of individuals monopolize reproduction and exhibit complex mating strategies. Here, we model social groups as age-structured units and investigate the genetic consequences of social structure under distinct mating strategies commonly found in mammals. Our results show that sociality maximizes genotypic diversity, which contradicts the belief that social groups are necessarily subject to strong genetic drift and at high risk of inbreeding depression. Social structure generates an excess of genotypic diversity. This is commonly observed in ecological studies but rarely reported in population genetic studies that ignore social structure. This heterozygosity excess, when detected, is often interpreted as a consequence of inbreeding avoidance mechanisms, but we show that it can occur even in the absence of such mechanisms. Many seemly contradictory results from ecology and population genetics can be reconciled by genetic models that include the complexities of social species. We find that such discrepancies can be explained by the intrinsic properties of social groups and by the sampling strategies of real populations. In particular, the number of social groups and the nature of the individuals that compose samples (e.g., nonreproductive and reproductive individuals) are key factors in generating outbreeding signatures. Sociality is an important component of population structure that needs to be revisited by ecologists and population geneticists alike.

Climate change and human colonization triggered habitat loss and fragmentation in Madagascar.

The relative effect of past climate fluctuations and anthropogenic activities on current biome distribution is subject to increasing attention, notably in biodiversity hot spots. In Madagascar, where humans arrived in the last ~4 to 5,000 years, the exact causes of the demise of large vertebrates that cohabited with humans are yet unclear. The prevailing narrative holds that Madagascar was covered with forest before human arrival and that the expansion of grasslands was the result of human-driven deforestation. However, recent studies have shown that vegetation and fauna structure substantially fluctuated during the Holocene. Here, we study the Holocene history of habitat fragmentation in the north of Madagascar using a population genetics approach. To do so, we infer the demographic history of two northern Madagascar neighbouring, congeneric and critically endangered forest dwelling lemur species-Propithecus tattersalli and Propithecus perrieri-using population genetic analyses. Our results highlight the necessity to consider population structure and changes in connectivity in demographic history inferences. We show that both species underwent demographic fluctuations which most likely occurred after the mid-Holocene transition. While mid-Holocene climate change probably triggered major demographic changes in the two lemur species range and connectivity, human settlements that expanded over the last four millennia in northern Madagascar likely played a role in the loss and fragmentation of the forest cover.

Demographic inference using genetic data from a single individual: Separating population size variation from population structure.

The rapid development of sequencing technologies represents new opportunities for population genetics research. It is expected that genomic data will increase our ability to reconstruct the history of populations. While this increase in genetic information will likely help biologists and anthropologists to reconstruct the demographic history of populations, it also represents new challenges. Recent work has shown that structured populations generate signals of population size change. As a consequence it is often difficult to determine whether demographic events such as expansions or contractions (bottlenecks) inferred from genetic data are real or due to the fact that populations are structured in nature. Given that few inferential methods allow us to account for that structure, and that genomic data will necessarily increase the precision of parameter estimates, it is important to develop new approaches. In the present study we analyze two demographic models. The first is a model of instantaneous population size change whereas the second is the classical symmetric island model. We (i) re-derive the distribution of coalescence times under the two models for a sample of size two, (ii) use a maximum likelihood approach to estimate the parameters of these models (iii) validate this estimation procedure under a wide array of parameter combinations, (iv) implement and validate a model rejection procedure by using a Kolmogorov-Smirnov test, and a model choice procedure based on the AIC, and (v) derive the explicit distribution for the number of differences between two non-recombining sequences. Altogether we show that it is possible to estimate parameters under several models and perform efficient model choice using genetic data from a single diploid individual.

Genetic data suggest a natural prehuman origin of open habitats in northern Madagascar and question the deforestation narrative in this region.

The impact of climate change and anthropogenic deforestation on biodiversity is of growing concern worldwide. Disentangling how past anthropogenic and natural factors contributed to current biome distribution is thus a crucial issue to understand their complex interactions on wider time scales and to improve predictions and conservation strategies. This is particularly important in biodiversity hotspots, such as Madagascar, dominated by large open habitats whose origins are increasingly debated. Although a dominant narrative argues that Madagascar was originally entirely covered by woodlands, which were destroyed by humans, a number of recent studies have suggested that past climatic fluctuations played a major role in shaping current biome distributions well before humans arrived. Here, we address the question of the origin of open habitats in the Daraina region in northern Madagascar, using a multiproxy approach combining population genetics modeling and remote-sensing analyses. We show that (i) contrary to most regions of Madagascar, the forest cover in Daraina remained remarkably stable over the past 60 y, and (ii) the golden-crowned sifaka (Propithecus tattersalli), a forest-dwelling lemur, underwent a strong population contraction before the arrival of the first humans, hence excluding an anthropogenic cause. Prehuman Holocene droughts may have led to a significant increase of grasslands and a reduction in the species' habitat. This contradicts the prevailing narrative that land cover changes are necessarily anthropogenic in Madagascar but does not preclude the later role played by humans in other regions in which recent lemur bottlenecks have been observed.

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.