Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

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.

A new in silico approach to investigate molecular aspects of factor IX missense causative mutations and their impact on the hemophilia B severity.

Factor IX (encoded by F9) is a protein in the coagulation process, where its lack or deficiency leads to hemophilia B. This condition has been much less studied than hemophilia A, especially in Latin America. We analyzed the structural and functional impact of 54 missense mutations (18 reported by us previously, and 36 other mutations from the Factor IX database) through molecular modeling approaches. To accomplish this task, we examine the electrostatic patterns, hydrophobicity/hydrophilicity, disulfide, and H-bond differences of the Factor IX structures harboring the missense mutations found, correlating them with their clinical effects. The 54 mutated sequences were modeled and their physicochemical features were determined and used as input in clusterization tools. The electrostatic pattern seems to influence in disease severity, especially for mutations investigated in epidermal growth factors 1 and 2 (EGF1/2) domains. The combined use of all physicochemical information improved the clustering of structures associated to similar phenotypes, especially for mutations from GLA and EGF1-2 domains. The effect of mutations in the disease phenotype severity seems to be a complex interplay of molecular features, each one contributing to different impacts. This highlights that previous studies and tools analyzing individually single features for single mutations are missing elements that fulfill the whole picture.

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.

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.

Reconstructing Native American population history.

The peopling of the Americas has been the subject of extensive genetic, archaeological and linguistic research; however, central questions remain unresolved. One contentious issue is whether the settlement occurred by means of a single migration or multiple streams of migration from Siberia. The pattern of dispersals within the Americas is also poorly understood. To address these questions at a higher resolution than was previously possible, we assembled data from 52 Native American and 17 Siberian groups genotyped at 364,470 single nucleotide polymorphisms. Here we show that Native Americans descend from at least three streams of Asian gene flow. Most descend entirely from a single ancestral population that we call 'First American'. However, speakers of Eskimo-Aleut languages from the Arctic inherit almost half their ancestry from a second stream of Asian gene flow, and the Na-Dene-speaking Chipewyan from Canada inherit roughly one-tenth of their ancestry from a third stream. We show that the initial peopling followed a southward expansion facilitated by the coast, with sequential population splits and little gene flow after divergence, especially in South America. A major exception is in Chibchan speakers on both sides of the Panama isthmus, who have ancestry from both North and South America.

Remembering the past - studies on evolution done by the genetics group at Universidade Federal do Rio Grande do Sul (UFRGS).

After a brief introduction about the factors that are involved in science development, and world and Brazilian evolutionary genetics, the studies developed in Porto Alegre in this area were reviewed. Four periods in the development of this group were distinguished: (a) Origins and first expansion (1949-1961); (b) Second expansion (1962-1988); (c) Third expansion (1989-2001); and (d) The last 15 years (2002-present). The international Porto Alegre Biological Evolution Workshops (PABEWs), with five biannual events from 2007 o 2015, were also mentioned. The final message stressed the importance of the maintenance of this and other Brazilian groups of research through adequate finance and recognition.

Molecular evolution and functional divergence of alcohol dehydrogenases in animals, fungi and plants.

Alcohol dehydrogenases belong to the large superfamily of medium-chain dehydrogenases/reductases, which occur throughout the biological world and are involved with many important metabolic routes. We considered the phylogeny of 190 ADH sequences of animals, fungi, and plants. Non-class III Caenorhabditis elegans ADHs were seen closely related to tetrameric fungal ADHs. ADH3 forms a sister group to amphibian, reptilian, avian and mammalian non-class III ADHs. In fishes, two main forms are identified: ADH1 and ADH3, whereas in amphibians there is a new ADH form (ADH8). ADH2 is found in Mammalia and Aves, and they formed a monophyletic group. Additionally, mammalian ADH4 seems to result from an ADH1 duplication, while in Fungi, ADH formed clusters based on types and genera. The plant ADH isoforms constitute a basal clade in relation to ADHs from animals. We identified amino acid residues responsible for functional divergence between ADH types in fungi, mammals, and fishes. In mammals, these differences occur mainly between ADH1/ADH4 and ADH3/ADH5, whereas functional divergence occurred in fungi between ADH1/ADH5, ADH5/ADH4, and ADH5/ADH3. In fishes, the forms also seem to be functionally divergent. The ADH family expansion exemplifies a neofunctionalization process where reiterative duplication events are related to new activities.

How strong was the bottleneck associated to the peopling of the Americas? New insights from multilocus sequence data.

In spite of many genetic studies that contributed for a deep knowledge about the peopling of the Americas, no consensus has emerged about important parameters such as the effective size of the Native Americans founder population. Previous estimates based on genomic datasets may have been biased by the use of admixed individuals from Latino populations, while other recent studies using samples from Native American individuals relied on approximated analytical approaches. In this study we use resequencing data for nine independent regions in a set of Native American and Siberian individuals and a full-likelihood approach based on isolation-with-migration scenarios accounting for recent flow between Asian and Native American populations. Our results suggest that, in agreement with previous studies, the effective size of the Native American population was small, most likely in the order of a few hundred individuals, with point estimates close to 250 individuals, even though credible intervals include a number as large as ~4,000 individuals. Recognizing the size of the genetic bottleneck during the peopling of the Americas is important for determining the extent of genetic markers needed to characterize Native American populations in genome-wide studies and to evaluate the adaptive potential of genetic variants in this population.

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.

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.

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.

Admixture in Latin America: geographic structure, phenotypic diversity and self-perception of ancestry based on 7,342 individuals.

The current genetic makeup of Latin America has been shaped by a history of extensive admixture between Africans, Europeans and Native Americans, a process taking place within the context of extensive geographic and social stratification. We estimated individual ancestry proportions in a sample of 7,342 subjects ascertained in five countries (Brazil, Chile, Colombia, México and Perú). These individuals were also characterized for a range of physical appearance traits and for self-perception of ancestry. The geographic distribution of admixture proportions in this sample reveals extensive population structure, illustrating the continuing impact of demographic history on the genetic diversity of Latin America. Significant ancestry effects were detected for most phenotypes studied. However, ancestry generally explains only a modest proportion of total phenotypic variation. Genetically estimated and self-perceived ancestry correlate significantly, but certain physical attributes have a strong impact on self-perception and bias self-perception of ancestry relative to genetically estimated ancestry.

The role of natural selection in human evolution - insights from Latin America.

A brief introduction considering Darwin's work, the evolutionary synthesis, and the scientific biological field around the 1970s and subsequently, with the molecular revolution, was followed by selected examples of recent investigations dealing with the selection-drift controversy. The studies surveyed included the comparison between essential genes in humans and mice, selection in Africa and Europe, and the possible reasons why females in humans remain healthy and productive after menopause, in contrast with what happens in the great apes. At the end, selected examples of investigations performed in Latin America, related to the action of selection for muscle performance, acetylation of xenobiotics, high altitude and tropical forest adaptations were considered. Despite dissenting views, the influence of positive selection in a considerable portion of the human genome cannot presently be dismissed.

Population Genetics of Jaguars (Panthera onca) in the Brazilian Pantanal: Molecular Evidence for Demographic Connectivity on a Regional Scale.

Habitat loss and fragmentation are important threats to carnivores worldwide, and can be especially intense for large predators. Jaguars have already been extirpated from over half of their original area of distribution, and few regions still maintain large populations. For these, detailed understanding is crucial for setting appropriate recovery targets in impacted areas. The Pantanal is among the best examples of a region with a large jaguar population in a healthy environment. Here, we analyzed 12 microsatellite loci to characterize genetic diversity and population structure of 52 jaguars sampled in 4 localities of the southern Pantanal, and compared them with prior studies of heavily fragmented populations of the Atlantic Forest. Although we observed some internal structure among the Pantanal localities, our results indicated that this area comprises a single population with high genetic variability. Moreover, our comparative analyses supported the hypothesis that the strong population structure observed in the Atlantic Forest derives from recent, anthropogenic fragmentation. We also observed significant but low levels of genetic differentiation between the Pantanal and Atlantic Forest populations, indicating recent connectivity between jaguars occurring in these biomes. Evidence for admixture between the Pantanal and a population on the western boundary of the Atlantic Forest corroborates the transitional nature of the latter area, where the jaguar population has already been extirpated. Our results can be used to understand jaguar population dynamics in a region that is less disturbed than the Atlantic forest, and to support the design of conservation strategies that maintain and restore natural connectivity among currently isolated areas.

High interpopulation homogeneity in Central Argentina as assessed by Ancestry Informative Markers (AIMs).

The population of Argentina has already been studied with regard to several genetic markers, but much more data are needed for the appropriate definition of its genetic profile. This study aimed at investigating the admixture patterns and genetic structure in Central Argentina, using biparental markers and comparing the results with those previously obtained by us with mitochondrial DNA (mtDNA) in the same samples. A total of 521 healthy unrelated individuals living in 13 villages of the Córdoba and San Luis provinces were tested. The individuals were genotyped for ten autosomal ancestry informative markers (AIMs). Allele frequencies were compared with those of African, European and Native American populations, chosen to represent parental contributions. The AIM estimates indicated a greater influence of the Native American ancestry as compared to previous studies in the same or other Argentinean regions, but smaller than that observed with the mtDNA tests. These differences can be explained, respectively, by different genetic contributions between rural and urban areas, and asymmetric gene flow occurred in the past. But a most unexpected finding was the marked interpopulation genetic homogeneity found in villages located in diverse geographic environments across a wide territory, suggesting considerable gene flow.

FGFR1 signaling is associated with the magnitude of morphological integration in human head shape.

OBJECTIVES: The head can be used as a model to study complex phenotypes controlled simultaneously by morphological integration (MI) due to common factors, and modular patterns caused by local factors affecting the development and functional demands of specific structures. The fibroblast growth factor and receptor system (FGF/FGFR) participates in cell communication and pattern formation in osseous tissues, among others, and there is compelling evidence from mouse model studies suggesting a role of the FGF/FGFR pathway as a covariance-generating signaling process in head development. Here we use human data to test if specific genetic variants of another gene of this pathway, the FGFR1 gene, can be associated with differences in the integration of the head. METHODS: We explored whether and how three specific variants on FGFR1, previously associated with human cephalic index, influence the pattern and level of head integration of one Native American and one admixed group from Mexico. MI, measured as the intensity of covariation among head traits, was assessed using data from three-dimensional head landmark coordinates taken on 176 individuals. RESULTS: Individuals carrying the derived allele of the rs4647905:G>C polymorphism present significantly greater levels of head MI, especially in facial structures and on the shape space where the modular portion of the covariation is explicitly removed. CONCLUSIONS: Since FGFR genes present nonconservative and tissue-specific splicing sites, they may have some effect on protein structure and performance likely involved in developmental processes responsible for the magnitude and pattern of MI in the human head.

Brief communication: 5-HTTLPR genetic diversity and mode of subsistence in Native Americans.

The relationship between the "individualism-collectivism" and the serotonin transporter functional polymorphism (5-HTTLPR), suggested in the previous reports, was tested in Native South Amerindian populations. A total of 170 individuals from 21 populations were genotyped for the 5-HTTLPR alleles. For comparative purposes, these populations were classified as individualistic (recent history of hunter-gathering) or collectivistic (agriculturalists). These two groups showed an almost identical S allele frequency (75 and 76%, respectively). The analysis of molecular variance showed no structural differences between them. Behavioral typologies like those suggested by JY Chiao and KD Blizinsky (Proc R Soc B 277 () 529-537) are always a simplification of complex phenomena and should be regarded with caution. In addition, classification of a whole nation in the individualist/collectivist dichotomy is controversial. The focus on modes of subsistence in preindustrial societies, as was tested here, may be a good alternative although the postulated association between the 5-HTTLPR S allele and the collectivist societies was not confirmed.

Demographic expansions in South America: enlightening a complex scenario with genetic and linguistic data.

Native Americans are characterized by specific and unique patterns of genetic and cultural/linguistic diversities, and this information has been used to understand patterns of geographic dispersion, and the relationship between these peoples. Particularly interesting are the Tupi and Je speaker dispersions. At present, a large number of individuals speak languages of these two stocks; for instance, Tupi-Guarani is one of the official languages in Paraguay, Bolivia, and the Mercosul economic block. Although the Tupi expansion can be compared in importance to the Bantu migration in Africa, little is known about this event relative to others. Equal and even deeper gaps exist concerning the Je-speakers' expansion. This study aims to elucidate some aspects of these successful expansions. To meet this purpose, we analyzed Native American mtDNA complete control region from nine different populations and included HVS-I sequences available in the literature, resulting in a total of 1,176 samples investigated. Evolutionary relationships were explored through median-joining networks and genetic/geographic/linguistic correlations with Mantel tests and spatial autocorrelation analyses. Both Tupi and Je showed general traces of ancient or more recent fission-fusion processes, but a very different pattern of demographic expansion. Tupi populations displayed a classical isolation-by-distance pattern, while Je groups presented an intricate and nonlinear mode of dispersion. We suggest that the collective memory and other cultural processes could be important factors influencing the fission-fusion events, which likely contributed to the genetic structure, evolution, and dispersion of Native American populations.

Fibroblast growth factor receptor 1 (FGFR1) variants and craniofacial variation in Amerindians and related populations.

OBJECTIVES: The polymorphic site rs4647905 of the FGFR1 gene was previously associated with a decrease in cephalic index (CI). Here, we evaluate the relationships between genotypes and cephalometric measurements and indices in one Mexican Native and two mestizo Mexican populations using two haplotype-tag SNPs (rs4647905 and rs3213849) that represent >85% of the FGFR1 variability, plus three other SNPs (rs2293971, rs2304000, and rs930828) situated nearby. In addition, we genotyped five South American natives, two European, one African, and one Siberian populations to evaluate their intra and intercontinental population diversity. METHODS: The five SNPs were tested and the craniofacial measurements and indices were collected using standardized procedures. Principal Component Analysis was used to verify individual/population comparisons. Associations were performed through the generalized linear model (GLM), coefficient of determination R(2) and linear regression tests. RESULTS: We found a tendency for a decrease in CI in individuals homozygous for allele rs4647905C, regardless of the population to which they belong, though the effect is more pronounced in mestizo. When the GLM analyses were performed using the absolute/linear cephalometric measurements, a statistically significant association was found between four SNPs and head length in the mestizo population. CONCLUSIONS: FGFR1 polymorphisms, especially rs4647905, can have an important role in the normal human skull variation, primarily due to their influence in head length, which would affect other cephalometric absolute/linear measures as well as indices like CI as a result of the pervasive nature of the morphological integration that characterizes the human skull.