New insights into the Tyrolean Iceman's origin and phenotype as inferred by whole-genome sequencing. Nat Commun. 2012;3:698. [PMID: 22426219 DOI: 10.1038/ncomms1701]

The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes

More than two hundred papers have reported genome-wide data from ancient humans. While the raw data for the vast majority are fully publicly available testifying to the commitment of the paleogenomics community to open data, formats for both raw data and meta-data differ. There is thus a need for uniform curation and a centralized, version-controlled compendium that researchers can download, analyze, and reference. Since 2019, we have been maintaining the Allen Ancient DNA Resource (AADR), which aims to provide an up-to-date, curated version of the world's published ancient human DNA data, represented at more than a million single nucleotide polymorphisms (SNPs) at which almost all ancient individuals have been assayed. The AADR has gone through six public releases at the time of writing and review of this manuscript, and crossed the threshold of >10,000 individuals with published genome-wide ancient DNA data at the end of 2022. This note is intended as a citable descriptor of the AADR.

Stable population structure in Europe since the Iron Age, despite high mobility

Ancient DNA research in the past decade has revealed that European population structure changed dramatically in the prehistoric period (14,000-3000 years before present, YBP), reflecting the widespread introduction of Neolithic farmer and Bronze Age Steppe ancestries. However, little is known about how population structure changed from the historical period onward (3000 YBP - present). To address this, we collected whole genomes from 204 individuals from Europe and the Mediterranean, many of which are the first historical period genomes from their region (e.g. Armenia and France). We found that most regions show remarkable inter-individual heterogeneity. At least 7% of historical individuals carry ancestry uncommon in the region where they were sampled, some indicating cross-Mediterranean contacts. Despite this high level of mobility, overall population structure across western Eurasia is relatively stable through the historical period up to the present, mirroring geography. We show that, under standard population genetics models with local panmixia, the observed level of dispersal would lead to a collapse of population structure. Persistent population structure thus suggests a lower effective migration rate than indicated by the observed dispersal. We hypothesize that this phenomenon can be explained by extensive transient dispersal arising from drastically improved transportation networks and the Roman Empire's mobilization of people for trade, labor, and military. This work highlights the utility of ancient DNA in elucidating finer scale human population dynamics in recent history.

Genetic Predisposition of Atherosclerotic Cardiovascular Disease in Ancient Human Remains

Background

Several computed tomographic studies have shown the presence of atherosclerosis in ancient human remains. However, while it is important to understand the development of atherosclerotic cardiovascular disease (ASCVD), genetic data concerning the prevalence of the disease-associated single nucleotide polymorphisms (SNPs) in our ancestors are scarce.

Objective

For a better understanding of the role of genetics in the evolution of ASCVD, we applied an enrichment capture sequencing approach to mummified human remains from different geographic regions and time periods.

Methods

Twenty-two mummified individuals were analyzed for their genetic predisposition of ASCVD. Next-generation sequencing methods were applied to ancient DNA (aDNA) samples, including a novel enrichment approach specifically designed to capture SNPs associated with ASCVD in genome-wide association studies of modern humans.

Findings

Five out of 22 ancient individuals passed all filter steps for calculating a weighted polygenic risk score (PRS) based on 87 SNPs in 56 genes. PRSs were correlated to scores obtained from contemporary people from around the world and cover their complete range. The genetic results of the ancient individuals reflect their phenotypic results, given that the only two mummies showing calcified atherosclerotic arterial plaques on computed tomography scans are the ones exhibiting the highest calculated PRSs.

Conclusions

These data show that alleles associated with ASCVD have been widespread for at least 5,000 years. Despite some limitations due to the nature of aDNA, our approach has the potential to lead to a better understanding of the interaction between environmental and genetic influences on the development of ASCVD.

Ancestry and kinship in a Late Antiquity-Early Middle Ages cemetery in the Eastern Italian Alps

In South Tyrol (Eastern Italian Alps), during Late Antiquity-Early Middle Ages, archeological records indicate cultural hybridization among alpine groups and peoples of various origin. Using paleogenomics, we reconstructed the ancestry of 20 individuals (4th-7th cent. AD) from a cemetery to analyze whether they had heterogeneous or homogeneous ancestry and to study their social organization. The results revealed a primary genetic ancestry from southern Europe and additional ancestries from south-western, western, and northern Europe, suggesting that cultural hybridization was accompanied by complex genetic admixture. Kinship analyses found no genetic relatedness between the only two individuals buried with grave goods. Instead, a father-son pair was discovered in one multiple grave, together with unrelated individuals and one possible non-local female. These genetic findings indicate the presence of a high social status familia, which is supported by the cultural materials and the proximity of the grave to the most sacred area of the church.

High-coverage genome of the Tyrolean Iceman reveals unusually high Anatolian farmer ancestry

The Tyrolean Iceman is known as one of the oldest human glacier mummies, directly dated to 3350-3120 calibrated BCE. A previously published low-coverage genome provided novel insights into European prehistory, despite high present-day DNA contamination. Here, we generate a high-coverage genome with low contamination (15.3×) to gain further insights into the genetic history and phenotype of this individual. Contrary to previous studies, we found no detectable Steppe-related ancestry in the Iceman. Instead, he retained the highest Anatolian-farmer-related ancestry among contemporaneous European populations, indicating a rather isolated Alpine population with limited gene flow from hunter-gatherer-ancestry-related populations. Phenotypic analysis revealed that the Iceman likely had darker skin than present-day Europeans and carried risk alleles associated with male-pattern baldness, type 2 diabetes, and obesity-related metabolic syndrome. These results corroborate phenotypic observations of the preserved mummified body, such as high pigmentation of his skin and the absence of hair on his head.

Bioarchaeological and paleogenomic profiling of the unusual Neolithic burial from Grotta di Pietra Sant'Angelo (Calabria, Italy)

The Neolithic burial of Grotta di Pietra Sant'Angelo (CS) represents a unique archaeological finding for the prehistory of Southern Italy. The unusual placement of the inhumation at a rather high altitude and far from inhabited areas, the lack of funerary equipment and the prone deposition of the body find limited similarities in coeval Italian sites. These elements have prompted wider questions on mortuary customs during the prehistory of Southern Italy. This atypical case requires an interdisciplinary approach aimed to build an integrated bioarchaeological profile of the individual. The paleopathological investigation of the skeletal remains revealed the presence of numerous markers that could be associated with craft activities, suggesting possible interpretations of the individual's lifestyle. CT analyses, carried out on the maxillary bones, showed the presence of a peculiar type of dental wear, but also a good density of the bone matrix. Biomolecular and micromorphological analyses of dental calculus highlight the presence of a rich Neolithic-like oral microbiome, the composition of which is consistent with the presence pathologies. Finally, paleogenomic data obtained from the individual were compared with ancient and modern Mediterranean populations, including unpublished high-resolution genome-wide data for 20 modern inhabitants of the nearby village of San Lorenzo Bellizzi, which provided interesting insights into the biodemographic landscape of the Neolithic in Southern Italy.

Revealing Joseon period People's single nucleotide polymorphism associated with lactase gene by ancient DNA analysis of human remains from archaeological sites in Korea

Lactase non-persistence (LNP), one of the causes of lactose intolerance, is related to lactase gene associated single nucleotide polymorphisms (SNPs). Since the frequency of LNP varies by ethnic group and country, the research to reveal the presence or absence of LNP for specific people has been conducted worldwide. However, in East Asia, the study of lactase gene associated SNPs have not been sufficiently examined so far using ancient human specimens from archaeological sites. In our study of Joseon period human remains (n=14), we successfully revealed genetic information of lactase gene associated SNPs (rs1679771596, rs41525747, rs4988236, rs4988235, rs41380347, rs869051967, rs145946881 and rs182549), further confirming that as for eight SNPs, the pre-modern Korean people had a lactase non-persistent genotype. Our report contributes to the establishment of LNP associated SNP analysis technique that can be useful in forthcoming studies on human bones and mummy samples from East Asian archaeological sites.

Methodological Changes in the Field of Paleogenetics

Paleogenetics has significantly changed since its inception almost forty years ago. Initially, molecular techniques available to the researchers offered minimal possibilities for ancient DNA analysis. The subsequent expansion of the scientific tool cabinet allowed for more remarkable achievements, combined has with the newfound popularity of this budding field of science. Finally, a breakthrough was made with the development of next-generation sequencing (NGS) technologies and the update of DNA isolation protocols, through which even very fragmented aDNA samples could be used to sequence whole genomes. In this paper, we review the achievements made thus far and compare the methodologies utilized in this field of science, discussing their benefits and challenges.

Morphological and Tissue Characterization with 3D Reconstruction of a 350-Year-Old Austrian Ardea purpurea Glacier Mummy

Glaciers are dwindling archives, releasing animal mummies preserved in the ice for centuries due to climate changes. As preservation varies, residual soft tissues may differently expand the biological information content of such mummies. DNA studies have proven the possibility of extracting and analyzing DNA preserved in skeletal residuals and sediments for hundreds or thousands of years. Paleoradiology is the method of choice as a non-destructive tool for analyzing mummies, including micro-computed tomography (micro-CT) and magnetic resonance imaging (MRI). Together with radiocarbon dating, histo-anatomical analyses, and DNA sequencing, these techniques were employed to identify a 350-year-old Austrian Ardea purpurea glacier mummy from the Ötztal Alps. Combining these techniques proved to be a robust methodological concept for collecting inaccessible information regarding the structural organization of the mummy. The variety of methodological approaches resulted in a distinct picture of the morphological patterns of the glacier animal mummy. The BLAST search in GenBank resulted in a 100% and 98.7% match in the cytb gene sequence with two entries of the species Purple heron (Ardea purpurea; Accession number KJ941160.1 and KJ190948.1) and a 98% match with the same species for the 16 s sequence (KJ190948.1), which was confirmed by the anatomic characteristics deduced from micro-CT and MRI.

Systemic lupus erythematosus and cardiovascular disease

The prognosis in systemic lupus erythematosus (SLE) has improved due to better treatment and care, but cardiovascular disease (CVD) still remains an important clinical problem, since the risk of CVD in SLE is much higher than among controls. Atherosclerosis is the main cause of CVD in the general population, and in SLE, increased atherosclerosis, especially the prevalence of atherosclerotic plaques, has been demonstrated. Atherosclerosis is an inflammatory condition, where immunity plays an important role. Interestingly, oxidized low-density lipoprotein, defective clearance of dead cells, and inflammation, with a pro-inflammatory T-cell profile are characteristics of both atherosclerosis and SLE. In addition to atherosclerosis as an underlying cause of CVD in SLE, there are also other non-mutually exclusive mechanisms, and the most important of these are antiphospholipid antibodies (aPL) leading to the antiphospholipid antibody syndrome with both arterial and venous thrombosis. aPL can cause direct pro-inflammatory and prothrombotic effects on endothelial and other cells and also interfere with the coagulation, for example, by inhibiting annexin A5 from its antithrombotic and protective effects. Antibodies against phosphorylcholine (anti-PC) and other small lipid-related epitopes, sometimes called natural antibodies, are negatively associated with CVD and atherosclerosis in SLE. Taken together, a combination of traditional risk factors such as hypertension and dyslipidemia, and nontraditional ones, especially aPL, inflammation, and low anti-PC are implicated in the increased risk of CVD in SLE. Close monitoring of both traditional risk factors and nontraditional ones, including treatment of disease manifestations, not lest renal disease in SLE, is warranted.

CONGA: Copy number variation genotyping in ancient genomes and low-coverage sequencing data

To date, ancient genome analyses have been largely confined to the study of single nucleotide polymorphisms (SNPs). Copy number variants (CNVs) are a major contributor of disease and of evolutionary adaptation, but identifying CNVs in ancient shotgun-sequenced genomes is hampered by typical low genome coverage (<1×) and short fragments (<80 bps), precluding standard CNV detection software to be effectively applied to ancient genomes. Here we present CONGA, tailored for genotyping CNVs at low coverage. Simulations and down-sampling experiments suggest that CONGA can genotype deletions >1 kbps with F-scores >0.75 at ≥1×, and distinguish between heterozygous and homozygous states. We used CONGA to genotype 10,002 outgroup-ascertained deletions across a heterogenous set of 71 ancient human genomes spanning the last 50,000 years, produced using variable experimental protocols. A fraction of these (21/71) display divergent deletion profiles unrelated to their population origin, but attributable to technical factors such as coverage and read length. The majority of the sample (50/71), despite originating from nine different laboratories and having coverages ranging from 0.44×-26× (median 4×) and average read lengths 52-121 bps (median 69), exhibit coherent deletion frequencies. Across these 50 genomes, inter-individual genetic diversity measured using SNPs and CONGA-genotyped deletions are highly correlated. CONGA-genotyped deletions also display purifying selection signatures, as expected. CONGA thus paves the way for systematic CNV analyses in ancient genomes, despite the technical challenges posed by low and variable genome coverage.

Ancient DNA analysis from epoxy resin Biodur®-embedded bones

For microscopic investigation, archaeological bone samples are often embedded in Biodur^® epoxy resin. This study wants to test whether it is possible to extract DNA suitable for PCR amplification from this sample type. For eight individuals a set of samples - each consisting of a Biodur-embedded femur sample, a native femur sample and a control sample of different anatomical origin - were submitted to organic DNA extraction. The extraction success was tested by autosomal short tandem repeat amplification. Seven out of eight Biodur-embedded femur samples revealed successful amplification results. If Biodur-embedded bone material exists from earlier microscopic investigations, our results encourage the use of this sample type as a source for genetic research.

Complex interactions between p.His558Arg and linked variants in the sodium voltage-gated channel alpha subunit 5 (Na V 1.5)

Common genetic polymorphisms may modify the phenotypic outcome when co-occurring with a disease-causing variant, and therefore understanding their modulating role in health and disease is of great importance. The polymorphic p.His558Arg variant of the sodium voltage-gated channel alpha subunit 5 (Na _V 1.5) encoded by the SCN5A gene is a case in point, as several studies have shown it can modify the clinical phenotype in a number of cardiac diseases. To evaluate the genetic backgrounds associated with this modulating effect, we reanalysed previous electrophysiological findings regarding the p.His558Arg variant and further assessed its patterns of genetic diversity in human populations. The Na _V 1.5 p.His558Arg variant was found to be in linkage disequilibrium with six other polymorphic variants that previously were also associated with cardiac traits in GWAS analyses. On account of this, incongruent reports that Arg558 allele can compensate, aggravate or have no effect on Na _V 1.5, likely might have arose due to a role of p.His558Arg depending on the additional linked variants. Altogether, these results indicate a major influence of the epistatic interactions between SCN5A variants, revealing also that phenotypic severity may depend on the polymorphic background associated to each individual genome.

Dairying, diseases and the evolution of lactase persistence in Europe

In European and many African, Middle Eastern and southern Asian populations, lactase persistence (LP) is the most strongly selected monogenic trait to have evolved over the past 10,000 years1. Although the selection of LP and the consumption of prehistoric milk must be linked, considerable uncertainty remains concerning their spatiotemporal configuration and specific interactions2,3. Here we provide detailed distributions of milk exploitation across Europe over the past 9,000 years using around 7,000 pottery fat residues from more than 550 archaeological sites. European milk use was widespread from the Neolithic period onwards but varied spatially and temporally in intensity. Notably, LP selection varying with levels of prehistoric milk exploitation is no better at explaining LP allele frequency trajectories than uniform selection since the Neolithic period. In the UK Biobank4,5 cohort of 500,000 contemporary Europeans, LP genotype was only weakly associated with milk consumption and did not show consistent associations with improved fitness or health indicators. This suggests that other reasons for the beneficial effects of LP should be considered for its rapid frequency increase. We propose that lactase non-persistent individuals consumed milk when it became available but, under conditions of famine and/or increased pathogen exposure, this was disadvantageous, driving LP selection in prehistoric Europe. Comparison of model likelihoods indicates that population fluctuations, settlement density and wild animal exploitation-proxies for these drivers-provide better explanations of LP selection than the extent of milk exploitation. These findings offer new perspectives on prehistoric milk exploitation and LP evolution.

Ancient Maltese genomes and the genetic geography of Neolithic Europe

Archaeological consideration of maritime connectivity has ranged from a biogeographical perspective that considers the sea as a barrier to a view of seaways as ancient highways that facilitate exchange. Our results illustrate the former. We report three Late Neolithic human genomes from the Mediterranean island of Malta that are markedly enriched for runs of homozygosity, indicating inbreeding in their ancestry and an effective population size of only hundreds, a striking illustration of maritime isolation in this agricultural society. In the Late Neolithic, communities across mainland Europe experienced a resurgence of hunter-gatherer ancestry, pointing toward the persistence of different ancestral strands that subsequently admixed. This is absent in the Maltese genomes, giving a further indication of their genomic insularity. Imputation of genome-wide genotypes in our new and 258 published ancient individuals allowed shared identity-by-descent segment analysis, giving a fine-grained genetic geography of Neolithic Europe. This highlights the differentiating effects of seafaring Mediterranean expansion and also island colonization, including that of Ireland, Britain, and Orkney. These maritime effects contrast profoundly with a lack of migratory barriers in the establishment of Central European farming populations from Anatolia and the Balkans.

Assessing temporal and geographic contacts across the Adriatic Sea through the analysis of genome-wide data from Southern Italy

Southern Italy was characterised by a complex prehistory that started with different Palaeolithic cultures, later followed by the Neolithization and the demic dispersal from the Pontic-Caspian Steppe during the Bronze Age. Archaeological and historical evidences point to a link between Southern Italians and the Balkans still present in modern times. To shed light on these dynamics, we analysed around 700 South Mediterranean genomes combined with informative ancient DNAs. Our findings revealed high affinities of South-Eastern Italians with modern Eastern Peloponnesians, and a closer affinity of ancient Greek genomes with those from specific regions of South Italy than modern Greek genomes. The higher similarity could be associated with a Bronze Age component ultimately originating from the Caucasus with high Iranian and Anatolian Neolithic ancestries. Furthermore, extremely differentiated allele frequencies among Northern and Southern Italy revealed putatively adapted SNPs in genes involved in alcohol metabolism, nevi features and immunological traits.

Bioarchaeological and palaeogenomic portrait of two Pompeians that died during the eruption of Vesuvius in 79 AD

The archaeological site of Pompeii is one of the 54 UNESCO World Heritage sites in Italy, thanks to its uniqueness: the town was completely destroyed and buried by a Vesuvius’ eruption in 79 AD. In this work, we present a multidisciplinary approach with bioarchaeological and palaeogenomic analyses of two Pompeian human remains from the Casa del Fabbro. We have been able to characterize the genetic profile of the first Pompeian’ genome, which has strong affinities with the surrounding central Italian population from the Roman Imperial Age. Our findings suggest that, despite the extensive connection between Rome and other Mediterranean populations, a noticeable degree of genetic homogeneity exists in the Italian peninsula at that time. Moreover, palaeopathological analyses identified the presence of spinal tuberculosis and we further investigated the presence of ancient DNA from Mycobacterium tuberculosis. In conclusion, our study demonstrates the power of a combined approach to investigate ancient humans and confirms the possibility to retrieve ancient DNA from Pompeii human remains. Our initial findings provide a foundation to promote an intensive and extensive paleogenetic analysis in order to reconstruct the genetic history of population from Pompeii, a unique archaeological site.

Serum Diamine Oxidase Values, Indicating Histamine Intolerance, Influence Lactose Tolerance Breath Test Results

Lactose intolerance (LIT) is one of the major causes of irritable bowel syndrome (IBS) spectrum complaints. Differences in inadequate lactose digestion are described as various LIT phenotypes with basically unknown pathophysiology. In LIT patients, we retrospectively assessed the effect of histamine intolerance (HIT) on expiratory hydrogen (H2) during H2 lactose breath tests. In a retrospective evaluation of charts from 402 LIT patients, 200 patients were identified as having only LIT. The other 202 LIT patients were found to additionally have diamine oxidase (DAO) values of <10 U/mL, which indicates histamine intolerance (HIT). To identify HIT, standardized questionnaires, low serum DAO values and responses to a histamine-reduced diet were used. Patients were separated into three diagnostic groups according to the result of H2 breath tests: (1) LIT, with an H2 increase of >20 parts per million (ppm), but a blood glucose (BG) increase of >20 mg/dL, (2) LIT with an H2 increase of 20 ppm in combination with a BG increase of <20 mg/dL, and (3) LIT with an exhaled H2 increase of <20 ppm and BG increase of <20 mg/dL. Pairwise comparison with the Kruskal Wallis test was used to compare the areas under the curve (AUC) of LIT and LIT with HIT patients. Exhaled H2 values were significantly higher in H2 > 20 ppm and BG < 20 mg/dL patients with LIT and HIT (p = 0.007). This diagnostic group also showed a significant higher number of patients (p = 0.012) and a significant higher number of patients with gastrointestinal (GI) symptoms during H2 breath tests (p < 0.001). Therefore, low serum DAO values, indicating HIT, influence results of lactose tolerance breath tests.

Genomes from Verteba cave suggest diversity within the Trypillians in Ukraine

The transition to agriculture occurred relatively late in Eastern Europe, leading researchers to debate whether it was a gradual, interactive process or a colonisation event. In the forest and forest-steppe regions of Ukraine, farming appeared during the fifth millennium BCE, associated with the Cucuteni-Trypillia cultural complex (CTCC, ~ 5000–3000 BCE). Across Europe, the Neolithisation process was highly variable across space and over time. Here, we investigate the population dynamics of early agriculturalists from the eastern forest-steppe region based on the analyses of 20 ancient genomes from the site of Verteba Cave (3935–825 cal BCE). Results reveal that the CTCC individuals’ ancestry is related to both western hunter-gatherers and Near Eastern farmers, has no local ancestry associated with Ukrainian Neolithic hunter-gatherers and has steppe ancestry. An Early Bronze Age individual has an ancestry profile related to the Yamnaya expansions but with 20% of ancestry related to the other Trypillian individuals, which suggests admixture between the Trypillians and the incoming populations carrying steppe-related ancestry. A Late Bronze Age individual dated to 980–825 cal BCE has a genetic profile indicating affinity to Beaker-related populations, detected close to 1000 years after the end of the Bell Beaker phenomenon during the third millennium BCE.

Genomic and dietary discontinuities during the Mesolithic and Neolithic in Sicily

Sicily is a key region for understanding the agricultural transition in the Mediterranean because of its central position. Here, we present genomic and stable isotopic data for 19 prehistoric Sicilians covering the Mesolithic to Bronze Age periods (10,700–4,100 yBP). We find that Early Mesolithic hunter-gatherers (HGs) from Sicily are a highly drifted lineage of the Early Holocene western European HGs, whereas Late Mesolithic HGs carry ∼20% ancestry related to northern and (south) eastern European HGs, indicating substantial gene flow. Early Neolithic farmers are genetically most similar to farmers from the Balkans and Greece, with only ∼7% of ancestry from local Mesolithic HGs. The genetic discontinuities during the Mesolithic and Early Neolithic match the changes in material culture and diet. Three outlying individuals dated to ∼8,000 yBP; however, suggest that hunter-gatherers interacted with incoming farmers at Grotta dell’Uzzo, resulting in a mixed economy and diet for a brief interlude at the Mesolithic-Neolithic transition.

•

Genetic transition between Early Mesolithic and Late Mesolithic hunter-gatherers

•

A near-complete genetic turnover during the Mesolithic-Neolithic transition

•

Exchange of subsistence practices between hunter-gatherers and early farmers

An integrative skeletal and paleogenomic analysis of stature variation suggests relatively reduced health for early European farmers

Human culture, biology, and health were shaped dramatically by the onset of agriculture ∼12,000 y B.P. This shift is hypothesized to have resulted in increased individual fitness and population growth as evidenced by archaeological and population genomic data alongside a decline in physiological health as inferred from skeletal remains. Here, we consider osteological and ancient DNA data from the same prehistoric individuals to study human stature variation as a proxy for health across a transition to agriculture. Specifically, we compared “predicted” genetic contributions to height from paleogenomic data and “achieved” adult osteological height estimated from long bone measurements for 167 individuals across Europe spanning the Upper Paleolithic to Iron Age (∼38,000 to 2,400 B.P.). We found that individuals from the Neolithic were shorter than expected (given their individual polygenic height scores) by an average of −3.82 cm relative to individuals from the Upper Paleolithic and Mesolithic (P = 0.040) and −2.21 cm shorter relative to post-Neolithic individuals (P = 0.068), with osteological vs. expected stature steadily increasing across the Copper (+1.95 cm relative to the Neolithic), Bronze (+2.70 cm), and Iron (+3.27 cm) Ages. These results were attenuated when we additionally accounted for genome-wide genetic ancestry variation: for example, with Neolithic individuals −2.82 cm shorter than expected on average relative to pre-Neolithic individuals (P = 0.120). We also incorporated observations of paleopathological indicators of nonspecific stress that can persist from childhood to adulthood in skeletal remains into our model. Overall, our work highlights the potential of integrating disparate datasets to explore proxies of health in prehistory.

A unified genealogy of modern and ancient genomes

The sequencing of modern and ancient genomes from around the world has revolutionized our understanding of human history and evolution. However, the problem of how best to characterize ancestral relationships from the totality of human genomic variation remains unsolved. Here, we address this challenge with nonparametric methods that enable us to infer a unified genealogy of modern and ancient humans. This compact representation of multiple datasets explores the challenges of missing and erroneous data and uses ancient samples to constrain and date relationships. We demonstrate the power of the method to recover relationships between individuals and populations as well as to identify descendants of ancient samples. Finally, we introduce a simple nonparametric estimator of the geographical location of ancestors that recapitulates key events in human history.

Mycobiome-Host Coevolution? The Mycobiome of Ancestral Human Populations Seems to Be Different and Less Diverse Than Those of Extant Native and Urban-Industrialized Populations

Few data exist on the human gut mycobiome in relation to lifestyle, ethnicity, and dietary habits. To understand the effect of these factors on the structure of the human gut mycobiome, we analyzed sequences belonging to two extinct pre-Columbian cultures inhabiting Puerto Rico (the Huecoid and Saladoid) and compared them to coprolite samples found in Mexico and Ötzi, the Iceman’s large intestine. Stool mycobiome samples from extant populations in Peru and urban cultures from the United States were also included. The ancient Puerto Rican cultures exhibited a lower fungal diversity in comparison to the extant populations. Dissimilarity distances showed that the Huecoid gut mycobiome resembled that from ancient Mexico. Fungal genera including Aspergillus spp., Penicillium spp., Rasamsonia spp., Byssochlamys spp., Talaromyces spp., Blastomyces spp., Monascus spp., and Penicilliopsis spp. were differentially abundant in the ancient and extant populations. Despite cultural differences, certain fungal taxa were present in all samples. These results suggest that culture and diet may impact the gut mycobiome and emphasize that modern lifestyles could be associated with the alteration of gut mycobiome diversity. The present study presents data on ancient and extant human gut mycobiomes in terms of lifestyle, ethnicity, and diet in the Americas.

The Genetic Origin of Daunians and the Pan-Mediterranean Southern Italian Iron Age Context

The geographical location and shape of Apulia, a narrow land stretching out in the sea at the South of Italy, made this region a Mediterranean crossroads connecting Western Europe and the Balkans. Such movements culminated at the beginning of the Iron Age with the Iapygian civilization which consisted of three cultures: Peucetians, Messapians, and Daunians. Among them, the Daunians left a peculiar cultural heritage, with one-of-a-kind stelae and pottery, but, despite the extensive archaeological literature, their origin has been lost to time. In order to shed light on this and to provide a genetic picture of Iron Age Southern Italy, we collected and sequenced human remains from three archaeological sites geographically located in Northern Apulia (the area historically inhabited by Daunians) and radiocarbon dated between 1157 and 275 calBCE. We find that Iron Age Apulian samples are still distant from the genetic variability of modern-day Apulians, they show a degree of genetic heterogeneity comparable with the cosmopolitan Republican and Imperial Roman civilization, even though a few kilometers and centuries separate them, and they are well inserted into the Iron Age Pan-Mediterranean genetic landscape. Our study provides for the first time a window on the genetic make-up of pre-Roman Apulia, whose increasing connectivity within the Mediterranean landscape, would have contributed to laying the foundation for modern genetic variability. In this light, the genetic profile of Daunians may be compatible with an at least partial autochthonous origin, with plausible contributions from the Balkan peninsula.

Predicting Physical Appearance from DNA Data-Towards Genomic Solutions

The idea of forensic DNA intelligence is to extract from genomic data any information that can help guide the investigation. The clues to the externally visible phenotype are of particular practical importance. The high heritability of the physical phenotype suggests that genetic data can be easily predicted, but this has only become possible with less polygenic traits. The forensic community has developed DNA-based predictive tools by employing a limited number of the most important markers analysed with targeted massive parallel sequencing. The complexity of the genetics of many other appearance phenotypes requires big data coupled with sophisticated machine learning methods to develop accurate genomic predictors. A significant challenge in developing universal genomic predictive methods will be the collection of sufficiently large data sets. These should be created using whole-genome sequencing technology to enable the identification of rare DNA variants implicated in phenotype determination. It is worth noting that the correctness of the forensic sketch generated from the DNA data depends on the inclusion of an age factor. This, however, can be predicted by analysing epigenetic data. An important limitation preventing whole-genome approaches from being commonly used in forensics is the slow progress in the development and implementation of high-throughput, low DNA input sequencing technologies. The example of palaeoanthropology suggests that such methods may possibly be developed in forensics.

Molecular Detection of Tick-Borne Pathogens in American Bison (Bison bison) at El Uno Ecological Reserve, Janos, Chihuahua, Mexico

American bison (Bison bison) is listed as near-threatened and in danger of extinction in Mexico. Recent studies have demonstrated the presence of several emerging pathogens at the Janos Biosphere Reserve (JBR), inhabited by one wild herd of American bison. Blood samples were collected from 26 American bison in the JBR. We tested for the presence of Anaplasma marginale, Babesia bigemina, B. bovis, Borrelia burgdorferi sensu lato, and Rickettsia rickettsii DNA using nested and semi-nested PCR protocols performing duplicates in two different laboratories. Results showed three animals (11.5%) positive for B. burgdorferi s. l., three more (11.5%) for Rickettsia rickettsii, and four (19.2%) for B. bovis. Two individuals were co-infected with B. burgdorferi s. l. and B. bovis. We found no animals positive for A. marginale and B. bigemina. This is the first report in America of R. rickettsii in American bison. American bison has been described as an important reservoir for pathogens of zoonotic and veterinary importance; thus, the presence of tick-borne pathogen DNA in the JBR American bison indicates the importance of continuous wildlife health surveys.

Whole-exome sequencing of the mummified remains of Cangrande della Scala (1291-1329 CE) indicates the first known case of late-onset Pompe disease

Mummified remains of relevant historical figures are nowadays an important source of information to retrace data concerning their private life and health, especially when historical archives are not available. Next-generation-sequencing was proved to be a valuable tool to unravel the characteristics of these individuals through their genetic heritage. Using the strictest criteria currently available for the validation of ancient DNA sequences, whole-genome and whole-exome sequencing were generated from the mummy remains of an Italian nobleman died almost 700 years ago, Cangrande della Scala. While its genome sequencing could not yield sufficient coverage for in depth investigation, exome sequencing could overcome the limitations of this approach to achieve significantly high coverage on coding regions, thus allowing to perform the first extensive exome analysis of a mummy genome. Similar to a standard "clinical exome analysis" conducted on modern DNA, an in-depth variant annotation, high-quality filtering and interpretation was performed, leading to the identification of a genotype associated with late-onset Pompe disease (glycogen storage disease type II). This genetic diagnosis was concordant with the limited clinical history available for Cangrande della Scala, who likely represents the earliest known case of this autosomal recessive metabolic disorder.

An etiology of human modernity

Following the refutation of the replacement hypothesis, which had proposed that a ‘superior’ hominin species arose in Africa and replaced all other humans existing at the time, the auto-domestication hypothesis remains the only viable explanation for the relatively abrupt change from robust to gracile humans in the Late Pleistocene. It invokes the incidental institution of the domestication syndrome in humans, most probably by newly introduced cultural practices. It also postulates that the induction of exograms compensated for the atrophy of the brain caused by domestication. This new explanation of the origins of modernity in humans elucidates practically all its many aspects, in stark contrast to the superseded replacement hypothesis, which explained virtually nothing. The first results of the domestication syndrome’s genetic exploration have become available in recent years, and they endorse the human self-domestication hypothesis.

The origin and legacy of the Etruscans through a 2000-year archeogenomic time transect

The origin, development, and legacy of the enigmatic Etruscan civilization from the central region of the Italian peninsula known as Etruria have been debated for centuries. Here we report a genomic time transect of 82 individuals spanning almost two millennia (800 BCE to 1000 CE) across Etruria and southern Italy. During the Iron Age, we detect a component of Indo-European–associated steppe ancestry and the lack of recent Anatolian-related admixture among the putative non–Indo-European–speaking Etruscans. Despite comprising diverse individuals of central European, northern African, and Near Eastern ancestry, the local gene pool is largely maintained across the first millennium BCE. This drastically changes during the Roman Imperial period where we report an abrupt population-wide shift to ~50% admixture with eastern Mediterranean ancestry. Last, we identify northern European components appearing in central Italy during the Early Middle Ages, which thus formed the genetic landscape of present-day Italian populations.

From old markers to next generation: reconstructing the history of the peopling of Sardinia

CONTEXT

For many years the Sardinian population has been the object of numerous studies because of its unique genetic structure. Despite the extreme abundance of papers, various aspects of the peopling and genetic structure of Sardinia still remain uncertain and sometimes controversial.

OBJECTIVE

We reviewed what has emerged from different studies, focussing on some still open questions, such as the origin of Sardinians, their relationship with the Corsican population, and the intra-regional genetic heterogeneity.

METHODS

The various issues have been addressed through the analysis of classical markers, molecular markers and, finally, genomic data through next generation sequencing.

RESULTS AND CONCLUSIONS

Although the most ancient human remains date back to the end of the Palaeolithic, Mesolithic populations brought founding lineages that left evident traces in the modern population. Then, with the Neolithic, the island underwent an important demographic expansion. Subsequently, isolation and genetic drift contributed to maintain a significant genetic heterogeneity, but preserving the overall homogeneity on a regional scale. At the same time, isolation and genetic drift contributed to differentiate Sardinia from Corsica, which saw an important gene flow from the mainland. However, the isolation did not prevent gene flow from the neighbouring populations whose contribution are still recognisable in the genome of Sardinians.

Paleogenomics of the prehistory of Europe: human migrations, domestication and disease

A substantial portion of ancient DNA research has been centred on understanding European populations' origin and evolution. A rchaeological evidence has already shown that the peopling of Europe involved an intricate pattern of demic and/or cultural diffusion since the Upper Palaeolithic, which became more evident during the Neolithic and Bronze Age periods. However, ancient DNA data has been crucial in determining if cultural changes occurred due to the movement of ideas or people. With the advent of next-generation sequencing and population-based paleogenomic research, ancient DNA studies have been directed not only at the study of continental human migrations, but also to the detailed analysis of particular archaeological sites, the processes of domestication, or the spread of disease during prehistoric times. With this vast paleogenomic effort added to a proper archaeological contextualisation of results, a deeper understanding of Europe's peopling is starting to emanate.

Through 40,000 years of human presence in Southern Europe: the Italian case study

The Italian Peninsula, a natural pier across the Mediterranean Sea, witnessed intricate population events since the very beginning of the human occupation in Europe. In the last few years, an increasing number of modern and ancient genomes from the area have been published by the international research community. This genomic perspective started unveiling the relevance of Italy to understand the post-Last Glacial Maximum (LGM) re-peopling of Europe, the earlier phase of the Neolithic westward migrations, and its linking role between Eastern and Western Mediterranean areas after the Iron Age. However, many open questions are still waiting for more data to be addressed in full. With this review, we summarize the current knowledge emerging from the available ancient Italian individuals and, by re-analysing them all at once, we try to shed light on the avenues future research in the area should cover. In particular, open questions concern (1) the fate of pre-Villabruna Europeans and to what extent their genomic components were absorbed by the post-LGM hunter-gatherers; (2) the role of Sicily and Sardinia before LGM; (3) to what degree the documented genetic structure within the Early Neolithic settlers can be described as two separate migrations; (4) what are the population events behind the marked presence of an Iranian Neolithic-like component in Bronze Age and Iron Age Italian and Southern European samples.

Reconstructing genetic histories and social organisation in Neolithic and Bronze Age Croatia

Ancient DNA studies have revealed how human migrations from the Neolithic to the Bronze Age transformed the social and genetic structure of European societies. Present-day Croatia lies at the heart of ancient migration routes through Europe, yet our knowledge about social and genetic processes here remains sparse. To shed light on these questions, we report new whole-genome data for 28 individuals dated to between ~ 4700 BCE-400 CE from two sites in present-day eastern Croatia. In the Middle Neolithic we evidence first cousin mating practices and strong genetic continuity from the Early Neolithic. In the Middle Bronze Age community that we studied, we find multiple closely related males suggesting a patrilocal social organisation. We also find in that community an unexpected genetic ancestry profile distinct from individuals found at contemporaneous sites in the region, due to the addition of hunter-gatherer-related ancestry. These findings support archaeological evidence for contacts with communities further north in the Carpathian Basin. Finally, an individual dated to Roman times exhibits an ancestry profile that is broadly present in the region today, adding an important data point to the substantial shift in ancestry that occurred in the region between the Bronze Age and today.

Novel mutations in the Myo5a gene cause a dilute coat color phenotype in mice

C57BL/6 laboratory mice usually show black coat color. We observed a dilute (gray) coat color phenotype in progenies of two C57BL/6 mice. This phenotype is inherited in an autosomal recessive mode. To uncover the molecular mechanism underlying this naturally occurring phenotypic variation, we performed whole-genome sequencing (25×) on 10 offspring of the two founder mice. The whole-genome DNA sequencing and additional RNA-Seq data reveal that Myo5a is the gene responsible for the coat color dilution in C57BL/6 mice, and novel mutations in the Myo5a gene are likely causal. We further performed reverse transcription-quantitative PCR, and showed increased expression of truncated Myo5a transcripts encoding dysfunctional proteins and decreased expression of Myo5a full-length transcripts encoding functional proteins in mutant individuals. The decrease in full-length messenger RNA abundance was accompanied by reduced Myo5a protein level and deficient melanosome transport, a potential mechanistic link between the Myo5a mutations and the dilute color phenotype. This study not only advances our understanding of the molecular mechanisms of pigmentation in mice, but also provides a typical case of deciphering the molecular basis of phenotypic variation in mice by genomic analyses and subsequent functional work.

[The Iceman : Life scenarios and pathological findings from 30 years of research on the glacier mummy "Ötzi"]

The comprehensive investigation of the excellently preserved mummy of Ötzi, the Iceman, and his equipment over the last 30 years has provided a wealth of information about the life and disease of this late Neolithic individual. This research has indicated that his origin was from a local southern Alpine population, that he grew up in the valleys of the Southern Alps, and that he had considerable local mobility. He had well-balanced nutrition with a mixed vegetable and animal diet. He was very mobile in the alpine terrain and of athletic constitution. The Iceman suffered from mild to moderate degenerative joint disease primarily of the right hip joint, slight spondylosis of the cervical and lumbar spine, a minor focal (premature) arteriosclerosis, lung anthracosis and possibly silicosis, previous pleuritic inflammation (possibly of post-specific origin), intestinal infections of the stomach by Helicobacter pylori and Trichuris trichiura worm infestation in the intestines, a mild osteomalacia of cancellous bone, and diverse pathologies of his teeth with dental caries and periodontitis, as well as hair anomalies. The presence of borreliosis is still under debate. As potential remedies, the Iceman carried some anthelmintic substances with him: a birch polypore and an anthelmintic fern. The numerous tattoos may also have had therapeutic effects. Finally, the last days of Ötzi could be reconstructed quite precisely: his gastrointestinal content indicates that the Iceman moved from Alpine heights to a lower location and then again up to the glacier region where he died. During this journey he encountered two attacks: the first, several days before his death, lead to a stabbing wound in his right hand; the second was an arrow hit that wounded the Iceman lethally at his left axilla by laceration of the subclavian artery.

Evidence of the interplay of genetics and culture in Ethiopia

The rich linguistic, ethnic and cultural diversity of Ethiopia provides an unprecedented opportunity to understand the level to which cultural factors correlate with-and shape-genetic structure in human populations. Using primarily new genetic variation data covering 1,214 Ethiopians representing 68 different ethnic groups, together with information on individuals' birthplaces, linguistic/religious practices and 31 cultural practices, we disentangle the effects of geographic distance, elevation, and social factors on the genetic structure of Ethiopians today. We provide evidence of associations between social behaviours and genetic differences among present-day peoples. We show that genetic similarity is broadly associated with linguistic affiliation, but also identify pronounced genetic similarity among groups from disparate language classifications that may in part be attributable to recent intermixing. We also illustrate how groups reporting the same culture traits are more genetically similar on average and show evidence of recent intermixing, suggesting that shared cultural traits may promote admixture. In addition to providing insights into the genetic structure and history of Ethiopia, we identify the most important cultural and geographic predictors of genetic differentiation and provide a resource for designing sampling protocols for future genetic studies involving Ethiopians.

rasterdiv-An Information Theory tailored R package for measuring ecosystem heterogeneity from space: To the origin and back

Ecosystem heterogeneity has been widely recognized as a key ecological indicator of several ecological functions, diversity patterns and change, metapopulation dynamics, population connectivity or gene flow.In this paper, we present a new R package-rasterdiv-to calculate heterogeneity indices based on remotely sensed data. We also provide an ecological application at the landscape scale and demonstrate its power in revealing potentially hidden heterogeneity patterns.The rasterdiv package allows calculating multiple indices, robustly rooted in Information Theory, and based on reproducible open-source algorithms.

Ancient genomes reveal structural shifts after the arrival of Steppe-related ancestry in the Italian Peninsula

Across Europe, the genetics of the Chalcolithic/Bronze Age transition is increasingly characterized in terms of an influx of Steppe-related ancestry. The effect of this major shift on the genetic structure of populations in the Italian Peninsula remains underexplored. Here, genome-wide shotgun data for 22 individuals from commingled cave and single burials in Northeastern and Central Italy dated between 3200 and 1500 BCE provide the first genomic characterization of Bronze Age individuals (n = 8; 0.001-1.2× coverage) from the central Italian Peninsula, filling a gap in the literature between 1950 and 1500 BCE. Our study confirms a diversity of ancestry components during the Chalcolithic and the arrival of Steppe-related ancestry in the central Italian Peninsula as early as 1600 BCE, with this ancestry component increasing through time. We detect close patrilineal kinship in the burial patterns of Chalcolithic commingled cave burials and a shift away from this in the Bronze Age (2200-900 BCE) along with lowered runs of homozygosity, which may reflect larger changes in population structure. Finally, we find no evidence that the arrival of Steppe-related ancestry in Central Italy directly led to changes in frequency of 115 phenotypes present in the dataset, rather that the post-Roman Imperial period had a stronger influence, particularly on the frequency of variants associated with protection against Hansen's disease (leprosy). Our study provides a closer look at local dynamics of demography and phenotypic shifts as they occurred as part of a broader phenomenon of widespread admixture during the Chalcolithic/Bronze Age transition.

The evolution of skin pigmentation-associated variation in West Eurasia

Some of the genes responsible for the evolution of light skin pigmentation in Europeans show signals of positive selection in present-day populations. Recently, genome-wide association studies have highlighted the highly polygenic nature of skin pigmentation. It is unclear whether selection has operated on all of these genetic variants or just a subset. By studying variation in over a thousand ancient genomes from West Eurasia covering 40,000 y, we are able to study both the aggregate behavior of pigmentation-associated variants and the evolutionary history of individual variants. We find that the evolution of light skin pigmentation in Europeans was driven by frequency changes in a relatively small fraction of the genetic variants that are associated with variation in the trait today.

Skin pigmentation is a classic example of a polygenic trait that has experienced directional selection in humans. Genome-wide association studies have identified well over a hundred pigmentation-associated loci, and genomic scans in present-day and ancient populations have identified selective sweeps for a small number of light pigmentation-associated alleles in Europeans. It is unclear whether selection has operated on all of the genetic variation associated with skin pigmentation as opposed to just a small number of large-effect variants. Here, we address this question using ancient DNA from 1,158 individuals from West Eurasia covering a period of 40,000 y combined with genome-wide association summary statistics from the UK Biobank. We find a robust signal of directional selection in ancient West Eurasians on 170 skin pigmentation-associated variants ascertained in the UK Biobank. However, we also show that this signal is driven by a limited number of large-effect variants. Consistent with this observation, we find that a polygenic selection test in present-day populations fails to detect selection with the full set of variants. Our data allow us to disentangle the effects of admixture and selection. Most notably, a large-effect variant at SLC24A5 was introduced to Western Europe by migrations of Neolithic farming populations but continued to be under selection post-admixture. This study shows that the response to selection for light skin pigmentation in West Eurasia was driven by a relatively small proportion of the variants that are associated with present-day phenotypic variation.

How Modern Mass Spectrometry Can Solve Ancient Questions: A Multi-Omics Study of the Stomach Content of the Oldest Human Ice Mummy, the 5300-Year-Old Iceman or Oetzi

In the past 40 years, mass spectrometry has seen a stunning development regarding increased sensitivity, resolution, and accuracy, especially for biomolecule analysis. These days without any doubt mass spectrometry is the most powerful analytical tool as a standalone technique or in conjunction with separation techniques such as high-performance liquid chromatography (HPLC), gas chromatography (GC), or capillary electrophoresis (CE). It is literally used to analyze any kind of small or large molecules ranging from basic elements to metabolites, pesticides, toxins, small or large molecule drugs, oligonucleotides, peptides, proteins, and many other molecule classes.Here, various modern mass spectrometry techniques such as LC-MS , GC-MS, ICP-MS, and elemental bio-imaging are briefly described how they were used for the first complex multi-omics study of the oldest human ice mummy, the 5300-year-old Iceman or Oetzi. The study comprised of mass spectrometry-driven proteomics (protein profiling and characterization), metabolomics, lipidomics, glycomics, and metallomics.

Lyme disease: "End of the debate?"

Lyme borreliosis is a tick-borne disease that is widespread throughout the northern hemisphere. Ixodes ricinus is present throughout metropolitan France, except for the Mediterranean region. The debate revolves around whether or not a chronic form of Lyme disease exists. This controversy is not limited to France but has been reported worldwide. In France, in 2019, 24 scientific societies representing the medical disciplines most involved in Lyme disease, including the Société Française de Rhumatologie (French Rheumatology Society [SFR]) and the Société de Pathologie Infectieuse de la Langue Française (French Infectious Disease Society-SPILF), published recommendations on the management of Lyme borreliosis following a submission to the Director General of Health. These recommendations conflict with those of the Haute Autorité de Santé (HAS), a multi-specialties independent group of physician, on a key point: whether to add a new nosological entity labeled as "persistent polymorphous signs and symptoms (or syndrome) possibly due to tick bite." The creation of this new syndrome risks should increase anchoring bias, leading to the attribution of all symptoms to a possible tick bite, without considering differential diagnoses. Lyme disease has been extensively studied. Erythema migrans is the primary clinical manifestation. In the presence of nonmetabolic, nonseptic monoarthritis involving the knee or radiculitis of a lower limb during the summer, Lyme disease should be suspected. Serologic testing for Lyme disease is reliable in the case of late forms such as chronic arthritis, while the detection of Borrelia DNA in synovial fluid by PCR is inconsistent. Sometimes, the serology can be misleading in early forms such as radiculitis. Treatment is based on doxycycline for 14 days in early forms (radiculitis), or 28 days in late forms (arthritis). Arthritis can persist or recur after antibiotic therapy. The prevalence of a diffuse polyalgia syndrome (fibromyalgia) following Lyme disease does not seem to differ much from that in the general population. It is not improved by prolonged antibiotic therapy, which is therefore not recommended.

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.

Giant virus-related sequences in the 5300-year-old Ötzi mummy metagenome

Giant viruses have brought new perspectives on the virosphere. They have been increasingly described in humans, including in several metagenomic studies. Here, we searched into the metagenome of the 5300-year-old Ötzi mummy for the presence of giant virus-related sequences using MG-Digger pipeline. We found 19 reads (0.00006% of the total read number) that best matched (mean ± standard deviation (range) for e-values of 5.0E-6 ± 1.4E-6 (6.0E-5-4.0E-10) and for amino acid identity of 69.9 ± 8.7% (46.4-84.9%) and most significantly with sequences from various giant viruses, including mostly mimiviruses. This expands current knowledge on the ubiquity and relationship with humans of giant viruses.

Population structure of indigenous inhabitants of Arabia

Modern day Saudi Arabia occupies the majority of historical Arabia, which may have contributed to ancient waves of migration out of Africa. This ancient history has left a lasting imprint in the genetics of the region, including the diverse set of tribes that call Saudi Arabia their home. How these tribes relate to each other and to the world's major populations remains an unanswered question. In an attempt to improve our understanding of the population structure of Saudi Arabia, we conducted genomic profiling of 957 unrelated individuals who self-identify with 28 large tribes in Saudi Arabia. Consistent with the tradition of intra-tribal unions, the subjects showed strong clustering along tribal lines with the distance between clusters correlating with their geographical proximities in Arabia. However, these individuals form a unique cluster when compared to the world's major populations. The ancient origin of these tribal affiliations is supported by analyses that revealed little evidence of ancestral origin from within the 28 tribes. Our results disclose a granular map of population structure and have important implications for future genetic studies into Mendelian and common diseases in the region.

Evolution of Human Brain Size-Associated NOTCH2NL Genes Proceeds toward Reduced Protein Levels

Ever since the availability of genomes from Neanderthals, Denisovans, and ancient humans, the field of evolutionary genomics has been searching for protein-coding variants that may hold clues to how our species evolved over the last ∼600,000 years. In this study, we identify such variants in the human-specific NOTCH2NL gene family, which were recently identified as possible contributors to the evolutionary expansion of the human brain. We find evidence for the existence of unique protein-coding NOTCH2NL variants in Neanderthals and Denisovans which could affect their ability to activate Notch signaling. Furthermore, in the Neanderthal and Denisovan genomes, we find unusual NOTCH2NL configurations, not found in any of the modern human genomes analyzed. Finally, genetic analysis of archaic and modern humans reveals ongoing adaptive evolution of modern human NOTCH2NL genes, identifying three structural variants acting complementary to drive our genome to produce a lower dosage of NOTCH2NL protein. Because copy-number variations of the 1q21.1 locus, encompassing NOTCH2NL genes, are associated with severe neurological disorders, this seemingly contradicting drive toward low levels of NOTCH2NL protein indicates that the optimal dosage of NOTCH2NL may have not yet been settled in the human population.

ContamLD: estimation of ancient nuclear DNA contamination using breakdown of linkage disequilibrium

We report a method called ContamLD for estimating autosomal ancient DNA (aDNA) contamination by measuring the breakdown of linkage disequilibrium in a sequenced individual due to the introduction of contaminant DNA. ContamLD leverages the idea that contaminants should have haplotypes uncorrelated to those of the studied individual. Using simulated data, we confirm that ContamLD accurately infers contamination rates with low standard errors: for example, less than 1.5% standard error in cases with less than 10% contamination and 500,000 sequences covering SNPs. This method is optimized for application to aDNA, taking advantage of characteristic aDNA damage patterns to provide calibrated contamination estimates, and is available at https://github.com/nathan-nakatsuka/ContamLD .