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.

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.

Analysis of a Protein Network Related to Copy Number Variations in Autism Spectrum Disorder.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder, with strong genetic influences as evidenced by its high heritability. Submicroscopic variations (ranging from one kilobase to several megabases) in DNA, called copy number variations (CNVs), have been associated with psychiatric diseases, including ASD. We aimed to identify CNVs in children diagnosed with idiopathic ASD. We used microarray-based comparative genomic hybridization analysis to detect the CNVs, and bioinformatic tools to evaluate their pathogenic potential, based on predicted functional aspects. Using combined cytogenetic and bioinformatic tools, we identified an autism network of genes/proteins related to the CNVs. Among the 40 children analyzed, we found 14 potentially pathogenic CNVs, including those previously associated with ASD (located at 16p11.2, 15q11.2, and 7p21 regions). We suggest that the most relevant biological process and functional attributes involve olfactory receptors. The CNV-related autism network comprised 90 proteins and 754 nodes and indicated the family of olfactory receptors as a significant pathway in ASD. Olfactory receptors were previously associated with neurologic diseases, and they are possibly related to cognition. This integrative analysis that combines cytogenetics and bioinformatics is a promising approach to understand complex conditions such as ASD.

Serotonin, behavior, and natural selection in New World monkeys.

Traits that undergo massive natural selection pressure, with multiple events of positive selection, are hard to find. Social behaviour, in social animals, is crucial for survival, and genetic networks involved in behaviour, such as those of serotonin (5-HT) and other neurotransmitters, must be the target of natural selection. Here, we used molecular analyses to search for signals of positive selection in the 5-HT system and found such signals in the M3-M4 intracellular domain of the 5-HT3A serotonin receptor subunit (HTR3A) in primates. We detected four amino acid sites with signs of putatively positive selection (398, 403, 432 and 416); the first three showed indications of being selected in New World monkeys (NWM, Platyrrhini), specifically in the Callitrichinae branch. Additionally, we searched for associations of these amino acid variants with social behavioural traits (i.e. sex-biased dispersal, dominance and social monogamy) using classical and Bayesian methods, and found statistically significant associations for unbiased sex dispersal (398L and 416S), unbiased sex dominance (416S) and social monogamy (416S), as well as significant positive correlation between female dispersal and 403G. Furthermore, we found putatively functional protein motifs determined by three selected sites, of which we highlight a ligand motif to GSK3 in the 416S variant, appearing only in Platyrrhini. 5-HT, 5-HT3A receptor and GSK3 are part of a network that participates in neurodevelopment and regulates behaviour, among other functions. We suggest that these genetic variations, together with those found in other neurotransmitter systems, must contribute to adaptive behaviours and consequently to fitness in NWMs.

Contrasting patterns of RUNX2 repeat variations are associated with palate shape in phyllostomid bats and New World primates.

Establishing the genetic basis that underlies craniofacial variability in natural populations is one of the main topics of evolutionary and developmental studies. One of the genes associated with mammal craniofacial variability is RUNX2, and in the present study we investigated the association between craniofacial length and width and RUNX2 across New World bats (Phyllostomidae) and primates (Catarrhini and Platyrrhini). Our results showed contrasting patterns of association between the glutamate/alanine ratios (Q/A ratio) and palate shape in these highly diverse groups. In phyllostomid bats, we found an association between shorter/broader faces and increase of the Q/A ratio. In New World monkeys (NWM) there was a positive correlation of increasing Q/A ratios to more elongated faces. Our findings reinforced the role of the Q/A ratio as a flexible genetic mechanism that would rapidly change the time of skull ossification throughout development. However, we propose a scenario in which the influence of this genetic adjustment system is indirect. The Q/A ratio would not lead to a specific phenotype, but throughout the history of a lineage, would act along with evolutionary constraints, as well as other genes, as a facilitator for adaptive morphological changes.

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.

A tale of agriculturalists and hunter-gatherers: Exploring the thrifty genotype hypothesis in native South Americans.

OBJECTIVES: To determine genetic differences between agriculturalist and hunter-gatherer southern Native American populations for selected metabolism-related markers and to test whether Neel's thrifty genotype hypothesis (TGH) could explain the genetic patterns observed in these populations. MATERIALS AND METHODS: 375 Native South American individuals from 17 populations were genotyped using six markers (APOE rs429358 and rs7412; APOA2 rs5082; CD36 rs3211883; TCF7L2 rs11196205; and IGF2BP2 rs11705701). Additionally, APOE genotypes from 39 individuals were obtained from the literature. AMOVA, main effects, and gene-gene interaction tests were performed. RESULTS: We observed differences in allele distribution patterns between agriculturalists and hunter-gatherers for some markers. For instance, between-groups component of genetic variance (FCT ) for APOE rs429358 showed strong differences in allelic distributions between hunter-gatherers and agriculturalists (p = 0.00196). Gene-gene interaction analysis indicated that the APOE E4/CD36 TT and APOE E4/IGF2BP2 A carrier combinations occur at a higher frequency in hunter-gatherers, but this combination is not replicated in archaic (Neanderthal and Denisovan) and ancient (Anzick, Saqqaq, Ust-Ishim, Mal'ta) hunter-gatherer individuals. DISCUSSION: A complex scenario explains the observed frequencies of the tested markers in hunter-gatherers. Different factors, such as pleotropic alleles, rainforest selective pressures, and population dynamics, may be collectively shaping the observed genetic patterns. We conclude that although TGH seems a plausible hypothesis to explain part of the data, other factors may be important in our tested populations.

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.

Origins and evolvability of the PAX family.

The paired box (PAX) family of transcription/developmental genes plays a key role in numerous stages of embryonic development, as well as in adult organogenesis. There is evidence linking the acquisition of a paired-like DNA binding domain (PD) to domestication of a Tc1/mariner transposon. Further duplication/deletion processes led to at least five paralogous metazoan protein groups, which can be classified into two supergroups, PAXB-like or PAXD-like, using ancestral defining structures; the PD plus an octapeptide motif (OP) and a paired-type homeobox DNA binding domain (PTHD), producing the PD-OP-PTHD structure characteristic of the PAXB-like group, whereas an additional domain, the paired-type homeodomain tail (PHT), is present in the PAXD-like group, producing a PD-OP-PTHD-PHT structure. We examined their patterns of distribution in various species, using both available data and new bioinformatic analyses, including vertebrate PAX genes and their shared and specific functions, as well as inter- and intraspecific variability of PAX in primates. These analyses revealed a relatively conserved PAX network, accompanied by specific changes that led to adaptive novelties. Therefore, both stability and evolvability shaped the molecular evolution of this key transcriptional network.

MSX1 and PAX9 investigation in monozygotic twins with variable expression of tooth agenesis.

Non-syndromic agenesis of permanent teeth is one of the most common anomalies in human development, a multifactorial characteristic caused by genetic and environmental factors. We describe a pair of monozygotic twins who showed second premolar and third molar agenesis, albeit with different expressions. We studied the DNA of two genes, paired domain box gene 9 (PAX9) and muscle segment homeodomain-homeobox1 (MSX1), encoding transcription factors that earlier studies found were involved in the manifestation of this condition. No specific causative mutation was found. However, we detected a C→T change in MSX1 exon 2 in both twins, suggesting that this polymorphism might be involved in the trait's expression.

The cognitive ability of extinct hominins: bringing down the hierarchy using genomic evidences.

BACKGROUND: The availability of the full genomes of Homo sapiens, Homo neanderthalensis, and Denisovans, as well as modern bioinformatic tools, are opening new possibilities for the understanding of the differences and similarities present in these taxa. METHODS: We searched for cognitive genes, examined their status in the genomes of these three entities. All substitutions present among them were retrieved. RESULTS: We found 93 nonsynonymous substitutions in 51 cognitive genes, in which the derived allele was present in archaic and modern humans and the ancestral allele in other nonhuman primates. CONCLUSIONS: The general picture obtained is of similarity in cognitive genes between extinct and extant humans.

Predicting Homo pigmentation phenotype through genomic data: from Neanderthal to James Watson.

BACKGROUND: Human pigmentation is regulated by several genes acting at different stages of melanin formation. Functional and association studies have elucidated the role of several of these genes in pigmentation phenotypes. Forensic and evolutionary studies can benefit from this knowledge. OBJECTIVES: To evaluate the reliability of the prediction of pigmentation phenotypes using a large database of genetic markers in individuals with known phenotypes; and from this try to predict the pigmentation phenotype of prehistoric Homo specimens and of contemporary individuals whose visible phenotypes are not known. METHODS: We compared predicted and observed phenotypic data through an analysis of 124 single nucleotide polymorphisms in 33 genic and seven intergenic regions of 30 subjects, five of them prehistoric, whose complete nuclear genomes are available in UCSC and PSU UCSC public databases. RESULTS: For the molecular predicted versus observed phenotypes, the percentage of agreement was as follows: freckles: 91; skin: 64; hair: 44; eyes: 36; total: 59; while the molecular predicted versus probable (no visible observation available; inferences based on ethnic population characteristics) it was, respectively, 83, 60, 42, 67, and 63. The difference between two sets is statistically nonsignificant (P = 0.75). CONCLUSION: To our knowledge, this is the first article to examine the effect of a large number of genetics markers for phenotype prediction. The approach could be useful for forensic applications, as well as for the determination of possible phenotypes of extinct prehistoric individuals.

A new subhaplogroup of native American Y-Chromosomes from the Andes.

The human Y chromosome contains highly informative markers for making historical inferences about the pre-Columbian peopling of Americas. However, the scarcity of these markers has limited its use in the inference of shared ancestry and past migrations relevant to the origin of the culturally and biologically diverse Native Americans. To identify new single nucleotide polymorphisms (SNPs) and increase the phylogenetic resolution of the major haplogroup Q found in the Americas, we have performed a search for new polymorphisms based on sequencing divergent Y chromosomes identified by microsatellite haplotype analysis. Using this approach, a new Y-SNP (SA01) has been identified in the Andean populations of South America, allowing for the detection of a new sublineage of Q1a3a. This sublineage displays a less complex phylogeographic network of associated microsatellites and more restricted geographic occurrence, and is given the designation Q1a3a4. This result indicates that our approach can be successfully used to identify sublineages of interest in a specific region that allow the investigation of particular histories of human populations.

Distribution of Y-chromosome Q lineages in Native Americans.

OBJECTIVES: This investigation was performed to identify and evaluate the distribution of all 15 Y-chromosome lineages belonging to the Q clade in a sample of natives from South America. METHODS: One hundred and forty-eight individuals from 20 Native American populations, as well as 24 Asian samples including Eskimos, were tested with 18 biallelic loci that can identify all currently known lineages of the Y-Chromosome Q clade. Sequencing was performed in part of the sample (∼180,000 nucleotides, which detected, for instance, several downstream markers related to the Q1a3a lineage). RESULTS: No new mutation was found and Q1a3a was consistently found in high frequencies in all populations, followed at a much lower frequency by Q1a3*, while Q1a3a derived-lineages are probably population/tribe/region-specific. CONCLUSION: The number of basal Y chromosome lineages in North America is apparently higher than in South America due probably to a bottleneck during the South American colonization and/or more recent Circum-Arctic gene flow.

Genetic variation among major human geographic groups supports a peculiar evolutionary trend in PAX9.

A total of 172 persons from nine South Amerindian, three African and one Eskimo populations were studied in relation to the Paired box gene 9 (PAX9) exon 3 (138 base pairs) as well as its 5'and 3'flanking intronic segments (232 bp and 220 bp, respectively) and integrated with the information available for the same genetic region from individuals of different geographical origins. Nine mutations were scored in exon 3 and six in its flanking regions; four of them are new South American tribe-specific singletons. Exon3 nucleotide diversity is several orders of magnitude higher than its intronic regions. Additionally, a set of variants in the PAX9 and 101 other genes related with dentition can define at least some dental morphological differences between Sub-Saharan Africans and non-Africans, probably associated with adaptations after the modern human exodus from Africa. Exon 3 of PAX9 could be a good molecular example of how evolvability works.