Archaic Hominin Introgression in Africa Contributes to Functional Salivary MUC7 Genetic Variation

Emergence of saliva protein genes in the secretory calcium-binding phosphoprotein (SCPP) locus and accelerated evolution in primates

Genes within the secretory calcium-binding phosphoprotein (SCPP) family evolved in conjunction with major evolutionary milestones: the formation of a calcified skeleton in vertebrates, the emergence of tooth enamel in fish, and the introduction of lactation in mammals. The SCPP gene family also contains genes expressed primarily and abundantly in human saliva. Here, we explored the evolution of the saliva-related SCPP genes by harnessing currently available genomic and transcriptomic resources. Our findings provide insights into the expansion and diversification of SCPP genes, notably identifying previously undocumented convergent gene duplications. In primate genomes, we found additional duplication and diversification events that affected genes coding for proteins secreted in saliva. These saliva-related SCPP genes exhibit signatures of positive selection in the primate lineage while the other genes in the same locus remain conserved. We found that regulatory shifts and gene turnover events facilitated the accelerated gain of salivary expression. Collectively, our results position the SCPP gene family as a hotbed of evolutionary innovation, suggesting the potential role of dietary and pathogenic pressures in the adaptive diversification of the saliva composition in primates, including humans.

A Catalog of Coding Sequence Variations in Salivary Proteins' Genes Occurring during Recent Human Evolution

Saliva houses over 2000 proteins and peptides with poorly clarified functions, including proline-rich proteins, statherin, P-B peptides, histatins, cystatins, and amylases. Their genes are poorly conserved across related species, reflecting an evolutionary adaptation. We searched the nucleotide substitutions fixed in these salivary proteins' gene loci in modern humans compared with ancient hominins. We mapped 3472 sequence variants/nucleotide substitutions in coding, noncoding, and 5'-3' untranslated regions. Despite most of the detected variations being within noncoding regions, the frequency of coding variations was far higher than the general rate found throughout the genome. Among the various missense substitutions, specific substitutions detected in PRB1 and PRB2 genes were responsible for the introduction/abrogation of consensus sequences recognized by convertase enzymes that cleave the protein precursors. Overall, these changes that occurred during the recent human evolution might have generated novel functional features and/or different expression ratios among the various components of the salivary proteome. This may have influenced the homeostasis of the oral cavity environment, possibly conditioning the eating habits of modern humans. However, fixed nucleotide changes in modern humans represented only 7.3% of all the substitutions reported in this study, and no signs of evolutionary pressure or adaptative introgression from archaic hominins were found on the tested genes.

The MUC19 gene in Denisovans, Neanderthals, and Modern Humans: An Evolutionary History of Recurrent Introgression and Natural Selection

All humans carry a small fraction of archaic ancestry across the genome, the legacy of gene flow from Neanderthals, Denisovans, and other hominids into the ancestors of modern humans. While the effects of Neanderthal ancestry on human fitness and health have been explored more thoroughly, there are fewer examples of adaptive introgression of Denisovan variants. Here, we study the gene MUC19, for which some modern humans carry a Denisovan-like haplotype. MUC19 is a mucin, a glycoprotein that forms gels with various biological functions, from lubrication to immunity. We find the diagnostic variants for the Denisovan-like MUC19 haplotype at high frequencies in admixed Latin American individuals among global population, and at highest frequency in 23 ancient Indigenous American individuals, all predating population admixture with Europeans and Africans. We find that some Neanderthals--Vindija and Chagyrskaya--carry the Denisovan-like MUC19 haplotype, and that it was likely introgressed into human populations through Neanderthal introgression rather than Denisovan introgression. Finally, we find that the Denisovan-like MUC19 haplotype carries a higher copy number of a 30 base-pair variable number tandem repeat relative to the Human-like haplotype, and that copy numbers of this repeat are exceedingly high in American populations. Our results suggest that the Denisovan-like MUC19 haplotype served as the raw genetic material for positive selection as American populations adapted to novel environments during their movement from Beringia into North and then South America.

A Population Genetic Perspective on Subsistence Systems in the Sahel/Savannah Belt of Africa and the Historical Role of Pastoralism

This review focuses on the Sahel/Savannah belt, a large region of Africa where two alternative subsistence systems (pastoralism and agriculture), nowadays, interact. It is a long-standing question whether the pastoralists became isolated here from other populations after cattle began to spread into Africa (~8 thousand years ago, kya) or, rather, began to merge with other populations, such as agropastoralists, after the domestication of sorghum and pearl millet (~5 kya) and with the subsequent spread of agriculture. If we look at lactase persistence, a trait closely associated with pastoral lifestyle, we see that its variants in current pastoralists distinguish them from their farmer neighbours. Most other (mostly neutral) genetic polymorphisms do not, however, indicate such clear differentiation between these groups; they suggest a common origin and/or an extensive gene flow. Genetic affinity and ecological symbiosis between the two subsistence systems can help us better understand the population history of this African region. In this review, we show that genomic datasets of modern Sahel/Savannah belt populations properly collected in local populations can complement the still insufficient archaeological research of this region, especially when dealing with the prehistory of mobile populations with perishable material culture and therefore precarious archaeological visibility.

MaLAdapt Reveals Novel Targets of Adaptive Introgression From Neanderthals and Denisovans in Worldwide Human Populations

Adaptive introgression (AI) facilitates local adaptation in a wide range of species. Many state-of-the-art methods detect AI with ad-hoc approaches that identify summary statistic outliers or intersect scans for positive selection with scans for introgressed genomic regions. Although widely used, approaches intersecting outliers are vulnerable to a high false-negative rate as the power of different methods varies, especially for complex introgression events. Moreover, population genetic processes unrelated to AI, such as background selection or heterosis, may create similar genomic signals to AI, compromising the reliability of methods that rely on neutral null distributions. In recent years, machine learning (ML) methods have been increasingly applied to population genetic questions. Here, we present a ML-based method called MaLAdapt for identifying AI loci from genome-wide sequencing data. Using an Extra-Trees Classifier algorithm, our method combines information from a large number of biologically meaningful summary statistics to capture a powerful composite signature of AI across the genome. In contrast to existing methods, MaLAdapt is especially well-powered to detect AI with mild beneficial effects, including selection on standing archaic variation, and is robust to non-AI selective sweeps, heterosis from deleterious mutations, and demographic misspecification. Furthermore, MaLAdapt outperforms existing methods for detecting AI based on the analysis of simulated data and the validation of empirical signals through visual inspection of haplotype patterns. We apply MaLAdapt to the 1000 Genomes Project human genomic data and discover novel AI candidate regions in non-African populations, including genes that are enriched in functionally important biological pathways regulating metabolism and immune responses.

Including diverse and admixed populations in genetic epidemiology research

The inclusion of ancestrally diverse participants in genetic studies can lead to new discoveries and is important to ensure equitable health care benefit from research advances. Here, members of the Ethical, Legal, Social, Implications (ELSI) committee of the International Genetic Epidemiology Society (IGES) offer perspectives on methods and analysis tools for the conduct of inclusive genetic epidemiology research, with a focus on admixed and ancestrally diverse populations in support of reproducible research practices. We emphasize the importance of distinguishing socially defined population categorizations from genetic ancestry in the design, analysis, reporting, and interpretation of genetic epidemiology research findings. Finally, we discuss the current state of genomic resources used in genetic association studies, functional interpretation, and clinical and public health translation of genomic findings with respect to diverse populations.

Role of salivary glycopatterns for oral microbiota associated with gastric cancer

Microbiota in the oral cavity plays an important role in maintaining human health. Our previous studies have revealed significant alterations of salivary glycopatterns in gastric cancer (GC) patients, but it is unclear whether these altered salivary glycopatterns can cause the dysbiosis of oral microbiota. In this study, the oral microbiome of healthy volunteers (HVs) and GC patients were detected. The neoglycoproteins were then synthesized according to the altered glycopatterns in GC patients and used to explore the effects of specific salivary glycopattern against oral microbiota. The results showed that five species were significantly increased (p < 0.05) while two species were significantly decreased (p < 0.01) in the saliva of GC patients compared with that of HVs. And the fucose-neoglycoproteins (30-100 μg/mL) could reduce the adhesion and toxicity of Aggregatibacter segnis (A. segnis) to oral cells (HOEC and CAL-27), change the glycan structures of lipopolysaccharide on the surface of A. segnis, and enhance the capacity of A. segnis to trigger innate immune responses. This study revealed that the changes of salivary protein glycopatterns in GC patients might contribute to the dysbiosis of oral microbiota, and had important implications in developing new carbohydrate drugs to maintain a balanced microbiota in the oral.

The history and evolution of the Denisovan-EPAS1 haplotype in Tibetans

Recent studies suggest that admixture with archaic hominins played an important role in facilitating biological adaptations to new environments. For example, interbreeding with Denisovans facilitated the adaptation to high-altitude environments on the Tibetan Plateau. Specifically, the EPAS1 gene, a transcription factor that regulates the response to hypoxia, exhibits strong signatures of both positive selection and introgression from Denisovans in Tibetan individuals. Interestingly, despite being geographically closer to the Denisova Cave, East Asian populations do not harbor as much Denisovan ancestry as populations from Melanesia. Recently, two studies have suggested two independent waves of Denisovan admixture into East Asians, one of which is shared with South Asians and Oceanians. Here, we leverage data from EPAS1 in 78 Tibetan individuals to interrogate which of these two introgression events introduced the EPAS1 beneficial sequence into the ancestral population of Tibetans, and we use the distribution of introgressed segment lengths at this locus to infer the timing of the introgression and selection event. We find that the introgression event unique to East Asians most likely introduced the beneficial haplotype into the ancestral population of Tibetans around 48,700 (16,000-59,500) y ago, and selection started around 9,000 (2,500-42,000) y ago. Our estimates suggest that one of the most convincing examples of adaptive introgression is in fact selection acting on standing archaic variation.

African genetic diversity and adaptation inform a precision medicine agenda

The deep evolutionary history of African populations, since the emergence of modern humans more than 300,000 years ago, has resulted in high genetic diversity and considerable population structure. Selected genetic variants have increased in frequency due to environmental adaptation, but recent exposures to novel pathogens and changes in lifestyle render some of them with properties leading to present health liabilities. The unique discoverability potential from African genomic studies promises invaluable contributions to understanding the genomic and molecular basis of health and disease. Globally, African populations are understudied, and precision medicine approaches are largely based on data from European and Asian-ancestry populations, which limits the transferability of findings to the continent of Africa. Africa needs innovative precision medicine solutions based on African data that use knowledge and implementation strategies aligned to its climatic, cultural, economic and genomic diversity.

Functional Specialization of Human Salivary Glands and Origins of Proteins Intrinsic to Human Saliva

Salivary proteins are essential for maintaining health in the oral cavity and proximal digestive tract, and they serve as potential diagnostic markers for monitoring human health and disease. However, their precise organ origins remain unclear. Through transcriptomic analysis of major adult and fetal salivary glands and integration with the saliva proteome, the blood plasma proteome, and transcriptomes of 28+ organs, we link human saliva proteins to their source, identify salivary-gland-specific genes, and uncover fetal- and adult-specific gene repertoires. Our results also provide insights into the degree of gene retention during gland maturation and suggest that functional diversity among adult gland types is driven by specific dosage combinations of hundreds of transcriptional regulators rather than by a few gland-specific factors. Finally, we demonstrate the heterogeneity of the human acinar cell lineage. Our results pave the way for future investigations into glandular biology and pathology, as well as saliva's use as a diagnostic fluid.

The Impact of Ancient Genome Studies in Archaeology

The study of ancient genomes has burgeoned at an incredible rate in the last decade. The result is a shift in archaeological narratives, bringing with it a fierce debate on the place of genetics in anthropological research. Archaeogenomics has challenged and scrutinized fundamental themes of anthropological research, including human origins, movement of ancient and modern populations, the role of social organization in shaping material culture, and the relationship between culture, language, and ancestry. Moreover, the discussion has inevitably invoked new debates on indigenous rights, ownership of ancient materials, inclusion in the scientific process, and even the meaning of what it is to be a human. We argue that the broad and seemingly daunting ethical, methodological, and theoretical challenges posed by archaeogenomics, in fact, represent the very cutting edge of social science research. Here, we provide a general review of the field by introducing the contemporary discussion points and summarizing methodological and ethical concerns, while highlighting the exciting possibilities of ancient genome studies in archaeology from an anthropological perspective.

Human and Nonhuman Primate Lineage-Specific Footprints in the Salivary Proteome

Proteins in saliva are needed for preprocessing food in the mouth, maintenance of tooth mineralization, and protection from microbial pathogens. Novel insights into human lineage-specific functions of salivary proteins and clues to their involvement in human disease can be gained through evolutionary studies, as recently shown for salivary amylase AMY1 and salivary agglutinin DMBT1/gp340. However, the entirety of proteins in saliva, the salivary proteome, has not yet been investigated from an evolutionary perspective. Here, we compared the proteomes of human saliva and the saliva of our closest extant evolutionary relatives, chimpanzees and gorillas, using macaques as an outgroup, with the aim to uncover features in saliva protein composition that are unique to each species. We found that humans produce a waterier saliva, containing less than half total protein than great apes and Old World monkeys. For all major salivary proteins in humans, we could identify counterparts in chimpanzee and gorilla saliva. However, we discovered unique protein profiles in saliva of humans that were distinct from those of nonhuman primates. These findings open up the possibility that dietary differences and pathogenic pressures may have shaped a distinct salivary proteome in the human lineage.

VolcanoFinder: Genomic scans for adaptive introgression

Recent research shows that introgression between closely-related species is an important source of adaptive alleles for a wide range of taxa. Typically, detection of adaptive introgression from genomic data relies on comparative analyses that require sequence data from both the recipient and the donor species. However, in many cases, the donor is unknown or the data is not currently available. Here, we introduce a genome-scan method—VolcanoFinder—to detect recent events of adaptive introgression using polymorphism data from the recipient species only. VolcanoFinder detects adaptive introgression sweeps from the pattern of excess intermediate-frequency polymorphism they produce in the flanking region of the genome, a pattern which appears as a volcano-shape in pairwise genetic diversity. Using coalescent theory, we derive analytical predictions for these patterns. Based on these results, we develop a composite-likelihood test to detect signatures of adaptive introgression relative to the genomic background. Simulation results show that VolcanoFinder has high statistical power to detect these signatures, even for older sweeps and for soft sweeps initiated by multiple migrant haplotypes. Finally, we implement VolcanoFinder to detect archaic introgression in European and sub-Saharan African human populations, and uncovered interesting candidates in both populations, such as TSHR in Europeans and TCHH-RPTN in Africans. We discuss their biological implications and provide guidelines for identifying and circumventing artifactual signals during empirical applications of VolcanoFinder.

The process by which beneficial alleles are introduced into a species from a closely-related species is termed adaptive introgression. We present an analytically-tractable model for the effects of adaptive introgression on non-adaptive genetic variation in the genomic region surrounding the beneficial allele. The result we describe is a characteristic volcano-shaped pattern of increased variability that arises around the positively-selected site, and we introduce an open-source method VolcanoFinder to detect this signal in genomic data. Importantly, VolcanoFinder is a population-genetic likelihood-based approach, rather than a comparative-genomic approach, and can therefore probe genomic variation data from a single population for footprints of adaptive introgression, even from a priori unknown and possibly extinct donor species.

Analysis of Haplotypic Variation and Deletion Polymorphisms Point to Multiple Archaic Introgression Events, Including from Altai Neanderthal Lineage

The time, extent, and genomic effect of the introgressions from archaic humans into ancestors of extant human populations remain some of the most exciting venues of population genetics research in the past decade. Several studies have shown population-specific signatures of introgression events from Neanderthals, Denisovans, and potentially other unknown hominin populations in different human groups. Moreover, it was shown that these introgression events may have contributed to phenotypic variation in extant humans, with biomedical and evolutionary consequences. In this study, we present a comprehensive analysis of the unusually divergent haplotypes in the Eurasian genomes and show that they can be traced back to multiple introgression events. In parallel, we document hundreds of deletion polymorphisms shared with Neanderthals. A locus-specific analysis of one such shared deletion suggests the existence of a direct introgression event from the Altai Neanderthal lineage into the ancestors of extant East Asian populations. Overall, our study is in agreement with the emergent notion that various Neanderthal populations contributed to extant human genetic variation in a population-specific manner.

Recovering signals of ghost archaic introgression in African populations

While introgression from Neanderthals and Denisovans has been documented in modern humans outside Africa, the contribution of archaic hominins to the genetic variation of present-day Africans remains poorly understood. We provide complementary lines of evidence for archaic introgression into four West African populations. Our analyses of site frequency spectra indicate that these populations derive 2 to 19% of their genetic ancestry from an archaic population that diverged before the split of Neanderthals and modern humans. Using a method that can identify segments of archaic ancestry without the need for reference archaic genomes, we built genome-wide maps of archaic ancestry in the Yoruba and the Mende populations. Analyses of these maps reveal segments of archaic ancestry at high frequency in these populations that represent potential targets of adaptive introgression. Our results reveal the substantial contribution of archaic ancestry in shaping the gene pool of present-day West African populations.

Archaic hominin introgression into modern human genomes

Ancient genomes from multiple Neanderthal and the Denisovan individuals, along with DNA sequence data from diverse contemporary human populations strongly support the prevalence of gene flow among different hominins. Recent studies now provide evidence for multiple gene flow events that leave genetic signatures in extant and ancient human populations. These events include older gene flow from an unknown hominin in Africa predating out-of-Africa migrations, and in the last 50,000-100,000 years, multiple gene flow events from Neanderthals into ancestral Eurasian human populations, and at least three distinct introgression events from a lineage close to Denisovans into ancestors of extant Southeast Asian and Oceanic populations. Some of these introgression events may have happened as late as 20,000 years before present and reshaped the way in which we think about human evolution. In this review, I aim to answer anthropologically relevant questions with regard to recent research on ancient hominin introgression in the human lineage. How have genomic data from archaic hominins changed our view of human evolution? Is there any doubt about whether introgression from ancient hominins to the ancestors of present-day humans occurred? What is the current view of human evolutionary history from the genomics perspective? What is the impact of introgression on human phenotypes?

Human saliva and model saliva at bulk to adsorbed phases - similarities and differences

Human saliva, a seemingly simple aqueous fluid, is, in fact, an extraordinarily complex biocolloid that is not fully understood, despite many decades of study. Salivary lubrication is widely believed to be a signature of good oral health and is also crucial for speech, food oral processing and swallowing. However, saliva has been often neglected in food colloid research, primarily due to its high intra- to inter-individual variability and altering material properties upon collection and storage, when used as an ex vivo research material. In the last few decades, colloid scientists have attempted designing model (i.e. 'saliva mimicking fluid') salivary formulations to understand saliva-food colloid interactions in an in vitro set up and its contribution on microstructural aspects, lubrication properties and sensory perception. In this Review, we critically examine the current state of knowledge on bulk and interfacial properties of model saliva in comparison to real human saliva and highlight how far such model salivary formulations can match the properties of real human saliva. Many, if not most, of these model saliva formulations share similarities with real human saliva in terms of biochemical compositions, including electrolytes, pH and concentrations of salivary proteins, such as α-amylase and highly glycosylated mucins. This, together with similarities between model and real saliva in terms of surface charge, has led to significant advancement in decoding various colloidal interactions (bridging, depletion) of charged emulsion droplets and associated sensory perception in the oral phase. However, model saliva represents significant dissimilarity to real saliva in terms of lubricating properties. Based on in-depth examination of properties of mucins derived from animal sources (e.g. pig gastric mucins (PGM) or bovine submaxillary mucin (BSM)), we can recommend that BSM is currently the most optimal commercially available mucin source when attempting to replicate saliva based on surface adsorption and lubrication properties. Even though purification via dialysis or chromatographic techniques may influence various physicochemical properties of BSM, such as structure and surface adsorption, the lubricating properties of model saliva formulations based on BSM are generally superior and more reliable than the PGM counterpart at orally relevant pH. Comparison of mucin-containing model saliva with ex vivo human salivary conditioning films suggests that mucin alone cannot replicate the lubricity of real human salivary pellicle. Mucin-based multi-layers containing mucin and oppositely charged polyelectrolytes may offer promising avenues in the future for engineering biomimetic salivary pellicle, however, this has not been explored in oral tribology experiments to date. Hence, there is a strong need for systematic studies with employment of model saliva formulations containing mucins with and without polycationic additives before a consensus on a standardized model salivary formulation can be achieved. Overall, this review provides the first comprehensive framework on simulating saliva for a particular bulk or surface property when doing food oral processing experiments.

An Evolutionary Perspective on the Impact of Genomic Copy Number Variation on Human Health

Copy number variants (CNVs), deletions and duplications of segments of DNA, account for at least five times more variable base pairs in humans than single-nucleotide variants. Several common CNVs were shown to change coding and regulatory sequences and thus dramatically affect adaptive phenotypes involving immunity, perception, metabolism, skin structure, among others. Some of these CNVs were also associated with susceptibility to cancer, infection, and metabolic disorders. These observations raise the possibility that CNVs are a primary contributor to human phenotypic variation and consequently evolve under selective pressures. Indeed, locus-specific haplotype-level analyses revealed signatures of natural selection on several CNVs. However, more traditional tests of selection which are often applied to single-nucleotide variation often have diminished statistical power when applied to CNVs because they often do not show strong linkage disequilibrium with nearby variants. Recombination-based formation mechanisms of CNVs lead to frequent recurrence and gene conversion events, breaking the linkage disequilibrium involving CNVs. Similar methodological challenges also prevent routine genome-wide association studies to adequately investigate the impact of CNVs on heritable human disease. Thus, we argue that the full relevance of CNVs to human health and evolution is yet to be elucidated. We further argue that a holistic investigation of formation mechanisms within an evolutionary framework would provide a powerful framework to understand the functional and biomedical impact of CNVs. In this paper, we review several cases where studies reveal diverse evolutionary histories and unexpected functional consequences of CNVs. We hope that this review will encourage further work on CNVs by both evolutionary and medical geneticists.

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.

Genomics of disease risk in globally diverse populations

Risk of disease is multifactorial and can be shaped by socio-economic, demographic, cultural, environmental and genetic factors. Our understanding of the genetic determinants of disease risk has greatly advanced with the advent of genome-wide association studies (GWAS), which detect associations between genetic variants and complex traits or diseases by comparing populations of cases and controls. However, much of this discovery has occurred through GWAS of individuals of European ancestry, with limited representation of other populations, including from Africa, The Americas, Asia and Oceania. Population demography, genetic drift and adaptation to environments over thousands of years have led globally to the diversification of populations. This global genomic diversity can provide new opportunities for discovery and translation into therapies, as well as a better understanding of population disease risk. Large-scale multi-ethnic and representative biobanks and population health resources provide unprecedented opportunities to understand the genetic determinants of disease on a global scale.

Ancient admixture from an extinct ape lineage into bonobos

Admixture is a recurrent phenomenon in humans and other great ape populations. Genetic information from extinct hominins allows us to study historical interactions with modern humans and discover adaptive functions of gene flow. Here, we investigate whole genomes from bonobo and chimpanzee populations for signatures of gene flow from unknown archaic populations, finding evidence for an ancient admixture event between bonobos and a divergent lineage. This result reveals a complex population history in our closest living relatives, probably several hundred thousand years ago. We reconstruct up to 4.8% of the genome of this 'ghost' ape, which represents genomic data of an extinct great ape population. Genes contained in archaic fragments might confer functional consequences for the immunity, behaviour and physiology of bonobos. Finally, comparing the landscapes of introgressed regions in humans and bonobos, we find that a recurrent depletion of introgression is rare, suggesting that genomic incompatibilities arose seldom in these lineages.

Whole-genome sequence analysis of a Pan African set of samples reveals archaic gene flow from an extinct basal population of modern humans into sub-Saharan populations

Background

Population demography and gene flow among African groups, as well as the putative archaic introgression of ancient hominins, have been poorly explored at the genome level.

Results

Here, we examine 15 African populations covering all major continental linguistic groups, ecosystems, and lifestyles within Africa through analysis of whole-genome sequence data of 21 individuals sequenced at deep coverage. We observe a remarkable correlation among genetic diversity and geographic distance, with the hunter-gatherer groups being more genetically differentiated and having larger effective population sizes throughout most modern-human history. Admixture signals are found between neighbor populations from both hunter-gatherer and agriculturalists groups, whereas North African individuals are closely related to Eurasian populations. Regarding archaic gene flow, we test six complex demographic models that consider recent admixture as well as archaic introgression. We identify the fingerprint of an archaic introgression event in the sub-Saharan populations included in the models (~ 4.0% in Khoisan, ~ 4.3% in Mbuti Pygmies, and ~ 5.8% in Mandenka) from an early divergent and currently extinct ghost modern human lineage.

Conclusion

The present study represents an in-depth genomic analysis of a Pan African set of individuals, which emphasizes their complex relationships and demographic history at population level.

Electronic supplementary material

The online version of this article (10.1186/s13059-019-1684-5) contains supplementary material, which is available to authorized users.

Identification of a nine-gene panel as a prognostic indicator for recurrence with muscle-invasive bladder cancer

BACKGROUND AND OBJECTIVES

Bladder cancer is one of the most common and highly recurrent cancers worldwide. Recurrence-associated genes may potentially predict cancer recurrence. We aimed to construct a recurrence-associated gene panel to improve the prognostic prediction of bladder cancer.

METHODS

Based on DNA sequencing and clinical data from the TCGA-BLCA project, we identified 10 potential driver genes significantly associated with recurrence of bladder cancer. We performed multivariable logistic regression analysis to construct an optimized recurrence prediction model with nine recurrence-associated genes (EME1, AKAP9, ZNF91, PARD3, STAG2, ZFP36L2, METTL3, POLR3B, and MUC7) and clinical information as the independent variables.

RESULTS

The area under the receiver operating characteristic (ROC) curve was 0.80 in this model, much higher than that of the baseline model (AUC = 0.73) and the same trend was also validated in its subset. Decision curve analysis also revealed that there is a significant net benefit gained by adding nine genes mutation to the baseline model. Furthermore, Kaplan-Meier survival analysis showed that eight out of the nine genes (excluding MUC7) had good effects on the overall prognosis of patients.

CONCLUSIONS

This nine-gene panel will most likely be a useful tool for prognostic evaluation and will facilitate the personalized management of patients with bladder cancer.

High-resolution quantitative proteomics applied to the study of the specific protein signature in the sputum and saliva of active tuberculosis patients and their infected and uninfected contacts

Our goal was to establish panels of protein biomarkers that are characteristic of patients with microbiologically confirmed pulmonary tuberculosis (TB) and their contacts, including latent TB-infected (LTBI) and uninfected patients. Since the first pathogen-host contact occurs in the oral and nasal passages the saliva and sputum were chosen as the biological fluids to be studied. Quantitative shotgun proteomics was performed using a LTQ-Orbitrap-Elite platform. For active TB patients, both fluids exhibited a specific accumulation of proteins that were related to complement activation, inflammation and modulation of immune response. In the saliva of TB patients, a decrease of in proteins related to glucose and lipid metabolism was detected. In contrast, the sputum of uninfected contacts presented a specific proteomic signature that was composed of proteins involved in the perception of bitter taste, defense against pathogens and innate immune response, suggesting that those are key events during the initial entry of the pathogen in the host. SIGNIFICANCE: This is the first study to compare the saliva and sputum from active TB patients and their contacts. Our findings strongly suggest that TB patients show not only an activation of processes that are related to complement activation and modulation of inflammation but also an imbalance in carbohydrate and lipid metabolism. In addition, those individuals who do not get infected after direct exposure to the pathogen display a typical proteomic signature in the sputum, which is a reflection of the secretion from the nasal and oral mucosa, the first immunological barriers that M. tuberculosis encounters in the host. Thus, this result indicates the importance of the processes related to the innate immune response in fighting the initial events of the infection.

Variation and Functional Impact of Neanderthal Ancestry in Western Asia

Neanderthals contributed genetic material to modern humans via multiple admixture events. Initial admixture events presumably occurred in Western Asia shortly after humans migrated out of Africa. Despite being a focal point of admixture, earlier studies indicate lower Neanderthal introgression rates in some Western Asian populations as compared with other Eurasian populations. To better understand the genome-wide and phenotypic impact of Neanderthal introgression in the region, we sequenced whole genomes of nine present-day Europeans, Africans, and the Western Asian Druze at high depth, and analyzed available whole genome data from various other populations, including 16 genomes from present-day Turkey. Our results confirmed previous observations that contemporary Western Asian populations, on an average, have lower levels of Neanderthal-introgressed DNA relative to other Eurasian populations. Modern Western Asians also show comparatively high variability in Neanderthal ancestry, which may be attributed to the complex demographic history of the region. We further replicated the previously described depletion of putatively functional sequences among Neanderthal-introgressed haplotypes. Still, we find dozens of common Neanderthal-introgressed haplotypes in the Turkish sample associated with human phenotypes, including anthropometric and metabolic traits, as well as the immune response. One of these haplotypes is unusually long and harbors variants that affect the expression of members of the CCR gene family and are associated with celiac disease. Overall, our results paint a complex first picture of the genomic impact of Neanderthal introgression in the Western Asian populations.

The Malaria-Protective Human Glycophorin Structural Variant DUP4 Shows Somatic Mosaicism and Association with Hemoglobin Levels

Glycophorin A and glycophorin B are red blood cell surface proteins and are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A-glycophorin B fusion gene and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40%. Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant and reveal somatic variation in copy number of the glycophorin B-glycophorin A fusion gene. By developing a simple, specific, PCR-based assay for DUP4, we show that the DUP4 variant reaches a frequency of 13% in the population of a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation and show that DUP4 is associated with a malarial-related phenotype in a longitudinally followed population.

Complex Haplotypes of GSTM1 Gene Deletions Harbor Signatures of a Selective Sweep in East Asian Populations

The deletion of the metabolizing Glutathione S-transferase Mu 1 (GSTM1) gene has been associated with multiple cancers, metabolic and autoimmune disorders, as well as drug response. It is unusually common, with allele frequency reaching up to 75% in some human populations. Such high allele frequency of a derived allele with apparent impact on an otherwise conserved gene is a rare phenomenon. To investigate the evolutionary history of this locus, we analyzed 310 genomes using population genetics tools. Our analysis revealed a surprising lack of linkage disequilibrium between the deletion and the flanking single nucleotide variants in this locus. Tests that measure extended homozygosity and rapid change in allele frequency revealed signatures of an incomplete sweep in the locus. Using empirical approaches, we identified the Tanuki haplogroup, which carries the GSTM1 deletion and is found in approximately 70% of East Asian chromosomes. This haplogroup has rapidly increased in frequency in East Asian populations, contributing to a high population differentiation among continental human groups. We showed that extended homozygosity and population differentiation for this haplogroup is incompatible with simulated neutral expectations in East Asian populations. In parallel, we found that the Tanuki haplogroup is significantly associated with the expression levels of other GSTM genes. Collectively, our results suggest that standing variation in this locus has likely undergone an incomplete sweep in East Asia with regulatory impact on multiple GSTM genes. Our study provides the necessary framework for further studies to elucidate the evolutionary reasons that maintain disease-susceptibility variants in the GSTM1 locus.

Glycan recognition at the saliva - oral microbiome interface

The mouth is a first critical interface where most potentially harmful substances or pathogens contact the host environment. Adaptive and innate immune defense mechanisms are established there to inactivate or eliminate pathogenic microbes that traverse the oral environment on the way to their target organs and tissues. Protein and glycoprotein components of saliva play a particularly important role in modulating the oral microbiota and helping with the clearance of pathogens. It has long been acknowledged that glycobiological and glycoimmunological aspects play a pivotal role in oral host-microbe, microbe-host, and microbe-microbe interactions in the mouth. In this review, we aim to delineate how glycan-mediated host defense mechanisms in the oral cavity support human health. We will describe the role of glycans attached to large molecular size salivary glycoproteins which act as a first line of primordial host defense in the human mouth. We will further discuss how glycan recognition contributes to both colonization and clearance of oral microbes.

Evolutionary and Medical Consequences of Archaic Introgression into Modern Human Genomes

The demographic history of anatomically modern humans (AMH) involves multiple migration events, population extinctions and genetic adaptations. As genome-wide data from complete genome sequencing becomes increasingly abundant and available even from extinct hominins, new insights of the evolutionary history of our species are discovered. It is currently known that AMH interbred with archaic hominins once they left the African continent. Current non-African human genomes carry fragments of archaic origin. This review focuses on the fitness consequences of archaic interbreeding in current human populations. We discuss new insights and challenges that researchers face when interpreting the potential impact of introgression on fitness and testing hypotheses about the role of selection within the context of health and disease.

Did Our Species Evolve in Subdivided Populations across Africa, and Why Does It Matter?

We challenge the view that our species, Homo sapiens, evolved within a single population and/or region of Africa. The chronology and physical diversity of Pleistocene human fossils suggest that morphologically varied populations pertaining to the H. sapiens clade lived throughout Africa. Similarly, the African archaeological record demonstrates the polycentric origin and persistence of regionally distinct Pleistocene material culture in a variety of paleoecological settings. Genetic studies also indicate that present-day population structure within Africa extends to deep times, paralleling a paleoenvironmental record of shifting and fractured habitable zones. We argue that these fields support an emerging view of a highly structured African prehistory that should be considered in human evolutionary inferences, prompting new interpretations, questions, and interdisciplinary research directions.

The view that Homo sapiens evolved from a single region/population within Africa has been given primacy in studies of human evolution.

However, developments across multiple fields show that relevant data are no longer consistent with this view.

We argue instead that Homo sapiens evolved within a set of interlinked groups living across Africa, whose connectivity changed through time.

Genetic models therefore need to incorporate a more complex view of ancient migration and divergence in Africa.

We summarize this new framework emphasizing population structure, outline how this changes our understanding of human evolution, and identify new research directions.

Something old, something borrowed: admixture and adaptation in human evolution

The sequencing of ancient DNA from archaic humans-Neanderthals and Denisovans-has revealed that modern and archaic humans interbred at least twice during the Pleistocene. The field of human paleogenomics has now turned its attention towards understanding the nature of this genetic legacy in the gene pool of present-day humans. What exactly did modern humans obtain from interbreeding with Neanderthals and Denisovans? Was the introgressed genetic material beneficial, neutral or maladaptive? Can differences in phenotypes among present-day human populations be explained by archaic human introgression? These questions are of prime importance for our understanding of recent human evolution, but will require careful computational modeling and extensive functional assays before they can be answered in full. Here, we review the recent literature characterizing introgressed DNA and the likely biological consequences for their modern human carriers. We focus particularly on archaic human haplotypes that were beneficial to modern humans as they expanded across the globe, and on ways to understand how populations harboring these haplotypes evolved over time.

Outstanding questions in the study of archaic hominin admixture

The complete sequencing of archaic and modern human genomes has revolutionized the study of human history and evolution. The application of paleogenomics has answered questions that were beyond the scope of archaeology alone-definitively proving admixture between archaic and modern humans. Despite the remarkable progress made in the study of archaic-modern human admixture, many outstanding questions remain. Here, we review some of these questions, which include how frequent archaic-modern human admixture was in history, to what degree drift and selection are responsible for the loss and retention of introgressed sequences in modern human genomes, and how surviving archaic sequences affect human phenotypes.

Recent advances in understanding evolution of the placenta: insights from transcriptomics

The mammalian placenta shows an extraordinary degree of variation in gross and fine structure, but this has been difficult to interpret in physiological terms. Transcriptomics offers a path to understanding how structure relates to function. This essay examines how studies of gene transcription can inform us about placental evolution in eutherian and marsupial mammals and more broadly about convergent evolution of viviparity and placentation in vertebrates. Thus far, the focus has been on the chorioallantoic placenta of eutherians at term, the reproductive strategies of eutherians and marsupials, and the decidual response of the uterus at implantation. Future work should address gene expression during early stages of placental development and endeavor to cover all major groups of mammals. Comparative studies across oviparous and viviparous vertebrates have centered on the chorioallantoic membrane and yolk sac. They point to the possibility of defining a set of genes that can be recruited to support commonalities in reproductive strategies. Further advances can be anticipated from single-cell transcriptomics if those techniques are applied to a range of placental structures and in species other than humans and mice.

VCFtoTree: a user-friendly tool to construct locus-specific alignments and phylogenies from thousands of anthropologically relevant genome sequences

BACKGROUND

Constructing alignments and phylogenies for a given locus from large genome sequencing studies with relevant outgroups allow novel evolutionary and anthropological insights. However, no user-friendly tool has been developed to integrate thousands of recently available and anthropologically relevant genome sequences to construct complete sequence alignments and phylogenies.

RESULTS

Here, we provide VCFtoTree, a user friendly tool with a graphical user interface that directly accesses online databases to download, parse and analyze genome variation data for regions of interest. Our pipeline combines popular sequence datasets and tree building algorithms with custom data parsing to generate accurate alignments and phylogenies using all the individuals from the 1000 Genomes Project, Neanderthal and Denisovan genomes, as well as reference genomes of Chimpanzee and Rhesus Macaque. It can also be applied to other phased human genomes, as well as genomes from other species. The output of our pipeline includes an alignment in FASTA format and a tree file in newick format.

CONCLUSION

VCFtoTree fulfills the increasing demand for constructing alignments and phylogenies for a given loci from thousands of available genomes. Our software provides a user friendly interface for a wider audience without prerequisite knowledge in programming. VCFtoTree can be accessed from https://github.com/duoduoo/VCFtoTree_3.0.0 .