MatchTope: A tool to predict the cross reactivity of peptides complexed with Major Histocompatibility Complex I.

The therapeutic targeting of the immune system, for example in vaccinology and cancer treatment, is a challenging task and the subject of active research. Several in silico tools used for predicting immunogenicity are based on the analysis of peptide sequences binding to the Major Histocompatibility Complex (pMHC). However, few of these bioinformatics tools take into account the pMHC three-dimensional structure. Here, we describe a new bioinformatics tool, MatchTope, developed for predicting peptide similarity, which can trigger cross-reactivity events, by computing and analyzing the electrostatic potentials of pMHC complexes. We validated MatchTope by using previously published data from in vitro assays. We thereby demonstrate the strength of MatchTope for similarity prediction between targets derived from several pathogens as well as for indicating possible cross responses between self and tumor peptides. Our results suggest that MatchTope can enhance and speed up future studies in the fields of vaccinology and cancer immunotherapy.

Variants p.Pro2063Ser and p.Arg324* co-segregate in type 3 von Willebrand disease patients from Southern Brazil.

INTRODUCTION: von Willebrand factor (VWF) is a multimeric plasma glycoprotein that plays an important role in haemostasis. von Willebrand disease (VWD) is an inherited heterogeneous bleeding disorder caused by either a quantitative or qualitative defect of VWF. Type 3 VWD, the most severe form of the disease, leads to complete quantitative VWF deficiency. AIM: The present study aims to investigate the molecular pathogenesis of type 3 VWD patients from Southern Brazil. METHODS: The VWF gene was sequenced in 26 cases clinically diagnosed with type 3 VWD by next-generation sequencing using Ion Torrent PGM. RESULTS: In 25 patients, we were able to identify both disease-causing variants. We identified 72 different variants: 31 intronic and 41 exonic. Five novel variants were found: c.6976+5G>T; c.6885_6886insC; c.3378C>T (p.Cys1126); c.3346_3347insCCA; and c.2503G>T (p.Glu835*). Variants p.Pro2063Ser and p.Arg324* co-segregated in 17 patients, 15 of them in homozygosity. CONCLUSION: Our results may contribute to the discussion on whether the variant p.Pro2063Ser is pathogenic or not. Finally, the presence of a common haplotype in patients bearing these two variants suggests a founder effect for this variant in our region.

Genomic insight into the origins and dispersal of the Brazilian coastal natives.

In the 15th century, ∼900,000 Native Americans, mostly Tupí speakers, lived on the Brazilian coast. By the end of the 18th century, the coastal native populations were declared extinct. The Tupí arrived on the east coast after leaving the Amazonian basin ∼2,000 y before present; however, there is no consensus on how this migration occurred: toward the northern Amazon and then directly to the Atlantic coast, or heading south into the continent and then migrating to the coast. Here we leveraged genomic data from one of the last remaining putative representatives of the Tupí coastal branch, a small, admixed, self-reported Tupiniquim community, as well as data of a Guaraní Mbyá native population from Southern Brazil and of three other native populations from the Amazonian region. We demonstrated that the Tupiniquim Native American ancestry is not related to any extant Brazilian Native American population already studied, and thus they could be considered the only living representatives of the extinct Tupí branch that used to settle the Atlantic Coast of Brazil. Furthermore, these data show evidence of a direct migration from Amazon to the Northeast Coast in pre-Columbian time, giving rise to the Tupí Coastal populations, and a single distinct migration southward that originated the Guaraní people from Brazil and Paraguay. This study elucidates the population dynamics and diversification of the Brazilian natives at a genomic level, which was made possible by recovering data from the Brazilian coastal population through the genomes of mestizo individuals.

Searching for ancient balanced polymorphisms shared between Neanderthals and Modern Humans.

Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

Searching for ancient balanced polymorphisms shared between Neanderthals and Modern Humans

Abstract Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

FOXP in Tetrapoda: Intrinsically Disordered Regions, Short Linear Motifs and their evolutionary significance.

The FOXP subfamily is probably the most extensively characterized subfamily of the forkhead superfamily, playing important roles in development and homeostasis in vertebrates. Intrinsically disorder protein regions (IDRs) are protein segments that exhibit multiple physical interactions and play critical roles in various biological processes, including regulation and signaling. IDRs in proteins may play an important role in the evolvability of genetic systems. In this study, we analyzed 77 orthologous FOXP genes/proteins from Tetrapoda, regarding protein disorder content and evolutionary rate. We also predicted the number and type of short linear motifs (SLIMs) in the IDRs. Similar levels of protein disorder (approximately 70%) were found for FOXP1, FOXP2, and FOXP4. However, for FOXP3, which is shorter in length and has a more specific function, the disordered content was lower (30%). Mammals showed higher protein disorders for FOXP1 and FOXP4 than non-mammals. Specific analyses related to linear motifs in the four genes showed also a clear differentiation between FOXPs in mammals and non-mammals. We predicted for the first time the role of IDRs and SLIMs in the FOXP gene family associated with possible adaptive novelties within Tetrapoda. For instance, we found gain and loss of important phosphorylation sites in the Homo sapiens FOXP2 IDR regions, with possible implication for the evolution of human speech.

FOXP in Tetrapoda: Intrinsically Disordered Regions, Short Linear Motifs and their evolutionary significance

Abstract The FOXP subfamily is probably the most extensively characterized subfamily of the forkhead superfamily, playing important roles in development and homeostasis in vertebrates. Intrinsically disorder protein regions (IDRs) are protein segments that exhibit multiple physical interactions and play critical roles in various biological processes, including regulation and signaling. IDRs in proteins may play an important role in the evolvability of genetic systems. In this study, we analyzed 77 orthologous FOXP genes/proteins from Tetrapoda, regarding protein disorder content and evolutionary rate. We also predicted the number and type of short linear motifs (SLIMs) in the IDRs. Similar levels of protein disorder (approximately 70%) were found for FOXP1, FOXP2, and FOXP4. However, for FOXP3, which is shorter in length and has a more specific function, the disordered content was lower (30%). Mammals showed higher protein disorders for FOXP1 and FOXP4 than non-mammals. Specific analyses related to linear motifs in the four genes showed also a clear differentiation between FOXPs in mammals and non-mammals. We predicted for the first time the role of IDRs and SLIMs in the FOXP gene family associated with possible adaptive novelties within Tetrapoda. For instance, we found gain and loss of important phosphorylation sites in the Homo sapiens FOXP2 IDR regions, with possible implication for the evolution of human speech.

Genetic signature of natural selection in first Americans.

When humans moved from Asia toward the Americas over 18,000 y ago and eventually peopled the New World they encountered a new environment with extreme climate conditions and distinct dietary resources. These environmental and dietary pressures may have led to instances of genetic adaptation with the potential to influence the phenotypic variation in extant Native American populations. An example of such an event is the evolution of the fatty acid desaturases (FADS) genes, which have been claimed to harbor signals of positive selection in Inuit populations due to adaptation to the cold Greenland Arctic climate and to a protein-rich diet. Because there was evidence of intercontinental variation in this genetic region, with indications of positive selection for its variants, we decided to compare the Inuit findings with other Native American data. Here, we use several lines of evidence to show that the signal of FADS-positive selection is not restricted to the Arctic but instead is broadly observed throughout the Americas. The shared signature of selection among populations living in such a diverse range of environments is likely due to a single and strong instance of local adaptation that took place in the common ancestral population before their entrance into the New World. These first Americans peopled the whole continent and spread this adaptive variant across a diverse set of environments.

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.

Variability of innate immune system genes in Native American populations-relationship with history and epidemiology.

OBJECTIVES: The immune system of a host, defending him/her against invading pathogens, has two main subsystems: innate immunity and acquired immunity. There are several evidences showing that Native American populations are immunologically different from non-Native populations. Our aim was to describe the variability of innate immune system genes in Native American populations. MATERIALS AND METHODS: We investigated heterozygozities and patterns of population differentiation (FST ) of 14 polymorphisms related to the innate immune response in five Native American populations (Aché, Guarani-Kaiowá, Guarani-Ñandeva, Kaingang, and Xavante) and the results were compared with the three major world population data (YRI, CEU, and CHB) available at the 1,000 genomes database. RESULTS: Mean heterozygosities ranged between 0.241 ± 0.057 (Aché) and 0.343 ± 0.033 (Kaingang), but no significant differences were observed (Friedman test, P = 0.197). Mean heterozygosities were also not significantly different when Amerindians were pooled and compared with the 1000 genomes populations (Friedman test, P = 0.506). When the Native American populations were grouped as Amerindians, a significantly higher FST value (0.194) was observed between the Amerindian and African populations. The Ewens-Watterson neutrality test showed that these markers are not under strong selective pressure. DISCUSSION: Native American populations present similar levels of heterozygosity as those of other continents, but are different from Africans in the frequency of polymorphisms of innate immune genes. This higher differentiation is probably due to demographic processes that occurred during the out-of-Africa event.

Genetic Variations in the TP53 Pathway in Native Americans Strongly Suggest Adaptation to the High Altitudes of the Andes.

The diversity of the five single nucleotide polymorphisms located in genes of the TP53 pathway (TP53, rs1042522; MDM2, rs2279744; MDM4, rs1563828; USP7, rs1529916; and LIF, rs929271) were studied in a total of 282 individuals belonging to Quechua, Aymara, Chivay, Cabanaconde, Yanke, Taquile, Amantani, Anapia, Uros, Guarani Ñandeva, and Guarani Kaiowá populations, characterized as Native American or as having a high level (> 90%) of Native American ancestry. In addition, published data pertaining to 100 persons from five other Native American populations (Surui, Karitiana, Maya, Pima, and Piapoco) were analyzed. The populations were classified as living in high altitude (≥ 2,500 m) or in lowlands (< 2,500 m). Our analyses revealed that alleles USP7-G, LIF-T, and MDM2-T showed significant evidence that they were selected for in relation to harsh environmental variables related to high altitudes. Our results show for the first time that alleles of classical TP53 network genes have been evolutionary co-opted for the successful human colonization of the Andes.

Genetic evidence for two founding populations of the Americas.

Genetic studies have consistently indicated a single common origin of Native American groups from Central and South America. However, some morphological studies have suggested a more complex picture, whereby the northeast Asian affinities of present-day Native Americans contrast with a distinctive morphology seen in some of the earliest American skeletons, which share traits with present-day Australasians (indigenous groups in Australia, Melanesia, and island Southeast Asia). Here we analyse genome-wide data to show that some Amazonian Native Americans descend partly from a Native American founding population that carried ancestry more closely related to indigenous Australians, New Guineans and Andaman Islanders than to any present-day Eurasians or Native Americans. This signature is not present to the same extent, or at all, in present-day Northern and Central Americans or in a ∼12,600-year-old Clovis-associated genome, suggesting a more diverse set of founding populations of the Americas than previously accepted.

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.

Interethnic admixture and the evolution of Latin American populations.

A general introduction to the origins and history of Latin American populations is followed by a systematic review of the data from molecular autosomal assessments of the ethnic/continental (European, African, Amerindian) ancestries for 24 Latin American countries or territories. The data surveyed are of varying quality but provide a general picture of the present constitution of these populations. A brief discussion about the applications of these results (admixture mapping) is also provided. Latin American populations can be viewed as natural experiments for the investigation of unique anthropological and epidemiological issues.

Evolutionary history of chordate PAX genes: dynamics of change in a complex gene family.

Paired box (PAX) genes are transcription factors that play important roles in embryonic development. Although the PAX gene family occurs in animals only, it is widely distributed. Among the vertebrates, its 9 genes appear to be the product of complete duplication of an original set of 4 genes, followed by an additional partial duplication. Although some studies of PAX genes have been conducted, no comprehensive survey of these genes across the entire taxonomic unit has yet been attempted. In this study, we conducted a detailed comparison of PAX sequences from 188 chordates, which revealed restricted variation. The absence of PAX4 and PAX8 among some species of reptiles and birds was notable; however, all 9 genes were present in all 74 mammalian genomes investigated. A search for signatures of selection indicated that all genes are subject to purifying selection, with a possible constraint relaxation in PAX4, PAX7, and PAX8. This result indicates asymmetric evolution of PAX family genes, which can be associated with the emergence of adaptive novelties in the chordate evolutionary trajectory.