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

Understanding the variant landscape, and genetic epidemiology of Multiple Endocrine Neoplasia in India

PURPOSE

Multiple Endocrine Neoplasia (MEN) is a group of familial cancer syndromes that encompasses several types of endocrine tumors differentiated by genetic mutations in RET, MEN1 and CDKN1B genes. Accurate diagnosis of MEN subtypes can thus be performed through genetic testing. However, MEN variants remain largely understudied in Indian populations. Additionally, few dedicated resources to understand these disorders currently exist.

METHODS

Using the gold-standard ACMG/AMP guidelines, we systematically classified variants reported across the three genes in the IndiGen dataset, and established the genetic epidemiology of MEN in the Indian population. We further classified ClinVar and Mastermind variants and compiled all into a database. Finally, we designed a multiplex primer panel for rapid variant identification.

RESULTS

We have established the genetic prevalence of MEN as the following: 1 in 1026 individuals is likely to be afflicted with MEN linked with pathogenic RET mutations. We have further created the MAPVar database containing 3280 ACMG-classified variants freely accessible at: https://clingen.igib.res.in/MAPVar/ . Finally, our NGS primer panel covers 33 exonic regions across two pools through 38 amplicons with a total amplified region of 65 kb.

CONCLUSION

Our work establishes that MEN is a prevalent disorder in India. The rare nature of Indian variants underscores the need of genomic and functional studies to establish a more comprehensive variant landscape. Additionally, our panel offers a means of cost-effective genetic testing, and the MAPVar database a ready reference to aid in a better understanding of variant pathogenicity in clinical as well as research settings.

Pharmacogenomic insights in psychiatric care: uncovering novel actionability, allele-specific CYP2D6 copy number variation, and phenoconversion in 15,000 patients

Pharmacogenomic testing has emerged as an aid in clinical decision making for psychiatric providers, but more data is needed regarding its utility in clinical practice and potential impact on patient care. In this cross-sectional study, we determined the real-world prevalence of pharmacogenomic actionability in patients receiving psychiatric care. Potential actionability was based on the prevalence of CYP2C19 and CYP2D6 phenotypes, including CYP2D6 allele-specific copy number variations (CNVs). Combined actionability additionally incorporated CYP2D6 phenoconversion and the novel CYP2C-TG haplotype in patients with available medication data. Across 15,000 patients receiving clinical pharmacogenomic testing, 65% had potentially actionable CYP2D6 and CYP2C19 phenotypes, and phenotype assignment was impacted by CYP2D6 allele-specific CNVs in 2% of all patients. Of 4114 patients with medication data, 42% had CYP2D6 phenoconversion from drug interactions and 20% carried a novel CYP2C haplotype potentially altering actionability. A total of 87% had some form of potential actionability from genetic findings and/or phenoconversion. Genetic variation detected via next-generation sequencing led to phenotype reassignment in 22% of individuals overall (2% in CYP2D6 and 20% in CYP2C19). Ultimately, pharmacogenomic testing using next-generation sequencing identified potential actionability in most patients receiving psychiatric care. Early pharmacogenomic testing may provide actionable insights to aid clinicians in drug prescribing to optimize psychiatric care.

9,000 years of genetic continuity in southernmost Africa demonstrated at Oakhurst rockshelter

Southern Africa has one of the longest records of fossil hominins and harbours the largest human genetic diversity in the world. Yet, despite its relevance for human origins and spread around the globe, the formation and processes of its gene pool in the past are still largely unknown. Here, we present a time transect of genome-wide sequences from nine individuals recovered from a single site in South Africa, Oakhurst Rockshelter. Spanning the whole Holocene, the ancient DNA of these individuals allows us to reconstruct the demographic trajectories of the indigenous San population and their ancestors during the last 10,000 years. We show that, in contrast to most regions around the world, the population history of southernmost Africa was not characterized by several waves of migration, replacement and admixture but by long-lasting genetic continuity from the early Holocene to the end of the Later Stone Age. Although the advent of pastoralism and farming substantially transformed the gene pool in most parts of southern Africa after 1,300 BP, we demonstrate using allele-frequency and identity-by-descent segment-based methods that the ‡Khomani San and Karretjiemense from South Africa still show direct signs of relatedness to the Oakhurst hunter-gatherers, a pattern obscured by recent, extensive non-Southern African admixture. Yet, some southern San in South Africa still preserve this ancient, Pleistocene-derived genetic signature, extending the period of genetic continuity until today.

An archaic HLA class I receptor allele diversifies natural killer cell-driven immunity in First Nations peoples of Oceania

Genetic variation in host immunity impacts the disproportionate burden of infectious diseases that can be experienced by First Nations peoples. Polymorphic human leukocyte antigen (HLA) class I and killer cell immunoglobulin-like receptors (KIRs) are key regulators of natural killer (NK) cells, which mediate early infection control. How this variation impacts their responses across populations is unclear. We show that HLA-A∗24:02 became the dominant ligand for inhibitory KIR3DL1 in First Nations peoples across Oceania, through positive natural selection. We identify KIR3DL1∗114, widespread across and unique to Oceania, as an allele lineage derived from archaic humans. KIR3DL1∗114+NK cells from First Nations Australian donors are inhibited through binding HLA-A∗24:02. The KIR3DL1∗114 lineage is defined by phenylalanine at residue 166. Structural and binding studies show phenylalanine 166 forms multiple unique contacts with HLA-peptide complexes, increasing both affinity and specificity. Accordingly, assessing immunogenetic variation and the functional implications for immunity are fundamental toward understanding population-based disease associations.

Genotyping sequence-resolved copy-number variations using pangenomes reveals paralog-specific global diversity and expression divergence of duplicated genes

Copy-number variable (CNV) genes are important in evolution and disease, yet sequence variation in CNV genes is a blindspot for large-scale studies. We present a method, ctyper, that leverages pangenomes to produce copy-number maps with allele-specific sequences containing locally phased variants of CNV genes from NGS reads. We extensively characterized accuracy and efficiency on a database of 3,351 CNV genes including HLA, SMN, and CYP2D6 as well as 212 non-CNV medically-relevant challenging genes. The genotypes capture 96.5% of underlying variants in new genomes, requiring 0.9 seconds per gene. Expression analysis of ctyper genotypes explains more variance than known eQTL variants. Comparing allele-specific expression quantified divergent expression on 7.94% of paralogs and tissue-specific biases on 4.7% of paralogs. We found reduced expression of SMN-1 converted from SMN-2, which potentially affects diagnosis of spinal muscular atrophy, and increased expression of a duplicative translocation of AMY2B. Overall, ctyper enables biobank-scale genotyping of CNV and challenging genes.

Towards the genomic sequence code of DNA fragility for machine learning

Genomic DNA breakages and the subsequent insertion and deletion mutations are important contributors to genome instability and linked diseases. Unlike the research in point mutations, the relationship between DNA sequence context and the propensity for strand breaks remains elusive. Here, by analyzing the differences and commonalities across myriads of genomic breakage datasets, we extract the sequence-linked rules and patterns behind DNA fragility. We show the overall deconvolution of the sequence influence into short-, mid- and long-range effects, and the stressor-dependent differences in defining the range and compositional effects on DNA fragility. We summarize and release our feature compendium as a library that can be seamlessly incorporated into genomic machine learning procedures, where DNA fragility is of concern, and train a generalized DNA fragility model on cancer-associated breakages. Structural variants (SVs) tend to stabilize regions in which they emerge, with the effect most pronounced for pathogenic SVs. In contrast, the effects of chromothripsis are seen across regions less prone to breakages. We find that viral integration may bring genome fragility, particularly for cancer-associated viruses. Overall, this work offers novel insights into the genomic sequence basis of DNA fragility and presents a powerful machine learning resource to further enhance our understanding of genome (in)stability and evolution.

Genetic ancestry in population pharmacogenomics unravels distinct geographical patterns related to drug toxicity

Genetic ancestry plays a major role in pharmacogenomics, and a deeper understanding of the genetic diversity among individuals holds immerse promise for reshaping personalized medicine. In this pivotal study, we have conducted a large-scale genomic analysis of 1,136 pharmacogenomic variants employing machine learning algorithms on 3,714 individuals from publicly available datasets to assess the risk proximity of experiencing drug-related adverse events. Our findings indicate that Admixed Americans and Europeans have demonstrated a higher risk of experiencing drug toxicity, whereas individuals with East Asian ancestry and, to a lesser extent, Oceanians displayed a lower risk proximity. Polygenic risk scores for drug-gene interactions did not necessarily follow similar assumptions, reflecting distinct genetic patterns and population-specific differences that vary depending on the drug class. Overall, our results provide evidence that genetic ancestry is a pivotal factor in population pharmacogenomics and should be further exploited to strengthen even more personalized drug therapy.

Reconstruction of the human amylase locus reveals ancient duplications seeding modern-day variation

Previous studies suggested that the copy number of the human salivary amylase gene, AMY1, correlates with starch-rich diets. However, evolutionary analyses are hampered by the absence of accurate, sequence-resolved haplotype variation maps. We identified 30 structurally distinct haplotypes at nucleotide resolution among 98 present-day humans, revealing that the coding sequences of AMY1 copies are evolving under negative selection. Genomic analyses of these haplotypes in archaic hominins and ancient human genomes suggest that a common three-copy haplotype, dating as far back as 800 KYA, has seeded rapidly evolving rearrangements through recurrent non-allelic homologous recombination. Additionally, haplotypes with more than three AMY1 copies have significantly increased in frequency among European farmers over the past 4,000 years, potentially as an adaptive response to increased starch digestion.

Microsatellite density landscapes illustrate short tandem repeats aggregation in the complete reference human genome

BACKGROUND

Microsatellites are increasingly realized to have biological significance in human genome and health in past decades, the assembled complete reference sequence of human genome T2T-CHM13 brought great help for a comprehensive study of short tandem repeats in the human genome.

RESULTS

Microsatellites density landscapes of all 24 chromosomes were built here for the first complete reference sequence of human genome T2T-CHM13. These landscapes showed that short tandem repeats (STRs) are prone to aggregate characteristically to form a large number of STRs density peaks. We classified 8,823 High Microsatellites Density Peaks (HMDPs), 35,257 Middle Microsatellites Density Peaks (MMDPs) and 199, 649 Low Microsatellites Density Peaks (LMDPs) on the 24 chromosomes; and also classified the motif types of every microsatellites density peak. These STRs density aggregation peaks are mainly composing of a single motif, and AT is the most dominant motif, followed by AATGG and CCATT motifs. And 514 genomic regions were characterized by microsatellite density feature in the full T2T-CHM13 genome.

CONCLUSIONS

These landscape maps exhibited that microsatellites aggregate in many genomic positions to form a large number of microsatellite density peaks with composing of mainly single motif type in the complete reference genome, indicating that the local microsatellites density varies enormously along the every chromosome of T2T-CHM13.

Large-scale selection of highly informative microhaplotypes for ancestry inference and population specific informativeness

Microhaplotypes (MHs) describe physically close genetic markers that are inherited together and are gaining prominence due to their efficiency in forensic, clinical, and population studies. They excel in kinship analysis, DNA mixture detection, and ancestry inference, offering advantages in precision over individual SNPs and STRs. In this study, a pipeline was developed to efficiently select highly informative MHs from large-scale genomic datasets. Over 120,000 MHs were identified from almost a million markers, which allow this non-independent information to be efficiently used for inference. The MHs were compared to SNPs in terms of their informativeness and performance of their subsets in ancestry inference and all the results consistently favored MHs. A method for ranking markers by specific population informativeness was also introduced, which showed improvement in the accuracy of Native American ancestry estimation, overcoming the challenges of its underrepresentation in datasets. In conclusion, this study presents a comprehensive way for selecting highly informative MHs for accurate ancestry inference. The proposed approach and the subsets selected by specific population informativeness offer valuable tools for improving ancestry inference accuracy, particularly for admixed populations as demonstrated for a Brazilian dataset.

Machine Learning Reveals the Diversity of Human 3D Chromatin Contact Patterns

Understanding variation in chromatin contact patterns across diverse humans is critical for interpreting noncoding variants and their effects on gene expression and phenotypes. However, experimental determination of chromatin contact patterns across large samples is prohibitively expensive. To overcome this challenge, we develop and validate a machine learning method to quantify the variation in 3D chromatin contacts at 2 kilobase resolution from genome sequence alone. We apply this approach to thousands of human genomes from the 1000 Genomes Project and the inferred hominin ancestral genome. While patterns of 3D contact divergence genome wide are qualitatively similar to patterns of sequence divergence, we find substantial differences in 3D divergence and sequence divergence in local 1 megabase genomic windows. In particular, we identify 392 windows with significantly greater 3D divergence than expected from sequence. Moreover, for 31% of genomic windows, a single individual has a rare divergent 3D contact map pattern. Using in silico mutagenesis, we find that most single nucleotide sequence changes do not result in changes to 3D chromatin contacts. However, in windows with substantial 3D divergence just one or a few variants can lead to divergent 3D chromatin contacts without the individuals carrying those variants having high sequence divergence. In summary, inferring 3D chromatin contact maps across human populations reveals variable contact patterns. We anticipate that these genetically diverse maps of 3D chromatin contact will provide a reference for future work on the function and evolution of 3D chromatin contact variation across human populations.

The Human Accelerated Region HAR202 Controls NPAS3 Expression in the Developing Forebrain Displaying Differential Enhancer Activity Between Modern and Archaic Human Sequences

It has been proposed that the phenotypic differences in cognitive abilities between humans and our closest living relatives, chimpanzees, are largely due to changes in the regulation of neurodevelopmental genes. We have previously found that the neurodevelopmental transcription factor gene NPAS3 accumulates the largest number of human accelerated regions (HARs), suggesting it may play some role in the phenotypic evolution of the human nervous system. In this work, we performed a comparative functional analysis of NPAS3-HAR202 using enhancer reporter assays in transgenic zebrafish and mice. We found that the Homo sapiens HAR202 ortholog failed to drive reporter expression to the zebrafish nervous system, in high contrast to the strong expression displayed by the rest of the vertebrate ortholog sequences tested. Remarkably, the HAR202 ortholog from archaic humans (Neanderthals/Denisovans) also displayed a pan-vertebrate expression pattern, despite the fact that archaic and modern humans have only one nucleotide substitution. Moreover, similar results were found when comparing enhancer activity in transgenic mice, where we observed a loss of activity of the modern human version in the mouse developing brain. To investigate the functional importance of HAR202, we generated mice lacking HAR202 and found a remarkable decrease of Npas3 expression in the forebrain during development. Our results place HAR202 as one of the very few examples of a neurodevelopmental transcriptional enhancer displaying functional evolution in the brain as a result of a fast molecular evolutionary process that specifically occurred in the human lineage.

An IGHG1 variant exhibits polarized prevalence and confers enhanced IgG1 antibody responses against life-threatening organisms

Evolutionary pressures sculpt population genetics, whereas immune adaptation fortifies humans against life-threatening organisms. How the evolution of selective genetic variation in adaptive immune receptors orchestrates the adaptation of human populations to contextual perturbations remains elusive. Here, we show that the G396R coding variant within the human immunoglobulin G1 (IgG1) heavy chain presents a concentrated prevalence in Southeast Asian populations. We uncovered a 190-kb genomic linkage disequilibrium block peaked in close proximity to this variant, suggestive of potential Darwinian selection. This variant confers heightened immune resilience against various pathogens and viper toxins in mice. Mechanistic studies involving severe acute respiratory syndrome coronavirus 2 infection and vaccinated individuals reveal that this variant enhances pathogen-specific IgG1+ memory B cell activation and antibody production. This G396R variant may have arisen on a Neanderthal haplotype background. These findings underscore the importance of an IGHG1 variant in reinforcing IgG1 antibody responses against life-threatening organisms, unraveling the intricate interplay between human evolution and immune adaptation.

Recurrent evolution and selection shape structural diversity at the amylase locus

The adoption of agriculture triggered a rapid shift towards starch-rich diets in human populations1. Amylase genes facilitate starch digestion, and increased amylase copy number has been observed in some modern human populations with high-starch intake2, although evidence of recent selection is lacking3,4. Here, using 94 long-read haplotype-resolved assemblies and short-read data from approximately 5,600 contemporary and ancient humans, we resolve the diversity and evolutionary history of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in agricultural populations than in fishing, hunting and pastoral populations. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each underwent multiple duplication/deletion events with mutation rates up to more than 10,000-fold the single-nucleotide polymorphism mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome-based approach, we infer structural haplotypes across thousands of humans identifying extensively duplicated haplotypes at higher frequency in modern agricultural populations. Leveraging 533 ancient human genomes, we find that duplication-containing haplotypes (with more gene copies than the ancestral haplotype) have rapidly increased in frequency over the past 12,000 years in West Eurasians, suggestive of positive selection. Together, our study highlights the potential effects of the agricultural revolution on human genomes and the importance of structural variation in human adaptation.

Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

Genetic disease risks of under-represented founder populations in New York City

The detection of founder pathogenic variants, those observed in high frequency only in a group of individuals with increased inter-relatedness, can help improve delivery of health care for that community. We identified 16 groups with shared ancestry, based on genomic segments that are shared through identity by descent (IBD), in New York City using the genomic data of 25,366 residents from the All Of Us Research Program and the Mount Sinai BioMe biobank. From these groups we defined 8 as founder populations, mostly communities currently under-represented in medical genomics research, such as Puerto Rican, Garifuna and Filipino/Pacific Islanders. The enrichment analysis of ClinVar pathogenic or likely pathogenic (P/LP) variants in each group identified 202 of these damaging variants across the 8 founder populations. We confirmed disease-causing variants previously reported to occur at increased frequencies in Ashkenazi Jewish and Puerto Rican genetic ancestry groups, but most of the damaging variants identified have not been previously associated with any such founder populations, and most of these founder populations have not been described to have increased prevalence of the associated rare disease. Twenty-five of 51 variants meeting Tier 2 clinical screening criteria (1/100 carrier frequency within these founder groups) have never previously been reported. We show how population structure studies can provide insights into rare diseases disproportionately affecting under-represented founder populations, delivering a health care benefit but also a potential source of stigmatization of these communities, who should be part of the decision-making about implementation into health care delivery.

Gene expansions contributing to human brain evolution

Genomic drivers of human-specific neurological traits remain largely undiscovered. Duplicated genes expanded uniquely in the human lineage likely contributed to brain evolution, including the increased complexity of synaptic connections between neurons and the dramatic expansion of the neocortex. Discovering duplicate genes is challenging because the similarity of paralogs makes them prone to sequence-assembly errors. To mitigate this issue, we analyzed a complete telomere-to-telomere human genome sequence (T2T-CHM13) and identified 213 duplicated gene families likely containing human-specific paralogs (>98% identity). Positing that genes important in universal human brain features should exist with at least one copy in all modern humans and exhibit expression in the brain, we narrowed in on 362 paralogs with at least one copy across thousands of ancestrally diverse genomes and present in human brain transcriptomes. Of these, 38 paralogs co-express in gene modules enriched for autism-associated genes and potentially contribute to human language and cognition. We narrowed in on 13 duplicate gene families with human-specific paralogs that are fixed among modern humans and show convincing brain expression patterns. Using long-read DNA sequencing revealed hidden variation across 200 modern humans of diverse ancestries, uncovering signatures of selection not previously identified, including possible balancing selection of CD8B. To understand the roles of duplicated genes in brain development, we generated zebrafish CRISPR "knockout" models of nine orthologs and transiently introduced mRNA-encoding paralogs, effectively "humanizing" the larvae. Morphometric, behavioral, and single-cell RNA-seq screening highlighted, for the first time, a possible role for GPR89B in dosage-mediated brain expansion and FRMPD2B function in altered synaptic signaling, both hallmark features of the human brain. Our holistic approach provides important insights into human brain evolution as well as a resource to the community for studying additional gene expansion drivers of human brain evolution.

A consensus-based classification workflow to determine genetically inferred ancestry from comprehensive genomic profiling of patients with solid tumors

Disparities in cancer diagnosis, treatment, and outcomes based on self-identified race and ethnicity (SIRE) are well documented, yet these variables have historically been excluded from clinical research. Without SIRE, genetic ancestry can be inferred using single-nucleotide polymorphisms (SNPs) detected from tumor DNA using comprehensive genomic profiling (CGP). However, factors inherent to CGP of tumor DNA increase the difficulty of identifying ancestry-informative SNPs, and current workflows for inferring genetic ancestry from CGP need improvements in key areas of the ancestry inference process. This study used genomic data from 4274 diverse reference subjects and CGP data from 491 patients with solid tumors and SIRE to develop and validate a workflow to obtain accurate genetically inferred ancestry (GIA) from CGP sequencing results. We use consensus-based classification to derive confident ancestral inferences from an expanded reference dataset covering eight world populations (African, Admixed American, Central Asian/Siberian, European, East Asian, Middle Eastern, Oceania, South Asian). Our GIA calls were highly concordant with SIRE (95%) and aligned well with reference populations of inferred ancestries. Further, our workflow could expand on SIRE by (i) detecting the ancestry of patients that usually lack appropriate racial categories, (ii) determining what patients have mixed ancestry, and (iii) resolving ancestries of patients in heterogeneous racial categories and who had missing SIRE. Accurate GIA provides needed information to enable ancestry-aware biomarker research, ensure the inclusion of underrepresented groups in clinical research, and increase the diverse representation of patient populations eligible for precision medicine therapies and trials.

Pervasive findings of directional selection realize the promise of ancient DNA to elucidate human adaptation

We present a method for detecting evidence of natural selection in ancient DNA time-series data that leverages an opportunity not utilized in previous scans: testing for a consistent trend in allele frequency change over time. By applying this to 8433 West Eurasians who lived over the past 14000 years and 6510 contemporary people, we find an order of magnitude more genome-wide significant signals than previous studies: 347 independent loci with >99% probability of selection. Previous work showed that classic hard sweeps driving advantageous mutations to fixation have been rare over the broad span of human evolution, but in the last ten millennia, many hundreds of alleles have been affected by strong directional selection. Discoveries include an increase from ~0% to ~20% in 4000 years for the major risk factor for celiac disease at HLA-DQB1; a rise from ~0% to ~8% in 6000 years of blood type B; and fluctuating selection at the TYK2 tuberculosis risk allele rising from ~2% to ~9% from ~5500 to ~3000 years ago before dropping to ~3%. We identify instances of coordinated selection on alleles affecting the same trait, with the polygenic score today predictive of body fat percentage decreasing by around a standard deviation over ten millennia, consistent with the "Thrifty Gene" hypothesis that a genetic predisposition to store energy during food scarcity became disadvantageous after farming. We also identify selection for combinations of alleles that are today associated with lighter skin color, lower risk for schizophrenia and bipolar disease, slower health decline, and increased measures related to cognitive performance (scores on intelligence tests, household income, and years of schooling). These traits are measured in modern industrialized societies, so what phenotypes were adaptive in the past is unclear. We estimate selection coefficients at 9.9 million variants, enabling study of how Darwinian forces couple to allelic effects and shape the genetic architecture of complex traits.

Long genetic and social isolation in Neanderthals before their extinction

Neanderthal genomes have been recovered from sites across Eurasia, painting an increasingly complex picture of their populations' structure that mostly indicates that late European Neanderthals belonged to a single metapopulation with no significant evidence of population structure. Here, we report the discovery of a late Neanderthal individual, nicknamed "Thorin," from Grotte Mandrin in Mediterranean France, and his genome. These dentognathic fossils, including a rare example of distomolars, are associated with a rich archeological record of Neanderthal final technological traditions in this region ∼50-42 thousand years ago. Thorin's genome reveals a relatively early divergence of ∼105 ka with other late Neanderthals. Thorin belonged to a population with a small group size that showed no genetic introgression with other known late European Neanderthals, revealing some 50 ka of genetic isolation of his lineage despite them living in neighboring regions. These results have important implications for resolving competing hypotheses about causes of the disappearance of the Neanderthals.

The Genomic and Cultural Diversity of the Inka Qhapaq Hucha Ceremony in Chile and Argentina

The South American archaeological record has ample evidence of the socio-cultural dynamism of human populations in the past. This has also been supported through the analysis of ancient genomes, by showing evidence of gene flow across the region. While the extent of these signals is yet to be tested, the growing number of ancient genomes allows for more fine-scaled hypotheses to be evaluated. In this study, we assessed the genetic diversity of individuals associated with the Inka ritual, Qhapaq hucha. As part of this ceremony, one or more individuals were buried with Inka and local-style offerings on mountain summits along the Andes, leaving a very distinctive record. Using paleogenomic tools, we analyzed three individuals: two newly generated genomes from El Plomo Mountain (Chile) and El Toro Mountain (Argentina), and a previously published genome from Argentina (Aconcagua Mountain). Our results reveal a complex demographic scenario with each of the individuals showing different genetic affinities. Furthermore, while two individuals showed genetic similarities with present-day and ancient populations from the southern region of the Inka empire, the third individual may have undertaken long-distance movement. The genetic diversity we observed between individuals from similar cultural contexts supports the highly diverse strategies Inka implemented while incorporating new territories. More broadly, this research contributes to our growing understanding of the population dynamics in the Andes by discussing the implications and temporality of population movements in the region.

Ancient Rapanui genomes reveal resilience and pre-European contact with the Americas

Rapa Nui (also known as Easter Island) is one of the most isolated inhabited places in the world. It has captured the imagination of many owing to its archaeological record, which includes iconic megalithic statues called moai1. Two prominent contentions have arisen from the extensive study of Rapa Nui. First, the history of the Rapanui has been presented as a warning tale of resource overexploitation that would have culminated in a major population collapse-the 'ecocide' theory2-4. Second, the possibility of trans-Pacific voyages to the Americas pre-dating European contact is still debated5-7. Here, to address these questions, we reconstructed the genomic history of the Rapanui on the basis of 15 ancient Rapanui individuals that we radiocarbon dated (1670-1950 CE) and whole-genome sequenced (0.4-25.6×). We find that these individuals are Polynesian in origin and most closely related to present-day Rapanui, a finding that will contribute to repatriation efforts. Through effective population size reconstructions and extensive population genetics simulations, we reject a scenario involving a severe population bottleneck during the 1600s, as proposed by the ecocide theory. Furthermore, the ancient and present-day Rapanui carry similar proportions of Native American admixture (about 10%). Using a Bayesian approach integrating genetic and radiocarbon dates, we estimate that this admixture event occurred about 1250-1430 CE.

Short Tandem Repeats in the era of next-generation sequencing: from historical loci to population databases

In this study, we explore the landscape of short tandem repeats (STRs) within the human genome through the lens of evolving technologies to detect genomic variations. STRs, which encompass approximately 3% of our genomic DNA, are crucial for understanding human genetic diversity, disease mechanisms, and evolutionary biology. The advent of high-throughput sequencing methods has revolutionized our ability to accurately map and analyze STRs, highlighting their significance in genetic disorders, forensic science, and population genetics. We review the current available methodologies for STR analysis, the challenges in interpreting STR variations across different populations, and the implications of STRs in medical genetics. Our findings underscore the urgent need for comprehensive STR databases that reflect the genetic diversity of global populations, facilitating the interpretation of STR data in clinical diagnostics, genetic research, and forensic applications. This work sets the stage for future studies aimed at harnessing STR variations to elucidate complex genetic traits and diseases, reinforcing the importance of integrating STRs into genetic research and clinical practice.

Association of genetic ancestry with molecular tumor profiles in colorectal cancer

BACKGROUND

There are known disparities in incidence and outcomes of colorectal cancer (CRC) by race and ethnicity. Some of these disparities may be mediated by molecular changes in tumors that occur at different rates across populations. Genetic ancestry is a measure complementary to race and ethnicity that can overcome missing data issues and better capture genetic similarity in admixed populations. We aimed to identify somatic mutations and tumor gene expression differences associated with both genetic ancestry and imputed race and ethnicity.

METHODS

Sequencing was performed with the Tempus xT NGS 648-gene panel and whole exome capture RNA-Seq for 8454 primarily late-stage CRC patients. Genetic ancestry proportions for five continental groups-Africa (AFR), American indigenous (AMR), East Asia (EAS), Europe (EUR), and South Asia (SAS)-were estimated using ancestry informative markers. To address data gaps, race and ethnicity categories were imputed, resulting in assignments for 952 Hispanic/Latino, 420 non-Hispanic (NH) Asian, 1061 NH Black, and 5763 NH White individuals. We assessed association of genetic ancestry proportions and imputed race and ethnicity categories with somatic mutations in relevant CRC genes and in 2608 expression profiles, as well as 1957 consensus molecular subtypes (CMS).

RESULTS

Increased AFR ancestry was associated with higher odds of somatic mutations in APC, KRAS, and PIK3CA and lower odds of BRAF mutations. Additionally, increased EAS ancestry was associated with lower odds of mutations in KRAS, EUR with higher odds in BRAF, and the Hispanic/Latino category with lower odds in BRAF. Greater AFR ancestry and the NH Black category were associated with higher rates of CMS3, while a higher proportion of Hispanic/Latino patients exhibited indeterminate CMS classifications.

CONCLUSIONS

Molecular differences in CRC tumor mutation frequencies and gene expression that may underlie observed differences by race and ethnicity were identified. The association of AFR ancestry with increased KRAS mutations aligns with higher CMS3 subtype rates in NH Black patients. The increase of indeterminate CMS in Hispanic/Latino patients suggests that subtype classification methods could benefit from enhanced patient diversity.

Frontal sinus size in South African Later Stone Age Holocene Khoe-San

Frontal size variation is comparatively poorly sampled among sub-Saharan African populations. This study assessed frontal sinus size in a sample of Khoe-San skeletal remains from South African Later Stone Age contexts. Volumes were determined from CT scans of 102 adult crania; individual sex could be estimated in 82 cases. Sinus volume is not sexually dimorphic in this sample. The lack of frontal sinus aplasia is concordant with the low incidences recorded for other sub-Saharan African and most other global populations save those that inhabit high latitudes. There is considerable variation in frontal sinus size among global populations, and the Khoe-San possess among the smallest. The Khoe-San have rather diminutive sinuses compared to sub-Saharan Bantu-speaking populations but resemble a northern African (Sudanese) population. Genetic studies indicate the earliest population divergence within Homo sapiens to have been between the Khoe-San and all other living groups, and that this likely occurred in Africa during the span of Marine Isotope Stages 8-6. There is scant information on frontal sinus development among Late Quaternary African fossils that are likely either closely related or attributable to Homo sapiens. Among these, the MIS 3 cranium from Hofmeyr, South Africa, exhibits distinct Khoe-San cranial affinities and despite its large size has a very small frontal sinus. This raises the possibility that the small frontal sinuses of the Holocene South African Khoe-San might be a feature retained from an earlier MIS 3 population.

Modelling the demographic history of human North African genomes points to a recent soft split divergence between populations

BACKGROUND

North African human populations present a complex demographic scenario due to the presence of an autochthonous genetic component and population substructure, plus extensive gene flow from the Middle East, Europe, and sub-Saharan Africa.

RESULTS

We conducted a comprehensive analysis of 364 genomes to construct detailed demographic models for the North African region, encompassing its two primary ethnic groups, the Arab and Amazigh populations. This was achieved through an Approximate Bayesian Computation with Deep Learning (ABC-DL) framework and a novel algorithm called Genetic Programming for Population Genetics (GP4PG). This innovative approach enabled us to effectively model intricate demographic scenarios, utilizing a subset of 16 whole genomes at > 30X coverage. The demographic model suggested by GP4PG exhibited a closer alignment with the observed data compared to the ABC-DL model. Both point to a back-to-Africa origin of North African individuals and a close relationship with Eurasian populations. Results support different origins for Amazigh and Arab populations, with Amazigh populations originating back in Epipaleolithic times, while GP4PG supports Arabization as the main source of Middle Eastern ancestry. The GP4PG model includes population substructure in surrounding populations (sub-Saharan Africa and Middle East) with continuous decaying gene flow after population split. Contrary to ABC-DL, the best GP4PG model does not require pulses of admixture from surrounding populations into North Africa pointing to soft splits as drivers of divergence in North Africa.

CONCLUSIONS

We have built a demographic model on North Africa that points to a back-to-Africa expansion and a differential origin between Arab and Amazigh populations.

Geographic EBV variants confound disease-specific variant interpretation and predict variable immune therapy responses

ABSTRACT

Epstein-Barr virus (EBV) is a potent carcinogen linked to hematologic and solid malignancies and causes significant global morbidity and mortality. Therapy using allogeneic EBV-specific lymphocytes shows promise in certain populations, but the impact of EBV genome variation on these strategies remains unexplored. To address this, we sequenced 217 EBV genomes, including hematologic malignancies from Guatemala, Peru, Malawi, and Taiwan, and analyzed them alongside 1307 publicly available EBV genomes from cancer, nonmalignant diseases, and healthy individuals across Africa, Asia, Europe, North America, and South America. These included, to our knowledge, the first natural killer (NK)/T-cell lymphoma (NKTCL) EBV genomes reported outside of East Asia. Our findings indicate that previously proposed EBV genome variants specific to certain cancer types are more closely tied to geographic origin than to cancer histology. This included variants previously reported to be specific to NKTCL but were prevalent in EBV genomes from other cancer types and healthy individuals in East Asia. After controlling for geographic region, we did identify multiple NKTCL-specific variants associated with a 7.8-fold to 21.9-fold increased risk. We also observed frequent variations in EBV genomes that affected peptide sequences previously reported to bind common major histocompatibility complex alleles. Finally, we found several nonsynonymous variants spanning the coding sequences of current vaccine targets BALF4, BKRF2, BLLF1, BXLF2, BZLF1, and BZLF2. These results highlight the need to consider geographic variation in EBV genomes when devising strategies for exploiting adaptive immune responses against EBV-related cancers, ensuring greater global effectiveness and equity in prevention and treatment.

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.

The impact of clinical genome sequencing in a global population with suspected rare genetic disease

There is mounting evidence of the value of clinical genome sequencing (cGS) in individuals with suspected rare genetic disease (RGD), but cGS performance and impact on clinical care in a diverse population drawn from both high-income countries (HICs) and low- and middle-income countries (LMICs) has not been investigated. The iHope program, a philanthropic cGS initiative, established a network of 24 clinical sites in eight countries through which it provided cGS to individuals with signs or symptoms of an RGD and constrained access to molecular testing. A total of 1,004 individuals (median age, 6.5 years; 53.5% male) with diverse ancestral backgrounds (51.8% non-majority European) were assessed from June 2016 to September 2021. The diagnostic yield of cGS was 41.4% (416/1,004), with individuals from LMIC sites 1.7 times more likely to receive a positive test result compared to HIC sites (LMIC 56.5% [195/345] vs. HIC 33.5% [221/659], OR 2.6, 95% CI 1.9-3.4, p < 0.0001). A change in diagnostic evaluation occurred in 76.9% (514/668) of individuals. Change of management, inclusive of specialty referrals, imaging and testing, therapeutic interventions, and palliative care, was reported in 41.4% (285/694) of individuals, which increased to 69.2% (480/694) when genetic counseling and avoidance of additional testing were also included. Individuals from LMIC sites were as likely as their HIC counterparts to experience a change in diagnostic evaluation (OR 6.1, 95% CI 1.1-∞, p = 0.05) and change of management (OR 0.9, 95% CI 0.5-1.3, p = 0.49). Increased access to genomic testing may support diagnostic equity and the reduction of global health care disparities.

An assessment of the genomic structural variation landscape in Sub-Saharan African populations

Structural variants are responsible for a large part of genomic variation between individuals and play a role in both common and rare diseases. Databases cataloguing structural variants notably do not represent the full spectrum of global diversity, particularly missing information from most African populations. To address this representation gap, we analysed 1,091 high-coverage African genomes, 545 of which are public data sets, and 546 which have been analysed for structural variants for the first time. Variants were called using five different tools and datasets merged and jointly called using SURVIVOR. We identified 67,795 structural variants throughout the genome, with 10,421 genes having at least one variant. Using a conservative overlap in merged data, 6,414 of the structural variants (9.5%) are novel compared to the Database of Genomic Variants. This study contributes to knowledge of the landscape of structural variant diversity in Africa and presents a reliable dataset for potential applications in population genetics and health-related research.

Forensic investigative genetic genealogy: expanding pedigree tracing and genetic inquiry in the genomic era

Genetic genealogy provides crucial insights into the complex biological relationships within contemporary and ancient human populations by analyzing shared alleles and chromosomal segments that are identical by descent to understand kinship, migration patterns, and population dynamics. Within forensic science, forensic investigative genetic genealogy (FIGG) has gained prominence by leveraging next-generation sequencing technologies and population-specific genomic resources, opening new investigative avenues. In this review, we synthesize current knowledge, underscore recent advancements, and discuss the growing role of FIGG in forensic genomics. FIGG has been pivotal in revitalizing dormant inquiries and offering new genetic leads in numerous cold cases. Its effectiveness relies on the extensive single-nucleotide polymorphism profiles contributed by individuals from diverse populations to specialized genomic databases. Advances in computational genomics and the growth of human genomic databases have spurred a profound shift in the application of genetic genealogy across forensics, anthropology, and ancient DNA studies. As the field progresses, FIGG is evolving from a nascent practice into a more sophisticated and specialized discipline, shaping the future of forensic investigations.

Multiple Human Population Movements and Cultural Dispersal Events Shaped the Landscape of Chinese Paternal Heritage

Large-scale genomic projects and ancient DNA innovations have ushered in a new paradigm for exploring human evolutionary history. However, the genetic legacy of spatiotemporally diverse ancient Eurasians within Chinese paternal lineages remains unresolved. Here, we report an integrated Y-chromosome genomic database encompassing 15,563 individuals from both modern and ancient Eurasians, including 919 newly reported individuals, to investigate the Chinese paternal genomic diversity. The high-resolution, time-stamped phylogeny reveals multiple diversification events and extensive expansions in the early and middle Neolithic. We identify four major ancient population movements, each associated with technological innovations that have shaped the Chinese paternal landscape. First, the expansion of early East Asians and millet farmers from the Yellow River Basin predominantly carrying O2/D subclades significantly influenced the formation of the Sino-Tibetan people and facilitated the permanent settlement of the Tibetan Plateau. Second, the dispersal of rice farmers from the Yangtze River Valley carrying O1 and certain O2 sublineages reshapes the genetic makeup of southern Han Chinese, as well as the Tai-Kadai, Austronesian, Hmong-Mien, and Austroasiatic people. Third, the Neolithic Siberian Q/C paternal lineages originated and proliferated among hunter-gatherers on the Mongolian Plateau and the Amur River Basin, leaving a significant imprint on the gene pools of northern China. Fourth, the J/G/R paternal lineages derived from western Eurasia, which were initially spread by Yamnaya-related steppe pastoralists, maintain their presence primarily in northwestern China. Overall, our research provides comprehensive genetic evidence elucidating the significant impact of interactions with culturally distinct ancient Eurasians on the patterns of paternal diversity in modern Chinese populations.

A harmonized public resource of deeply sequenced diverse human genomes

Underrepresented populations are often excluded from genomic studies owing in part to a lack of resources supporting their analyses. The 1000 Genomes Project (1kGP) and Human Genome Diversity Project (HGDP), which have recently been sequenced to high coverage, are valuable genomic resources because of the global diversity they capture and their open data sharing policies. Here, we harmonized a high-quality set of 4094 whole genomes from 80 populations in the HGDP and 1kGP with data from the Genome Aggregation Database (gnomAD) and identified over 153 million high-quality SNVs, indels, and SVs. We performed a detailed ancestry analysis of this cohort, characterizing population structure and patterns of admixture across populations, analyzing site frequency spectra, and measuring variant counts at global and subcontinental levels. We also show substantial added value from this data set compared with the prior versions of the component resources, typically combined via liftOver and variant intersection; for example, we catalog millions of new genetic variants, mostly rare, compared with previous releases. In addition to unrestricted individual-level public release, we provide detailed tutorials for conducting many of the most common quality-control steps and analyses with these data in a scalable cloud-computing environment and publicly release this new phased joint callset for use as a haplotype resource in phasing and imputation pipelines. This jointly called reference panel will serve as a key resource to support research of diverse ancestry populations.

Human Y chromosome haplogroup L1-M22 traces Neolithic expansion in West Asia and supports the Elamite and Dravidian connection

West and South Asian populations profoundly influenced Eurasian genetic and cultural diversity. We investigate the genetic history of the Y chromosome haplogroup L1-M22, which, while prevalent in these regions, lacks in-depth study. Robust Bayesian analyses of 165 high-coverage Y chromosomes favor a West Asian origin for L1-M22 ∼20.6 thousand years ago (kya). Moreover, this haplogroup parallels the genome-wide genetic ancestry of hunter-gatherers from the Iranian Plateau and the Caucasus. We characterized two L1-M22 harboring population groups during the Early Holocene. One expanded with the West Asian Neolithic transition. The other moved to South Asia ∼8-6 kya but showed no expansion. This group likely participated in the spread of Dravidian languages. These South Asian L1-M22 lineages expanded ∼4-3 kya, coinciding with the Steppe ancestry introduction. Our findings advance the current understanding of Eurasian historical dynamics, emphasizing L1-M22's West Asian origin, associated population movements, and possible linguistic impacts.

Ancient genomes revealed the complex human interactions of the ancient western Tibetans

The western Tibetan Plateau is the crossroad between the Tibetan Plateau, Central Asia, and South Asia, and it is a potential human migration pathway connecting these regions. However, the population history of the western Tibetan Plateau remains largely unexplored due to the lack of ancient genomes covering a long-time interval from this area. Here, we reported genome-wide data of 65 individuals dated to 3,500-300 years before present (BP) in the Ngari prefecture. The ancient western Tibetan Plateau populations share the majority of their genetic components with the southern Tibetan Plateau populations and have maintained genetic continuity since 3,500 BP while maintaining interactions with populations within and outside the Tibetan Plateau. Within the Tibetan Plateau, the ancient western Tibetan Plateau populations were influenced by the additional expansion from the south to the southwest plateau before 1,800 BP. Outside the Tibetan Plateau, the western Tibetan Plateau populations interacted with both South and Central Asian populations at least 2,000 years ago, and the South Asian-related genetic influence, despite being very limited, was from the Indus Valley Civilization (IVC) migrants in Central Asia instead of the IVC populations from the Indus Valley. In light of the new genetic data, our study revealed the complex population interconnections across and within the Tibetan Plateau.

Global diversity, recurrent evolution, and recent selection on amylase structural haplotypes in humans

The adoption of agriculture, first documented ~12,000 years ago in the Fertile Crescent, triggered a rapid shift toward starch-rich diets in human populations. Amylase genes facilitate starch digestion and increased salivary amylase copy number has been observed in some modern human populations with high starch intake, though evidence of recent selection is lacking. Here, using 52 long-read diploid assemblies and short read data from ~5,600 contemporary and ancient humans, we resolve the diversity, evolutionary history, and selective impact of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in populations with agricultural subsistence compared to fishing, hunting, and pastoral groups. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each exhibit multiple duplications/deletions with mutation rates >10,000-fold the SNP mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome graph-based approach to infer structural haplotypes across thousands of humans, we identify extensively duplicated haplotypes present at higher frequencies in modern day populations with traditionally agricultural diets. Leveraging 533 ancient human genomes we find that duplication-containing haplotypes (i.e. haplotypes with more amylase gene copies than the ancestral haplotype) have increased in frequency more than seven-fold over the last 12,000 years providing evidence for recent selection in West Eurasians. Together, our study highlights the potential impacts of the agricultural revolution on human genomes and the importance of long-read sequencing in identifying signatures of selection at structurally complex loci.

A disease-associated gene desert directs macrophage inflammation through ETS2

Increasing rates of autoimmune and inflammatory disease present a burgeoning threat to human health1. This is compounded by the limited efficacy of available treatments1 and high failure rates during drug development2, highlighting an urgent need to better understand disease mechanisms. Here we show how functional genomics could address this challenge. By investigating an intergenic haplotype on chr21q22-which has been independently linked to inflammatory bowel disease, ankylosing spondylitis, primary sclerosing cholangitis and Takayasu's arteritis3-6-we identify that the causal gene, ETS2, is a central regulator of human inflammatory macrophages and delineate the shared disease mechanism that amplifies ETS2 expression. Genes regulated by ETS2 were prominently expressed in diseased tissues and more enriched for inflammatory bowel disease GWAS hits than most previously described pathways. Overexpressing ETS2 in resting macrophages reproduced the inflammatory state observed in chr21q22-associated diseases, with upregulation of multiple drug targets, including TNF and IL-23. Using a database of cellular signatures7, we identified drugs that might modulate this pathway and validated the potent anti-inflammatory activity of one class of small molecules in vitro and ex vivo. Together, this illustrates the power of functional genomics, applied directly in primary human cells, to identify immune-mediated disease mechanisms and potential therapeutic opportunities.

Ancient genomes reveal insights into ritual life at Chichén Itzá

The ancient city of Chichén Itzá in Yucatán, Mexico, was one of the largest and most influential Maya settlements during the Late and Terminal Classic periods (AD 600-1000) and it remains one of the most intensively studied archaeological sites in Mesoamerica1-4. However, many questions about the social and cultural use of its ceremonial spaces, as well as its population's genetic ties to other Mesoamerican groups, remain unanswered2. Here we present genome-wide data obtained from 64 subadult individuals dating to around AD 500-900 that were found in a subterranean mass burial near the Sacred Cenote (sinkhole) in the ceremonial centre of Chichén Itzá. Genetic analyses showed that all analysed individuals were male and several individuals were closely related, including two pairs of monozygotic twins. Twins feature prominently in Mayan and broader Mesoamerican mythology, where they embody qualities of duality among deities and heroes5, but until now they had not been identified in ancient Mayan mortuary contexts. Genetic comparison to present-day people in the region shows genetic continuity with the ancient inhabitants of Chichén Itzá, except at certain genetic loci related to human immunity, including the human leukocyte antigen complex, suggesting signals of adaptation due to infectious diseases introduced to the region during the colonial period.

Evolutionary patterns of DNA base composition at polymorphic sites highlight the role of the environment in shaping barley and rice genomes

Insights into changes in genome base composition underlying crop domestication can be gained by using comparative genomics. With this approach, previous studies have reported that crop genomes during domestication accumulate more nucleotides adenine (A) and thymine (T) (termed as [AT]-increase) across polymorphic sites. However, the potential influence of the environment or its factors, for example, solar ultraviolet (UV) radiation and temperature, on the [AT]-increase has not been well elucidated. Here, we investigated the [AT]-increase in barley (Hordeum vulgare L.) and rice (Oryza sativa L.) and the association with natural environments, where accessions are distributed. With 12,798,376 and 2,861,535 single-nucleotide polymorphisms from 368 barley and 1375 rice accessions, respectively, we discovered that [AT] increases from wild accessions to improved cultivars, and genomic regions with larger [AT]-increase tend to have higher UV-related motif frequencies, suggesting solar UV radiation as a potential factor in driving genome variation. To link [AT] change with the geographic distribution, we gathered georeferenced accessions and examined their local environments. Interestingly, negative correlations between [AT] and environmental factors were observed (r = -0.39 ∼ -0.75) and modern accessions with higher [AT] values, as compared with wild relatives, are from the environments with lower solar UV radiation or lower temperature. With [AT] and environmental factors as phenotypes, genome-wide association mapping identified three candidate genes that have the potential to contribute to [AT] variation under the effect of environmental conditions. Our findings provide genomic and environmental insights into evolutionary pattern of DNA base composition and underlying mechanisms.

A survey of BWT variants for string collections

MOTIVATION

In recent years, the focus of bioinformatics research has moved from individual sequences to collections of sequences. Given the fundamental role of the Burrows-Wheeler Transform (BWT) in string processing, a number of dedicated tools have been developed for computing the BWT of string collections. While the focus has been on improving efficiency, both in space and time, the exact definition of the BWT employed has not been at the center of attention. As we show in this paper, the different tools in use often compute non-equivalent BWT variants: the resulting transforms can differ from each other significantly, including the number r of runs, a central parameter of the BWT. Moreover, with many tools, the transform depends on the input order of the collection. In other words, on the same dataset, the same tool may output different transforms if the dataset is given in a different order.

RESULTS

We studied 18 dedicated tools for computing the BWT of string collections and were able to identify 6 different BWT variants computed by these tools. We review the differences between these BWT variants, both from a theoretical and from a practical point of view, comparing them on 8 real-life biological datasets with different characteristics. We find that the differences can be extensive, depending on the datasets, and are largest on collections of many similar short sequences. The parameter r, the number of runs of the BWT, also shows notable variation between the different BWT variants; on our datasets, it varied by a multiplicative factor of up to 4.2.

AVAILABILITY

Source code and scripts to replicate the results and download the data used in the article are available at https://github.com/davidecenzato/BWT-variants-for-string-collections.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Abrupt onset of intensive human occupation 44,000 years ago on the threshold of Sahul

Archaeological evidence attests multiple early dispersals of Homo sapiens out of Africa, but genetic evidence points to the primacy of a single dispersal 70-40 ka. Laili in Timor-Leste is on the southern dispersal route between Eurasia and Australasia and has the earliest record of human occupation in the eastern Wallacean archipelago. New evidence from the site shows that, unusually in the region, sediment accumulated in the shelter without human occupation, in the window 59-54 ka. This was followed by an abrupt onset of intensive human habitation beginning ~44 ka. The initial occupation is distinctive from overlying layers in the aquatic focus of faunal exploitation, while it has similarities in material culture to other early Homo sapiens sites in Wallacea. We suggest that the intensive early occupation at Laili represents a colonisation phase, which may have overwhelmed previous human dispersals in this part of the world.

Neandertal ancestry through time: Insights from genomes of ancient and present-day humans

Gene flow from Neandertals has shaped the landscape of genetic and phenotypic variation in modern humans. We identify the location and size of introgressed Neandertal ancestry segments in more than 300 genomes spanning the last 50,000 years. We study how Neandertal ancestry is shared among individuals to infer the time and duration of the Neandertal gene flow. We find the correlation of Neandertal segment locations across individuals and their divergence to sequenced Neandertals, both support a model of single major Neandertal gene flow. Our catalog of introgressed segments through time confirms that most natural selection-positive and negative-on Neandertal ancestry variants occurred immediately after the gene flow, and provides new insights into how the contact with Neandertals shaped human origins and adaptation.

Comprehensive landscape of non-CODIS STRs in global populations provides new insights into challenging DNA profiles

The worldwide implementation of short tandem repeats (STR) profiles in forensic genetics necessitated establishing and expanding the CODIS core loci set to facilitated efficient data management and exchange. Currently, the mainstay CODIS STRs are adopted in most general-purpose forensic kits. However, relying solely on these loci failed to yield satisfactory results for challenging tasks, such as bio-geographical ancestry inference, complex DNA mixture profile interpretation, and distant kinship analysis. In this context, non-CODIS STRs are potent supplements to enhance the systematic discriminating power, particularly when combined with the high-throughput next-generation sequencing (NGS) technique. Nevertheless, comprehensive evaluation on non-CODIS STRs in diverse populations was scarce, hindering their further application in routine caseworks. To address this gap, we investigated genetic variations of 178 historically available non-CODIS STRs from ethnolinguistically different worldwide populations and studied their characteristics and forensic potentials via high-coverage whole genome sequencing (WGS) data. Initially, we delineated the genomic properties of these non-CODIS markers through sequence searching, repeat structure scanning, and manual inspection. Subsequent population genetics analysis suggested that these non-CODIS STRs had comparable polymorphism levels and forensic utility to CODIS STRs. Furthermore, we constructed a theoretical next-generation sequencing (NGS) panel comprising 108 STRs (20 CODIS STRs and 88 non-CODIS STRs), and evaluated its performance in inferring bio-geographical ancestry origins, deconvoluting complex DNA mixtures, and differentiating distant kinships using real and simulated datasets. Our findings demonstrated that incorporating supplementary non-CODIS STRs enabled the extrapolation of multidimensional information from a single STR profile, thereby facilitating the analysis of challenging forensic tasks. In conclusion, this study presents an extensive genomic landscape of forensic non-CODIS STRs among global populations, and emphasized the imperative inclusion of additional polymorphic non-CODIS STRs in future NGS-based forensic systems.

Understanding the genomic heterogeneity of North African Imazighen: from broad to microgeographical perspectives

The strategic location of North Africa has led to cultural and demographic shifts, shaping its genetic structure. Historical migrations brought different genetic components that are evident in present-day North African genomes, along with autochthonous components. The Imazighen (plural of Amazigh) are believed to be the descendants of autochthonous North Africans and speak various Amazigh languages, which belong to the Afro-Asiatic language family. However, the arrival of different human groups, especially during the Arab conquest, caused cultural and linguistic changes in local populations, increasing their heterogeneity. We aim to characterize the genetic structure of the region, using the largest Amazigh dataset to date and other reference samples. Our findings indicate microgeographical genetic heterogeneity among Amazigh populations, modeled by various admixture waves and different effective population sizes. A first admixture wave is detected group-wide around the twelfth century, whereas a second wave appears in some Amazigh groups around the nineteenth century. These events involved populations with higher genetic ancestry from south of the Sahara compared to the current North Africans. A plausible explanation would be the historical trans-Saharan slave trade, which lasted from the Roman times to the nineteenth century. Furthermore, our investigation shows that assortative mating in North Africa has been rare.

Polymorphic transposable elements provide new insights on high-altitude adaptation in the Tibetan Plateau

Several studies demonstrated that populations living in the Tibetan plateau are genetically and physiologically adapted to high-altitude conditions, showing genomic signatures ascribable to the action of natural selection. However, so far most of them relied solely on inferences drawn from the analysis of coding variants and point mutations. To fill this gap, we focused on the possible role of polymorphic transposable elements in influencing the adaptation of Tibetan and Sherpa highlanders. To do so, we compared high-altitude and middle/low-lander individuals of East Asian ancestry by performing in silico analyses and differentiation tests on 118 modern and ancient samples. We detected several transposable elements associated with high altitude, which map genes involved in cardiovascular, hematological, chem-dependent and respiratory conditions, suggesting that metabolic and signaling pathways taking part in these functions are disproportionately impacted by the effect of environmental stressors in high-altitude individuals. To our knowledge, our study is the first hinting to a possible role of transposable elements in the adaptation of Tibetan and Sherpa highlanders.

Positive selection in the genomes of two Papua New Guinean populations at distinct altitude levels

Highlanders and lowlanders of Papua New Guinea have faced distinct environmental stress, such as hypoxia and environment-specific pathogen exposure, respectively. In this study, we explored the top genomics regions and the candidate driver SNPs for selection in these two populations using newly sequenced whole-genomes of 54 highlanders and 74 lowlanders. We identified two candidate SNPs under selection - one in highlanders, associated with red blood cell traits and another in lowlanders, which is associated with white blood cell count - both potentially influencing the heart rate of Papua New Guineans in opposite directions. We also observed four candidate driver SNPs that exhibit linkage disequilibrium with an introgressed haplotype, highlighting the need to explore the possibility of adaptive introgression within these populations. This study reveals that the signatures of positive selection in highlanders and lowlanders of Papua New Guinea align closely with the challenges they face, which are specific to their environments.

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

A noncoding regulatory variant in IKZF1 increases acute lymphoblastic leukemia risk in Hispanic/Latino children

Hispanic/Latino children have the highest risk of acute lymphoblastic leukemia (ALL) in the US compared to other racial/ethnic groups, yet the basis of this remains incompletely understood. Through genetic fine-mapping analyses, we identified a new independent childhood ALL risk signal near IKZF1 in self-reported Hispanic/Latino individuals, but not in non-Hispanic White individuals, with an effect size of ∼1.44 (95% confidence interval = 1.33-1.55) and a risk allele frequency of ∼18% in Hispanic/Latino populations and <0.5% in European populations. This risk allele was positively associated with Indigenous American ancestry, showed evidence of selection in human history, and was associated with reduced IKZF1 expression. We identified a putative causal variant in a downstream enhancer that is most active in pro-B cells and interacts with the IKZF1 promoter. This variant disrupts IKZF1 autoregulation at this enhancer and results in reduced enhancer activity in B cell progenitors. Our study reveals a genetic basis for the increased ALL risk in Hispanic/Latino children.

Shape variation in modern human upper premolars

Morphological variation in modern human dentition is still an open field of study. The understanding of dental shape and metrics is relevant for the advancement of human biology and evolution and is thus of interest in the fields of dental anthropology, as well as human anatomy and medicine. Of concern is also the variation of the inner aspects of the crown which can be investigated using the tools and methods of virtual anthropology. In this study, we explored inter- and intra-population morphometric variation of modern humans' upper third and fourth premolars (P3s and P4s, respectively) considering both the inner and outer aspects of the crown, and discrete traits. We worked by means of geometric morphometrics on 3D image data from a geographically balanced sample of human populations from five continents, to analyse the shape of the dentinal crown, and the crown outline in 78 P3s and 76 P4s from 85 individuals. For the study of dental traits, we referred to the Arizona State University Dental Anthropology System integrated with more recent classification systems. The 3D shape variation of upper premolar crowns varied between short and mesio-distally broad, and tall and mesio-distally narrow. The observed shape variation was independent from the geographical origin of the populations, and resulted in extensive overlap. We noted a high pairwise correlation (r1 = 0.83) between upper P3s and P4s. We did not find any significant geographic differences in the analysed non-metric traits. Our outcomes thus suggest that geographical provenance does not play a determinant role in the shaping of the dental crown, whose genesis is under strict genetic control.