The Genome of the Endangered Dryas Monkey Provides New Insights into the Evolutionary History of the Vervets

The phylogenetic relationship and demographic history of rhesus macaques (Macaca mulatta) in subtropical and temperate regions, China

Pleistocene climatic oscillations exerted significant influences on the genetic structure and demography of rhesus macaque (Macaca mulatta) in eastern China. However, the evolutionary history of rhesus macaques in subtropical and temperate China remained unclear and/or controversial. Herein, we analyzed the autosomes, mitochondrial genomes, and Y-chromosomes from 84 individuals of Chinese rhesus macaque. The results revealed that (1) all individuals were clustered into pan-west and pan-east genetic groups, which exhibited Shaanxi Province as the northernmost region of western dispersal route of rhesus macaques in China; (2) in subtropical and temperate China, rhesus macaques were divided into four lineages (TH, DB, HS, and QL), and their divergence times corresponded to the Penultimate Glaciation (300-130 kya) and Last Glaciation (70-10 kya), respectively; (3) the individuals from Mt. Taihangshan (TH) are closely related to individuals from Mt. Dabashan (DB) in the autosomal tree, rather than individuals from Mt. Huangshan (HS) as indicated by the mitogenome tree, which supports the hypothesis that the ancestral rhesus macaques radiated into Mt. Taihangshan from Mt. Huangshan via Mt. Dabashan; and (4) the demographic scenario of the four lineages showed the ancestral rhesus macaques bottleneck and expansion corresponding to the suitable habitat reduction and expansion, which confirmed they had experienced northward recolonization and southward retreat events from Mt. Huangshan area via Northern China Plain to Northernmost China along with Pleistocene glacial cycles. This study provides a new insight into understanding how Pleistocene glaciation has influenced faunal diversity in subtropical and temperate China, especially for those exhibiting differential patterns of sex dispersal.

Complex Evolutionary History With Extensive Ancestral Gene Flow in an African Primate Radiation

Understanding the drivers of speciation is fundamental in evolutionary biology, and recent studies highlight hybridization as an important evolutionary force. Using whole-genome sequencing data from 22 species of guenons (tribe Cercopithecini), one of the world's largest primate radiations, we show that rampant gene flow characterizes their evolutionary history and identify ancient hybridization across deeply divergent lineages that differ in ecology, morphology, and karyotypes. Some hybridization events resulted in mitochondrial introgression between distant lineages, likely facilitated by cointrogression of coadapted nuclear variants. Although the genomic landscapes of introgression were largely lineage specific, we found that genes with immune functions were overrepresented in introgressing regions, in line with adaptive introgression, whereas genes involved in pigmentation and morphology may contribute to reproductive isolation. In line with reports from other systems that hybridization might facilitate diversification, we find that some of the most species-rich guenon clades are of admixed origin. This study provides important insights into the prevalence, role, and outcomes of ancestral hybridization in a large mammalian radiation.

Phylogenomics Reveals High Levels of Incomplete Lineage Sorting at the Ancestral Nodes of the Macaque Radiation

The genus Macaca includes 23 species assigned into 4 to 7 groups. It exhibits the largest geographic range and represents the most successful example of adaptive radiation of nonhuman primates. However, intrageneric phylogenetic relationships among species remain controversial and have not been resolved so far. In this study, we conducted a phylogenomic analysis on 16 newly generated and 8 published macaque genomes. We found strong evidence supporting the division of this genus into 7 species groups. Incomplete lineage sorting (ILS) was the primary factor contributing to the discordance observed among gene trees; however, we also found evidence of hybridization events, specifically between the ancestral arctoides/sinica and silenus/nigra lineages that resulted in the hybrid formation of the fascicularis/mulatta group. Combined with fossil data, our phylogenomic data were used to establish a scenario for macaque radiation. These findings provide insights into ILS and potential ancient introgression events that were involved in the radiation of macaques, which will lead to a better understanding of the rapid speciation occurring in nonhuman primates.

A global catalog of whole-genome diversity from 233 primate species

The rich diversity of morphology and behavior displayed across primate species provides an informative context in which to study the impact of genomic diversity on fundamental biological processes. Analysis of that diversity provides insight into long-standing questions in evolutionary and conservation biology and is urgent given severe threats these species are facing. Here, we present high-coverage whole-genome data from 233 primate species representing 86% of genera and all 16 families. This dataset was used, together with fossil calibration, to create a nuclear DNA phylogeny and to reassess evolutionary divergence times among primate clades. We found within-species genetic diversity across families and geographic regions to be associated with climate and sociality, but not with extinction risk. Furthermore, mutation rates differ across species, potentially influenced by effective population sizes. Lastly, we identified extensive recurrence of missense mutations previously thought to be human specific. This study will open a wide range of research avenues for future primate genomic research.

Genomics, Population Divergence, and Historical Demography of the World's Largest and Endangered Butterfly, The Queen Alexandra's Birdwing

The world's largest butterfly is the microendemic Papua New Guinean Ornithoptera alexandrae. Despite years of conservation efforts to protect its habitat and breed this up-to-28-cm butterfly, this species still figures as endangered in the IUCN Red List and is only known from two allopatric populations occupying a total of only ∼140 km². Here we aim at assembling reference genomes for this species to investigate its genomic diversity, historical demography and determine whether the population is structured, which could provide guidance for conservation programs attempting to (inter)breed the two populations. Using a combination of long and short DNA reads and RNA sequencing, we assembled six reference genomes of the tribe Troidini, with four annotated genomes of O. alexandrae and two genomes of related species Ornithoptera priamus and Troides oblongomaculatus. We estimated the genomic diversity of the three species, and we proposed scenarios for the historical population demography using two polymorphism-based methods taking into account the characteristics of low-polymorphic invertebrates. Indeed, chromosome-scale assemblies reveal very low levels of nuclear heterozygosity across Troidini, which appears to be exceptionally low for O. alexandrae (lower than 0.01%). Demographic analyses demonstrate low and steadily declining Ne throughout O. alexandrae history, with a divergence into two distinct populations about 10,000 years ago. These results suggest that O. alexandrae distribution has been microendemic for a long time. It should also make local conservation programs aware of the genomic divergence of the two populations, which should not be ignored if any attempt is made to cross the two populations.

Current advances in primate genomics: novel approaches for understanding evolution and disease

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

Genomic insights into inter- and intraspecific mating system shifts in Primulina

The mating system shift from outcrossing to selfing is one of the most frequent evolutionary trends in flowering plants. However, the genomic consequences of this shift remain poorly understood. Specifically, the relative importance of the demographic and genetic processes causing changes in genetic variation and selection efficacy associated with the evolution of selfing is unclear. Here we sequenced the genomes of two Primulina species with contrasting mating systems: P. eburnea (outcrossing) versus P. tabacum (outcrossing, mixed-mating and selfing populations). Whole-genome resequencing data were used to investigate the genomic consequences of mating system shifts within and between species. We found that highly selfing populations of P. tabacum display loss of genetic diversity, increased deleterious mutations, higher genomic burden and fewer adaptive substitutions. However, compared with outcrossing populations, mixed-mating populations did not display loss of genetic diversity and accumulation of genetic load. We find no evidence of population bottlenecks associated with the shift to selfing, which suggests that the genetic effects of selfing on N_e and possibly linked selection, rather than demographic history, are the primary drivers of diversity reduction in highly selfing populations. Our results highlight the importance of distinguishing the relative contribution of mating system and demography on the genomic consequences associated with mating system evolution in plants.

Evidence of selection in the uncoupling protein 1 gene region suggests local adaptation to solar irradiance in savannah monkeys (Chlorocebus spp.)

In the last 300 thousand years, the genus Chlorocebus expanded from equatorial Africa into the southernmost latitudes of the continent, where colder climate was a probable driver of natural selection. We investigated population-level genetic variation in the mitochondrial uncoupling protein 1 (UCP1) gene region-implicated in non-shivering thermogenesis (NST)-in 73 wild savannah monkeys from three taxa representing this southern expansion (Chlorocebus pygerythrus hilgerti, Chlorocebus cynosuros and Chlorocebus pygerythrus pygerythrus) ranging from Kenya to South Africa. We found 17 single nucleotide polymorphisms with extended haplotype homozygosity consistent with positive selective sweeps, 10 of which show no significant linkage disequilibrium with each other. Phylogenetic generalized least-squares modelling with ecological covariates suggest that most derived allele frequencies are significantly associated with solar irradiance and winter precipitation, rather than overall low temperatures. This selection and association with irradiance is demonstrated by a relatively isolated population in the southern coastal belt of South Africa. We suggest that sunbathing behaviours common to savannah monkeys, in combination with the strength of solar irradiance, may mediate adaptations to thermal stress via NST among savannah monkeys. The variants we discovered all lie in non-coding regions, some with previously documented regulatory functions, calling for further validation and research.

Rocks and clocks revised: New promises and challenges in dating the primate tree of life

In recent years, multiple technological and methodological advances have increased our ability to estimate phylogenies, leading to more accurate dating of the primate tree of life. Here we provide an overview of the limitations and potentials of some of these advancements and discuss how dated phylogenies provide the crucial temporal scale required to understand primate evolution. First, we review new methods, such as the total-evidence dating approach, that promise a better integration between the fossil record and molecular data. We then explore how the ever-increasing availability of genomic-level data for more primate species can impact our ability to accurately estimate timetrees. Finally, we discuss more recent applications of mutation rates to date divergence times. We highlight example studies that have applied these approaches to estimate divergence dates within primates. Our goal is to provide a critical overview of these new developments and explore the promises and challenges of their application in evolutionary anthropology.

Genetic load: genomic estimates and applications in non-model animals

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Divergence and introgression in small apes, the genus Hylobates, revealed by reduced representation sequencing

Gibbons of the genus Hylobates, which inhabit Southeast Asia, show great diversity and comprise seven to nine species. Natural hybridisation has been observed in several species contact zones, but the history and extent of hybridisation and introgression in possibly historical and the current contact zones remain unclear. To uncover Hylobates species phylogeny and the extent of introgression in their evolution, genotyping by random amplicon sequencing-direct (GRAS-Di) was applied to 47 gibbons, representing seven Hylobates species/subspecies and two outgroup gibbon species. Over 200,000 autosomal single-nucleotide variant sites were identified. The autosomal phylogeny supported that divergence from the mainland species began ~3.5 million years ago, and subsequently occurred among the Sundaic island species. Significant introgression signals were detected between H. lar and H. pileatus, H. lar and H. agilis and H. albibarbis and H. muelleri, which all are parapatric and form ongoing hybrid zones. Furthermore, the introgression signals were detected in every analysed individual of these species, indicating a relatively long history of hybridisation, which might have affected the entire gene pool. By contrast, signals of introgression were either not detected or doubtful in other species pairs living on different islands, indicating the rarity of hybridisation and introgression, even though the Sundaic islands were connected during the Pliocene and Pleistocene glacial events.

Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi

Inbreeding can lead to a loss of heterozygosity in a population and when combined with genetic drift may reduce the adaptive potential of a species. However, there is uncertainty about whether resequencing data can provide accurate and consistent inbreeding estimates. Here, we performed an in-depth inbreeding analysis for hihi (Notiomystis cincta), an endemic and nationally vulnerable passerine bird of Aotearoa New Zealand. We first focused on subsampling variants from a reference genome male, and found that low-density data sets tend to miss runs of homozygosity (ROH) in some places and overestimate ROH length in others, resulting in contrasting homozygosity landscapes. Low-coverage resequencing and 50 K SNP array densities can yield comparable inbreeding results to high-coverage resequencing approaches, but the results for all data sets are highly dependent on the software settings employed. Second, we extended our analysis to 10 hihi where low-coverage whole genome resequencing, RAD-seq and SNP array genotypes are available. We inferred ROH and individual inbreeding to evaluate the relative effects of sequencing depth versus SNP density on estimating inbreeding coefficients and found that high rates of missingness downwardly bias both the number and length of ROH. In summary, when using genomic data to evaluate inbreeding, studies must consider that ROH estimates are heavily dependent on analysis parameters, data set density and individual sequencing depth.

A natural history of Chlorocebus dryas from camera traps in Lomami National Park and its buffer zone, Democratic Republic of the Congo, with notes on the species status of Cercopithecus salongo

The natural history and taxonomic status of two central African primates, Cercopithecus dryas and Cercopithecus salongo have long been in question. Recent studies confirmed that C. dryas is a basal member of the savanna monkey clade, and that it prefers dense undergrowth in lowland rainforest. While these studies advanced our knowledge of this enigmatic species, key aspects of its natural history remain poorly documented. Furthermore, the lack of a field study that documents pelage patterns of both sexes and different age classes of C. dryas has led to a disagreement over the validity of C. salongo as a sister taxon to C. dryas. Using the results of two multi-strata camera trap surveys in Lomami National Park (Democratic Republic of the Congo) and its buffer zone, we conducted a third survey in the understory of degraded forest to accumulate videos of C. dryas/salongo. We used these videos to test the hypothesis that C. dryas and C. salongo are synonymous, and to assess the species' group composition, density, behavior and vocalizations. Camera traps revealed an ontogenetic change in pelage pattern that supports the view that C. salongo is the adult of C. dryas. Videos revealed that adult males develop a blue perineum and scrotum, and a red subcaudal patch, similar to other savanna monkeys. We provide a preliminary assessment of C. dryas' group composition, density, behavior, and vocalizations. This long-overlooked monkey is an exceptional member of the Chlorocebus clade, and all aspects of its biology require further investigation.

Morphological analysis of new Dryas Monkey specimens from the Central Congo Basin: Taxonomic considerations and an emended diagnosis

OBJECTIVES

The little known guenon Cercopithecus dryas has a controversial taxonomic history with some recognizing two taxa (C. dryas and C. salongo) instead of one. New adult specimens from the TL2 region of the central Congo Basin allow further assessment of C. dryas morphology and, along with CT scans of the juvenile holotype, provide ontogenetically stable comparisons across all C. dryas and "C. salongo" specimens for the first time.

MATERIALS AND METHODS

The skins and skulls of two newly acquired C. dryas specimens, male YPM MAM 16890 and female YPM MAM 17066, were compared to previously described C. dryas and "C. salongo" specimens, along with a broader guenon comparative sample (cranial sample n = 146, dental sample n = 102). Qualitative and quantitative assessments were made on the basis of commonly noted pelage features as well as craniodental characters in the form of shape ratios and multivariate discriminant analyses.

RESULTS

All C. dryas specimens, including the TL2 adults, are comparatively small in overall cranial size, have relatively small I¹ s, and display tall molar cusps; these osteological characters, along with pelage features, are shared with known "C. salongo" specimens. Discriminant analyses of dental features separate C. dryas/salongo specimens from all other guenons.

DISCUSSION

In addition to pelage-based evidence, direct osteological evidence suggests "C. salongo" is a junior synonym of C. dryas. Combined with molecular analyses suggesting C. dryas is most closely related to Chlorocebus spp., we emend the species diagnosis and support its transfer to Chlorocebus or possibly a new genus to reflect its distinctiveness.

Using local knowledge and camera traps to investigate occurrence and habitat preference of an Endangered primate: the endemic dryas monkey in the Democratic Republic of the Congo

The Endangered dryas monkey Cercopithecus dryas, endemic to the Democratic Republic of the Congo, is one of Africa's most enigmatic primates. The discovery of a dryas monkey killed by a hunter in the buffer zone of Lomami National Park in 2014 prompted field research on the species’ distribution, habitat use and stratum preference. We used local knowledge to determine the distribution of this species and to select sites for camera-trap surveys in Lomami National Park and its buffer zone. We employed a multi-strata (0–29 m) camera-trap placement technique to determine habitat use at Camp Bartho in Lomami National Park and Bafundo Forest in the Park's buffer zone. We confirmed the occurrence of the dryas monkey at seven locations over a total area of 3,453 km2, in both the Park and its buffer zone. Dryas monkeys were detected most frequently (2.22 events/100 trap-days) in disturbed areas of Bafundo Forest and less in mature forest in Camp Bartho (0.82 events/100 trap-days). Dryas monkeys appear to prefer structurally complex understories and forest edges. We found that camera traps at 2–10 m above ground over at least 365 trap-days are required to determine if the species is present. We recommend utilizing local knowledge and using this species-specific camera-trap method in other areas of the central Congo basin to determine the wider distribution of the dryas monkey.

Population Genomics Reveals Incipient Speciation, Introgression, and Adaptation in the African Mona Monkey (Cercopithecus mona)

Guenons (tribe Cercopithecini) are the most widely distributed nonhuman primate in the tropical forest belt of Africa and show considerable phenotypic, taxonomic, and ecological diversity. However, genomic information for most species within this group is still lacking. Here, we present a high-quality de novo genome (total 2.90 Gb, contig N50 equal to 22.7 Mb) of the mona monkey (Cercopithecus mona), together with genome resequencing data of 13 individuals sampled across Nigeria. Our results showed differentiation between populations from East and West of the Niger River ∼84 ka and potential ancient introgression in the East population from other mona group species. The PTPRK, FRAS1, BNC2, and EDN3 genes related to pigmentation displayed signals of introgression in the East population. Genomic scans suggest that immunity genes such as AKT3 and IL13 (possibly involved in simian immunodeficiency virus defense), and G6PD, a gene involved in malaria resistance, are under positive natural selection. Our study gives insights into differentiation, natural selection, and introgression in guenons.

Genetic Diversity, Inbreeding Level, and Genetic Load in Endangered Snub-Nosed Monkeys (Rhinopithecus)

The snub-nosed monkey genus (Rhinopithecus) comprises five closely related species (R. avunculus, R. bieti, R. brelichi, R. roxellana, and R. strykeri). All are among the world's rarest and most endangered primates. However, the genomic impact associated with their population decline remains unknown. We analyzed population genomic data of all five snub-nosed monkey species to assess their genetic diversity, inbreeding level, and genetic load. For R. roxellana, R. bieti, and R. strykeri, population size is positively correlated with genetic diversity and negatively correlated with levels of inbreeding. Other species, however, which possess small population sizes, such as R. brelichi and R. avunculus, show high levels of genetic diversity and low levels of genomic inbreeding. Similarly, in the three populations of R. roxellana, the Shennongjia population, which possesses the lowest population size, displays a higher level of genetic diversity and lower level of genomic inbreeding. These findings suggest that although R. brelichi and R. avunculus and the Shennongjia population might be at risk, it possess significant genetic diversity and could thus help strengthen their long-term survival potential. Intriguingly, R. roxellana with large population size possess high genetic diversity and low level of genetic load, but they show the highest recent inbreeding level compared with the other snub-nosed monkeys. This suggests that, despite its large population size, R. roxellana has likely been experiencing recent inbreeding, which has not yet affected its mutational load and fitness. Analyses of homozygous-derived deleterious mutations identified in all snub-nosed monkey species indicate that these mutations are affecting immune, especially in smaller population sizes, indicating that the long-term consequences of inbreeding may be resulting in an overall reduction of immune capability in the snub-nosed monkeys, which could provide a dramatic effect on their long-term survival prospects. Altogether, our study provides valuable information concerning the genomic impact of population decline of the snub-nosed monkeys. We revealed multiple counterintuitive and unexpected patterns of genetic diversity in small and large population, which will be essential for conservation management of these endangered species.

The Diversity of Primates: From Biomedicine to Conservation Genomics

Until now, the field of primate genomics has focused on two major themes: understanding human evolution and advancing biomedical research. We propose that it is now time for a third theme to receive attention: conservation genomics. As a result of anthropogenic effects, the majority of primate species have become threatened with extinction. A more robust primate conservation genomics will allow for genetically informed population management. Thanks to a steady decline in the cost of sequencing, it has now become feasible to sequence whole primate genomes at the population level. Furthermore, technological advances in noninvasive genomic methods have made it possible to acquire genome-scale data from noninvasive biomaterials. Here, we review recent advances in the analysis of primate diversity, with a focus on genomic data sets across the radiation.

Fifty Years of Research on European Mink Mustela lutreola L., 1761 Genetics: Where Are We Now in Studies on One of the Most Endangered Mammals?

The purpose of this review is to present the current state of knowledge about the genetics of European mink Mustela lutreola L., 1761, which is one of the most endangered mammalian species in the world. This article provides a comprehensive description of the studies undertaken over the last 50 years in terms of cytogenetics, molecular genetics, genomics (including mitogenomics), population genetics of wild populations and captive stocks, phylogenetics, phylogeography, and applied genetics (including identification by genetic methods, molecular ecology, and conservation genetics). An extensive and up-to-date review and critical analysis of the available specialist literature on the topic is provided, with special reference to conservation genetics. Unresolved issues are also described, such as the standard karyotype, systematic position, and whole-genome sequencing, and hotly debated issues are addressed, like the origin of the Southwestern population of the European mink and management approaches of the most distinct populations of the species. Finally, the most urgent directions of future research, based on the research questions arising from completed studies and the implementation of conservation measures to save and restore M. lutreola populations, are outlined. The importance of the popularization of research topics related to European mink genetics among scientists is highlighted.

An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony

High density single nucleotide polymorphism (SNP) arrays allow large numbers of individuals to be rapidly and cost-effectively genotyped at large numbers of genetic markers. However, despite being widely used in studies of humans and domesticated plants and animals, SNP arrays are lacking for most wild organisms. We developed a custom 85K Affymetrix Axiom array for an intensively studied pinniped, the Antarctic fur seal (Arctocephalus gazella). SNPs were discovered from a combination of genomic and transcriptomic resources and filtered according to strict criteria. Out of a total of 85,359 SNPs tiled on the array, 75,601 (88.6%) successfully converted and were polymorphic in 270 animals from a breeding colony at Bird Island in South Georgia. Evidence was found for inbreeding, with three genomic inbreeding coefficients being strongly intercorrelated and the proportion of the genome in runs of homozygosity being non-zero in all individuals. Furthermore, analysis of genomic relatedness coefficients identified previously unknown first-degree relatives and multiple second-degree relatives among a sample of ostensibly unrelated individuals. Such “cryptic relatedness” within fur seal breeding colonies may increase the likelihood of consanguineous matings and could therefore have implications for understanding fitness variation and mate choice. Finally, we demonstrate the cross-amplification potential of the array in three related pinniped species. Overall, our SNP array will facilitate future studies of Antarctic fur seals and has the potential to serve as a more general resource for the wider pinniped research community.