From molecular variations to behavioral adaptations: Unveiling adaptive epistasis in primate oxytocin system.

OBJECTIVE: Our primary objective was to investigate the variability of oxytocin (OT) and the GAMEN binding motif within the LNPEP oxytocinase in primates. MATERIALS AND METHODS: We sequenced the LNPEP segment encompassing the GAMEN motif in 34 Platyrrhini species, with 21 of them also sequenced for the OT gene. Our dataset was supplemented with primate sequences of LNPEP, OT, and the oxytocin receptor (OTR) sourced from public databases. Evolutionary analysis and coevolution predictions were made followed by the macroevolution analysis of relevant amino acids associated with phenotypic traits, such as mating systems, parental care, and litter size. To account for phylogenetic structure, we utilized two distinct statistical tests. Additionally, we calculated binding energies focusing on the interaction between Callithtrix jacchus VAMEN and Pro8OT. RESULTS: We identified two novel motifs (AAMEN and VAMEN), challenging the current knowledge of motif conservation in placental mammals. Coevolution analysis demonstrated a correlation between GAMEN, AAMEN, and VAMEN and their corresponding OTs and OTRs. Callithrix jacchus exhibited a higher binding energy between VAMEN and Pro8OT than orthologous molecules found in humans (GAMEN and Leu8OT). DISCUSSION: The coevolution of AAMEN and VAMEN with their corresponding OTs and OTRs suggests a functional relationship that could have contributed to specific reproductive and adaptive behaviors, including paternal care, social monogamy, and twin births, prominent traits in Cebidae species, such as marmosets and tamarins. Our findings underscore the coevolution of taxon-specific amino acids among the three studied molecules, shedding light on the oxytocinergic system as an adaptive epistatic repertoire in primates.

Insights into the evolutionary history of the most skilled tool-handling platyrrhini monkey: Sapajus libidinosus from the Serra da Capivara National Park.

Sapajus libidinosus members of the Pedra Furada group, living in the Serra da Capivara National Park, use stone tools in a wider variety of behaviors than any other living animal, except humans. To rescue the evolutionary history of the Caatinga S. libidinosus and identify factors that may have contributed to the emergence and maintenance of their tool-use culture, we conducted fieldwork seasons to obtain biological samples of these capuchin monkeys. UsingCYTBsequences, we show a discrete but constant population growth from the beginning of the Holocene to the present, overlapping the emergence of the Caatinga biome. Our habitat suitability reconstruction reports the presence of plants whose hard fruits, seeds, or roots are processed by capuchins using tools. TheS. libidinosusindividuals in the Caatinga were capable of dynamically developing and maintaining their autochthonous culture thanks to: a) cognitive capacity to generate and execute innovation under selective pressure; b) tolerance favoring learning and cultural inheritance; c) an unknown genetic repertoire that underpins the adaptive traits; d) a high degree of terrestriality; e) presence and abundance of natural resources, which makes some places "hot spots" for innovation, and cultural diversification within a relatively short time.

Exploring the diversity of AVPR2 in Primates and its evolutionary implications.

The current study focuses on the investigation of AVPR2 (VTR2C) protein-coupled receptor variants specific to different primate taxa. AVPR2 is activated by the neurohormone AVP, which modulates physiological processes, including water homeostasis. Our findings reveal positive selection at three AVPR2 sites at positions 190, 250, and 346. Variation at position 250 is associated with human Congenital Nephrogenic Diabetes Insipidus (cNDI), a condition characterized by excessive water loss. Other 13 functional sites with potential adaptive relevance include positions 185, 202, 204, and 252 associated with cNDI. We identified SH3-binding motifs in AVPR2's ICL3 and N-terminus domains, with some losses observed in clades of Cercopithecidae, Callitrichinae, and Atelidae. SH3-binding motifs are crucial in regulating cellular physiology, indicating that the differences may be adaptive. Co-evolution was found between AVPR2 residues and those in the AVP signal peptide/Neurophysin-2 and AQP2, other molecules in the same signaling cascade. No significant correlation was found between these Primates' taxon-specific variants and the bioclimatic variables of the areas where they live. Distinct co-evolving amino acid sequences in functional sites were found in Platyrrhini and Catarrhini, which may have adaptive implications involving glucocorticoid hormones, suggesting varied selective pressures. Further studies are required to confirm these results.

Molecular evolutionary insights from PRLR in mammals.

Prolactin (PRL) is a pleiotropic neurohormone secreted by the mammalian pituitary gland into the blood, thus reaching many tissues and organs beyond the brain. PRL binds to its receptor, PRLR, eliciting a molecular signaling cascade. This system modulates essential mammalian behaviors and promotes notable modifications in the reproductive female tissues and organs. Here, we explore how the intracellular domain of PRLR (PRLR-ICD) modulates the expression of the PRLR gene. Despite differences in the reproductive strategies between eutherian and metatherian mammals, there is no clear distinction between PRLR-ICD functional motifs. However, we found selection signatures that showed differences between groups, with many conserved functional elements strongly maintained through purifying selection across the class Mammalia. We observed a few residues under relaxed selection, the levels of which were more pronounced in Eutheria and particularly striking in primates (Simiiformes), which could represent a pre-adaptive genetic element protected from purifying selection. Alternative, new motifs, such as YLDP (318-321) and others with residues Y283 and Y290, may already be functional. These motifs would have been co-opted in primates as part of a complex genetic repertoire related to some derived adaptive phenotypes, but these changes would have no impact on the primordial functions that characterize the mammals as a whole and that are related to the PRL-PRLR system.

ACE2 diversity in placental mammals reveals the evolutionary strategy of SARS-CoV-2.

The recent emergence of SARS-CoV-2 is responsible for the current pandemic of COVID-19, which uses the human membrane protein ACE2 as a gateway to host-cell infection. We performed a comparative genomic analysis of 70 ACE2 placental mammal orthologues to identify variations and contribute to the understanding of evolutionary dynamics behind this successful adaptation to infect humans. Our results reveal that 4% of the ACE2 sites are under positive selection, all located in the catalytic domain, suggesting possibly taxon-specific adaptations related to the ACE2 function, such as cardiovascular physiology. Considering all variable sites, we selected 30 of them located at the critical ACE2 binding sites to the SARS-CoV-like viruses for analysis in more detail. Our results reveal a relatively high diversity of ACE2 between placental mammal species, while showing no polymorphism within human populations, at least considering the 30 inter-species variable sites. A perfect scenario for natural selection favored this opportunistic new coronavirus in its trajectory of infecting humans. We suggest that SARS-CoV-2 became a specialist coronavirus for human hosts. Differences in the rate of infection and mortality could be related to the innate immune responses, other unknown genetic factors, as well as non-biological factors.

AVPR1b variation and the emergence of adaptive phenotypes in Platyrrhini primates.

Platyrrhini (New World monkeys, NWm) are a group of primates characterized by behavioral and reproductive traits that are otherwise uncommon among primates, including social monogamy, direct paternal care, and twin births. As a consequence, the study of Platyrrhine primates is an invaluable tool for the discovery of the genetic repertoire underlying these taxon-specific traits. Recently, high conservation of vasopressin (AVP) sequence, in contrast with high variability of oxytocin (OXT), has been described in NWm. AVP and OXT functions are possible due to interaction with their receptors: AVPR1a, AVPR1b, AVPR2, and OXTR; and the variability in this system is associated with the traits mentioned above. Understanding the variability in the receptors is thus fundamental to understand the function and evolution of the system as a whole. Here we describe the variability of AVPR1b coding region in 20 NWm species, which is well-known to influence behavioral traits such as aggression, anxiety, and stress control in placental mammals. Our results indicate that 4% of AVPR1b sites may be under positive selection and a significant number of sites under relaxed selective constraint. Considering the known role of AVPR1b, we suggest that some of the changes described here for the Platyrrhini may be a part of the genetic repertoire connected with the complex network of neuroendocrine mechanisms of AVP-OXT system in the modulation of the HPA axis. Thus, these changes may have promoted the emergence of social behaviors such as direct paternal care in socially monogamous species that are also characterized by small body size and twin births.

Evidence for Association Between OXTR Gene and ASD Clinical Phenotypes.

Autism spectrum disorder (ASD) is an early-onset neurodevelopmental disorder characterized by impairments in social behaviors and communication. Oxytocin and its signaling pathway are related to a range of human behaviors, from facial expression recognition to aggressive behaviors, and have been suggested as involved in the etiology of ASD. Our aim was to evaluate the influence of two polymorphisms (rs1042778, rs53576) at the oxytocin receptor gene (OXTR) on ASD diagnosis and on specific ASD-related clinical symptoms (seizures, panic, and aggressive behaviors). We also assessed if these SNPs could be related to changes in OXTR availability and functionality using a bioinformatic approach. The sample was composed by 209 probands with ASD and their biological parents. Family-based approach and logistic regression models were used to investigated the outcomes. We observed that panic and aggressive behaviors were nominally associated with presence of rs1042778 T allele (P = 0.019/Pcorr = 0.114; P = 0.046/Pcorr = 0.276 respectively). Also, in the family-based analysis, a trend towards association with ASD susceptibility was observed for rs1042778 (G allele) (P = 0.066). In a bioinformatic approach, we demonstrated that rs1042778 G allele is determinant for the binding of the transcription factor MAZ, suggesting that when the T allele is present, the absence of MAZ binding might be associated with lower transcription levels of the OXTR gene. The overall findings suggest that the OXTR gene may play a role in ASD diagnosis and some of its clinical phenotypes, supported by previous animal and clinical studies. Further investigations are necessary to replicate our findings and fully understand the effects of the oxytocin pathway on ASD.

Oxytocin and arginine vasopressin systems in the domestication process.

Domestication is of unquestionable importance to the technological revolution that has given rise to modern human societies. In this study, we analyzed the DNA and protein sequences of six genes of the oxytocin and arginine vasopressin systems (OXT-OXTR; AVP-AVPR1a, AVPR1b and AVPR2) in 40 placental mammals. These systems play an important role in the control of physiology and behavior. According to our analyses, neutrality does not explain the pattern of molecular evolution found in some of these genes. We observed specific sites under positive selection in AVPR1b (ω = 1.429, p = 0.001) and AVPR2 (ω= 1.49, p = 0.001), suggesting that they could be involved in behavior and physiological changes, including those related to the domestication process. Furthermore, AVPR1a, which plays a role in social behavior, is under relaxed selective constraint in domesticated species. These results provide new insights into the nature of the domestication process and its impact on the OXT-AVP system.

Oxytocin and arginine vasopressin systems in the domestication process

Abstract Domestication is of unquestionable importance to the technological revolution that has given rise to modern human societies. In this study, we analyzed the DNA and protein sequences of six genes of the oxytocin and arginine vasopressin systems (OXT-OXTR; AVP-AVPR1a, AVPR1b and AVPR2) in 40 placental mammals. These systems play an important role in the control of physiology and behavior. According to our analyses, neutrality does not explain the pattern of molecular evolution found in some of these genes. We observed specific sites under positive selection in AVPR1b (ω = 1.429, p = 0.001) and AVPR2 (ω= 1.49, p = 0.001), suggesting that they could be involved in behavior and physiological changes, including those related to the domestication process. Furthermore, AVPR1a, which plays a role in social behavior, is under relaxed selective constraint in domesticated species. These results provide new insights into the nature of the domestication process and its impact on the OXT-AVP system.

Functional New World monkey oxytocin forms elicit an altered signaling profile and promotes parental care in rats.

The neurohormone oxytocin is a key player in the modulation of reproductive and social behavioral traits, such as parental care. Recently, a correlation between different forms of oxytocin and behavioral phenotypes has been described in the New World Monkeys (NWMs). Here, we demonstrate that, compared with the Leu8OXT found in most placental mammals, the Cebidae Pro8OXT and Saguinus Val3Pro8OXT taxon-specific variants act as equi-efficacious agonists for the Gq-dependent pathway but are weaker agonists for the ß-arrestin engagement and subsequent endocytosis toward the oxytocin receptor (OXTR). Upon interaction with the AVPR1a, Pro8OXT and the common Leu8OXT yielded similar signaling profiles, being equally efficacious on Gq and ß-arrestin, while Val3Pro8OXT showed reduced relative efficacy toward ß-arrestin. Intranasal treatment with either of the variants increased maternal behavior and also promoted unusual paternal care in rats, as measured by pup-retrieval tests. We therefore suggest that Val3Pro8OXT and Pro8OXT are functional variants, which might have been evolutionarily co-opted as an essential part of the adaptive genetic repertoire that allowed the emergence of taxon-specific complex social behaviors, such as intense parental care in the Cebidae and the genus Saguinus.

Progesterone Response Element Variation in the OXTR Promoter Region and Paternal Care in New World Monkeys.

Paternal care is a complex social behavior common in primate species with socially monogamous mating systems and twin births. Evolutionary causes and consequences of such behavior are not well understood, nor are their neuroendocrine and genetic bases. However, the neuropeptide oxytocin (OXT) and its receptor (OXTR) are associated with parental care in mammalian lineages. Here we investigated the interspecific variation in the number of progesterone response elements (PREs) in the OXTR promoter region of 32 primate species, correlating genetic data with behavior, social systems, and ecological/life-history parameters, while controlling for phylogeny. We verified that PREs are only present in New World monkeys and that PRE number is significantly correlated with the presence of paternal care in this branch. We suggest that PRE number could be an essential part of the genetic repertoire that allowed the emergence of taxon-specific complex social behaviors, such as paternal care in marmosets and tamarins.

Oxytocin and arginine vasopressin receptor evolution: implications for adaptive novelties in placental mammals.

Oxytocin receptor (OXTR) and arginine vasopressin receptors (AVPR1a, AVPR1b, and AVPR2) are paralogous genes that emerged through duplication events; along the evolutionary timeline, owing to speciation, numerous orthologues emerged as well. In order to elucidate the evolutionary forces that shaped these four genes in placental mammals and to reveal specific aspects of their protein structures, 35 species were selected. Specifically, we investigated their molecular evolutionary history and intrinsic protein disorder content, and identified the presence of short linear interaction motifs. OXTR seems to be under evolutionary constraint in placental mammals, whereas AVPR1a, AVPR1b, and AVPR2 exhibit higher evolutionary rates, suggesting that they have been under relaxed or experienced positive selection. In addition, we describe here, for the first time, that the OXTR, AVPR1a, AVPR1b, and AVPR2 mammalian orthologues preserve their disorder content, while this condition varies among the paralogues. Finally, our results reveal the presence of short linear interaction motifs, indicating possible functional adaptations related to physiological and/or behavioral taxa-specific traits.

Evolutionary pattern in the OXT-OXTR system in primates: coevolution and positive selection footprints.

Oxytocin is a nonapeptide involved in a wide range of physiologic and behavioral functions. Until recently, it was believed that an unmodified oxytocin sequence was present in all placental mammals. This study analyzed oxytocin (OXT) in 29 primate species and the oxytocin receptor (OXTR) in 21 of these species. We report here three novel OXT forms in the New World monkeys, as well as a more extensive distribution of a previously described variant (Leu8Pro). In structural terms, these OXTs share the same three low-energy conformations in solution during molecular dynamic simulations, with subtle differences in their side chains. A consistent signal of positive selection was detected in the Cebidae family, and OXT position 8 showed a statistically significant (P = 0.013) correlation with litter size. Several OXTR changes were identified, some of them promoting gain or loss of putative phosphorylation sites, with possible consequences for receptor internalization and desensitization. OXTR amino acid sites are under positive selection, and intramolecular and intermolecular coevolutionary processes with OXT were also detected. We suggest that some New World monkey OXT-OXTR forms can be correlated to male parental care through the increase of cross-reactivity with its correlated vasopressin system.

Implications of the admixture process in skin color molecular assessment.

The understanding of the complex genotype-phenotype architecture of human pigmentation has clear implications for the evolutionary history of humans, as well as for medical and forensic practices. Although dozens of genes have previously been associated with human skin color, knowledge about this trait remains incomplete. In particular, studies focusing on populations outside the European-North American axis are rare, and, until now, admixed populations have seldom been considered. The present study was designed to help fill this gap. Our objective was to evaluate possible associations of 18 single nucleotide polymorphisms (SNPs), located within nine genes, and one pseudogene with the Melanin Index (MI) in two admixed Brazilian populations (Gaucho, N = 352; Baiano, N = 148) with different histories of geographic and ethnic colonization. Of the total sample, four markers were found to be significantly associated with skin color, but only two (SLC24A5 rs1426654, and SLC45A2 rs16891982) were consistently associated with MI in both samples (Gaucho and Baiano). Therefore, only these 2 SNPs should be preliminarily considered to have forensic significance because they consistently showed the association independently of the admixture level of the populations studied. We do not discard that the other two markers (HERC2 rs1129038 and TYR rs1126809) might be also relevant to admixed samples, but additional studies are necessary to confirm the real importance of these markers for skin pigmentation. Finally, our study shows associations of some SNPs with MI in a modern Brazilian admixed sample, with possible applications in forensic genetics. Some classical genetic markers in Euro-North American populations are not associated with MI in our sample. Our results point out the relevance of considering population differences in selecting an appropriate set of SNPs as phenotype predictors in forensic practice.