Archaeological, osteological and genetic evidence suggests that Neanderthals lived in small groups1,2; however, less is known about whether these groups were part of isolated communities or belonged to larger, well-connected populations3. The dense concentration of broadly contemporaneous Neanderthal sites in the Meuse Basin, Belgium4, provides a rare opportunity to study regional populations at high resolution. Here we generated genetic data from 27 Neanderthals who lived less than approximately 52,500 years ago from ten archaeological sites in Belgium and France, including a high-coverage genome from a 45,000-year-old individual from Goyet, Belgium. We show that most of these individuals are more closely related to one another than to other contemporaneous late Neanderthals in Europe. Further, some of these individuals carry DNA from a Neanderthal lineage predating the split of late Neanderthals. Although these Neanderthals overlapped temporally with early modern humans in northwestern Europe from around 47,000 years ago, we find no evidence of recent gene flow from modern humans. They also do not show the genetic signatures of mating among close relatives found in Altai Neanderthals, suggesting that they lived in larger or better-connected groups. Moreover, genetic load did not accumulate over time, arguing against progressive genetic deterioration as a driver of Neanderthal extinction.

Abstract The Teshik-Tash 1 child whose remains were found in Uzbekistan represents the southeastern-most extent of the known Neandertal range, providing an important link with the better studied Caucasus and Altai Mountain ranges. However, due to poor DNA preservation, studying the genetics of Teshik Tash 1 has remained elusive. Here we present analyses of the nuclear DNA from the Teshik-Tash 1, from extracts that are highly contaminated with present-day human DNA. To achieve this, we developed a new computational method, admixslug , that jointly models contamination and population relationships, in order to infer the relationship of a target individual from which only low-quality nuclear DNA is available, to high-quality archaic human genomes. After validating admixslug , we show that Teshik-Tash 1 is genetically more similar to later Neandertals from Western Eurasia than to older Neandertals from the Altai Mountains. We estimate that Teshik-Tash 1 split from the Western Eurasian lineage between 80,000 and 100,000 years ago. Despite the geographical proximity of Teshik-Tash 1 to the Denisovan range, we find no evidence for Denisovan ancestry in his genome. Our results demonstrate that admixslug enables the study of archaic human specimens in cases where DNA preservation was previously considered too poor for population genetic analyses.

Denisovans, an extinct sister group of Neandertals who lived in Eastern Eurasia during the Middle and Late Pleistocene, are known only from a handful of skeletal remains and limited genetic data, including the high-coverage genome of a woman who lived ∼65,000 years ago. Here, we present a second high-quality Denisovan genome, reconstructed from a molar found at Denisova Cave. It belonged to a man who lived ∼200,000 years ago in a small Denisovan group. This group mixed with early Neandertals and was then replaced by Denisovans who had mixed with later Neandertals. We show that in addition Denisovans received gene flow from hominins that diverged before the split of the ancestors of Denisovans and modern humans. The two Denisovan genomes allow us to disentangle Denisovan ancestry in present-day humans revealing contributions from at least three distinct Denisovan groups. In particular, Oceanians and South Asians independently inherited DNA from a deeply diverged Denisovan population which was likely isolated in South Asia. This supports an early migration of the ancestors of Oceanians through South Asia followed by the later arrival of the ancestors of present-day South Asians. East Asians do not share this Denisovan component in their genomes, suggesting that their ancestors arrived independently, perhaps by a northerly route. Finally, the two high-quality Denisovan genomes allow us to refine the catalogue of genetic changes that arose on the Denisovan lineage, some of which were contributed to present-day humans. ### Competing Interest Statement The authors have declared no competing interest. Max Planck Society, https://ror.org/01hhn8329 European Research Council, 694707 Russian Science Foundation, 24-18-00069

Denisovan ancestry in present-day humans is heterogeneously distributed and comes from genetically distinct Denisovan groups. Understanding the origin of this heterogeneity could provide insights into the early population history of modern humans in Eurasia. However, population movements and admixture after the initial dispersals of modern humans have obscured the origin of this heterogeneity. To address this, we investigated how Denisovan ancestry in early modern humans relates to that in present-day humans. We found that varying levels of Denisovan ancestry in Eurasians were shaped by admixture between diverse early modern human lineages. In particular, ancient Japanese individuals from the Jomon period have the least Denisovan ancestry among individuals from Eastern Eurasia, providing evidence for an ancient East Asian lineage with little to no Denisovan ancestry. By contrast, the earliest mainland East Asians harbor the most Denisovan ancestry in Eastern Eurasia, including ancestry from multiple divergent Denisovan groups, already before the last glacial maximum. Together with evidence of Denisovan ancestry sharing among ancient and present-day genomes, these patterns show how Denisovan ancestry in Eurasia spread through gene flow from early East Asians. Our study provides the first systematic investigation of Denisovan admixture across time and extends our understanding of human population history in Eurasia.

Disruption of gut barrier function and intestinal immune cell homeostasis are increasingly considered critical players in pathogenesis of extra-intestinal inflammatory diseases, including multiple sclerosis (MS) and its prototypical animal model, the experimental autoimmune encephalomyelitis (EAE). Breakdown of epithelial barriers increases intestinal permeability and systemic dissemination of microbiota-derived molecules. However, whether the gut-vascular barrier (GVB) is altered during EAE has not been reported. Here, we demonstrate that endothelial cell proliferation and vessel permeability increase before EAE clinical onset, leading to vascular remodeling and expansion of intestinal villi capillary bed during disease symptomatic phase in an antigen-independent manner. Concomitant to onset of angiogenesis observed prior to neurological symptoms, we identify an increase of intestinal perivascular immune cells characterized by the surface marker lymphatic vessel endothelial hyaluronic acid receptor 1 (LYVE-1). LYVE-1+ is expressed more frequently on B cells that show high levels of CD73 and have proangiogenic properties. B cell depletion was sufficient to mitigate enteric blood endothelial cell proliferation following immunization for EAE. In conclusion, we propose that altered intestinal vasculature driven by a specialized LYVE-1+ B cell subset promotes angiogenesis and that loss of GVB function is implicated in EAE development and autoimmunity.

Modern humans arrived in Europe more than 45,000 years ago, overlapping at least 5,000 years with Neanderthals1, 2, 3-4. Limited genomic data from these early modern humans have shown that at least two genetically distinct groups inhabited Europe, represented by Zlat & yacute; k & uring;& ncaron;, Czechia3 and Bacho Kiro, Bulgaria2. Here we deepen our understanding of early modern humans by analysing one high-coverage genome and five low-coverage genomes from approximately 45,000-year-old remains from Ilsenh & ouml;hle in Ranis, Germany4, and a further high-coverage genome from Zlat & yacute; k & uring;& ncaron;. We show that distant familial relationships link the Ranis and Zlat & yacute; k & uring;& ncaron; individuals and that they were part of the same small, isolated population that represents the deepest known split from the Out-of-Africa lineage. Ranis genomes harbour Neanderthal segments that originate from a single admixture event shared with all non-Africans that we date to approximately 45,000-49,000 years ago. This implies that ancestors of all non-Africans sequenced so far resided in a common population at this time, and further suggests that modern human remains older than 50,000 years from outside Africa represent different non-African populations.

Antibody-producing plasma cells fuel humoral immune responses. They also contribute to autoimmune diseases such as systemic lupus erythematosus or IgA nephropathy. Interleukin-6 and the tumor necrosis factor (TNF) family ligands BAFF (B cell-activating factor) and APRIL (a proliferation-inducing ligand) participate in plasma cell survival. BAFF binds to three receptors, BAFFR (BAFF receptor), TACI (transmembrane activator and CAML interactor), and BCMA (B cell maturation antigen), while APRIL binds to TACI, BCMA, and proteoglycans. However, which ligand–receptor pair(s) are required to maintain plasma cells in different body locations remains unknown. Here, by combining mouse genetic and pharmacological approaches, we found that plasma cells required BCMA and/or TACI but not BAFFR. BCMA responded exclusively to APRIL, while TACI responded to both BAFF and APRIL, identifying three self-sufficient ligand–receptor pairs for plasma cell maintenance: BAFF–TACI, APRIL–TACI, and APRIL–BCMA. Together, these actors accounted for 90% of circulating antibodies. In BAFF-ko mice, the reduction of plasma cells upon APRIL inhibition indicated that APRIL could function in the absence of BAFF–APRIL heteromers. No evidence was found that in the absence of BCMA and TACI, binding of APRIL to proteoglycans would help maintain plasma cells. IL-6, alone or together with BAFF and APRIL, supported mainly splenic plasmablasts and plasma cells and contributed to circulating IgG but not IgA levels. In conclusion, survival factors for plasma cells can vary with body location and with the antibody isotype that plasma cells produce. To efficiently target plasma cells, in particular IgA-producing ones, dual inhibition of BAFF and APRIL is required.

Gene flow from Neanderthals has shaped genetic and phenotypic variation in modern humans. We generated a catalog of Neanderthal ancestry segments in more than 300 genomes spanning the past 50,000 years. We examined how Neanderthal ancestry is shared among individuals over time. Our analysis revealed that the vast majority of Neanderthal gene flow is attributable to a single, shared extended period of gene flow that occurred between 50,500 to 43,500 years ago, as evidenced by ancestry correlation, colocalization of Neanderthal segments across individuals, and divergence from the sequenced Neanderthals. Most natural selection—positive and negative—on Neanderthal variants occurred rapidly after the gene flow. Our findings provide new insights into how contact with Neanderthals shaped modern human origins and adaptation.

Genetic kinship of ancient individuals can provide insights into their culture and social hierarchy, and is relevant for downstream genetic analyses. However, estimating relatedness from ancient DNA is difficult due to low-coverage, ascertainment bias, or contamination from various sources. Here, we present KIN, a method to estimate the relatedness of a pair of individuals from the identical-by-descent segments they share. KIN accurately classifies up to 3rd-degree relatives using at least 0.05x sequence coverage and differentiates siblings from parent-child pairs. It incorporates additional models to adjust for contamination and detect inbreeding, which improves classification accuracy.

The vasculature is a key regulator of leukocyte trafficking into the central nervous system (CNS) during inflammatory diseases including multiple sclerosis (MS). However, the impact of endothelial‐derived factors on CNS immune responses remains unknown. Bioactive lipids, in particular oxysterols downstream of Cholesterol‐25‐hydroxylase (Ch25h), promote neuroinflammation but their functions in the CNS are not well‐understood. Using floxed‐reporter Ch25h knock‐in mice, we trace Ch25h expression to CNS endothelial cells (ECs) and myeloid cells and demonstrate that Ch25h ablation specifically from ECs attenuates experimental autoimmune encephalomyelitis (EAE). Mechanistically, inflamed Ch25h‐deficient CNS ECs display altered lipid metabolism favoring polymorphonuclear myeloid‐derived suppressor cell (PMN‐MDSC) expansion, which suppresses encephalitogenic T lymphocyte proliferation. Additionally, endothelial Ch25h‐deficiency combined with immature neutrophil mobilization into the blood circulation nearly completely protects mice from EAE. Our findings reveal a central role for CNS endothelial Ch25h in promoting neuroinflammation by inhibiting the expansion of immunosuppressive myeloid cell populations.

BACKGROUND:Multiple sclerosis (MS) is a chronic disabling disease of the central nervous system (CNS) commonly affecting young adults. There is increasing evidence that environmental factors are important in the development and course of MS. The metabolic syndrome (MetS) which comprises dyslipidemia has been associated with a worse outcome in MS disease. Furthermore, the lipid-lowering drug class of statins has been proposed to improve MS disease course. However, cholesterol is also rate-limiting for myelin biogenesis and promotes remyelination in MS animal models. Thus, the impact of circulating blood cholesterol levels during the disease remains debated and controversial.METHODS:We assessed the role of circulating cholesterol on the murine model of MS, the experimental autoimmune encephalomyelitis (EAE) disease using two different approaches: (1) the mouse model of familial hypercholesterolemia induced by low-density lipoprotein receptor (LDLr) deficiency, and (2) the use of the monoclonal anti-PCSK9 neutralizing antibody alirocumab, which reduces LDLr degradation and consequently lowers blood levels of cholesterol.RESULTS:Elevated blood cholesterol levels induced by LDLr deficiency did not worsen clinical symptoms of mice during EAE. In addition, we observed that the anti-PCSK9 antibody alirocumab did not influence EAE disease course, nor modulate the immune response in EAE.CONCLUSIONS:These findings suggest that blood cholesterol level has no direct role in neuro-inflammatory diseases and that the previously shown protective effects of statins in MS are not related to circulating cholesterol.

Microbiota-derived metabolites are important molecules connecting the gut to the brain. Over the last decade, several studies have highlighted the importance of gut-derived metabolites in the development of multiple sclerosis (MS). Indeed, microbiota-derived metabolites modulate the immune system and affect demyelination. Here, we discuss the current knowledge about microbiota-derived metabolites implications in MS and in different mouse models of neuroinflammation. We focus on the main families of microbial metabolites that play a role during neuroinflammation. A better understanding of the role of those metabolites may lead to new therapeutical avenues to treat neuroinflammatory diseases targeting the gut–brain axis.

Principal component analysis (PCA) and F-statistics sensu Patterson are two of the most widely used population genetic tools to study human genetic variation. Here, I derive explicit connections between the two approaches and show that these two methods are closely related. F-statistics have a simple geometrical interpretation in the context of PCA, and orthogonal projections are a key concept to establish this link. I show that for any pair of populations, any population that is admixed as determined by an F3-statistic will lie inside a circle on a PCA plot. Furthermore, the F4-statistic is closely related to an angle measurement, and will be zero if the differences between pairs of populations intersect at a right angle in PCA space. I illustrate my results on two examples, one of Western Eurasian, and one of global human diversity. In both examples, I find that the first few PCs are sufficient to approximate most F-statistics, and that PCA plots are effective at predicting F-statistics. Thus, while F-statistics are commonly understood in terms of discrete populations, the geometric perspective illustrates that they can be viewed in a framework of populations that vary in a more continuous manner. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Genomic analyses of Neanderthals have previously provided insights into their population history and relationship to modern humans 1 – 8 , but the social organization of Neanderthal communities remains poorly understood. Here we present genetic data for 13 Neanderthals from two Middle Palaeolithic sites in the Altai Mountains of southern Siberia: 11 from Chagyrskaya Cave 9 , 10 and 2 from Okladnikov Cave 11 —making this one of the largest genetic studies of a Neanderthal population to date. We used hybridization capture to obtain genome-wide nuclear data, as well as mitochondrial and Y-chromosome sequences. Some Chagyrskaya individuals were closely related, including a father–daughter pair and a pair of second-degree relatives, indicating that at least some of the individuals lived at the same time. Up to one-third of these individuals’ genomes had long segments of homozygosity, suggesting that the Chagyrskaya Neanderthals were part of a small community. In addition, the Y-chromosome diversity is an order of magnitude lower than the mitochondrial diversity, a pattern that we found is best explained by female migration between communities. Thus, the genetic data presented here provide a detailed documentation of the social organization of an isolated Neanderthal community at the easternmost extent of their known range.

The Massim, a cultural region encompassing the southeastern tip of mainland Papua New Guinea (PNG) and nearby islands, is famous for the Kula trading network, in which different valuables circulate in different directions among the islands. To explore Massim genetic history, we generated genome-wide data from across the region and found variable levels of Papuan-related (indigenous) ancestry, including a novel ancestry associated with Rossel Island, and Austronesian-related ancestry that arrived later. We find genetic evidence for different patterns of contact across PNG, the Massim, and the Bismarck and Solomon Archipelagoes for Austronesian-related vs. Papuan-related ancestry, with more geographic restriction for the latter. Moreover, Kula-practicing groups share more genetic similarity than do other groups, and this likely predates the origin of Kula, suggesting that high between-group contact facilitated the formation of Kula. Our study provides the first comprehensive genome-wide assessment of Massim inhabitants and new insights into the fascinating Kula system.

Proteins associated with the spindle apparatus, a cytoskeletal structure that ensures the proper segregation of chromosomes during cell division, experienced an unusual number of amino acid substitutions in modern humans after the split from the ancestors of Neandertals and Denisovans. Here, we analyze the history of these substitutions and show that some of the genes in which they occur may have been targets of positive selection. We also find that the two changes in the kinetochore scaffold 1 (KNL1) protein, previously believed to be specific to modern humans, were present in some Neandertals. We show that the KNL1 gene of these Neandertals shared a common ancestor with present-day Africans about 200,000 years ago due to gene flow from the ancestors (or relatives) of modern humans into Neandertals. Subsequently, some non-Africans inherited this modern human-like gene variant from Neandertals, but none inherited the ancestral gene variants. These results add to the growing evidence of early contacts between modern humans and archaic groups in Eurasia and illustrate the intricate relationships among these groups.

Neandertal DNA makes up 2-3% of the genomes of all non-African individuals. The patterns of Neandertal ancestry in modern humans have been used to estimate that this is the result of gene flow that occurred during the expansion of modern humans into Eurasia, but the precise dates of this event remain largely unknown. Here, we introduce an extended admixture pulse model that allows joint estimation of the timing and duration of gene flow. This model leads to simple expressions for both the admixture segment distribution and the decay curve of ancestry linkage disequilibrium, and we show that these two statistics are closely related. In simulations, we find that estimates of the mean time of admixture are largely robust to details in gene flow models, but that the duration of the gene flow can only be recovered if gene flow is very recent and the exact recombination map is known. These results imply that gene flow from Neandertals into modern humans could have happened over hundreds of generations. Ancient genomes from the time around the admixture event are thus likely required to resolve the question when, where, and for how long humans and Neandertals interacted.