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.

Reconstruction of ancient microbial genomes from the human gut.

Loss of gut microbial diversity1-6 in industrial populations is associated with chronic diseases7, underscoring the importance of studying our ancestral gut microbiome. However, relatively little is known about the composition of pre-industrial gut microbiomes. Here we performed a large-scale de novo assembly of microbial genomes from palaeofaeces. From eight authenticated human palaeofaeces samples (1,000-2,000 years old) with well-preserved DNA from southwestern USA and Mexico, we reconstructed 498 medium- and high-quality microbial genomes. Among the 181 genomes with the strongest evidence of being ancient and of human gut origin, 39% represent previously undescribed species-level genome bins. Tip dating suggests an approximate diversification timeline for the key human symbiont Methanobrevibacter smithii. In comparison to 789 present-day human gut microbiome samples from eight countries, the palaeofaeces samples are more similar to non-industrialized than industrialized human gut microbiomes. Functional profiling of the palaeofaeces samples reveals a markedly lower abundance of antibiotic-resistance and mucin-degrading genes, as well as enrichment of mobile genetic elements relative to industrial gut microbiomes. This study facilitates the discovery and characterization of previously undescribed gut microorganisms from ancient microbiomes and the investigation of the evolutionary history of the human gut microbiota through genome reconstruction from palaeofaeces.

Stone Age Yersinia pestis genomes shed light on the early evolution, diversity, and ecology of plague.

The bacterial pathogen Yersinia pestis gave rise to devastating outbreaks throughout human history, and ancient DNA evidence has shown it afflicted human populations as far back as the Neolithic. Y. pestis genomes recovered from the Eurasian Late Neolithic/Early Bronze Age (LNBA) period have uncovered key evolutionary steps that led to its emergence from a Yersinia pseudotuberculosis-like progenitor; however, the number of reconstructed LNBA genomes are too few to explore its diversity during this critical period of development. Here, we present 17 Y. pestis genomes dating to 5,000 to 2,500 y BP from a wide geographic expanse across Eurasia. This increased dataset enabled us to explore correlations between temporal, geographical, and genetic distance. Our results suggest a nonflea-adapted and potentially extinct single lineage that persisted over millennia without significant parallel diversification, accompanied by rapid dispersal across continents throughout this period, a trend not observed in other pathogens for which ancient genomes are available. A stepwise pattern of gene loss provides further clues on its early evolution and potential adaptation. We also discover the presence of the flea-adapted form of Y. pestis in Bronze Age Iberia, previously only identified in in the Caucasus and the Volga regions, suggesting a much wider geographic spread of this form of Y. pestis. Together, these data reveal the dynamic nature of plague's formative years in terms of its early evolution and ecology.

Reconstructing the genetic history of late Neanderthals.

Although it has previously been shown that Neanderthals contributed DNA to modern humans, not much is known about the genetic diversity of Neanderthals or the relationship between late Neanderthal populations at the time at which their last interactions with early modern humans occurred and before they eventually disappeared. Our ability to retrieve DNA from a larger number of Neanderthal individuals has been limited by poor preservation of endogenous DNA and contamination of Neanderthal skeletal remains by large amounts of microbial and present-day human DNA. Here we use hypochlorite treatment of as little as 9 mg of bone or tooth powder to generate between 1- and 2.7-fold genomic coverage of five Neanderthals who lived around 39,000 to 47,000 years ago (that is, late Neanderthals), thereby doubling the number of Neanderthals for which genome sequences are available. Genetic similarity among late Neanderthals is well predicted by their geographical location, and comparison to the genome of an older Neanderthal from the Caucasus indicates that a population turnover is likely to have occurred, either in the Caucasus or throughout Europe, towards the end of Neanderthal history. We find that the bulk of Neanderthal gene flow into early modern humans originated from one or more source populations that diverged from the Neanderthals that were studied here at least 70,000 years ago, but after they split from a previously sequenced Neanderthal from Siberia around 150,000 years ago. Although four of the Neanderthals studied here post-date the putative arrival of early modern humans into Europe, we do not detect any recent gene flow from early modern humans in their ancestry.

Ancient DNA from Guam and the peopling of the Pacific.

Humans reached the Mariana Islands in the western Pacific by ∼3,500 y ago, contemporaneous with or even earlier than the initial peopling of Polynesia. They crossed more than 2,000 km of open ocean to get there, whereas voyages of similar length did not occur anywhere else until more than 2,000 y later. Yet, the settlement of Polynesia has received far more attention than the settlement of the Marianas. There is uncertainty over both the origin of the first colonizers of the Marianas (with different lines of evidence suggesting variously the Philippines, Indonesia, New Guinea, or the Bismarck Archipelago) as well as what, if any, relationship they might have had with the first colonizers of Polynesia. To address these questions, we obtained ancient DNA data from two skeletons from the Ritidian Beach Cave Site in northern Guam, dating to ∼2,200 y ago. Analyses of complete mitochondrial DNA genome sequences and genome-wide SNP data strongly support ancestry from the Philippines, in agreement with some interpretations of the linguistic and archaeological evidence, but in contradiction to results based on computer simulations of sea voyaging. We also find a close link between the ancient Guam skeletons and early Lapita individuals from Vanuatu and Tonga, suggesting that the Marianas and Polynesia were colonized from the same source population, and raising the possibility that the Marianas played a role in the eventual settlement of Polynesia.

Combining machine learning and optimization for the operational patient-bed assignment problem.

Assigning inpatients to hospital beds impacts patient satisfaction and the workload of nurses and doctors. The assignment is subject to unknown inpatient arrivals, in particular for emergency patients. Hospitals, therefore, need to deal with uncertainty on actual bed requirements and potential shortage situations as bed capacities are limited. This paper develops a model and solution approach for solving the patient bed-assignment problem that is based on a machine learning (ML) approach to forecasting emergency patients. First, it contributes by improving the anticipation of emergency patients using ML approaches, incorporating weather data, time and dates, important local and regional events, as well as current and historical occupancy levels. Drawing on real-life data from a large case hospital, we were able to improve forecasting accuracy for emergency inpatient arrivals. We achieved up to 17% better root mean square error (RMSE) when using ML methods compared to a baseline approach relying on averages for historical arrival rates. We further show that the ML methods outperform time series forecasts. Second, we develop a new hyper-heuristic for solving real-life problem instances based on the pilot method and a specialized greedy look-ahead (GLA) heuristic. When applying the hyper-heuristic in test sets we were able to increase the objective function by up to 5.3% in comparison to the benchmark approach in [40]. A benchmark with a Genetic Algorithm shows also the superiority of the hyper-heuristic. Third, the combination of ML for emergency patient admission forecasting with advanced optimization through the hyper-heuristic allowed us to obtain an improvement of up to 3.3% on a real-life problem.

Contrasting Paternal and Maternal Genetic Histories of Thai and Lao Populations.

The human demographic history of Mainland Southeast Asia (MSEA) has not been well studied; in particular, there have been very few sequence-based studies of variation in the male-specific portions of the Y chromosome (MSY). Here, we report new MSY sequences of ∼2.3 mB from 914 males and combine these with previous data for a total of 928 MSY sequences belonging to 59 populations from Thailand and Laos who speak languages belonging to three major Mainland Southeast Asia families: Austroasiatic, Tai-Kadai, and Sino-Tibetan. Among the 92 MSY haplogroups, two main MSY lineages (O1b1a1a* [O-M95*] and O2a* [O-M324*]) contribute substantially to the paternal genetic makeup of Thailand and Laos. We also analyze complete mitochondrial DNA genome sequences published previously from the same groups and find contrasting pattern of male and female genetic variation and demographic expansions, especially for the hill tribes, Mon, and some major Thai groups. In particular, we detect an effect of postmarital residence pattern on genetic diversity in patrilocal versus matrilocal groups. Additionally, both male and female demographic expansions were observed during the early Mesolithic (∼10 ka), with two later major male-specific expansions during the Neolithic period (∼4-5 ka) and the Bronze/Iron Age (∼2.0-2.5 ka). These two later expansions are characteristic of the modern Austroasiatic and Tai-Kadai groups, respectively, consistent with recent ancient DNA studies. We simulate MSY data based on three demographic models (continuous migration, demic diffusion, and cultural diffusion) of major Thai groups and find different results from mitochondrial DNA simulations, supporting contrasting male and female genetic histories.

Comparison of extraction methods for recovering ancient microbial DNA from paleofeces.

OBJECTIVES: Paleofeces are valuable to archeologists and evolutionary biologists for their potential to yield health, dietary, and host information. As a rich source of preserved biomolecules from host-associated microorganisms, they can also provide insights into the recent evolution and changing ecology of the gut microbiome. However, there is currently no standard method for DNA extraction from paleofeces, which combine the dual challenges of complex biological composition and degraded DNA. Due to the scarcity and relatively poor preservation of paleofeces when compared with other archeological remains, it is important to use efficient methods that maximize ancient DNA (aDNA) recovery while also minimizing downstream taxonomic biases. METHODS: In this study, we use shotgun metagenomics to systematically compare the performance of five DNA extraction methods on a set of well-preserved human and dog paleofeces from Mexico (~1,300 BP). RESULTS: Our results show that all tested DNA extraction methods yield a consistent microbial taxonomic profile, but that methods optimized for ancient samples recover significantly more DNA. CONCLUSIONS: These results show promise for future studies that seek to explore the evolution of the human gut microbiome by comparing aDNA data with those generated in modern studies.

Cultural Innovations Influence Patterns of Genetic Diversity in Northwestern Amazonia.

Human populations often exhibit contrasting patterns of genetic diversity in the mtDNA and the nonrecombining portion of the Y-chromosome (NRY), which reflect sex-specific cultural behaviors and population histories. Here, we sequenced 2.3 Mb of the NRY from 284 individuals representing more than 30 Native American groups from Northwestern Amazonia (NWA) and compared these data to previously generated mtDNA genomes from the same groups, to investigate the impact of cultural practices on genetic diversity and gain new insights about NWA population history. Relevant cultural practices in NWA include postmarital residential rules and linguistic exogamy, a marital practice in which men are required to marry women speaking a different language. We identified 2,969 SNPs in the NRY sequences, only 925 of which were previously described. The NRY and mtDNA data showed different sex-specific demographic histories: female effective population size has been larger than that of males through time, which might reflect larger variance in male reproductive success. Both markers show an increase in lineage diversification beginning ∼5,000 years ago, which may reflect the intensification of agriculture, technological innovations, and the expansion of regional trade networks documented in the archaeological evidence. Furthermore, we find similar excesses of NRY versus mtDNA between-population divergence at both the local and continental scale, suggesting long-term stability of female versus male migration. We also find evidence of the impact of sociocultural practices on diversity patterns. Finally, our study highlights the importance of analyzing high-resolution mtDNA and NRY sequences to reconstruct demographic history, since this can differ considerably between sexes.

Age-Related and Heteroplasmy-Related Variation in Human mtDNA Copy Number.

The mitochondrial (mt) genome is present in many copies in human cells, and intra-individual variation in mtDNA sequences is known as heteroplasmy. Recent studies found that heteroplasmies are highly tissue-specific, site-specific, and allele-specific, however the functional implications have not been explored. This study investigates variation in mtDNA copy numbers (mtCN) in 12 different tissues obtained at autopsy from 152 individuals (ranging in age from 3 days to 96 years). Three different methods to estimate mtCN were compared: shotgun sequencing (in 4 tissues), capture-enriched sequencing (in 12 tissues) and droplet digital PCR (ddPCR, in 2 tissues). The highest precision in mtCN estimation was achieved using shotgun sequencing data. However, capture-enrichment data provide reliable estimates of relative (albeit not absolute) mtCNs. Comparisons of mtCN from different tissues of the same individual revealed that mtCNs in different tissues are, with few exceptions, uncorrelated. Hence, each tissue of an individual seems to regulate mtCN in a tissue-related rather than an individual-dependent manner. Skeletal muscle (SM) samples showed an age-related decrease in mtCN that was especially pronounced in males, while there was an age-related increase in mtCN for liver (LIV) samples. MtCN in SM samples was significantly negatively correlated with both the total number of heteroplasmic sites and with minor allele frequency (MAF) at two heteroplasmic sites, 408 and 16327. Heteroplasmies at both sites are highly specific for SM, accumulate with aging and are part of functional elements that regulate mtDNA replication. These data support the hypothesis that selection acting on these heteroplasmic sites is reducing mtCN in SM of older individuals.

Genetic structure and sex-biased gene flow in the history of southern African populations.

OBJECTIVES: We investigated the genetic history of southern African populations with a special focus on their paternal history. We reexamined previous claims that the Y-chromosome haplogroup E1b1b (E-M293) was brought to southern Africa by pastoralists from eastern Africa, and investigated patterns of sex-biased gene flow in southern Africa. MATERIALS AND METHODS: We analyzed previously published complete mtDNA genome sequences and ∼900 kb of NRY sequences from 23 populations from Namibia, Botswana, and Zambia, as well as haplogroup frequencies from a large sample of southern African populations and 23 newly genotyped Y-linked STR loci for samples assigned to haplogroup E1b1b. RESULTS: Our results support an eastern African origin for Y-chromosome haplogroup E1b1b (E-M293); however, its current distribution in southern Africa is not strongly associated with pastoralism, suggesting more complex demographic events and/or changes in subsistence practices in this region. The Bantu expansion in southern Africa had a notable genetic impact and was probably a rapid, male-dominated expansion. Our finding of a significant increase in the intensity of the sex-biased gene flow from north to south may reflect changes in the social dynamics between Khoisan and Bantu groups over time. CONCLUSIONS: Our study shows that the population history of southern Africa has been complex, with different immigrating groups mixing to different degrees with the autochthonous populations. The Bantu expansion led to heavily sex-biased admixture as a result of interactions between Khoisan females and Bantu males, with a geographic gradient which may reflect changes in the social dynamics between Khoisan and Bantu groups over time.

Direct radiocarbon dating and DNA analysis of the Darra-i-Kur (Afghanistan) human temporal bone.

The temporal bone discovered in the 1960s from the Darra-i-Kur cave in Afghanistan is often cited as one of the very few Pleistocene human fossils from Central Asia. Here we report the first direct radiocarbon date for the specimen and the genetic analyses of DNA extracted and sequenced from two areas of the bone. The new radiocarbon determination places the find to ∼4500 cal BP (∼2500 BCE) contradicting an assumed Palaeolithic age of ∼30,000 years, as originally suggested. The DNA retrieved from the specimen originates from a male individual who carried mitochondrial DNA of the modern human type. The petrous part yielded more endogenous ancient DNA molecules than the squamous part of the same bone. Molecular dating of the Darra-i-Kur mitochondrial DNA sequence corroborates the radiocarbon date and suggests that the specimen is younger than previously thought. Taken together, the results consolidate the fact that the human bone is not associated with the Pleistocene-age deposits of Darra-i-Kur; instead it is intrusive, possibly re-deposited from upper levels dating to much later periods (Neolithic). Despite its Holocene age, the Darra-i-Kur specimen is, so far, the first and only ancient human from Afghanistan whose DNA has been sequenced.

Forming B-B Bonds by the Controlled Reduction of a Tetraaryl-diborane(6).

Dimeric aryl(hydro)boranes can provide suitable platforms for the synthesis of boron-containing graphene flakes through reductive B-B coupling. Two-electron reduction of 1,2:1,2-bis(4,4'-di-tert-butyl-2,2'-biphenylylene)diborane(6) (4) with LiNaph/THF establishes a B-B σ bond but can be accompanied by substituent redistribution. In the singly rearranged product, Li2[6], only one 1,2-phenyl shift has occurred. The doubly ring-contracted product, Li2[7], consists of two 9H-9-borafluorenyl moieties that are linked via their boron atoms. When the amount of LiNaph/THF is increased to 4 equiv, Li2[6] is subsequently observed as the dominant species. Addition of 11 equiv of LiNaph/THF results in over-reduction with hydride elimination to afford the doubly boron-doped dibenzo[g,p]chrysene Li2[1]. In contrast, excess KC8 reduces 4 to the corresponding dihydro-dibenzo[g,p]chrysene, K2[5], with a trans-HB-BH core. Hydride abstraction from K2[5] with 1 equiv of 4 leads to K[8], in which the central B-B bond is bridged by a single hydrogen atom. K[8] is also obtained upon treatment of 4 with 1 equiv of KC8. All products have been characterized by multinuclear NMR spectroscopy and X-ray crystallography.

Refining the Y chromosome phylogeny with southern African sequences.

The recent availability of large-scale sequence data for the human Y chromosome has revolutionized analyses of and insights gained from this non-recombining, paternally inherited chromosome. However, the studies to date focus on Eurasian variation, and hence the diversity of early-diverging branches found in Africa has not been adequately documented. Here, we analyze over 900 kb of Y chromosome sequence obtained from 547 individuals from southern African Khoisan- and Bantu-speaking populations, identifying 232 new sequences from basal haplogroups A and B. We identify new clades in the phylogeny, an older age for the root, and substantially older ages for some individual haplogroups. Furthermore, while haplogroup B2a is traditionally associated with the spread of Bantu speakers, we find that it probably also existed in Khoisan groups before the arrival of Bantu speakers. Finally, there is pronounced variation in branch length between major haplogroups; in particular, haplogroups associated with Bantu speakers have significantly longer branches. Technical artifacts cannot explain this branch length variation, which instead likely reflects aspects of the demographic history of Bantu speakers, such as recent population expansion and an older average paternal age. The influence of demographic factors on branch length variation has broader implications both for the human Y phylogeny and for similar analyses of other species.

A preorganized ditopic borane as highly efficient one- or two-electron trap.

Reduction of the bis(9-borafluorenyl)methane 1 with excess lithium furnishes the red dianion salt Li2[1]. The corresponding dark green monoanion radical Li[1] is accessible through the comproportionation reaction between 1 and Li2[1]. EPR spectroscopy on Li[1] reveals hyperfine coupling of the unpaired electron to two magnetically equivalent boron nuclei (a((11)B) = 5.1 ± 0.1 G, a((10)B) = 1.7 ± 0.2 G). Further coupling is observed to the unique B-CH-B bridgehead proton (a((1)H) = 7.2 ± 0.2 G) and to eight aromatic protons (a((1)H) = 1.4 ± 0.1 G). According to X-ray crystallography, the B···B distances continuously decrease along the sequence 1 → [1](•-) → [1](2-) with values of 2.534(2), 2.166(4), and 1.906(3) Å, respectively. Protonation of Li2[1] leads to the cyclic borohydride species Li[1H] featuring a B-H-B two-electron-three-center bond. This result strongly indicates a nucleophilic character of the boron atoms; the reaction can also be viewed as rare example of the protonation of an element-element σ bond. According to NMR spectroscopy, EPR spectroscopy, and quantum-chemical calculations, [1](2-) represents a closed-shell singlet without any spin contamination. Detailed wave function analyses of [1](•-) and [1](2-) reveal strongly localized interactions of the two boron pz-type orbitals, with small delocalized contributions of the 9-borafluorenyl π systems. Overall, our results provide evidence for a direct B-B one-electron and two-electron bonding interaction in [1](•-) and [1](2-), respectively.

Extensive structural rearrangements upon reduction of 9H-9-borafluorene.

Common wisdom has it that organoboranes are readily oxidized. Described herein is that also their reduction can result in remarkable chemistry. Treatment of dimeric 9H-9-borafluorene with Li metal in toluene yields two strikingly different classes of compounds. One part of the sample reacts in a way similar to B2H6, thus affording an aryl(hydro)borane cluster reminiscent of the [B3H8](-) anion. The other part furnishes a dianionic boron-doped graphene flake devoid of hydrogen substituents at the boron centers and featuring a central B=B bond. A change in the solvent to THF allows an isolation of this dibenzo[g,p]chrysene analogue in good yields.

Confirmed by X-ray crystallography: the B⋠B one-electron σ bond.

Is one electron sufficient to bring about significant σ bonding between two atoms? The chemist's view on the chemical bond is usually tied to the concept of shared electron pairs, and not too much experimental evidence exists to challenge this firm belief. Whilst species with the unusual one-electron σ-bonding motif between homonuclear atoms have so far been identified mainly by spectroscopic evidence, we present herein the first crystallographic characterization, augmented by a detailed quantum-chemical validation, for a radical anion featuring a B⋅B one-electron-two-center σ bond.

Metagenomic and paleopathological analyses of a historic documented collection explore ancient dental calculus as a diagnostic tool.

Dental calculus is a microbial biofilm that contains biomolecules from oral commensals and pathogens, including those potentially related to cause of death (CoD). To assess the utility of calculus as a diagnostically informative substrate, in conjunction with paleopathological analysis, calculus samples from 39 individuals in the Smithsonian Institution's Robert J. Terry Collection with CoDs of either syphilis or tuberculosis were assessed via shotgun metagenomic sequencing for the presence of Treponema pallidum subsp. pallidum and Mycobacterium tuberculosis complex (MTBC) DNA. Paleopathological analysis revealed that frequencies of skeletal lesions associated with these diseases were partially inconsistent with diagnostic criteria. Although recovery of T. p. pallidum DNA from individuals with a syphilis CoD was elusive, MTBC DNA was identified in at least one individual with a tuberculosis CoD. The authenticity of MTBC DNA was confirmed using targeted quantitative PCR assays, MTBC genome enrichment, and in silico bioinformatic analyses; however, the lineage of the MTBC strain present could not be determined. Overall, our study highlights the utility of dental calculus for molecular detection of tuberculosis in the archaeological record and underscores the effect of museum preparation techniques and extensive handling on pathogen DNA preservation in skeletal collections.

Investigating demic versus cultural diffusion and sex bias in the spread of Austronesian languages in Vietnam.

Austronesian (AN) is the second-largest language family in the world, particularly widespread in Island Southeast Asia (ISEA) and Oceania. In Mainland Southeast Asia (MSEA), groups speaking these languages are concentrated in the highlands of Vietnam. However, our knowledge of the spread of AN-speaking populations in MSEA remains limited; in particular, it is not clear if AN languages were spread by demic or cultural diffusion. In this study, we present and analyze new data consisting of complete mitogenomes from 369 individuals and 847 Y-chromosomal single nucleotide polymorphisms (SNPs) from 170 individuals from all five Vietnamese Austronesian groups (VN-AN) and five neighboring Vietnamese Austroasiatic groups (VN-AA). We found genetic signals consistent with matrilocality in some, but not all, of the VN-AN groups. Population affinity analyses indicated connections between the AN-speaking Giarai and certain Taiwanese AN groups (Rukai, Paiwan, and Bunun). However, overall, there were closer genetic affinities between VN-AN groups and neighboring VN-AA groups, suggesting language shifts. Our study provides insights into the genetic structure of AN-speaking communities in MSEA, characterized by some contact with Taiwan and language shift in neighboring groups, indicating that the expansion of AN speakers in MSEA was a combination of cultural and demic diffusion.

Natural products from reconstructed bacterial genomes of the Middle and Upper Paleolithic.

Major advances over the past decade in the field of ancient DNA are providing access to past paleogenomic diversity, but the diverse functions and biosynthetic capabilities of this growing paleome remain largely elusive. We investigated the dental calculus of 12 Neanderthals and 52 anatomically modern humans ranging from 100,000 years ago to the present and reconstructed 459 bacterial metagenome-assembled genomes. We identified a biosynthetic gene cluster shared by seven Middle and Upper Paleolithic individuals that allows for the heterologous production of a class of previously unknown metabolites that we name "paleofurans." This paleobiotechnological approach demonstrates that viable biosynthetic machinery can be produced from the preserved genetic material of ancient organisms, allowing access to natural products from the Pleistocene and providing a promising area for natural product exploration.