Illuminating the function of the orphan transporter, SLC22A10, in humans and other primates.

SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17ß-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.

Curcumin alleviates myocardial inflammation, apoptosis, and oxidative stress induced by acute pulmonary embolism by regulating microRNA-145-5P/insulin receptor substrate 1 axis.

The treatment of patients with acute pulmonary embolism (APE) is extremely challenging due to the complex clinical presentation and prognosis of APE related to the patient's hemodynamic status and insufficient arterial blood flow and right ventricular overload. Protective efficacy against cardiovascular diseases of curcumin, a common natural polyphenolic compound, which has antithrombotic properties and reduces platelet accumulation in the circulation by inhibiting thromboxane synthesis has been demonstrated. However, the direct effect of curcumin on APE has rarely been studied. Therefore, the present study aimed to investigate the therapeutic potential of curcumin in APE and associated myocardial injury to provide new insights into curcumin as a promising competitive new target for the treatment of APE. A suspension of 12 mg/kg microspheres was injected intravenously into rats. An APE rat model was built. Before modeling, intragastric 100 mg/kg curcumin was given, and/or lentiviral plasmid vector targeting microRNA-145-5p or insulin receptor substrate 1 (IRS1) was injected. Pulmonary artery pressure was measured to assess right ventricular systolic pressure (RVSP). Hematoxylin and eosin (H&E) staining was performed on liver tissues and myocardial tissues of APE rats. TUNEL (terminal deoxynucleotidyl transferase biotin-dUTP nick end labeling) staining and immunohistochemical (IHC) staining were conducted to measure apoptosis and CyPA-CD147 expression in the myocardium, respectively. Inflammatory indices interleukin-1beta (IL-1ß), interleukin-6 (IL-6) and tumor necrosis factor alpha (TNF-α) were measured by ELISA in cardiac tissues. RT-qPCR and Western blot were performed to determine the expression levels of related genes. In addition, by dual luciferase reporter assay and RIP assay, the relationship between microRNA-145-5p and insulin receptor substrate 1 (IRS1) was confirmed. In results: curcumin improved APE-induced myocardial injury, reduced myocardial tissue edema, and thrombus volume. It attenuated APE-induced myocardial inflammation and apoptosis, as well as reduced lung injury and pulmonary artery pressure. Curcumin promoted microRNA-145-5p expression in APE rat myocardium. MicroRNA-145-5p overexpression protected against APE-induced myocardial injury, and microRNA-145-5p silencing abolished the beneficial effects of curcumin in APE-induced myocardial injury. IRS1 was targeted by microRNA-145-5p. IRS1 silencing attenuated APE-induced myocardial injury, and enhanced therapeutic effect of curcumin on myocardial injury in APE rats. In conclusion, curcumin alleviates myocardial inflammation, apoptosis, and oxidative stress induced by APE by regulating microRNA-145-5p/IRS1 axis.

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.

Diverse cytomegalovirus US11 antagonism and MHC-A evasion strategies reveal a tit-for-tat coevolutionary arms race in hominids.

Recurrent, ancient arms races between viruses and hosts have shaped both host immunological defense strategies as well as viral countermeasures. One such battle is waged by the glycoprotein US11 encoded by the persisting human cytomegalovirus. US11 mediates degradation of major histocompatibility class I (MHC-I) molecules to prevent CD8+ T-cell activation. Here, we studied the consequences of the arms race between US11 and primate MHC-A proteins, leading us to uncover a tit-for-tat coevolution and its impact on MHC-A diversification. We found that US11 spurred MHC-A adaptation to evade viral antagonism: In an ancestor of great apes, the MHC-A A2 lineage acquired a Pro184Ala mutation, which confers resistance against the ancestral US11 targeting strategy. In response, US11 deployed a unique low-complexity region (LCR), which exploits the MHC-I peptide loading complex to target the MHC-A2 peptide-binding groove. In addition, the global spread of the human HLA-A*02 allelic family prompted US11 to employ a superior LCR strategy with an optimally fitting peptide mimetic that specifically antagonizes HLA-A*02. Thus, despite cytomegaloviruses low pathogenic potential, the increasing commitment of US11 to MHC-A has significantly promoted diversification of MHC-A in hominids.

Distilling functional variations for human UGT2B4 upstream region based on selection signals and implications for phenotypes of Neanderthal and Denisovan.

Our previous work identified one region upstream human UGT2B4 (UDP glucuronosyltransferase family 2 member B4) which is associated with breast cancer and under balancing selection. However, the distribution, functional variation and molecular mechanism underlying breast cancer and balancing selection remain unclear. In current study, the two haplotypes with deep divergence are described by analyzing 1000 genomes project data and observed to be with high frequencies in all human populations. Through population genetics analysis and genome annotation, the potential functional region is identified and verified by reporter gene assay. Further mutagenesis indicates that the functional mutations are rs66862535 and rs68096061. Both SNPs can alter the interaction efficiency of transcription factor POU2F1 (POU class 2 homeobox 1). Through chromosome conformation capture, it is identified that the enhancer containing these two SNPs can interact with UGT2B4 promoter. Expression quantitative trait loci analysis indicates that UGT2B4 expression is dependent on the genotype of this locus. The common haplotype in human is lost in four genomes of archaic hominins, which suggests that Neanderthal and Denisovan should present relatively lower UGT2B4 expression and further higher steroid hormone level. This study provides new insight into the contribution of ancient population structure to human phenotypes.

Secondary structure construction and molecular identification of rRNA 18S V4 region E23-5-E23-6 of parasitic lice of Hominidae.

The parasitic lice of Hominidae are a class of blood-sucking insects, having a large fragment expansion region in ribosome 18S V4 region. In this study, the value of the E23-5-E23-6 stem-loop structure in the insertion region for molecular identification of lice were explored through motif analysis and secondary structure construction. Five pubic lice samples from China were morphologically identified, and primers for the rRNA 18S V4 region were designed for molecular identification. The V4 sequence of the parasitic lice of Hominidae was retrieved from GenBank for sequence analysis. The five samples were identified as pubic lice based on V4 region, which were of the same specie but geographically different from Australian strains in Genbank, with the identity of 99.06-99.46%. Compared with the human lice, both the chimpanzee lice and pubic lice had large indel fragments in the V4 region. Comparison results showed that Muscle and MAFFT had better alignment and phylogeny results than Clustal. The large expansion region, comprising E23-5 and E23-6, was located between E23-4 and E23-7. The V4 secondary structure showed that the stem-loop structures of the lice parasitizing on pubic area, human, and chimpanzee were different in the E23-5 and E23-6, which could effectively distinguish the three parasitic lice and divide the human lice into five genotypes. This is suitable not only for the identification of three lice species in higher taxonomic ranks but also for genotype identification of human lice in lower taxonomic ranks. The difference between the stem-loop structure is more intuitive than that between the primary sequences.

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.

Ancestral APOBEC3B Nuclear Localization Is Maintained in Humans and Apes and Altered in Most Other Old World Primate Species.

APOBEC3B is an innate immune effector enzyme capable of introducing mutations in viral genomes through DNA cytosine-to-uracil editing. Recent studies have shown that gamma-herpesviruses, such as Epstein-Barr virus (EBV), have evolved a potent APOBEC3B neutralization mechanism to protect lytic viral DNA replication intermediates in the nuclear compartment. APOBEC3B is additionally unique as the only human DNA deaminase family member that is constitutively nuclear. Nuclear localization has therefore been inferred to be essential for innate antiviral function. Here, we combine evolutionary, molecular, and cell biology approaches to address whether nuclear localization is a conserved feature of APOBEC3B in primates. Despite the relatively recent emergence of APOBEC3B approximately 30 to 40 million years ago (MYA) in Old World primates by genetic recombination (after the split from the New World monkey lineage 40 to 50 MYA), we find that the hallmark nuclear localization of APOBEC3B shows variability. For instance, although human and several nonhuman primate APOBEC3B enzymes are predominantly nuclear, rhesus macaque and other Old World primate APOBEC3B proteins are clearly cytoplasmic or cell wide. A series of human/rhesus macaque chimeras and mutants combined to map localization determinants to the N-terminal half of the protein with residues 15, 19, and 24 proving critical. Ancestral APOBEC3B reconstructed from present-day primate species also shows strong nuclear localization. Together, these results indicate that the ancestral nuclear localization of APOBEC3B is maintained in present-day human and ape proteins, but nuclear localization is not conserved in all Old World monkey species despite a need for antiviral functions in the nuclear compartment. IMPORTANCE APOBEC3 enzymes are single-stranded DNA cytosine-to-uracil deaminases with beneficial roles in antiviral immunity and detrimental roles in cancer mutagenesis. Regarding viral infection, all seven human APOBEC3 enzymes have overlapping roles in restricting virus types that require DNA for replication, including EBV, HIV, human papillomavirus (HPV), and human T-cell leukemia virus (HTLV). Regarding cancer, at least two APOBEC3 enzymes, APOBEC3B and APOBEC3A, are prominent sources of mutation capable of influencing clinical outcomes. Here, we combine evolutionary, molecular, and cell biology approaches to characterize primate APOBEC3B enzymes. We show that nuclear localization is an ancestral property of APOBEC3B that is maintained in present-day human and ape enzymes, but not conserved in other nonhuman primates. This partial mechanistic conservation indicates that APOBEC3B is important for limiting the replication of DNA-based viruses in the nuclear compartment. Understanding these pathogen-host interactions may contribute to the development of future antiviral and antitumor therapies.

Molecular diversity and phenotypic pleiotropy of ancient genomic regulatory loci derived from human endogenous retrovirus type H (HERVH) promoter LTR7 and HERVK promoter LTR5_Hs and their contemporary impacts on pathophysiology of Modern Humans.

Timelines of population-level effects of viruses on humans varied from the evolutionary scale of million years to contemporary spread of viral infections. Correspondingly, these events are exemplified by: (i) emergence of human endogenous retroviruses (HERVs) from ancient germline infections leading to stable integration of viral genomes into human chromosomes; and (ii) wide-spread viral infections reaching a global pandemic state such as the COVID-19 pandemic. Despite significant efforts, understanding of HERV's roles in governance of genomic regulatory networks, their impacts on primate evolution and development of human-specific physiological and pathological phenotypic traits remains limited. Remarkably, present analyses revealed that expression of a dominant majority of genes (1696 of 1944 genes; 87%) constituting high-confidence down-steam regulatory targets of defined HERV loci was significantly altered in cells infected with the SARS-CoV-2 coronavirus, a pathogen causing the global COVID-19 pandemic. This study focused on defined sub-sets of DNA sequences derived from HERVs that are expressed at specific stages of human preimplantation embryogenesis and exert regulatory actions essential for self-renewal and pluripotency. Evolutionary histories of LTR7/HERVH and LTR5_Hs/HERVK were charted based on evidence of the earliest presence and expansion of highly conserved (HC) LTR sequences. Sequence conservation analyses of most recent releases 17 primate species' genomes revealed that LTR7/HERVH have entered germlines of primates in Africa after the separation of the New World Monkey lineage, while LTR5_Hs/HERVK successfully colonized primates' germlines after the segregation of Gibbons' species. Subsequently, both LTR7 and LTR5_Hs undergo a marked ~ fourfold-fivefold expansion in genomes of Great Apes. Timelines of quantitative expansion of both LTR7 and LTR5_Hs loci during evolution of Great Apes appear to replicate the consensus evolutionary sequence of increasing cognitive and behavioral complexities of non-human primates, which seems particularly striking for LTR7 loci and 11 distinct LTR7 subfamilies. Consistent with previous reports, identified in this study, 351 human-specific (HS) insertions of LTR7 (175 loci) and LTR5_Hs (176 loci) regulatory sequences have been linked to genes implicated in establishment and maintenance of naïve and primed pluripotent states and preimplantation embryogenesis phenotypes. Unexpectedly, HS-LTRs manifest regulatory connectivity to genes encoding markers of 12 distinct cells' populations of fetal gonads, as well as genes implicated in physiology and pathology of human spermatogenesis, including Y-linked spermatogenic failure, oligo- and azoospermia. Granular interrogations of genes linked with 11 distinct LTR7 subfamilies revealed that mammalian offspring survival (MOS) genes seem to remain one of consistent regulatory targets throughout ~ 30 MYA of the divergent evolution of LTR7 loci. Differential GSEA of MOS versus non-MOS genes identified clearly discernable dominant enrichment patterns of phenotypic traits affected by MOS genes linked with LTR7 (562 MOS genes) and LTR5_Hs (126 MOS genes) regulatory loci across the large panel of genomics and proteomics databases reflecting a broad spectrum of human physiological and pathological traits. GSEA of LTR7-linked MOS genes identified more than 2200 significantly enriched records of human common and rare diseases and gene signatures of 466 significantly enriched records of Human Phenotype Ontology traits, including Autosomal Dominant (92 genes) and Autosomal Recessive (93 genes) Inheritance. LTR7 regulatory elements appear linked with genes implicated in functional and morphological features of central nervous system, including synaptic transmission and protein-protein interactions at synapses, as well as gene signatures differentially regulated in cells of distinct neurodevelopmental stages and morphologically diverse cell types residing and functioning in human brain. These include Neural Stem/Precursor cells, Radial Glia cells, Bergman Glia cells, Pyramidal cells, Tanycytes, Immature neurons, Interneurons, Trigeminal neurons, GABAergic neurons, and Glutamatergic neurons. GSEA of LTR7-linked genes identified significantly enriched gene sets encoding markers of more than 80 specialized types of neurons and markers of 521 human brain regions, most prominently, subiculum and dentate gyrus. Identification and characterization of 1944 genes comprising high-confidence down-steam regulatory targets of LTR7 and/or LTR5_Hs loci validated and extended these observations by documenting marked enrichments for genes implicated in neoplasm metastasis, intellectual disability, autism, multiple cancer types, Alzheimer's, schizophrenia, and other brain disorders. Overall, genes representing down-stream regulatory targets of ancient retroviral LTRs exert the apparently cooperative and exceedingly broad phenotypic impacts on human physiology and pathology. This is exemplified by altered expression of 93% high-confidence LTR targets in cells infected by contemporary viruses, revealing a convergence of virus-inflicted aberrations on genomic regulatory circuitry governed by ancient retroviral LTR elements and interference with human cells' differentiation programs.

Genome Analyses of Ten New Ape Adenoviruses with Similarity to Human Mastadenovirus C.

The adenoviruses (AdVs) isolated from humans are taxonomically grouped in seven different species in the Mastadenovirus genus (HAdV-A through G). AdVs isolated from apes are often included in one of the human AdV species. Here we describe the sequence analyses of ten new AdVs that are related to the HAdV-C species and that were isolated from healthy western lowland gorillas, bonobos, chimpanzees, and orangutans kept in Dutch zoos. We analyzed these viruses and compared their genome sequences to those of human- and ape-derived AdV sequences in the NCBI GenBank database. Our data demonstrated that the ape-derived viruses clustering to HAdV-C are markedly distinct from the human HAdV-C species in the size and nucleotide composition (%GC) of their genome, differ in the amino-acid sequence of AdV proteins, and have longer RGD-loops in their penton-base proteins. The viruses form three well-separated clades (the human, the gorilla, and the combined group of the bonobo and chimpanzee viruses), and we propose that these should each be given species-level ranks. The Ad-lumc005 AdV isolated from orangutans was found to be very similar to the gorilla AdVs, and bootstrap inference provided evidence of recombination between the orangutan AdV and the gorilla AdVs. This suggests that this virus may not be a genuine orangutan AdV but may have been transferred from a gorilla to an orangutan host.

Human reproduction is regulated by retrotransposons derived from ancient Hominidae-specific viral infections.

Germ cells are essential to pass DNA from one generation to the next. In human reproduction, germ cell development begins with the specification of primordial germ cells (PGCs) and a failure to specify PGCs leads to human infertility. Recent studies have revealed that the transcription factor network required for PGC specification has diverged in mammals, and this has a significant impact on our understanding of human reproduction. Here, we reveal that the Hominidae-specific Transposable Elements (TEs) LTR5Hs, may serve as TEENhancers (TE Embedded eNhancers) to facilitate PGC specification. LTR5Hs TEENhancers become transcriptionally active during PGC specification both in vivo and in vitro with epigenetic reprogramming leading to increased chromatin accessibility, localized DNA demethylation, enrichment of H3K27ac, and occupation of key hPGC transcription factors. Inactivation of LTR5Hs TEENhancers with KRAB mediated CRISPRi has a significant impact on germ cell specification. In summary, our data reveals the essential role of Hominidae-specific LTR5Hs TEENhancers in human germ cell development.

Positive Selection and Enhancer Evolution Shaped Lifespan and Body Mass in Great Apes.

Within primates, the great apes are outliers both in terms of body size and lifespan, since they include the largest and longest-lived species in the order. Yet, the molecular bases underlying such features are poorly understood. Here, we leveraged an integrated approach to investigate multiple sources of molecular variation across primates, focusing on over 10,000 genes, including approximately 1,500 previously associated with lifespan, and additional approximately 9,000 for which an association with longevity has never been suggested. We analyzed dN/dS rates, positive selection, gene expression (RNA-seq), and gene regulation (ChIP-seq). By analyzing the correlation between dN/dS, maximum lifespan, and body mass, we identified 276 genes whose rate of evolution positively correlates with maximum lifespan in primates. Further, we identified five genes, important for tumor suppression, adaptive immunity, metastasis, and inflammation, under positive selection exclusively in the great ape lineage. RNA-seq data, generated from the liver of six species representing all the primate lineages, revealed that 8% of approximately 1,500 genes previously associated with longevity are differentially expressed in apes relative to other primates. Importantly, by integrating RNA-seq with ChIP-seq for H3K27ac (which marks active enhancers), we show that the differentially expressed longevity genes are significantly more likely than expected to be located near a novel "ape-specific" enhancer. Moreover, these particular ape-specific enhancers are enriched for young transposable elements, and specifically SINE-Vntr-Alus. In summary, we demonstrate that multiple evolutionary forces have contributed to the evolution of lifespan and body size in primates.

Growth of Biological Complexity from Prokaryotes to Hominids Reflected in the Human Genome.

The growth of complexity in evolution is a most intriguing phenomenon. Using gene phylostratigraphy, we showed this growth (as reflected in regulatory mechanisms) in the human genome, tracing the path from prokaryotes to hominids. Generally, the different regulatory gene families expanded at different times, yet only up to the Euteleostomi (bony vertebrates). The only exception was the expansion of transcription factors (TF) in placentals; however, we argue that this was not related to increase in general complexity. Surprisingly, although TF originated in the Prokaryota while chromatin appeared only in the Eukaryota, the expansion of epigenetic factors predated the expansion of TF. Signaling receptors, tumor suppressors, oncogenes, and aging- and disease-associated genes (indicating vulnerabilities in terms of complex organization and strongly enrichment in regulatory genes) also expanded only up to the Euteleostomi. The complexity-related gene properties (protein size, number of alternative splicing mRNA, length of untranslated mRNA, number of biological processes per gene, number of disordered regions in a protein, and density of TF-TF interactions) rose in multicellular organisms and declined after the Euteleostomi, and possibly earlier. At the same time, the speed of protein sequence evolution sharply increased in the genes that originated after the Euteleostomi. Thus, several lines of evidence indicate that molecular mechanisms of complexity growth were changing with time, and in the phyletic lineage leading to humans, the most salient shift occurred after the basic vertebrate body plan was fixed with bony skeleton. The obtained results can be useful for evolutionary medicine.

Interspecies variation in hominid gut microbiota controls host gene regulation.

The gut microbiome exhibits extreme compositional variation between hominid hosts. However, it is unclear how this variation impacts host physiology across species and whether this effect can be mediated through microbial regulation of host gene expression in interacting epithelial cells. Here, we characterize the transcriptional response of human colonic epithelial cells in vitro to live microbial communities extracted from humans, chimpanzees, gorillas, and orangutans. We find that most host genes exhibit a conserved response, whereby they respond similarly to the four hominid microbiomes. However, hundreds of host genes exhibit a divergent response, whereby they respond only to microbiomes from specific host species. Such genes are associated with intestinal diseases in humans, including inflammatory bowel disease and Crohn's disease. Last, we find that inflammation-associated microbial species regulate the expression of host genes previously associated with inflammatory bowel disease, suggesting health-related consequences for species-specific host-microbiome interactions across hominids.

Shaping modern human skull through epigenetic, transcriptional and post-transcriptional regulation of the RUNX2 master bone gene.

RUNX2 encodes the master bone transcription factor driving skeletal development in vertebrates, and playing a specific role in craniofacial and skull morphogenesis. The anatomically modern human (AMH) features sequence changes in the RUNX2 locus compared with archaic hominins' species. We aimed to understand how these changes may have contributed to human skull globularization occurred in recent evolution. We compared in silico AMH and archaic hominins' genomes, and used mesenchymal stromal cells isolated from skull sutures of craniosynostosis patients for in vitro functional assays. We detected 459 and 470 nucleotide changes in noncoding regions of the AMH RUNX2 locus, compared with the Neandertal and Denisovan genomes, respectively. Three nucleotide changes in the proximal promoter were predicted to alter the binding of the zinc finger protein Znf263 and long-distance interactions with other cis-regulatory regions. By surface plasmon resonance, we selected nucleotide substitutions in the 3'UTRs able to affect miRNA binding affinity. Specifically, miR-3150a-3p and miR-6785-5p expression inversely correlated with RUNX2 expression during in vitro osteogenic differentiation. The expression of two long non-coding RNAs, AL096865.1 and RUNX2-AS1, within the same locus, was modulated during in vitro osteogenic differentiation and correlated with the expression of specific RUNX2 isoforms. Our data suggest that RUNX2 may have undergone adaptive phenotypic evolution caused by epigenetic and post-transcriptional regulatory mechanisms, which may explain the delayed suture fusion leading to the present-day globular skull shape.

Reintroduction of the archaic variant of NOVA1 in cortical organoids alters neurodevelopment.

The evolutionarily conserved splicing regulator neuro-oncological ventral antigen 1 (NOVA1) plays a key role in neural development and function. NOVA1 also includes a protein-coding difference between the modern human genome and Neanderthal and Denisovan genomes. To investigate the functional importance of an amino acid change in humans, we reintroduced the archaic allele into human induced pluripotent cells using genome editing and then followed their neural development through cortical organoids. This modification promoted slower development and higher surface complexity in cortical organoids with the archaic version of NOVA1 Moreover, levels of synaptic markers and synaptic protein coassociations correlated with altered electrophysiological properties in organoids expressing the archaic variant. Our results suggest that the human-specific substitution in NOVA1, which is exclusive to modern humans since divergence from Neanderthals, may have had functional consequences for our species' evolution.

Deleterious variants in genes regulating mammalian reproduction in Neanderthals, Denisovans and extant humans.

STUDY QUESTION: Were Neanderthals and Denisovans (referred here also as extinct hominidae) carrying deleterious variants in genes regulating reproduction? SUMMARY ANSWER: The majority of extinct hominidae analyzed here, presented a considerable number of deleterious variants per individual in proteins regulating different aspects of reproduction, including gonad and uterine function, and gametogenesis. WHAT IS KNOWN ALREADY: Neanderthals, Denisovans and extant humans were interfertile and hybridized while occupying geographically overlapping areas in Europe and Asia. This is evidenced by the small archaic genome component (average ∼2%) present in non-African extant humans. STUDY DESIGN, SIZE, DURATION: The genome of eight extinct hominidae, together with five human genome databases, plus 44 mothers and 48 fathers (fertile controls), were screened to look for deleterious variants in 1734 protein-coding genes regulating reproduction. PARTICIPANTS/MATERIALS, SETTING, METHODS: Ancient DNA from six Neanderthals and two Denisovans dated between ∼82 000 and 43 000 calibrated years was retrieved from the public European Nucleotide Archive. The hominins analyzed include Altai, Vindija 33.15, 33.19, 33.25 and 33.26, El Sidron 1253, Denisova 3 and 11. Their DNA was analyzed using the CLC Genomics Workbench 12, by mapping overlapping paired-end reads (Illumina, FASTQ files) to the human genome assembly GRCh37 (hg19) (Vindija 33.19, 33.25, 33.26, Denisova 3 and Denisova 11) or by analyzing BAM files (Altai, El Sidron 1253 and Vindija 33.15) (human genome reference, GRCh37 (hg19)). Non-synonymous reproductive variants were classified as deleterious or tolerated (PolyPhen-2 and SIFT analyses) and were compared to deleterious variants obtained from extant human genome databases (Genome Aggregation Database (GnomAD), 1000 Genomes, the Haplotype Map (HapMap), Single Nucleotide Polymorphism Database (dbSNPs)) across different populations. A genetic intersection between extant or extinct DNA variants and other genetic disorders was evaluated by annotating the obtained variants with the Clinical Variant (ClinVar) database. MAIN RESULTS AND THE ROLE OF CHANCE: Among the eight extinct hominidae analyzed, a total of 9650 non-synonymous variants (only coverage ≥20 reads included; frameshift mutations were excluded) in 1734 reproductive protein-coding genes were found, 24% of which were classified as deleterious. The majority (73%) of the deleterious alleles present in extant humans that are shared between extant humans and extinct hominidae were found to be rare (<1%) in extant human populations. A set of 8044 variants were found uniquely in extinct hominidae. At the single-gene level, no extinct individual was found to be homozygous for deleterious variants in genes necessary for gamete recognition and fusion, and no higher chance of embryo-lethality (calculated by Mendelian Genetics) was found upon simulated mating between extant human and extinct hominidae compared to extant human-extant human. However, three of the eight extinct hominidae were found to be homozygous for 48-69 deleterious variants in 55 genes controlling ovarian and uterine functions, or oogenesis (AKAP1, BUB1B, CCDC141, CDC73, DUSP6, ESR1, ESR2, PATL2, PSMC3IP, SEMA3A, WT1 and WNT4). Moreover, we report the distribution of nine Neanderthal variants in genes associated with a human fertility phenotype found in extant human populations, one of which has been associated with polycystic ovarian syndrome and primary congenital glaucoma. LIMITATIONS, REASONS FOR CAUTION: While analyzing archaic DNA, stringent filtering criteria were adopted to screen for deleterious variants in Neanderthals and Denisovans, which could result in missing a number of variants. Such restraints preserve the potential for detection of additional deleterious variants in reproductive proteins in extinct hominidae. WIDER IMPLICATIONS OF THE FINDINGS: This study provides a comprehensive overview of putatively deleterious variants in extant human populations and extinct individuals occurring in 1734 protein-coding genes controlling reproduction and provides the fundaments for future functional studies of extinct variants in human reproduction. STUDY FUNDING/COMPETING INTEREST(S): This study was supported by the Department of Biological Science and by the Office of Research and Sponsored Programs at the University of Tulsa (Faculty Research Grant and Faculty Research Summer Fellowship) to M.A. and the University of Tulsa, Tulsa Undergraduate Research Challenge (TURC) program to E.L.; no conflict of interest to declare. TRIAL REGISTRATION NUMBER: N/A.

Multiple selective sweeps of ancient polymorphisms in and around LTα located in the MHC class III region on chromosome 6.

BACKGROUND: Lymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and near LTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassing LTα and the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations. RESULTS: Here, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating > 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or near LTα occurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes. CONCLUSIONS: Overall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

Soluble PD-L1 generated by endogenous retroelement exaptation is a receptor antagonist.

Immune regulation is a finely balanced process of positive and negative signals. PD-L1 and its receptor PD-1 are critical regulators of autoimmune, antiviral and antitumoural T cell responses. Although the function of its predominant membrane-bound form is well established, the source and biological activity of soluble PD-L1 (sPD-L1) remain incompletely understood. Here, we show that sPD-L1 in human healthy tissues and tumours is produced by exaptation of an intronic LINE-2A (L2A) endogenous retroelement in the CD274 gene, encoding PD-L1, which causes omission of the transmembrane domain and the regulatory sequence in the canonical 3' untranslated region. The alternatively spliced CD274-L2A transcript forms the major source of sPD-L1 and is highly conserved in hominids, but lost in mice and a few related species. Importantly, CD274-L2A-encoded sPD-L1 lacks measurable T cell inhibitory activity. Instead, it functions as a receptor antagonist, blocking the inhibitory activity of PD-L1 bound on cellular or exosomal membranes.

The repertoire of family A-peptide GPCRs in archaic hominins.

Given the importance of G-protein coupled receptors in the regulation of many physiological functions, deciphering the relationships between genotype and phenotype in past and present hominin GPCRs is of main interest to understand the evolutionary process that contributed to the present-day variability in human traits and health. Here, we carefully examined the publicly available genomic and protein sequence databases of the archaic hominins (Neanderthal and Denisova) to draw up the catalog of coding variations in GPCRs for peptide ligands, in comparison with living humans. We then searched in the literature the functional changes, phenotypes and risk of disease possibly associated with the detected variants. Our survey suggests that Neanderthal and Denisovan hominins were likely prone to lower risk of obesity, to enhanced platelet aggregation in response to thrombin, to better response to infection, to less anxiety and aggressiveness and to favorable sociability. While some archaic variants were likely advantageous in the past, they might be responsible for maladaptive disorders today in the context of modern life and/or specific regional distribution. For example, an archaic haplotype in the neuromedin receptor 2 is susceptible to confer risk of diabetic nephropathy in type 1 diabetes in present-day Europeans. Paying attention to the pharmacological properties of some of the archaic variants described in this study may be helpful to understand the variability of therapeutic efficacy between individuals or ethnic groups.