Tracing the peopling of the world through genomics

BaTwa populations from Zambia retain ancestry of past hunter-gatherer groups

Sub-equatorial Africa is today inhabited predominantly by Bantu-speaking groups of Western African descent who brought agriculture to the Luangwa valley in eastern Zambia ~2000 years ago. Before their arrival the area was inhabited by hunter-gatherers, who in many cases were subsequently replaced, displaced or assimilated. In Zambia, we know little about the genetic affinities of these hunter-gatherers. We examine ancestry of two isolated communities in Zambia, known as BaTwa and possible descendants of recent hunter-gatherers. We genotype over two million genome-wide SNPs from two BaTwa populations (total of 80 individuals) and from three comparative farming populations to: (i) determine if the BaTwa carry genetic links to past hunter-gatherer-groups, and (ii) characterise the genetic affinities of past Zambian hunter-gatherer-groups. The BaTwa populations do harbour a hunter-gatherer-like genetic ancestry and Western African ancestry. The hunter-gatherer component is a unique local signature, intermediate between current-day Khoe-San ancestry from southern Africa and central African rainforest hunter-gatherer ancestry.

Beyond the Human Genome Project: The Age of Complete Human Genome Sequences and Pangenome References

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

Recurrent gene flow between Neanderthals and modern humans over the past 200,000 years

Although it is well known that the ancestors of modern humans and Neanderthals admixed, the effects of gene flow on the Neanderthal genome are not well understood. We develop methods to estimate the amount of human-introgressed sequences in Neanderthals and apply it to whole-genome sequence data from 2000 modern humans and three Neanderthals. We estimate that Neanderthals have 2.5 to 3.7% human ancestry, and we leverage human-introgressed sequences in Neanderthals to revise estimates of Neanderthal ancestry in modern humans, show that Neanderthal population sizes were significantly smaller than previously estimated, and identify two distinct waves of modern human gene flow into Neanderthals. Our data provide insights into the genetic legacy of recurrent gene flow between modern humans and Neanderthals.

Biases in ARG-Based Inference of Historical Population Size in Populations Experiencing Selection

Inferring the demographic history of populations provides fundamental insights into species dynamics and is essential for developing a null model to accurately study selective processes. However, background selection and selective sweeps can produce genomic signatures at linked sites that mimic or mask signals associated with historical population size change. While the theoretical biases introduced by the linked effects of selection have been well established, it is unclear whether ancestral recombination graph (ARG)-based approaches to demographic inference in typical empirical analyses are susceptible to misinference due to these effects. To address this, we developed highly realistic forward simulations of human and Drosophila melanogaster populations, including empirically estimated variability of gene density, mutation rates, recombination rates, purifying, and positive selection, across different historical demographic scenarios, to broadly assess the impact of selection on demographic inference using a genealogy-based approach. Our results indicate that the linked effects of selection minimally impact demographic inference for human populations, although it could cause misinference in populations with similar genome architecture and population parameters experiencing more frequent recurrent sweeps. We found that accurate demographic inference of D. melanogaster populations by ARG-based methods is compromised by the presence of pervasive background selection alone, leading to spurious inferences of recent population expansion, which may be further worsened by recurrent sweeps, depending on the proportion and strength of beneficial mutations. Caution and additional testing with species-specific simulations are needed when inferring population history with non-human populations using ARG-based approaches to avoid misinference due to the linked effects of selection.

Widespread horse-based mobility arose around 2200 BCE in Eurasia

Horses revolutionized human history with fast mobility1. However, the timeline between their domestication and their widespread integration as a means of transport remains contentious2-4. Here we assemble a collection of 475 ancient horse genomes to assess the period when these animals were first reshaped by human agency in Eurasia. We find that reproductive control of the modern domestic lineage emerged around 2200 BCE, through close-kin mating and shortened generation times. Reproductive control emerged following a severe domestication bottleneck starting no earlier than approximately 2700 BCE, and coincided with a sudden expansion across Eurasia that ultimately resulted in the replacement of nearly every local horse lineage. This expansion marked the rise of widespread horse-based mobility in human history, which refutes the commonly held narrative of large horse herds accompanying the massive migration of steppe peoples across Europe around 3000 BCE and earlier3,5. Finally, we detect significantly shortened generation times at Botai around 3500 BCE, a settlement from central Asia associated with corrals and a subsistence economy centred on horses6,7. This supports local horse husbandry before the rise of modern domestic bloodlines.

Increased loss-of-function filaggrin gene mutation prevalence in atopic dermatitis patients across northern latitudes indicates genetic fitness: A systematic review and meta-analysis

Loss-of-function (LoF) mutations in the filaggrin gene (FLG) constitute the strongest genetic risk for atopic dermatitis (AD). A latitude-dependent difference in the prevalence of LoF FLG mutations was systematically evaluated. A systematic review and meta-analysis were performed to estimate the prevalence of LoF FLG mutations in AD patients and the general population by geography and ethnicity. Risk of bias was assessed by Newcastle-Ottawa Scale and Jadad score. StatsDirect, version 3 software was used to calculate all outcomes. PubMed and EMBASE were searched until 9th December 2021. Studies were included if they contained data on the prevalence of LoF FLG mutations in AD patients or from the general population or associations between AD and LoF FLG mutations and were authored in English. Overall, 248 studies and 229 310 AD patients and individuals of the general population were included in the quantitative analysis. The prevalence of LoF FLG mutations was 19.1% (95% CI, 17.3-21.0) in AD patients and 5.8% (95% CI, 5.3-6.2) in the general population. There was a significant positive association between AD and LoF FLG mutations in all latitudes in the Northern hemisphere, but not in all ethnicities. The prevalence of LoF FLG mutations became gradually more prevalent in populations residing farther north of the Equator but was negligible in Middle Easterners and absent in most African populations. FLG LoF mutations are common and tend to increase with northern latitude, suggesting potential clinical implications for future AD management. The existence of possible genetic fitness from FLG LoF mutations remains unknown.

Environmental conditions associated with initial northern expansion of anatomically modern humans

The ability of our ancestors to switch food sources and to migrate to more favourable environments enabled the rapid global expansion of anatomically modern humans beyond Africa as early as 120,000 years ago. Whether this versatility was largely the result of environmentally determined processes or was instead dominated by cultural drivers, social structures, and interactions among different groups, is unclear. We develop a statistical approach that combines both archaeological and genetic data to infer the more-likely initial expansion routes in northern Eurasia and the Americas. We then quantify the main differences in past environmental conditions between the more-likely routes and other potential (less-likely) routes of expansion. We establish that, even though cultural drivers remain plausible at finer scales, the emergent migration corridors were predominantly constrained by a combination of regional environmental conditions, including the presence of a forest-grassland ecotone, changes in temperature and precipitation, and proximity to rivers.

Shaping immunity: The influence of natural selection on population immune diversity

Humans exhibit considerable variability in their immune responses to the same immune challenges. Such variation is widespread and affects individual and population-level susceptibility to infectious diseases and immune disorders. Although the factors influencing immune response diversity are partially understood, what mechanisms lead to the wide range of immune traits in healthy individuals remain largely unexplained. Here, we discuss the role that natural selection has played in driving phenotypic differences in immune responses across populations and present-day susceptibility to immune-related disorders. Further, we touch on future directions in the field of immunogenomics, highlighting the value of expanding this work to human populations globally, the utility of modeling the immune response as a dynamic process, and the importance of considering the potential polygenic nature of natural selection. Identifying loci acted upon by evolution may further pinpoint variants critically involved in disease etiology, and designing studies to capture these effects will enrich our understanding of the genetic contributions to immunity and immune dysregulation.

Ghosts of extinct apes: genomic insights into African hominid evolution

We are accustomed to regular announcements of new hominin fossils. There are now some 6000 hominin fossils, and up to 31 species. However, where are the announcements of African ape fossils? The answer is that there are almost none. Our knowledge of African ape evolution is based entirely on genomic analyses, which show that extant diversity is very young. This contrasts with the extensive and deep diversity of hominins known from fossils. Does this difference point to low and late diversification of ape lineages, or high rates of extinction? The comparative evolutionary dynamics of African hominids are central to interpreting living ape adaptations, as well as understanding the patterns of hominin evolution and the nature of the last common ancestor.

Physiography, foraging mobility, and the first peopling of Sahul

The route and speed of migration into Sahul by Homo sapiens remain a major research question in archaeology. Here, we introduce an approach which models the impact of the physical environment on human mobility by combining time-evolving landscapes with Lévy walk foraging patterns, this latter accounting for a combination of short-distance steps and occasional longer moves that hunter-gatherers likely utilised for efficient exploration of new environments. Our results suggest a wave of dispersal radiating across Sahul following riverine corridors and coastlines. Estimated migration speeds, based on archaeological sites and predicted travelled distances, fall within previously reported range from Sahul and other regions. From our mechanistic movement simulations, we then analyse the likelihood of archaeological sites and highlight areas in Australia that hold archaeological potential. Our approach complements existing methods and provides interesting perspectives on the Pleistocene archaeology of Sahul that could be applied to other regions around the world.

Tumor phylogeography reveals block-shaped spatial heterogeneity and the mode of evolution in Hepatocellular Carcinoma

Solid tumors are complex ecosystems with heterogeneous 3D structures, but the spatial intra-tumor heterogeneity (sITH) at the macroscopic (i.e., whole tumor) level is under-explored. Using a phylogeographic approach, we sequence genomes and transcriptomes from 235 spatially informed sectors across 13 hepatocellular carcinomas (HCC), generating one of the largest datasets for studying sITH. We find that tumor heterogeneity in HCC segregates into spatially variegated blocks with large genotypic and phenotypic differences. By dissecting the transcriptomic heterogeneity, we discover that 30% of patients had a "spatially competing distribution" (SCD), where different spatial blocks have distinct transcriptomic subtypes co-existing within a tumor, capturing the critical transition period in disease progression. Interestingly, the tumor regions with more advanced transcriptomic subtypes (e.g., higher cell cycle) often take clonal dominance with a wider geographic range, rejecting neutral evolution for SCD patients. Extending the statistical tests for detecting natural selection to many non-SCD patients reveal varying levels of selective signal across different tumors, implying that many evolutionary forces including natural selection and geographic isolation can influence the overall pattern of sITH. Taken together, tumor phylogeography unravels a dynamic landscape of sITH, pinpointing important evolutionary and clinical consequences of spatial heterogeneity in cancer.

Sex-biased gene regulation varies across human populations as a result of adaptive evolution

OBJECTIVES

Studies of human sexual dimorphism and gender disparities in health focus on ostensibly universal molecular sex differences, such as sex chromosomes and circulating hormone levels, while ignoring the extraordinary diversity in biology, behavior, and culture acquired by different human populations over their unique evolutionary histories.

MATERIALS AND METHODS

Using RNA-Seq data and whole genome sequences from 1000G and HGDP, we investigate variation in sex-biased gene expression across 11 human populations and test whether population-level variation in sex-biased expression may have resulted from adaptive evolution in regions containing sex-specific regulatory variants.

RESULTS

We find that sex-biased gene expression in humans is highly variable, mostly population-specific, and demonstrates between population reversals. Expression quantitative trait locus mapping reveals sex-specific regulatory regions with evidence of recent positive natural selection, suggesting that variation in sex-biased expression may have evolved as an adaptive response to ancestral environments experienced by human populations.

DISCUSSION

These results indicate that sex-biased gene expression is more flexible than previously thought and is not generally shared among human populations. Instead, molecular phenotypes associated with sex depend on complex interactions between population-specific molecular evolution and physiological responses to contemporary socioecologies.

Reconstructing colonization dynamics to establish how human activities transformed island biodiversity

Drivers and dynamics of initial human migrations across individual islands and archipelagos are poorly understood, hampering assessments of subsequent modification of island biodiversity. We developed and tested a new statistical-simulation approach for reconstructing the pattern and pace of human migration across islands at high spatiotemporal resolutions. Using Polynesian colonisation of New Zealand as an example, we show that process-explicit models, informed by archaeological records and spatiotemporal reconstructions of past climates and environments, can provide new and important insights into the patterns and mechanisms of arrival and establishment of people on islands. We find that colonisation of New Zealand required there to have been a single founding population of approximately 500 people, arriving between 1233 and 1257 AD, settling multiple areas, and expanding rapidly over both North and South Islands. These verified spatiotemporal reconstructions of colonisation dynamics provide new opportunities to explore more extensively the potential ecological impacts of human colonisation on New Zealand's native biota and ecosystems.

An indel introduced by Neanderthal introgression, rs3835124:ATTTATT > ATT, might contribute to prostate cancer risk by regulating PDK1 expression

INTRODUCTION

Prostate cancer is one of the most common cancer types in males and rs12621278:A > G has been suggested to be associated with this disease by previous genome-wide association studies. One thousand genomes project data analysis indicated that rs12621278:A > G is within two long-core haplotypes. However, the origin, causal variant(s), and molecular function of these haplotypes were remaining unclear.

MATERIALS AND METHODS

Population genetics analysis and functional genomics work was performed for this locus.

RESULTS

Phylogeny analysis verified that the rare haplotype is derived from Neanderthal introgression. Genome annotation suggested that three genetic variants in the core haplotypes, rs116108611:G > A, rs139972066:AAAAAAAA > AAAAAAAAA, and rs3835124:ATTTATT > ATT, are located in functional regions. Luciferase assay indicated that rs139972066:AAAAAAAA > AAAAAAAAA and rs116108611:G > A are not able to alter ITGA6 (integrin alpha 6) and ITGA6 antisense RNA 1 expression, respectively. In contrast, rs3835124:ATTTATT > ATT can significantly influence PDK1 (pyruvate dehydrogenase kinase 1) expression, which was verified by expression quantitative trait locus analysis. This genetic variant can alter transcription factor cut like homeobox 1 interaction efficiency. The introgressed haplotype was observed to be subject to positive selection in East Asian populations. The molecular function of the haplotype suggested that Neanderthal should be with lower PDK1 expression and further different energy homeostasis from modern human.

CONCLUSION

This study provided new insight into the contribution of Neanderthal introgression to human phenotypes.

Equine herpesvirus 4 infected domestic horses associated with Sintashta spoke-wheeled chariots around 4,000 years ago

Equine viral outbreaks have disrupted the socio-economic life of past human societies up until the late 19th century and continue to be of major concern to the horse industry today. With a seroprevalence of 60-80 per cent, equine herpesvirus 4 (EHV-4) is the most common horse pathogen on the planet. Yet, its evolutionary history remains understudied. Here, we screen the sequenced data of 264 archaeological horse remains to detect the presence of EHV-4. We recover the first ancient EHV-4 genome with 4.2× average depth-of-coverage from a specimen excavated in the Southeastern Urals and dated to the Early Bronze Age period, approximately 3,900 years ago. The recovery of an EHV-4 virus outside the upper respiratory tract not only points to an animal particularly infected but also highlights the importance of post-cranial bones in pathogen characterisation. Bayesian phylogenetic reconstruction provides a minimal time estimate for EHV-4 diversification to around 4,000 years ago, a time when modern domestic horses spread across the Central Asian steppes together with spoke-wheeled Sintashta chariots, or earlier. The analyses also considerably revise the diversification time of the two EHV-4 subclades from the 16th century based solely on modern data to nearly a thousand years ago. Our study paves the way for a robust reconstruction of the history of non-human pathogens and their impact on animal health.

A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics

The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.

An evolutionary perspective on complex neuropsychiatric disease

The forces of evolution-mutation, selection, migration, and genetic drift-shape the genetic architecture of human traits, including the genetic architecture of complex neuropsychiatric illnesses. Studying these illnesses in populations that are diverse in genetic ancestry, historical demography, and cultural history can reveal how evolutionary forces have guided adaptation over time and place. A fundamental truth of shared human biology is that an allele responsible for a disease in anyone, anywhere, reveals a gene critical to the normal biology underlying that condition in everyone, everywhere. Understanding the genetic causes of neuropsychiatric disease in the widest possible range of human populations thus yields the greatest possible range of insight into genes critical to human brain development. In this perspective, we explore some of the relationships between genes, adaptation, and history that can be illuminated by an evolutionary perspective on studies of complex neuropsychiatric disease in diverse populations.

Speciation-by-Extinction

Extinction is a dominant force shaping patterns of biodiversity through time; however its role as a catalyst of speciation through its interaction with intraspecific variation has been overlooked. Here, we synthesize ideas alluded to by Darwin and others into the model of "speciation-by-extinction" in which speciation results from the extinction of intermediate populations within a single geographically variable species. We explore the properties and distinguishing features of speciation-by-extinction with respect to other established speciation models. We demonstrate its plausibility by showing that the experimental extinction of populations within variable species can result in speciation. The prerequisites for speciation-by-extinction, geographically structured intraspecific variation and local extinction, are ubiquitous in nature. We propose that speciation-by-extinction may be a prevalent, but underappreciated, speciation mechanism.

Undiscovered bird extinctions obscure the true magnitude of human-driven extinction waves

Birds are among the best-studied animal groups, but their prehistoric diversity is poorly known due to low fossilization potential. Hence, while many human-driven bird extinctions (i.e., extinctions caused directly by human activities such as hunting, as well as indirectly through human-associated impacts such as land use change, fire, and the introduction of invasive species) have been recorded, the true number is likely much larger. Here, by combining recorded extinctions with model estimates based on the completeness of the fossil record, we suggest that at least ~1300-1500 bird species (~12% of the total) have gone extinct since the Late Pleistocene, with 55% of these extinctions undiscovered (not yet discovered or left no trace). We estimate that the Pacific accounts for 61% of total bird extinctions. Bird extinction rate varied through time with an intense episode ~1300 CE, which likely represents the largest human-driven vertebrate extinction wave ever, and a rate 80 (60-95) times the background extinction rate. Thus, humans have already driven more than one in nine bird species to extinction, with likely severe, and potentially irreversible, ecological and evolutionary consequences.

Whole genomes from Angola and Mozambique inform about the origins and dispersals of major African migrations

As the continent of origin for our species, Africa harbours the highest levels of diversity anywhere on Earth. However, many regions of Africa remain under-sampled genetically. Here we present 350 whole genomes from Angola and Mozambique belonging to ten Bantu ethnolinguistic groups, enabling the construction of a reference variation catalogue including 2.9 million novel SNPs. We investigate the emergence of Bantu speaker population structure, admixture involving migrations across sub-Saharan Africa and model the demographic histories of Angolan and Mozambican Bantu speakers. Our results bring together concordant views from genomics, archaeology, and linguistics to paint an updated view of the complexity of the Bantu Expansion. Moreover, we generate reference panels that better represents the diversity of African populations involved in the trans-Atlantic slave trade, improving imputation accuracy in African Americans and Brazilians. We anticipate that our collection of genomes will form the foundation for future African genomic healthcare initiatives.

Haplotype-based inference of recent effective population size in modern and ancient DNA samples

Individuals sharing recent ancestors are likely to co-inherit large identical-by-descent (IBD) genomic regions. The distribution of these IBD segments in a population may be used to reconstruct past demographic events such as effective population size variation, but accurate IBD detection is difficult in ancient DNA data and in underrepresented populations with limited reference data. In this work, we introduce an accurate method for inferring effective population size variation during the past ~2000 years in both modern and ancient DNA data, called HapNe. HapNe infers recent population size fluctuations using either IBD sharing (HapNe-IBD) or linkage disequilibrium (HapNe-LD), which does not require phasing and can be computed in low coverage data, including data sets with heterogeneous sampling times. HapNe shows improved accuracy in a range of simulated demographic scenarios compared to currently available methods for IBD-based and LD-based inference of recent effective population size, while requiring fewer computational resources. We apply HapNe to several modern populations from the 1,000 Genomes Project, the UK Biobank, the Allen Ancient DNA Resource, and recently published samples from Iron Age Britain, detecting multiple instances of recent effective population size variation across these groups.

Indigenous Australian genomes show deep structure and rich novel variation

The Indigenous peoples of Australia have a rich linguistic and cultural history. How this relates to genetic diversity remains largely unknown because of their limited engagement with genomic studies. Here we analyse the genomes of 159 individuals from four remote Indigenous communities, including people who speak a language (Tiwi) not from the most widespread family (Pama-Nyungan). This large collection of Indigenous Australian genomes was made possible by careful community engagement and consultation. We observe exceptionally strong population structure across Australia, driven by divergence times between communities of 26,000-35,000 years ago and long-term low but stable effective population sizes. This demographic history, including early divergence from Papua New Guinean (47,000 years ago) and Eurasian groups1, has generated the highest proportion of previously undescribed genetic variation seen outside Africa and the most extended homozygosity compared with global samples. A substantial proportion of this variation is not observed in global reference panels or clinical datasets, and variation with predicted functional consequence is more likely to be homozygous than in other populations, with consequent implications for medical genomics2. Our results show that Indigenous Australians are not a single homogeneous genetic group and their genetic relationship with the peoples of New Guinea is not uniform. These patterns imply that the full breadth of Indigenous Australian genetic diversity remains uncharacterized, potentially limiting genomic medicine and equitable healthcare for Indigenous Australians.

Revolutionizing Biological Science: The Synergy of Genomics in Health, Bioinformatics, Agriculture, and Artificial Intelligence

With climate emergency, COVID-19, and the rise of planetary health scholarship, the binary of human and ecosystem health has been deeply challenged. The interdependence of human and nonhuman animal health is increasingly acknowledged and paving the way for new frontiers in integrative biology. The convergence of genomics in health, bioinformatics, agriculture, and artificial intelligence (AI) has ushered in a new era of possibilities and applications. However, the sheer volume of genomic/multiomics big data generated also presents formidable sociotechnical challenges in extracting meaningful biological, planetary health and ecological insights. Over the past few years, AI-guided bioinformatics has emerged as a powerful tool for managing, analyzing, and interpreting complex biological datasets. The advances in AI, particularly in machine learning and deep learning, have been transforming the fields of genomics, planetary health, and agriculture. This article aims to unpack and explore the formidable range of possibilities and challenges that result from such transdisciplinary integration, and emphasizes its radically transformative potential for human and ecosystem health. The integration of these disciplines is also driving significant advancements in precision medicine and personalized health care. This presents an unprecedented opportunity to deepen our understanding of complex biological systems and advance the well-being of all life in planetary ecosystems. Notwithstanding in mind its sociotechnical, ethical, and critical policy challenges, the integration of genomics, multiomics, planetary health, and agriculture with AI-guided bioinformatics opens up vast opportunities for transnational collaborative efforts, data sharing, analysis, valorization, and interdisciplinary innovations in life sciences and integrative biology.

Archaic Introgression Shaped Human Circadian Traits

When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely, Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultraviolet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology and whether archaic introgression adaptively contributed to human chronotypes remain unknown. Here, we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, and RORC) and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among expression quantitative trait loci for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency. These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.

Ethnopsychopharmacology: Clinical and scientific writing pearls

The concept of ethnopsychopharmacology aims to predict or explain the pharmacologic response to psychiatric medications based on the influence of biologic and nonbiologic factors. Interactions involving these factors are complex and influence patient outcomes in health care. Pharmacists and other clinicians working in patient care environments, research, or medical education should engage in lifelong learning to enhance ethnopsychopharmacologic knowledge gaps, which ultimately may improve and individualize care across diverse populations. Through two cases, this paper provides pearls on how biogeographical ancestry and cytochrome P450 status may influence pharmacotherapy selection, dosing, or response. A third scenario highlights a publication, like many other published works, with deficiencies in how data on ancestry, race, and ethnicity are collected or reported. Current recommendations on the use of inclusive language in scientific writing are reviewed, with attention to specific examples.

Getting it right: Teaching undergraduate biology to undermine racial essentialism

How we teach human genetics matters for social equity. The biology curriculum appears to be a crucial locus of intervention for either reinforcing or undermining students' racial essentialist views. The Mendelian genetic models dominating textbooks, particularly in combination with racially inflected language sometimes used when teaching about monogenic disorders, can increase middle and high school students' racial essentialism and opposition to policies to increase equity. These findings are of particular concern given the increasing spread of racist misinformation online and the misappropriation of human genomics research by white supremacists, who take advantage of low levels of genetics literacy in the general public. Encouragingly, however, teaching updated information about the geographical distribution of human genetic variation and the complex, multifactorial basis of most human traits, reduces students' endorsement of racial essentialism. The genetics curriculum is therefore a key tool in combating misinformation and scientific racism. Here, we describe a framework and example teaching materials for teaching students key concepts in genetics, human evolutionary history, and human phenotypic variation at the undergraduate level. This framework can be flexibly applied in biology and anthropology classes and adjusted based on time availability. Our goal is to provide undergraduate-level instructors with varying levels of expertise with a set of evidence-informed tools for teaching human genetics to combat scientific racism, including an evolving set of instructional resources, as well as learning goals and pedagogical approaches. Resources can be found at https://noto.li/YIlhZ5. Additionally, we hope to generate conversation about integrating modern genetics into the undergraduate curriculum, in light of recent findings about the risks and opportunities associated with teaching genetics.

Nuclear genetic diversity of head lice sheds light on human dispersal around the world

The human louse, Pediculus humanus, is an obligate blood-sucking ectoparasite that has coevolved with humans for millennia. Given the intimate relationship between this parasite and the human host, the study of human lice has the potential to shed light on aspects of human evolution that are difficult to interpret using other biological evidence. In this study, we analyzed the genetic variation in 274 human lice from 25 geographic sites around the world by using nuclear microsatellite loci and female-inherited mitochondrial DNA sequences. Nuclear genetic diversity analysis revealed the presence of two distinct genetic clusters I and II, which are subdivided into subclusters: Ia-Ib and IIa-IIb, respectively. Among these samples, we observed the presence of the two most common louse mitochondrial haplogroups: A and B that were found in both nuclear Clusters I and II. Evidence of nuclear admixture was uncommon (12%) and was predominate in the New World potentially mirroring the history of colonization in the Americas. These findings were supported by novel DIYABC simulations that were built using both host and parasite data to define parameters and models suggesting that admixture between cI and cII was very recent. This pattern could also be the result of a reproductive barrier between these two nuclear genetic clusters. In addition to providing new evolutionary knowledge about this human parasite, our study could guide the development of new analyses in other host-parasite systems.

Efficient Bayesian inference under the multispecies coalescent with migration

Analyses of genome sequence data have revealed pervasive interspecific gene flow and enriched our understanding of the role of gene flow in speciation and adaptation. Inference of gene flow using genomic data requires powerful statistical methods. Yet current likelihood-based methods involve heavy computation and are feasible for small datasets only. Here, we implement the multispecies-coalescent-with-migration model in the Bayesian program bpp, which can be used to test for gene flow and estimate migration rates, as well as species divergence times and population sizes. We develop Markov chain Monte Carlo algorithms for efficient sampling from the posterior, enabling the analysis of genome-scale datasets with thousands of loci. Implementation of both introgression and migration models in the same program allows us to test whether gene flow occurred continuously over time or in pulses. Analyses of genomic data from Anopheles mosquitoes demonstrate rich information in typical genomic datasets about the mode and rate of gene flow.

Past human expansions shaped the spatial pattern of Neanderthal ancestry

The worldwide expansion of modern humans (Homo sapiens) started before the extinction of Neanderthals (Homo neanderthalensis). Both species coexisted and interbred, leading to slightly higher introgression in East Asians than in Europeans. This distinct ancestry level has been argued to result from selection, but range expansions of modern humans could provide an alternative explanation. This hypothesis would lead to spatial introgression gradients, increasing with distance from the expansion source. We investigate the presence of Neanderthal introgression gradients after past human expansions by analyzing Eurasian paleogenomes. We show that the out-of-Africa expansion resulted in spatial gradients of Neanderthal ancestry that persisted through time. While keeping the same gradient orientation, the expansion of early Neolithic farmers contributed decisively to reducing the Neanderthal introgression in European populations compared to Asian populations. This is because Neolithic farmers carried less Neanderthal DNA than preceding Paleolithic hunter-gatherers. This study shows that inferences about past human population dynamics can be made from the spatiotemporal variation in archaic introgression.

Demographic modeling of admixed Latin American populations from whole genomes

Demographic models of Latin American populations often fail to fully capture their complex evolutionary history, which has been shaped by both recent admixture and deeper-in-time demographic events. To address this gap, we used high-coverage whole-genome data from Indigenous American ancestries in present-day Mexico and existing genomes from across Latin America to infer multiple demographic models that capture the impact of different timescales on genetic diversity. Our approach, which combines analyses of allele frequencies and ancestry tract length distributions, represents a significant improvement over current models in predicting patterns of genetic variation in admixed Latin American populations. We jointly modeled the contribution of European, African, East Asian, and Indigenous American ancestries into present-day Latin American populations. We infer that the ancestors of Indigenous Americans and East Asians diverged ∼30 thousand years ago, and we characterize genetic contributions of recent migrations from East and Southeast Asia to Peru and Mexico. Our inferred demographic histories are consistent across different genomic regions and annotations, suggesting that our inferences are robust to the potential effects of linked selection. In conjunction with published distributions of fitness effects for new nonsynonymous mutations in humans, we show in large-scale simulations that our models recover important features of both neutral and deleterious variation. By providing a more realistic framework for understanding the evolutionary history of Latin American populations, our models can help address the historical under-representation of admixed groups in genomics research and can be a valuable resource for future studies of populations with complex admixture and demographic histories.

The complete and fully-phased diploid genome of a male Han Chinese

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

Assessing HLA imputation accuracy in a West African population

The Human Leukocyte Antigen (HLA) region plays an important role in autoimmune and infectious diseases. HLA is a highly polymorphic region and thus difficult to impute. We, therefore, sought to evaluate HLA imputation accuracy, specifically in a West African population, since they are understudied and are known to harbor high genetic diversity. The study sets were selected from 315 Gambian individuals within the Gambian Genome Variation Project (GGVP) Whole Genome Sequence datasets. Two different arrays, Illumina Omni 2.5 and Human Hereditary and Health in Africa (H3Africa), were assessed for the appropriateness of their markers, and these were used to test several imputation panels and tools. The reference panels were chosen from the 1000 Genomes (1kg-All), 1000 Genomes African (1kg-Afr), 1000 Genomes Gambian (1kg-Gwd), H3Africa, and the HLA Multi-ethnic datasets. HLA-A, HLA-B, and HLA-C alleles were imputed using HIBAG, SNP2HLA, CookHLA, and Minimac4, and concordance rate was used as an assessment metric. The best performing tool was found to be HIBAG, with a concordance rate of 0.84, while the best performing reference panel was the H3Africa panel, with a concordance rate of 0.62. Minimac4 (0.75) was shown to increase HLA-B allele imputation accuracy compared to HIBAG (0.71), SNP2HLA (0.51) and CookHLA (0.17). The H3Africa and Illumina Omni 2.5 array performances were comparable, showing that genotyping arrays have less influence on HLA imputation in West African populations. The findings show that using a larger population-specific reference panel and the HIBAG tool improves the accuracy of HLA imputation in a West African population.

An Evolutionary Model for the Ancient Origins of Polycystic Ovary Syndrome

Polycystic ovary syndrome (PCOS) is a common endocrinopathy of reproductive-aged women, characterized by hyperandrogenism, oligo-anovulation and insulin resistance and closely linked with preferential abdominal fat accumulation. As an ancestral primate trait, PCOS was likely further selected in humans when scarcity of food in hunter-gatherers of the late Pleistocene additionally programmed for enhanced fat storage to meet the metabolic demands of reproduction in later life. As an evolutionary model for PCOS, healthy normal-weight women with hyperandrogenic PCOS have subcutaneous (SC) abdominal adipose stem cells that favor fat storage through exaggerated lipid accumulation during development to adipocytes in vitro. In turn, fat storage is counterbalanced by reduced insulin sensitivity and preferential accumulation of highly lipolytic intra-abdominal fat in vivo. This metabolic adaptation in PCOS balances energy storage with glucose availability and fatty acid oxidation for optimal energy use during reproduction; its accompanying oligo-anovulation allowed PCOS women from antiquity sufficient time and strength for childrearing of fewer offspring with a greater likelihood of childhood survival. Heritable PCOS characteristics are affected by today's contemporary environment through epigenetic events that predispose women to lipotoxicity, with excess weight gain and pregnancy complications, calling for an emphasis on preventive healthcare to optimize the long-term, endocrine-metabolic health of PCOS women in today's obesogenic environment.

Archaic Introgression Shaped Human Circadian Traits

Introduction

When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultra-violet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology, and whether archaic introgression adaptively contributed to human chronotypes remains unknown.

Results

Here we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, RORC), and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among eQTLs for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency.

Conclusions

These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.

Genomic inference of a severe human bottleneck during the Early to Middle Pleistocene transition

Population size history is essential for studying human evolution. However, ancient population size history during the Pleistocene is notoriously difficult to unravel. In this study, we developed a fast infinitesimal time coalescent process (FitCoal) to circumvent this difficulty and calculated the composite likelihood for present-day human genomic sequences of 3154 individuals. Results showed that human ancestors went through a severe population bottleneck with about 1280 breeding individuals between around 930,000 and 813,000 years ago. The bottleneck lasted for about 117,000 years and brought human ancestors close to extinction. This bottleneck is congruent with a substantial chronological gap in the available African and Eurasian fossil record. Our results provide new insights into our ancestry and suggest a coincident speciation event.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

A global view of the genetic basis of Alzheimer disease

The risk of Alzheimer disease (AD) increases with age, family history and informative genetic variants. Sadly, there is still no cure or means of prevention. As in other complex diseases, uncovering genetic causes of AD could identify underlying pathological mechanisms and lead to potential treatments. Rare, autosomal dominant forms of AD occur in middle age as a result of highly penetrant genetic mutations, but the most common form of AD occurs later in life. Large-scale, genome-wide analyses indicate that 70 or more genes or loci contribute to AD. One of the major factors limiting progress is that most genetic data have been obtained from non-Hispanic white individuals in Europe and North America, preventing the development of personalized approaches to AD in individuals of other ethnicities. Fortunately, emerging genetic data from other regions - including Africa, Asia, India and South America - are now providing information on the disease from a broader range of ethnicities. Here, we summarize the current knowledge on AD genetics in populations across the world. We predominantly focus on replicated genetic discoveries but also include studies in ethnic groups where replication might not be feasible. We attempt to identify gaps that need to be addressed to achieve a complete picture of the genetic and molecular factors that drive AD in individuals across the globe.

On whole-genome demography of world's ethnic groups and individual genomic identity

All current categorizations of human population, such as ethnicity, ancestry and race, are based on various selections and combinations of complex and dynamic common characteristics, that are mostly societal and cultural in nature, perceived by the members within or from outside of the categorized group. During the last decade, a massive amount of a new type of characteristics, that are exclusively genomic in nature, became available that allows us to analyze the inherited whole-genome demographics of extant human, especially in the fields such as human genetics, health sciences and medical practices (e.g., 1,2,3), where such health-related characteristics can be related to whole-genome-based categorization. Here we show the feasibility of deriving such whole-genome-based categorization. We observe that, within the available genomic data at present, (a) the study populations form about 14 genomic groups, each consisting of multiple ethnic groups; and (b), at an individual level, approximately 99.8%, on average, of the whole autosomal-genome contents are identical between any two individuals regardless of their genomic or ethnic groups.

The Increase in Hemoglobin Concentration With Altitude Differs Between World Regions and Is Less in Children Than in Adults

To compensate for decreased oxygen partial pressure, high-altitude residents increase hemoglobin concentrations [Hb]. The elevation varies between world regions, posing problems in defining cutoff values for anemia or polycythemia. The currently used altitude adjustments (World Health Organization [WHO]), however, do not account for regional differences. Data from The Demographic and Health Survey (DHS) Program were analyzed from 32 countries harboring >4% of residents at altitudes above 1000 m. [Hb]-increase, (ΔHb/km altitude) was calculated by linear regression analysis. Tables show 95% reference intervals (RIs) for different altitude ranges, world regions, and age groups. The prevalence of anemia and polycythemia was calculated using regressions in comparison to WHO adjustments. The most pronounced Δ[Hb]/km was found in East Africans and South Americans while [Hb] increased least in South/South-East Asia. In African regions and Middle East, [Hb] was decreased in some altitude regions showing inconsistent changes in different age groups. Of note, in all regions, the Δ[Hb]/km was lower in children than in adults, and in the Middle East, it was even negative. Overall, the Δ[Hb]/km from our analysis differed from the region-independent adjustments currently suggested by the WHO resulting in a lower anemia prevalence at very high altitudes. The distinct patterns of Δ[Hb] with altitude in residents from different world regions imply that one single, region-independent correction factor for altitude is not be applicable for diagnosing abnormal [Hb]. Therefore, we provide regression coefficients and reference-tables that are specific for world regions and altitude ranges to improve diagnosing abnormal [Hb].

Spatiotemporal fluctuations of population structure in the Americas revealed by a meta-analysis of the first decade of archaeogenomes

Objectives

Since 2010, genome-wide data from hundreds of ancient Native Americans have contributed to the understanding of Americas' prehistory. However, these samples have never been studied as a single dataset, and distinct relationships among themselves and with present-day populations may have never come to light. Here, we reassess genomic diversity and population structure of 223 ancient Native Americans published between 2010 and 2019.

Materials and Methods

The genomic data from ancient Americas was merged with a worldwide reference panel of 278 present-day genomes from the Simons Genome Diversity Project and then analyzed through ADMIXTURE, D-statistics, PCA, t-SNE, and UMAP.

Results

We find largely similar population structures in ancient and present-day Americas. However, the population structure of contemporary Native Americans, traced here to at least 10,000 years before present, is noticeably less diverse than their ancient counterparts, a possible outcome of the European contact. Additionally, in the past there were greater levels of population structure in North than in South America, except for ancient Brazil, which harbors comparatively high degrees of structure. Moreover, we find a component of genetic ancestry in the ancient dataset that is closely related to that of present-day Oceanic populations but does not correspond to the previously reported Australasian signal. Lastly, we report an expansion of the Ancient Beringian ancestry, previously reported for only one sample.

Discussion

Overall, our findings support a complex scenario for the settlement of the Americas, accommodating the occurrence of founder effects and the emergence of ancestral mixing events at the regional level.

Modern Japanese ancestry-derived variants reveal the formation process of the current Japanese regional gradations

Modern Japanese people have two major ancestral populations: indigenous Jomon hunter-gatherers and continental East Asian farmers. To determine the formation process of the current Japanese population, we developed a detection method for variants derived from ancestral populations using a summary statistic, the ancestry marker index (AMI). We applied AMI to modern Japanese population samples and identified 208,648 single nucleotide polymorphisms (SNPs) that were likely derived from the Jomon people (Jomon-derived variants). Analysis of Jomon-derived variants in 10,842 modern Japanese individuals recruited from all over Japan revealed that the admixture proportions of the Jomon people varied between prefectures, probably owing to the prehistoric population size difference. The estimated allele frequencies of genome-wide SNPs in the ancestral populations of the modern Japanese suggested their adaptive phenotypic characteristics to their respective livelihoods. Based on our findings, we propose a formation model for the genotypic and phenotypic gradations of the current Japanese archipelago populations.

Whole-genome phylogeography of the blue-faced honeyeater (Entomyzon cyanotis) and discovery and characterization of a neo-Z chromosome

Whole-genome surveys of genetic diversity and geographic variation often yield unexpected discoveries of novel structural variation, which long-read DNA sequencing can help clarify. Here, we report on whole-genome phylogeography of a bird exhibiting classic vicariant geographies across Australia and New Guinea, the blue-faced honeyeater (Entomyzon cyanotis), and the discovery and characterization of a novel neo-Z chromosome by long-read sequencing. Using short-read genome-wide SNPs, we inferred population divergence events within E. cyanotis across the Carpentarian and other biogeographic barriers during the Pleistocene (~0.3-1.7 Ma). Evidence for introgression between nonsister populations supports a hypothesis of reticulate evolution around a triad of dynamic barriers around Pleistocene Lake Carpentaria between Australia and New Guinea. During this phylogeographic survey, we discovered a large (134 Mbp) neo-Z chromosome and we explored its diversity, divergence and introgression landscape. We show that, as in some sylvioid passerine birds, a fusion occurred between chromosome 5 and the Z chromosome to form a neo-Z chromosome; and in E. cyanotis, the ancestral pseudoautosomal region (PAR) appears nonrecombinant between Z and W, along with most of the fused chromosome 5. The added recombination-suppressed portion of the neo-Z (~37.2 Mbp) displays reduced diversity and faster population genetic differentiation compared with the ancestral-Z. Yet, the new PAR (~17.4 Mbp) shows elevated diversity and reduced differentiation compared to autosomes, potentially resulting from introgression. In our case, long-read sequencing helped clarify the genomic landscape of population divergence on autosomes and sex chromosomes in a species where prior knowledge of genome structure was still incomplete.

Advancements and Challenges in Ancient DNA Research: Bridging the Global North-South Divide

Ancient DNA (aDNA) research first began in 1984 and ever since has greatly expanded our understanding of evolution and migration. Today, aDNA analysis is used to solve various puzzles about the origin of mankind, migration patterns, and the spread of infectious diseases. The incredible findings ranging from identifying the new branches within the human family to studying the genomes of extinct flora and fauna have caught the world by surprise in recent times. However, a closer look at these published results points out a clear Global North and Global South divide. Therefore, through this research, we aim to emphasize encouraging better collaborative opportunities and technology transfer to support researchers in the Global South. Further, the present research also focuses on expanding the scope of the ongoing conversation in the field of aDNA by reporting relevant literature published around the world and discussing the advancements and challenges in the field.

The evolution of human skin pigmentation: A changing medley of vitamins, genetic variability, and UV radiation during human expansion

This review examines putative, yet likely critical evolutionary pressures contributing to human skin pigmentation and subsequently, depigmentation phenotypes. To achieve this, it provides a synthesis of ideas that frame contemporary thinking, without limiting the narrative to pigmentation genes alone. It examines how geography and hence the quality and quantity of UV exposure, pigmentation genes, diet-related genes, vitamins, anti-oxidant nutrients, and cultural practices intersect and interact to facilitate the evolution of human skin color. The article has a strong focus on the vitamin D-folate evolutionary model, with updates on the latest biophysical research findings to support this paradigm. This model is examined within a broad canvas that takes human expansion out of Africa and genetic architecture into account. A thorough discourse on the biology of melanization is provided (includes relationship to BH4 and DNA damage repair), with the relevance of this to the UV sensitivity of folate and UV photosynthesis of vitamin D explained in detail, including the relevance of these vitamins to reproductive success. It explores whether we might be able to predict vitamin-related gene polymorphisms that pivot metabolism to the prevailing UVR exposome within the vitamin D-folate evolutionary hypothesis context. This is discussed in terms of a primary adaptive phenotype (pigmentation/depigmentation), a secondary adaptive phenotype (flexible metabolic phenotype based on vitamin-related gene polymorphism profile), and a tertiary adaptive strategy (dietary anti-oxidants to support the secondary adaptive phenotype). Finally, alternative evolutionary models for pigmentation are discussed, as are challenges to future research in this area.

Assessing HLA imputation accuracy in a West African population

The Human Leukocyte Antigen (HLA) region plays an important role in autoimmune and infectious diseases. HLA is a highly polymorphic region and thus difficult to impute. We therefore sought to evaluate HLA imputation accuracy, specifically in a West African population, since they are understudied and are known to harbor high genetic diversity. The study sets were selected from Gambian individuals within the Gambian Genome Variation Project (GGVP) Whole Genome Sequence datasets. Two different arrays, Illumina Omni 2.5 and Human Hereditary and Health in Africa (H3Africa), were assessed for the appropriateness of their markers, and these were used to test several imputation panels and tools. The reference panels were chosen from the 1000 Genomes dataset (1kg-All), 1000 Genomes African dataset (1kg-Afr), 1000 Genomes Gambian dataset (1kg-Gwd), H3Africa dataset and the HLA Multi-ethnic dataset. HLA-A, HLA-B and HLA-C alleles were imputed using HIBAG, SNP2HLA, CookHLA and Minimac4, and concordance rate was used as an assessment metric. Overall, the best performing tool was found to be HIBAG, with a concordance rate of 0.84, while the best performing reference panel was the H3Africa panel with a concordance rate of 0.62. Minimac4 (0.75) was shown to increase HLA-B allele imputation accuracy compared to HIBAG (0.71), SNP2HLA (0.51) and CookHLA (0.17). The H3Africa and Illumina Omni 2.5 array performances were comparable, showing that genotyping arrays have less influence on HLA imputation in West African populations. The findings show that using a larger population-specific reference panel and the HIBAG tool improves the accuracy of HLA imputation in West African populations.

Author Summary

For studies that associate a particular HLA type to a phenotypic trait for instance HIV susceptibility or control, genotype imputation remains the main method for acquiring a larger sample size. Genotype imputation, process of inferring unobserved genotypes, is a statistical technique and thus deals with probabilities. Also, the HLA region is highly variable and therefore difficult to impute. In view of this, it is important to assess HLA imputation accuracy especially in African populations. This is because the African genome has high diversity, and such studies have hardly been conducted in African populations. This work highlights that using HIBAG imputation tool and a larger population-specific reference panel increases HLA imputation accuracy in an African population.

A review of ancestrality and admixture in Latin America and the caribbean focusing on native American and African descendant populations

Genomics can reveal essential features about the demographic evolution of a population that may not be apparent from historical elements. In recent years, there has been a significant increase in the number of studies applying genomic epidemiological approaches to understand the genetic structure and diversity of human populations in the context of demographic history and for implementing precision medicine. These efforts have traditionally been applied predominantly to populations of European origin. More recently, initiatives in the United States and Africa are including more diverse populations, establishing new horizons for research in human populations with African and/or Native ancestries. Still, even in the most recent projects, the under-representation of genomic data from Latin America and the Caribbean (LAC) is remarkable. In addition, because the region presents the most recent global miscegenation, genomics data from LAC may add relevant information to understand population admixture better. Admixture in LAC started during the colonial period, in the 15th century, with intense miscegenation between European settlers, mainly from Portugal and Spain, with local indigenous and sub-Saharan Africans brought through the slave trade. Since, there are descendants of formerly enslaved and Native American populations in the LAC territory; they are considered vulnerable populations because of their history and current living conditions. In this context, studying LAC Native American and African descendant populations is important for several reasons. First, studying human populations from different origins makes it possible to understand the diversity of the human genome better. Second, it also has an immediate application to these populations, such as empowering communities with the knowledge of their ancestral origins. Furthermore, because knowledge of the population genomic structure is an essential requirement for implementing genomic medicine and precision health practices, population genomics studies may ensure that these communities have access to genomic information for risk assessment, prevention, and the delivery of optimized treatment; thus, helping to reduce inequalities in the Western Hemisphere. Hoping to set the stage for future studies, we review different aspects related to genetic and genomic research in vulnerable populations from LAC countries.

ACL microtrauma: healing through nutrition, modified sports training, and increased recovery time

PURPOSE

Sports injuries among youth and adolescent athletes are a growing concern, particularly at the knee. Based on our current understanding of microtrauma and anterior cruciate ligament (ACL) healing characteristics, this clinical commentary describes a comprehensive plan to better manage ACL microtrauma and mitigate the likelihood of progression to a non-contact macrotraumatic ACL rupture.

METHODS

Medical literature related to non-contact ACL injuries among youth and adolescent athletes, collagen and ACL extracellular matrix metabolism, ACL microtrauma and sudden failure, and concerns related to current sports training were reviewed and synthesized into a comprehensive intervention plan.

RESULTS

With consideration for biopsychosocial model health factors, proper nutrition and modified sports training with increased recovery time, a comprehensive primary ACL injury prevention plan is described for the purpose of better managing ACL microtrauma, thereby reducing the incidence of non-contact macrotraumatic ACL rupture among youth and adolescent athletes.

CONCLUSION

Preventing non-contact ACL injuries may require greater consideration for reducing accumulated ACL microtrauma. Proper nutrition including glycine-rich collagen peptides, or gelatin-vitamin C supplementation in combination with healthy sleep, and adjusted sports training periodization with increased recovery time may improve ACL extracellular matrix collagen deposition homeostasis, decreasing sudden non-contact ACL rupture incidence likelihood in youth and adolescent athletes. Successful implementation will require compliance from athletes, parents, coaches, the sports medicine healthcare team, and event organizers. Studies are needed to confirm the efficacy of these concepts.

LEVEL OF EVIDENCE

V.

Disease-causing human viruses: novelty and legacy

About 270 viruses are known to infect humans. Some of these viruses have been known for centuries, whereas others have recently emerged. During their evolutionary history, humans have moved out of Africa to populate the world. In historical times, human migrations resulted in the displacement of large numbers of people. All these events determined the movement and dispersal of human-infecting viruses. Technological advances have resulted in the characterization of the genetic variability of human viruses, both in extant and in archaeological samples. Field studies investigated the diversity of viruses hosted by other animals. In turn, these advances provided insight into the evolutionary history of human viruses back in time and defined the key events through which they originated and spread.

Mosaic patterns of selection in genomic regions associated with diverse human traits

Natural selection shapes the genetic architecture of many human traits. However, the prevalence of different modes of selection on genomic regions associated with variation in traits remains poorly understood. To address this, we developed an efficient computational framework to calculate positive and negative enrichment of different evolutionary measures among regions associated with complex traits. We applied the framework to summary statistics from >900 genome-wide association studies (GWASs) and 11 evolutionary measures of sequence constraint, population differentiation, and allele age while accounting for linkage disequilibrium, allele frequency, and other potential confounders. We demonstrate that this framework yields consistent results across GWASs with variable sample sizes, numbers of trait-associated SNPs, and analytical approaches. The resulting evolutionary atlas maps diverse signatures of selection on genomic regions associated with complex human traits on an unprecedented scale. We detected positive enrichment for sequence conservation among trait-associated regions for the majority of traits (>77% of 290 high power GWASs), which included reproductive traits. Many traits also exhibited substantial positive enrichment for population differentiation, especially among hair, skin, and pigmentation traits. In contrast, we detected widespread negative enrichment for signatures of balancing selection (51% of GWASs) and absence of enrichment for evolutionary signals in regions associated with late-onset Alzheimer's disease. These results support a pervasive role for negative selection on regions of the human genome that contribute to variation in complex traits, but also demonstrate that diverse modes of evolution are likely to have shaped trait-associated loci. This atlas of evolutionary signatures across the diversity of available GWASs will enable exploration of the relationship between the genetic architecture and evolutionary processes in the human genome.