Comparative primate genomics: emerging patterns of genome content and dynamics

Analysis of the aging-related biomarker in a nonhuman primate model using multilayer omics

BACKGROUND

Aging is a prominent risk factor for diverse diseases; therefore, an in-depth understanding of its physiological mechanisms is required. Nonhuman primates, which share the closest genetic relationship with humans, serve as an ideal model for exploring the complex aging process. However, the potential of the nonhuman primate animal model in the screening of human aging markers is still not fully exploited. Multiomics analysis of nonhuman primate peripheral blood offers a promising approach to evaluate new therapies and biomarkers. This study explores aging-related biomarker through multilayer omics, including transcriptomics (mRNA, lncRNA, and circRNA) and proteomics (serum and serum-derived exosomes) in rhesus monkeys (Macaca mulatta).

RESULTS

Our findings reveal that, unlike mRNAs and circRNAs, highly expressed lncRNAs are abundant during the key aging period and are associated with cancer pathways. Comparative analysis highlighted exosomal proteins contain more types of proteins than serum proteins, indicating that serum-derived exosomes primarily regulate aging through metabolic pathways. Finally, eight candidate aging biomarkers were identified, which may serve as blood-based indicators for detecting age-related brain changes.

CONCLUSIONS

Our results provide a comprehensive understanding of nonhuman primate blood transcriptomes and proteomes, offering novel insights into the aging mechanisms for preventing or treating age-related diseases.

Speeding up Glioblastoma Cancer Research: Highlighting the Zebrafish Xenograft Model

Glioblastoma multiforme (GBM) is a very aggressive and lethal primary brain cancer in adults. The multifaceted nature of GBM pathogenesis, rising from complex interactions between cells and the tumor microenvironment (TME), has posed great treatment challenges. Despite significant scientific efforts, the prognosis for GBM remains very poor, even after intensive treatment with surgery, radiation, and chemotherapy. Efficient GBM management still requires the invention of innovative treatment strategies. There is a strong necessity to complete cancer in vitro studies and in vivo studies to properly evaluate the mechanisms of tumor progression within the complex TME. In recent years, the animal models used to study GBM tumors have evolved, achieving highly invasive GBM models able to provide key information on the molecular mechanisms of GBM onset. At present, the most commonly used animal models in GBM research are represented by mammalian models, such as mouse and canine ones. However, the latter present several limitations, such as high cost and time-consuming management, making them inappropriate for large-scale anticancer drug evaluation. In recent years, the zebrafish (Danio rerio) model has emerged as a valuable tool for studying GBM. It has shown great promise in preclinical studies due to numerous advantages, such as its small size, its ability to generate a large cohort of genetically identical offspring, and its rapid development, permitting more time- and cost-effective management and high-throughput drug screening when compared to mammalian models. Moreover, due to its transparent nature in early developmental stages and genetic and anatomical similarities with humans, it allows for translatable brain cancer research and related genetic screening and drug discovery. For this reason, the aim of the present review is to highlight the potential of relevant transgenic and xenograft zebrafish models and to compare them to the traditionally used animal models in GBM research.

Comparative genomics of macaques and integrated insights into genetic variation and population history

The crab-eating macaques ( Macaca fascicularis ) and rhesus macaques ( M. mulatta ) are widely studied nonhuman primates in biomedical and evolutionary research. Despite their significance, the current understanding of the complex genomic structure in macaques and the differences between species requires substantial improvement. Here, we present a complete genome assembly of a crab-eating macaque and 20 haplotype-resolved macaque assemblies to investigate the complex regions and major genomic differences between species. Segmental duplication in macaques is ∼42% lower, while centromeres are ∼3.7 times longer than those in humans. The characterization of ∼2 Mbp fixed genetic variants and ∼240 Mbp complex loci highlights potential associations with metabolic differences between the two macaque species (e.g., CYP2C76 and EHBP1L1 ). Additionally, hundreds of alternative splicing differences show post-transcriptional regulation divergence between these two species (e.g., PNPO ). We also characterize 91 large-scale genomic differences between macaques and humans at a single-base-pair resolution and highlight their impact on gene regulation in primate evolution (e.g., FOLH1 and PIEZO2 ). Finally, population genetics recapitulates macaque speciation and selective sweeps, highlighting potential genetic basis of reproduction and tail phenotype differences (e.g., STAB1 , SEMA3F , and HOXD13 ). In summary, the integrated analysis of genetic variation and population genetics in macaques greatly enhances our comprehension of lineage-specific phenotypes, adaptation, and primate evolution, thereby improving their biomedical applications in human diseases.

Adaptive functions of structural variants in human brain development

Quantifying the structural variants (SVs) in nonhuman primates could provide a niche to clarify the genetic backgrounds underlying human-specific traits, but such resource is largely lacking. Here, we report an accurate SV map in a population of 562 rhesus macaques, verified by in-house benchmarks of eight macaque genomes with long-read sequencing and another one with genome assembly. This map indicates stronger selective constrains on inversions at regulatory regions, suggesting a strategy for prioritizing them with the most important functions. Accordingly, we identified 75 human-specific inversions and prioritized them. The top-ranked inversions have substantially shaped the human transcriptome, through their dual effects of reconfiguring the ancestral genomic architecture and introducing regional mutation hotspots at the inverted regions. As a proof of concept, we linked APCDD1, located on one of these inversions and down-regulated specifically in humans, to neuronal maturation and cognitive ability. We thus highlight inversions in shaping the human uniqueness in brain development.

Spatiotemporal transcriptomic changes of human ovarian aging and the regulatory role of FOXP1

Limited understanding exists regarding how aging impacts the cellular and molecular aspects of the human ovary. This study combines single-cell RNA sequencing and spatial transcriptomics to systematically characterize human ovarian aging. Spatiotemporal molecular signatures of the eight types of ovarian cells during aging are observed. An analysis of age-associated changes in gene expression reveals that DNA damage response may be a key biological pathway in oocyte aging. Three granulosa cells subtypes and five theca and stromal cells subtypes, as well as their spatiotemporal transcriptomics changes during aging, are identified. FOXP1 emerges as a regulator of ovarian aging, declining with age and inhibiting CDKN1A transcription. Silencing FOXP1 results in premature ovarian insufficiency in mice. These findings offer a comprehensive understanding of spatiotemporal variability in human ovarian aging, aiding the prioritization of potential diagnostic biomarkers and therapeutic strategies.

Epigenetic MLH1 silencing concurs with mismatch repair deficiency in sporadic, naturally occurring colorectal cancer in rhesus macaques

BACKGROUND

Naturally occurring colorectal cancers (CRC) in rhesus macaques share many features with their human counterparts and are useful models for cancer immunotherapy; but mechanistic data are lacking regarding the comparative molecular pathogenesis of these cancers.

METHODS

We conducted state-of-the-art imaging including CT and PET, clinical assessments, and pathological review of 24 rhesus macaques with naturally occurring CRC. Additionally, we molecularly characterized these tumors utilizing immunohistochemistry (IHC), microsatellite instability assays, DNAseq, transcriptomics, and developed a DNA methylation-specific qPCR assay for MLH1, CACNA1G, CDKN2A, CRABP1, and NEUROG1, human markers for CpG island methylator phenotype (CIMP). We furthermore employed Monte-Carlo simulations to in-silico model alterations in DNA topology in transcription-factor binding site-rich promoter regions upon experimentally demonstrated DNA methylation.

RESULTS

Similar cancer histology, progression patterns, and co-morbidities could be observed in rhesus as reported for human CRC patients. IHC identified loss of MLH1 and PMS2 in all cases, with functional microsatellite instability. DNA sequencing revealed the close genetic relatedness to human CRCs, including a similar mutational signature, chromosomal instability, and functionally-relevant mutations affecting KRAS (G12D), TP53 (R175H, R273*), APC, AMER1, ALK, and ARID1A. Interestingly, MLH1 mutations were rarely identified on a somatic or germline level. Transcriptomics not only corroborated the similarities of rhesus and human CRCs, but also demonstrated the significant downregulation of MLH1 but not MSH2, MSH6, or PMS2 in rhesus CRCs. Methylation-specific qPCR suggested CIMP-positivity in 9/16 rhesus CRCs, but all 16/16 exhibited significant MLH1 promoter hypermethylation. DNA hypermethylation was modelled to affect DNA topology, particularly propeller twist and roll profiles. Modelling the DNA topology of a transcription factor binding motif (TFAP2A) in the MLH1 promoter that overlapped with a methylation-specific probe, we observed significant differences in DNA topology upon experimentally shown DNA methylation. This suggests a role of transcription factor binding interference in epigenetic silencing of MLH1 in rhesus CRCs.

CONCLUSIONS

These data indicate that epigenetic silencing suppresses MLH1 transcription, induces the loss of MLH1 protein, abrogates mismatch repair, and drives genomic instability in naturally occurring CRC in rhesus macaques. We consider this spontaneous, uninduced CRC in immunocompetent, treatment-naïve rhesus macaques to be a uniquely informative model for human CRC.

Structurally divergent and recurrently mutated regions of primate genomes

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

First whole-genome sequence and assembly of the Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps), a critically endangered species, using Oxford Nanopore Technologies

The Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps) is currently considered one of the most endangered primates in the world and is classified as critically endangered [International union for conservation of nature (IUCN)]. It faces multiple threats, the most significant one being habitat loss due to deforestation in western Ecuador. Genomic tools are keys for the management of endangered species, but this requires a reference genome, which until now was unavailable for A. f. fusciceps. The present study reports the first whole-genome sequence and assembly of A. f. fusciceps generated using Oxford Nanopore long reads. DNA was extracted from a subadult male, and libraries were prepared for sequencing following the Ligation Sequencing Kit SQK-LSK112 workflow. Sequencing was performed using a MinION Mk1C sequencer. The sequencing reads were processed to generate a genome assembly. Two different assemblers were used to obtain draft genomes using raw reads, of which the Flye assembly was found to be superior. The final assembly has a total length of 2.63 Gb and contains 3,861 contigs, with an N50 of 7,560,531 bp. The assembly was analyzed for annotation completeness based on primate ortholog prediction using a high-resolution database, and was found to be 84.3% complete, with a low number of duplicated genes indicating a precise assembly. The annotation of the assembly predicted 31,417 protein-coding genes, comparable with other mammal assemblies. A reference genome for this critically endangered species will allow researchers to gain insight into the genetics of its populations and thus aid conservation and management efforts of this vulnerable species.

On the genetic basis of tail-loss evolution in humans and apes

The loss of the tail is among the most notable anatomical changes to have occurred along the evolutionary lineage leading to humans and to the 'anthropomorphous apes'1-3, with a proposed role in contributing to human bipedalism4-6. Yet, the genetic mechanism that facilitated tail-loss evolution in hominoids remains unknown. Here we present evidence that an individual insertion of an Alu element in the genome of the hominoid ancestor may have contributed to tail-loss evolution. We demonstrate that this Alu element-inserted into an intron of the TBXT gene7-9-pairs with a neighbouring ancestral Alu element encoded in the reverse genomic orientation and leads to a hominoid-specific alternative splicing event. To study the effect of this splicing event, we generated multiple mouse models that express both full-length and exon-skipped isoforms of Tbxt, mimicking the expression pattern of its hominoid orthologue TBXT. Mice expressing both Tbxt isoforms exhibit a complete absence of the tail or a shortened tail depending on the relative abundance of Tbxt isoforms expressed at the embryonic tail bud. These results support the notion that the exon-skipped transcript is sufficient to induce a tail-loss phenotype. Moreover, mice expressing the exon-skipped Tbxt isoform develop neural tube defects, a condition that affects approximately 1 in 1,000 neonates in humans10. Thus, tail-loss evolution may have been associated with an adaptive cost of the potential for neural tube defects, which continue to affect human health today.

The Impact and Effects of Host Immunogenetics on Infectious Disease Studies Using Non-Human Primates in Biomedical Research

Understanding infectious disease pathogenesis and evaluating novel candidate treatment interventions for human use frequently requires prior or parallel analysis in animal model systems. While rodent species are frequently applied in such studies, there are situations where non-human primate (NHP) species are advantageous or required. These include studies of animals that are anatomically more akin to humans, where there is a need to interrogate the complexity of more advanced biological systems or simply reflect susceptibility to a specific infectious agent. The contribution of different arms of the immune response may be addressed in a variety of NHP species or subspecies in specific physiological compartments. Such studies provide insights into immune repertoires not always possible from human studies. However, genetic variation in outbred NHP models may confound, or significantly impact the outcome of a particular study. Thus, host factors need to be considered when undertaking such studies. Considerable knowledge of the impact of host immunogenetics on infection dynamics was elucidated from HIV/SIV research. NHP models are now important for studies of emerging infections. They have contributed to delineating the pathogenesis of SARS-CoV-2/COVID-19, which identified differences in outcomes attributable to the selected NHP host. Moreover, their use was crucial in evaluating the immunogenicity and efficacy of vaccines against COVID-19 and establishing putative correlates of vaccine protection. More broadly, neglected or highly pathogenic emerging or re-emergent viruses may be studied in selected NHPs. These studies characterise protective immune responses following infection or the administration of candidate immunogens which may be central to the accelerated licensing of new vaccines. Here, we review selected aspects of host immunogenetics, specifically MHC background and TRIM5 polymorphism as exemplars of adaptive and innate immunity, in commonly used Old and New World host species. Understanding this variation within and between NHP species will ensure that this valuable laboratory source is used most effectively to combat established and emerging virus infections and improve human health worldwide.

Resilience, aging, and response to radiation exposure (RARRE) in nonhuman primates: a resource review

The Wake Forest nonhuman primate (NHP) Radiation Late Effects Cohort (RLEC) is a unique and irreplaceable population of aging NHP radiation survivors which serves the nation's need to understand the late effects of radiation exposure. Over the past 16 years, Wake Forest has evaluated > 250 previously irradiated rhesus macaques (Macaca mulatta) that were exposed to single total body irradiation (IR) doses of 1.14-8.5 Gy or to partial body exposures of up to 10 Gy (5% bone marrow sparing) or 10.75 Gy (whole thorax). Though primarily used to examine IR effects on disease-specific processes or to develop radiation countermeasures, this resource provides insights on resilience across physiologic systems and its relationship with biological aging. Exposure to IR has well documented deleterious effects on health, but the late effects of IR are highly variable. Some animals exhibit multimorbidity and accumulated health deficits, whereas others remain relatively resilient years after exposure to total body IR. This provides an opportunity to evaluate biological aging at the nexus of resilient/vulnerable responses to a stressor. Consideration of inter-individual differences in response to this stressor can inform individualized strategies to manage late effects of radiation exposure, and provide insight into mechanisms underlying systemic resilience and aging. The utility of this cohort for age-related research questions was summarized at the 2022 Trans-NIH Geroscience Interest Group's Workshop on Animal Models for Geroscience. We present a brief review of radiation injury and its relationship to aging and resilience in NHPs with a focus on the RLEC.

CEGA: a method for inferring natural selection by comparative population genomic analysis across species

We developed maximum likelihood method for detecting positive selection or balancing selection using multilocus or genomic polymorphism and divergence data from two species. The method is especially useful for investigating natural selection in noncoding regions. Simulations demonstrate that the method outperforms existing methods in detecting both positive and balancing selection. We apply the method to population genomic data from human and chimpanzee. The list of genes identified under selection in the noncoding regions is prominently enriched in pathways related to the brain and nervous system. Therefore, our method will serve as a useful tool for comparative population genomic analysis.

Human-specific genetics: new tools to explore the molecular and cellular basis of human evolution

Our ancestors acquired morphological, cognitive and metabolic modifications that enabled humans to colonize diverse habitats, develop extraordinary technologies and reshape the biosphere. Understanding the genetic, developmental and molecular bases for these changes will provide insights into how we became human. Connecting human-specific genetic changes to species differences has been challenging owing to an abundance of low-effect size genetic changes, limited descriptions of phenotypic differences across development at the level of cell types and lack of experimental models. Emerging approaches for single-cell sequencing, genetic manipulation and stem cell culture now support descriptive and functional studies in defined cell types with a human or ape genetic background. In this Review, we describe how the sequencing of genomes from modern and archaic hominins, great apes and other primates is revealing human-specific genetic changes and how new molecular and cellular approaches - including cell atlases and organoids - are enabling exploration of the candidate causal factors that underlie human-specific traits.

Genome sequencing and application of Taiwanese macaque Macaca cyclopis

Formosan macaque (Macaca cyclopis) is the only non-human primate in Taiwan Island. We performed de novo hybrid assembly for M. cyclopis using Illumina paired-end short reads, mate-pair reads and Nanopore long reads and obtained 5065 contigs with a N50 of 2.66 megabases. M. cyclopis contigs >  = 10 kb were assigned to chromosomes using Indian rhesus macaque (Macaca mulatta mulatta) genome assembly Mmul_10 as reference, resulting in a draft of M. cyclopis genome of 2,846,042,475 bases, distributed in 21 chromosomes. The draft genome contains 23,462 transcriptional origins (genes), capable of expressing 716,231 exons in 59,484 transcripts. Genome-based phylogenetic study using the assembled M. cyclopis genome together with genomes of four other macaque species, human, orangutan and chimpanzee showed similar result as previously reported. However, the M. cyclopis species was found to diverge from Chinese M. mulatta lasiota about 1.8 million years ago. Fossil gene analysis detected the presence of gap and pol endogenous viral elements of simian retrovirus in all macaques tested, including M. fascicularis, M. m. mulatta and M. cyclopis. However, M. cyclopis showed ~ 2 times less in number and more uniform in chromosomal locations. The constrain in foreign genome disturbance, presumably due to geographical isolation, should be able to simplify genomics-related investigations, making M. cyclopis an ideal primate species for medical research.

Phylogenomic analyses provide insights into primate evolution

Comparative analysis of primate genomes within a phylogenetic context is essential for understanding the evolution of human genetic architecture and primate diversity. We present such a study of 50 primate species spanning 38 genera and 14 families, including 27 genomes first reported here, with many from previously less well represented groups, the New World monkeys and the Strepsirrhini. Our analyses reveal heterogeneous rates of genomic rearrangement and gene evolution across primate lineages. Thousands of genes under positive selection in different lineages play roles in the nervous, skeletal, and digestive systems and may have contributed to primate innovations and adaptations. Our study reveals that many key genomic innovations occurred in the Simiiformes ancestral node and may have had an impact on the adaptive radiation of the Simiiformes and human evolution.

A revised perspective on the evolution of the lateral frontal cortex in primates

Detailed neuroscientific data from macaque monkeys have been essential in advancing understanding of human frontal cortex function, particularly for regions of frontal cortex without homologs in other model species. However, precise transfer of this knowledge for direct use in human applications requires an understanding of monkey to hominid homologies, particularly whether and how sulci and cytoarchitectonic regions in the frontal cortex of macaques relate to those in hominids. We combine sulcal pattern analysis with resting-state functional magnetic resonance imaging and cytoarchitectonic analysis to show that old-world monkey brains have the same principles of organization as hominid brains, with the notable exception of sulci in the frontopolar cortex. This essential comparative framework provides insights into primate brain evolution and a key tool to drive translation from invasive research in monkeys to human applications.

Structurally divergent and recurrently mutated regions of primate genomes

To better understand the pattern of primate genome structural variation, we sequenced and assembled using multiple long-read sequencing technologies the genomes of eight nonhuman primate species, including New World monkeys (owl monkey and marmoset), Old World monkey (macaque), Asian apes (orangutan and gibbon), and African ape lineages (gorilla, bonobo, and chimpanzee). Compared to the human genome, we identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. Across 50 million years of primate evolution, we estimate that 819.47 Mbp or ~27% of the genome has been affected by SVs based on analysis of these primate lineages. We identify 1,607 structurally divergent regions (SDRs) wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (CARDs, ABCD7, OLAH) and new lineage-specific genes are generated (e.g., CKAP2, NEK5) and have become targets of rapid chromosomal diversification and positive selection (e.g., RGPDs). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species for the first time.

Primate protein-ligand interfaces exhibit significant conservation and unveil human-specific evolutionary drivers

Despite the vast phenotypic differences observed across primates, their protein products are largely similar to each other at the sequence level. We hypothesized that, since proteins accomplish all their functions via interactions with other molecules, alterations in the sites that participate in these interactions may be of critical importance. To uncover the extent to which these sites evolve across primates, we built a structurally-derived dataset of ~4,200 one-to-one orthologous sequence groups across 18 primate species, consisting of ~68,000 ligand-binding sites that interact with DNA, RNA, small molecules, ions, or peptides. Using this dataset, we identify functionally important patterns of conservation and variation within the amino acid residues that facilitate protein-ligand interactions across the primate phylogeny. We uncover that interaction sites are significantly more conserved than other sites, and that sites binding DNA and RNA further exhibit the lowest levels of variation. We also show that the subset of ligand-binding sites that do vary are enriched in components of gene regulatory pathways and uncover several instances of human-specific ligand-binding site changes within transcription factors. Altogether, our results suggest that ligand-binding sites have experienced selective pressure in primates and propose that variation in these sites may have an outsized effect on phenotypic variation in primates through pleiotropic effects on gene regulation.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

Nonhuman primate genetic models for the study of rare diseases

Pre-clinical research and development relies heavily upon translationally valid models of disease. A major difficulty in understanding the biology of, and developing treatments for, rare disease is the lack of animal models. It is important that these models not only recapitulate the presentation of the disease in humans, but also that they share functionally equivalent underlying genetic causes. Nonhuman primates share physiological, anatomical, and behavioral similarities with humans resulting from close evolutionary relationships and high genetic homology. As the post-genomic era develops and next generation sequencing allows for the resequencing and screening of large populations of research animals, naturally occurring genetic variation in nonhuman primates with clinically relevant phenotypes is regularly emerging. Here we review nonhuman primate models of multiple rare genetic diseases with a focus on the similarities and differences in manifestation and etiologies across species. We discuss how these models are being developed and how they can offer new tools and opportunities for researchers interested in exploring novel therapeutics for these and other genetic diseases. Modeling human genetic diseases in translationally relevant nonhuman primates presents new prospects for development of therapeutics and a better understanding of rare diseases. The post-genomic era offers the opportunity for the discovery and further development of more models like those discussed here.

Current advances in primate genomics: novel approaches for understanding evolution and disease

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

Omics and Remote Homology Integration to Decipher Protein Functionality

In the recent years, several "omics" technologies based on specific biomolecules (from DNA, RNA, proteins, or metabolites) have won growing importance in the scientific field. Despite each omics possess their own laboratorial protocols, they share a background of bioinformatic tools for data integration and analysis. A recent subset of bioinformatic tools, based on available templates or remote homology protocols, allow computational fast and high-accuracy prediction of protein structures. The quickly predict of actually unsolved protein structures, together with late omics findings allow a boost of scientific advances in multiple fields such as cancer, longevity, immunity, mitochondrial function, toxicology, drug design, biosensors, and recombinant protein engineering. In this chapter, we assessed methodological approaches for the integration of omics and remote homology inferences to decipher protein functionality, opening the door to the next era of biological knowledge.

Real or fake? Measuring the impact of protein annotation errors on estimates of domain gain and loss events

Protein annotation errors can have significant consequences in a wide range of fields, ranging from protein structure and function prediction to biomedical research, drug discovery, and biotechnology. By comparing the domains of different proteins, scientists can identify common domains, classify proteins based on their domain architecture, and highlight proteins that have evolved differently in one or more species or clades. However, genome-wide identification of different protein domain architectures involves a complex error-prone pipeline that includes genome sequencing, prediction of gene exon/intron structures, and inference of protein sequences and domain annotations. Here we developed an automated fact-checking approach to distinguish true domain loss/gain events from false events caused by errors that occur during the annotation process. Using genome-wide ortholog sets and taking advantage of the high-quality human and Saccharomyces cerevisiae genome annotations, we analyzed the domain gain and loss events in the predicted proteomes of 9 non-human primates (NHP) and 20 non-S. cerevisiae fungi (NSF) as annotated in the Uniprot and Interpro databases. Our approach allowed us to quantify the impact of errors on estimates of protein domain gains and losses, and we show that domain losses are over-estimated ten-fold and three-fold in the NHP and NSF proteins respectively. This is in line with previous studies of gene-level losses, where issues with genome sequencing or gene annotation led to genes being falsely inferred as absent. In addition, we show that insistent protein domain annotations are a major factor contributing to the false events. For the first time, to our knowledge, we show that domain gains are also over-estimated by three-fold and two-fold respectively in NHP and NSF proteins. Based on our more accurate estimates, we infer that true domain losses and gains in NHP with respect to humans are observed at similar rates, while domain gains in the more divergent NSF are observed twice as frequently as domain losses with respect to S. cerevisiae. This study highlights the need to critically examine the scientific validity of protein annotations, and represents a significant step toward scalable computational fact-checking methods that may 1 day mitigate the propagation of wrong information in protein databases.

“Natural Laboratory Complex” for novel primate neuroscience

We propose novel strategies for primate experimentation that are ethically valuable and pragmatically useful for cognitive neuroscience and neuropsychiatric research. Specifically, we propose Natural Laboratory Complex or Natural Labs, which are a combination of indoor-outdoor structures for studying free moving and socially housed primates in natural or naturalistic environment. We contend that Natural Labs are pivotal to improve primate welfare, and at the same time to implement longitudinal and socio-ecological studies of primate brain and behavior. Currently emerging advanced technologies and social systems (including recent COVID-19 induced “remote” infrastructures) can speed-up cognitive neuroscience approaches in freely behaving animals. Experimental approaches in natural(istic) settings are not in competition with conventional approaches of laboratory investigations, and could establish several benefits at the ethical, experimental, and economic levels.

Effects of Anesthesia on Cerebral Blood Flow and Functional Connectivity of Nonhuman Primates

Nonhuman primates (NHPs) are the closest living relatives of humans and play a critical and unique role in neuroscience research and pharmaceutical development. General anesthesia is usually required in neuroimaging studies of NHPs to keep the animal from stress and motion. However, the adverse effects of anesthesia on cerebral physiology and neural activity are pronounced and can compromise the data collection and interpretation. Functional connectivity is frequently examined using resting-state functional MRI (rsfMRI) to assess the functional abnormality in the animal brain under anesthesia. The fMRI signal can be dramatically suppressed by most anesthetics in a dose-dependent manner. In addition, rsfMRI studies may be further compromised by inter-subject variations when the sample size is small (as seen in most neuroscience studies of NHPs). Therefore, proper use of anesthesia is strongly demanded to ensure steady and consistent physiology maintained during rsfMRI data collection of each subject. The aim of this review is to summarize typical anesthesia used in rsfMRI scans of NHPs and the effects of anesthetics on cerebral physiology and functional connectivity. Moreover, the protocols with optimal rsfMRI data acquisition and anesthesia procedures for functional connectivity study of macaque monkeys are introduced.

Total-Body Irradiation Is Associated With Increased Incidence of Mesenchymal Neoplasia in a Radiation Late Effects Cohort of Rhesus Macaques (Macaca mulatta)

PURPOSE

Cancer is a severe delayed effect of acute radiation exposure. Total-body irradiation has been associated with an increased risk of solid cancer and leukemia in Japanese atomic bomb survivors, and secondary malignancies, such as sarcoma, are a serious consequence of cancer radiation therapy. The radiation late effects cohort (RLEC) of rhesus macaques (Macaca mulatta) is a unique resource of more than 200 animals for studying the long-term consequences of total-body irradiation in an animal model that closely resembles humans at the genetic and physiologic levels.

METHODS AND MATERIALS

Using clinical records, clinical imaging, histopathology, and immunohistochemistry, this retrospective study characterized the incidence of neoplasia in the RLEC.

RESULTS

Since 2007, 61 neoplasms in 44 of 239 irradiated animals were documented (18.4% of the irradiated population). Only 1 neoplasm was diagnosed among the 51 nonirradiated controls of the RLEC (2.0%). The most common malignancies in the RLEC were sarcomas (38.3% of diagnoses), which are rare neoplasms in nonirradiated macaques. The most common sarcomas included malignant nerve sheath tumors and malignant glomus tumors. Carcinomas were less common (19.7% of diagnoses), and consisted primarily of renal cell and hepatocellular carcinomas. Neoplasia occurred in most major body systems, with the skin and subcutis being the most common site (40%). RNA analysis showed similarities in transcriptional profiles between RLEC and human malignant nerve sheath tumors.

CONCLUSIONS

This study indicates that total-body irradiation is associated with an increased incidence of neoplasia years following irradiation, at more than double the incidence described in aging, nonirradiated animals, and promotes tumor histotypes that are rarely observed in nonirradiated, aging rhesus macaques.

Society of Toxicologic Pathology Neuropathology Interest Group Article: Neuropathologic Findings in Nonhuman Primates Associated With Administration of Biomolecule-Based Test Articles

The increasing specificity of novel druggable targets coupled with the complexity of emerging therapeutic modalities for treating human diseases has created a growing need for nonhuman primates (NHPs) as models for translational drug discovery and nonclinical safety assessment. In particular, NHPs are critical for investigating potential unexpected/undesired on-target and off-target liabilities associated with administration of candidate biotherapeutics (nucleic acids, proteins, viral gene therapy vectors, etc.) to treat nervous system disorders. Nervous system findings unique to or overrepresented in NHPs administered biomolecule-based ("biologic") test articles include mononuclear cell infiltration in most neural tissues for all biomolecule classes as well as neuronal necrosis with glial cell proliferation in sensory ganglia for certain viral vectors. Such test article-related findings in NHPs often must be differentiated from procedural effects (e.g., local parenchymal or meningeal reactions associated with an injection site or implanted catheter to administer a test article directly into the central nervous system) or spontaneous background findings (e.g., neuronal autophagy in sensory ganglia).

Virtual ChIP-seq: predicting transcription factor binding by learning from the transcriptome

Existing methods for computational prediction of transcription factor (TF) binding sites evaluate genomic regions with similarity to known TF sequence preferences. Most TF binding sites, however, do not resemble known TF sequence motifs, and many TFs are not sequence-specific. We developed Virtual ChIP-seq, which predicts binding of individual TFs in new cell types, integrating learned associations with gene expression and binding, TF binding sites from other cell types, and chromatin accessibility data in the new cell type. This approach outperforms methods that predict TF binding solely based on sequence preference, predicting binding for 36 TFs (MCC>0.3).

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

Rocks and clocks revised: New promises and challenges in dating the primate tree of life

In recent years, multiple technological and methodological advances have increased our ability to estimate phylogenies, leading to more accurate dating of the primate tree of life. Here we provide an overview of the limitations and potentials of some of these advancements and discuss how dated phylogenies provide the crucial temporal scale required to understand primate evolution. First, we review new methods, such as the total-evidence dating approach, that promise a better integration between the fossil record and molecular data. We then explore how the ever-increasing availability of genomic-level data for more primate species can impact our ability to accurately estimate timetrees. Finally, we discuss more recent applications of mutation rates to date divergence times. We highlight example studies that have applied these approaches to estimate divergence dates within primates. Our goal is to provide a critical overview of these new developments and explore the promises and challenges of their application in evolutionary anthropology.

Animal Models of Neurodegenerative Disease: Recent Advances in Fly Highlight Innovative Approaches to Drug Discovery

Neurodegenerative diseases represent a formidable challenge to global health. As advances in other areas of medicine grant healthy living into later decades of life, aging diseases such as Alzheimer's disease (AD) and other neurodegenerative disorders can diminish the quality of these additional years, owed largely to the lack of efficacious treatments and the absence of durable cures. Alzheimer's disease prevalence is predicted to more than double in the next 30 years, affecting nearly 15 million Americans, with AD-associated costs exceeding $1 billion by 2050. Delaying onset of AD and other neurodegenerative diseases is critical to improving the quality of life for patients and reducing the burden of disease on caregivers and healthcare systems. Significant progress has been made to model disease pathogenesis and identify points of therapeutic intervention. While some researchers have contributed to our understanding of the proteins and pathways that drive biological dysfunction in disease using in vitro and in vivo models, others have provided mathematical, biophysical, and computational technologies to identify potential therapeutic compounds using in silico modeling. The most exciting phase of the drug discovery process is now: by applying a target-directed approach that leverages the strengths of multiple techniques and validates lead hits using Drosophila as an animal model of disease, we are on the fast-track to identifying novel therapeutics to restore health to those impacted by neurodegenerative disease.

Comparative molecular genomic analyses of a spontaneous rhesus macaque model of mismatch repair-deficient colorectal cancer

Colorectal cancer (CRC) remains the third most common cancer in the US with 15% of cases displaying Microsatellite Instability (MSI) secondary to Lynch Syndrome (LS) or somatic hypermethylation of the MLH1 promoter. A cohort of rhesus macaques from our institution developed spontaneous mismatch repair deficient (MMRd) CRC with a notable fraction harboring a pathogenic germline mutation in MLH1 (c.1029C Collapse

Key Words Collapse

MESH Headings

• Animals
• Colonic Neoplasms/genetics
• Colorectal Neoplasms/genetics
• Colorectal Neoplasms/pathology
• Colorectal Neoplasms, Hereditary Nonpolyposis/genetics
• Colorectal Neoplasms, Hereditary Nonpolyposis/pathology
• DNA Methylation/genetics
• DNA Mismatch Repair/genetics
• Genomics
• Humans
• Macaca mulatta/genetics
• Mice
• Microsatellite Instability
• MutL Protein Homolog 1/genetics

Collapse

Grants

• P30 CA016672 NCI NIH HHS
• R24 OD011173 NIH HHS
• Feinberg Family Foundation
• MDACC Institutional Research Grant Program
• Cattleman for Cancer Research
• National Institutes of Health
• National Cancer Institute

Collapse

Affiliation(s)

Nejla Ozirmak Lermi Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America School of Health Professions, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Stanton B. Gray Michale E. Keeling Center for Comparative Medicine and Research, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Charles M. Bowen Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Laura Reyes-Uribe Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Beth K. Dray Charles River Laboratories, Ashland, Ohio, United States of America
Nan Deng Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
R. Alan Harris Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Muthuswamy Raveendran Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Fernando Benavides Department of Epigenetics and Molecular Carcinogenesis, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Carolyn L. Hodo Michale E. Keeling Center for Comparative Medicine and Research, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Melissa W. Taggart Department of Pathology, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Karen Colbert Maresso Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Krishna M. Sinha Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Jeffrey Rogers Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Eduardo Vilar Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America Clinical Cancer Genetics Program, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America

Collapse

Flexible auditory training, psychophysics, and enrichment of common marmosets with an automated, touchscreen-based system

Devising new and more efficient protocols to analyze the phenotypes of non-human primates, as well as their complex nervous systems, is rapidly becoming of paramount importance. This is because with genome-editing techniques, recently adopted to non-human primates, new animal models for fundamental and translational research have been established. One aspect in particular, namely cognitive hearing, has been difficult to assess compared to visual cognition. To address this, we devised autonomous, standardized, and unsupervised training and testing of auditory capabilities of common marmosets with a cage-based standalone, wireless system. All marmosets tested voluntarily operated the device on a daily basis and went from naïve to experienced at their own pace and with ease. Through a series of experiments, here we show, that animals autonomously learn to associate sounds with images; to flexibly discriminate sounds, and to detect sounds of varying loudness. The developed platform and training principles combine in-cage training of common marmosets for cognitive and psychoacoustic assessment with an enriched environment that does not rely on dietary restriction or social separation, in compliance with the 3Rs principle.

The authors present a cage-based stand-alone platform for autonomous, standardized, and unsupervised training and testing of visuo-auditory-cued behaviours of common marmosets. The experiments do not require dietary restriction or social separation.

Modeling genetic diseases in nonhuman primates through embryonic and germline modification: Considerations and challenges

Genetic modification of the embryo or germ line of nonhuman primates is envisioned as a method to develop improved models of human disease, yet the promise of such animal models remains unfulfilled. Here, we discuss current methods and their limitations for producing nonhuman primate genetic models that faithfully genocopy and phenocopy human disease. We reflect on how to ethically maximize the translational relevance of such models in the search for new therapeutic strategies to treat human disease.

Importance of Nonhuman Primates as a Model System for Gene Therapy Development in Ophthalmology

Gene therapy is a treatment concept that uses, in most cases, viral vectors to deliver a therapeutic transgene to target cells. Although the idea of gene therapy dates back over 50 years ago, due to the complexity of the treatment concept, it took until the last decade for the responsible agencies like FDA and EMA to recommend the first gene therapy products for clinical use. The development of these therapies relies on molecular engineering of specifically designed vectors and models to test the effectiveness and safety of the treatment. Despite an increasing effort to find effective surrogates, animal models are still irreplaceable in gene therapy development. Rodents are important for exploring pathways and disease mechanisms and identifying potential treatment targets. However, only the primate eye resembles the human eye to a degree where most structures are nearly identical. Some research questions can therefore only be answered using a nonhuman primate (NHP) model. In this review, we want to summarize these key features and highlight the importance of the NHP model for gene therapy development in ophthalmology.

Social pressure drives "conversational rules" in great apes

In the last decade, two hypotheses, one on the evolution of animal vocal communication in general and the other on the origins of human language, have gained ground. The first hypothesis argues that the complexity of communication co-evolved with the complexity of sociality. Species forming larger groups with complex social networks have more elaborate vocal repertoires. The second hypothesis posits that the core of communication is represented not only by what can be expressed by an isolated caller, but also by the way that vocal interactions are structured, language being above all a social act. Primitive forms of conversational rules based on a vocal turn-taking principle are thought to exist in primates. To support and bring together these hypotheses, more comparative studies of socially diverse species at different levels of the primate phylogeny are needed. However, the majority of available studies focus on monkeys, primates that are distant from the human lineage. Great apes represent excellent candidates for such comparative studies because of their phylogenetic proximity to humans and their varied social lives. We propose that studying vocal turn-taking in apes could address several major gaps regarding the social relevance of vocal turn-taking and the evolutionary trajectory of this behaviour among anthropoids. Indeed, how the social structure of a species may influence the vocal interaction patterns observed among group members remains an open question. We gathered data from the literature as well as original unpublished data (where absent in the literature) on four great ape species: chimpanzees Pan troglodytes, bonobos Pan paniscus, western lowland gorillas Gorilla gorilla gorilla and Bornean orang-utans Pongo pygmaeus. We found no clear-cut relationship between classical social complexity metrics (e.g. number of group members, interaction rates) and vocal complexity parameters (e.g. repertoire size, call rates). Nevertheless, the nature of the society (i.e. group composition, diversity and valence of social bonds) and the type of vocal interaction patterns (isolated calling, call overlap, turn-taking-based vocal exchanges) do appear to be related. Isolated calling is the main vocal pattern found in the species with the smallest social networks (orang-utan), while the other species show vocal interactions that are structured according to temporal rules. A high proportion of overlapping vocalisations is found in the most competitive species (chimpanzee), while vocal turn-taking predominates in more tolerant bonobos and gorillas. Also, preferentially interacting individuals and call types used to interact are not randomly distributed. Vocal overlap ('chorusing') and vocal exchange ('conversing') appear as possible social strategies used to advertise/strengthen social bonds. Our analyses highlight that: (i) vocal turn-taking is also observed in non-human great apes, revealing universal rules for conversing that may be deeply rooted in the primate lineage; (ii) vocal interaction patterns match the species' social lifestyle; (iii) although limited to four species here, adopting a targeted comparative approach could help to identify the multiple and subtle factors underlying social and vocal complexity. We believe that vocal interaction patterns form the basis of a promising field of investigation that may ultimately improve our understanding of the socially driven evolution of communication.

A population-level invasion by transposable elements triggers genome expansion in a fungal pathogen

Genome evolution is driven by the activity of transposable elements (TEs). The spread of TEs can have deleterious effects including the destabilization of genome integrity and expansions. However, the precise triggers of genome expansions remain poorly understood because genome size evolution is typically investigated only among deeply divergent lineages. Here, we use a large population genomics dataset of 284 individuals from populations across the globe of Zymoseptoria tritici, a major fungal wheat pathogen. We built a robust map of genome-wide TE insertions and deletions to track a total of 2456 polymorphic loci within the species. We show that purifying selection substantially depressed TE frequencies in most populations, but some rare TEs have recently risen in frequency and likely confer benefits. We found that specific TE families have undergone a substantial genome-wide expansion from the pathogen's center of origin to more recently founded populations. The most dramatic increase in TE insertions occurred between a pair of North American populations collected in the same field at an interval of 25 years. We find that both genome-wide counts of TE insertions and genome size have increased with colonization bottlenecks. Hence, the demographic history likely played a major role in shaping genome evolution within the species. We show that both the activation of specific TEs and relaxed purifying selection underpin this incipient expansion of the genome. Our study establishes a model to recapitulate TE-driven genome evolution over deeper evolutionary timescales.

Schistosoma mansoni and soil transmtted helminths in olive baboons and potential zoonosis

Zoonotic pathogens are among the most important causes of ill health all over the world. The presence of these pathogens in free ranging baboons may have significant implications for humans. In Kenya, the encroachment of wildlife habitats has led to increased interaction between humans and wildlife especially non-human primates. The current study therefore aimed at investigating any possible zoonotic gastrointestinal helminths of olive baboons (Papio anubis) at the human-wildlife interface in two park borders and a ranch in Kenya, namely, Tsavo West National Park, Tana River Primate Reserve and Mutara Ranch, Laikipia, Kenya. One hundred and forty-seven baboons were used in the study. They were trapped in the wild, sampled for stool marked and then released back to the wild. Gastrointestinal (GIT) helminths identified were Strongyloides, Oesophagostomum, Enterobius spp and Trichuris Trichiura from all the three sites while Schistosoma mansoni was only detected from Tsavo baboons and with very low incidence (2.1%). The prevalence of these parasites varied among the sites but significant difference in prevalence was only noted in Strongyloides and Oesophagostomum (p < 0.05) among the three sites. This therefore implies that even with control measures instituted on the human population, baboons will always be a source of zoonotic GIT helminths especially S. mansoni even if the incidence are low. There is need to put in place measures aiming to reduce their interactions with humans and also try to control these infections in the baboons.

Evolutionary and functional lessons from human-specific amino acid substitution matrices

Human genetic variation in coding regions is fundamental to the study of protein structure and function. Most methods for interpreting missense variants consider substitution measures derived from homologous proteins across different species. In this study, we introduce human-specific amino acid (AA) substitution matrices that are based on genetic variations in the modern human population. We analyzed the frequencies of >4.8M single nucleotide variants (SNVs) at codon and AA resolution and compiled human-centric substitution matrices that are fundamentally different from classic cross-species matrices (e.g. BLOSUM, PAM). Our matrices are asymmetric, with some AA replacements showing significant directional preference. Moreover, these AA matrices are only partly predicted by nucleotide substitution rates. We further test the utility of our matrices in exposing functional signals of experimentally-validated protein annotations. A significant reduction in AA transition frequencies was observed across nine post-translational modification (PTM) types and four ion-binding sites. Our results propose a purifying selection signal in the human proteome across a diverse set of functional protein annotations and provide an empirical baseline for interpreting human genetic variation in coding regions.

Variational inference using approximate likelihood under the coalescent with recombination

Coalescent methods are proven and powerful tools for population genetics, phylogenetics, epidemiology, and other fields. A promising avenue for the analysis of large genomic alignments, which are increasingly common, is coalescent hidden Markov model (coalHMM) methods, but these methods have lacked general usability and flexibility. We introduce a novel method for automatically learning a coalHMM and inferring the posterior distributions of evolutionary parameters using black-box variational inference, with the transition rates between local genealogies derived empirically by simulation. This derivation enables our method to work directly with three or four taxa and through a divide-and-conquer approach with more taxa. Using a simulated data set resembling a human–chimp–gorilla scenario, we show that our method has comparable or better accuracy to previous coalHMM methods. Both species divergence times and population sizes were accurately inferred. The method also infers local genealogies, and we report on their accuracy. Furthermore, we discuss a potential direction for scaling the method to larger data sets through a divide-and-conquer approach. This accuracy means our method is useful now, and by deriving transition rates by simulation, it is flexible enough to enable future implementations of various population models.

A comprehensive analysis of gorilla-specific LINE-1 retrotransposons

BACKGROUND

Long interspersed element-1 (LINE-1 or L1) is the most abundant retrotransposons in the primate genome. They have approximately 520,000 copies and make up ~ 17% of the primate genome. Full-length L1s can mobilize to a new genomic location using their enzymatic machinery. Gorilla is the second closest species to humans after the chimpanzee, and human-gorilla split 7-12 million years ago. The gorilla genome provides an opportunity to explore primate origins and evolution.

OBJECTIVE

L1s have contributed to genome diversity and variations during primate evolution. This study aimed to identify gorilla-specific L1s using a more recent version of the gorilla reference genome (Mar. 2016 GSMRT3/gorGor5).

METHODS

We collected gorilla-specific L1 candidates through computational analysis and manual inspection. L1Xplorer was used to identify whether full-length gorilla-specific L1s were intact. In addition, to determine the level of sequence conservation between intact fulllength gorilla-specific L1s, two ORFs of intact L1s were aligned with the L1PA2 consensus sequence.

RESULTS

2002 gorilla-specific L1 candidates were identified through computational analysis. Among them, we manually inspected 1,883 gorilla-specific L1s, among which most of them belong to the L1PA2 subfamily and 12 were intact L1s that could influence genomic variations in the gorilla genome. Interestingly, the 12 intact full-length gorilla-specific L1s have 14 highly conserved nonsynonymous mutations, including 6 mutations and 8 mutations in ORF1 and ORF2, respectively. In comparison to the intact full-length chimpanzee-specific L1s and human-specific hot-L1s, two of these in ORF1 (L256F and E293G) were shown as gorilla-specific nonsynonymous mutations.

CONCLUSION

The gorilla-specific L1s may have had significantly affected the gorilla genome to compose a genome different form that of other primates during primate evolution.

MiREDiBase, a manually curated database of validated and putative editing events in microRNAs

MicroRNAs (miRNAs) are regulatory small non-coding RNAs that function as translational repressors. MiRNAs are involved in most cellular processes, and their expression and function are presided by several factors. Amongst, miRNA editing is an epitranscriptional modification that alters the original nucleotide sequence of selected miRNAs, possibly influencing their biogenesis and target-binding ability. A-to-I and C-to-U RNA editing are recognized as the canonical types, with the A-to-I type being the predominant one. Albeit some bioinformatics resources have been implemented to collect RNA editing data, it still lacks a comprehensive resource explicitly dedicated to miRNA editing. Here, we present MiREDiBase, a manually curated catalog of editing events in miRNAs. The current version includes 3,059 unique validated and putative editing sites from 626 pre-miRNAs in humans and three primates. Editing events in mature human miRNAs are supplied with miRNA-target predictions and enrichment analysis, while minimum free energy structures are inferred for edited pre-miRNAs. MiREDiBase represents a valuable tool for cell biology and biomedical research and will be continuously updated and expanded at https://ncrnaome.osumc.edu/miredibase .

Transcriptome analysis of non human primate-induced pluripotent stem cell-derived cardiomyocytes in 2D monolayer culture vs. 3D engineered heart tissue

AIMS

Stem cell therapy has shown promise for treating myocardial infarction via re-muscularization and paracrine signalling in both small and large animals. Non-human primates (NHPs), such as rhesus macaques (Macaca mulatta), are primarily utilized in preclinical trials due to their similarity to humans, both genetically and physiologically. Currently, induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) are delivered into the infarcted myocardium by either direct cell injection or an engineered tissue patch. Although both approaches have advantages in terms of sample preparation, cell-host interaction, and engraftment, how the iPSC-CMs respond to ischaemic conditions in the infarcted heart under these two different delivery approaches remains unclear. Here, we aim to gain a better understanding of the effects of hypoxia on iPSC-CMs at the transcriptome level.

METHODS AND RESULTS

NHP iPSC-CMs in both monolayer culture (2D) and engineered heart tissue (EHT) (3D) format were exposed to hypoxic conditions to serve as surrogates of direct cell injection and tissue implantation in vivo, respectively. Outcomes were compared at the transcriptome level. We found the 3D EHT model was more sensitive to ischaemic conditions and similar to the native in vivo myocardium in terms of cell-extracellular matrix/cell-cell interactions, energy metabolism, and paracrine signalling.

CONCLUSION

By exposing NHP iPSC-CMs to different culture conditions, transcriptome profiling improves our understanding of the mechanism of ischaemic injury.

Transcriptomic analysis of primate placentas and novel rhesus trophoblast cell lines informs investigations of human placentation

BACKGROUND

Proper placentation, including trophoblast differentiation and function, is essential for the health and well-being of both the mother and baby throughout pregnancy. Placental abnormalities that occur during the early stages of development are thought to contribute to preeclampsia and other placenta-related pregnancy complications. However, relatively little is known about these stages in humans due to obvious ethical and technical limitations. Rhesus macaques are considered an ideal surrogate for studying human placentation, but the unclear translatability of known human placental markers and lack of accessible rhesus trophoblast cell lines can impede the use of this animal model.

RESULTS

Here, we performed a cross-species transcriptomic comparison of human and rhesus placenta and determined that while the majority of human placental marker genes (HPGs) were similarly expressed, 952 differentially expressed genes (DEGs) were identified between the two species. Functional enrichment analysis of the 447 human-upregulated DEGs, including ADAM12, ERVW-1, KISS1, LGALS13, PAPPA2, PGF, and SIGLEC6, revealed over-representation of genes implicated in preeclampsia and other pregnancy disorders. Additionally, to enable in vitro functional studies of early placentation, we generated and thoroughly characterized two highly pure first trimester telomerase (TERT) immortalized rhesus trophoblast cell lines (iRP-D26 and iRP-D28A) that retained crucial features of isolated primary trophoblasts.

CONCLUSIONS

Overall, our findings help elucidate the molecular translatability between human and rhesus placenta and reveal notable expression differences in several HPGs and genes implicated in pregnancy complications that should be considered when using the rhesus animal model to study normal and pathological human placentation.

Self-assembled peptide and protein nanostructures for anti-cancer therapy: Targeted delivery, stimuli-responsive devices and immunotherapy

Self-assembled peptides and proteins possess tremendous potential as targeted drug delivery systems and key applications of these well-defined nanostructures reside in anti-cancer therapy. Peptides and proteins can self-assemble into nanostructures of diverse sizes and shapes in response to changing environmental conditions such as pH, temperature, ionic strength, as well as host and guest molecular interactions; their countless benefits include good biocompatibility and high loading capacity for hydrophobic and hydrophilic drugs. These self-assembled nanomaterials can be adorned with functional moieties to specifically target tumor cells. Stimuli-responsive features can also be incorporated with respect to the tumor microenvironment. This review sheds light on the growing interest in self-assembled peptides and proteins and their burgeoning applications in cancer treatment and immunotherapy.

Ticks, Hair Loss, and Non-Clinging Babies: A Novel Tick-Based Hypothesis for the Evolutionary Divergence of Humans and Chimpanzees

Human straight-legged bipedalism represents one of the earliest events in the evolutionary split between humans (Homo spp.) and chimpanzees (Pan spp.), although its selective basis is a mystery. A carrying-related hypothesis has recently been proposed in which hair loss within the hominin lineage resulted in the inability of babies to cling to their mothers, requiring mothers to walk upright to carry their babies. However, a question remains for this model: what drove the hair loss that resulted in upright walking? Observers since Darwin have suggested that hair loss in humans may represent an evolutionary strategy for defence against ticks. The aim of this review is to propose and evaluate a novel tick-based evolutionary hypothesis wherein forest fragmentation in hominin paleoenvironments created conditions that were favourable for tick proliferation, selecting for hair loss in hominins and grooming behaviour in chimpanzees as divergent anti-tick strategies. It is argued that these divergent anti-tick strategies resulted in different methods for carrying babies, driving the locomotor divergence of humans and chimpanzees.

Mitogenomic phylogeny of Callithrix with special focus on human transferred taxa

Background

Callithrix marmosets are a relatively young primate radiation, whose phylogeny is not yet fully resolved. These primates are naturally para- and allopatric, but three species with highly invasive potential have been introduced into the southeastern Brazilian Atlantic Forest by the pet trade. There, these species hybridize with each other and endangered, native congeners. We aimed here to reconstruct a robust Callithrix phylogeny and divergence time estimates, and identify the biogeographic origins of autochthonous and allochthonous Callithrix mitogenome lineages. We sequenced 49 mitogenomes from four species (C. aurita, C. geoffroyi, C. jacchus, C. penicillata) and anthropogenic hybrids (C. aurita x Callithrix sp., C. penicillata x C. jacchus, Callithrix sp. x Callithrix sp., C. penicillata x C. geoffroyi) via Sanger and whole genome sequencing. We combined these data with previously published Callithrix mitogenomes to analyze five Callithrix species in total.

Results

We report the complete sequence and organization of the C. aurita mitogenome. Phylogenetic analyses showed that C. aurita was the first to diverge within Callithrix 3.54 million years ago (Ma), while C. jacchus and C. penicillata lineages diverged most recently 0.5 Ma as sister clades. MtDNA clades of C. aurita, C. geoffroyi, and C. penicillata show intraspecific geographic structure, but C. penicillata clades appear polyphyletic. Hybrids, which were identified by phenotype, possessed mainly C. penicillata or C. jacchus mtDNA haplotypes. The biogeographic origins of mtDNA haplotypes from hybrid and allochthonous Callithrix were broadly distributed across natural Callithrix ranges. Our phylogenetic results also evidence introgression of C. jacchus mtDNA into C. aurita.

Conclusion

Our robust Callithrix mitogenome phylogeny shows C. aurita lineages as basal and C. jacchus lineages among the most recent within Callithrix. We provide the first evidence that parental mtDNA lineages of anthropogenic hybrid and allochthonous marmosets are broadly distributed inside and outside of the Atlantic Forest. We also show evidence of cryptic hybridization between allochthonous Callithrix and autochthonous C. aurita. Our results encouragingly show that further development of genomic resources will allow to more clearly elucidate Callithrix evolutionary relationships and understand the dynamics of Callithrix anthropogenic introductions into the Brazilian Atlantic Forest.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07533-1.

Population Genomics Reveals Incipient Speciation, Introgression, and Adaptation in the African Mona Monkey (Cercopithecus mona)

Guenons (tribe Cercopithecini) are the most widely distributed nonhuman primate in the tropical forest belt of Africa and show considerable phenotypic, taxonomic, and ecological diversity. However, genomic information for most species within this group is still lacking. Here, we present a high-quality de novo genome (total 2.90 Gb, contig N50 equal to 22.7 Mb) of the mona monkey (Cercopithecus mona), together with genome resequencing data of 13 individuals sampled across Nigeria. Our results showed differentiation between populations from East and West of the Niger River ∼84 ka and potential ancient introgression in the East population from other mona group species. The PTPRK, FRAS1, BNC2, and EDN3 genes related to pigmentation displayed signals of introgression in the East population. Genomic scans suggest that immunity genes such as AKT3 and IL13 (possibly involved in simian immunodeficiency virus defense), and G6PD, a gene involved in malaria resistance, are under positive natural selection. Our study gives insights into differentiation, natural selection, and introgression in guenons.

Evidence for the loss and recovery of SLAMF9 during human evolution: implications on Dollo's law

Signaling lymphocyte activation molecule family member 9 (SLAMF9) is a cell surface protein of the CD2/SLAM family of leukocyte surface receptors. It is conserved throughout mammals and has roles in the initiation of inflammatory responses and regulation of plasmacytoid dendritic cell function. Through comparison of reference sequences encoding SLAMF9 in human, mouse, and primate sequences, we have determined that the SLAMF9 gene underwent successive mutation events, resulting in the loss of the protein and subsequent recovery of a less stable version. The mutations included a single base pair deletion in the second exon and a change in the splice acceptor site of that same exon. These changes would have had the effect of creating and later repairing a frameshift in the coding sequence. These events took place since the divergence of the human lineage from the chimpanzee-human last common ancestor and represent the first known case of the functional loss and recovery of a gene within the human lineage.