"Outlaw" mutations in quasispecies of SARS-CoV-2 inhibit replication.

The evolution of SARS-CoV-2, the agent of COVID-19, has been remarkable for its high mutation potential, leading to the appearance of variants. Some mutations have never appeared in the published genomes, which represent consensus, or bona fide genomes. Here we tested the hypothesis that mutations that did not appear in consensus genomes were, in fact, as frequent as the mutations that appeared during the various epidemic episodes, but were not expressed because lethal. To identify these mutations, we analysed the genomes of 90 nasopharyngeal samples and the quasispecies determined by next-generation sequencing. Mutations observed in the quasispecies and not in the consensus genomes were considered to be lethal, what we called "outlaw" mutations. Among these mutations, we analysed the 21 most frequent. Eight of these "outlaws" were in the RNA polymerase and we were able to use a structural biology model and molecular dynamics simulations to demonstrate the functional incapacity of these mutated RNA polymerases. Three other mutations affected the spike, a major protein involved in the pathogenesis of COVID-19. Overall, by analysing the SARS-CoV-2 quasispecies obtained during sequencing, this method made it possible to identify "outlaws," showing areas that could potentially become the target of treatments.

Mechanisms of acquisition of the vanA operon among vancomycin-resistant Staphylococcus aureus genomes: The tip of the iceberg?

BACKGROUND: Vancomycin is frequently used as a last line of defence against infections due to multidrug-resistant Staphylococcus aureus (S. aureus). A recent finding described the acquisition of vancomycin-resistant S. aureus strains by the integration of an enterococcal plasmid containing the vanA operon into the S. aureus chromosome via homologous recombination involving a specific integration site called locus L2. METHODS: To characterise all mechanisms of acquisition of vanA, this study analysed the 15 706 S. aureus genomes to look for vanA and described its genetic environment. RESULTS: A complete vanA operon was found in 25 S. aureus strains isolated from 12 patients, including nine co-isolated with vancomycin-resistant Enterococcus strains. VanA was found within transposon Tn1546-like elements on 17 plasmids and eight chromosomes. VanA might be acquired through conjugation of enterococcal and staphylococcal plasmids, transposition of Tn1546 carrying vanA and plasmid integration into the chromosome. Further, L2 was detected in 2087 genomes (13.3%) of S. aureus strains across different continents. Six potential chromosomal hotspots for integration of the entire vanA-containing enterococcal plasmid were identified by homologous recombination via L2. CONCLUSIONS: These findings suggest that the recently described scenario in a New York patient could be reproduced anywhere. Surveillance of this possibility is mandatory, especially in patients with vancomycin-resistant Enterococcus infection or colonisation.

Role of SARS-CoV-2 mutations in the evolution of the COVID-19 pandemic.

OBJECTIVES: The SARS-CoV-2 pandemic and large-scale genomic surveillance provided an exceptional opportunity to analyze mutations that appeared over three years in viral genomes. Here we studied mutations and their epidemic consequences for SARS-CoV-2 genomes from our center. METHODS: We analyzed 61,397 SARS-CoV-2 genomes we sequenced from respiratory samples for genomic surveillance. Mutations frequencies were calculated using Nextclade, Microsoft Excel, and an in-house Python script. RESULTS: A total of 22,225 nucleotide mutations were identified, 220 (1.0%) being each at the root of ≥836 genomes, classifying mutations as 'hyperfertile'. Two seeded the European pandemic: P323L in RNA polymerase, associated with an increased mutation rate, and D614G in spike that improved fitness. Most 'hyperfertile' mutations occurred in areas not predicted with increased virulence. Their mean number was 8±6 (0-22) per 1000 nucleotides per gene. They were 3.7-times more frequent in accessory than informational genes (13.8 versus 3.7/1000 nucleotides). Particularly, they were 4.1-times more frequent in ORF8 than in the RNA polymerase gene. Interestingly, stop codons were present in 97 positions, almost only in accessory genes, including ORF8 (21/100 codons). CONCLUSIONS: most 'hyperfertile' mutations did not predict emergence of a new epidemic, and some were stop codons indicating the existence of so-named 'non-virulence' genes.

The IL-17 pathway intertwines with neurotrophin and TLR/IL-1R pathways since its domain shuffling origin.

The IL-17 pathway displays remarkably diverse functional modes between different subphyla, classes, and even orders, yet its driving factors remains elusive. Here, we demonstrate that the IL-17 pathway originated through domain shuffling between a Toll-like receptor (TLR)/IL-1R pathway and a neurotrophin-RTK (receptor-tyrosine-kinase) pathway (a Trunk-Torso pathway). Unlike other new pathways that evolve independently, the IL-17 pathway remains intertwined with its donor pathways throughout later evolution. This intertwining not only influenced the gains and losses of domains and components in the pathway but also drove the diversification of the pathway's functional modes among animal lineages. For instance, we reveal that the crustacean female sex hormone, a neurotrophin inducing sex differentiation, could interact with IL-17Rs and thus be classified as true IL-17s. Additionally, the insect prothoracicotropic hormone, a neurotrophin initiating ecdysis in Drosophila by binding to Torso, could bind to IL-17Rs in other insects. Furthermore, IL-17R and TLR/IL-1R pathways maintain crosstalk in amphioxus and zebrafish. Moreover, the loss of the Death domain in the pathway adaptor connection to IκB kinase and stress-activated protein kinase (CIKSs) dramatically reduced their abilities to activate nuclear factor-kappaB (NF-κB) and activator protein 1 (AP-1) in amphioxus and zebrafish. Reinstating this Death domain not only enhanced NF-κB/AP-1 activation but also strengthened anti-bacterial immunity in zebrafish larvae. This could explain why the mammalian IL-17 pathway, whose CIKS also lacks Death, is considered a weak signaling activator, relying on synergies with other pathways. Our findings provide insights into the functional diversity of the IL-17 pathway and unveil evolutionary principles that could govern the pathway and be used to redesign and manipulate it.

Resistance-associated mutations to the anti-SARS-CoV-2 agent nirmatrelvir: Selection not induction.

Mutations associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resistance to antiprotease nirmatrelvir were reported. We aimed to detect them in SARS-CoV-2 genomes and quasispecies retrieved in our institute before drug availability in January 2022 and to analyze the impact of mutations on protease (3CLpro) structure. We sought for 38 3CLpro nirmatrelvir resistance mutations in a set of 62 673 SARS-CoV-2 genomes obtained in our institute from respiratory samples collected between 2020 and 2023 and for these mutations in SARS-CoV-2 quasispecies for 90 samples collected in 2020, using Python. SARS-CoV-2 protease with major mutation E166V was generated with Swiss Pdb Viewer and Molegro Molecular Viewer. We detected 22 (58%) of the resistance-associated mutations in 417 (0.67%) of the genomes analyzed; 325 (78%) of these genomes had been obtained from samples collected in 2020-2021. APOBEC signatures were found for 12/22 mutations. We also detected among viral quasispecies from 90 samples some minority reads harboring any of 15 nirmatrelvir resistance mutations, including E166V. Also, we predicted that E166V has a very limited effect on 3CLpro structure but may prevent drug attachment. Thus, we evidenced that mutations associated with nirmatrelvir resistance pre-existed in SARS-CoV-2 before drug availability. These findings further warrant SARS-CoV-2 genomic surveillance and SARS-CoV-2 quasispecies characterization.

From viral democratic genomes to viral wild bunch of quasispecies.

The tremendous majority of RNA genomes from pathogenic viruses analyzed and deposited in databases are consensus or "democratic" genomes. They represent the genomes most frequently found in the clinical samples of patients but do not account for the huge genetic diversity of coexisting genomes, which is better described as quasispecies. A viral quasispecies is defined as the dynamic distribution of nonidentical but closely related mutants, variants, recombinant, or reassortant viral genomes. Viral quasispecies have collective behavior and dynamics and are the subject of internal interactions that comprise interference, complementation, or cooperation. In the setting of SARS-CoV-2 infection, intrahost SARS-CoV-2 genetic diversity was recently notably reported for immunocompromised, chronically infected patients, for patients treated with monoclonal antibodies targeting the viral spike protein, and for different body compartments of a single patient. A question that deserves attention is whether such diversity is generated postinfection from a clonal genome in response to selection pressure or is already present at the time of infection as a quasispecies. In the present review, we summarize the data supporting that hosts are infected by a "wild bunch" of viruses rather than by multiple virions sharing the same genome. Each virion in the "wild bunch" may have different virulence and tissue tropisms. As the number of viruses replicated during host infections is huge, a viral quasispecies at any time of infection is wide and is also influenced by host-specific selection pressure after infection, which accounts for the difficulty in deciphering and predicting the appearance of more fit variants and the evolution of epidemics of novel RNA viruses.

Insights into RAG Evolution from the Identification of "Missing Link" Family A RAGL Transposons.

A series of "molecular domestication" events are thought to have converted an invertebrate RAG-like (RAGL) transposase into the RAG1-RAG2 (RAG) recombinase, a critical enzyme for adaptive immunity in jawed vertebrates. The timing and order of these events are not well understood, in part because of a dearth of information regarding the invertebrate RAGL-A transposon family. In contrast to the abundant and divergent RAGL-B transposon family, RAGL-A most closely resembles RAG and is represented by a single orphan RAG1-like (RAG1L) gene in the genome of the hemichordate Ptychodera flava (PflRAG1L-A). Here, we provide evidence for the existence of complete RAGL-A transposons in the genomes of P. flava and several echinoderms. The predicted RAG1L-A and RAG2L-A proteins encoded by these transposons intermingle sequence features of jawed vertebrate RAG and RAGL-B transposases, leading to a prediction of DNA binding, catalytic, and transposition activities that are a hybrid of RAG and RAGL-B. Similarly, the terminal inverted repeats (TIRs) of the RAGL-A transposons combine features of both RAGL-B transposon TIRs and RAG recombination signal sequences. Unlike all previously described RAG2L proteins, RAG2L-A proteins contain an acidic hinge region, which we demonstrate is capable of efficiently inhibiting RAG-mediated transposition. Our findings provide evidence for a critical intermediate in RAG evolution and argue that certain adaptations thought to be specific to jawed vertebrates (e.g. the RAG2 acidic hinge) actually arose in invertebrates, thereby focusing attention on other adaptations as the pivotal steps in the completion of RAG domestication in jawed vertebrates.

The return of the "Mistigri" (virus adaptative gain by gene loss) through the SARS-CoV-2 XBB.1.5 chimera that predominated in 2023.

Severe acute respiratory syndrome coronavirus 2 XBB.1.5 is the first recombinant lineage to predominate at the country and global scales. Very interestingly, like the Marseille-4B subvariant (or B.1.160) and the pandemic variant B.1.1.7 (or Alpha) previously, it has its ORF8 gene inactivated by a stop codon. We aimed here to study the distribution of stop codons in ORF8 of XBB.1.5 and non-XBB.1.5 genomes. We identified that a stop codon was present at 89 (74%) ORF8 codons in ≥1 of 15 222 404 genomes available in GISAID. The mean proportion of genomes with a stop codon per codon was 0.11% (range, 0%-7.8%). In addition, a stop codon was detected at 15 (12%) codons in at least 1000 genomes. These 15 codons are notably located on seven stem-loop hairpin regions and in the signal peptide region for the case of the XBB.1.5 lineage (codon 8). Thus, it is very likely that stop codons in ORF8 gene contributed on at least three occasions and independently during the pandemic to the evolutionary success of a lineage that became transiently predominant. Such association of gene loss with evolutionary success, which suits the recently described Mistigri rule, is an important biological phenomenon very unknown in virology while largely described in cellular organisms.

Integration of the immune memory into the pathogen-driven MHC polymorphism hypothesis.

Major histocompatibility complex (MHC) genes (referred to as human leukocyte antigen or HLA in humans) are a key component of vertebrate immune systems, coding for proteins which present antigens to T-cells. These genes are outstanding in their degree of polymorphism, with important consequences for human and animal health. The polymorphism is thought to arise from selection pressures imposed by pathogens on MHC allomorphs, which differ in their antigen-binding capacity. However, the existing theory has not considered MHC selection in relation to the formation of immune memory. In this paper, we argue that this omission limits our understanding of the evolution of MHC polymorphism and its role in disease. We review recent evidence that has emerged from the massive research effort related to the SARS-CoV-2 pandemics, and which provides new evidence for the role of MHC in shaping immune memory. We then discuss why the inclusion of immune memory within the existing theory may have non-trivial consequence for our understanding of the evolution of MHC polymorphism. Finally, we will argue that neglecting immune memory hinders our interpretation of empirical findings, and postulate that future studies focusing on pathogen-driven MHC selection would benefit from stratifying the available data according to the history of infection (and vaccination, if relevant).

Insights into RAG evolution from the identification of "missing link" family A RAGL transposons.

A series of "molecular domestication" events are thought to have converted an invertebrate RAG-like (RAGL) transposase into the RAG1-RAG2 (RAG) recombinase, a critical enzyme for adaptive immunity in jawed vertebrates. The timing and order of these events is not well understood, in part because of a dearth of information regarding the invertebrate RAGL-A transposon family. In contrast to the abundant and divergent RAGL-B transposon family, RAGL-A most closely resembles RAG and is represented by a single orphan RAG1-like (RAG1L) gene in the genome of the hemichordate Ptychodera flava (PflRAG1L-A). Here, we provide evidence for the existence of complete RAGL-A transposons in the genomes of P. flava and several echinoderms. The predicted RAG1L-A and RAG2L-A proteins encoded by these transposons intermingle sequence features of jawed vertebrate RAG and RAGL-B transposases, leading to a prediction of DNA binding, catalytic, and transposition activities that are a hybrid of RAG and RAGL-B. Similarly, the terminal inverted repeats (TIRs) of the RAGL-A transposons combine features of both RAGL-B transposon TIRs and RAG recombination signal sequences. Unlike all previously described RAG2L proteins, PflRAG2L-A and echinoderm RAG2L-A contain an acidic hinge region, which we demonstrate is capable of efficiently inhibiting RAG-mediated transposition. Our findings provide evidence for a critical intermediate in RAG evolution and argue that certain adaptations thought to be specific to jawed vertebrates (e.g., the RAG2 acidic hinge) actually arose in invertebrates, thereby focusing attention on other adaptations as the pivotal steps in the completion of RAG domestication in jawed vertebrates.

Metallo-Beta-Lactamase-like Encoding Genes in Candidate Phyla Radiation: Widespread and Highly Divergent Proteins with Potential Multifunctionality.

The Candidate Phyla Radiation (CPR) was found to harbor a vast repertoire of genes encoding for enzymes with potential antibiotic resistance activity. Among these, as many as 3349 genes were predicted in silico to contain a metallo-beta-lactamase-like (MBL-like) fold. These proteins were subject to an in silico functional characterization by comparing their protein profiles (presence/absence of conserved protein domains) to other MBLs, including 24 already expressed in vitro, along with those of the beta-lactamase database (BLDB) (n = 761). The sequence similarity network (SSN) was then used to predict the functional clusters of CPR MBL-like sequences. Our findings showed that CPR MBL-like sequences were longer and more diverse than bacterial MBL sequences, with a high content of functional domains. Most CPR MBL-like sequences did not show any SSN connectivity with expressed MBLs, indicating the presence of many potential, yet unidentified, functions in CPR. In conclusion, CPR was shown to have many protein functions and a large sequence variability of MBL-like folds, exceeding all known MBLs. Further experimental and evolutionary studies of this superfamily of hydrolyzing enzymes are necessary to illustrate their functional annotation, origin, and expansion for adaptation or specialization within a given niche or compared to a specific substrate.

Origin, Diversity, and Multiple Roles of Enzymes with Metallo-ß-Lactamase Fold from Different Organisms.

ß-lactamase enzymes have generated significant interest due to their ability to confer resistance to the most commonly used family of antibiotics in human medicine. Among these enzymes, the class B ß-lactamases are members of a superfamily of metallo-ß-lactamase (MßL) fold proteins which are characterised by conserved motifs (i.e., HxHxDH) and are not only limited to bacteria. Indeed, as the result of several barriers, including low sequence similarity, default protein annotation, or untested enzymatic activity, MßL fold proteins have long been unexplored in other organisms. However, thanks to search approaches which are more sensitive compared to classical Blast analysis, such as the use of common ancestors to identify distant homologous sequences, we are now able to highlight their presence in different organisms including Bacteria, Archaea, Nanoarchaeota, Asgard, Humans, Giant viruses, and Candidate Phyla Radiation (CPR). These MßL fold proteins are multifunctional enzymes with diverse enzymatic or non-enzymatic activities of which, at least thirteen activities have been reported such as ß-lactamase, ribonuclease, nuclease, glyoxalase, lactonase, phytase, ascorbic acid degradation, anti-cancer drug degradation, or membrane transport. In this review, we (i) discuss the existence of MßL fold enzymes in the different domains of life, (ii) present more suitable approaches to better investigating their homologous sequences in unsuspected sources, and (iii) report described MßL fold enzymes with demonstrated enzymatic or non-enzymatic activities.

The emergence, spread and vanishing of a French SARS-CoV-2 variant exemplifies the fate of RNA virus epidemics and obeys the Mistigri rule.

The nature and dynamics of mutations associated with the emergence, spread, and vanishing of SARS-CoV-2 variants causing successive waves are complex. We determined the kinetics of the most common French variant ("Marseille-4") for 10 months since its onset in July 2020. Here, we analyzed and classified into subvariants and lineages 7453 genomes obtained by next-generation sequencing. We identified two subvariants, Marseille-4A, which contains 22 different lineages of at least 50 genomes, and Marseille-4B. Their average lifetime was 4.1 ± 1.4 months, during which 4.1 ± 2.6 mutations accumulated. Growth rate was 0.079 ± 0.045, varying from 0.010 to 0.173. Most of the lineages exhibited a bell-shaped distribution. Several beneficial mutations at unpredicted sites initiated a new outbreak, while the accumulation of other mutations resulted in more viral heterogenicity, increased diversity and vanishing of the lineages. Marseille-4B emerged when the other Marseille-4 lineages vanished. Its ORF8 gene was knocked out by a stop codon, as reported in SARS-CoV-2 of mink and in the Alpha variant. This subvariant was associated with increased hospitalization and death rates, suggesting that ORF8 is a nonvirulence gene. We speculate that the observed heterogenicity of a lineage may predict the end of the outbreak.

Is it time to switch to a formulation other than the live attenuated poliovirus vaccine to prevent poliomyelitis?

The polioviruses (PVs) are mainly transmitted by direct contact with an infected person through the fecal-oral route and respiratory secretions (or more rarely via contaminated water or food) and have a primary tropism for the gut. After their replication in the gut, in rare cases (far less than 1% of the infected individuals), PVs can spread to the central nervous system leading to flaccid paralysis, which can result in respiratory paralysis and death. By the middle of the 20th century, every year the wild polioviruses (WPVs) are supposed to have killed or paralyzed over half a million people. The introduction of the oral poliovirus vaccines (OPVs) through mass vaccination campaigns (combined with better application of hygiene measures), was a success story which enabled the World Health Organization (WHO) to set the global eradication of poliomyelitis as an objective. However this strategy of viral eradication has its limits as the majority of poliomyelitis cases today arise in individuals infected with circulating vaccine-derived polioviruses (cVDPVs) which regain pathogenicity following reversion or recombination. In recent years (between January 2018 and May 2023), the WHO recorded 8.8 times more cases of polio which were linked to the attenuated OPV vaccines (3,442 polio cases after reversion or recombination events) than cases linked to a WPV (390 cases). Recent knowledge of the evolution of RNA viruses and the exchange of genetic material among biological entities of the intestinal microbiota, call for a reassessment of the polio eradication vaccine strategies.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications.

BACKGROUND: Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates. RESULTS: We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements. CONCLUSIONS: In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

Species and population genomic differentiation in Pocillopora corals (Cnidaria, Hexacorallia).

Correctly delimiting species and populations is a prerequisite for studies of connectivity, adaptation and conservation. Genomic data are particularly useful to test species differentiation for organisms with few informative morphological characters or low discrimination of cytoplasmic markers, as in Scleractinians. Here we applied Restriction site Associated DNA sequencing (RAD-sequencing) to the study of species differentiation and genetic structure in populations of Pocillopora spp. from Oman and French Polynesia, with the objectives to test species hypotheses, and to study the genetic structure among sampling sites within species. We focused here on coral colonies morphologically similar to P. acuta (damicornis type ß). We tested the impact of different filtering strategies on the stability of the results. The main genetic differentiation was observed between samples from Oman and French Polynesia. These samples corresponded to different previously defined primary species hypotheses (PSH), i.e., PSHs 12 and 13 in Oman, and PSH 5 in French Polynesia. In Oman, we did not observe any clear differentiation between the two putative species PSH 12 and 13, nor between sampling sites. In French Polynesia, where a single species hypothesis was studied, there was no differentiation between sites. Our analyses allowed the identification of clonal lineages in Oman and French Polynesia. The impact of clonality on genetic diversity is discussed in light of individual-based simulations.

RadA, a Key Gene of the Circadian Rhythm of Escherichia coli.

Circadian rhythms are present in almost all living organisms, and their activity relies on molecular clocks. In prokaryotes, a functional molecular clock has been defined only in cyanobacteria. Here, we investigated the presence of circadian rhythms in non-cyanobacterial prokaryotes. The bioinformatic approach was used to identify a homologue of KaiC (circadian gene in cyanobacteria) in Escherichia coli. Then, strains of E. coli (wild type and mutants) were grown on blood agar, and sampling was made every 3 h for 24 h at constant conditions. Gene expression was determined by qRT-PCR, and the rhythmicity was analyzed using the Cosinor model. We identified RadA as a KaiC homologue in E. coli. Expression of radA showed a circadian rhythm persisting at least 3 days, with a peak in the morning. The circadian expression of other E. coli genes was also observed. Gene circadian oscillations were lost in radA mutants of E. coli. This study provides evidence of molecular clock gene expression in E. coli with a circadian rhythm. Such a finding paves the way for new perspectives in antibacterial treatment.

COVID-19 Pandemic: Escape of Pathogenic Variants and MHC Evolution.

We propose a new hypothesis that explains the maintenance and evolution of MHC polymorphism. It is based on two phenomena: the constitution of the repertoire of naive T lymphocytes and the evolution of the pathogen and its impact on the immune memory of T lymphocytes. Concerning the latter, pathogen evolution will have a different impact on reinfection depending on the MHC allomorph. If a mutation occurs in a given region, in the case of MHC allotypes, which do not recognize the peptide in this region, the mutation will have no impact on the memory repertoire. In the case where the MHC allomorph binds to the ancestral peptides and not to the mutated peptide, that individual will have a higher chance of being reinfected. This difference in fitness will lead to a variation of the allele frequency in the next generation. Data from the SARS-CoV-2 pandemic already support a significant part of this hypothesis and following up on these data may enable it to be confirmed. This hypothesis could explain why some individuals after vaccination respond less well than others to variants and leads to predict the probability of reinfection after a first infection depending upon the variant and the HLA allomorph.

'Cannibalism' of exogenous DNA sequences: The ancestral form of adaptive immunity which entails recognition of danger.

Adaptive immunity is a sophisticated form of immune response capable of retaining the molecular memory of a very great diversity of target antigens (epitopes) as non-self. It is capable of reactivating itself upon a second encounter with an immunoglobulin or T-cell receptor antigen-binding site with a known epitope that had previously primed the host immune system. It has long been considered that adaptive immunity is a highly evolved form of non-self recognition that appeared quite late in speciation and complemented a more generalist response called innate immunity. Innate immunity offers a relatively non-specific defense (although mediated by sensors that could specifically recognize virus or bacteria compounds) and which does not retain a memory of the danger. But this notion of recent acquisition of adaptive immunity is challenged by the fact that another form of specific recognition mechanisms already existed in prokaryotes that may be able to specifically auto-protect against external danger. This recognition mechanism can be considered a primitive form of specific (adaptive) non-self recognition. It is based on the fact that many archaea and bacteria use a genome editing system that confers the ability to appropriate viral DNA sequences allowing prokaryotes to prevent host damage through a mechanism very similar to adaptive immunity. This is indistinctly called, 'endogenization of foreign DNA' or 'viral DNA predation' or, more pictorially 'DNA cannibalism'. For several years evidence has been accumulating, highlighting the crucial role of endogenization of foreign DNA in the fundamental processes related to adaptive immunity and leading to a change in the dogma that adaptive immunity appeared late in speciation.

HLA-H*02:07 Is a Membrane-Bound Ligand of Denisovan Origin That Protects against Lysis by Activated Immune Effectors.

The biological relevance of genes initially categorized as "pseudogenes" is slowly emerging, notably in innate immunity. In the HLA region on chromosome 6, HLA-H is one such pseudogene; yet, it is transcribed, and its variation is associated with immune properties. Furthermore, two HLA-H alleles, H*02:07 and H*02:14, putatively encode a complete, membrane-bound HLA protein. Here we thus hypothesized that HLA-H contributes to immune homeostasis similarly to tolerogenic molecules HLA-G, -E, and -F. We tested if HLA-H*02:07 encodes a membrane-bound protein that can inhibit the cytotoxicity of effector cells. We used an HLA-null human erythroblast cell line transduced with HLA-H*02:07 cDNA to demonstrate that HLA-H*02:07 encodes a membrane-bound protein. Additionally, using a cytotoxicity assay, our results support that K562 HLA-H*02:07 inhibits human effector IL-2-activated PBMCs and human IL-2-independent NK92-MI cell line activity. Finally, through in silico genotyping of the Denisovan genome and haplotypic association with Denisovan-derived HLA-A*11, we also show that H*02:07 is of archaic origin. Hence, admixture with archaic humans brought a functional HLA-H allele into modern European and Asian populations.