Maja P. Greminger

Marked Population Structure and Recent Migration in the Critically Endangered Sumatran Orangutan (Pongo abelii)

A multitude of factors influence how natural populations are genetically structured, including dispersal barriers, inhomogeneous habitats, and social organization. Such population subdivision is of special concern in endangered species, as it may lead to reduced adaptive potential and inbreeding in local subpopulations, thus increasing the risk of future extinctions. With only 6600 animals left in the wild, Sumatran orangutans (Pongo abelii) are among the most endangered, but also most enigmatic, great ape species. In order to infer the fine-scale population structure and connectivity of Sumatran orangutans, we analyzed the most comprehensive set of samples to date, including mitochondrial hyper-variable region I haplotypes for 123 individuals and genotypes of 27 autosomal microsatellite markers for 109 individuals. For both mitochondrial and autosomal markers, we found a pronounced population structure, caused by major rivers, mountain ridges, and the Toba caldera. We found that genetic diversity and corresponding long-term effective population size estimates vary strongly among sampling regions for mitochondrial DNA, but show remarkable similarity for autosomal markers, hinting at male-driven long-distance gene flow. In support of this, we identified several individuals that were most likely sired by males originating from other genetic clusters. Our results highlight the effect of natural barriers in shaping the genetic structure of great ape populations, but also point toward important dispersal corridors on northern Sumatra that allow for genetic exchange.

Development of SNP markers identifying European wildcats, domestic cats, and their admixed progeny

Introgression can be an important evolutionary force but it can also lead to species extinction and as such is a crucial issue for species conservation. However, introgression is difficult to detect, morphologically as well as genetically. Hybridization with domestic cats (Felis silvestris catus) is a major concern for the conservation of European wildcats (Felis s. silvestris). The available morphologic and genetic markers for the two Felis subspecies are not sufficient to reliably detect hybrids beyond first generation. Here we present a single nucleotide polymorphism (SNP) based approach that allows the identification of introgressed individuals. Using high‐throughput sequencing of reduced representation libraries we developed a diagnostic marker set containing 48 SNPs (Fst > 0.8) which allows the identification of wildcats, domestic cats, their hybrids and backcrosses. This allows assessing introgression rate in natural wildcat populations and is key for a better understanding of hybridization processes.

Parentage-based pedigree reconstruction reveals female matrilineal clusters and male-biased dispersal in nongregarious Asian great apes, the Bornean orang-utans (Pongo pygmaeus)

Philopatry and sex‐biased dispersal have a strong influence on population genetic structure, so the study of species dispersal patterns and evolutionary mechanisms shaping them are of great interest. Particularly nongregarious mammalian species present an underexplored field of study: despite their lower levels of sociality compared to group‐living species, interactions among individuals do occur, providing opportunities for cryptic kin selection. Among the least gregarious primates are orang‐utans (genus: Pongo), in which preferential associations among females have nevertheless been observed, but for which the presence of kin structures was so far unresolved because of the equivocal results of previous genetic studies. To clarify relatedness and dispersal patterns in orang‐utans, we examined the largest longitudinal set of individuals with combined genetic, spatial and behavioural data. We found that males had significantly higher mitochondrial DNA (mtDNA) variation and more unique haplotypes, thus underscoring their different maternal ancestries compared to females. Moreover, pedigree reconstruction based on 24 highly polymorphic microsatellite markers and mtDNA haplotypes demonstrated the presence of three matrilineal clusters of generally highly related females with substantially overlapping ranges. In orang‐utans and possibly other nongregarious species, comparing average biparental relatedness (r) of males and females to infer sex‐biased dispersal is extremely problematic. This is because the opportunistic sampling regime frequently employed in nongregarious species, combined with overlapping space use of distinct matrilineal clusters, leads to a strong downward bias when mtDNA lineage membership is ignored. Thus, in nongregarious species, correct inferences of dispersal can only be achieved by combining several genetic approaches with detailed spatial information.

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

BackgroundHigh-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.ResultsWe 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.ConclusionsOur 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.

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.

A multiplex-system to target 16 male-specific and 15 autosomal genetic markers for orang-utans (genus: Pongo)

Genetic studies of dispersal on local spatial and short temporal scales require a large number of autosomal microsatellites. However, the study of dispersal over large spatial scales and the resolution of deep evolutionary histories require marker systems that are preferentially inherited through the male or female line. Addressing such questions in endangered orang-utans (genus: Pongo) bears significant relevance to species conservation, as habitat destruction and fragmentation pose a significant threat to the whole genus. Here, we report 16 male-specific markers (nine human-derived microsatellites, six single nucleotide and one insertion-deletion polymorphisms), and 15 novel Pongo-derived autosomal microsatellite loci. All 31 markers can be amplified in four multiplex polymerase chain reactions even in DNA derived from faecal material. The markers can be applied to studying a wide range of important questions in this genus, such as conservation genetics, social structure, phylogeny and phylogeography.

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

Significance The rate of genomic adaptation is determined by the rate of environmental change, the availability of beneficial mutations, and the efficiency of positive selection. The relative importance of these factors has been actively discussed. We address the questions using whole genome sequences of great apes, which have very different population sizes whereas their genomic architectures are highly similar. We infer that the impact of selection on the genomic diversity of a species increases with the effective population size, most likely due to the differential influx rate of beneficial mutations. This explanation is, among other possibilities, expected if adaptive evolution is limited by the waiting time for new favorable mutations in 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.

Development of polymorphic microsatellite markers for the dung fly (Sepsis cynipsea).

The polyandrous fly Sepsis cynipsea has been used extensively in studies of sexual selection and local adaptation. We isolated and characterized 11 novel microsatellite markers for S. cynipsea from a genomic library and screened 32 flies for polymorphism. All microsatellite markers show high allelic diversity with an average of 9.64 alleles per locus. Two microsatellites were found likely to be X‐linked. These novel markers will significantly advance studies of sexual selection and evolutionary genetics of S. cynipsea and related species, especially given the low numbers of markers currently available in this family.