Past Connectivity but Recent Inbreeding in Cross River Gorillas Determined Using Whole Genomes from Single Hairs

Estimating realized relatedness in free-ranging macaques by inferring identity-by-descent segments

Biological relatedness is a key consideration in studies of behavior, population structure, and trait evolution. Except for parent-offspring dyads, pedigrees capture relatedness imperfectly. The number and length of identical-by-descent DNA segments (IBD) yield the most precise relatedness estimates. Here, we leverage different methods for estimating IBD segments from low-depth whole genome resequencing data to demonstrate the feasibility and value of resolving fine-scaled gradients of relatedness in free-living animals. Using primarily 4 to 6× depth data from a rhesus macaque (Macaca mulatta) population with long-term pedigree data, we show that we can infer the number and length of IBD segments across the genome with high accuracy even at 0.5× sequencing depth. In line with expectations based on simulation, the resulting estimates demonstrate substantial variation in genetic relatedness within kin classes, leading to overlapping distributions between kin classes. By comparing the IBD-based estimates with pedigree and short tandem repeat-based methods, we show that IBD estimates are more reliable and provide more detailed information on kinship. The inferred IBD segments also identify cryptic genetic relatives not represented in the pedigree and reveal elevated recombination rates in females relative to males, which enables the majority of close maternal and paternal kin to be distinguished with genotype data alone. Our findings represent a breakthrough in the ability to study the predictors and consequences of genetic relatedness in natural populations, contributing to our understanding of a fundamental component of population structure in the wild.

Ecological and anthropogenic effects on the genomic diversity of lemurs in Madagascar

Ecological variation and anthropogenic landscape modification have had key roles in the diversification and extinction of mammals in Madagascar. Lemurs represent a radiation with more than 100 species, constituting roughly one-fifth of the primate order. Almost all species of lemurs are threatened with extinction, but little is known about their genetic diversity and demographic history. Here, we analyse high-coverage genome-wide resequencing data from 162 unique individuals comprising 50 species of Lemuriformes, including multiple individuals from most species. Genomic diversity varies widely across the infraorder and yet is broadly consistent among individuals within species. We show widespread introgression in multiple genera and generally high levels of genomic diversity likely resulting from allele sharing that occurred during periods of connectivity and fragmentation during climatic shifts. We find distinct patterns of demographic history in lemurs across the ecogeographic regions of Madagascar within the last million years. Within the past 2,000 years, lemurs underwent major declines in effective population size that corresponded to the timing of human population expansion in Madagascar. In multiple regions of the island, we identified chronological trajectories of inbreeding that are consistent across genera and species, suggesting localized effects of human activity. Our results show how the extraordinary diversity of these long-neglected, endangered primates has been influenced by ecological and anthropogenic factors.

Comparing Genetic Ne Reconstructions Over Time With Long-Time Wolf Monitoring Data in Two Populations

Many methods are now available to calculate N e , but their performance varies depending on assumptions. Although simulated data are useful to discover certain types of bias, real empirical data supported by detailed known population histories allow us to discern how well methods perform with actual messy and complex data. Here, we focus on two genomic data sets of grey wolf populations for which population size changes of the past 40-120 years are well documented. We use this background to explore in what detail we can retrieve the known population history from these populations, in the light of pitfalls relating to population history, sampling design and the change in the spatial scale at which N e is estimated as we go further back in time. The Scandinavian wolf population was founded in the early 1980s from a few individuals and has gradually expanded up to 510 wolves. Although the founder event of the Scandinavian population was detected by GONE, the founding effective population size was strongly overestimated when the most recent samples were used, but less so when older samples were considered. Nevertheless, the present-day N e corresponds to theoretical expectations. The western Great Lakes wolf population of Minnesota is the only population in the contiguous United States that persisted throughout the 20th century, surviving intense persecution. We found a good concordance between the estimated N e and trends in census size data, but the reconstruction of N e clearly highlights the difficulty of interpreting results in spatially structured populations that underwent demographic fluctuations.

Ghosts of extinct apes: genomic insights into African hominid evolution

We are accustomed to regular announcements of new hominin fossils. There are now some 6000 hominin fossils, and up to 31 species. However, where are the announcements of African ape fossils? The answer is that there are almost none. Our knowledge of African ape evolution is based entirely on genomic analyses, which show that extant diversity is very young. This contrasts with the extensive and deep diversity of hominins known from fossils. Does this difference point to low and late diversification of ape lineages, or high rates of extinction? The comparative evolutionary dynamics of African hominids are central to interpreting living ape adaptations, as well as understanding the patterns of hominin evolution and the nature of the last common ancestor.

Comparative genomic analyses provide new insights into evolutionary history and conservation genomics of gorillas

Genome sequencing is a powerful tool to understand species evolutionary history, uncover genes under selection, which could be informative of local adaptation, and infer measures of genetic diversity, inbreeding and mutational load that could be used to inform conservation efforts. Gorillas, critically endangered primates, have received considerable attention and with the recently sequenced Bwindi mountain gorilla population, genomic data is now available from all gorilla subspecies and both mountain gorilla populations. Here, we reanalysed this rich dataset with a focus on evolutionary history, local adaptation and genomic parameters relevant for conservation. We estimate a recent split between western and eastern gorillas of 150,000-180,000 years ago, with gene flow around 20,000 years ago, primarily between the Cross River and Grauer's gorilla subspecies. This gene flow event likely obscures evolutionary relationships within eastern gorillas: after excluding putatively introgressed genomic regions, we uncover a sister relationship between Virunga mountain gorillas and Grauer's gorillas to the exclusion of Bwindi mountain gorillas. This makes mountain gorillas paraphyletic. Eastern gorillas are less genetically diverse and more inbred than western gorillas, yet we detected lower genetic load in the eastern species. Analyses of indels fit remarkably well with differences in genetic diversity across gorilla taxa as recovered with nucleotide diversity measures. We also identified genes under selection and unique gene variants specific for each gorilla subspecies, encoding, among others, traits involved in immunity, diet, muscular development, hair morphology and behavior. The presence of this functional variation suggests that the subspecies may be locally adapted. In conclusion, using extensive genomic resources we provide a comprehensive overview of gorilla genomic diversity, including a so-far understudied Bwindi mountain gorilla population, identify putative genes involved in local adaptation, and detect population-specific gene flow across gorilla species.

An empirical test of the estimation of historical effective population size using Drosophila melanogaster

The availability of a large number of high-density markers (SNPs) allows the estimation of historical effective population size (Ne ) from linkage disequilibrium between loci. A recent refinement of methods to estimate historical Ne from the recent past has been shown to be rather accurate with simulation data. The method has also been applied to real data for numerous species. However, the simulation data cannot encompass all the complexities of real genomes, and the performance of any estimation method with real data is always uncertain, as the true demography of the populations is not known. Here, we carried out an experimental design with Drosophila melanogaster to test the method with real data following a known demographic history. We used a population maintained in the laboratory with a constant census size of about 2800 individuals and subjected the population to a drastic decline to a size of 100 individuals. After a few generations, the population was expanded back to the previous size and after a few further generations again expanded to twice the initial size. Estimates of historical Ne were obtained with the software GONE both for autosomal and X chromosomes from samples of 17 individuals sequenced for the whole genome. Estimates of the historical effective size were able to infer the patterns of changes that occurred in the populations showing generally good performance of the method. We discuss the limitations of the method and the application of the software carried out so far.