Primate phylogenetic relationships and divergence dates inferred from complete mitochondrial genomes

Conservation of a Chromosome 8 Inversion and Exon Mutations Confirm Common Gulonolactone Oxidase Gene Evolution Among Primates, Including H. Neanderthalensis

Ascorbic acid functions as an antioxidant and facilitates other biochemical processes such as collagen triple helix formation, and iron uptake by cells. Animals which endogenously produce ascorbic acid have a functional gulonolactone oxidase gene (GULO); however, humans have a GULO pseudogene (GULOP) and depend on dietary ascorbic acid. In this study, the conservation of GULOP sequences in the primate haplorhini suborder were investigated and compared to the GULO sequences belonging to the primates strepsirrhini suborder. Phylogenetic analysis suggested that the conserved GULOP exons in the haplorhini primates experienced a high rate of mutations following the haplorhini/strepsirrhini divergence. This high mutation rate has decreased during the evolution of the haplorhini primates. Additionally, indels of the haplorhini GULOP sequences were conserved across the suborder. A separate analysis for GULO sequences and well-conserved GULOP sequences focusing on placental mammals identified an in-frame GULO sequence in the Brazilian guinea pig, and a potential GULOP sequence in the pika. Similar to haplorhini primates, the guinea pig and lagomorph species have experienced a high substitution rate when compared to the mammals used in this study. A shared synteny to examine the conservation of local genes near GULO/GULOP identified a conserved inversion around the GULO/GULOP locus between the haplorhini and strepsirrhini primates. Fischer's exact test did not support an association between GULOP and the chromosomal inversion. Mauve alignment showed that the inversion of the length of the syntenic block that the GULO/GULOP genes belonged to was variable. However, there were frequent rearrangements around ~ 2 million base pairs adjacent to GULOP involving the KIF13B and MSRA genes. These data may suggest that genes acquiring deleterious mutations in the coding sequence may respond to these deleterious mutations with rapid substitution rates.

Synchronized Expansion and Contraction of Olfactory, Vomeronasal, and Taste Receptor Gene Families in Hystricomorph Rodents

Chemical senses, including olfaction, pheromones, and taste, are crucial for the survival of most animals. There has long been a debate about whether different types of senses might influence each other. For instance, primates with a strong sense of vision are thought to have weakened olfactory abilities, although the oversimplified trade-off theory is now being questioned. It is uncertain whether such interactions between different chemical senses occur during evolution. To address this question, we examined four receptor gene families related to olfaction, pheromones, and taste: olfactory receptor (OR), vomeronasal receptor type 1 and type 2 (V1R and V2R), and bitter taste receptor (T2R) genes in Hystricomorpha, which is morphologically and ecologically the most diverse group of rodents. We also sequenced and assembled the genome of the grasscutter, Thryonomys swinderianus. By examining 16 available genome assemblies alongside the grasscutter genome, we identified orthologous gene groups among hystricomorph rodents for these gene families to separate the gene gain and loss events in each phylogenetic branch of the Hystricomorpha evolutionary tree. Our analysis revealed that the expansion or contraction of the four gene families occurred synchronously, indicating that when one chemical sense develops or deteriorates, the others follow suit. The results also showed that V1R/V2R genes underwent the fastest evolution, followed by OR genes, and T2R genes were the most evolutionarily stable. This variation likely reflects the difference in ligands of V1R/V2Rs, ORs, and T2Rs: species-specific pheromones, environment-based scents, and toxic substances common to many animals, respectively.

Structurally divergent and recurrently mutated regions of primate genomes

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

Generative modeling of biological shapes and images using a probabilistic α-shape sampler

Understanding morphological variation is an important task in many areas of computational biology. Recent studies have focused on developing computational tools for the task of sub-image selection which aims at identifying structural features that best describe the variation between classes of shapes. A major part in assessing the utility of these approaches is to demonstrate their performance on both simulated and real datasets. However, when creating a model for shape statistics, real data can be difficult to access and the sample sizes for these data are often small due to them being expensive to collect. Meanwhile, the current landscape of generative models for shapes has been mostly limited to approaches that use black-box inference-making it difficult to systematically assess the power and calibration of sub-image models. In this paper, we introduce the α -shape sampler: a probabilistic framework for generating realistic 2D and 3D shapes based on probability distributions which can be learned from real data. We demonstrate our framework using proof-of-concept examples and in two real applications in biology where we generate (i) 2D images of healthy and septic neutrophils and (ii) 3D computed tomography (CT) scans of primate mandibular molars. The α -shape sampler R package is open-source and can be downloaded at https://github.com/lcrawlab/ashapesampler.

Pan-primate studies of age and sex

Age and sex have a profound effect on cytosine methylation levels in humans and many other species. Here we analyzed DNA methylation profiles of 2400 tissues derived from 37 primate species including 11 haplorhine species (baboons, marmosets, vervets, rhesus macaque, chimpanzees, gorillas, orangutan, humans) and 26 strepsirrhine species (suborders Lemuriformes and Lorisiformes). From these we present here, pan-primate epigenetic clocks which are highly accurate for all primates including humans (age correlation R = 0.98). We also carried out in-depth analysis of baboon DNA methylation profiles and generated five epigenetic clocks for baboons (Olive-yellow baboon hybrid), one of which, the pan-tissue epigenetic clock, was trained on seven tissue types (fetal cerebral cortex, adult cerebral cortex, cerebellum, adipose, heart, liver, and skeletal muscle) with ages ranging from late fetal life to 22.8 years of age. Using the primate data, we characterize the effect of age and sex on individual cytosines in highly conserved regions. We identify 11 sex-related CpGs on autosomes near genes (POU3F2, CDYL, MYCL, FBXL4, ZC3H10, ZXDC, RRAS, FAM217A, RBM39, GRIA2, UHRF2). Low overlap can be observed between age- and sex-related CpGs. Overall, this study advances our understanding of conserved age- and sex-related epigenetic changes in primates, and provides biomarkers of aging for all primates.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

Deciphering the mandibular shape variation in a group of Malagasy primates using Fourier outline analysis

OBJECTIVES

Among living Malagasy primates, the family Lemuridae has previously been recognized as presenting a higher mandibular morphological variation than other families. We conducted a quantitative analysis of mandibular size and shape within the five genera (Lemur, Eulemur, Hapalemur, Prolemur, and Varecia) associated with a set of covariables that could explain this variation.

MATERIALS AND METHODS

We used Fourier outline analysis on the left hemimandible of 182 specimens covering the Lemuridae family. The influence of the phylogeny but also seven covariables (genus, diet, sex, sexual behavior, mating system, ecoregion, and forest type) on mandibular variation was examined using multivariate statistics and model selection.

RESULTS

Our results indicate that the high level of morphological variation within the family, associated with a phylogenetic effect and differences in diet, is due to a strong distinction between the genera Prolemur and Hapalemur and the other genera of the family. A second analysis, correcting this strong effect, indicates that mandibular shape variation is influenced not only by the phylogeny and the diet but by a combination of all the covariables.

DISCUSSION

The analysis of morphological variation is a powerful tool with major applications, both for the estimation of biological diversity and for the understanding of the fundamental parameters of species' ecology. Our work indicates that, if mandibular shape variation is mainly driven by dietary adaptation, other variables describing ecology and habitat should be considered and taken into account for an integrative understanding of species resources and the establishment of conservation measures.

Exploring objective feature sets in constructing the evolution relationship of animal genome sequences

BACKGROUND

Exploring evolution regularities of genome sequences and constructing more objective species evolution relationships at the genomic level are high-profile topics. Based on the evolution mechanism of genome sequences proposed in our previous research, we found that only the 8-mers containing CG or TA dinucleotides correlate directly with the evolution of genome sequences, and the relative frequency rather than the actual frequency of these 8-mers is more suitable to characterize the evolution of genome sequences.

RESULT

Therefore, two types of feature sets were obtained, they are the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers. The evolution relationships of mammals and reptiles were constructed by the relative frequency set of CG1 + CG2 8-mers, and two types of evolution relationships of insects were constructed by the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers respectively. Through comparison and analysis, we found that evolution relationships are consistent with the known conclusions. According to the evolution mechanism, we considered that the evolution relationship constructed by CG1 + CG2 8-mers reflects the evolution state of genome sequences in current time, and the evolution relationship constructed by TA1 + TA2 8-mers reflects the evolution state in the early stage.

CONCLUSION

Our study provides objective feature sets in constructing evolution relationships at the genomic level.

Museum DNA reveals a new, potentially extinct species of rinkhals (Serpentes: Elapidae: Hemachatus) from the Eastern Highlands of Zimbabwe

Genetic information plays a pivotal role in species recognition and delimitation, but rare or extinct animals can be difficult to obtain genetic samples from. While natural history wet collections have proven invaluable in the description of novel species, the use of these historical samples in genetic studies has been greatly impeded by DNA degradation, especially because of formalin-fixation prior to preservation. Here, we use recently developed museum genomics approaches to determine the status of an isolated population of the elapid snake genus Hemachatus from Zimbabwe. We used multiple digestion phases followed by single strand sequencing library construction and hybridisation capture to obtain 12S and 16S rDNA sequences from a poorly preserved tissue sample of this population. Phylogenetic and morphological analyses in an integrated taxonomic framework demonstrate that the Zimbabwean rinkhals population represents an old and highly distinct lineage, which we describe as a new species, Hemachatus nyangensis sp. nov. Our phylogenetic dating analysis is compatible with venom spitting having evolved in response to the threat posed by early hominins, although more data are required for a robust test of this hypothesis. This description demonstrates the power of museum genomics in revealing rare or even extinct species: Hemachatus from Zimbabwe are only known from a small area of the Eastern Highlands known for high endemism. No living specimens have been seen since the 1980s, most likely due to dramatic land-use changes in the Eastern Highlands, suggesting that the species could be extinct. In view of its recognition as a highly distinct lineage, urgent action is required to determine whether any populations survive, and to safeguard remaining habitat.

Allo: Accurate allocation of multi-mapped reads enables regulatory element analysis at repeats

Transposable elements (TEs) and other repetitive regions have been shown to contain gene regulatory elements, including transcription factor binding sites. Unfortunately, regulatory elements harbored by repeats have proven difficult to characterize using short-read sequencing assays such as ChIP-seq or ATAC-seq. Most regulatory genomics analysis pipelines discard "multi-mapped" reads that align equally well to multiple genomic locations. Since multi-mapped reads arise predominantly from repeats, current analysis pipelines fail to detect a substantial portion of regulatory events that occur in repetitive regions. To address this shortcoming, we developed Allo, a new approach to allocate multi-mapped reads in an efficient, accurate, and user-friendly manner. Allo combines probabilistic mapping of multi-mapped reads with a convolutional neural network that recognizes the read distribution features of potential peaks, offering enhanced accuracy in multi-mapping read assignment. Allo also provides read-level output in the form of a corrected alignment file, making it compatible with existing regulatory genomics analysis pipelines and downstream peak-finders. In a demonstration application on CTCF ChIP-seq data, we show that Allo results in the discovery of thousands of new CTCF peaks. Many of these peaks contain the expected cognate motif and/or serve as TAD boundaries. We additionally apply Allo to a diverse collection of ENCODE ChIP-seq datasets, resulting in multiple previously unidentified interactions between transcription factors and repetitive element families. Finally, we show that Allo may be particularly effective in identifying ChIP-seq peaks in younger TEs, which hold evolutionary significance due to their emergence during human evolution from primates.

Domestication in dry-cured meat Penicillium fungi: Convergent specific phenotypes and horizontal gene transfers without strong genetic subdivision

Some fungi have been domesticated for food production, with genetic differentiation between populations from food and wild environments, and food populations often acquiring beneficial traits through horizontal gene transfers (HGTs). Studying their adaptation to human-made substrates is of fundamental and applied importance for understanding adaptation processes and for further strain improvement. We studied here the population structures and phenotypes of two distantly related Penicillium species used for dry-cured meat production, P. nalgiovense, the most common species in the dry-cured meat food industry, and P. salamii, used locally by farms. Both species displayed low genetic diversity, lacking differentiation between strains isolated from dry-cured meat and those from other environments. Nevertheless, the strains collected from dry-cured meat within each species displayed slower proteolysis and lipolysis than their wild conspecifics, and those of P. nalgiovense were whiter. Phenotypically, the non-dry-cured meat strains were more similar to their sister species than to their conspecific dry-cured meat strains, indicating an evolution of specific phenotypes in dry-cured meat strains. A comparison of available Penicillium genomes from various environments revealed HGTs, particularly between P. nalgiovense and P. salamii (representing almost 1.5 Mb of cumulative length). HGTs additionally involved P. biforme, also found in dry-cured meat products. We further detected positive selection based on amino acid changes. Our findings suggest that selection by humans has shaped the P. salamii and P. nalgiovense populations used for dry-cured meat production, which constitutes domestication. Several genetic and phenotypic changes were similar in P. salamii, P. nalgiovense and P. biforme, indicating convergent adaptation to the same human-made environment. Our findings have implications for fundamental knowledge on adaptation and for the food industry: the discovery of different phenotypes and of two mating types paves the way for strain improvement by conventional breeding, to elucidate the genomic bases of beneficial phenotypes and to generate diversity.

Identification of 3 novel herpesviruses in prosimians with lymphoproliferative disease

Although many studies have characterized catarrhine and platyrrhine primate herpesviruses, little is known about herpesviruses in prosimians. We aimed to identify and characterize herpesviruses in prosimians with proliferative lymphocytic disease. DNA was extracted from tissues of 9 gray mouse lemurs (Microcebus murinus) and 3 pygmy slow lorises (Nycticebus pygmaeus) with lymphoproliferative lesions, and we performed nested PCR and sequencing for detection of herpesviruses and polyomaviruses. We identified 3 novel herpesviruses and performed phylogenetic analyses to characterize their relationship with other herpesviruses. A gray mouse lemur herpesvirus clustered with other primate herpesviruses within the subfamily Betaherpesvirinae, just basal to the genus Cytomegalovirus. The other gray mouse lemur herpesvirus and the pygmy slow loris herpesvirus clustered within the subfamily Gammaherpesvirinae, although the relationships within the subfamily were less resolved. Quantitative PCR assays were developed for the 2 new gray mouse lemur viruses, providing specific, faster, less expensive, and quantitative detection tools. Further studies are needed to elucidate the relationship between the presence of these viruses and the severity or presence of lymphoproliferative lesions in prosimians.

Comparative Genomic Analysis Identifies Great-Ape-Specific Structural Variants and Their Evolutionary Relevance

During the origin of great apes about 14 million years ago, a series of phenotypic innovations emerged, such as the increased body size, the enlarged brain volume, the improved cognitive skill, and the diversified diet. Yet, the genomic basis of these evolutionary changes remains unclear. Utilizing the high-quality genome assemblies of great apes (including human), gibbon, and macaque, we conducted comparative genome analyses and identified 15,885 great ape-specific structural variants (GSSVs), including eight coding GSSVs resulting in the creation of novel proteins (e.g., ACAN and CMYA5). Functional annotations of the GSSV-related genes revealed the enrichment of genes involved in development and morphogenesis, especially neurogenesis and neural network formation, suggesting the potential role of GSSVs in shaping the great ape-shared traits. Further dissection of the brain-related GSSVs shows great ape-specific changes of enhancer activities and gene expression in the brain, involving a group of GSSV-regulated genes (such as NOL3) that potentially contribute to the altered brain development and function in great apes. The presented data highlight the evolutionary role of structural variants in the phenotypic innovations during the origin of the great ape lineage.

Multiple routes of communication within the amygdala-mPFC network: A comparative approach in humans and macaques

The network formed by the amygdala (AMG) and the medial Prefrontal Cortex (mPFC), at the interface between our internal and external environment, has been shown to support some important aspects of behavioral adaptation. Whether and how the anatomo-functional organization of this network evolved across primates remains unclear. Here, we compared AMG nuclei morphological characteristics and their functional connectivity with the mPFC in humans and macaques to identify potential homologies and differences between these species. Based on selected studies, we highlight two subsystems within the AMG-mPFC circuits, likely involved in distinct temporal dynamics of integration during behavioral adaptation. We also show that whereas the mPFC displays a large expansion but a preserved intrinsic anatomo-functional organization, the AMG displays a volume reduction and morphological changes related to specific nuclei. We discuss potential commonalities and differences in the dialogue between AMG nuclei and mPFC in humans and macaques based on available data.

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.

Comprehension of own and other species' alarm calls in sooty mangabey vocal development

Abstract

Primates understand the meaning of their own and other species' alarm calls, but little is known about how they acquire such knowledge. Here, we combined direct behavioural observations with playback experiments to investigate two key processes underlying vocal development: comprehension and usage. Especifically, we studied the development of con- and heterospecific alarm call recognition in free-ranging sooty mangabeys, Cercocebus atys, across three age groups: young juveniles (1-2y), old juveniles (3-4y) and adults (> 5y). We observed that, during natural predator encounters, juveniles alarm called to a significantly wider range of species than adults, with evidence of refinement during the first four years of life. In the experiments, we exposed subjects to leopard, eagle and snake alarm calls given by other group members or sympatric Diana monkeys. We found that young juveniles' locomotor and vocal responses were least appropriate and that they engaged in more social referencing (look at adults when hearing an alarm call) than older individuals, suggesting that vocal competence is obtained via social learning. In conclusion, our results suggest that alarm call comprehension is socially learned during the juvenile stage, with comprehension preceding appropriate usage but no difference between learning their own or other species' alarm calls.

Significance statement

Under natural conditions, animals do not just interact with members of their own species, but usually operate in a network of associated species. However, ontogenetic research on primate communication frequently ignores this significant element. We studied the development of con- and heterospecific alarm call recognition in wild sooty mangabeys. We found that communicative competence was acquired during the juvenile stages, with alarm call comprehension learning preceding appropriate vocal usage and with no clear difference in learning of con- and heterospecific signals. We also found that, during early stages of life, social referencing, a proactive form of social learning, was key in the acquisition of competent alarm call behaviour. Our results show that primates equally learn to interpret alarm calls from their own and other species during their early stages of life and that this learning process is refined as the animals mature.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00265-023-03318-6.

A revised perspective on the evolution of the lateral frontal cortex in primates

Detailed neuroscientific data from macaque monkeys have been essential in advancing understanding of human frontal cortex function, particularly for regions of frontal cortex without homologs in other model species. However, precise transfer of this knowledge for direct use in human applications requires an understanding of monkey to hominid homologies, particularly whether and how sulci and cytoarchitectonic regions in the frontal cortex of macaques relate to those in hominids. We combine sulcal pattern analysis with resting-state functional magnetic resonance imaging and cytoarchitectonic analysis to show that old-world monkey brains have the same principles of organization as hominid brains, with the notable exception of sulci in the frontopolar cortex. This essential comparative framework provides insights into primate brain evolution and a key tool to drive translation from invasive research in monkeys to human applications.

Call combinations and compositional processing in wild chimpanzees

Through syntax, i.e., the combination of words into larger phrases, language can express a limitless number of messages. Data in great apes, our closest-living relatives, are central to the reconstruction of syntax's phylogenetic origins, yet are currently lacking. Here, we provide evidence for syntactic-like structuring in chimpanzee communication. Chimpanzees produce "alarm-huus" when surprised and "waa-barks" when potentially recruiting conspecifics during aggression or hunting. Anecdotal data suggested chimpanzees combine these calls specifically when encountering snakes. Using snake presentations, we confirm call combinations are produced when individuals encounter snakes and find that more individuals join the caller after hearing the combination. To test the meaning-bearing nature of the call combination, we use playbacks of artificially-constructed call combinations and both independent calls. Chimpanzees react most strongly to call combinations, showing longer looking responses, compared with both independent calls. We propose the "alarm-huu + waa-bark" represents a compositional syntactic-like structure, where the meaning of the call combination is derived from the meaning of its parts. Our work suggests that compositional structures may not have evolved de novo in the human lineage, but that the cognitive building-blocks facilitating syntax may have been present in our last common ancestor with chimpanzees.

NGSEP 4: Efficient and accurate identification of orthogroups and whole-genome alignment

Whole-genome alignment allows researchers to understand the genomic structure and variation among genomes. Approaches based on direct pairwise comparisons of DNA sequences require large computational capacities. As a consequence, pipelines combining tools for orthologous gene identification and synteny have been developed. In this manuscript, we present the latest functionalities implemented in NGSEP 4, to identify orthogroups and perform whole genome alignments. NGSEP implements functionalities for identification of clusters of homologus genes, synteny analysis and whole genome alignment. Our results showed that the NGSEP algorithm for orthogroups identification has competitive accuracy and efficiency in comparison to commonly used tools. The implementation also includes a visualization of the whole genome alignment based on synteny of the orthogroups that were identified, and a reconstruction of the pangenome based on frequencies of the orthogroups among the genomes. NGSEP 4 also includes a new graphical user interface based on the JavaFX technology. We expect that these new developments will be very useful for several studies in evolutionary biology and population genomics.

Non-structural genes of novel lemur adenoviruses reveal codivergence of virus and host

Adenoviruses (AdVs) are important human and animal pathogens and are frequently used as vectors for gene therapy and vaccine delivery. Surprisingly, there are only scant data regarding primate AdV origin and evolution, especially in the most basal primate hosts. We detect and sequence AdVs from faeces of two Madagascan lemur species. Complete genome sequence analyses define a new AdV species with a particularly large gene encoding a protein of unknown function in the early gene region 3. Unexpectedly, the new AdV species is not most similar to human or other simian AdVs but to bat adenovirus C. Genome characterisation shows signals of virus-host codivergence in non-structural genes, which show lower diversity than structural genes. Outside a lemur species mixing zone, recombination less frequently separates structural genes, as in human adenovirus C. The evolutionary history of lemur AdVs likely involves both a host switch and codivergence with the lemur hosts.

Molecular and cellular mechanisms of human cortical connectivity

Comparative studies of the cerebral cortex have identified various human and primate-specific changes in both local and long-range connectivity, which are thought to underlie our advanced cognitive capabilities. These changes are likely mediated by the divergence of spatiotemporal regulation of gene expression, which is particularly prominent in the prenatal and early postnatal human and non-human primate cerebral cortex. In this review, we describe recent advances in characterizing human and primate genetic and cellular innovations including identification of novel species-specific, especially human-specific, genes, gene expression patterns, and cell types. Finally, we highlight three recent studies linking these molecular changes to reorganization of cortical connectivity.

Structurally divergent and recurrently mutated regions of primate genomes

To better understand the pattern of primate genome structural variation, we sequenced and assembled using multiple long-read sequencing technologies the genomes of eight nonhuman primate species, including New World monkeys (owl monkey and marmoset), Old World monkey (macaque), Asian apes (orangutan and gibbon), and African ape lineages (gorilla, bonobo, and chimpanzee). Compared to the human genome, we identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. Across 50 million years of primate evolution, we estimate that 819.47 Mbp or ~27% of the genome has been affected by SVs based on analysis of these primate lineages. We identify 1,607 structurally divergent regions (SDRs) wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (CARDs, ABCD7, OLAH) and new lineage-specific genes are generated (e.g., CKAP2, NEK5) and have become targets of rapid chromosomal diversification and positive selection (e.g., RGPDs). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species for the first time.

Nonhuman primate genetic models for the study of rare diseases

Pre-clinical research and development relies heavily upon translationally valid models of disease. A major difficulty in understanding the biology of, and developing treatments for, rare disease is the lack of animal models. It is important that these models not only recapitulate the presentation of the disease in humans, but also that they share functionally equivalent underlying genetic causes. Nonhuman primates share physiological, anatomical, and behavioral similarities with humans resulting from close evolutionary relationships and high genetic homology. As the post-genomic era develops and next generation sequencing allows for the resequencing and screening of large populations of research animals, naturally occurring genetic variation in nonhuman primates with clinically relevant phenotypes is regularly emerging. Here we review nonhuman primate models of multiple rare genetic diseases with a focus on the similarities and differences in manifestation and etiologies across species. We discuss how these models are being developed and how they can offer new tools and opportunities for researchers interested in exploring novel therapeutics for these and other genetic diseases. Modeling human genetic diseases in translationally relevant nonhuman primates presents new prospects for development of therapeutics and a better understanding of rare diseases. The post-genomic era offers the opportunity for the discovery and further development of more models like those discussed here.

NANOGP1, a tandem duplicate of NANOG, exhibits partial functional conservation in human naïve pluripotent stem cells

Gene duplication events can drive evolution by providing genetic material for new gene functions, and they create opportunities for diverse developmental strategies to emerge between species. To study the contribution of duplicated genes to human early development, we examined the evolution and function of NANOGP1, a tandem duplicate of the transcription factor NANOG. We found that NANOGP1 and NANOG have overlapping but distinct expression profiles, with high NANOGP1 expression restricted to early epiblast cells and naïve-state pluripotent stem cells. Sequence analysis and epitope-tagging revealed that NANOGP1 is protein coding with an intact homeobox domain. The duplication that created NANOGP1 occurred earlier in primate evolution than previously thought and has been retained only in great apes, whereas Old World monkeys have disabled the gene in different ways, including homeodomain point mutations. NANOGP1 is a strong inducer of naïve pluripotency; however, unlike NANOG, it is not required to maintain the undifferentiated status of human naïve pluripotent cells. By retaining expression, sequence and partial functional conservation with its ancestral copy, NANOGP1 exemplifies how gene duplication and subfunctionalisation can contribute to transcription factor activity in human pluripotency and development.

CD-MAWS: An Alignment-Free Phylogeny Estimation Method Using Cosine Distance on Minimal Absent Word Sets

Multiple sequence alignment has been the traditional and well established approach of sequence analysis and comparison, though it is time and memory consuming. As the scale of sequencing data is increasing day by day, the importance of faster yet accurate alignment-free methods is on the rise. Several alignment-free sequence analysis methods have been established in the literature in recent years, which extract numerical features from genomic data to analyze sequences and also to estimate phylogenetic relationship among genes and species. Minimal Absent Word (MAW) is an effective concept for representing characteristics of a sequence in an alignment-free manner. In this study, we present CD-MAWS, a distance measure based on cosine of the angle between composition vectors constructed using minimal absent words, for sequence analysis in a computationally inexpensive manner. We have benchmarked CD-MAWS using several AFProject datasets, such as Fish mtDNA, E.coli, Plants, Shigella and Yersinia datasets, and found it to perform quite well. Applied on several other biological datasets such as mammal mtDNA, bacterial genomes and viral genomes, CD-MAWS resolved phylogenetic relationships similar to or better than state-of-the-art alignment-free methods such as Mash, Skmer, Co-phylog and kSNP3.

Evolutionary analysis of endogenous intronic retroviruses in primates reveals an enrichment in transcription binding sites associated with key regulatory processes

Background

Endogenous retroviruses (ERVs) are the result of the integration of retroviruses into host DNA following germline infection. Endogenous retroviruses are made up of three main genes: gag, pol, and env, each of which encodes viral proteins that can be conserved or not. ERVs have been observed in a wide range of vertebrate genomes and their functions are associated with viral silencing and gene regulation.

Results

In this work, we studied the evolutionary history of endogenous retroviruses associated with five human genes (INPP5B, DET1, PSMA1, USH2A, and MACROD2), which are located within intron sections. To verify the retroviral origin of the candidates, several approaches were used to detect and locate ERV elements. Both orthologous and paralogous genes were identified by Ensembl and then analyzed for ERV presence using RetroTector. A phylogenetic tree was reconstructed to identify the minimum time point of ERV acquisition. From that search, we detected ERVs throughout the primate lineage and in some other groups. Also, we identified the minimum origin of the ERVs from the parvorder Catarrhini to the Homininae subfamily.

Conclusions

With the data collected, and by observing the transcription factors annotated inside ERVs, we propose that these elements play a relevant role in gene expression regulation and they probably possess important features for tumorigenesis control.

Ancestral APOBEC3B Nuclear Localization Is Maintained in Humans and Apes and Altered in Most Other Old World Primate Species

APOBEC3B is an innate immune effector enzyme capable of introducing mutations in viral genomes through DNA cytosine-to-uracil editing. Recent studies have shown that gamma-herpesviruses, such as Epstein-Barr virus (EBV), have evolved a potent APOBEC3B neutralization mechanism to protect lytic viral DNA replication intermediates in the nuclear compartment. APOBEC3B is additionally unique as the only human DNA deaminase family member that is constitutively nuclear. Nuclear localization has therefore been inferred to be essential for innate antiviral function. Here, we combine evolutionary, molecular, and cell biology approaches to address whether nuclear localization is a conserved feature of APOBEC3B in primates. Despite the relatively recent emergence of APOBEC3B approximately 30 to 40 million years ago (MYA) in Old World primates by genetic recombination (after the split from the New World monkey lineage 40 to 50 MYA), we find that the hallmark nuclear localization of APOBEC3B shows variability. For instance, although human and several nonhuman primate APOBEC3B enzymes are predominantly nuclear, rhesus macaque and other Old World primate APOBEC3B proteins are clearly cytoplasmic or cell wide. A series of human/rhesus macaque chimeras and mutants combined to map localization determinants to the N-terminal half of the protein with residues 15, 19, and 24 proving critical. Ancestral APOBEC3B reconstructed from present-day primate species also shows strong nuclear localization. Together, these results indicate that the ancestral nuclear localization of APOBEC3B is maintained in present-day human and ape proteins, but nuclear localization is not conserved in all Old World monkey species despite a need for antiviral functions in the nuclear compartment. IMPORTANCE APOBEC3 enzymes are single-stranded DNA cytosine-to-uracil deaminases with beneficial roles in antiviral immunity and detrimental roles in cancer mutagenesis. Regarding viral infection, all seven human APOBEC3 enzymes have overlapping roles in restricting virus types that require DNA for replication, including EBV, HIV, human papillomavirus (HPV), and human T-cell leukemia virus (HTLV). Regarding cancer, at least two APOBEC3 enzymes, APOBEC3B and APOBEC3A, are prominent sources of mutation capable of influencing clinical outcomes. Here, we combine evolutionary, molecular, and cell biology approaches to characterize primate APOBEC3B enzymes. We show that nuclear localization is an ancestral property of APOBEC3B that is maintained in present-day human and ape enzymes, but not conserved in other nonhuman primates. This partial mechanistic conservation indicates that APOBEC3B is important for limiting the replication of DNA-based viruses in the nuclear compartment. Understanding these pathogen-host interactions may contribute to the development of future antiviral and antitumor therapies.

Phylogeography and evolutionary lineage diversity in the small-eared greater galago, Otolemur garnettii (Primates: Galagidae)

Assessing the true lineage diversity in elusive nocturnal organisms is particularly challenging due to their subtle phenotypic variation in diagnostic traits. The cryptic small-eared greater galago (Otolemur garnettii) offers a great opportunity to test if currently recognized subspecies, suggested by discontinuities in coat colour pattern and geographic barriers, represent distinct evolutionary lineages. To answer this question, we conducted the first population-level phylogeographic study of the species, sampling wild specimens from across almost its entire latitudinal range, including the Zanzibar Archipelago. We applied five species-delimitation algorithms to investigate the genetic diversity and distribution pattern of mitochondrial DNA across the geographic range of three out of four subspecies. Our results suggest that far-northern populations of O. g. lasiotis potentially represent an independently evolving lineage, but populations assigned to O. g. garnettii from Zanzibar Island and of O. g panganiensis from mainland Tanzania do not constitute two independent lineages. A dated phylogeny suggests that this northern clade diverged from all remaining samples approximately 4 Mya. Such old divergence age is in line with the split between many galagid species. This northern lineage could potentially represent an incipient species; however, there is not yet enough evidence to support a new taxonomic status for this unique mitochondrial group.

The earliest hylobatid from the Late Miocene of China

Yuanmoupithecus xiaoyuan, a small catarrhine from the Late Miocene of Yunnan in southern China, was initially suggested to be related to Miocene proconsuloids or dendropithecoids from East Africa, but subsequent reports indicated that it might be more closely related to hylobatids. Here, detailed comparisons of the material, including seven newly discovered teeth and a partial lower face of a juvenile individual, provide crucial evidence to help establish its phylogenetic relationships. Yuanmoupithecus exhibits a suite of synapomorphies that support a close phylogenetic relationship with extant hylobatids. Furthermore, based on the retention of several primitive features of the dentition, Yuanmoupithecus can be shown to be the sister taxon of crown hylobatids. The contention that Kapi ramnagarensis from the Middle Miocene of India might represent an earlier species of hylobatid is not supported here. Instead, Kapi is inferred to be a specialized pliopithecoid more closely related to Krishnapithecus krishnaii from the Late Miocene of India. Currently then, Yuanmoupithecus represents the earliest known definitively identified hylobatid and the only member of the clade predating the Pleistocene. It extends the fossil record of hylobatids back to 7-8 Ma and fills a critical gap in the evolutionary history of hominoids that has up until now remained elusive. Even so, molecular estimates of a divergence date of hylobatids from other hominoids at about 17-22 Ma signifies that there is still a substantial gap in the fossil record of more than 10 million years that needs to be filled in order to document the biogeographic origins and early evolution of hylobatids.

Genomic signatures of high-altitude adaptation and chromosomal polymorphism in geladas

Primates have adapted to numerous environments and lifestyles, but very few species are native to high elevations. Here, we investigated high-altitude adaptations in the gelada (Theropithecus gelada), a monkey endemic to the Ethiopian Plateau. We examined genome-wide variation in conjunction with measurements of hematological and morphological traits. Our new gelada reference genome is highly intact and assembled at chromosome-length levels. Unexpectedly, we identified a chromosomal polymorphism in geladas that could potentially contribute to reproductive barriers between populations. Compared to baboons at low altitude, we found that high-altitude geladas exhibit significantly expanded chest circumferences, potentially allowing for greater lung surface area for increased oxygen diffusion. We identified gelada-specific amino acid substitutions in the alpha-chain subunit of adult hemoglobin but found that gelada hemoglobin does not exhibit markedly altered oxygenation properties compared to lowland primates. We also found that geladas at high altitude do not exhibit elevated blood hemoglobin concentrations, in contrast to the normal acclimatization response to hypoxia in lowland primates. The absence of altitude-related polycythemia suggests that geladas are able to sustain adequate tissue-oxygen delivery despite environmental hypoxia. Finally, we identified numerous genes and genomic regions exhibiting accelerated rates of evolution, as well as gene families exhibiting expansions in the gelada lineage, potentially reflecting altitude-related selection. Our findings lend insight into putative mechanisms of high-altitude adaptation while suggesting promising avenues for functional hypoxia 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.

The Island of Female Power? Intersexual Dominance Relationships in the Lemurs of Madagascar

The extant primates of Madagascar (Lemuriformes) represent the endpoints of an adaptive radiation following a single colonization event more than 50 million years ago. They have since evolved a diversity of life history traits, ecological adaptations and social systems that rivals that of all other living primates combined. Their social systems are characterized by a unique combination of traits, including the ability of adult females to dominate adult males. In fact, there is no other group of mammals in which female dominance is so widespread. Yet, recent research has indicated that there is more interspecific variation in lemur intersexual relationships than previously acknowledged. Here, we therefore review and summarize the relevant literature, quantifying the extent of sex-bias in intersexual dominance relations documented in observational and experimental studies in captivity and the wild. Female dominance is often, but not always, implemented by spontaneous male submission in the absence of female aggression and linked to female sexual maturation. We connect the available evidence to the hypotheses that have been proposed to explain the evolution of female dominance among lemurs. The occurrence of female dominance in all lemur families and the interspecific variation in its extent indicate that it has evolved soon after lemurs colonized Madagascar – presumably in response to particular ecological challenges – and that it has since been reduced in magnitude independently in some taxa. Our study contributes important comparative information on sex roles from an independent primate radiation and provides general insights into the conditions, opportunities and obstacles in the evolution of female-biased power.

Molecular characterization of matrix metalloproteinase gene family across primates

Deregulation of matrix metalloproteinases (MMPs) contributes considerably to cancers, psychiatric disorders, macular degeneration and bone diseases. The use of humans in the development of MMPs as prognostic biomarkers and therapeutic targets is complicated by many factors, while primate models can be useful alternatives for this purpose. Here, we performed genome-enabled identification of putative MMPs across primate species, and comprehensively investigated the genes. Phylogenetic topology of the MMP family showed each type formulates a distinct clade, and was further clustered to classes, largely agreeing with classification based on biochemical properties and domain organization. Across primates, the excess of candidate sites of positive selection was detected for MMP-19, in addition to 1-3 sites in MMP-8, MMP-10 and MMP-26. MMP-26 showed Ka/Ks value above 1 between human and chimpanzee copies. We observed two copies of MMP-19 in the old-world monkey genomes, suggesting gene duplication at the early stage of or prior to the emergence of the lineage. Furin-activatable MMPs demonstrate the most variable properties regarding Domain organization and gene structure. During human aging, MMP-11 showed gradually decreased expression in testis, so as MMP-2, MMP-14, MMP15 and MMP-28 in ovary, while MMP-7 and MMP-21 showed elevated expression, implying their distinct roles in different reproductive organs. Co-expression clusters were formed among human MMPs both within and across classes, and expression correlation was observed in MMP genes across primates. Our results illuminate the utilization of MMPs for the discovery of prognostic biomarkers and therapeutic targets for aging-related diseases and carry new messages on MMP classification.

Recurrent evolution of an inhibitor of ESCRT-dependent virus budding and LINE-1 retrotransposition in primates

Most antiviral proteins recognize specific features of viruses. In contrast, the recently described antiviral factor retroCHMP3 interferes with the "host endosomal complexes required for transport" (ESCRT) pathway to inhibit the budding of enveloped viruses. RetroCHMP3 arose independently on multiple occasions via duplication and truncation of the gene encoding the ESCRT-III factor CHMP3. However, since the ESCRT pathway is essential for cellular membrane fission reactions, ESCRT inhibition is potentially cytotoxic. This raises fundamental questions about how hosts can repurpose core cellular functions into antiviral functions without incurring a fitness cost due to excess cellular toxicity. We reveal the evolutionary process of detoxification for retroCHMP3 in New World monkeys using a combination of ancestral reconstructions, cytotoxicity, and virus release assays. A duplicated, full-length copy of retroCHMP3 in the ancestors of New World monkeys provides modest inhibition of virus budding while exhibiting subtle cytotoxicity. Ancient retroCHMP3 then accumulated mutations that reduced cytotoxicity but preserved virus inhibition before a truncating stop codon arose in the more recent ancestors of squirrel monkeys, resulting in potent inhibition. In species where full-length copies of retroCHMP3 still exist, their artificial truncation generated potent virus-budding inhibitors with little cytotoxicity, revealing the potential for future antiviral defenses in modern species. In addition, we discovered that retroCHMP3 restricts LINE-1 retrotransposition, revealing how different challenges to genome integrity might explain multiple independent origins of retroCHMP3 in different species to converge on new immune functions.

A Comprehensive Atlas of Immunological Differences Between Humans, Mice, and Non-Human Primates

Animal models are an integral part of the drug development and evaluation process. However, they are unsurprisingly imperfect reflections of humans, and the extent and nature of many immunological differences are unknown. With the rise of targeted and biological therapeutics, it is increasingly important that we understand the molecular differences in the immunological behavior of humans and model organisms. However, very few antibodies are raised against non-human primate antigens, and databases of cross-reactivity between species are incomplete. Thus, we screened 332 antibodies in five immune cell populations in blood from humans and four non-human primate species generating a comprehensive cross-reactivity catalog that includes cell type-specificity. We used this catalog to create large mass cytometry universal cross-species phenotyping and signaling panels for humans, along with three of the model organisms most similar to humans: rhesus and cynomolgus macaques and African green monkeys; and one of the mammalian models most widely used in drug development: C57BL/6 mice. As a proof-of-principle, we measured immune cell signaling responses across all five species to an array of 15 stimuli using mass cytometry. We found numerous instances of different cellular phenotypes and immune signaling events occurring within and between species, and detailed three examples (double-positive T cell frequency and signaling; granulocyte response to Bacillus anthracis antigen; and B cell subsets). We also explore the correlation of herpes simian B virus serostatus on the immune profile. Antibody panels and the full dataset generated are available online as a resource to enable future studies comparing immune responses across species during the evaluation of therapeutics.

Identification of multiple TAR DNA binding protein retropseudogene lineages during the evolution of primates

The TAR DNA Binding Protein (TARDBP) gene has become relevant after the discovery of its several pathogenic mutations. The lack of evolutionary history is in contrast to the amount of studies found in the literature. This study investigated the evolutionary dynamics associated with the retrotransposition of the TARDBP gene in primates. We identified novel retropseudogenes that likely originated in the ancestors of anthropoids, catarrhines, and lemuriformes, i.e. the strepsirrhine clade that inhabit Madagascar. We also found species-specific retropseudogenes in the Philippine tarsier, Bolivian squirrel monkey, capuchin monkey and vervet. The identification of a retropseudocopy of the TARDBP gene overlapping a lncRNA that is potentially expressed opens a new avenue to investigate TARDBP gene regulation, especially in the context of TARDBP associated pathologies.

Comparative morphology and physiology of the vocal production apparatus and the brain in the extant primates

Objective data mainly from the comparative anatomy of various organs related to human speech and language is considered to unearth clues about the mechanisms behind language development. The two organs of the larynx and hyoid bone are considered to have evolved towards suitable positions and forms in preparation for the occurrence of the large repertoire of vocalization necessary for human speech. However, some researchers have asserted that there is no significant difference of these organs between humans and non-human primates. Speech production is dependent on the voluntary control of the respiratory, laryngeal, and vocal tract musculature. Such control is fully present in humans but only partially so in non-human primates, which appear to be able to voluntarily control only supralaryngeal articulators. Both humans and non-human primates have direct cortical innervation of motor neurons controlling the supralaryngeal vocal tract but only human appear to have direct cortical innervation of motor neurons controlling the larynx. In this review, we investigate the comparative morphology and function of the wide range of components involved in vocal production, including the larynx, the hyoid bone, the tongue, and the vocal brain. We would like to emphasize the importance of the tongue in the primary development of human speech and language. It is now time to reconsider the possibility of the tongue playing a definitive role in the emergence of human speech.

Monkeying around with venom: an increased resistance to α-neurotoxins supports an evolutionary arms race between Afro-Asian primates and sympatric cobras

Background

Snakes and primates have a multi-layered coevolutionary history as predators, prey, and competitors with each other. Previous work has explored the Snake Detection Theory (SDT), which focuses on the role of snakes as predators of primates and argues that snakes have exerted a selection pressure for the origin of primates’ visual systems, a trait that sets primates apart from other mammals. However, primates also attack and kill snakes and so snakes must simultaneously avoid primates. This factor has been recently highlighted in regard to the movement of hominins into new geographic ranges potentially exerting a selection pressure leading to the evolution of spitting in cobras on three independent occasions.

Results

Here, we provide further evidence of coevolution between primates and snakes, whereby through frequent encounters and reciprocal antagonism with large, diurnally active neurotoxic elapid snakes, Afro-Asian primates have evolved an increased resistance to α-neurotoxins, which are toxins that target the nicotinic acetylcholine receptors. In contrast, such resistance is not found in Lemuriformes in Madagascar, where venomous snakes are absent, or in Platyrrhini in the Americas, where encounters with neurotoxic elapids are unlikely since they are relatively small, fossorial, and nocturnal. Within the Afro-Asian primates, the increased resistance toward the neurotoxins was significantly amplified in the last common ancestor of chimpanzees, gorillas, and humans (clade Homininae). Comparative testing of venoms from Afro-Asian and American elapid snakes revealed an increase in α-neurotoxin resistance across Afro-Asian primates, which was likely selected against cobra venoms. Through structure-activity studies using native and mutant mimotopes of the α-1 nAChR receptor orthosteric site (loop C), we identified the specific amino acids responsible for conferring this increased level of resistance in hominine primates to the α-neurotoxins in cobra venom.

Conclusion

We have discovered a pattern of primate susceptibility toward α-neurotoxins that supports the theory of a reciprocal coevolutionary arms-race between venomous snakes and primates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01195-x.

Comparative anatomy of the carotid canal in the Miocene small-bodied catarrhine Pliobates cataloniae

The small-bodied Miocene catarrhine Pliobates cataloniae (11.6 Ma, Spain) displays a mosaic of catarrhine symplesiomorphies and hominoid synapomorphies that hinders deciphering its phylogenetic relationships. Based on cladistic analyses, it has been interpreted as a stem hominoid or as a pliopithecoid. Intriguingly, the carotid canal orientation of Pliobates was originally described as hylobatid-like. The variation in carotid canal morphology among anthropoid clades shown in previous studies suggests that this structure might be phylogenetically informative. However, its potential for phylogenetic reconstruction among extinct catarrhines remains largely unexplored. Here we quantify the orientation, proportions, and course of the carotid canal in Pliobates, extant anthropoids and other Miocene catarrhines (Epipliopithecus, Victoriapithecus, and Ekembo) using three-dimensional morphometric techniques. We also compute phylogenetic signal and reconstruct the ancestral carotid canal course for main anthropoid clades. Our results reveal that carotid canal morphology embeds strong phylogenetic signal but mostly discriminates between platyrrhines and catarrhines, with an extensive overlap among extant catarrhine families. The analyzed extinct taxa display a quite similar carotid canal morphology more closely resembling that of extant catarrhines. Nevertheless, our results for Pliobates highlight some differences compared with the pliopithecid Epipliopithecus, which displays a somewhat more platyrrhine-like morphology. In contrast, Pliobates appears as derived toward the modern catarrhine condition as the stem cercopithecid Victoriapithecus and the stem hominoid Ekembo, which more closely resemble one another. Moreover, Pliobates appears somewhat derived toward the reconstructed ancestral hominoid morphotype, being more similar than other Miocene catarrhines to the condition of great apes and the hylobatid Symphalangus. Overall, our results rule out previously noted similarities in carotid canal morphology between Pliobates and hylobatids, but do not show particular similarities with pliopithecoids either-as opposed to extant and other extinct catarrhines. Additional analyses will be required to clarify the phylogenetic relationships of Pliobates, particularly given its dental similarities with dendropithecids.

Ten millennia of hepatitis B virus evolution

[Figure: see text].

Evolutionary History of Alpha Satellite DNA Repeats Dispersed within Human Genome Euchromatin

Major human alpha satellite DNA repeats are preferentially assembled within (peri)centromeric regions but are also dispersed within euchromatin in the form of clustered or short single repeat arrays. To study the evolutionary history of single euchromatic human alpha satellite repeats (ARs), we analyzed their orthologous loci across the primate genomes. The continuous insertion of euchromatic ARs throughout the evolutionary history of primates starting with the ancestors of Simiformes (45-60 Ma) and continuing up to the ancestors of Homo is revealed. Once inserted, the euchromatic ARs were stably transmitted to the descendant species, some exhibiting copy number variation, whereas their sequence divergence followed the species phylogeny. Many euchromatic ARs have sequence characteristics of (peri)centromeric alpha repeats suggesting heterochromatin as a source of dispersed euchromatic ARs. The majority of euchromatic ARs are inserted in the vicinity of other repetitive elements such as L1, Alu, and ERV or are embedded within them. Irrespective of the insertion context, each AR insertion seems to be unique and once inserted, ARs do not seem to be subsequently spread to new genomic locations. In spite of association with (retro)transposable elements, there is no indication that such elements play a role in ARs proliferation. The presence of short duplications at most of ARs insertion sites suggests site-directed recombination between homologous motifs in ARs and in the target genomic sequence, probably mediated by extrachromosomal circular DNA, as a mechanism of spreading within euchromatin.

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.

Phylogenomics of African radiation of Praomyini (Muridae: Murinae) rodents: First fully resolved phylogeny, evolutionary history and delimitation of extant genera

The tribe Praomyini is a diversified group including 64 species and eight extant rodent genera. They live in a broad spectrum of habitats across whole sub-Saharan Africa. Members of this tribe are often very abundant, they have a key ecological role in ecosystems, they are hosts of many potentially pathogenic microorganisms and comprise numerous agricultural pests. Although this tribe is well supported by both molecular and morphological data, its intergeneric relationships and the species contents of several genera are not yet fully resolved. Recent molecular data suggest that at least three genera in current sense are paraphyletic. However, in these studies the species sampling was sparse and the resolution of relationships among genera was poor, probably due to a fast radiation of the tribe dated to the Miocene and insufficient amount of genetic data. Here we used genomic scale data (395 nuclear loci = 610,965 bp long alignment and mitogenomes = 14,745 bp) and produced the first fully resolved species tree containing most major lineages of the Praomyini tribe (i.e. all but one currently delimited genera and major intrageneric clades). Results of a fossil-based divergence dating analysis suggest that the radiation started during the Messinian stage (ca. 7 Ma) and was likely linked to a fragmentation of the pan-African Miocene forest. Some lineages remained in the rain forests, while many others adapted to a broad spectrum of new open lowland and montane habitats that appeared at the beginning of Pliocene. Our analyses clearly confirmed the presence of three polyphyletic genera (Praomys, Myomyscus and Mastomys). We review current knowledge of these three genera and suggest corresponding taxonomic changes. To keep genera monophyletic, we propose taxonomic re-arrangements and delimit four new genera. Furthermore, we discovered a new highly divergent genetic lineage of Praomyini in southwestern Ethiopia, which is described as a new species and genus.

Evolutionary history of orangutan plasmodia revealed by phylogenetic analysis of complete mtDNA genomes and new biogeographical divergence dating calibration models

During the past 15 years, researchers have shown a renewed interest in the study of the Plasmodium parasites that infect orangutans. Most recently, studies examined the phylogenetic relationships and divergence dates of these parasites in orangutans using complete mitochondrial DNA genomes. Questions regarding the dating of these parasites, however, remain. In the present study, we provide a new calibration model for dating the origins of Plasmodium parasites in orangutans using a modified date range for the origin of macaques in Asia. Our Bayesian phylogenetic analyses of complete Plasmodium sp. mitochondrial DNA genomes inferred two clades of plasmodia in orangutans (Pongo 1 and Pongo 2), and that these clades likely represent the previously identified species Plasmodium pitheci and Plasmodium silvaticum. However, we cannot identify which Pongo clade is representative of the morphologically described species. The most recent common ancestor of both Pongo sp. plasmodia, Plasmodium. hylobati, and Plasmodium. inui dates to 3-3.16 million years ago (mya) (95% highest posterior density [HPD]: 2.09-4.08 mya). The Pongo 1 parasite diversified 0.33-0.36 mya (95% HPD: 0.12-0.63), while the Pongo 2 parasite diversified 1.15-1.22 mya (95% HPD: 0.63-1.82 mya). It now seems likely that the monkey Plasmodium (P. inui) is the result of a host switch event from the Pongo 2 parasite to sympatric monkeys, or P. hylobati. Our new estimates for the divergence of orangutan malaria parasites, and subsequent diversification, are all several hundred thousand years later than previous Bayesian estimates.

The germline mutational process in rhesus macaque and its implications for phylogenetic dating

BACKGROUND

Understanding the rate and pattern of germline mutations is of fundamental importance for understanding evolutionary processes.

RESULTS

Here we analyzed 19 parent-offspring trios of rhesus macaques (Macaca mulatta) at high sequencing coverage of ∼76× per individual and estimated a mean rate of 0.77 × 10-8de novo mutations per site per generation (95% CI: 0.69 × 10-8 to 0.85 × 10-8). By phasing 50% of the mutations to parental origins, we found that the mutation rate is positively correlated with the paternal age. The paternal lineage contributed a mean of 81% of the de novo mutations, with a trend of an increasing male contribution for older fathers. Approximately 3.5% of de novo mutations were shared between siblings, with no parental bias, suggesting that they arose from early development (postzygotic) stages. Finally, the divergence times between closely related primates calculated on the basis of the yearly mutation rate of rhesus macaque generally reconcile with divergence estimated with molecular clock methods, except for the Cercopithecoidea/Hominoidea molecular divergence dated at 58 Mya using our new estimate of the yearly mutation rate.

CONCLUSIONS

When compared to the traditional molecular clock methods, new estimated rates from pedigree samples can provide insights into the evolution of well-studied groups such as primates.

Glutaredoxin: Discovery, redox defense and much more

Glutaredoxin, Grx, is a small protein containing an active site cysteine pair and was discovered in 1976 by Arne Holmgren. The Grx system, comprised of Grx, glutathione, glutathione reductase, and NADPH, was first described as an electron donor for Ribonucleotide Reductase but, from the first discovery in E.coli, the Grx family has impressively grown, particularly in the last two decades. Several isoforms have been described in different organisms (from bacteria to humans) and with different functions.

The unique characteristic of Grxs is their ability to catalyse glutathione-dependent redox regulation via glutathionylation, the conjugation of glutathione to a substrate, and its reverse reaction, deglutathionylation. Grxs have also recently been enrolled in iron sulphur cluster formation. These functions have been implied in various physiological and pathological conditions, from immune defense to neurodegeneration and cancer development thus making Grx a possible drug target.

This review aims to give an overview on Grxs, starting by a phylogenetic analysis of vertebrate Grxs, followed by an analysis of the mechanisms of action, the specific characteristics of the different human isoforms and a discussion on aspects related to human physiology and diseases.

Divergence-time estimates for hominins provide insight into encephalization and body mass trends in human evolution

Quantifying speciation times during human evolution is fundamental as it provides a timescale to test for the correlation between key evolutionary transitions and extrinsic factors such as climatic or environmental change. Here, we applied a total evidence dating approach to a hominin phylogeny to estimate divergence times under different topological hypotheses. The time-scaled phylogenies were subsequently used to perform ancestral state reconstructions of body mass and phylogenetic encephalization quotient (PEQ). Our divergence-time estimates are consistent with other recent studies that analysed extant species. We show that the origin of the genus Homo probably occurred between 4.30 and 2.56 million years ago. The ancestral state reconstructions show a general trend towards a smaller body mass before the emergence of Homo, followed by a trend towards a greater body mass. PEQ estimations display a general trend of gradual but accelerating encephalization evolution. The obtained results provide a rigorous temporal framework for human evolution.