Primate Lineage Research Articles

Heritability in the Rhesus Macaque (Macaca mulatta) Vertebral Column.

The vertebral column plays a central role in primate locomotion and positional behavior. Understanding its evolution, therefore, has the potential to clarify evolutionary processes that have occurred in the primate lineage as well as the specific behaviors of extinct primates. However, to understand primate vertebral anatomy, it is important to determine how much of this anatomy is heritable and how much develops as a response to environmental factors during life. We estimated heritability for vertebral counts as well as typical cervical, thoracic, and lumbar elements from 210 individuals from the pedigreed Cayo Santiago Macaca mulatta skeletal collection. We found moderate heritability of vertebral counts (h2 = 0.216-0.326), but with strong heritability of the type of variation (e.g., a tendency to meristic or homeotic change) in the vertebral count (h2 = 0.599), suggesting a possible explanation for high variability in vertebral numbers among the hominoids in particular. The moderate heritability of vertebral count also suggests that vertebral count is an unsuitable metric for estimating the ancestral state for some taxa. We found strong heritability in the morphology of cervical and upper lumbar zygapophyseal facets (h2 = 0.548-0.550) and the thoracic spinous processes (h2 = 0.609-0.761), including high heritability of the spinous process angle in the upper thoracic and upper lumbar elements (h2 = 0.649-0.752). We suggest these are related to maintaining stability in the cervical and lumbar regions, and reducing motion in the thoracic region, respectively. We propose that spinous processes may contain greater phylogenetic information, whereas transverse processes may contain greater information of function 'in life'. We also found important size effects, suggesting that size is the most heritable component of overall form and largely responsible for intertrait differences. This suggests that it is inappropriate to indiscriminately remove size effects from morphological comparisons.

Directionality theory and mortality patterns across the primate lineage

Empirical studies of aging in primates show that local selective forces rather than phylogenetic history determine the exceptional nature of human longevity (Bronikowski et al., Science 331:1325–1328, 2011). This article proposes an evolutionary rationale for this pattern of primate mortality by invoking the parameter, Life-Table Entropy, a measure of the uncertainty in the life span of a randomly chosen newborn. Life-table entropy is positively correlated with maximal life span, that is, the mean life span of a species living under favourable conditions.The logic which underlies the exceptional nature of human longevity derives from the terrestrial life-history of humans – a singularity within the primate lineage; and the concomitant ecological constraints—the hunter-gatherer, agricultural, and industrial modes of subsistence, that have defined human evolutionary history. The effect of these ecological constraints on the evolution of life span is encoded in the Entropic Principle of Longevity: life-table entropy increases in equilibrium species, populations evolving in environments with stable, renewable resources; and decreases in opportunistic species, populations subject to fluctuating resource endowments.The Entropic Principle of Longevity is a derivative of Directionality Theory, an analytic study of the evolutionary process of variation and selection based on Evolutionary Entropy, a statistical measure of the uncertainty in the age of the mother of a randomly chosen newborn. Evolutionary entropy is the organizing concept of The Entropic Principle of Evolution: Evolutionary Entropy increases in equilibrium species and decreases in opportunistic species.

CellCover Captures Neural Stem Cell Progression in Mammalian Neocortical Development.

Definition of cell classes across the tissues of living organisms is central in the analysis of growing atlases of single-cell RNA sequencing (scRNA-seq) data across biomedicine. Marker genes for cell classes are most often defined by differential expression (DE) methods that serially assess individual genes across landscapes of diverse cells. This serial approach has been extremely useful, but is limited because it ignores possible redundancy or complementarity across genes that can only be captured by analyzing multiple genes simultaneously. We aim to identify discriminating panels of genes. To efficiently explore the vast space of possible marker panels, leverage the large number of cells often sequenced, and overcome zero-inflation in scRNA-seq data, we propose viewing gene panel selection as a variation of the "minimal set-covering problem" in combinatorial optimization. We show that this new method, CellCover, captures cell-class-specific signals in the developing mouse neocortex that are distinct from those defined by DE methods. Transfer learning experiments across mouse, primate, and human data demonstrate that CellCover identifies markers of conserved cell classes in neurogenesis, as well as temporal progression in both progenitors and neurons. Exploring markers of human outer radial glia (oRG, or basal RG) across mammals, we show that transcriptomic elements of this key cell type in the expansion of the human cortex appeared in gliogenic precursors of the rodent before the full program emerged in the primate lineage. We have assembled the public datasets we use in this report at NeMO analytics where the expression of individual genes {NeMO Individual Genes} and marker gene panels can be freely explored {NeMO: Telley 3 Sets Covering Panels}, {NeMO: Telley 12 Sets Covering Panels}, and {NeMO: Sorted Brain Cell Covering Panels}. CellCover is available in {CellCover R} and {CellCover Python}.

Independent expansion, selection, and hypervariability of the TBC1D3 gene family in humans.

TBC1D3 is a primate-specific gene family that has expanded in the human lineage and has been implicated in neuronal progenitor proliferation and expansion of the frontal cortex. The gene family and its expression have been challenging to investigate because it is embedded in high-identity and highly variable segmental duplications. We sequenced and assembled the gene family using long-read sequencing data from 34 humans and 11 nonhuman primate species. Our analysis shows that this particular gene family has independently duplicated in at least five primate lineages, and the duplicated loci are enriched at sites of large-scale chromosomal rearrangements on Chromosome 17. We find that all human copy-number variation maps to two distinct clusters located at Chromosome 17q12 and that humans are highly structurally variable at this locus, differing by as many as 20 copies and ∼1 Mbp in length depending on haplotypes. We also show evidence of positive selection, as well as a significant change in the predicted human TBC1D3 protein sequence. Last, we find that, despite multiple duplications, human TBC1D3 expression is limited to a subset of copies and, most notably, from a single paralog group: TBC1D3-CDKL These observations may help explain why a gene potentially important in cortical development can be so variable in the human population.

The primate-specific presence of interferon regulatory factor-5 pseudogene 1.

Interferon regulatory factor 5 (IRF5) is a key transcription factor in inflammatory and immune responses, with its dysregulation linked to autoimmune diseases. Using bioinformatic approaches, including Basic Local Alignment Search Tool (BLAST) for sequence similarity searches, BLAST-Like Alignment Tool (BLAT) for genome-wide alignments, and several phylogenetics software, such as Multiple Alignment using Fast Fourier Transform (MAFFT), for phylogenetic analyses, we characterized the structure, origin, and evolutionary history of the human IRF5 pseudogene 1 (IRF5P1). Our analyses reveal that IRF5P1 is a chimeric processed pseudogene containing sequences derived from multiple sources, including IRF5-like sequences from disparate organisms. We find that IRF5P1 is specific to higher primates, likely originating through an ancient retroviral integration event approximately 60 million years ago. Interestingly, IRF5P1 resides within the triple QxxK/R motif-containing (TRIQK) gene, and its antisense strand is predominantly expressed as part of the TRIQK pre-messenger RNA (mRNA). Analysis of publicly available RNA-seq data suggests potential expression of antisense IRF5P1 RNA. We hypothesize that this antisense RNA may regulate IRF5 expression through complementary binding to IRF5 mRNA, with human genetic variants potentially modulating this interaction. The conservation of IRF5P1 in the primate lineage suggests its positive effects on primate evolution and innate immunity. This study highlights the importance of investigating pseudogenes and their potential regulatory roles in shaping lineage-specific immune adaptations.

Neutral evolution of snoRNA Host Gene long non-coding RNA affects cell fate control

A fundamental challenge in molecular biology is to understand how evolving genomes can acquire new functions. Actively transcribed, non-coding parts of the genome provide a potential platform for the development of new functional sequences, but their biological and evolutionary roles remain largely unexplored. Here, we show that a set of neutrally evolving long non-coding RNAs (lncRNAs) whose introns encode small nucleolar RNAs (snoRNA Host Genes, SNHGs) are highly expressed in skin and dysregulated in inflammatory conditions. Using SNHG7 and human epidermal keratinocytes as a model, we describe a mechanism by which these lncRNAs can increase self-renewal and inhibit differentiation. The activity of SNHG7 lncRNA has been recently acquired in the primate lineage and depends on a short sequence required for microRNA binding. Taken together, our results highlight the importance of understanding the role of fast-evolving transcripts in normal and diseased epithelia, and show how poorly conserved, actively transcribed non-coding sequences can participate in the evolution of genomic functionality.

A prominent vertical occipital white matter fasciculus unique to primate brains

Vision in humans and other primates enlists parallel processing streams in the dorsal and ventral visual cortex, known to support spatial and object processing, respectively. These streams are bridged, however, by a prominent white matter tract, the vertical occipital fasciculus (VOF), identified in both classical neuroanatomy and recent diffusion-weighted magnetic resonance imaging (dMRI) studies. Understanding the evolution of the VOF may shed light on its origin, function, and role in visually guided behaviors. To this end, we acquired high-resolution dMRI data from the brains of select mammalian species, including anthropoid and strepsirrhine primates, a tree shrew, rodents, and carnivores. In each species, we attempted to delineate the VOF after first locating the optic radiations in the occipital white matter. In all primate species examined, the optic radiation was flanked laterally by a prominent and coherent white matter fasciculus recognizable as the VOF. By contrast, the equivalent analysis applied to four non-primate species from the same superorder as primates (tree shrew, ground squirrel, paca, and rat) failed to reveal white matter tracts in the equivalent location. Clear evidence for a VOF was also absent in two larger carnivore species (ferret and fox). Although we cannot rule out the existence of minor or differently organized homologous fiber pathways in the non-primate species, the results suggest that the VOF has greatly expanded, or possibly emerged, in the primate lineage. This adaptation likely facilitated the evolution of unique visually guided behaviors in primates, with direct impacts on manual object manipulation, social interactions, and arboreal locomotion.

Benefits but not the dual functions of submissive signals differ between two Malagasy primates.

Many animals use formalized signals to communicate dominance relationships. In some primates, such as macaques, the function of such signals varies with dominance style. Despotic species produce unidirectional submission signals that have a dual function: in conflict contexts, they signal a willingness to withdraw, whereas in peaceful contexts, they indicate the agreement to subordination. More despotic species produce these calls to a lesser extent than less despotic species. Here, we investigated whether the use of unidirectional submission signals is also related to dominance style in two lemur species and whether signalling subordination stabilizes social relationships at the group level. Ring-tailed lemurs (Lemur catta) exhibit a more despotic dominance hierarchy than Verreaux's sifakas (Propithecus verreauxi). We observed social interactions in 75 dyads of Verreaux's sifakas and 118 dyads of ring-tailed lemurs. Both species used unidirectional submissive calls that have a dual function, potentially suggesting convergent evolution of the function of these signals in independent primate lineages. However, signalling subordination did not stabilize social relationships at the group level in both species. Additionally, subordination occurred more frequently in dyads of the more despotic ring-tailed lemurs than in Verreaux's sifakas, indicating opposite patterns to macaques in the coevolution of social traits with dominance style. This article is part of the theme issue 'The power of sound: unravelling how acoustic communication shapes group dynamics'.

Coevolution of non-homologous end joining efficiency and encephalization.

Double-strand breaks (DSB), the most difficult to repair DNA damage, are mainly repaired by non-homologous end-joining (NHEJ) or homologous recombination (HR). Previous studies seem to indicate that primates, and particularly humans, have a better NHEJ system. A distinctive feature of the primate lineage (beside longevity) is encephalization, i.e., the expansion of the brain relative to body mass (BM). Using existing transcriptome data from 34 mammalian species, we investigated the possible correlations between the expression of genes involved in NHEJ and encephalization, BM, and longevity. The same was done also for genes involved in the HR pathway. We found that, while HR gene expression is better correlated with longevity, NHEJ gene expression is strongly (and better) correlated with encephalization. Since the brain is composed of postmitotic cells, DSB repair should be mainly performed by NHEJ in this organ. Therefore, we interpret the correlation we found as an indication that NHEJ efficiency coevolved with encephalization.

Multimodal communication and audience directedness in the greeting behaviour of semi-captive African savannah elephants

Many species communicate by combining signals into multimodal combinations. Elephants live in multi-level societies where individuals regularly separate and reunite. Upon reunion, elephants often engage in elaborate greeting rituals, where they use vocalisations and body acts produced with different body parts and of various sensory modalities (e.g., audible, tactile). However, whether these body acts represent communicative gestures and whether elephants combine vocalisations and gestures during greeting is still unknown. Here we use separation-reunion events to explore the greeting behaviour of semi-captive elephants (Loxodonta africana). We investigate whether elephants use silent-visual, audible, and tactile gestures directing them at their audience based on their state of visual attention and how they combine these gestures with vocalisations during greeting. We show that elephants select gesture modality appropriately according to their audience’s visual attention, suggesting evidence of first-order intentional communicative use. We further show that elephants integrate vocalisations and gestures into different combinations and orders. The most frequent combination consists of rumble vocalisations with ear-flapping gestures, used most often between females. By showing that a species evolutionarily distant to our own primate lineage shows sensitivity to their audience’s visual attention in their gesturing and combines gestures with vocalisations, our study advances our understanding of the emergence of first-order intentionality and multimodal communication across taxa.

Insertion-Deletion Events Are Depleted in Protein Regions with Predicted Secondary Structure.

A fundamental goal in evolutionary biology and population genetics is to understand how selection shapes the fate of new mutations. Here, we test the null hypothesis that insertion-deletion (indel) events in protein-coding regions occur randomly with respect to secondary structures. We identified indels across 11,444 sequence alignments in mouse, rat, human, chimp, and dog genomes and then quantified their overlap with four different types of secondary structure-alpha helices, beta strands, protein bends, and protein turns-predicted by deep-learning methods of AlphaFold2. Indels overlapped secondary structures 54% as much as expected and were especially underrepresented over beta strands, which tend to form internal, stable regions of proteins. In contrast, indels were enriched by 155% over regions without any predicted secondary structures. These skews were stronger in the rodent lineages compared to the primate lineages, consistent with population genetic theory predicting that natural selection will be more efficient in species with larger effective population sizes. Nonsynonymous substitutions were also less common in regions of protein secondary structure, although not as strongly reduced as in indels. In a complementary analysis of thousands of human genomes, we showed that indels overlapping secondary structure segregated at significantly lower frequency than indels outside of secondary structure. Taken together, our study shows that indels are selected against if they overlap secondary structure, presumably because they disrupt the tertiary structure and function of a protein.

Changes in movement patterns in relation to sun conditions and spatial scales in wild western gorillas.

For most primates living in tropical forests, food resources occur in patchworks of different habitats that vary seasonally in quality and quantity. Efficient navigation (i.e., spatial memory-based orientation) towards profitable food patches should enhance their foraging success. The mechanisms underpinning primate navigating ability remain nonetheless mostly unknown. Using GPS long-term tracking (596days) of one group of wild western lowland gorillas (Gorilla gorilla gorilla), we investigated their ability to navigate at long distances, and tested for how the sun was used to navigate at any scale by improving landmark visibility and/or by acting as a compass. Long episodic movements ending at a distant swamp, a unique place in the home range where gorillas could find mineral-rich aquatic plants, were straighter and faster than their everyday foraging movements relying on spatial memory. This suggests intentional targeting of the swamp based on long-distance navigation skills, which can thus be efficient over a couple of kilometres. Interestingly, for both long-distance movements towards the swamp and everyday foraging movements, gorillas moved straighter under sunlight conditions even under a dense vegetation cover. By contrast, movement straightness was not markedly different when the sun elevation was low (the sun azimuth then being potentially usable as a compass) or high (so providing no directional information) and the sky was clear or overcast. This suggests that gorillas navigate their home range by relying on visual place recognition but do not use the sun azimuth as a compass. Like humans, who rely heavily on vision to navigate, gorillas should benefit from better lighting to help them identify landmarks as they move through shady forests. This study uncovers a neglected aspect of primate navigation. Spatial memory and vision might have played an important role in the evolutionary success of diurnal primate lineages.

Development and evolution of the primate neocortex from a progenitor cell perspective.

The generation of neurons in the developing neocortex is a major determinant of neocortex size. Crucially, the increase in cortical neuron numbers in the primate lineage, notably in the upper-layer neurons, contributes to increased cognitive abilities. Here, we review major evolutionary changes affecting the apical progenitors in the ventricular zone and focus on the key germinal zone constituting the foundation of neocortical neurogenesis in primates, the outer subventricular zone (OSVZ). We summarize characteristic features of the OSVZ and its key stem cell type, the basal (or outer) radial glia. Next, we concentrate on primate-specific and human-specific genes, expressed in OSVZ-progenitors, the ability of which to amplify these progenitors by targeting the regulation of the cell cycle ultimately underlies the evolutionary increase in upper-layer neurons. Finally, we address likely differences in neocortical development between present-day humans and Neanderthals that are based on human-specific amino acid substitutions in proteins operating in cortical progenitors.

Tempo and mode of gene expression evolution in the brain across primates.

Primate evolution has led to a remarkable diversity of behavioral specializations and pronounced brain size variation among species (Barton, 2012; DeCasien and Higham, 2019; Powell et al., 2017). Gene expression provides a promising opportunity for studying the molecular basis of brain evolution, but it has been explored in very few primate species to date (e.g. Khaitovich et al., 2005; Khrameeva et al., 2020; Ma et al., 2022; Somel et al., 2009). To understand the landscape of gene expression evolution across the primate lineage, we generated and analyzed RNA-seq data from four brain regions in an unprecedented eighteen species. Here, we show a remarkable level of variation in gene expression among hominid species, including humans and chimpanzees, despite their relatively recent divergence time from other primates. We found that individual genes display a wide range of expression dynamics across evolutionary time reflective of the diverse selection pressures acting on genes within primate brain tissue. Using our samples that represent a 190-fold difference in primate brain size, we identified genes with variation in expression most correlated with brain size. Our study extensively broadens the phylogenetic context of what is known about the molecular evolution of the brain across primates and identifies novel candidate genes for the study of genetic regulation of brain evolution.

Deciphering the Role of Rapidly Evolving Conserved Elements in Primate Brain Development and Exploring Their Potential Involvement in Alzheimer's Disease.

Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.

Examining infantile facial features and their influence on caretaking behaviors in free-ranging Japanese macaques (Macaca fuscata).

Facial features of immature individuals play a pivotal role in eliciting caretaking behaviors in humans. It has been posited that non-human animals share particular infantile facial features with humans, which can elicit caregivers' attention and caretaking behaviors. Nevertheless, the empirical examination of this hypothesis is extremely limited. In this study, we investigated infantile facial features in Japanese macaques (Macaca fuscata), their developmental processes, and their correlation with caretaking and infant behaviors, based on 470 facial photographs from one free-ranging group. We measured the size of facial parts and evaluated these features using non-contact procedures with the animals. The results indicated that, although some partial species differences were observed, the infantile facial features in Japanese macaques were broadly consistent with those previously observed in humans and great apes. Furthermore, half of the infant subjects displayed non-linear developmental trajectories of infantile faces, similar to those suggested in humans. However, unlike previous studies in humans, infantile faces were not significantly associated with maternal or non-maternal caretaking behaviors, nor were their developmental changes correlated with infant behavioral development. These findings indicate that while many aspects of infantile facial features are shared among particular primates, humans may have evolved a uniquely elevated preference for selecting such features among the primate lineage.

Generation of de novo miRNAs from template switching during DNA replication.

The mechanisms generating novel genes and genetic information are poorly known, even for microRNA (miRNA) genes with an extremely constrained design. All miRNA primary transcripts need to fold into a stem-loop structure to yield short gene products ([Formula: see text]22 nt) that bind and repress their mRNA targets. While a substantial number of miRNA genes are ancient and highly conserved, short secondary structures coding for entirely novel miRNA genes have been shown to emerge in a lineage-specific manner. Template switching is a DNA-replication-related mutation mechanism that can introduce complex changes and generate perfect base pairing for entire hairpin structures in a single event. Here, we show that the template-switching mutations (TSMs) have participated in the emergence of over 6,000 suitable hairpin structures in the primate lineage to yield at least 18 new human miRNA genes, that is 26% of the miRNAs inferred to have arisen since the origin of primates. While the mechanism appears random, the TSM-generated miRNAs are enriched in introns where they can be expressed with their host genes. The high frequency of TSM events provides raw material for evolution. Being orders of magnitude faster than other mechanisms proposed for de novo creation of genes, TSM-generated miRNAs enable near-instant rewiring of genetic information and rapid adaptation to changing environments.

ELOA3: A primate-specific RNA polymerase II elongation factor encoded by a tandem repeat gene cluster.

The biological role of the repetitive DNA sequences in the human genome remains an outstanding question. Recent long-read human genome assemblies have allowed us to identify a function for one of these repetitive regions. We have uncovered a tandem array of conserved primate-specific retrogenes encoding the protein Elongin A3 (ELOA3), a homolog of the RNA polymerase II (RNAPII) elongation factor Elongin A (ELOA). Our genomic analysis shows that the ELOA3 gene cluster is conserved among primates and the number of ELOA3 gene repeats is variable in the human population and across primate species. Moreover, the gene cluster has undergone concerted evolution and homogenization within primates. Our biochemical studies show that ELOA3 functions as a promoter-associated RNAPII pause-release elongation factor with distinct biochemical and functional features from its ancestral homolog, ELOA. We propose that the ELOA3 gene cluster has evolved to fulfil a transcriptional regulatory function unique to the primate lineage that can be targeted to regulate cellular hyperproliferation.

Phylogenetic comparative analysis of the cerebello-cerebral system in 34 species highlights primate-general expansion of cerebellar crura I-II

The reciprocal connections between the cerebellum and the cerebrum have been suggested to simultaneously play a role in brain size increase and to support a broad array of brain functions in primates. The cerebello-cerebral system has undergone marked functionally relevant reorganization. In particular, the lateral cerebellar lobules crura I-II (the ansiform) have been suggested to be expanded in hominoids. Here, we manually segmented 63 cerebella (34 primate species; 9 infraorders) and 30 ansiforms (13 species; 8 infraorders) to understand how their volumes have evolved over the primate lineage. Together, our analyses support proportional cerebellar-cerebral scaling, whereas ansiforms have expanded faster than the cerebellum and cerebrum. We did not find different scaling between strepsirrhines and haplorhines, nor between apes and non-apes. In sum, our study shows primate-general structural reorganization of the ansiform, relative to the cerebello-cerebral system, which is relevant for specialized brain functions in an evolutionary context.

Development of a Novel L-Asparaginase with Enhanced Properties Using Ancestral Sequence Reconstruction

The overall survival of children with acute lymphoblastic leukemia (ALL) now exceeds 90% due to the precise disease risk stratification and intensification of multi-agent chemotherapy regimens. L-asparaginase (L-ASNase) is an indispensable biotherapeutic enzyme used in treatment of pediatric ALL. Discontinuation of L-ASNase treatment due to immunological or adverse reactions leads to significantly inferior disease-free survival (73% ±7% vs 92% ± 2%; P <0.01, Gupta et al. JCO 2020). Clinical L-ASNases are bacterial in origin, derived from either Escherichia coli or Erwinia chrysanthemi and are therefore highly immunogenic. Additionally, emerging data suggests that the concurrent glutaminase activity seen in bacterial asparaginases exacerbates its liver and pancreatic toxicity, especially in adults thereby limiting its extensive use beyond children. Therefore, finding alternative sources with preferable biochemical features remains paramount. The inferior catalytic properties of human L-ASNase restricts its functionality in ALL treatment. However, characterization has shown that guinea pig (GP) L-ASNase demonstrates significant anti-tumor properties and improved enzyme kinetics despite a 69.8% amino acid sequence identity with human L-ASNase, as compared to 30% for bacterial L-ASNases. Additionally, guinea pig L-ASNase has no glutaminase activity. Thus, to overcome the major deficiencies of bacterial L-ASNase and develop a more humanized form of L-ASNase we performed ancestral sequence reconstruction (ASR). ASR is a protein drug discovery / optimization platform that predicts ancient DNA and protein sequences from extant sequences. This approach permits higher-resolution mapping potentially highlighting functional resides. The objective therefore is to exploit functional diversity refined by natural selection using ASR to identify therapeutic variants with superior properties. Recent work on guinea pig (GP) L-ASNase (Schalk et al. J Biol Chem 2014) identified a low micromolar (μM) Michaelis constant (Km) at 57.7 ± 6.4 μM essential for hydrolyzing the extra-cellular pool of asparagine ~50 μM in human blood; thus, equipping us with the ideal mammalian L-ASNase ortholog to serve as a template to resurrect ancient ASNase sequences. Ancestral L-ASNase sequences were generated utilizing 54 extant L- asparaginase sequences, aligned using MUSCLE and an evolutionary tree inferred using MrBayes. A total of 53 ancestral nodes were identified, and the ten ancestral L-ASNase variants spanning the ancient primate and guinea pig (GP) lineage were resurrected. The sequence identity of these ancestral L-ASNase variants ranged from 81 to 98% when compared to human L-ASNase. E. coli codon optimized complementary DNA (cDNA) sequences were synthesized and subcloned intoan expression vector. Correctly cloned plasmids were transformed into E. coli BL21 (DE3) cells for protein expression. An-ASNase candidates were isolated by Ni 2+ affinity chromatography and then further purified by size exclusion chromatography. Asparaginase activity of the An-ASNase candidates were determined via a modified Nessler's reagent assay by using a continuous spectroscopic enzyme-coupled assay. Three of the ancestral L-ASNase candidates, An-88, An-104 and An-107 demonstrated excellent asparaginase activity when tested at an enzyme concentration of 0.1 mg/mL and an excess substrate concentration, with activity within a comparable range to the clinically relevant E. coli and Erwinia asparaginases. An-88 had 81% similarity whereas both An-104 and An-107 ASNases had 88% identity to human L-ASNase, which is in stark contrast to the only 30% similarity seen with bacterial L-ASNases. Detailed enzymatic characterization of these three lead drug candidates is currently being performed. Preliminary cytotoxicity data testing the ancestral L-ASNase candidates against a T-cell ALL cell line CCRF-CEM showed an expected dose response, with An-107 showing the highest cytotoxicity. Our findings demonstrate that ASR is an efficient protein drug design/optimization platform that we have now extended beyond our previous work with recombinant coagulation factors (Zakas et al. Nat Biotech 2017, Knight et al. Blood Adv 2021) to therapeutic enzymes; thereby, enabling the creation of novel recombinant L-ASNases with translational implications due to improved enzymatic specificity and reduced toxicity.