Genetic insights into the social organization of Neanderthals.

Genomic analyses of Neanderthals have previously provided insights into their population history and relationship to modern humans1-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 Cave9,10 and 2 from Okladnikov Cave11-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.

Genome of a middle Holocene hunter-gatherer from Wallacea.

Much remains unknown about the population history of early modern humans in southeast Asia, where the archaeological record is sparse and the tropical climate is inimical to the preservation of ancient human DNA1. So far, only two low-coverage pre-Neolithic human genomes have been sequenced from this region. Both are from mainland Hòabìnhian hunter-gatherer sites: Pha Faen in Laos, dated to 7939-7751 calibrated years before present (yr cal BP; present taken as AD 1950), and Gua Cha in Malaysia (4.4-4.2 kyr cal BP)1. Here we report, to our knowledge, the first ancient human genome from Wallacea, the oceanic island zone between the Sunda Shelf (comprising mainland southeast Asia and the continental islands of western Indonesia) and Pleistocene Sahul (Australia-New Guinea). We extracted DNA from the petrous bone of a young female hunter-gatherer buried 7.3-7.2 kyr cal BP at the limestone cave of Leang Panninge2 in South Sulawesi, Indonesia. Genetic analyses show that this pre-Neolithic forager, who is associated with the 'Toalean' technocomplex3,4, shares most genetic drift and morphological similarities with present-day Papuan and Indigenous Australian groups, yet represents a previously unknown divergent human lineage that branched off around the time of the split between these populations approximately 37,000 years ago5. We also describe Denisovan and deep Asian-related ancestries in the Leang Panninge genome, and infer their large-scale displacement from the region today.

Initial Upper Palaeolithic humans in Europe had recent Neanderthal ancestry.

Modern humans appeared in Europe by at least 45,000 years ago1-5, but the extent of their interactions with Neanderthals, who disappeared by about 40,000 years ago6, and their relationship to the broader expansion of modern humans outside Africa are poorly understood. Here we present genome-wide data from three individuals dated to between 45,930 and 42,580 years ago from Bacho Kiro Cave, Bulgaria1,2. They are the earliest Late Pleistocene modern humans known to have been recovered in Europe so far, and were found in association with an Initial Upper Palaeolithic artefact assemblage. Unlike two previously studied individuals of similar ages from Romania7 and Siberia8 who did not contribute detectably to later populations, these individuals are more closely related to present-day and ancient populations in East Asia and the Americas than to later west Eurasian populations. This indicates that they belonged to a modern human migration into Europe that was not previously known from the genetic record, and provides evidence that there was at least some continuity between the earliest modern humans in Europe and later people in Eurasia. Moreover, we find that all three individuals had Neanderthal ancestors a few generations back in their family history, confirming that the first European modern humans mixed with Neanderthals and suggesting that such mixing could have been common.

Unique and redundant roles of mouse BCMA, TACI, BAFF, APRIL, and IL-6 in supporting antibody-producing cells in different tissues.

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.

Endothelial cell-derived oxysterol ablation attenuates experimental autoimmune encephalomyelitis.

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.

Perivascular B cells link intestinal angiogenesis to immunity and to the gut-brain axis during neuroinflammation.

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.

Assessing Human Genome-wide Variation in the Massim Region of Papua New Guinea and Implications for the Kula Trading Tradition.

The Massim, a cultural region that includes the southeastern tip of mainland Papua New Guinea (PNG) and nearby PNG offshore islands, is renowned for a trading network called Kula, in which different valuable items circulate in different directions among some of the islands. Although the Massim has been a focus of anthropological investigation since the pioneering work of Malinowski in 1922, the genetic background of its inhabitants remains relatively unexplored. To characterize the Massim genomically, we generated genome-wide SNP data from 192 individuals from 15 groups spanning the entire region. Analyzing these together with comparative data, we found that all Massim individuals have variable Papuan-related (indigenous) and Austronesian-related (arriving ∼3,000 years ago) ancestries. Individuals from Rossel Island in southern Massim, speaking an isolate Papuan language, have the highest amount of a distinct Papuan ancestry. We also investigated the recent contact via sharing of identical by descent (IBD) genomic segments and found that Austronesian-related IBD tracts are widely distributed geographically, but Papuan-related tracts are shared exclusively between the PNG mainland and Massim, and between the Bismarck and Solomon Archipelagoes. Moreover, the Kula-practicing groups of the Massim show higher IBD sharing among themselves than do groups that do not participate in Kula. This higher sharing predates the formation of Kula, suggesting that extensive contact between these groups since the Austronesian settlement may have 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.

A high-coverage Neandertal genome from Chagyrskaya Cave.

We sequenced the genome of a Neandertal from Chagyrskaya Cave in the Altai Mountains, Russia, to 27-fold genomic coverage. We show that this Neandertal was a female and that she was more related to Neandertals in western Eurasia [Prüfer et al., Science 358, 655-658 (2017); Hajdinjak et al., Nature 555, 652-656 (2018)] than to Neandertals who lived earlier in Denisova Cave [Prüfer et al., Nature 505, 43-49 (2014)], which is located about 100 km away. About 12.9% of the Chagyrskaya genome is spanned by homozygous regions that are between 2.5 and 10 centiMorgans (cM) long. This is consistent with the fact that Siberian Neandertals lived in relatively isolated populations of less than 60 individuals. In contrast, a Neandertal from Europe, a Denisovan from the Altai Mountains, and ancient modern humans seem to have lived in populations of larger sizes. The availability of three Neandertal genomes of high quality allows a view of genetic features that were unique to Neandertals and that are likely to have been at high frequency among them. We find that genes highly expressed in the striatum in the basal ganglia of the brain carry more amino-acid-changing substitutions than genes expressed elsewhere in the brain, suggesting that the striatum may have evolved unique functions in Neandertals.

An Extended Admixture Pulse Model Reveals the Limitations to Human-Neandertal Introgression Dating.

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.

Lowering blood cholesterol does not affect neuroinflammation in experimental autoimmune encephalomyelitis.

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.

How Microbiota-Derived Metabolites Link the Gut to the Brain during Neuroinflammation.

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.

Genetic Landscapes Reveal How Human Genetic Diversity Aligns with Geography.

Geographic patterns in human genetic diversity carry footprints of population history and provide insights for genetic medicine and its application across human populations. Summarizing and visually representing these patterns of diversity has been a persistent goal for human geneticists, and has revealed that genetic differentiation is frequently correlated with geographic distance. However, most analytical methods to represent population structure do not incorporate geography directly, and it must be considered post hoc alongside a visual summary of the genetic structure. Here, we estimate "effective migration" surfaces to visualize how human genetic diversity is geographically structured. The results reveal local patterns of differentiation in detail and emphasize that while genetic similarity generally decays with geographic distance, the relationship is often subtly distorted. Overall, the visualizations provide a new perspective on genetics and geography in humans and insight to the geographic distribution of human genetic variation.

Mitonuclear interactions and introgression genomics of macaque monkeys (Macaca) highlight the influence of behaviour on genome evolution.

In most macaques, females are philopatric and males migrate from their natal ranges, which results in pronounced divergence of mitochondrial genomes within and among species. We therefore predicted that some nuclear genes would have to acquire compensatory mutations to preserve compatibility with diverged interaction partners from the mitochondria. We additionally expected that these sex-differences would have distinctive effects on gene flow in the X and autosomes. Using new genomic data from 29 individuals from eight species of Southeast Asian macaque, we identified evidence of natural selection associated with mitonuclear interactions, including extreme outliers of interspecies differentiation and metrics of positive selection, low intraspecies polymorphism and atypically long runs of homozygosity associated with nuclear-encoded genes that interact with mitochondria-encoded genes. In one individual with introgressed mitochondria, we detected a small but significant enrichment of autosomal introgression blocks from the source species of her mitochondria that contained genes which interact with mitochondria-encoded loci. Our analyses also demonstrate that sex-specific demography sculpts genetic exchange across multiple species boundaries. These findings show that behaviour can have profound but indirect effects on genome evolution by influencing how interacting components of different genomic compartments (mitochondria, the autosomes and the sex chromosomes) move through time and space.

Human local adaptation of the TRPM8 cold receptor along a latitudinal cline.

Ambient temperature is a critical environmental factor for all living organisms. It was likely an important selective force as modern humans recently colonized temperate and cold Eurasian environments. Nevertheless, as of yet we have limited evidence of local adaptation to ambient temperature in populations from those environments. To shed light on this question, we exploit the fact that humans are a cosmopolitan species that inhabit territories under a wide range of temperatures. Focusing on cold perception-which is central to thermoregulation and survival in cold environments-we show evidence of recent local adaptation on TRPM8. This gene encodes for a cation channel that is, to date, the only temperature receptor known to mediate an endogenous response to moderate cold. The upstream variant rs10166942 shows extreme population differentiation, with frequencies that range from 5% in Nigeria to 88% in Finland (placing this SNP in the 0.02% tail of the FST empirical distribution). When all populations are jointly analyzed, allele frequencies correlate with latitude and temperature beyond what can be explained by shared ancestry and population substructure. Using a Bayesian approach, we infer that the allele originated and evolved neutrally in Africa, while positive selection raised its frequency to different degrees in Eurasian populations, resulting in allele frequencies that follow a latitudinal cline. We infer strong positive selection, in agreement with ancient DNA showing high frequency of the allele in Europe 3,000 to 8,000 years ago. rs10166942 is important phenotypically because its ancestral allele is protective of migraine. This debilitating disorder varies in prevalence across human populations, with highest prevalence in individuals of European descent-precisely the population with the highest frequency of rs10166942 derived allele. We thus hypothesize that local adaptation on previously neutral standing variation may have contributed to the genetic differences that exist in the prevalence of migraine among human populations today.

Altitude adaptation in Tibetans caused by introgression of Denisovan-like DNA.

As modern humans migrated out of Africa, they encountered many new environmental conditions, including greater temperature extremes, different pathogens and higher altitudes. These diverse environments are likely to have acted as agents of natural selection and to have led to local adaptations. One of the most celebrated examples in humans is the adaptation of Tibetans to the hypoxic environment of the high-altitude Tibetan plateau. A hypoxia pathway gene, EPAS1, was previously identified as having the most extreme signature of positive selection in Tibetans, and was shown to be associated with differences in haemoglobin concentration at high altitude. Re-sequencing the region around EPAS1 in 40 Tibetan and 40 Han individuals, we find that this gene has a highly unusual haplotype structure that can only be convincingly explained by introgression of DNA from Denisovan or Denisovan-related individuals into humans. Scanning a larger set of worldwide populations, we find that the selected haplotype is only found in Denisovans and in Tibetans, and at very low frequency among Han Chinese. Furthermore, the length of the haplotype, and the fact that it is not found in any other populations, makes it unlikely that the haplotype sharing between Tibetans and Denisovans was caused by incomplete ancestral lineage sorting rather than introgression. Our findings illustrate that admixture with other hominin species has provided genetic variation that helped humans to adapt to new environments.

Estimating the Ages of Selection Signals from Different Epochs in Human History.

Genetic variation harbors signatures of natural selection driven by selective pressures that are often unknown. Estimating the ages of selection signals may allow reconstructing the history of environmental changes that shaped human phenotypes and diseases. We have developed an approximate Bayesian computation (ABC) approach to estimate allele ages under a model of selection on new mutations and under demographic models appropriate for human populations. We have applied it to two resequencing data sets: An ultra-high depth data set from a relatively small sample of unrelated individuals and a lower depth data set in a larger sample with transmission information. In addition to evaluating the accuracy of our method based on simulations, for each SNP, we assessed the consistency between the posterior probabilities estimated by the ABC approach and the ancient DNA record, finding good agreement between the two types of data and methods. Applying this ABC approach to data for eight single nucleotide polymorphisms (SNPs), we were able to rule out an onset of selection prior to the dispersal out-of-Africa for three of them and more recent than the spread of agriculture for an additional three SNPs.

A Selective Sweep on a Deleterious Mutation in CPT1A in Arctic Populations.

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6-23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.

Selection on a variant associated with improved viral clearance drives local, adaptive pseudogenization of interferon lambda 4 (IFNL4).

Interferon lambda 4 gene (IFNL4) encodes IFN-λ4, a new member of the IFN-λ family with antiviral activity. In humans IFNL4 open reading frame is truncated by a polymorphic frame-shift insertion that eliminates IFN-λ4 and turns IFNL4 into a polymorphic pseudogene. Functional IFN-λ4 has antiviral activity but the elimination of IFN-λ4 through pseudogenization is strongly associated with improved clearance of hepatitis C virus (HCV) infection. We show that functional IFN-λ4 is conserved and evolutionarily constrained in mammals and thus functionally relevant. However, the pseudogene has reached moderately high frequency in Africa, America, and Europe, and near fixation in East Asia. In fact, the pseudogenizing variant is among the 0.8% most differentiated SNPs between Africa and East Asia genome-wide. Its raise in frequency is associated with additional evidence of positive selection, which is strongest in East Asia, where this variant falls in the 0.5% tail of SNPs with strongest signatures of recent positive selection genome-wide. Using a new Approximate Bayesian Computation (ABC) approach we infer that the pseudogenizing allele appeared just before the out-of-Africa migration and was immediately targeted by moderate positive selection; selection subsequently strengthened in European and Asian populations resulting in the high frequency observed today. This provides evidence for a changing adaptive process that, by favoring IFN-λ4 inactivation, has shaped present-day phenotypic diversity and susceptibility to disease.

Phylogenomics at the tips: inferring lineages and their demographic history in a tropical lizard, Carlia amax.

High-throughput sequencing approaches offer opportunities to better understand the evolutionary processes driving diversification, particularly in nonmodel organisms. In particular, the 100-1000's of loci that can now be sequenced are providing unprecedented power in population, speciation and phylogenetic studies. Here, we apply an exon capture approach to generate >99% complete sequence and SNP data across >2000 loci from a tropical skink, Carlia amax, and exploit these data to identify divergent lineages and infer their relationships and demographic histories. This is especially relevant to low-dispersal tropical taxa that often have cryptic diversity and spatially dynamic histories. For C. amax, clustering of nuclear SNPs and coalescent-based species delimitation analyses identify four divergent lineages, one fewer than predicted based on geographically coherent mtDNA clades (>9.4% sequence divergence). Three of these lineages are widespread and parapatric on the mainland, whereas the most divergent is restricted to islands off the northeast Northern Territory. Tests for population expansion reject an equilibrium isolation-by-distance model for two of the three widespread lineages and infer refugial expansion sources in the relatively mesic northeast Top End and northwest Kimberley. The latter is already recognized as a hotspot of endemism, but our results also suggest that a stronger focus on the northeast Top End, and adjacent islands is warranted. More generally, our results show how genome-reduction methods such as exon capture can yield insights into the pattern and dynamics of biodiversity across complex landscapes with as yet poorly understood biogeographic history and how exon data can link between population and phylogenetic questions.

Distinguishing between selective sweeps from standing variation and from a de novo mutation.

An outstanding question in human genetics has been the degree to which adaptation occurs from standing genetic variation or from de novo mutations. Here, we combine several common statistics used to detect selection in an Approximate Bayesian Computation (ABC) framework, with the goal of discriminating between models of selection and providing estimates of the age of selected alleles and the selection coefficients acting on them. We use simulations to assess the power and accuracy of our method and apply it to seven of the strongest sweeps currently known in humans. We identify two genes, ASPM and PSCA, that are most likely affected by selection on standing variation; and we find three genes, ADH1B, LCT, and EDAR, in which the adaptive alleles seem to have swept from a new mutation. We also confirm evidence of selection for one further gene, TRPV6. In one gene, G6PD, neither neutral models nor models of selective sweeps fit the data, presumably because this locus has been subject to balancing selection.