Biallelic human SHARPIN loss of function induces autoinflammation and immunodeficiency.

The linear ubiquitin assembly complex (LUBAC) consists of HOIP, HOIL-1 and SHARPIN and is essential for proper immune responses. Individuals with HOIP and HOIL-1 deficiencies present with severe immunodeficiency, autoinflammation and glycogen storage disease. In mice, the loss of Sharpin leads to severe dermatitis due to excessive keratinocyte cell death. Here, we report two individuals with SHARPIN deficiency who manifest autoinflammatory symptoms but unexpectedly no dermatological problems. Fibroblasts and B cells from these individuals showed attenuated canonical NF-κB responses and a propensity for cell death mediated by TNF superfamily members. Both SHARPIN-deficient and HOIP-deficient individuals showed a substantial reduction of secondary lymphoid germinal center B cell development. Treatment of one SHARPIN-deficient individual with anti-TNF therapies led to complete clinical and transcriptomic resolution of autoinflammation. These findings underscore the critical function of the LUBAC as a gatekeeper for cell death-mediated immune dysregulation in humans.

Antibody Lineages with Vaccine-Induced Antigen-Binding Hotspots Develop Broad HIV Neutralization.

The vaccine-mediated elicitation of antibodies (Abs) capable of neutralizing diverse HIV-1 strains has been a long-standing goal. To understand how broadly neutralizing antibodies (bNAbs) can be elicited, we identified, characterized, and tracked five neutralizing Ab lineages targeting the HIV-1-fusion peptide (FP) in vaccinated macaques over time. Genetic and structural analyses revealed two of these lineages to belong to a reproducible class capable of neutralizing up to 59% of 208 diverse viral strains. B cell analysis indicated each of the five lineages to have been initiated and expanded by FP-carrier priming, with envelope (Env)-trimer boosts inducing cross-reactive neutralization. These Abs had binding-energy hotspots focused on FP, whereas several FP-directed Abs induced by immunization with Env trimer-only were less FP-focused and less broadly neutralizing. Priming with a conserved subregion, such as FP, can thus induce Abs with binding-energy hotspots coincident with the target subregion and capable of broad neutralization.

Vaccine-Induced Antibodies that Neutralize Group 1 and Group 2 Influenza A Viruses.

Antibodies capable of neutralizing divergent influenza A viruses could form the basis of a universal vaccine. Here, from subjects enrolled in an H5N1 DNA/MIV-prime-boost influenza vaccine trial, we sorted hemagglutinin cross-reactive memory B cells and identified three antibody classes, each capable of neutralizing diverse subtypes of group 1 and group 2 influenza A viruses. Co-crystal structures with hemagglutinin revealed that each class utilized characteristic germline genes and convergent sequence motifs to recognize overlapping epitopes in the hemagglutinin stem. All six analyzed subjects had sequences from at least one multidonor class, and-in half the subjects-multidonor-class sequences were recovered from >40% of cross-reactive B cells. By contrast, these multidonor-class sequences were rare in published antibody datasets. Vaccination with a divergent hemagglutinin can thus increase the frequency of B cells encoding broad influenza A-neutralizing antibodies. We propose the sequence signature-quantified prevalence of these B cells as a metric to guide universal influenza A immunization strategies.

Maturation Pathway from Germline to Broad HIV-1 Neutralizer of a CD4-Mimic Antibody.

Antibodies with ontogenies from VH1-2 or VH1-46-germline genes dominate the broadly neutralizing response against the CD4-binding site (CD4bs) on HIV-1. Here, we define with longitudinal sampling from time-of-infection the development of a VH1-46-derived antibody lineage that matured to neutralize 90% of HIV-1 isolates. Structures of lineage antibodies CH235 (week 41 from time-of-infection, 18% breadth), CH235.9 (week 152, 77%), and CH235.12 (week 323, 90%) demonstrated the maturing epitope to focus on the conformationally invariant portion of the CD4bs. Similarities between CH235 lineage and five unrelated CD4bs lineages in epitope focusing, length-of-time to develop breadth, and extraordinary level of somatic hypermutation suggested commonalities in maturation among all CD4bs antibodies. Fortunately, the required CH235-lineage hypermutation appeared substantially guided by the intrinsic mutability of the VH1-46 gene, which closely resembled VH1-2. We integrated our CH235-lineage findings with a second broadly neutralizing lineage and HIV-1 co-evolution to suggest a vaccination strategy for inducing both lineages.

Maturation and Diversity of the VRC01-Antibody Lineage over 15 Years of Chronic HIV-1 Infection.

HIV-1-neutralizing antibodies develop in most HIV-1-infected individuals, although highly effective antibodies are generally observed only after years of chronic infection. Here, we characterize the rate of maturation and extent of diversity for the lineage that produced the broadly neutralizing antibody VRC01 through longitudinal sampling of peripheral B cell transcripts over 15 years and co-crystal structures of lineage members. Next-generation sequencing identified VRC01-lineage transcripts, which encompassed diverse antibodies organized into distinct phylogenetic clades. Prevalent clades maintained characteristic features of antigen recognition, though each evolved binding loops and disulfides that formed distinct recognition surfaces. Over the course of the study period, VRC01-lineage clades showed continuous evolution, with rates of ∼2 substitutions per 100 nucleotides per year, comparable to that of HIV-1 evolution. This high rate of antibody evolution provides a mechanism by which antibody lineages can achieve extraordinary diversity and, over years of chronic infection, develop effective HIV-1 neutralization.

Structural Repertoire of HIV-1-Neutralizing Antibodies Targeting the CD4 Supersite in 14 Donors.

The site on the HIV-1 gp120 glycoprotein that binds the CD4 receptor is recognized by broadly reactive antibodies, several of which neutralize over 90% of HIV-1 strains. To understand how antibodies achieve such neutralization, we isolated CD4-binding-site (CD4bs) antibodies and analyzed 16 co-crystal structures -8 determined here- of CD4bs antibodies from 14 donors. The 16 antibodies segregated by recognition mode and developmental ontogeny into two types: CDR H3-dominated and VH-gene-restricted. Both could achieve greater than 80% neutralization breadth, and both could develop in the same donor. Although paratope chemistries differed, all 16 gp120-CD4bs antibody complexes showed geometric similarity, with antibody-neutralization breadth correlating with antibody-angle of approach relative to the most effective antibody of each type. The repertoire for effective recognition of the CD4 supersite thus comprises antibodies with distinct paratopes arrayed about two optimal geometric orientations, one achieved by CDR H3 ontogenies and the other achieved by VH-gene-restricted ontogenies.

The genome of the colonial hydroid Hydractinia reveals that their stem cells use a toolkit of evolutionarily shared genes with all animals.

Hydractinia is a colonial marine hydroid that shows remarkable biological properties, including the capacity to regenerate its entire body throughout its lifetime, a process made possible by its adult migratory stem cells, known as i-cells. Here, we provide an in-depth characterization of the genomic structure and gene content of two Hydractinia species, Hydractinia symbiolongicarpus and Hydractinia echinata, placing them in a comparative evolutionary framework with other cnidarian genomes. We also generated and annotated a single-cell transcriptomic atlas for adult male H. symbiolongicarpus and identified cell-type markers for all major cell types, including key i-cell markers. Orthology analyses based on the markers revealed that Hydractinia's i-cells are highly enriched in genes that are widely shared amongst animals, a striking finding given that Hydractinia has a higher proportion of phylum-specific genes than any of the other 41 animals in our orthology analysis. These results indicate that Hydractinia's stem cells and early progenitor cells may use a toolkit shared with all animals, making it a promising model organism for future exploration of stem cell biology and regenerative medicine. The genomic and transcriptomic resources for Hydractinia presented here will enable further studies of their regenerative capacity, colonial morphology, and ability to distinguish self from nonself.

A Neutralizing Antibody Recognizing Primarily N-Linked Glycan Targets the Silent Face of the HIV Envelope.

Virtually the entire surface of the HIV-1-envelope trimer is recognized by neutralizing antibodies, except for a highly glycosylated region at the center of the "silent face" on the gp120 subunit. From an HIV-1-infected donor, #74, we identified antibody VRC-PG05, which neutralized 27% of HIV-1 strains. The crystal structure of the antigen-binding fragment of VRC-PG05 in complex with gp120 revealed an epitope comprised primarily of N-linked glycans from N262, N295, and N448 at the silent face center. Somatic hypermutation occurred preferentially at antibody residues that interacted with these glycans, suggesting somatic development of glycan recognition. Resistance to VRC-PG05 in donor #74 involved shifting of glycan-N448 to N446 or mutation of glycan-proximal residue E293. HIV-1 neutralization can thus be achieved at the silent face center by glycan-recognizing antibody; along with other known epitopes, the VRC-PG05 epitope completes coverage by neutralizing antibody of all major exposed regions of the prefusion closed trimer.

Mutations that prevent caspase cleavage of RIPK1 cause autoinflammatory disease.

RIPK1 is a key regulator of innate immune signalling pathways. To ensure an optimal inflammatory response, RIPK1 is regulated post-translationally by well-characterized ubiquitylation and phosphorylation events, as well as by caspase-8-mediated cleavage1-7. The physiological relevance of this cleavage event remains unclear, although it is thought to inhibit activation of RIPK3 and necroptosis8. Here we show that the heterozygous missense mutations D324N, D324H and D324Y prevent caspase cleavage of RIPK1 in humans and result in an early-onset periodic fever syndrome and severe intermittent lymphadenopathy-a condition we term 'cleavage-resistant RIPK1-induced autoinflammatory syndrome'. To define the mechanism for this disease, we generated a cleavage-resistant Ripk1D325A mutant mouse strain. Whereas Ripk1-/- mice died postnatally from systemic inflammation, Ripk1D325A/D325A mice died during embryogenesis. Embryonic lethality was completely prevented by the combined loss of Casp8 and Ripk3, but not by loss of Ripk3 or Mlkl alone. Loss of RIPK1 kinase activity also prevented Ripk1D325A/D325A embryonic lethality, although the mice died before weaning from multi-organ inflammation in a RIPK3-dependent manner. Consistently, Ripk1D325A/D325A and Ripk1D325A/+ cells were hypersensitive to RIPK3-dependent TNF-induced apoptosis and necroptosis. Heterozygous Ripk1D325A/+ mice were viable and grossly normal, but were hyper-responsive to inflammatory stimuli in vivo. Our results demonstrate the importance of caspase-mediated RIPK1 cleavage during embryonic development and show that caspase cleavage of RIPK1 not only inhibits necroptosis but also maintains inflammatory homeostasis throughout life.

Telomere-to-telomere assembly of a complete human X chromosome.

After two decades of improvements, the current human reference genome (GRCh38) is the most accurate and complete vertebrate genome ever produced. However, no single chromosome has been finished end to end, and hundreds of unresolved gaps persist1,2. Here we present a human genome assembly that surpasses the continuity of GRCh382, along with a gapless, telomere-to-telomere assembly of a human chromosome. This was enabled by high-coverage, ultra-long-read nanopore sequencing of the complete hydatidiform mole CHM13 genome, combined with complementary technologies for quality improvement and validation. Focusing our efforts on the human X chromosome3, we reconstructed the centromeric satellite DNA array (approximately 3.1 Mb) and closed the 29 remaining gaps in the current reference, including new sequences from the human pseudoautosomal regions and from cancer-testis ampliconic gene families (CT-X and GAGE). These sequences will be integrated into future human reference genome releases. In addition, the complete chromosome X, combined with the ultra-long nanopore data, allowed us to map methylation patterns across complex tandem repeats and satellite arrays. Our results demonstrate that finishing the entire human genome is now within reach, and the data presented here will facilitate ongoing efforts to complete the other human chromosomes.

A family of unusual immunoglobulin superfamily genes in an invertebrate histocompatibility complex.

Most colonial marine invertebrates are capable of allorecognition, the ability to distinguish between themselves and conspecifics. One long-standing question is whether invertebrate allorecognition genes are homologous to vertebrate histocompatibility genes. In the cnidarian Hydractinia symbiolongicarpus, allorecognition is controlled by at least two genes, Allorecognition 1 (Alr1) and Allorecognition 2 (Alr2), which encode highly polymorphic cell-surface proteins that serve as markers of self. Here, we show that Alr1 and Alr2 are part of a family of 41 Alr genes, all of which reside in a single genomic interval called the Allorecognition Complex (ARC). Using sensitive homology searches and highly accurate structural predictions, we demonstrate that the Alr proteins are members of the immunoglobulin superfamily (IgSF) with V-set and I-set Ig domains unlike any previously identified in animals. Specifically, their primary amino acid sequences lack many of the motifs considered diagnostic for V-set and I-set domains, yet they adopt secondary and tertiary structures nearly identical to canonical Ig domains. Thus, the V-set domain, which played a central role in the evolution of vertebrate adaptive immunity, was present in the last common ancestor of cnidarians and bilaterians. Unexpectedly, several Alr proteins also have immunoreceptor tyrosine-based activation motifs and immunoreceptor tyrosine-based inhibitory motifs in their cytoplasmic tails, suggesting they could participate in pathways homologous to those that regulate immunity in humans and flies. This work expands our definition of the IgSF with the addition of a family of unusual members, several of which play a role in invertebrate histocompatibility.

Somatic Mutations in UBA1 and Severe Adult-Onset Autoinflammatory Disease.

BACKGROUND: Adult-onset inflammatory syndromes often manifest with overlapping clinical features. Variants in ubiquitin-related genes, previously implicated in autoinflammatory disease, may define new disorders. METHODS: We analyzed peripheral-blood exome sequence data independent of clinical phenotype and inheritance pattern to identify deleterious mutations in ubiquitin-related genes. Sanger sequencing, immunoblotting, immunohistochemical testing, flow cytometry, and transcriptome and cytokine profiling were performed. CRISPR-Cas9-edited zebrafish were used as an in vivo model to assess gene function. RESULTS: We identified 25 men with somatic mutations affecting methionine-41 (p.Met41) in UBA1, the major E1 enzyme that initiates ubiquitylation. (The gene UBA1 lies on the X chromosome.) In such patients, an often fatal, treatment-refractory inflammatory syndrome develops in late adulthood, with fevers, cytopenias, characteristic vacuoles in myeloid and erythroid precursor cells, dysplastic bone marrow, neutrophilic cutaneous and pulmonary inflammation, chondritis, and vasculitis. Most of these 25 patients met clinical criteria for an inflammatory syndrome (relapsing polychondritis, Sweet's syndrome, polyarteritis nodosa, or giant-cell arteritis) or a hematologic condition (myelodysplastic syndrome or multiple myeloma) or both. Mutations were found in more than half the hematopoietic stem cells, including peripheral-blood myeloid cells but not lymphocytes or fibroblasts. Mutations affecting p.Met41 resulted in loss of the canonical cytoplasmic isoform of UBA1 and in expression of a novel, catalytically impaired isoform initiated at p.Met67. Mutant peripheral-blood cells showed decreased ubiquitylation and activated innate immune pathways. Knockout of the cytoplasmic UBA1 isoform homologue in zebrafish caused systemic inflammation. CONCLUSIONS: Using a genotype-driven approach, we identified a disorder that connects seemingly unrelated adult-onset inflammatory syndromes. We named this disorder the VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome. (Funded by the NIH Intramural Research Programs and the EU Horizon 2020 Research and Innovation Program.).

Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants.

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

HLA and autoantibodies define scleroderma subtypes and risk in African and European Americans and suggest a role for molecular mimicry.

Systemic sclerosis (SSc) is a clinically heterogeneous autoimmune disease characterized by mutually exclusive autoantibodies directed against distinct nuclear antigens. We examined HLA associations in SSc and its autoantibody subsets in a large, newly recruited African American (AA) cohort and among European Americans (EA). In the AA population, the African ancestry-predominant HLA-DRB1*08:04 and HLA-DRB1*11:02 alleles were associated with overall SSc risk, and the HLA-DRB1*08:04 allele was strongly associated with the severe antifibrillarin (AFA) antibody subset of SSc (odds ratio = 7.4). These African ancestry-predominant alleles may help explain the increased frequency and severity of SSc among the AA population. In the EA population, the HLA-DPB1*13:01 and HLA-DRB1*07:01 alleles were more strongly associated with antitopoisomerase (ATA) and anticentromere antibody-positive subsets of SSc, respectively, than with overall SSc risk, emphasizing the importance of HLA in defining autoantibody subtypes. The association of the HLA-DPB1*13:01 allele with the ATA+ subset of SSc in both AA and EA patients demonstrated a transancestry effect. A direct correlation between SSc prevalence and HLA-DPB1*13:01 allele frequency in multiple populations was observed (r = 0.98, P = 3 × 10-6). Conditional analysis in the autoantibody subsets of SSc revealed several associated amino acid residues, mostly in the peptide-binding groove of the class II HLA molecules. Using HLA α/ß allelic heterodimers, we bioinformatically predicted immunodominant peptides of topoisomerase 1, fibrillarin, and centromere protein A and discovered that they are homologous to viral protein sequences from the Mimiviridae and Phycodnaviridae families. Taken together, these data suggest a possible link between HLA alleles, autoantibodies, and environmental triggers in the pathogenesis of SSc.

Exome sequencing identifies genetic variants in anophthalmia and microphthalmia.

Anophthalmia and microphthalmia (A/M) are rare birth defects affecting up to 2 per 10,000 live births. These conditions are manifested by the absence of an eye or reduced eye volumes within the orbit leading to vision loss. Although clinical case series suggest a strong genetic component in A/M, few systematic investigations have been conducted on potential genetic contributions owing to low population prevalence. To overcome this challenge, we utilized DNA samples and data collected as part of the National Birth Defects Prevention Study (NBDPS). The NBDPS employed multi-center ascertainment of infants affected by A/M. We performed exome sequencing on 67 family trios and identified numerous genes affected by rare deleterious nonsense and missense variants in this cohort, including de novo variants. We identified 9 nonsense changes and 86 missense variants that are absent from the reference human population (Genome Aggregation Database), and we suggest that these are high priority candidate genes for A/M. We also performed literature curation, single cell transcriptome comparisons, and molecular pathway analysis on the candidate genes and performed protein structure modeling to determine the potential pathogenic variant consequences on PAX6 in this disease.

Multidonor analysis reveals structural elements, genetic determinants, and maturation pathway for HIV-1 neutralization by VRC01-class antibodies.

Antibodies of the VRC01 class neutralize HIV-1, arise in diverse HIV-1-infected donors, and are potential templates for an effective HIV-1 vaccine. However, the stochastic processes that generate repertoires in each individual of >10(12) antibodies make elicitation of specific antibodies uncertain. Here we determine the ontogeny of the VRC01 class by crystallography and next-generation sequencing. Despite antibody-sequence differences exceeding 50%, antibody-gp120 cocrystal structures reveal VRC01-class recognition to be remarkably similar. B cell transcripts indicate that VRC01-class antibodies require few specific genetic elements, suggesting that naive-B cells with VRC01-class features are generated regularly by recombination. Virtually all of these fail to mature, however, with only a few-likely one-ancestor B cell expanding to form a VRC01-class lineage in each donor. Developmental similarities in multiple donors thus reveal the generation of VRC01-class antibodies to be reproducible in principle, thereby providing a framework for attempts to elicit similar antibodies in the general population.

Evaluation of Recipients of Positive and Negative Secondary Findings Evaluations in a Hybrid CLIA-Research Sequencing Pilot.

While consensus regarding the return of secondary genomic findings in the clinical setting has been reached, debate about such findings in the research setting remains. We developed a hybrid, research-clinical translational genomics process for research exome data coupled with a CLIA-validated secondary findings analysis. Eleven intramural investigators from ten institutes at the National Institutes of Health piloted this process. Nearly 1,200 individuals were sequenced and 14 secondary findings were identified in 18 participants. Positive secondary findings were returned by a genetic counselor following a standardized protocol, including referrals for specialty follow-up care for the secondary finding local to the participants. Interviews were undertaken with 13 participants 4 months after receipt of a positive report. These participants reported minimal psychologic distress within a process to assimilate their results. Of the 13, 9 reported accessing the recommended health care services. A sample of 107 participants who received a negative findings report were surveyed 4 months after receiving it. They demonstrated good understanding of the negative secondary findings result and most expressed reassurance (64%) from that report. However, a notable minority (up to 17%) expressed confusion regarding the distinction of primary from secondary findings. This pilot shows it is feasible to couple CLIA-compliant secondary findings to research sequencing with minimal harms. Participants managed the surprise of a secondary finding with most following recommended follow up, yet some with negative findings conflated secondary and primary findings. Additional work is needed to understand barriers to follow-up care and help participants distinguish secondary from primary findings.

Exome sequencing of child-parent trios with bladder exstrophy: Findings in 26 children.

Bladder exstrophy (BE) is a rare, lower ventral midline defect with the bladder and part of the urethra exposed. The etiology of BE is unknown but thought to be influenced by genetic variation with more recent studies suggesting a role for rare variants. As such, we conducted paired-end exome sequencing in 26 child/mother/father trios. Three children had rare (allele frequency ≤ 0.0001 in several public databases) inherited variants in TSPAN4, one with a loss-of-function variant and two with missense variants. Two children had loss-of-function variants in TUBE1. Four children had rare missense or nonsense variants (one per child) in WNT3, CRKL, MYH9, or LZTR1, genes previously associated with BE. We detected 17 de novo missense variants in 13 children and three de novo loss-of-function variants (AKR1C2, PRRX1, PPM1D) in three children (one per child). We also detected rare compound heterozygous loss-of-function variants in PLCH2 and CLEC4M and rare inherited missense or loss-of-function variants in additional genes applying autosomal recessive (three genes) and X-linked recessive inheritance models (13 genes). Variants in two genes identified may implicate disruption in cell migration (TUBE1) and adhesion (TSPAN4) processes, mechanisms proposed for BE, and provide additional evidence for rare variants in the development of this defect.

Biallelic Mutations in DNAJC12 Cause Hyperphenylalaninemia, Dystonia, and Intellectual Disability.

Phenylketonuria (PKU, phenylalanine hydroxylase deficiency), an inborn error of metabolism, can be detected through newborn screening for hyperphenylalaninemia (HPA). Most individuals with HPA harbor mutations in the gene encoding phenylalanine hydroxylase (PAH), and a small proportion (2%) exhibit tetrahydrobiopterin (BH4) deficiency with additional neurotransmitter (dopamine and serotonin) deficiency. Here we report six individuals from four unrelated families with HPA who exhibited progressive neurodevelopmental delay, dystonia, and a unique profile of neurotransmitter deficiencies without mutations in PAH or BH4 metabolism disorder-related genes. In these six affected individuals, whole-exome sequencing (WES) identified biallelic mutations in DNAJC12, which encodes a heat shock co-chaperone family member that interacts with phenylalanine, tyrosine, and tryptophan hydroxylases catalyzing the BH4-activated conversion of phenylalanine into tyrosine, tyrosine into L-dopa (the precursor of dopamine), and tryptophan into 5-hydroxytryptophan (the precursor of serotonin), respectively. DNAJC12 was undetectable in fibroblasts from the individuals with null mutations. PAH enzyme activity was reduced in the presence of DNAJC12 mutations. Early treatment with BH4 and/or neurotransmitter precursors had dramatic beneficial effects and resulted in the prevention of neurodevelopmental delay in the one individual treated before symptom onset. Thus, DNAJC12 deficiency is a preventable and treatable cause of intellectual disability that should be considered in the early differential diagnosis when screening results are positive for HPA. Sequencing of DNAJC12 may resolve any uncertainty and should be considered in all children with unresolved HPA.

The complete genome sequence of a Neanderthal from the Altai Mountains.

We present a high-quality genome sequence of a Neanderthal woman from Siberia. We show that her parents were related at the level of half-siblings and that mating among close relatives was common among her recent ancestors. We also sequenced the genome of a Neanderthal from the Caucasus to low coverage. An analysis of the relationships and population history of available archaic genomes and 25 present-day human genomes shows that several gene flow events occurred among Neanderthals, Denisovans and early modern humans, possibly including gene flow into Denisovans from an unknown archaic group. Thus, interbreeding, albeit of low magnitude, occurred among many hominin groups in the Late Pleistocene. In addition, the high-quality Neanderthal genome allows us to establish a definitive list of substitutions that became fixed in modern humans after their separation from the ancestors of Neanderthals and Denisovans.