Functional characterization of archaic-specific variants in mitonuclear genes: insights from comparative analysis in S. cerevisiae.

Neanderthal and Denisovan hybridisation with modern humans has generated a non-random genomic distribution of introgressed regions, the result of drift and selection dynamics. Cross-species genomic incompatibility and more efficient removal of slightly deleterious archaic variants have been proposed as selection-based processes involved in the post-hybridisation purge of archaic introgressed regions. Both scenarios require the presence of functionally different alleles across Homo species onto which selection operated differently according to which populations hosted them, but only a few of these variants have been pinpointed so far. In order to identify functionally divergent archaic variants removed in humans, we focused on mitonuclear genes, which are underrepresented in the genomic landscape of archaic humans. We searched for non-synonymous, fixed, archaic-derived variants present in mitonuclear genes, rare or absent in human populations. We then compared the functional impact of archaic and human variants in the model organism Saccharomyces cerevisiae. Notably, a variant within the mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) gene exhibited a significant decrease in respiratory activity and a substantial reduction of Cox2 levels, a proxy for mitochondrial protein biosynthesis, coupled with the accumulation of the YARS2 protein precursor and a lower amount of mature enzyme. Our work suggests that this variant is associated with mitochondrial functionality impairment, thus contributing to the purging of archaic introgression in YARS2. While different molecular mechanisms may have impacted other mitonuclear genes, our approach can be extended to the functional screening of mitonuclear genetic variants present across species and populations.

The first Miocene fossils from coastal woodlands in the southern East African Rift.

The Miocene was a key time in the evolution of African ecosystems witnessing the origin of the African apes and the isolation of eastern coastal forests through an expanding arid corridor. Until recently, however, Miocene sites from the southeastern regions of the continent were unknown. Here, we report the first Miocene fossil teeth from the shoulders of the Urema Rift in Gorongosa National Park, Mozambique. We provide the first 1) radiometric ages of the Mazamba Formation, 2) reconstructions of paleovegetation in the region based on pedogenic carbonates and fossil wood, and 3) descriptions of fossil teeth. Gorongosa is unique in the East African Rift in combining marine invertebrates, marine vertebrates, reptiles, terrestrial mammals, and fossil woods in coastal paleoenvironments. The Gorongosa fossil sites offer the first evidence of woodlands and forests on the coastal margins of southeastern Africa during the Miocene, and an exceptional assemblage of fossils including new species.

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes.

Y chromosome markers can shed light on male-specific population dynamics but for many species no such markers have been discovered and are available yet, despite the potential for recovering Y-linked loci from available genome sequences. Here, we investigated how effective available bioinformatic tools are in recovering informative Y chromosome microsatellites from whole genome sequence data. In order to do so, we initially explored a large dataset of whole genome sequences comprising individuals at various coverages belonging to different species of baboons (genus: Papio) using Y chromosome references belonging to the same genus and more distantly related species (Macaca mulatta). We then further tested this approach by recovering Y-STRs from available Theropithecus gelada genomes using Papio and Macaca Y chromosome as reference sequences. Identified loci were validated in silico by a) comparing within-species relationships of Y chromosome lineages and b) genotyping male individuals in available pedigrees. Each STR was selected not to extend in its variable region beyond 100 base pairs, so that loci can be developed for PCR-based genotyping of non-invasive DNA samples. In addition to assembling a first set of Papio and Theropithecus Y-specific microsatellite markers, we released TYpeSTeR, an easy-to-use script to identify and genotype Y chromosome STRs using population genomic data which can be modulated according to available male reference genomes and genomic data, making it widely applicable across taxa.

Genetic heritage of the Baphuthi highlights an over-ethnicized notion of "Bushman" in the Maloti-Drakensberg, southern Africa.

Using contemporary people as proxies for ancient communities is a contentious but necessary practice in anthropology. In southern Africa, the distinction between the Cape KhoeSan and eastern KhoeSan remains unclear, as ethnicity labels have been changed through time and most communities were decimated if not extirpated. The eastern KhoeSan may have had genetic distinctions from neighboring communities who speak Bantu languages and KhoeSan further away; alternatively, the identity may not have been tied to any notion of biology, instead denoting communities with a nomadic "lifeway" distinct from African agro-pastoralism. The Baphuthi of the 1800s in the Maloti-Drakensberg, southern Africa had a substantial KhoeSan constituency and a lifeway of nomadism, cattle raiding, and horticulture. Baphuthi heritage could provide insights into the history of the eastern KhoeSan. We examine genetic affinities of 23 Baphuthi to discern whether the narrative of KhoeSan descent reflects distinct genetic ancestry. Genome-wide SNP data (Illumina GSA) were merged with 52 global populations, for 160,000 SNPs. Genetic analyses show no support for a unique eastern KhoeSan ancestry distinct from other KhoeSan or southern Bantu speakers. The Baphuthi have strong affinities with early-arriving southern Bantu-speaking (Nguni) communities, as the later-arriving non-Nguni show strong evidence of recent African admixture possibly related to late-Iron Age migrations. The references to communities as "San" and "Bushman" in historic literature has often been misconstrued as notions of ethnic/biological distinctions. The terms may have reflected ambiguous references to non-sedentary polities instead, as seems to be the case for the eastern "Bushman" heritage of the Baphuthi.

The diverse genetic origins of a Classical period Greek army.

Trade and colonization caused an unprecedented increase in Mediterranean human mobility in the first millennium BCE. Often seen as a dividing force, warfare is in fact another catalyst of culture contact. We provide insight into the demographic dynamics of ancient warfare by reporting genome-wide data from fifth-century soldiers who fought for the army of the Greek Sicilian colony of Himera, along with representatives of the civilian population, nearby indigenous settlements, and 96 present-day individuals from Italy and Greece. Unlike the rest of the sample, many soldiers had ancestral origins in northern Europe, the Steppe, and the Caucasus. Integrating genetic, archaeological, isotopic, and historical data, these results illustrate the significant role mercenaries played in ancient Greek armies and highlight how participation in war contributed to continental-scale human mobility in the Classical world.

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.

Genomic variation in baboons from central Mozambique unveils complex evolutionary relationships with other Papio species.

BACKGROUND: Gorongosa National Park in Mozambique hosts a large population of baboons, numbering over 200 troops. Gorongosa baboons have been tentatively identified as part of Papio ursinus on the basis of previous limited morphological analysis and a handful of mitochondrial DNA sequences. However, a recent morphological and morphometric analysis of Gorongosa baboons pinpointed the occurrence of several traits intermediate between P. ursinus and P. cynocephalus, leaving open the possibility of past and/or ongoing gene flow in the baboon population of Gorongosa National Park. In order to investigate the evolutionary history of baboons in Gorongosa, we generated high and low coverage whole genome sequence data of Gorongosa baboons and compared it to available Papio genomes. RESULTS: We confirmed that P. ursinus is the species closest to Gorongosa baboons. However, the Gorongosa baboon genomes share more derived alleles with P. cynocephalus than P. ursinus does, but no recent gene flow between P. ursinus and P. cynocephalus was detected when available Papio genomes were analyzed. Our results, based on the analysis of autosomal, mitochondrial and Y chromosome data, suggest complex, possibly male-biased, gene flow between Gorongosa baboons and P. cynocephalus, hinting to direct or indirect contributions from baboons belonging to the "northern" Papio clade, and signal the presence of population structure within P. ursinus. CONCLUSIONS: The analysis of genome data generated from baboon samples collected in central Mozambique highlighted a complex set of evolutionary relationships with other baboons. Our results provided new insights in the population dynamics that have shaped baboon diversity.

Evaluating the Impact of Sex-Biased Genetic Admixture in the Americas through the Analysis of Haplotype Data.

A general imbalance in the proportion of disembarked males and females in the Americas has been documented during the Trans-Atlantic Slave Trade and the Colonial Era and, although less prominent, more recently. This imbalance may have left a signature on the genomes of modern-day populations characterised by high levels of admixture. The analysis of the uniparental systems and the evaluation of continental proportion ratio of autosomal and X chromosomes revealed a general sex imbalance towards males for European and females for African and Indigenous American ancestries. However, the consistency and degree of this imbalance are variable, suggesting that other factors, such as cultural and social practices, may have played a role in shaping it. Moreover, very few investigations have evaluated the sex imbalance using haplotype data, containing more critical information than genotypes. Here, we analysed genome-wide data for more than 5000 admixed American individuals to assess the presence, direction and magnitude of sex-biased admixture in the Americas. For this purpose, we applied two haplotype-based approaches, ELAI and NNLS, and we compared them with a genotype-based method, ADMIXTURE. In doing so, besides a general agreement between methods, we unravelled that the post-colonial admixture dynamics show higher complexity than previously described.

Continental-scale genomic analysis suggests shared post-admixture adaptation in the Americas.

American populations are one of the most interesting examples of recently admixed groups, where ancestral components from three major continental human groups (Africans, Eurasians and Native Americans) have admixed within the last 15 generations. Recently, several genetic surveys focusing on thousands of individuals shed light on the geography, chronology and relevance of these events. However, even though gene flow could drive adaptive evolution, it is unclear whether and how natural selection acted on the resulting genetic variation in the Americas. In this study, we analysed the patterns of local ancestry of genomic fragments in genome-wide data for ~ 6000 admixed individuals from 10 American countries. In doing so, we identified regions characterized by a divergent ancestry profile (DAP), in which a significant over or under ancestral representation is evident. Our results highlighted a series of genomic regions with DAPs associated with immune system response and relevant medical traits, with the longest DAP region encompassing the human leukocyte antigen locus. Furthermore, we found that DAP regions are enriched in genes linked to cancer-related traits and autoimmune diseases. Then, analysing the biological impact of these regions, we showed that natural selection could have acted preferentially towards variants located in coding and non-coding transcripts and characterized by a high deleteriousness score. Taken together, our analyses suggest that shared patterns of post admixture adaptation occurred at a continental scale in the Americas, affecting more often functional and impactful genomic variants.

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.

A Worldwide Map of Human Structural Variants.

Genomic variation extends from single nucleotide variants to large chromosomal rearrangements, but the extent of structural variation in Homo sapiens is still unclear. Almarri et al. provide a worldwide catalogue of structural variants present in human populations. Most of the reported variation is novel, with some variants being inherited from Neanderthals and Denisovans. Drift and selection shaped the distribution of these variants with some suggested to have functional implications.

The Genomic Impact of European Colonization of the Americas.

The human genetic diversity of the Americas has been affected by several events of gene flow that have continued since the colonial era and the Atlantic slave trade. Moreover, multiple waves of migration followed by local admixture occurred in the last two centuries, the impact of which has been largely unexplored. Here, we compiled a genome-wide dataset of ∼12,000 individuals from twelve American countries and ∼6,000 individuals from worldwide populations and applied haplotype-based methods to investigate how historical movements from outside the New World affected (1) the genetic structure, (2) the admixture profile, (3) the demographic history, and (4) sex-biased gene-flow dynamics of the Americas. We revealed a high degree of complexity underlying the genetic contribution of European and African populations in North and South America, from both geographic and temporal perspectives, identifying previously unreported sources related to Italy, the Middle East, and to specific regions of Africa.

Searching for archaic contribution in Africa.

Context: Africa's role in the narrative of human evolution is indisputably emphasised in the emergence of Homo sapiens. However, once humans dispersed beyond Africa, the history of those who stayed remains vastly under-studied, lacking the proper attention the birthplace of both modern and archaic humans deserves. The sequencing of Neanderthal and Denisovan genomes has elucidated evidence of admixture between archaic and modern humans outside of Africa, but has not aided efforts in answering whether archaic admixture happened within Africa. Objectives: This article reviews the state of research for archaic introgression in African populations and discusses recent insights into this topic. Methods: Gathering published sources and recently released preprints, this review reports on the different methods developed for detecting archaic introgression. Particularly it discusses how relevant these are when implemented on African populations and what findings these studies have shown so far. Results: Methods for detecting archaic introgression have been predominantly developed and implemented on non-African populations. Recent preprints present new methods considering African populations. While a number of studies using these methods suggest archaic introgression in Africa, without an African archaic genome to validate these results, such findings remain as putative archaic introgression. Conclusion: In light of the caveats with implementing current archaic introgression detection methods in Africa, we recommend future studies to concentrate on unravelling the complicated demographic history of Africa through means of ancient DNA where possible and through more focused efforts to sequence modern DNA from more representative populations across the African continent.

A missing piece of the Papio puzzle: Gorongosa baboon phenostructure and intrageneric relationships.

Most authors recognize six baboon species: hamadryas (Papio hamadryas), Guinea (Papio papio), olive (Papio anubis), yellow (Papio cynocephalus), chacma (Papio ursinus), and Kinda (Papio kindae). However, there is still debate regarding the taxonomic status, phylogenetic relationships, and the amount of gene flow occurring between species. Here, we present ongoing research on baboon morphological diversity in Gorongosa National Park (GNP), located in central Mozambique, south of the Zambezi River, at the southern end of the East African Rift System. The park exhibits outstanding ecological diversity and hosts more than 200 baboon troops. Gorongosa National Park baboons have previously been classified as chacma baboons (P. ursinus). In accordance with this, two mtDNA samples from the park have been placed in the same mtDNA clade as the northern chacma baboons. However, GNP baboons exhibit morphological features common in yellow baboons (e.g., yellow fur color), suggesting that parapatric gene flow between chacma and yellow baboons might have occurred in the past or could be ongoing. We investigated the phenostructure of the Gorongosa baboons using two approaches: 1) description of external phenotypic features, such as coloration and body size, and 2) 3D geometric morphometric analysis of 43 craniofacial landmarks on 11 specimens from Gorongosa compared to a pan-African sample of 352 baboons. The results show that Gorongosa baboons exhibit a mosaic of features shared with southern P. cynocephalus and P. ursinus griseipes. The GNP baboon phenotype fits within a geographic clinal pattern of replacing allotaxa. We put forward the hypothesis of either past and/or ongoing hybridization between the gray-footed chacma and southern yellow baboons in Gorongosa or an isolation-by-distance scenario in which the GNP baboons are geographically and morphologically intermediate. These two scenarios are not mutually exclusive. We highlight the potential of baboons as a useful model to understand speciation and hybridization in early human evolution.

The evolutionary history of Southern Africa.

The genomic variability of Southern African groups is characterized by an exceptional degree of diversity, which is the result of long-term local evolutionary history, migrations and gene-flow. Over the last few years several investigations have identified and described signatures related to these processes, revealing how ancient and more recent events have shaped the structure and ancestry composition of local populations. Here we discuss recent insights into the genetic history of the Southernmost part of the African continent provided by the analysis of modern and ancient genomes. Future work is expected to clarify the population dynamics associated with the emergence of Homo sapiens across Africa and the details of the process of dispersion and admixture associated with the arrival of Bantu-speaking groups in the region.

Stuck in fragments: Population genetics of the Endangered collared brown lemur Eulemur collaris in the Malagasy littoral forest.

OBJECTIVES: The Endangered collared brown lemur (Eulemur collaris) is the largest primate living in the littoral forest of southeastern Madagascar, a top priority habitat for biodiversity conservation on the island. Because this lemur is a key seed-disperser, an evaluation of the structure and connectivity of the populations surviving in the forest fragments is urgently needed to guide conservation plans. MATERIALS AND METHODS: Genetic variability at autosomal microsatellites and mitochondrial DNA was investigated in a total of 49 collared brown lemurs sampled by non-invasive methods in three littoral forest fragments and in the nearby lowland humid forest. RESULTS: The overall genetic diversity of E. collaris in the southeastern coastal region of Madagascar was lower than in other populations, as well as in other lemur species. The population appears highly structured, with less variable and more inbred groups inhabiting the littoral forest fragments compared to the inland area. Major barriers to gene flow were identified isolating littoral forest fragments from each other and from the inland lowland humid forest. DISCUSSION: Medium to long-term drift and scarce gene flow is the scenario that best explains the current genetic distribution. Habitat discontinuities such as rivers and grassland between forest fragments played a major role in structuring the population. A common history of size contraction is pointed out by several genetic estimators, indicating a possible ecological crisis triggered around 1,300 years ago. The adoption of strategies aimed at facilitating gene flow and population growth appears crucial to delay further loss of genetic diversity.

Complex Ancient Genetic Structure and Cultural Transitions in Southern African Populations.

The characterization of the structure of southern African populations has been the subject of numerous genetic, medical, linguistic, archaeological, and anthropological investigations. Current diversity in the subcontinent is the result of complex events of genetic admixture and cultural contact between early inhabitants and migrants that arrived in the region over the last 2000 years. Here, we analyze 1856 individuals from 91 populations, comprising novel and published genotype data, to characterize the genetic ancestry profiles of 631 individuals from 51 southern African populations. Combining both local ancestry and allele frequency based analyses, we identify a tripartite, ancient, Khoesan-related genetic structure. This structure correlates neither with linguistic affiliation nor subsistence strategy, but with geography, revealing the importance of isolation-by-distance dynamics in the area. Fine-mapping of these components in southern African populations reveals admixture and cultural reversion involving several Khoesan groups, and highlights that Bantu speakers and Coloured individuals have different mixtures of these ancient ancestries.

The Simons Genome Diversity Project: 300 genomes from 142 diverse populations.

Here we report the Simons Genome Diversity Project data set: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioural modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that of other non-Africans.

Low AMY1 Gene Copy Number Is Associated with Increased Body Mass Index in Prepubertal Boys.

BACKGROUND: Genome-wide association studies have identified more than 60 single nucleotide polymorphisms associated with Body Mass Index (BMI). Additional genetic variants, such as copy number variations (CNV), have also been investigated in relation to BMI. Recently, the highly polymorphic CNV in the salivary amylase (AMY1) gene, encoding an enzyme implicated in the first step of starch digestion, has been associated with obesity in adults and children. We assessed the potential association between AMY1 copy number and a wide range of BMI in a population of Italian school-children. METHODS: 744 children (354 boys, 390 girls, mean age (±SD): 8.4±1.4years) underwent anthropometric assessments (height, weight) and collection of saliva samples for DNA extraction. AMY1 copies were evaluated by quantitative PCR. RESULTS: A significant increase of BMI z-score by decreasing AMY1 copy number was observed in boys (ß: -0.117, p = 0.033), but not in girls. Similarly, waist circumference (ß: -0.155, p = 0.003, adjusted for age) was negatively influenced by AMY1 copy number in boys. Boys with 8 or more AMY1 copy numbers presented a significant lower BMI z-score (p = 0.04) and waist circumference (p = 0.01) when compared to boys with less than 8 copy numbers. CONCLUSIONS: In this pediatric-only, population-based study, a lower AMY1 copy number emerged to be associated with increased BMI in boys. These data confirm previous findings from adult studies and support a potential role of a higher copy number of the salivary AMY1 gene in protecting from excess weight gain.