Estimation of divergence times for major lineages of primate species

Isochronous singing in 3 crested gibbon species (Nomascus spp.)

The search for common characteristics between the musical abilities of humans and other animal species is still taking its first steps. One of the most promising aspects from a comparative point of view is the analysis of rhythmic components, which are crucial features of human communicative performance but also well-identifiable patterns in the vocal displays of other species. Therefore, the study of rhythm is becoming essential to understand the mechanisms of singing behavior and the evolution of human communication. Recent findings provided evidence that particular rhythmic structures occur in human music and some singing animal species, such as birds and rock hyraxes, but only 2 species of nonhuman primates have been investigated so far (Indri indri and Hylobates lar). Therefore, our study aims to consistently broaden the list of species studied regarding the presence of rhythmic categories. We investigated the temporal organization in the singing of 3 species of crested gibbons (Nomascus gabriellae, Nomascus leucogenys, and Nomascus siki) and found that the most prominent rhythmic category was isochrony. Moreover, we found slight variation in songs' tempo among species, with N. gabriellae and N. siki singing with a temporal pattern involving a gradually increasing tempo (a musical accelerando), and N. leucogenys with a more regular pattern. Here, we show how the prominence of a peak at the isochrony establishes itself as a shared characteristic in the small apes considered so far.

The variation and evolution of complete human centromeres

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

Spontaneous playful teasing in four great ape species

Joking draws on complex cognitive abilities: understanding social norms, theory of mind, anticipating others' responses and appreciating the violation of others' expectations. Playful teasing, which is present in preverbal infants, shares many of these cognitive features. There is some evidence that great apes can tease in structurally similar ways, but no systematic study exists. We developed a coding system to identify playful teasing and applied it to video of zoo-housed great apes. All four species engaged in intentionally provocative behaviour, frequently accompanied by characteristics of play. We found playful teasing to be characterized by attention-getting, one-sidedness, response looking, repetition and elaboration/escalation. It takes place mainly in relaxed contexts, has a wide variety of forms, and differs from play in several ways (e.g. asymmetry, low rates of play signals like the playface and absence of movement-final 'holds' characteristic of intentional gestures). As playful teasing is present in all extant great ape genera, it is likely that the cognitive prerequisites for joking evolved in the hominoid lineage at least 13 million years ago.

Hypoconulid loss in cercopithecins: Functional and developmental considerations

Cercopithecins differ from papionins in lacking a M3 hypoconulid. Although this loss may be related to dietary differences, the functional and developmental ramifications of hypoconulid loss are currently unclear. The following makes use of dental topographic analysis to quantify shape variation in a sample of cercopithecin M3s, as well as in a sample of Macaca, which has a hypoconulid. To help understand the consequences of hypoconulid loss, Macaca M3s were virtually cropped to remove the hypoconulid and were also subjected to dental topographic analysis. The patterning cascade model and the inhibitory cascade model attempt to explain variation in cusp pattern and molar proportions, respectively. These models have both previously been used to explain patterns of variation in cercopithecines, but have not been examined in the context of hypoconulid loss. For example, previous work suggests that earlier developing cusps impact the development of later developing cusps (i.e., the hypoconulid) and that cercopithecines do not conform to the predictions of the inhibitory cascade model in that the size of the molars is not linear moving distally. Results of the current study suggest that the loss of the hypoconulid is associated with a reduction in dental topography among cercopithecins, which is potentially related to diet, although the connection to diet is not necessarily clear. Results also suggest that the loss of the hypoconulid can be explained by the patterning cascade model, and that hypoconulid loss explains the apparent lack of support for the inhibitory cascade model among cercopithecines. These findings highlight the importance of a holistic approach to studying variation in molar proportions and developmental models.

Genetic Diversity and Sequence Conservation of Peptide-Binding Regions of MHC Class I Genes in Pig, Cattle, Chimpanzee, and Human

Comparative analyses of MHC gene diversity and evolution across different species could offer valuable insights into the evolution of MHC genes. Intra- and inter-species sequence diversity and conservation of 12 classical major histocompatibility complex (MHC) class I genes from cattle, chimpanzees, pigs, and humans was analyzed using 20 representative allelic groups for each gene. The combined analysis of paralogous loci for each species revealed that intra-locus amino-acid sequence variations in the peptide-binding region (PBR) of MHC I genes did not differ significantly between species, ranging from 8.44% for SLA to 10.75% for BoLA class I genes. In contrast, intraspecies differences in the non-PBRs of these paralogous genes were more pronounced, varying from 4.59% for SLA to 16.89% for HLA. Interestingly, the Shannon diversity index and rate of nonsynonymous substitutions for PBR were significantly higher in SLA and BoLA than those in Patr and HLA. Analysis of peptide-binding pockets across all analyzed MHC class I genes of the four species indicated that pockets A and E showed the lowest and highest diversity, respectively. The estimated divergence times suggest that primate and artiodactyl MHC class I genes diverged 60.41 Mya, and BoLA and SLA genes diverged 35.34 Mya. These results offer new insights into the conservation and diversity of MHC class I genes in various mammalian species.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Transcript Isoform Diversity of Ampliconic Genes on the Y Chromosome of Great Apes

Y chromosomal ampliconic genes (YAGs) are important for male fertility, as they encode proteins functioning in spermatogenesis. The variation in copy number and expression levels of these multicopy gene families has been studied in great apes; however, the diversity of splicing variants remains unexplored. Here, we deciphered the sequences of polyadenylated transcripts of all nine YAG families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) from testis samples of six great ape species (human, chimpanzee, bonobo, gorilla, Bornean orangutan, and Sumatran orangutan). To achieve this, we enriched YAG transcripts with capture probe hybridization and sequenced them with long (Pacific Biosciences) reads. Our analysis of this data set resulted in several findings. First, we observed evolutionarily conserved alternative splicing patterns for most YAG families except for BPY2 and PRY. Second, our results suggest that BPY2 transcripts and proteins originate from separate genomic regions in bonobo versus human, which is possibly facilitated by acquiring new promoters. Third, our analysis indicates that the PRY gene family, having the highest representation of noncoding transcripts, has been undergoing pseudogenization. Fourth, we have not detected signatures of selection in the five YAG families shared among great apes, even though we identified many species-specific protein-coding transcripts. Fifth, we predicted consensus disorder regions across most gene families and species, which could be used for future investigations of male infertility. Overall, our work illuminates the YAG isoform landscape and provides a genomic resource for future functional studies focusing on infertility phenotypes in humans and critically endangered great apes.

The promise of great apes as model organisms for understanding the downstream consequences of early life experiences

Early life experiences have a significant influence on adult health and aging processes in humans. Despite widespread interest in the evolutionary roots of this phenomenon, very little research on this topic has been conducted in humans' closest living relatives, the great apes. The longitudinal data sets that are now available on wild and captive great ape populations hold great promise to clarify the nature, evolutionary function, and mechanisms underlying these connections in species which share key human life history characteristics. Here, we explain features of great ape life history and socioecologies that make them of particular interest for this topic, as well as those that may limit their utility as comparative models; outline the ways in which available data are complementary to and extend the kinds of data that are available for humans; and review what is currently known about the connections among early life experiences, social behavior, and adult physiology and biological fitness in our closest living relatives. We conclude by highlighting key next steps for this emerging area of research.

TP53 germline pathogenic variants in modern humans were likely originated during recent human history

TP53 is crucial for maintaining genome stability and preventing oncogenesis. Germline pathogenic variation in TP53 damages its function, causing genome instability and increased cancer risk. Despite extensive study in TP53, the evolutionary origin of the human TP53 germline pathogenic variants remains largely unclear. In this study, we applied phylogenetic and archaeological approaches to identify the evolutionary origin of TP53 germline pathogenic variants in modern humans. In the phylogenic analysis, we searched 406 human TP53 germline pathogenic variants in 99 vertebrates distributed in eight clades of Primate, Euarchontoglires, Laurasiatheria, Afrotheria, Mammal, Aves, Sarcopterygii and Fish, but we observed no direct evidence for the cross-species conservation as the origin; in the archaeological analysis, we searched the variants in 5031 ancient human genomes dated between 45045 and 100 years before present, and identified 45 pathogenic variants in 62 ancient humans dated mostly within the last 8000 years; we also identified 6 pathogenic variants in 3 Neanderthals dated 44000 to 38515 years before present and 1 Denisovan dated 158 550 years before present. Our study reveals that TP53 germline pathogenic variants in modern humans were likely originated in recent human history and partially inherited from the extinct Neanderthals and Denisovans.

Principal component and linear discriminant analyses for the classification of hominoid primate specimens based on bone shape data

In this study, we tested the hypothesis that machine learning methods can accurately classify extant primates based on triquetrum shape data. We then used this classification tool to observe the affinities between extant primates and fossil hominoids. We assessed the discrimination accuracy for an unsupervised and supervised learning pipeline, i.e. with principal component analysis (PCA) and linear discriminant analysis (LDA) feature extraction, when tasked with the classification of extant primates. The trained algorithm is used to classify a sample of known fossil hominoids. For the visualization, PCA and uniform manifold approximation and projection (UMAP) are used. The results show that the discriminant function correctly classified the extant specimens with an F1-score of 0.90 for both PCA and LDA. In addition, the classification of fossil hominoids reflects taxonomy and locomotor behaviour reported in literature. This classification based on shape data using PCA and LDA is a powerful tool that can discriminate between the triquetrum shape of extant primates with high accuracy and quantitatively compare fossil and extant morphology. It can be used to support taxonomic differentiation and aid the further interpretation of fossil remains. Further testing is necessary by including other bones and more species and specimens per species extinct primates.

Comparison of antiviral responses in two bat species reveals conserved and divergent innate immune pathways

Bats host a range of disease-causing viruses without displaying clinical symptoms. The mechanisms behind this are a continuous source of interest. Here, we studied the antiviral response in the Egyptian fruit bat and Kuhl's pipistrelle, representing two subordinal clades. We profiled the antiviral response in fibroblasts using RNA sequencing and compared bat with primate and rodent responses. Both bats upregulate similar genes; however, a subset of these genes is transcriptionally divergent between them. These divergent genes also evolve rapidly in sequence, have specific promoter architectures, and are associated with programs underlying tolerance and resistance. Finally, we characterized antiviral genes that expanded in bats, with duplicates diverging in sequence and expression. Our study reveals a largely conserved antiviral program across bats and points to a set of genes that rapidly evolve through multiple mechanisms. These can contribute to bat adaptation to viral infection and provide directions to understanding the mechanisms behind it.

Comparative landscape of genetic dependencies in human and chimpanzee stem cells

Comparative studies of great apes provide a window into our evolutionary past, but the extent and identity of cellular differences that emerged during hominin evolution remain largely unexplored. We established a comparative loss-of-function approach to evaluate whether human cells exhibit distinct genetic dependencies. By performing genome-wide CRISPR interference screens in human and chimpanzee pluripotent stem cells, we identified 75 genes with species-specific effects on cellular proliferation. These genes comprised coherent processes, including cell-cycle progression and lysosomal signaling, which we determined to be human-derived by comparison with orangutan cells. Human-specific robustness to CDK2 and CCNE1 depletion persisted in neural progenitor cells and cerebral organoids, supporting the G1-phase length hypothesis as a potential evolutionary mechanism in human brain expansion. Our findings demonstrate that evolutionary changes in human cells reshaped the landscape of essential genes and establish a platform for systematically uncovering latent cellular and molecular differences between species.

Comparative landscape of genetic dependencies in human and chimpanzee stem cells

Comparative studies of great apes provide a window into our evolutionary past, but the extent and identity of cellular differences that emerged during hominin evolution remain largely unexplored. We established a comparative loss-of-function approach to evaluate whether changes in human cells alter requirements for essential genes. By performing genome-wide CRISPR interference screens in human and chimpanzee pluripotent stem cells, we identified 75 genes with species-specific effects on cellular proliferation. These genes comprised coherent processes, including cell cycle progression and lysosomal signaling, which we determined to be human-derived by comparison with orangutan cells. Human-specific robustness to CDK2 and CCNE1 depletion persisted in neural progenitor cells, providing support for the G1-phase length hypothesis as a potential evolutionary mechanism in human brain expansion. Our findings demonstrate that evolutionary changes in human cells can reshape the landscape of essential genes and establish a platform for systematically uncovering latent cellular and molecular differences between species.

Transcript Isoform Diversity of Ampliconic Genes on the Y Chromosome of Great Apes

Y-chromosomal Ampliconic Genes (YAGs) are important for male fertility, as they encode proteins functioning in spermatogenesis. The variation in copy number and expression levels of these multicopy gene families has been recently studied in great apes, however, the diversity of splicing variants remains unexplored. Here we deciphered the sequences of polyadenylated transcripts of all nine YAG families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) from testis samples of six great ape species (human, chimpanzee, bonobo, gorilla, Bornean orangutan, and Sumatran orangutan). To achieve this, we enriched YAG transcripts with capture-probe hybridization and sequenced them with long (Pacific Biosciences) reads. Our analysis of this dataset resulted in several findings. First, we uncovered a high diversity of YAG transcripts across great apes. Second, we observed evolutionarily conserved alternative splicing patterns for most YAG families except for BPY2 and PRY. Our results suggest that BPY2 transcripts and predicted proteins in several great ape species (bonobo and the two orangutans) have independent evolutionary origins and are not homologous to human reference transcripts and proteins. In contrast, our results suggest that the PRY gene family, having the highest representation of transcripts without open reading frames, has been undergoing pseudogenization. Third, even though we have identified many species-specific protein-coding YAG transcripts, we have not detected any signatures of positive selection. Overall, our work illuminates the YAG isoform landscape and its evolutionary history, and provides a genomic resource for future functional studies focusing on infertility phenotypes in humans and critically endangered great apes.

Puberty initiates a unique stage of social learning and development prior to adulthood: Insights from studies of adolescence in wild chimpanzees

In humans, puberty initiates a period of rapid growth, change, and formative neurobehavioral development. Brain and behavior changes during this maturational window contribute to opportunities for social learning. Here we provide new insights into adolescence as a unique period of social learning and development by describing field studies of our closest living relatives, chimpanzees. Like humans, chimpanzees have a multiyear juvenile life stage between weaning and puberty onset followed by a multiyear adolescent life stage after pubertal onset but prior to socially-recognized adulthood. As they develop increasing autonomy from caregivers, adolescent chimpanzees explore and develop many different types of social relationships with a wide range of individuals in a highly flexible social environment. We describe how adolescent social motivations and experiences differ from those of juveniles and adults and expose adolescents to high levels of uncertainty, risk, and vulnerability, as well as opportunities for adaptive social learning. We discuss how these adolescent learning experiences may be shaped by early life and in turn shape varied adult social outcomes. We outline how future chimpanzee field research can contribute in new ways to a more integrative interdisciplinary understanding of adolescence as a developmental window of adaptive social learning and resilience.

Long-read assembly of major histocompatibility complex and killer cell immunoglobulin-like receptor genome regions in cynomolgus macaque

BACKGROUND

The major histocompatibility complex (MHC) and the killer cell immunoglobulin-like receptors (KIR) are key regulators of immune responses. The cynomolgus macaque, an Old World monkey species, can be applied as an important preclinical model for studying human diseases, including coronavirus disease 2019 (COVID-19). Several MHC-KIR combinations have been associated with either a poor or good prognosis. Therefore, macaques with a well-characterized immunogenetic profile may improve drug evaluation and speed up vaccine development. At present, a complete overview of the MHC and KIR haplotype organizations in cynomolgus macaques is lacking, and characterization by conventional techniques is hampered by the extensive expansion of the macaque MHC-B region that complicates the discrimination between genes and alleles.

METHODS

We assembled complete MHC and KIR genomic regions of cynomolgus macaque using third-generation long-read sequencing approach. We identified functional Mafa-B loci at the transcriptome level using locus-specific amplification in a cohort of 33 Vietnamese cynomolgus macaques.

RESULTS

This is the first physical mapping of complete MHC and KIR gene regions in a Vietnamese cynomolgus macaque. Furthermore, we identified four functional Mafa-B loci (B2, B3, B5, and B6) and showed that alleles of the Mafa-I*01, -B*056, -B*034, and -B*001 functional lineages, respectively, are highly frequent in the Vietnamese cynomolgus macaque population.

CONCLUSION

The insights into the MHC and KIR haplotype organizations and the level of diversity may refine the selection of animals with specific genetic markers for future medical research.

Dynamic evolution of Mhc haplotypes in cynomolgus macaques of different geographic origins

The major histocompatibility complex (MHC) plays a key role in immune defense, and the Mhc genes of cynomolgus macaque display a high degree of polymorphism. Based on their geographic distribution, different populations of cynomolgus macaques are recognized. Here we present the characterization of the Mhc class I and II repertoire of a large pedigreed group of cynomolgus macaques originating from the mainland north of the isthmus of Kra (N = 42). Segregation analyses resulted in the definition of 81 unreported Mafa-A/B/DRB/DQ/DP haplotypes, which include 32 previously unknown DRB regions. In addition, we report 13 newly defined Mafa-A/B/DRB/DQ/DP haplotypes in a group of cynomolgus macaques originating from the mainland south of the isthmus of Kra/Maritime Southeast Asia (N = 16). A relatively high level of sharing of Mafa-A (51%) and Mafa-B (40%) lineage groups is observed between the populations native to the north and the south of isthmus of Kra. At the allelic level, however, the Mafa-A/B haplotypes seem to be characteristic of a population. An overall comparison of all currently known data revealed that each geographic population has its own specific combinations of Mhc class I and II haplotypes. This illustrates the dynamic evolution of the cynomolgus macaque Mhc region, which was most likely generated by recombination and maintained by selection due to the differential pathogenic pressures encountered in different geographic areas.

Impact of natural selection on global patterns of genetic variation and association with clinical phenotypes at genes involved in SARS-CoV-2 infection

Human genomic diversity has been shaped by both ancient and ongoing challenges from viruses. The current coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a devastating impact on population health. However, genetic diversity and evolutionary forces impacting host genes related to SARS-CoV-2 infection are not well understood. We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (angiotensin converting enzyme 2 [ACE2], transmembrane protease serine 2 [TMPRSS2], dipeptidyl peptidase 4 [DPP4], and lymphocyte antigen 6 complex locus E [LY6E]). We analyzed data from 2,012 ethnically diverse Africans and 15,977 individuals of European and African ancestry with electronic health records and integrated with global data from the 1000 Genomes Project. At ACE2, we identified 41 nonsynonymous variants that were rare in most populations, several of which impact protein function. However, three nonsynonymous variants (rs138390800, rs147311723, and rs145437639) were common among central African hunter-gatherers from Cameroon (minor allele frequency 0.083 to 0.164) and are on haplotypes that exhibit signatures of positive selection. We identify signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage compared with the chimpanzee genome. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19. Our study provides insights into global variation at host genes related to SARS-CoV-2 infection, which have been shaped by natural selection in some populations, possibly due to prior viral infections.

Cryptic Species Exist in Vietnamella sinensis Hsu, 1936 (Insecta: Ephemeroptera) from Studies of Complete Mitochondrial Genomes

Ephemeroptera (Insecta: Pterygota) are widely distributed all over the world with more than 3500 species. During the last decade, the phylogenetic relationships within Ephemeroptera have been a hot topic of research, especially regarding the phylogenetic relationships among Vietnamellidae. In this study, three mitochondrial genomes from three populations of Vienamella sinensis collected from Tonglu (V. sinensis TL), Chun’an (V. sinensis CN), and Qingyuan (V. sinensis QY) in Zhejiang Province, China were compared to discuss the potential existence of cryptic species. We also established their phylogenetic relationship by combining the mt genomes of 69 Ephemeroptera downloaded from NCBI. The mt genomes of V. sinensis TL, V. sinensis CN, and V. sinensis QY showed the same gene arrangement with lengths of 15,674 bp, 15,674 bp, and 15,610 bp, respectively. Comprehensive analyses of these three mt genomes revealed significant differences in mt genome organization, genetic distance, and divergence time. Our results showed that the specimens collected from Chun’an and Tonglu in Zhejiang Province, China belonged to V. sinensis, and the specimens collected from Qingyuan, Zhejiang Province, China were a cryptic species of V. sinensis. In maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees, the monophyly of the family Vietnamellidae was supported and Vietnamellidae has a close relationship with Ephemerellidae.

A Comprehensive Atlas of Immunological Differences Between Humans, Mice, and Non-Human Primates

Animal models are an integral part of the drug development and evaluation process. However, they are unsurprisingly imperfect reflections of humans, and the extent and nature of many immunological differences are unknown. With the rise of targeted and biological therapeutics, it is increasingly important that we understand the molecular differences in the immunological behavior of humans and model organisms. However, very few antibodies are raised against non-human primate antigens, and databases of cross-reactivity between species are incomplete. Thus, we screened 332 antibodies in five immune cell populations in blood from humans and four non-human primate species generating a comprehensive cross-reactivity catalog that includes cell type-specificity. We used this catalog to create large mass cytometry universal cross-species phenotyping and signaling panels for humans, along with three of the model organisms most similar to humans: rhesus and cynomolgus macaques and African green monkeys; and one of the mammalian models most widely used in drug development: C57BL/6 mice. As a proof-of-principle, we measured immune cell signaling responses across all five species to an array of 15 stimuli using mass cytometry. We found numerous instances of different cellular phenotypes and immune signaling events occurring within and between species, and detailed three examples (double-positive T cell frequency and signaling; granulocyte response to Bacillus anthracis antigen; and B cell subsets). We also explore the correlation of herpes simian B virus serostatus on the immune profile. Antibody panels and the full dataset generated are available online as a resource to enable future studies comparing immune responses across species during the evaluation of therapeutics.

Bifidobacterium castoris strains isolated from wild mice show evidence of frequent host switching and diverse carbohydrate metabolism potential

Members of the gut microbiota genus Bifidobacterium are widely distributed human and animal symbionts believed to exert beneficial effects on their hosts. However, in-depth genomic analyses of animal-associated species and strains are somewhat lacking, particularly in wild animal populations. Here, to examine patterns of host specificity and carbohydrate metabolism capacity, we sequenced whole genomes of Bifidobacterium isolated from wild-caught small mammals from two European countries (UK and Lithuania). Members of Bifidobacterium castoris, Bifidobacterium animalis and Bifodobacterium pseudolongum were detected in wild mice (Apodemus sylvaticus, Apodemus agrarius and Apodemus flavicollis), but not voles or shrews. B. castoris constituted the most commonly recovered Bifidobacterium (78% of all isolates), with the majority of strains only detected in a single population, although populations frequently harboured multiple co-circulating strains. Phylogenetic analysis revealed that the mouse-associated B. castoris clades were not specific to a particular location or host species, and their distribution across the host phylogeny was consistent with regular host shifts rather than host-microbe codiversification. Functional analysis, including in vitro growth assays, suggested that mouse-derived B. castoris strains encoded an extensive arsenal of carbohydrate-active enzymes, including putative novel glycosyl hydrolases such as chitosanases, along with genes encoding putative exopolysaccharides, some of which may have been acquired via horizontal gene transfer. Overall, these results provide a rare genome-level analysis of host specificity and genomic capacity among important gut symbionts of wild animals, and reveal that Bifidobacterium has a labile relationship with its host over evolutionary time scales.

The evolving story of apolipoprotein L1 nephropathy: the end of the beginning

Genetic coding variants in APOL1, which encodes apolipoprotein L1 (APOL1), were identified in 2010 and are relatively common among individuals of sub-Saharan African ancestry. Approximately 13% of African Americans carry two APOL1 risk alleles. These variants, termed G1 and G2, are a frequent cause of kidney disease — termed APOL1 nephropathy — that typically manifests as focal segmental glomerulosclerosis and the clinical syndrome of hypertension and arterionephrosclerosis. Cell culture studies suggest that APOL1 variants cause cell dysfunction through several processes, including alterations in cation channel activity, inflammasome activation, increased endoplasmic reticulum stress, activation of protein kinase R, mitochondrial dysfunction and disruption of APOL1 ubiquitinylation. Risk of APOL1 nephropathy is mostly confined to individuals with two APOL1 risk variants. However, only a minority of individuals with two APOL1 risk alleles develop kidney disease, suggesting the need for a ‘second hit’. The best recognized factor responsible for this ‘second hit’ is a chronic viral infection, particularly HIV-1, resulting in interferon-mediated activation of the APOL1 promoter, although most individuals with APOL1 nephropathy do not have an obvious cofactor. Current therapies for APOL1 nephropathies are not adequate to halt progression of chronic kidney disease, and new targeted molecular therapies are in clinical trials.

This Review summarizes current understanding of the role of APOL1 variants in kidney disease. The authors discuss the genetics, protein structure and biological functions of APOL1 variants and provide an overview of promising therapeutic strategies.

In contrast to other APOL family members, which are primarily intracellular, APOL1 contains a unique secretory signal peptide, resulting in its secretion into plasma.

APOL1 renal risk alleles provide protection from African human trypanosomiasis but are a risk factor for progressive kidney disease in those carrying two risk alleles.

APOL1 risk allele frequency is ~35% in the African American population in the United States, with ~13% of individuals having two risk alleles; the highest allele frequencies are found in West African populations and their descendants.

Cell and mouse models implicate endolysosomal and mitochondrial dysfunction, altered ion channel activity, altered autophagy, and activation of protein kinase R in the pathogenesis of APOL1-associated kidney disease; however, the relevance of these injury pathways to human disease has not been resolved.

APOL1 kidney disease tends to be progressive, and current standard therapies are generally ineffective; targeted therapeutic strategies hold the most promise.

The Genomic Organization of the LILR Region Remained Largely Conserved Throughout Primate Evolution: Implications for Health And Disease

The genes of the leukocyte immunoglobulin-like receptor (LILR) family map to the leukocyte receptor complex (LRC) on chromosome 19, and consist of both activating and inhibiting entities. These receptors are often involved in regulating immune responses, and are considered to play a role in health and disease. The human LILR region and evolutionary equivalents in some rodent and bird species have been thoroughly characterized. In non-human primates, the LILR region is annotated, but a thorough comparison between humans and non-human primates has not yet been documented. Therefore, it was decided to undertake a comprehensive comparison of the human and non-human primate LILR region at the genomic level. During primate evolution the organization of the LILR region remained largely conserved. One major exception, however, is provided by the common marmoset, a New World monkey species, which seems to feature a substantial contraction of the number of LILR genes in both the centromeric and the telomeric region. Furthermore, genomic analysis revealed that the killer-cell immunoglobulin-like receptor gene KIR3DX1, which maps in the LILR region, features one copy in humans and great ape species. A second copy, which might have been introduced by a duplication event, was observed in the lesser apes, and in Old and New World monkey species. The highly conserved gene organization allowed us to standardize the LILR gene nomenclature for non-human primate species, and implies that most of the receptors encoded by these genes likely fulfill highly preserved functions.

Factors Regulating the Activity of LINE1 Retrotransposons

LINE-1 (L1) is a class of autonomous mobile genetic elements that form somatic mosaicisms in various tissues of the organism. The activity of L1 retrotransposons is strictly controlled by many factors in somatic and germ cells at all stages of ontogenesis. Alteration of L1 activity was noted in a number of diseases: in neuropsychiatric and autoimmune diseases, as well as in various forms of cancer. Altered activity of L1 retrotransposons for some pathologies is associated with epigenetic changes and defects in the genes involved in their repression. This review discusses the molecular genetic mechanisms of the retrotransposition and regulation of the activity of L1 elements. The contribution of various factors controlling the expression and distribution of L1 elements in the genome occurs at all stages of the retrotransposition. The regulation of L1 elements at the transcriptional, post-transcriptional and integration into the genome stages is described in detail. Finally, this review also focuses on the evolutionary aspects of L1 accumulation and their interplay with the host regulation system.

Impact of natural selection on global patterns of genetic variation, and association with clinical phenotypes, at genes involved in SARS-CoV-2 infection

We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (ACE2, TMPRSS2, DPP4, and LY6E). We analyzed novel data from 2,012 ethnically diverse Africans and 15,997 individuals of European and African ancestry with electronic health records, and integrated with global data from the 1000GP. At ACE2, we identified 41 non-synonymous variants that were rare in most populations, several of which impact protein function. However, three non-synonymous variants were common among Central African hunter-gatherers from Cameroon and are on haplotypes that exhibit signatures of positive selection. We identify strong signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19.

Impact of natural selection on global patterns of genetic variation, and association with clinical phenotypes, at genes involved in SARS-CoV-2 infection

We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection ( ACE2, TMPRSS2, DPP4 , and LY6E ). We analyzed novel data from 2,012 ethnically diverse Africans and 15,997 individuals of European and African ancestry with electronic health records, and integrated with global data from the 1000GP. At ACE2 , we identified 41 non-synonymous variants that were rare in most populations, several of which impact protein function. However, three non-synonymous variants were common among Central African hunter-gatherers from Cameroon and are on haplotypes that exhibit signatures of positive selection. We identify strong signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2 , we identified 13 amino acid changes that are adaptive and specific to the human lineage. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19.

Ampliconic Genes on the Great Ape Y Chromosomes: Rapid Evolution of Copy Number but Conservation of Expression Levels

Multicopy ampliconic gene families on the Y chromosome play an important role in spermatogenesis. Thus, studying their genetic variation in endangered great ape species is critical. We estimated the sizes (copy number) of nine Y ampliconic gene families in population samples of chimpanzee, bonobo, and orangutan with droplet digital polymerase chain reaction, combined these estimates with published data for human and gorilla, and produced genome-wide testis gene expression data for great apes. Analyzing this comprehensive data set within an evolutionary framework, we, first, found high inter- and intraspecific variation in gene family size, with larger families exhibiting higher variation as compared with smaller families, a pattern consistent with random genetic drift. Second, for four gene families, we observed significant interspecific size differences, sometimes even between sister species—chimpanzee and bonobo. Third, despite substantial variation in copy number, Y ampliconic gene families’ expression levels did not differ significantly among species, suggesting dosage regulation. Fourth, for three gene families, size was positively correlated with gene expression levels across species, suggesting that, given sufficient evolutionary time, copy number influences gene expression. Our results indicate high variability in size but conservation in gene expression levels in Y ampliconic gene families, significantly advancing our understanding of Y-chromosome evolution in great apes.

Does play shape hand use skill in rats?

Hand use is a widespread act in many vertebrate lineages and subserves behaviors including locomotion, predation, feeding, nest construction, and grooming. In order to determine whether hand use is similarly used in social behavior, the present paper describes hand use in the social play of rats. In the course of rough and tumble play sessions, rats are found to make as many as twenty different movements a minute with each hand for the purposes of manipulating a partner into a subordinate position or defending against a partner's attack. The hand movements comprise signaling movements of touching, offensive manipulating of a partner to control a play engagement, and defensive hand movements directed toward blocking, pushing and pulling to parry an attack. For signaling, attack and defense, hand movements have a structure that is similar to the structure of hand movements used for other purposes including eating, but in their contact points on an opponent, they are tailored for partner control. Given the time devoted to play by rats, play likely features the social rat behavior with the most extensive use of hand movements. This extensive use of hand movements for social play is discussed in relation to the ubiquity of hand use in adaptive behavior, the evolution of hand use in the play of mammals, and in relation to extending the multifunctional theory of the purposes of play to include the education of skilled hand movements for various adult functions including as feeding.

Two Human Monoclonal HLA-Reactive Antibodies Cross-React with Mamu-B*008, a Rhesus Macaque MHC Allotype Associated with Control of Simian Immunodeficiency Virus Replication

MHC class I molecules play an important role in adaptive immune responses against intracellular pathogens. These molecules are highly polymorphic, and many allotypes have been characterized. In a transplantation setting, a mismatch between MHC allotypes may initiate an alloimmune response. Rhesus macaques (Macaca mulatta, Mamu) are valuable as a preclinical model species in transplantation research as well as to evaluate the safety and efficacy of vaccine candidates. In both lines of research, the availability of nonhuman primate MHC-reactive mAbs may enable in vitro monitoring and detection of presence of particular Mamu molecules. In this study, we screened a collection of thoroughly characterized HLA class I-specific human mAbs for cross-reactivity with rhesus macaque MHC class I allotypes. Two mAbs, OK4F9 and OK4F10, recognize an epitope that is defined by isoleucine (I) at amino acid position 142 that is present on the Indian rhesus macaque Mamu-B*008:01 allotype, which is an allotype known to be associated with elite control of SIV replication. The reactive pattern of a third mAb, MUS4H4, is more complex and includes an epitope shared on Mamu-A2*05:01 and -B*001:01-encoded Ags. This is the first description, to our knowledge, of human HLA-reactive mAbs that can recognize Mamu allotypes, and these can be useful tools for in vitro monitoring the presence of the relevant allelic products. Moreover, OK4F9 and OK4F10 can be powerful mAbs for application in SIV-related research.

Linguistic laws of brevity: conformity in Indri indri

Vocal and gestural sequences of several primates have been found to conform to two general principles of information compression: the compensation between the duration of a construct and that of its components (Menzerath-Altmann law) and an inverse relationship between signal duration and its occurrence (Zipf's law of abbreviation). Even though Zipf's law of brevity has been proposed as a universal in animal communication, evidence on non-human primate vocal behavior conformity to linguistic laws is still debated, and information on strepsirrhine primates is lacking. We analyzed the vocal behavior of the unique singing lemur species (Indri indri) to assess whether the song of the species shows evidence for compression. As roars have a chaotic structure that impedes the recognition of each individual utterance, and long notes are usually given by males, we focused on the core part of the song (i.e., the descending phrases, composed of two-six units). Our results indicate that indris' songs conform to Zipf's and Menzerath-Altmann linguistic laws. Indeed, shorter phrases are more likely to be included in the song, and units' duration decrease at the increase of the size of the phrases. We also found that, despite a sexual dimorphism in the duration of both units and phrases, these laws characterize sequences of both males and females. Overall, we provide the first evidence for a trade-off between signal duration and occurrence in the vocal behavior of a strepsirrhine species, suggesting that selective pressures for vocal compression are more ancestral than previously assumed within primates.

The immune-modulating pregnancy-specific glycoproteins evolve rapidly and their presence correlates with hemochorial placentation in primates

Background

Pregnancy-specific glycoprotein (PSG) genes belong to the carcinoembryonic antigen (CEA) gene family, within the immunoglobulin gene superfamily. In humans, 10 PSG genes encode closely related secreted glycoproteins. They are exclusively expressed in fetal syncytiotrophoblast cells and represent the most abundant fetal proteins in the maternal blood. In recent years, a role in modulation of the maternal immune system possibly to avoid rejection of the semiallogeneic fetus and to facilitate access of trophoblast cells to maternal resources via the blood system has been suggested. Alternatively, they could serve as soluble pathogen decoy receptors like other members of the CEA family. Despite their clearly different domain organization, similar functional properties have also been observed for murine and bat PSG. As these species share a hemochorial type of placentation and a seemingly convergent formation of PSG genes during evolution, we hypothesized that hemochorial placentae support the evolution of PSG gene families.

Results

To strengthen this hypothesis, we have analyzed PSG genes in 57 primate species which exhibit hemochorial or epitheliochorial placentation. In nearly all analyzed apes some 10 PSG genes each could be retrieved from genomic databases, while 6 to 24 PSG genes were found in Old World monkey genomes. Surprisingly, only 1 to 7 PSG genes could be identified in New World monkeys. Interestingly, no PSG genes were found in more distantly related primates with epitheliochorial placentae like lemurs and lorises. The exons encoding the putative receptor-binding domains exhibit strong selection for diversification in most primate PSG as revealed by rapid loss of orthologous relationship during evolution and high ratios of nonsynonymous and synonymous mutations.

Conclusion

The distribution of trophoblast-specific PSGs in primates and their pattern of selection supports the hypothesis that PSG are still evolving to optimize fetal-maternal or putative pathogen interactions in mammals with intimate contact of fetal cells with the immune system of the mother like in hemochorial placentation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07413-8.

Dynamic evolution of great ape Y chromosomes

The male-specific Y chromosome harbors genes important for sperm production. Because Y is repetitive, its DNA sequence was deciphered for only a few species, and its evolution remains elusive. Here we compared the Y chromosomes of great apes (human, chimpanzee, bonobo, gorilla, and orangutan) and found that many of their repetitive sequences and multicopy genes were likely already present in their common ancestor. Y repeats had increased intrachromosomal contacts, which might facilitate preservation of genes and gene regulatory elements. Chimpanzee and bonobo, experiencing high sperm competition, underwent many DNA changes and gene losses on the Y. Our research is significant for understanding the role of the Y chromosome in reproduction of nonhuman great apes, all of which are endangered.

The mammalian male-specific Y chromosome plays a critical role in sex determination and male fertility. However, because of its repetitive and haploid nature, it is frequently absent from genome assemblies and remains enigmatic. The Y chromosomes of great apes represent a particular puzzle: their gene content is more similar between human and gorilla than between human and chimpanzee, even though human and chimpanzee share a more recent common ancestor. To solve this puzzle, here we constructed a dataset including Ys from all extant great ape genera. We generated assemblies of bonobo and orangutan Ys from short and long sequencing reads and aligned them with the publicly available human, chimpanzee, and gorilla Y assemblies. Analyzing this dataset, we found that the genus Pan, which includes chimpanzee and bonobo, experienced accelerated substitution rates. Pan also exhibited elevated gene death rates. These observations are consistent with high levels of sperm competition in Pan. Furthermore, we inferred that the great ape common ancestor already possessed multicopy sequences homologous to most human and chimpanzee palindromes. Nonetheless, each species also acquired distinct ampliconic sequences. We also detected increased chromatin contacts between and within palindromes (from Hi-C data), likely facilitating gene conversion and structural rearrangements. Our results highlight the dynamic mode of Y chromosome evolution and open avenues for studies of male-specific dispersal in endangered great ape species.

A bigger brain for a more complex environment

The environment increased complexity required more neural functions to develop in the hominin brains, and the hominins adapted to the complexity by developing a bigger brain with a greater interconnection between its parts. Thus, complex environments drove the growth of the brain. In about two million years during hominin evolution, the brain increased three folds in size, one of the largest and most complex amongst mammals, relative to body size. The size increase has led to anatomical reorganization and complex neuronal interactions in a relatively small skull. At birth, the human brain is only about 20% of its adult size. That facilitates the passage through the birth canal. Therefore, the human brain, especially cortex, develops postnatally in a rich stimulating environment with continuous brain wiring and rewiring and insertion of billions of new neurons. One of the consequence is that in the newborn brain, neuroplasticity is always turned "on" and it remains active throughout life, which gave humans the ability to adapt to complex and often hostile environments, integrate external experiences, solve problems, elaborate abstract ideas and innovative technologies, store a lot of information. Besides, hominins acquired unique abilities as music, language, and intense social cooperation. Overwhelming ecological, social, and cultural challenges have made the human brain so unique. From these events, as well as the molecular genetic changes that took place in those million years, under the pressure of natural selection, derive the distinctive cognitive abilities that have led us to complex social organizations and made our species successful.

Comparative genetics of the major histocompatibility complex in humans and nonhuman primates

The major histocompatibility complex (MHC) is one of the most gene-dense regions of the mammalian genome. Multiple genes within the human MHC (HLA) show extensive polymorphism, and currently, more than 26,000 alleles divided over 39 different genes are known. Nonhuman primate (NHP) species are grouped into great and lesser apes and Old and New World monkeys, and their MHC is studied mostly because of their important role as animal models in preclinical research or in connection with conservation biology purposes. The evolutionary equivalents of many of the HLA genes are present in NHP species, and these genes may also show abundant levels of polymorphism. This review is intended to provide a comprehensive comparison relating to the organization and polymorphism of human and NHP MHC regions.

Non-Human Primate Models of Dengue Virus Infection: A Comparison of Viremia Levels and Antibody Responses during Primary and Secondary Infection among Old World and New World Monkeys

Due to the global burden of dengue disease, a vaccine is urgently needed. One of the key points in vaccine development is the development of a robust and reliable animal model of dengue virus infection. Characteristics including the ability to sustain viral replication, demonstration of clinical signs, and immune response that resemble those of human dengue virus infection are vital in animal models. Preclinical studies in vaccine development usually include parameters such as safety evaluation, induction of viremia and antigenemia, immunogenicity, and vaccine effectiveness. Although mice have been used as a model, non-human primates have an advantage over mice because of their relative similarity to humans in their genetic composition and immune responses. This review compares the viremia kinetics and antibody responses of cynomolgus macaques (Macaca fasicularis), common marmosets (Callithrix jacchus), and tamarins (Saguinus midas and Saguinus labitus) and summarize the perspectives and the usefulness along with challenges in dengue vaccine development.

The draft genome of mandrill (Mandrillus sphinx): An Old World monkey

Mandrill (Mandrillus sphinx) is a primate species, which belongs to the Old World monkey (Cercopithecidae) family. It is closely related to human, serving as a model for human health related research. However, the genetic studies on and genomic resources of mandrill are limited, especially in comparison to other primate species. Here we produced 284 Gb data, providing 96-fold coverage (considering the estimated genome size of 2.9 Gb), to construct a reference genome for the mandrill. The assembled draft genome was 2.79 Gb with contig N50 of 20.48 Kb and scaffold N50 of 3.56 Mb. We annotated the mandrill genome to find 43.83% repeat elements, as well as 21,906 protein-coding genes. The draft genome was of good quality with 98% gene annotation coverage by Benchmarking Universal Single-Copy Orthologs (BUSCO). Based on comparative genomic analyses of  the Major Histocompatibility Complex (MHC) of the immune system in mandrill and human, we found that 17 genes in the mandrill that have been associated with disease phenotypes in human such as Lung cancer, cranial volume and asthma, barbored amino acids changing mutations. Gene family analyses revealed expansion of several genes, and several genes associated with stress environmental adaptation and innate immunity responses exhibited signatures of positive selection. In summary, we established the first draft genome of  the mandrill of value for studies on evolution and human health.

Taxonomic variation in the supraorbital region of catarrhine primates

OBJECTIVES

This study aimed to test the taxonomic utility of the catarrhine supraorbital region using 3D geometric morphometrics, with the aim of establishing its potential use in elucidating the position of more debated hominin groups.

MATERIALS AND METHODS

230 3D coordinates were used to record the supraorbital morphology of two datasets: one containing 460 non-hominin catarrhine primates from species and subspecies of Gorilla, Pan, Papio, and Macaca; and the other containing 55 Pleistocene hominins from Homo, Australopithecus, and Paranthropus. Principal component analyses in tangent, form, and allometry-free shape space were used to assess differentiation of taxa, with biological distinctiveness of taxa being established using step-wise discriminant analysis with subsampling.

RESULTS

Results indicated that the recorded supraorbital morphology could be used to separate non-hominin catarrhine primate genera, species, and subspecies, although accuracy was found to decrease with decreasing Linnaean rank. In addition, analyses in tangent space were found to produce the highest accuracy when classifying primates of known taxonomy. Biological distinctiveness of the middle and later Homo species was comparable to or higher than that of the non-hominin primates, and relatively lower for the earlier groups of Homo.

DISCUSSION

This study indicates that the supraorbital region preserves taxonomic information that can be used to delineate between closely related groups, both within hominins and wider catarrhine primates. Therefore, this region may be used to provide insight when assessing the taxonomic affiliation of disputed hominin specimens.

Cerebral blood flow rates in recent great apes are greater than in Australopithecus species that had equal or larger brains

Brain metabolic rate (MR) is linked mainly to the cost of synaptic activity, so may be a better correlate of cognitive ability than brain size alone. Among primates, the sizes of arterial foramina in recent and fossil skulls can be used to evaluate brain blood flow rate, which is proportional to brain MR. We use this approach to calculate flow rate in the internal carotid arteries (Q˙ICA), which supply most of the primate cerebrum. Q˙ICA is up to two times higher in recent gorillas, chimpanzees and orangutans compared with 3-million-year-old australopithecine human relatives, which had equal or larger brains. The scaling relationships between Q˙ICA and brain volume (V_br) show exponents of 1.03 across 44 species of living haplorhine primates and 1.41 across 12 species of fossil hominins. Thus, the evolutionary trajectory for brain perfusion is much steeper among ancestral hominins than would be predicted from living primates. Between 4.4-million-year-old Ardipithecus and Homo sapiens, V_br increased 4.7-fold, but Q˙ICA increased 9.3-fold, indicating an approximate doubling of metabolic intensity of brain tissue. By contrast, Q˙ICA is proportional to V_br among haplorhine primates, suggesting a constant volume-specific brain MR.

High Satellite Repeat Turnover in Great Apes Studied with Short- and Long-Read Technologies

Satellite repeats are a structural component of centromeres and telomeres, and in some instances, their divergence is known to drive speciation. Due to their highly repetitive nature, satellite sequences have been understudied and underrepresented in genome assemblies. To investigate their turnover in great apes, we studied satellite repeats of unit sizes up to 50 bp in human, chimpanzee, bonobo, gorilla, and Sumatran and Bornean orangutans, using unassembled short and long sequencing reads. The density of satellite repeats, as identified from accurate short reads (Illumina), varied greatly among great ape genomes. These were dominated by a handful of abundant repeated motifs, frequently shared among species, which formed two groups: 1) the (AATGG)n repeat (critical for heat shock response) and its derivatives; and 2) subtelomeric 32-mers involved in telomeric metabolism. Using the densities of abundant repeats, individuals could be classified into species. However, clustering did not reproduce the accepted species phylogeny, suggesting rapid repeat evolution. Several abundant repeats were enriched in males versus females; using Y chromosome assemblies or Fluorescent In Situ Hybridization, we validated their location on the Y. Finally, applying a novel computational tool, we identified many satellite repeats completely embedded within long Oxford Nanopore and Pacific Biosciences reads. Such repeats were up to 59 kb in length and consisted of perfect repeats interspersed with other similar sequences. Our results based on sequencing reads generated with three different technologies provide the first detailed characterization of great ape satellite repeats, and open new avenues for exploring their functions.

Analysis of differentially methylated regions in great apes and extinct hominids provides support for the evolutionary hypothesis of schizophrenia

INTRODUCTION

The persistence of schizophrenia in human populations separated by geography and time led to the evolutionary hypothesis that proposes schizophrenia as a by-product of the higher cognitive abilities of modern humans. To explore this hypothesis, we used here an evolutionary epigenetics approach building on differentially methylated regions (DMRs) of the genome.

METHODS

We implemented a polygenic enrichment testing pipeline using the summary statistics of genome-wide association studies (GWAS) of schizophrenia and 12 other phenotypes. We investigated the enrichment of association of these traits across genomic regions with variable methylation between modern humans and great apes (orangutans, chimpanzees and gorillas; great ape DMRs) and between modern humans and recently extinct hominids (Neanderthals and Denisovans; hominid DMRs).

RESULTS

Regions that are hypo-methylated in humans compared to great apes show enrichment of association with schizophrenia only if the major histocompatibility complex (MHC) region is included. With the MHC region removed from the analysis, only a modest enrichment for SNPs of low effect persists. The INRICH pipeline confirms this finding after rigorous permutation and bootstrapping procedures.

CONCLUSION

The analyses of regions with differential methylation changes in humans and great apes do not provide compelling evidence of enrichment of association with schizophrenia, in contrast to our previous findings on more recent methylation differences between modern humans, Neanderthals and Denisovans. Our results further support the evolutionary hypothesis of schizophrenia and indicate that the origin of some of the genetic susceptibility factors of schizophrenia may lie in recent human evolution.

A comprehensive analysis of the Baboon-specific full-length LINE-1 retrotransposons

BACKGROUND

Long interspersed elements-1 (LINE-1s or L1s) and Alu elements are most successful retrotransposons that have generated genetic diversity and genomic fluidity in the primate genome. They account for ~ 27.7% of the primate genome. Interestingly, a previous study has shown that the retrotransposition rate of Alu elements is nine times higher in baboons than in humans.

OBJECTIVE

The expansion of Alu copies could be dependent on the activity of L1-encoded proteins. Thus, we aimed to investigate full-length baboon-specific L1s and characterize structurally and functionally intact baboon-specific L1s (ORF1p/ORF2p and ORF2p only) that could induce trans-mobilization of Alu elements in the baboon genome.

RESULTS

A total of 673 baboon-specific L1 candidates (> 4 kb) were identified through the comparative genomic analysis. Applying the baboon-specific correction value obtained from the experimental validation, it demonstrated that approximately 446 baboon-specific L1s (> 4 kb) were present in the baboon reference genome (papAnu2). In addition, we observed phylogenetic relationship of the baboon-specific L1s through the neighbor-joining method and they diverged from the L1PA6 consensus sequence. Finally, we identified 36 full-length baboon-specific L1s that were intact both ORF1p and ORF2p.

CONCLUSION

The number of baboon-specific full-length L1s is fewer than the number of human-specific full-length L1s. Therefore, there is possibility that the "L1 master gene" or "L1 source gene" is more abundant in the baboon genome, or that in trans retrotransposition activity of baboon-specific L1s is relatively stronger than in the other genomes.

Uncovering the inferior fronto-occipital fascicle and its topological organization in non-human primates: the missing connection for language evolution

Whether brain networks underlying the multimodal processing of language in humans are present in non-human primates is an unresolved question in primate evolution. Conceptual awareness in humans, which is the backbone of verbal and non-verbal semantic elaboration, involves intracerebral connectivity via the inferior fronto-occipital fascicle (IFOF). While non-human primates can communicate through visual information channels, there has been no formal demonstration that they possess a functional homologue of the human IFOF. Therefore, we undertook a post-mortem diffusion MRI tractography study in conjunction with Klingler micro-dissection to search for IFOF fiber tracts in brain of Old-World (vervet) monkeys. We found clear and concordant evidence from both techniques for the existence of bilateral fiber tracts connecting the frontal and occipital lobes. These tracts closely resembled the human IFOF with respect to trajectory, topological organization, and cortical terminal fields. Moreover, these fibers are clearly distinct from other bundles previously described in this region of monkey brain, i.e., the inferior longitudinal and uncinate fascicles, and the external and extreme capsules. This demonstration of an IFOF in brain of a species that diverged from the human lineage some 22 millions years ago enhances our comprehension about the evolution of language and social behavior.

Predominant patterns of splicing evolution on human, chimpanzee and macaque evolutionary lineages

Although splicing is widespread and evolves rapidly among species, the mechanisms driving this evolution, as well as its functional implications, are not yet fully understood. We analyzed the evolution of splicing patterns based on transcriptome data from five tissues of humans, chimpanzees, rhesus macaques and mice. In total, 1526 exons and exon sets from 1236 genes showed significant splicing differences among primates. More than 60% of these differences represent constitutive-to-alternative exon transitions while an additional 25% represent changes in exon inclusion frequency. These two dominant evolutionary patterns have contrasting conservation, regulation and functional features. The sum of these features indicates that, despite their prevalence, constitutive-to-alternative exon transitions do not substantially contribute to long-term functional transcriptome changes. Conversely, changes in exon inclusion frequency appear to be functionally relevant, especially for changes taking place in the brain on the human evolutionary lineage.

snRPN controls the ability of neurons to regenerate axons

BACKGROUND

Retinal ganglion cells (RGCs) of mammals lose the ability to regenerate injured axons during postnatal maturation, but little is known about the underlying molecular mechanisms.

OBJECTIVE

It remains of particular importance to understand the mechanisms of axonal regeneration to develop new therapeutic approaches for nerve injuries.

METHODS

Retinas from newborn to adult monkeys (Callithrix jacchus)1 were obtained immediately after death and cultured in vitro. Growths of axons were monitored using microscopy and time-lapse video cinematography. Immunohistochemistry, Western blotting, qRT-PCR, and genomics were performed to characterize molecules associated with axonal regeneration and growth. A genomic screen was performed by using retinal explants versus native and non-regenerative explants obtained from eye cadavers on the day of birth, and hybridizing the mRNA with cross-reacting cDNA on conventional human microarrays. Followed the genomic screen, siRNA experiments were conducted to identify the functional involvement of identified candidates.

RESULTS

Neuron-specific human ribonucleoprotein N (snRPN) was found to be a potential regulator of impaired axonal regeneration during neuronal maturation in these animals. In particular, up-regulation of snRPN was observed during retinal maturation, coinciding with a decline in regenerative ability. Axon regeneration was reactivated in snRPN-knockout retinal ex vivo explants of adult monkey.

CONCLUSION

These results suggest that coordinated snRPN-driven activities within the neuron-specific ribonucleoprotein complex regulate the regenerative ability of RGCs in primates, thereby highlighting a potential new role for snRPN within neurons and the possibility of novel postinjury therapies.

Trypanosoma cruzi genetic diversity: Something new for something known about Chagas disease manifestations, serodiagnosis and drug sensitivity

The genetic diversity of Trypanosoma cruzi, the protozoan agent of Chagas disease, is widely recognized. At present, T. cruzi is partitioned into seven discrete typing units (DTUs), TcI-TcVI and Tcbat. This article reviews the present knowledge on the parasite population structure, the evolutionary relationships among DTUs and their distinct, but not exclusive ecological and epidemiological associations. Different models for the origin of hybrid DTUs are examined, which agree that genetic exchange among T. cruzi populations is frequent and has contributed to the present parasite population structure. The geographic distribution of the prevalent DTUs in humans from the southern United States to Argentina is here presented and the circumstantial evidence of a possible association between T. cruzi genotype and Chagas disease manifestations is discussed. The available information suggests that parasite strains detected in patients, regardless of the clinical presentation, reflect the principal DTU circulating in the domestic transmission cycles of a particular region. In contrast, in several orally transmitted outbreaks, sylvatic strains are implicated. As a consequence of the genotypic and phenotypic differences of T. cruzi strains and the differential geographic distribution of DTUs in humans, regional variations in the sensitivity of the serological tests are verified. The natural resistance to benznidazole and nifurtimox, verified in vivo and in vitro for some parasite stocks, is not associated with any particular DTU, and does not explain the marked difference in the anti-parasitic efficacy of both drugs in the acute and chronic phases of Chagas disease. Throughout this review, it is emphasized that the interplay between parasite and host genetics should have an important role in the definition of Chagas disease pathogenesis, anti-T. cruzi immune response and chemotherapy outcome and should be considered in future investigations.

Comprehensive skin microbiome analysis reveals the uniqueness of human skin and evidence for phylosymbiosis within the class Mammalia

Skin forms a critical protective barrier between a mammal and its external environment. Baseline data on the mammalian skin microbiome elucidates which microorganisms are found on healthy skin and provides insight into mammalian evolutionary history. To our knowledge, this study represents the largest existing mammalian skin microbiome survey. Our findings demonstrate that human skin is distinct, not only from other Primates, but from all 10 mammalian orders sampled. Identifying significant similarities between branching of mammalian phylogenetic trees and relatedness trees for their corresponding microbial communities raises the possibility that mammals have experienced coevolution between skin microbiota and their corresponding host species.

Skin is the largest organ of the body and represents the primary physical barrier between mammals and their external environment, yet the factors that govern skin microbial community composition among mammals are poorly understood. The objective of this research was to generate a skin microbiota baseline for members of the class Mammalia, testing the effects of host species, geographic location, body region, and biological sex. Skin from the back, torso, and inner thighs of 177 nonhuman mammals was sampled, representing individuals from 38 species and 10 mammalian orders. Animals were sampled from farms, zoos, households, and the wild. The DNA extracts from all skin swabs were amplified by PCR and sequenced, targeting the V3-V4 regions of bacterial and archaeal 16S rRNA genes. Previously published skin microbiome data from 20 human participants, sampled and sequenced using an identical protocol to the nonhuman mammals, were included to make this a comprehensive analysis. Human skin microbial communities were distinct and significantly less diverse than all other sampled mammalian orders. The factor most strongly associated with microbial community data for all samples was whether the host was a human. Within nonhuman samples, host taxonomic order was the most significant factor influencing skin microbiota, followed by the geographic location of the habitat. By comparing the congruence between host phylogeny and microbial community dendrograms, we observed that Artiodactyla (even-toed ungulates) and Perissodactyla (odd-toed ungulates) had significant congruence, providing evidence of phylosymbiosis between skin microbial communities and their hosts.

Why Cross-Cultural Psychology Is Incomplete Without Comparative and Developmental Perspectives

We argue that comparing adult behavior and cognition across cultures is insufficient to capture the multifaceted complexity of cultural variation. We champion a multidisciplinary perspective that draws on biological and psychological theory and methods. We provide examples for ways in which cross-cultural, developmental, and comparative studies might be combined to unravel the interplay between universal species-typical behaviors and behavioral variation across groups and, at the same time, to explain uniquely human cultural diversity by identifying the unique and universal patterns of human behavior and cognition in early childhood that create, structure, and maintain variation across groups. Such a perspective adds depth to explanations of cultural variation and universality and firmly roots accounts of human culture in a broader, biological framework. We believe that, therefore, the field of cross-cultural psychology may benefit from combining efforts with comparative and developmental psychologists.

A 200K SNP chip reveals a novel Pacific salmon louse genotype linked to differential efficacy of emamectin benzoate

Antiparasitic drugs such as emamectin benzoate (EMB) are relied upon to reduce the parasite load, particularly of the sea louse Lepeophtheirus salmonis, on farmed salmon. The decline in EMB treatment efficacy for this purpose is an important issue for salmon producers around the world, and particularly for those in the Atlantic Ocean where widespread EMB tolerance in sea lice is recognized as a significant problem. Salmon farms in the Northeast Pacific Ocean have not historically experienced the same issues with treatment efficacy, possibly due to the relatively large population of endemic salmonid hosts that serve to both redistribute surviving lice and dilute populations potentially under selection by introducing naïve lice to farms. Frequent migration of lice among farmed and wild hosts should limit the effect of farm-specific selection pressures on changes to the overall allele frequencies of sea lice in the Pacific Ocean. A previous study using microsatellites examined L. salmonis oncorhynchi from 10 Pacific locations from wild and farmed hosts and found no population structure. Recently however, a farm population of sea lice was detected where EMB bioassay exposure tolerance was abnormally elevated. In response, we have developed a Pacific louse draft genome that complements the previously-released Atlantic louse sequence. These genomes were combined with whole-genome re-sequencing data to design a highly sensitive 201,279 marker SNP array applicable for both subspecies (90,827 validated Pacific loci; 153,569 validated Atlantic loci). Notably, kmer spectrum analysis of the re-sequenced samples indicated that Pacific lice exhibit a large within-individual heterozygosity rate (average of 1 in every 72 bases) that is markedly higher than that of Atlantic individuals (1 in every 173 bases). The SNP chip was used to produce a high-density map for Atlantic sea louse linkage group 5 that was previously shown to be associated with EMB tolerance in Atlantic lice. Additionally, 478 Pacific louse samples from farmed and wild hosts obtained between 2005 and 2014 were also genotyped on the array. Clustering analysis allowed us to detect the apparent emergence of an otherwise rare genotype at a high frequency among the lice collected from two farms in 2013 that had reported elevated EMB tolerance. This genotype was not observed in louse samples collected from the same farm in 2010, nor in any lice sampled from other locations prior to 2013. However, this genotype was detected at low frequencies in louse samples from farms in two locations reporting elevated EMB tolerance in 2014. These results suggest that a rare genotype present in Pacific lice may be locally expanded in farms after EMB treatment. Supporting this hypothesis, 437 SNPs associated with this genotype were found to be in a region of linkage group 5 that overlaps the region associated with EMB resistance in Atlantic lice. Finally, five of the top diagnostic SNPs within this region were used to screen lice that had been subjected to an EMB survival assay, revealing a significant association between these SNPs and EMB treatment outcome. To our knowledge this work is the first report to identify a genetic link to altered EMB efficacy in L. salmonis in the Pacific Ocean.

Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics

During the last decade, the vast omics field has revolutionized biological research, especially the genomics, transcriptomics and proteomics branches, as technological tools become available to the field researcher and allow difficult question-driven studies to be addressed. Parasitology has greatly benefited from next generation sequencing (NGS) projects, which have resulted in a broadened comprehension of basic parasite molecular biology, ecology and epidemiology. Malariology is one example where application of this technology has greatly contributed to a better understanding of Plasmodium spp. biology and host-parasite interactions. Among the several parasite species that cause human malaria, the neglected Plasmodium vivax presents great research challenges, as in vitro culturing is not yet feasible and functional assays are heavily limited. Therefore, there are gaps in our P. vivax biology knowledge that affect decisions for control policies aiming to eradicate vivax malaria in the near future. In this review, we provide a snapshot of key discoveries already achieved in P. vivax sequencing projects, focusing on developments, hurdles, and limitations currently faced by the research community, as well as perspectives on future vivax malaria research.

DDX5/p68 associated lncRNA LOC284454 is differentially expressed in human cancers and modulates gene expression

Long non-coding RNAs (lncRNAs) are emerging as important players in regulation of gene expression in higher eukaryotes. DDX5/p68 RNA helicase protein which is involved in splicing of precursor mRNAs also interacts with lncRNAs like, SRA and mrhl, to modulate gene expression. We performed RIP-seq analysis in HEK293T cells to identify the complete repertoire of DDX5/p68 interacting transcripts including 73 single exonic (SE) lncRNAs. The LOC284454 lncRNA is the second top hit of the list of SE lncRNAs which we have characterized in detail for its molecular features and cellular functions. The RNA is located in the same primary transcript harboring miR-23a∼27a∼24-2 cluster. LOC284454 is a stable, nuclear restricted and chromatin associated lncRNA. The sequence is conserved only in primates among 26 different species and is expressed in multiple human tissues. Expression of LOC284454 is significantly reduced in breast, prostate, uterus and kidney cancer and also in breast cancer cell lines (MCF7 and T47D). Global gene expression studies upon loss and gain of function of LOC284454 revealed perturbation of genes related to cancer-related pathways. Focal adhesion and cell migration pathway genes are downregulated under overexpression condition, and these genes are significantly upregulated in breast cancer cell lines as well as breast cancer tissue samples suggesting a functional role of LOC284454 lncRNA in breast cancer pathobiology.