Genetic Adaptation and Neandertal Admixture Shaped the Immune System of Human Populations.

Humans differ in the outcome that follows exposure to life-threatening pathogens, yet the extent of population differences in immune responses and their genetic and evolutionary determinants remain undefined. Here, we characterized, using RNA sequencing, the transcriptional response of primary monocytes from Africans and Europeans to bacterial and viral stimuli-ligands activating Toll-like receptor pathways (TLR1/2, TLR4, and TLR7/8) and influenza virus-and mapped expression quantitative trait loci (eQTLs). We identify numerous cis-eQTLs that contribute to the marked differences in immune responses detected within and between populations and a strong trans-eQTL hotspot at TLR1 that decreases expression of pro-inflammatory genes in Europeans only. We find that immune-responsive regulatory variants are enriched in population-specific signals of natural selection and show that admixture with Neandertals introduced regulatory variants into European genomes, affecting preferentially responses to viral challenges. Together, our study uncovers evolutionarily important determinants of differences in host immune responsiveness between human populations.

Denisovan admixture facilitated environmental adaptation in Papua New Guinean populations.

Neandertals and Denisovans, having inhabited distinct regions in Eurasia and possibly Oceania for over 200,000 y, experienced ample time to adapt to diverse environmental challenges these regions presented. Among present-day human populations, Papua New Guineans (PNG) stand out as one of the few carrying substantial amounts of both Neandertal and Denisovan DNA, a result of past admixture events with these archaic human groups. This study investigates the distribution of introgressed Denisovan and Neandertal DNA within two distinct PNG populations, residing in the highlands of Mt Wilhelm and the lowlands of Daru Island. These locations exhibit unique environmental features, some of which may parallel the challenges that archaic humans once confronted and adapted to. Our results show that PNG highlanders carry higher levels of Denisovan DNA compared to PNG lowlanders. Among the Denisovan-like haplotypes with higher frequencies in highlander populations, those exhibiting the greatest frequency difference compared to lowlander populations also demonstrate more pronounced differences in population frequencies than frequency-matched nonarchaic variants. Two of the five most highly differentiated of those haplotypes reside in genomic areas linked to brain development genes. Conversely, Denisovan-like haplotypes more frequent in lowlanders overlap with genes associated with immune response processes. Our findings suggest that Denisovan DNA has provided genetic variation associated with brain biology and immune response to PNG genomes, some of which might have facilitated adaptive processes to environmental challenges.

Detection of Neanderthal Adaptively Introgressed Genetic Variants That Modulate Reporter Gene Expression in Human Immune Cells.

Although some variation introgressed from Neanderthals has undergone selective sweeps, little is known about its functional significance. We used a Massively Parallel Reporter Assay (MPRA) to assay 5,353 high-frequency introgressed variants for their ability to modulate the gene expression within 170 bp of endogenous sequence. We identified 2,548 variants in active putative cis-regulatory elements (CREs) and 292 expression-modulating variants (emVars). These emVars are predicted to alter the binding motifs of important immune transcription factors, are enriched for associations with neutrophil and white blood cell count, and are associated with the expression of genes that function in innate immune pathways including inflammatory response and antiviral defense. We combined the MPRA data with other data sets to identify strong candidates to be driver variants of positive selection including an emVar that may contribute to protection against severe COVID-19 response. We endogenously deleted two CREs containing expression-modulation variants linked to immune function, rs11624425 and rs80317430, identifying their primary genic targets as ELMSAN1, and PAN2 and STAT2, respectively, three genes differentially expressed during influenza infection. Overall, we present the first database of experimentally identified expression-modulating Neanderthal-introgressed alleles contributing to potential immune response in modern humans.

A novel compound heterozygous leptin receptor mutation causes more severe obesity than in Leprdb/db mice.

The leptin receptor (Lepr) pathway is important for food intake regulation, energy expenditure, and body weight. Mutations in leptin and the Lepr have been shown to cause early-onset severe obesity in mice and humans. In studies with C57BL/6NCrl mice, we found a mouse with extreme obesity. To identify a putative spontaneous new form of monogenic obesity, we performed backcross studies with this mouse followed by a quantitative trait locus (QTL) analysis and sequencing of the selected chromosomal QTL region. We thereby identified a novel Lepr mutation (C57BL/6N-LeprL536Hfs*6-1NKB), which is located at chromosome 4, exon 11 within the CRH2-leptin-binding site. Compared with C57BL/6N mice, LeprL536Hfs*6 develop early onset obesity and their body weight exceeds that of Leprdb/db mice at an age of 30 weeks. Similar to Leprdb/db mice, the LeprL536Hfs*6 model is characterized by hyperphagia, obesity, lower energy expenditure and activity, hyperglycemia, and hyperinsulinemia compared with C57BL/6N mice. Crossing Leprdb/wt with LeprL536Hfs*6/wt mice results in compound heterozygous LeprL536Hfs*6/db mice, which develop even higher body weight and fat mass than both homozygous Leprdb/db and LeprL536Hfs*6 mice. Compound heterozygous Lepr deficiency affecting functionally different regions of the Lepr causes more severe obesity than the parental homozygous mutations.

The Contribution of Neanderthals to Phenotypic Variation in Modern Humans.

Assessing the genetic contribution of Neanderthals to non-disease phenotypes in modern humans has been difficult because of the absence of large cohorts for which common phenotype information is available. Using baseline phenotypes collected for 112,000 individuals by the UK Biobank, we can now elaborate on previous findings that identified associations between signatures of positive selection on Neanderthal DNA and various modern human traits but not any specific phenotypic consequences. Here, we show that Neanderthal DNA affects skin tone and hair color, height, sleeping patterns, mood, and smoking status in present-day Europeans. Interestingly, multiple Neanderthal alleles at different loci contribute to skin and hair color in present-day Europeans, and these Neanderthal alleles contribute to both lighter and darker skin tones and hair color, suggesting that Neanderthals themselves were most likely variable in these traits.

Detecting ancient positive selection in humans using extended lineage sorting.

Natural selection that affected modern humans early in their evolution has likely shaped some of the traits that set present-day humans apart from their closest extinct and living relatives. The ability to detect ancient natural selection in the human genome could provide insights into the molecular basis for these human-specific traits. Here, we introduce a method for detecting ancient selective sweeps by scanning for extended genomic regions where our closest extinct relatives, Neandertals and Denisovans, fall outside of the present-day human variation. Regions that are unusually long indicate the presence of lineages that reached fixation in the human population faster than expected under neutral evolution. Using simulations, we show that the method is able to detect ancient events of positive selection and that it can differentiate those from background selection. Applying our method to the 1000 Genomes data set, we find evidence for ancient selective sweeps favoring regulatory changes and present a list of genomic regions that are predicted to underlie positively selected human specific traits.

The genomic landscape of Neanderthal ancestry in present-day humans.

Genomic studies have shown that Neanderthals interbred with modern humans, and that non-Africans today are the products of this mixture. The antiquity of Neanderthal gene flow into modern humans means that genomic regions that derive from Neanderthals in any one human today are usually less than a hundred kilobases in size. However, Neanderthal haplotypes are also distinctive enough that several studies have been able to detect Neanderthal ancestry at specific loci. We systematically infer Neanderthal haplotypes in the genomes of 1,004 present-day humans. Regions that harbour a high frequency of Neanderthal alleles are enriched for genes affecting keratin filaments, suggesting that Neanderthal alleles may have helped modern humans to adapt to non-African environments. We identify multiple Neanderthal-derived alleles that confer risk for disease, suggesting that Neanderthal alleles continue to shape human biology. An unexpected finding is that regions with reduced Neanderthal ancestry are enriched in genes, implying selection to remove genetic material derived from Neanderthals. Genes that are more highly expressed in testes than in any other tissue are especially reduced in Neanderthal ancestry, and there is an approximately fivefold reduction of Neanderthal ancestry on the X chromosome, which is known from studies of diverse species to be especially dense in male hybrid sterility genes. These results suggest that part of the explanation for genomic regions of reduced Neanderthal ancestry is Neanderthal alleles that caused decreased fertility in males when moved to a modern human genetic background.

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.

Introgression of Neandertal- and Denisovan-like Haplotypes Contributes to Adaptive Variation in Human Toll-like Receptors.

Pathogens and the diseases they cause have been among the most important selective forces experienced by humans during their evolutionary history. Although adaptive alleles generally arise by mutation, introgression can also be a valuable source of beneficial alleles. Archaic humans, who lived in Europe and Western Asia for more than 200,000 years, were probably well adapted to this environment and its local pathogens. It is therefore conceivable that modern humans entering Europe and Western Asia who admixed with them obtained a substantial immune advantage from the introgression of archaic alleles. Here we document a cluster of three Toll-like receptors (TLR6-TLR1-TLR10) in modern humans that carries three distinct archaic haplotypes, indicating repeated introgression from archaic humans. Two of these haplotypes are most similar to the Neandertal genome, and the third haplotype is most similar to the Denisovan genome. The Toll-like receptors are key components of innate immunity and provide an important first line of immune defense against bacteria, fungi, and parasites. The unusually high allele frequencies and unexpected levels of population differentiation indicate that there has been local positive selection on multiple haplotypes at this locus. We show that the introgressed alleles have clear functional effects in modern humans; archaic-like alleles underlie differences in the expression of the TLR genes and are associated with increased [corrected] microbial resistance and increased allergic disease in large cohorts. This provides strong evidence for recurrent adaptive introgression at the TLR6-TLR1-TLR10 locus, resulting in differences in disease phenotypes in modern humans.

Erratum: Introgression of Neandertal- and Denisovan-like Haplotypes Contributes to Adaptive Variation in Human Toll-like Receptors.

[This corrects the article DOI: 10.1016/j.ajhg.2015.11.015.].

Palaeoproteomic evidence identifies archaic hominins associated with the Châtelperronian at the Grotte du Renne.

In Western Europe, the Middle to Upper Paleolithic transition is associated with the disappearance of Neandertals and the spread of anatomically modern humans (AMHs). Current chronological, behavioral, and biological models of this transitional period hinge on the Châtelperronian technocomplex. At the site of the Grotte du Renne, Arcy-sur-Cure, morphological Neandertal specimens are not directly dated but are contextually associated with the Châtelperronian, which contains bone points and beads. The association between Neandertals and this "transitional" assemblage has been controversial because of the lack either of a direct hominin radiocarbon date or of molecular confirmation of the Neandertal affiliation. Here we provide further evidence for a Neandertal-Châtelperronian association at the Grotte du Renne through biomolecular and chronological analysis. We identified 28 additional hominin specimens through zooarchaeology by mass spectrometry (ZooMS) screening of morphologically uninformative bone specimens from Châtelperronian layers at the Grotte du Renne. Next, we obtain an ancient hominin bone proteome through liquid chromatography-MS/MS analysis and error-tolerant amino acid sequence analysis. Analysis of this palaeoproteome allows us to provide phylogenetic and physiological information on these ancient hominin specimens. We distinguish Late Pleistocene clades within the genus Homo based on ancient protein evidence through the identification of an archaic-derived amino acid sequence for the collagen type X, alpha-1 (COL10α1) protein. We support this by obtaining ancient mtDNA sequences, which indicate a Neandertal ancestry for these specimens. Direct accelerator mass spectometry radiocarbon dating and Bayesian modeling confirm that the hominin specimens date to the Châtelperronian at the Grotte du Renne.

Using the Neandertal genome to study the evolution of small insertions and deletions in modern humans.

BACKGROUND: Small insertions and deletions occur in humans at a lower rate compared to nucleotide changes, but evolve under more constraint than nucleotide changes. While the evolution of insertions and deletions have been investigated using ape outgroups, the now available genome of a Neandertal can shed light on the evolution of indels in more recent times. RESULTS: We used the Neandertal genome together with several primate outgroup genomes to differentiate between human insertion/deletion changes that likely occurred before the split from Neandertals and those that likely arose later. Changes that pre-date the split from Neandertals show a smaller proportion of deletions than those that occurred later. The presence of a Neandertal-shared allele in Europeans or Asians but the absence in Africans was used to detect putatively introgressed indels in Europeans and Asians. A larger proportion of these variants reside in intergenic regions compared to other modern human variants, and some variants are linked to SNPs that have been associated with traits in modern humans. CONCLUSIONS: Our results are in agreement with earlier results that suggested that deletions evolve under more constraint than insertions. When considering Neandertal introgressed variants, we find some evidence that negative selection affected these variants more than other variants segregating in modern humans. Among introgressed variants we also identify indels that may influence the phenotype of their carriers. In particular an introgressed deletion associated with a decrease in the time to menarche may constitute an example of a former Neandertal-specific trait contributing to modern human phenotypic diversity.

ABAEnrichment: an R package to test for gene set expression enrichment in the adult and developing human brain.

We present ABAEnrichment, an R package that tests for expression enrichment in specific brain regions at different developmental stages using expression information gathered from multiple regions of the adult and developing human brain, together with ontologically organized structural information about the brain, both provided by the Allen Brain Atlas. We validate ABAEnrichment by successfully recovering the origin of gene sets identified in specific brain cell-types and developmental stages. AVAILABILITY AND IMPLEMENTATION: ABAEnrichment was implemented as an R package and is available under GPL (≥ 2) from the Bioconductor website (http://bioconductor.org/packages/3.3/bioc/html/ABAEnrichment.html). CONTACTS: steffi_grote@eva.mpg.de, kelso@eva.mpg.de or michael_dannemann@eva.mpg.deSupplementary information: Supplementary data are available at Bioinformatics online.

Humanized Foxp2 accelerates learning by enhancing transitions from declarative to procedural performance.

The acquisition of language and speech is uniquely human, but how genetic changes might have adapted the nervous system to this capacity is not well understood. Two human-specific amino acid substitutions in the transcription factor forkhead box P2 (FOXP2) are outstanding mechanistic candidates, as they could have been positively selected during human evolution and as FOXP2 is the sole gene to date firmly linked to speech and language development. When these two substitutions are introduced into the endogenous Foxp2 gene of mice (Foxp2(hum)), cortico-basal ganglia circuits are specifically affected. Here we demonstrate marked effects of this humanization of Foxp2 on learning and striatal neuroplasticity. Foxp2(hum/hum) mice learn stimulus-response associations faster than their WT littermates in situations in which declarative (i.e., place-based) and procedural (i.e., response-based) forms of learning could compete during transitions toward proceduralization of action sequences. Striatal districts known to be differently related to these two modes of learning are affected differently in the Foxp2(hum/hum) mice, as judged by measures of dopamine levels, gene expression patterns, and synaptic plasticity, including an NMDA receptor-dependent form of long-term depression. These findings raise the possibility that the humanized Foxp2 phenotype reflects a different tuning of corticostriatal systems involved in declarative and procedural learning, a capacity potentially contributing to adapting the human brain for speech and language acquisition.

Long-Term Balancing Selection in LAD1 Maintains a Missense Trans-Species Polymorphism in Humans, Chimpanzees, and Bonobos.

Balancing selection maintains advantageous genetic and phenotypic diversity in populations. When selection acts for long evolutionary periods selected polymorphisms may survive species splits and segregate in present-day populations of different species. Here, we investigate the role of long-term balancing selection in the evolution of protein-coding sequences in the Homo-Pan clade. We sequenced the exome of 20 humans, 20 chimpanzees, and 20 bonobos and detected eight coding trans-species polymorphisms (trSNPs) that are shared among the three species and have segregated for approximately 14 My of independent evolution. Although the majority of these trSNPs were found in three genes of the major histocompatibility locus cluster, we also uncovered one coding trSNP (rs12088790) in the gene LAD1. All these trSNPs show clustering of sequences by allele rather than by species and also exhibit other signatures of long-term balancing selection, such as segregating at intermediate frequency and lying in a locus with high genetic diversity. Here, we focus on the trSNP in LAD1, a gene that encodes for Ladinin-1, a collagenous anchoring filament protein of basement membrane that is responsible for maintaining cohesion at the dermal-epidermal junction; the gene is also an autoantigen responsible for linear IgA disease. This trSNP results in a missense change (Leucine257Proline) and, besides altering the protein sequence, is associated with changes in gene expression of LAD1.

Positive selection in the genomes of two Papua New Guinean populations at distinct altitude levels.

Highlanders and lowlanders of Papua New Guinea have faced distinct environmental stress, such as hypoxia and environment-specific pathogen exposure, respectively. In this study, we explored the top genomics regions and the candidate driver SNPs for selection in these two populations using newly sequenced whole-genomes of 54 highlanders and 74 lowlanders. We identified two candidate SNPs under selection - one in highlanders, associated with red blood cell traits and another in lowlanders, which is associated with white blood cell count - both potentially influencing the heart rate of Papua New Guineans in opposite directions. We also observed four candidate driver SNPs that exhibit linkage disequilibrium with an introgressed haplotype, highlighting the need to explore the possibility of adaptive introgression within these populations. This study reveals that the signatures of positive selection in highlanders and lowlanders of Papua New Guinea align closely with the challenges they face, which are specific to their environments.

High-throughput sequencing of the melanoma genome.

Next-generation sequencing technologies are now common for whole-genome, whole-exome and whole-transcriptome sequencing (RNA-seq) of tumors to identify point mutations, structural or copy number alterations and changes in gene expression. A substantial number of studies have already been performed for melanoma. One study analysed eight melanoma cell lines with RNA-Seq technology and identified 11 novel melanoma gene fusions. Whole-exome sequencing of seven melanoma cell lines identified overlapping gain of function mutations in MAP2K1 (MEK1) and MAP2K2 (MEK2) genes. Integrative sequencing of cutaneous melanoma metastases using different sequencing platforms revealed a new somatic point mutation in HRAS and a structural rearrangement affecting CDKN2C (a CDK4 inhibitor). These latter sequencing-based discoveries may be used to motivate the inclusion of the affected patients into clinical trials with specific signalling pathway inhibitors. Taken together, we are at the beginning of an era with new sequencing technologies providing a more comprehensive view of cancer mutational landscapes and hereby a better understanding of their pathogenesis. This will also open interesting perspectives for new treatment approaches and clinical trial designs.

Long-range regulatory effects of Neandertal DNA in modern humans.

The admixture between modern humans and Neandertals has resulted in ∼2% of the genomes of present-day non-Africans being composed of Neandertal DNA. Introgressed Neandertal DNA has been demonstrated to significantly affect the transcriptomic landscape in people today and via this molecular mechanism influence phenotype variation as well. However, little is known about how much of that regulatory impact is mediated through long-range regulatory effects that have been shown to explain ∼20% of expression variation. Here we identified 60 transcription factors (TFs) with their top cis-eQTL SNP in GTEx being of Neandertal ancestry and predicted long-range Neandertal DNA-induced regulatory effects by screening for the predicted target genes of those TFs. We show that the TFs form a significantly connected protein-protein interaction network. Among them are JUN and PRDM5, two brain-expressed TFs that have their predicted target genes enriched in regions devoid of Neandertal DNA. Archaic cis-eQTLs for the 60 TFs include multiple candidates for local adaptation, some of which show significant allele frequency increases over the last ∼10,000 years. A large proportion of the cis-eQTL-associated archaic SNPs have additional associations with various immune traits, schizophrenia, blood cell type composition and anthropometric measures. Finally, we demonstrate that our results are consistent with those of Neandertal DNA-associated empirical trans-eQTLs. Our results suggest that Neandertal DNA significantly influences regulatory networks, that its regulatory reach goes beyond the 40% of genomic sequence it still covers in present-day non-Africans and that via the investigated mechanism Neandertal DNA influences the phenotypic variation in people today.

'maskBAD'--a package to detect and remove Affymetrix probes with binding affinity differences.

BACKGROUND: Hybridization differences caused by target sequence differences can be a confounding factor in analyzing gene expression on microarrays, lead to false positives and reduce power to detect real expression differences. We prepared an R Bioconductor compatible package to detect, characterize and remove such probes in Affymetrix 3'IVT and exon-based arrays on the basis of correlation of signal intensities from probes within probe sets. RESULTS: Using completely mouse genomes we determined type 1 (false negatives) and type 2 (false positives) errors with high accuracy and we show that our method routinely outperforms previous methods. When detecting 76.2% of known SNP/indels in mouse expression data, we obtain at most 5.5% false positives. At the same level of false positives, best previous method detected 72.6%. We also show that probes with differing binding affinity both hinder differential expression detection and introduce artifacts in cancer-healthy tissue comparison. CONCLUSIONS: Detection and removal of such probes should be a routine step in Affymetrix data preprocessing. We prepared a user friendly R package, compatible with Bioconductor, that allows the filtering and improving of data from Affymetrix microarrays experiments.

Annotation of primate miRNAs by high throughput sequencing of small RNA libraries.

BACKGROUND: In addition to genome sequencing, accurate functional annotation of genomes is required in order to carry out comparative and evolutionary analyses between species. Among primates, the human genome is the most extensively annotated. Human miRNA gene annotation is based on multiple lines of evidence including evidence for expression as well as prediction of the characteristic hairpin structure. In contrast, most miRNA genes in non-human primates are annotated based on homology without any expression evidence. We have sequenced small-RNA libraries from chimpanzee, gorilla, orangutan and rhesus macaque from multiple individuals and tissues. Using patterns of miRNA expression in conjunction with a model of miRNA biogenesis we used these high-throughput sequencing data to identify novel miRNAs in non-human primates. RESULTS: We predicted 47 new miRNAs in chimpanzee, 240 in gorilla, 55 in orangutan and 47 in rhesus macaque. The algorithm we used was able to predict 64% of the previously known miRNAs in chimpanzee, 94% in gorilla, 61% in orangutan and 71% in rhesus macaque. We therefore added evidence for expression in between one and five tissues to miRNAs that were previously annotated based only on homology to human miRNAs. We increased from 60 to 175 the number miRNAs that are located in orthologous regions in humans and the four non-human primate species studied here. CONCLUSIONS: In this study we provide expression evidence for homology-based annotated miRNAs and predict de novo miRNAs in four non-human primate species. We increased the number of annotated miRNA genes and provided evidence for their expression in four non-human primates. Similar approaches using different individuals and tissues would improve annotation in non-human primates and allow for further comparative studies in the future.