Inference of population structure using dense haplotype data

Phylogenomics and population genomics of Nothobranchius in lowland Tanzania: species delimitation and comparative genetic structure

Annual killifishes of the genus Nothobranchius are widespread across East Africa, with a particularly high biodiversity in lowland Tanzania. While they are typically found in ephemeral pools, the pools vary greatly in size, connectivity and inundation patterns. It was previously suggested that main river channels formed significant barriers to Nothobranchius dispersal. Here, we study the distribution of genetic lineages in an equatorial part of their range where main river channels that may act as barriers occur and closely related lineages frequently coexist in secondary contact zones. We used single-nucleotide polymorphism (SNP) dataset from double-digest restriction site-associated DNA (ddRAD) sequencing to investigate how genetic diversity is structured in Nothobranchius species from the coastal lowlands of Tanzania. Our analyses resolved some uncertain phylogenetic relationships within the N. melanospilus and N. guentheri species groups and placed N. flammicomantis outside the Coastal clade. Rather than a shared intraspecific genetic diversity pattern across four coexisting and widely distributed species, we found highly diverse patterns of intra-specific genetic structure among N. eggersi, N. janpapi, N. melanospilus and N. ocellatus. Populations of Nothobranchius species from the humid coastal lowlands of Tanzania are therefore structured, but not constrained by barriers formed by river channels or by basins - in contrast to Nothobranchius species from the dry part of their distribution. Some of the genetic relationships determined call for a re-evaluation of taxonomic delimitations.

Urbanization and genetic homogenization in the medieval Low Countries revealed through a ten-century paleogenomic study of the city of Sint-Truiden

BACKGROUND

Processes shaping the formation of the present-day population structure in highly urbanized Northern Europe are still poorly understood. Gaps remain in our understanding of when and how currently observable regional differences emerged and what impact city growth, migration, and disease pandemics during and after the Middle Ages had on these processes.

RESULTS

We perform low-coverage sequencing of the genomes of 338 individuals spanning the eighth to the eighteenth centuries in the city of Sint-Truiden in Flanders, in the northern part of Belgium. The early/high medieval Sint-Truiden population was more heterogeneous, having received migrants from Scotland or Ireland, and displayed less genetic relatedness than observed today between individuals in present-day Flanders. We find differences in gene variants associated with high vitamin D blood levels between individuals with Gaulish or Germanic ancestry. Although we find evidence of a Yersinia pestis infection in 5 of the 58 late medieval burials, we were unable to detect a major population-scale impact of the second plague pandemic on genetic diversity or on the elevated differentiation of immunity genes.

CONCLUSIONS

This study reveals that the genetic homogenization process in a medieval city population in the Low Countries was protracted for centuries. Over time, the Sint-Truiden population became more similar to the current population of the surrounding Limburg province, likely as a result of reduced long-distance migration after the high medieval period, and the continuous process of local admixture of Germanic and Gaulish ancestries which formed the genetic cline observable today in the Low Countries.

Population structure of Helicobacter pylori and antibiotic resistance-associated variants in a high-risk area of gastric cancer

The increasing antibiotic resistance of Helicobacter pylori has had a serious impact on gastric cancer prevention. Our study aimed to profile the genomic characteristics and explore variants associated with resistance in H. pylori strains from a high-risk area of gastric cancer in China. We isolated 153 strains from a community-based cohort and assessed their susceptibility to six antibiotics by MIC Test Strip and genomic characteristics by whole-genome sequencing. Phylogenetic analysis identified the strains as an independent cluster within H. pylori East Asian population (hpEastAsia). HefA, an efflux pump gene, showed the highest differentiation in the Linqu strains compared with the other Chinese strains. Bacterial genome-wide association study (GWAS) identified 86 resistance variants covering 44 genes. Novel resistance variants were found in lon and babA for metronidazole, HP1168 for clarithromycin, hcpC for levofloxacin, and sabA for rifamycin. Two newly identified hefA mutations (R229K and A283V) showed significant associations with metronidazole (P = 0.012) and tetracycline (P = 0.044) resistance, respectively. HefA mutations and GWAS variants were integrated with the significant literature-reported mutations to optimize the prediction models for metronidazole, levofloxacin, clarithromycin, and tetracycline resistance with area under the receiver operating characteristic curves of 0.82-0.93. Double-antibiotic resistance models were established for clinical applicability. Furthermore, hefA expression may play a potential mediating role in the associations between mutations and resistance. This study identified genetic independence in the representative H. pylori strains from a high-risk area of gastric cancer. Optimized resistance prediction panels, including novel hefA mutations and GWAS variants, may provide preliminary guidance for localized precise treatment and helpful experiences for the similar high-risk populations.IMPORTANCEHelicobacter pylori is a remarkable pathogen due to its virulence in gastric cancer and high genetic plasticity. Linqu County in China, a high-risk area of gastric cancer, faces serious antibiotic resistance issues and necessitates genomic profiling of local H. pylori strains. Phylogenetic analysis revealed the Linqu strains as a relatively independent cluster within the hpEastAsia population. Novel antibiotic resistance-associated hefA mutations and variants from our bacterial genome-wide association study in the Linqu strains were optimized to improve the prediction performances for single antibiotic and double-drug combination resistance compared with traditional literature-reported mutations. This study identified relative genetic independence and high differentiation in the representative H. pylori strains from a population with high risk of gastric cancer and high prevalence of antibiotic resistance. The optimized panels with novel variants improve antibiotic resistance prediction models compared with literature-reported mutations, providing guidance for localized precise treatment and suggesting prevention strategies for similar high-risk populations.

Molecular epidemiology of Ascaris lumbricoides following multiple rounds of community-wide treatment

Control and elimination of the parasite Ascaris lumbricoides relies on mass drug administration (MDA) using a limited number of anti-helminthics. Whilst these programs have reduced the infection intensity and prevalence within many endemic regions, patterns of transmission remain poorly understood. Reinfection commonly occurs following cessation of treatment due to the absence of acquired immunity post infection. Here, we utilise genomic data to understand parasite transmission within and between households in a community and the genomic impact of repeated MDA. We sequenced 54 whole-genomes from Ascaris worms obtained from individuals in a longitudinal cohort epidemiological study of transmission and drug treatment extending over 6 years. We found that fine-scale population structure exists in spatially distinct clusters of infected individuals with reinfection occurring within or between geographically close households. This observation helps inform the policy for future control in low prevalence settings suggesting more targeted treatment of infection hotspots. We found evidence of positive selection acting on members of gene families previously implicated in reduced drug efficacy but detected no impactful variants. As efforts to eliminate A. lumbricoides intensify, our study provides a foundation for genomic surveillance to help identify both who infects whom and the impact of repeated drug treatment.

Unique demographic history and population substructure among the Coorgs of Southern India

The fascinating genetic architecture of today's Indian population is the result of thousands of years of population mixing and eventual isolation. The Coorgs are one such small and religiously/socioculturally homogeneous community in Karnataka, India, whose origins and demographic history are much debated due to their stark sociocultural contrast with surrounding populations. Here, we analyzed Coorgs using both autosomal (n = 70) and uniparental markers (n = 144). Our analyses suggest population substructure among Coorgs and showed significant population drift in Coorg3 in both allele frequency and haplotype-based analysis methods. Further sharing of haplotype and identity by descent suggests a shared genetic history of Coorg1 with the Palliyar population, and founder event analysis clearly indicates that the founder event in Coorg1 was around 40 GBP (Generations Before Present). The demographic models based on fastGlobeTrotter and Moments highlighted the recent admixture of Coorg3 with the northwest Indian Sikh Jatt population (~23 GBP); and also showed that Coorg2 was formed by mixing Coorg1 and Coorg3 at ~11 GBP, explaining their current sociocultural homogeneity. F-statistics-based admixture graph models suggest an as yet unknown lineage in Coorg3. mtDNA analysis revealed about 40% South Asia-specific mitochondrial lineages in Coorgs; while Y chromosome analysis revealed a predominance of Eurasian, Middle Eastern, and Indian-specific haplogroups, suggesting male-mediated migration and eventual assimilation with native females. These insights into ancient and diverse genealogies among Coorgs not only explain their unique status in the Indian diaspora but also encourage further research to identify unknown migrations to the Indian subcontinent and thus further unravel its unique demography.

Factors Beyond Karstification Have Shaped the Population Structure of a Surface-Dwelling Minnow (Phoxinus lumaireul) Able to Disperse Underground

The Dinaric Karst, a biodiversity hotspot, features complex surface and subterranean hydrological networks that influence aquatic species distribution. This study investigates how karst hydrology shapes the genetic structure of the surface-dwelling minnow Phoxinus lumaireul, examining both large-scale and small-scale population patterns. Using mitochondrial DNA and genome-wide single nucleotide polymorphism (SNP) data of 827 specimens of P. lumaireul, three hypotheses were tested: (1) karst underground water connections facilitate genetic connectivity within and across river systems, whereas non-karst rivers exhibit genetic connectivity mostly within the same system; (2) historical and occasional hydrological connections have shaped present-day population structure, leaving genetic signatures of relatedness where no contemporary hydrological links exist; and (3) genomic approaches provide additional insights into biologically relevant connections that may not be captured by classical tracing tests. The large-scale analyses confirmed three main genetic groups (1a-c), whose structure was likely shaped by Pleistocene glaciations and associated microrefugia rather than by karst hydrology. Small-scale structure analyses revealed that while karst hydrology facilitated gene flow within specific areas, connectivity was uneven and influenced by local hydrological dynamics and historical admixture events. Furthermore, some underground pathways identified by classical tracing tests lacked evidence of genetic connectivity, underscoring the limitations of traditional methods and the added value of genomic data in indirectly detecting biologically relevant hydrological connections. These findings highlight the influence of both historical processes and contemporary karst hydrology on P. lumaireul populations, emphasizing their vulnerability in karst ecosystems and the need for targeted conservation efforts.

Exploring Mitonuclear Discordance: Ghost Introgression From an Ancient Extinction Lineage in the Odorrana swinhoana Complex

Mitonuclear discordance, the incongruence between mitochondrial DNA (mtDNA) and nuclear DNA (nuDNA), is a well-documented phenomenon with various potential explanations. One emerging hypothesis, ghost introgression, refers to the genetic contribution of an ancient, extinct or unsampled lineage and can now be tested using modern genomic data and demographic models. In this study, we investigated the evolutionary history of the Odorrana swinhoana complex (Anura: Ranidae), which includes O. swinhoana, O. utsunomiyaorum and an unidentified population with highly divergent mtDNA. While mitochondrial phylogeny suggested this population as a basal lineage, nuclear data from ddRADseq revealed it as a mixture of the most derived O. swinhoana nuclear sequences combined with ancient mtDNA. Demographic modelling further supported ghost introgression, as all models incorporating a ghost population outperformed those without it. These findings suggest that an eastward expansion of western O. swinhoana replaced an ancient Odorrana lineage, leaving only its mtDNA and fragments of its nuclear genome in the hybrid population. Our results provide one of the first documented cases of ghost introgression in amphibians and highlight its potential as a widespread evolutionary process. This study also underscores the risks of relying solely on mtDNA for phylogenetic reconstruction and species delimitation.

Leveraging Whole Genomes, Mitochondrial DNA and Haploblocks to Decipher Complex Demographic Histories: An Example From a Broadly Admixed Arctic Fish

The study of phylogeography has transitioned from mitochondrial haplotypes to genome-wide analyses, borrowing from population genomics methods along the way. Whole-genome sequencing allows the study of both mitochondrial and nuclear DNA and provides the density of markers to investigate recombination along the genome. This level of resolution could unravel complex histories of admixture between lineages, which are commonly observed in species evolving in recently deglaciated habitats. In this study, we sequenced 1120 Arctic Char genomes from 33 populations across Canada and Greenland to characterise patterns of genetic variation and diversity, and how they are shaped by hybridisation between the Arctic and Atlantic glacial lineages. Mitochondrial genomes across the study area were predominantly of Arctic origin, except in Greenland, where we observed some Atlantic descent. Through admixture analyses and demographic inferences on nuclear markers, we identified that all Canadian populations under the 66th parallel showed introgression from the Atlantic lineage, leading to higher genetic diversity. By scanning the genome using local principal component analyses, we identified putative large low-recombining haploblocks as local ancestry tracts from either lineage. Since haplotypes might retain different signatures of postglacial histories by sheltering sequences from recombination, we attempted to infer origins of recolonisation using whole genomes vs. ancestry tracts for the Arctic lineage. Despite limitations, we unveiled clues suggesting a complex postglacial history in Arctic Char. Overall, our study demonstrates that, even at low depth, making the most of whole-genome sequencing by analysing several genomic compartments provides a versatile and powerful way to address phylogeographic dynamics.

Genomic origin of Citrus reticulata "Unshiu"

Satsuma mandarin (Citrus reticulata "Unshiu") is a global cultivar with superior fruit characteristics and ranking among the top citrus cultivars in terms of production. It is also a key contributor to citrus breeding. However, the lack of high-quality genome makes the origin of Satsuma mandarin has long been a matter of debate. Here, we assembled a gap-free, high-quality genome of Satsuma mandarin. Meanwhile, we collected and sequenced 15 indigenous citrus varieties in Zhejiang Province, 12 Satsuma mandarins, 21 citrus hybrids related to Satsuma mandarin, 10 modern citrus varieties, and 7 other mandarins. Through high-resolution genome analysis, we inferred that Satsuma mandarin originated from a cross between C. reticulata "Ruju" × C. reticulata "Bendiguang" and proposed that Satsuma mandarin most probably originated in East area in Zhejiang Province of China, where the two parents-like cultivars are still found in a sympatric region to date. These results provide new insights into the origin model of Satsuma mandarin. The spread of mandarin is also discussed, which probably associated with the culture exchange and trade activities between Japan and China from Tang Dynasty and afterwards.

Submerged Corridors of Ancient Gene Flow in an Island Amphibian

Many island archipelagos sit on shallow continental shelves, and during the Pleistocene, these islands were often connected as global sea levels dropped following glaciation. Given a continental shelf only 30-60 m below sea level, the terrestrial biota of the Seychelles Archipelago likely dispersed amongst now isolated islands during the Pleistocene. Hypogeophis rostratus is an egg-laying, direct-developing caecilian amphibian found on 10 islands in the granitic Seychelles. Despite the seemingly limited dispersal abilities of this salt-intolerant amphibian, its distribution on multiple islands suggests likely historic dispersal across now submerged continental shelf corridors. We tested for the genetic signature of these historic corridors using fine-scale genomic data (ddRADseq). We found that genomic clusters often did not correspond to islands in the archipelago and that isolation-by-distance patterns were more consistent with gene flow across a continuous landscape than with isolated island populations. Using effective migration surfaces and ancestral range expansion prediction, we found support for contemporary populations originating near the large southern island of Mahé and dispersing to northern islands via the isolated Frégate island, with additional historic migration across the flat expanse of the Seychelles bank. Collectively, our results suggest that biogeographic patterns can retain signals from Pleistocene 'palaeo-islands' and that present-day islands can be thought of as hosting bottlenecks or transient refugia rather than discrete genetic units. Thus, the signatures of gene flow associated with palaeo-islands may be stronger than the isolating effects of contemporary islands in terrestrial species distributed on continental shelf islands.

Population genomics and molecular epidemiology of wheat powdery mildew in Europe

Agricultural diseases are a major threat to sustainable food production. Yet, for many pathogens we know exceptionally little about their epidemiological and population dynamics, and this knowledge gap is slowing the development of efficient control strategies. Here we study the population genomics and molecular epidemiology of wheat powdery mildew, a disease caused by the biotrophic fungus Blumeria graminis forma specialis tritici (Bgt). We sampled Bgt across two consecutive years, 2022 and 2023, and compiled a genomic dataset of 415 Bgt isolates from 22 countries in Europe and surrounding regions. We identified a single epidemic unit in the north of Europe, consisting of a highly homogeneous population. Conversely, the south of Europe hosts smaller local populations which are less interconnected. In addition, we show that the population structure can be largely predicted by the prevalent wind patterns. We identified several loci that were under selection in the recent past, including fungicide targets and avirulence genes. Some of these loci are common between populations, while others are not, suggesting different local selective pressures. We reconstructed the evolutionary history of one of these loci, AvrPm17, coding for an effector recognized by the wheat receptor Pm17. We found evidence for a soft sweep on standing genetic variation. Multiple AvrPm17 haplotypes, which can partially escape recognition by Pm17, spread rapidly throughout the continent upon its introduction in the early 2000s. We also identified a new virulent variant, which emerged more recently and can evade Pm17 resistance altogether. Overall, we highlight the potential of genomic surveillance in resolving the evolutionary and epidemiological dynamics of agricultural pathogens, as well as in guiding control strategies.

Standing genetic variation and introgression shape the cryptic radiation of Aquilegia in the mountains of Southwest China

Cryptic diversity in evolutionary radiation offers an excellent system for investigating the intricacies of evolutionary progress. Understanding the evolution of cryptic diversity is imperative for unraveling the hidden complexities of biodiversity. However, empirical evidence elucidating the mechanisms behind cryptic radiation remains limited, particularly in plants. Here, we focus on a monophyletic group of Aquilegia species mainly distributed in the mountains of Southwest China, one of the world's biodiversity hotspots. Using whole-genome resequencing of 158 individuals from 23 natural populations, we identify three to four paraphyletic lineages within each morphological species. Our findings reveal that 39 out of 43 detected instances of introgression occurred post-lineage formation. Identifying shared genomic regions indicates that the divergence of fixed singletons in lineages from morphological species A. kansuensis and A. rockii predates lineage formation, supporting a scenario where incomplete lineage sorting of standing variation contributes to morphological parallelism. Furthermore, strong positive correlations among genomic differentiation, divergence, and introgression suggest that standing variations and introgression from non-sister lineages contribute to the rapid genetic divergence. Our study illuminates the important roles of standing variations and introgression in plant cryptic radiation, advancing our understanding of the complex mechanisms behind the evolution of biodiversity in recent radiation events.

Inferring past demography and genetic adaptation in Spain using the GCAT cohort

Located in the southwestern corner of Europe, the Iberian Peninsula is separated from the rest of the continent by the Pyrenees Mountains and from Africa by the Strait of Gibraltar. This geographical position may have conditioned distinct selective pressures compared to the rest of Europe and influenced differential patterns of gene flow. In this work, we analyse 704 whole-genome sequences from the GCAT reference panel to quantify gene flow into Spain from various historical sources and identify the top signatures of positive (adaptive) selection. While we found no clear evidence of a 16th-century admixture event putatively related to the French diaspora during the Wars of Religion, we detected signals of North African admixture matching the Muslim period and the subsequent Christian Reconquista. Notably, besides finding that well-known candidate genes previously described in Eurasians also seem to be adaptive in Spain, we discovered novel top candidates for positive selection putatively associated with immunity and diet (UBL7, SMYD1, VAC14 and FDFT1). Finally, local ancestry deviation analysis revealed that the MHCIII genomic region underwent post-admixture selection following the post-Neolithic admixture with Steppe ancestry.

SURFBAT: a surrogate family based association test building on large imputation reference panels

Genotype-phenotype association tests are typically adjusted for population stratification using principal components that are estimated genome-wide. This lacks resolution when analyzing populations with fine structure and/or individuals with fine levels of admixture. This can affect power and precision, and is a particularly relevant consideration when control individuals are recruited using geographic selection criteria. Such is the case in France where we have recently created reference panels of individuals anchored to different geographic regions. To make correct comparisons against case groups, who would likely be gathered from large urban areas, new methods are needed. We present SURFBAT (a surrogate family based association test), which performs an approximation of the transmission-disequilibrium test. Our method hinges on the application of genotype imputation algorithms to match similar haplotypes between the case and control groups. This permits us to approximate local ancestry informed posterior probabilities of un-transmitted parental alleles of each case individual. This is achieved by assuming haplotypes from the imputation panel are well-matched for ancestry with the case individuals. When the first haplotype of an individual from the imputation panel matches that of a case individual, it is assumed that the second haplotype of the same reference individual can be used as a locally ancestry matched control haplotype and to approximately impute un-transmitted parental alleles. SURFBAT provides an association test that is inherently robust to fine-scale population stratification and opens up the possibility of efficiently using large imputation reference panels as control groups for association testing. In contrast to other methods for association testing that incorporate local-ancestry inference, SURFBAT does not require a set of ancestry groups to be defined, nor for local ancestry to be explicitly estimated. We demonstrate the interest of our tool on simulated datasets, as well as on a real-data example for a group of case individuals affected by Brugada syndrome.

Allopatric Speciation and Interspecific Gene Flow Driven by Niche Conservatism of Diploderma Tree Lizards in Taiwan

Allopatric speciation is a widely accepted hypothesis for species distributed across geographic barriers. Meanwhile, niche conservatism, the tendency of species to retain their ancestral ecological traits, helps reinforce genetic differentiation by stabilising species distributions over time and reducing the role of competition in shaping range boundaries. In contrast, hybridisation can occur at the edges of distribution after secondary contact following climatic or geological events, leading to a reduction in genetic divergence between divergent lineages. In this study, we investigated the role of geographic barriers, niche conservatism and gene flow in the speciation history of Diploderma species in Taiwan, where geographically distinct taxa share similar environmental preferences. By using ddRAD-seq data, seven distinct genetic clusters were identified with two putatively new cryptic species in D. brevipes and D. polygonatum. Most sister species pairs share similar climatic niches based on niche equivalency and similarity tests. We further detected significant historical gene flow between lineages of D. brevipes and D. polygonatum, where secondary contact might have occurred because of palaeoclimate changes and historical demographic expansion. Our results demonstrate that niche conservatism does not always act in concert to strengthen the result of allopatric speciation; instead, it may also lead to gene flow between divergent lineages following secondary contact. On the other hand, postdivergence gene flow may be a creating force generating phenotypic diversity in sexually selected traits in our study system. The underestimated species diversity of Diploderma in Taiwan requires further taxonomic work in the future.

The Genetic Variability of Present-Day Bulgarians Captures Ancient and Recent Ancestral Contributions

OBJECTIVES

Thanks to its pivotal crossroad position, Bulgaria played a fundamental key role during all the migration processes that interested the continent through time. While the genetic variability of the country has been deeply investigated using uniparental markers, previous genome-wide autosomal-based surveys mainly consisted of wider-range analyzes on Europe and the whole Balkan Peninsula. Here, we specifically focused on the Bulgarian population to recapitulate the main patterns of genomic variation and the major events that shaped the present-day genetic landscape.

METHODS

A total of 112 samples from seven highly representative areas of present-day Bulgaria were collected and genotyped for approximately 720 K genome-wide SNPs, and integrated with previously generated genomic data from wide modern and ancient reference panels to explore fine-scale relationship patterns and detail ancestral contributions.

RESULTS

In addition to the combination of ancient ancestries related to the early Mesolithic hunter-gatherers, Neolithic farmers, and Bronze Age Steppe pastoralists, both haplotype-based analyzes on modern populations and the comparisons with ancient genomes suggest the contribution of population processes that have occurred after the Roman rule and during the Medieval period in shaping the current Bulgarian genetic pool.

CONCLUSIONS

Our results align with previous evidence highlighting the impact that some historical events may have had not only in contributing to the ethnical and socio-cultural richness of present-day populations, but also in participating in the formation of the current genomic landscape. By providing new data from modern highly-representative samples, this study integrates further research to provide a comprehensive overview of the genetic history of Bulgaria.

The genetic footprint of the European Roma diaspora: evidence from the Balkans to the Iberian Peninsula

The Roma people have a complex demographic history shaped by their recent dispersal from a South Asian origin into Europe, accompanied by continuous population bottlenecks and gene flow. After settling in the Balkans around 1,000 years ago, the Roma gradually dispersed across Europe, and approximately 500 years ago, they established in the Iberian Peninsula what is now one of the largest Roma populations in Western Europe. Focusing specifically on the Iberian Roma, we conducted the most comprehensive genome-wide analysis of European Roma populations to date. Using allele frequency and haplotype-based methods, we analysed 181 individuals to investigate their genetic diversity, social dynamics, and migration histories at both continental and local scales. Our findings demonstrate significant gene flow from populations encountered during the Roma's dispersal and confirm their South Asian origins. We show that, between the 14th and 19th centuries, the Roma spread westward from the Balkans in various waves, with multiple admixture events. Furthermore, our findings refute previous hypotheses of a North African dispersal route into Iberia and genetic connections to Jewish populations. The Iberian Roma exhibit ten times greater genetic differentiation compared to non-Roma Iberians, indicating significant regional substructure. Additionally, we provide the first genetic evidence of assortative mating within Roma groups, highlighting distinct mating patterns and suggesting a gradual shift towards increased integration with non-Roma individuals. This study significantly enhances our understanding of how demographic history and complex genetic structure have shaped the genetic diversity of Roma populations, while also highlighting the influence of their evolving social dynamics.

A geographic history of human genetic ancestry

Describing the distribution of genetic variation across individuals is a fundamental goal of population genetics. We present a method that capitalizes on the rich genealogical information encoded in genomic tree sequences to infer the geographic locations of the shared ancestors of a sample of sequenced individuals. We used this method to infer the geographic history of genetic ancestry of a set of human genomes sampled from Europe, Asia, and Africa, accurately recovering major population movements on those continents. Our findings demonstrate the importance of defining the spatiotemporal context of genetic ancestry when describing human genetic variation and caution against the oversimplified interpretations of genetic data prevalent in contemporary discussions of race and ancestry.

Protocol for reconstructing ancestral genomes from present-day samples by applying local ancestry inference

The genome of the most recent common ancestor is generally not available but can greatly facilitate the inference of demographic history and the detection of local adaptations. Here, we present a protocol for applying local ancestry inference in present-day samples to reconstruct ancestral genomes. We describe steps for estimating haplotypes, inferring local ancestry, and assembling ancestral haplotypes. This protocol describes the analytic steps of reconstructing ancestral genomes using the example data of the Miao and She target populations. For complete details on the use and execution of this protocol, please refer to Gao et al.1.

Sparse haplotype-based fine-scale local ancestry inference at scale reveals recent selection on immune responses

Increasingly efficient methods for inferring the ancestral origin of genome regions are needed to gain insights into genetic function and history as biobanks grow in scale. Here we describe two near-linear time algorithms to learn ancestry harnessing the strengths of a Positional Burrows-Wheeler Transform. SparsePainter is a faster, sparse replacement of previous model-based 'chromosome painting' algorithms to identify recently shared haplotypes, whilst PBWTpaint uses further approximations to obtain lightning-fast estimation optimized for genome-wide relatedness estimation. The computational efficiency gains of these tools for fine-scale local ancestry inference offer the possibility to analyse large-scale genomic datasets using different approaches. Application to the UK Biobank shows that haplotypes better represent ancestries than principal components, whilst linkage-disequilibrium of ancestry identifies signals of recent changes to population-specific selection for many genomic regions associated with immune responses, suggesting avenues for understanding the pathogen-immune system interplay on a historical timescale.

Convergent evolution of complex adaptive traits modulates angiogenesis in high-altitude Andean and Himalayan human populations

Convergent adaptations represent paradigmatic examples of the capacity of natural selection to influence organisms' biology. However, the possibility to investigate the genetic determinants underpinning convergent complex adaptive traits has been offered only recently by methods for inferring polygenic adaptations from genomic data. Relying on this approach, we demonstrate how high-altitude Andean human groups experienced pervasive selective events at angiogenic pathways, which resemble those previously attested for Himalayan populations despite partial convergence at the single-gene level was observed. This provides additional evidence for the drivers of convergent evolution of enhanced blood perfusion in populations exposed to hypobaric hypoxia for thousands of years.

Global Population Structure, Virulence Factors and Antibiotic Resistance of Helicobacter pylori: A Pooled Analysis of 4067 Isolates From 76 Countries

BACKGROUND

Helicobacter pylori (H. pylori) is a common pathogen that has co-evolved with the human host for approximately 100,000 years; however, our understanding of its population structure remains limited. Furthermore, the detailed characteristics of its virulence factors and antibiotic resistance for H. pylori are not yet fully elucidated.

METHODS

In this study, we curated a global genome dataset of 4067 H. pylori isolates from 76 countries and explored H. pylori characteristics, including population genetic structure, virulence factors, and antibiotic resistance. We used three approaches (fineSTRUCTURE, ADMIXTURE, and DAPC) to infer the population structure of H. pylori. We investigated the virulence of each isolate by calling genotypes of cagA and vacA and evaluated the correlations of virulence factors with subpopulation. For antibiotic resistance, we identified mutations to determine the genotypic antibiotic resistance. Then we estimated the prevalence of genotypic antibiotic resistance grouped by geographical location, subpopulation, and study period.

RESULT

We identified 21 subpopulations in 4067 H. pylori isolates, including 20 previously reported subpopulations and a novel subpopulation hspEuropeIsrael, and found that the population structure of H. pylori was geographically restricted. The novel subpopulation hspEuropeIsrael had a higher proportion of less virulent cagA and vacA genotypes compared to other subpopulations. After evaluating the rates of H. pylori genotypic resistance to four antibiotics, we found that the prevalence of genotypic resistance to amoxicillin and metronidazole was > 15% across all five continents. Genotypic resistance to levofloxacin was > 15% on all continents except for Oceania. Additionally, the genotypic resistance rate to clarithromycin was > 15% in Asia, Europe, and Oceania. A trend of increased genotypic resistance over time was observed in several continents during subgroup analyses. Furthermore, we constructed a comprehensive database for H. pylori, named Helicobacter Pylori Encyclopedia for Research (HELPER, http://ccra.njmu.edu.cn/helper).

CONCLUSION

Our results provide a detailed characterization of H. pylori and extend previous schemas. HELPER serves as an informative and comprehensive database that will be a valuable resource for researchers and lay the foundation for future studies on H. pylori.

Origin and de novo domestication of sweet orange

Sweet orange is cultivated worldwide but suffers from various devastating diseases because of its monogenetic background. The elucidation of the origin of a crop facilitates the domestication of new crops that may better cope with new challenges. Here we collected and sequenced 226 citrus accessions and assembled telomere-to-telomere phased diploid genomes of sweet orange and sour orange. On the basis of a high-resolution haplotype-resolved genome analysis, we inferred that sweet orange originated from a sour orange × mandarin cross and confirmed this model using artificial hybridization experiments. We identified defense-related metabolites that potently inhibited the growth of multiple industrially important pathogenic bacteria. We introduced diversity to sweet orange, which showed wide segregation in fruit flavor and disease resistance and produced canker-resistant sweet orange by selecting defense-related metabolites. Our findings elucidate the origin of sweet orange and de novo domesticated disease-resistant sweet oranges, illuminating a strategy for the rapid domestication of perennial crops.

PANE: fast and reliable ancestral reconstruction on ancient genotype data with non-negative least square and principal component analysis

The history of human populations has been strongly shaped by admixture events, contributing to patterns of observed genetic diversity across populations. In this study, we introduce the Principal component Ancestry proportions using NNLS Estimation (PANE) method that leverages principal component analysis and non-negative least squares to assess the ancestral compositions of admixed individuals given a large set of populations. Our results show its ability to reliably estimate ancestry across several scenarios, even those with a significant proportion of missing genotypes, in a fraction of the time required when using other tools.

Fine-scale population structure and widespread conservation of genetic effect sizes between human groups across traits

Understanding genetic differences between populations is essential for avoiding confounding in genome-wide association studies and improving polygenic score (PGS) portability. We developed a statistical pipeline to infer fine-scale Ancestry Components and applied it to UK Biobank data. Ancestry Components identify population structure not captured by widely used principal components, improving stratification correction for geographically correlated traits. To estimate the similarity of genetic effect sizes between groups, we developed ANCHOR, which estimates changes in the predictive power of an existing PGS in distinct local ancestry segments. ANCHOR infers highly similar (estimated correlation 0.98 ± 0.07) effect sizes between UK Biobank participants of African and European ancestry for 47 of 53 quantitative phenotypes, suggesting that gene-environment and gene-gene interactions do not play major roles in poor cross-ancestry PGS transferability for these traits in the United Kingdom, and providing optimism that shared causal mutations operate similarly in different populations.

Epigenetic variation in light of population genetic practice

The evolutionary impact of epigenetic variation depends on its transgenerational stability and source - whether genetically determined, environmentally induced, or due to spontaneous, genotype-independent mutations. Here, we evaluate current approaches for investigating an independent role of epigenetics in evolution, pinpointing methodological challenges. We further identify opportunities arising from integrating epigenetic data with population genetic analyses in natural populations. Efforts to advance data quality, study design, and statistical treatment are encouraged to consolidate our understanding of the source of heritable epigenetic variation, quantify its autonomous potential for evolution, and enrich population genetic analyses with an additional layer of information.

Diversity in women and their vaginal microbiota

Women's health is essential to global societal and economic wellbeing, yet health disparities remain prevalent. The vaginal microbiota plays a critical role in health, with research indicating that reduced levels of core bacteria, such as lactobacilli, are associated with conditions like bacterial vaginosis (BV) and increased infection susceptibility. Lower levels of vaginal lactobacilli are reported more frequently in women of African and Latin American descent compared with women of European and Asian descent. However, geographical and other study inclusion and analysis biases influence current research. This opinion highlights the need for a more comprehensive understanding of a 'healthy' vaginal microbiome. It underscores efforts to broaden global research on microbiome diversity in socially relevant contexts, avoiding inappropriate applications of terms such as race and ethnicity.

Successful Traceability of Wildlife Samples Contributes to Wildlife Conservation: A Case Study of Tracing the Snub-Nosed Monkey (Rhinopithecus spp.)

Rapid and effective methods for tracing the geographic origin of wildlife samples are essential for tackling the illegal wildlife trade. Traditional morphological categorization methods are often inadequate as relying on the mitochondrial COXI barcode is insufficient for determining geographic populations. To address these limitations, we developed a bioinformatics-based pipeline for the rapid identification of traceable nuclear genome loci. This pipeline has been applied to the whole-genome sequence (WGS) data of China's flagship species, the snub-nosed monkey (Rhinopithecus spp.). These species are known for sex-biased dispersal and hybrid speciation, which complicates genealogy tracing. Using phylogenetic principles, we employed the Robinson and Foulds (RF) distance and scanned over 1,850,726 population-specific loci, identifying five pairs that can trace genealogy origins rapidly and cost-effectively using PCR. Additionally, we found that relying only on mitochondrial genetic information is insufficient for rapid and accurate traceability to subspecies-level geographic populations. Our pipeline efficiently identifies loci and traces the geographic origin of snub-nosed monkey individuals, providing a valuable tool for species preservation and combating the wildlife trade. This approach can be extended to other species, aiding in the conservation of endangered wildlife and tracing criminal evidence.

Population genomics and connectivity of Vazella pourtalesii sponge grounds of the northwest Atlantic with conservation implications of deep sea vulnerable marine ecosystems

Sponges are key ecosystem engineers that shape, structure and enhance the biodiversity of marine benthic communities globally. Sponge aggregations and reefs are recognized as vulnerable marine ecosystems (or VMEs) due to their susceptibility to damage from bottom-contact fishing gears. Ensuring their long-term sustainability, preservation, and ecosystem functions requires the implementation of sound scientific conservation tools. Here, the genetic diversity, structure, and connectivity of the deep-sea glass sponge, Vazella pourtalesii (Schmidt, 1870), was investigated using 1,102 neutral SNPs obtained in RADseq. This species is distributed across the northwest Atlantic from Florida, USA to Nova Scotia, Canada and we sequenced samples covering this full distribution and provided evidence of strong genetic structure with two distinct clusters: Florida together with the Carolina Shelves and the Scotian Shelf. We estimated moderate levels of diversity with low migration across large distances (> 1000 kms) and high connectivity at smaller scales (< 300 kms). Further, fishing pressure on genetic diversity was evaluated, within two Sponge Conservation Areas (SCAs) on the Scotian Shelf. Those areas have different disturbance histories, and cumulative fishing pressure. Slightly lower levels of genetic diversity were found inside the SCAs, and yet they encompassed a high proportion of the diversity observed within the Scotian Shelf. We provide baseline data for future monitoring of the SCAs, discussing our findings in the light of existing area-based management tools.

Demographic history and genetic variation of the Armenian population

We introduce a sizable (n = 34) whole-genome dataset on Armenians, a population inhabiting the region in West Asia known as the Armenian highlands. Equipped with this genetic data, we conducted a whole-genome study of Armenians and deciphered their fine-scale population structure and complex demographic history. We demonstrated that the Armenian populations from western, central, and eastern parts of the highlands are relatively homogeneous. The Sasun, a population in the south that had been argued to have received a major genetic contribution from Assyrians, was instead shown to have derived its slightly divergent genetic profile from a bottleneck that occurred in the recent past. We also investigated the debated question on the genetic origin of Armenians and failed to find any significant support for historical suggestions by Herodotus of their Balkan-related ancestry. We checked the degree of continuity of modern Armenians with ancient inhabitants of the eastern Armenian highlands and detected a genetic input into the region from a source linked to Neolithic Levantine Farmers at some point after the Early Bronze Age. Additionally, we cataloged an abundance of new mutations unique to the population, including a missense mutation predicted to cause familial Mediterranean fever, an autoinflammatory disorder highly prevalent in Armenians. Thus, we highlight the importance of further genetic and medical studies of this population.

Leveraging haplotype information in heritability estimation and polygenic prediction

Polygenic prediction has yet to make a major clinical breakthrough in precision medicine and psychiatry, where the application of polygenic risk scores is expected to improve clinical decision-making. Most widely used approaches for estimating polygenic risk scores are based on summary statistics from external large-scale genome-wide association studies, which rely on assumptions of matching data distributions. This may hinder the impact of polygenic risk scores in modern diverse populations due to small differences in genetic architectures. Reference-free estimators of polygenic scores are instead based on genomic best linear unbiased predictions and model the population of interest directly. We introduce a framework, named hapla, with a novel algorithm for clustering haplotypes in phased genotype data to estimate heritability and perform reference-free polygenic prediction in complex traits. We utilize inferred haplotype clusters to compute accurate heritability estimates and polygenic scores in a simulation study and the iPSYCH2012 case-cohort for depression disorders and schizophrenia. We demonstrate that our haplotype-based approach robustly outperforms standard genotype-based approaches, which can help pave the way for polygenic risk scores in the future of precision medicine and psychiatry.

Genomics Reveal Population Structure and Intergeneric Hybridization in an Endangered South American Bird: Implications for Management and Conservation

Genomics is an invaluable tool for conservation, particularly for endangered species impacted by wildlife trafficking. This study uses genomic data to provide new insights to aid conservation and management of endangered species, using as a case study the Yellow cardinal (Gubernatrix cristata), a bird endemic to southern South America severely affected by illegal trade and the transformation of its natural habitat. We explore population structure within the Yellow cardinal, delimiting management units and describing connectivity among them. Additionally, we develop and assess the accuracy of a panel of 189 informative SNPs, and demonstrate how these can reliably assign confiscated individuals to one of the management units established. Lastly, we assess hybridization between the Yellow cardinal and the Diuca finch (Diuca diuca), which is reported to occur in regions of sympatry. We confirm that hybridization occurs, although it is not as common as previously thought, and that hybrids might be fertile, as we found evidence of backcrossing with Yellow cardinals. We discuss the implications of this introgression for the evolution and conservation of Yellow cardinals. Our study provides new, valuable information that can guide conservation efforts, comprising a test case for the use of genomics in combating illegal trafficking, with potential application beyond the case of the Yellow cardinal.

High-resolution genomic history of early medieval Europe

Many known and unknown historical events have remained below detection thresholds of genetic studies because subtle ancestry changes are challenging to reconstruct. Methods based on shared haplotypes1,2 and rare variants3,4 improve power but are not explicitly temporal and have not been possible to adopt in unbiased ancestry models. Here we develop Twigstats, an approach of time-stratified ancestry analysis that can improve statistical power by an order of magnitude by focusing on coalescences in recent times, while remaining unbiased by population-specific drift. We apply this framework to 1,556 available ancient whole genomes from Europe in the historical period. We are able to model individual-level ancestry using preceding genomes to provide high resolution. During the first half of the first millennium CE, we observe at least two different streams of Scandinavian-related ancestry expanding across western, central and eastern Europe. By contrast, during the second half of the first millennium CE, ancestry patterns suggest the regional disappearance or substantial admixture of these ancestries. In Scandinavia, we document a major ancestry influx by approximately 800 CE, when a large proportion of Viking Age individuals carried ancestry from groups related to central Europe not seen in individuals from the early Iron Age. Our findings suggest that time-stratified ancestry analysis can provide a higher-resolution lens for genetic history.

Genomic Architecture Underlying the Striking Colour Variation in the Presence of Gene Flow for the Guinan Toad-Headed Lizard

How divergence occurs between closely related organisms in the absence of geographic barriers to gene flow stands as one of the long-standing questions in evolutionary biology. Previous studies suggested that the interplay between selection, gene flow and recombination strongly affected the process of divergence with gene flow. However, the extent to which these forces interact to drive divergence remains largely ambiguous. Guinan toad-headed lizards (Phrynocephalus guinanensis) in the Mugetan Desert exhibit striking colour differences from lizards outside the desert and provide an excellent model to address this question. Through extensive sampling and whole genome sequencing, we obtained genotypes for 191 samples from 14 populations inside and outside the desert. Despite the colour differences, continuous and asymmetric gene flow was detected across the desert border. More importantly, 273 highly diverged regions (HDRs) were identified between them, accounting only for 0.47% of the genome but widely distributed across 20 (out of the total 24) chromosomes. Strong signatures of selection were identified in HDRs, and local recombination rates were repressed. Furthermore, five HDRs exhibited significantly higher divergence, which contained key genes associated with crucial functions in animal coloration, including pteridine and melanocyte pigmentation. Genes related to retinal cells and steroid hormones were identified in other HDRs, which might have also contributed to the formation of colour variation in the presence of gene flow. This study provided novel insights into the understanding of the evolutionary mechanisms of genetic divergence in the presence of gene flow.

Continental influx and pervasive matrilocality in Iron Age Britain

Roman writers found the relative empowerment of Celtic women remarkable1. In southern Britain, the Late Iron Age Durotriges tribe often buried women with substantial grave goods2. Here we analyse 57 ancient genomes from Durotrigian burial sites and find an extended kin group centred around a single maternal lineage, with unrelated (presumably inward migrating) burials being predominantly male. Such a matrilocal pattern is undescribed in European prehistory, but when we compare mitochondrial haplotype variation among European archaeological sites spanning six millennia, British Iron Age cemeteries stand out as having marked reductions in diversity driven by the presence of dominant matrilines. Patterns of haplotype sharing reveal that British Iron Age populations form fine-grained geographical clusters with southern links extending across the channel to the continent. Indeed, whereas most of Britain shows majority genomic continuity from the Early Bronze Age to the Iron Age, this is markedly reduced in a southern coastal core region with persistent cross-channel cultural exchange3. This southern core has evidence of population influx in the Middle Bronze Age but also during the Iron Age. This is asynchronous with the rest of the island and points towards a staged, geographically granular absorption of continental influence, possibly including the acquisition of Celtic languages.

Human migration from the Levant and Arabia into Yemen since Last Glacial Maximum

While a broad consensus about the first successful migration modern humans out of Africa seems established, the peopling of Arabia remains somewhat enigmatic. Identifying the ancestral populations that contributed to the gene pool of the current populations inhabiting Arabia and the impact of their contributions remains a challenging task. We investigate the genetic makeup of the current Yemeni population using 46 whole genomes and 169 genotype arrays derived from Yemeni individuals from all geographic regions across Yemen and 351 genotype arrays derived from neighboring populations providing regional context. Principal Component Analysis shows stratification between Yemen districts but also with respect to nearby populations: Yemeni, other Arabian and Bedouin samples form a continuum towards the populations of the Levant, whereas East Africa and India appear strongly differentiated. This finding is further supported by higher Principal Components, admixture and haplogroup analyses, and F-statistics. Moreover, two-reference linkage disequilibrium decay estimates are most significant for Yemeni admixture from an ancient northern influx (up to 5220BP from Palestine) and East Africa (750BP). We show that the initial gene flow into the Yemeni populations of today came from the rest of Arabia and the Levant, and a less substantial and more recent genetic impact into coastal Yemen from East Africa, particularly.

Reconstructing the population history of the Nicobarese

The Nicobarese are the major tribal groups in the Nicobar district, situated south of the Andaman group of Islands. Linguistic phylogeny suggests that the linguistic ancestors of the Nicobarese settled the Nicobar archipelago in the early Holocene. So far, genetic research on them is low-resolution and restricted to the haploid DNA markers. Therefore, in the present analysis, we have used the high-resolution biparental (1554 published and 5 newly genotyped Nicobarese individuals) and uniparental genetic markers and looked at the genetic association of Nicobarese with the South and Southeast Asian populations. We report a common ancestral component shared among the Austroasiatic of South and Southeast Asia. Our analyses have suggested that the Nicobarese peoples retain this ancestral Austroasiatic predominant component in their genomes in the highest proportion. On the Southeast Asian mainland, the Htin Mal, who speak an Austroasiatic language of the Khmuic branch, represent a population that has preserved their ethnic distinctness from other groups over time and consequently shown the highest drift with the Nicobarese. The analysis based on haplotypes indicated a significant level of genomic segment sharing across linguistic groups, indicating an ancient broader distribution of Austroasiatic populations in Southeast Asia. Based on the temporal analyses of haploid DNA, it is suggested that the forebears of the Nicobarese people may have arrived on the Nicobar Islands in the last 5000 YBP. Therefore, among the modern populations, the Nicobarese peoples and the Htin Mal language community represent good genetic proxies for ancient Austroasiatics.

Ancient genomics support deep divergence between Eastern and Western Mediterranean Indo-European languages

The Indo-European languages are among the most widely spoken in the world, yet their early diversification remains contentious1-5. It is widely accepted that the spread of this language family across Europe from the 5th millennium BP correlates with the expansion and diversification of steppe-related genetic ancestry from the onset of the Bronze Age6,7. However, multiple steppe-derived populations co-existed in Europe during this period, and it remains unclear how these populations diverged and which provided the demographic channels for the ancestral forms of the Italic, Celtic, Greek, and Armenian languages8,9. To investigate the ancestral histories of Indo-European-speaking groups in Southern Europe, we sequenced genomes from 314 ancient individuals from the Mediterranean and surrounding regions, spanning from 5,200 BP to 2,100 BP, and co-analysed these with published genome data. We additionally conducted strontium isotope analyses on 224 of these individuals. We find a deep east-west divide of steppe ancestry in Southern Europe during the Bronze Age. Specifically, we show that the arrival of steppe ancestry in Spain, France, and Italy was mediated by Bell Beaker (BB) populations of Western Europe, likely contributing to the emergence of the Italic and Celtic languages. In contrast, Armenian and Greek populations acquired steppe ancestry directly from Yamnaya groups of Eastern Europe. These results are consistent with the linguistic Italo-Celtic10,11 and Graeco-Armenian1,12,13 hypotheses accounting for the origins of most Mediterranean Indo-European languages of Classical Antiquity. Our findings thus align with specific linguistic divergence models for the Indo-European language family while contradicting others. This underlines the power of ancient DNA in uncovering prehistoric diversifications of human populations and language communities.

Estimation of genetic admixture proportions via haplotypes

Estimation of ancestral admixture is essential for creating personal genealogies, studying human history, and conducting genome-wide association studies (GWAS). The following three primary methods exist for estimating admixture coefficients. The frequentist approach directly maximizes the binomial loglikelihood. The Bayesian approach adds a reasonable prior and samples the posterior distribution. Finally, the nonparametric approach decomposes the genotype matrix algebraically. Each approach scales successfully to datasets with a million individuals and a million single nucleotide polymorphisms (SNPs). Despite their variety, all current approaches assume independence between SNPs. To achieve independence requires performing LD (linkage disequilibrium) filtering before analysis. Unfortunately, this tactic loses valuable information and usually retains many SNPs still in LD. The present paper explores the option of explicitly incorporating haplotypes in ancestry estimation. Our program, HaploADMIXTURE, operates on adjacent SNP pairs and jointly estimates their haplotype frequencies along with admixture coefficients. This more complex strategy takes advantage of the rich information available in haplotypes and ultimately yields better admixture estimates and better clustering of real populations in curated datasets.

Adjusting for principal components can induce collider bias in genome-wide association studies

Principal component analysis (PCA) is widely used to control for population structure in genome-wide association studies (GWAS). Top principal components (PCs) typically reflect population structure, but challenges arise in deciding how many PCs are needed and ensuring that PCs do not capture other artifacts such as regions with atypical linkage disequilibrium (LD). In response to the latter, many groups suggest performing LD pruning or excluding known high LD regions prior to PCA. However, these suggestions are not universally implemented and the implications for GWAS are not fully understood, especially in the context of admixed populations. In this paper, we investigate the impact of pre-processing and the number of PCs included in GWAS models in African American samples from the Women's Health Initiative SNP Health Association Resource and two Trans-Omics for Precision Medicine Whole Genome Sequencing Project contributing studies (Jackson Heart Study and Genetic Epidemiology of Chronic Obstructive Pulmonary Disease Study). In all three samples, we find the first PC is highly correlated with genome-wide ancestry whereas later PCs often capture local genomic features. The pattern of which, and how many, genetic variants are highly correlated with individual PCs differs from what has been observed in prior studies focused on European populations and leads to distinct downstream consequences: adjusting for such PCs yields biased effect size estimates and elevated rates of spurious associations due to the phenomenon of collider bias. Excluding high LD regions identified in previous studies does not resolve these issues. LD pruning proves more effective, but the optimal choice of thresholds varies across datasets. Altogether, our work highlights unique issues that arise when using PCA to control for ancestral heterogeneity in admixed populations and demonstrates the importance of careful pre-processing and diagnostics to ensure that PCs capturing multiple local genomic features are not included in GWAS models.

Computational Genomics and Its Applications to Anthropological Questions

The advent of affordable genome sequencing and the development of new computational tools have established a new era of genomic knowledge. Sequenced human genomes number in the tens of thousands, including thousands of ancient human genomes. The abundance of data has been met with new analysis tools that can be used to understand populations' demographic and evolutionary histories. Thus, a variety of computational methods now exist that can be leveraged to answer anthropological questions. This includes novel likelihood and Bayesian methods, machine learning techniques, and a vast array of population simulators. These computational tools provide powerful insights gained from genomic datasets, although they are generally inaccessible to those with less computational experience. Here, we outline the theoretical workings behind computational genomics methods, limitations and other considerations when applying these computational methods, and examples of how computational methods have already been applied to anthropological questions. We hope this review will empower other anthropologists to utilize these powerful tools in their own research.

Evolutionary genomics of Leishmania braziliensis across the neotropical realm

The Neotropical realm, one of the most biodiverse regions on Earth, houses a broad range of zoonoses that pose serious public health threats. Protozoan parasites of the Leishmania (Viannia) braziliensis clade cause zoonotic leishmaniasis in Latin America with clinical symptoms ranging from simple cutaneous to destructive, disfiguring mucosal lesions. We present the first comprehensive genome-wide continental study including 257 cultivated isolates representing most of the geographical distribution of this major human pathogen. The L. braziliensis clade is genetically highly heterogeneous, consisting of divergent parasite groups that are associated with different environments and vary greatly in diversity. Apart from several small ecologically isolated groups with little diversity, our sampling identifies two major parasite groups, one associated with the Amazon and the other with the Atlantic Forest biomes. These groups show different recombination histories, as suggested by high levels of heterozygosity and effective population sizes in the Amazonian group in contrast to high levels of linkage and clonality in the Atlantic group. We argue that these differences are linked to strong eco-epidemiological differences between the two regions. In contrast to geographically focused studies, our study provides a broad understanding of the molecular epidemiology of zoonotic parasites circulating in tropical America.

AncestryGrapher toolkit: Python command-line pipelines to visualize global- and local- ancestry inferences from the RFMIX version 2 software

SUMMARY

Admixture is a fundamental process that has shaped levels and patterns of genetic variation in human populations. RFMIX version 2 (RFMIX2) utilizes a robust modeling approach to identify the genetic ancestries in admixed populations. However, this software does not have a built-in method to visually summarize the results of analyses. Here, we introduce the AncestryGrapher toolkit, which converts the numerical output of RFMIX2 into graphical representations of global and local ancestry (i.e. the per-individual ancestry components and the genetic ancestry along chromosomes, respectively).

RESULTS

To demonstrate the utility of our methods, we applied the AncestryGrapher toolkit to visualize the global and local ancestry of individuals in the North African Mozabite Berber population from the Human Genome Diversity Panel. Our results showed that the Mozabite Berbers derived their ancestry from the Middle East, Europe, and sub-Saharan Africa (global ancestry). We also found that the population origin of ancestry varied considerably along chromosomes (local ancestry). For example, we observed variance in local ancestry in the genomic region on Chromosome 2 containing the regulatory sequence in the MCM6 gene associated with lactase persistence, a human trait tied to the cultural development of adult milk consumption. Overall, the AncestryGrapher toolkit facilitates the exploration, interpretation, and reporting of ancestry patterns in human populations.

AVAILABILITY AND IMPLEMENTATION

The AncestryGrapher toolkit is free and open source on https://github.com/alisi1989/RFmix2-Pipeline-to-plot.

An ancient ecospecies of Helicobacter pylori

Helicobacter pylori disturbs the stomach lining during long-term colonization of its human host, with sequelae including ulcers and gastric cancer1,2. Numerous H. pylori virulence factors have been identified, showing extensive geographic variation1. Here we identify a 'Hardy' ecospecies of H. pylori that shares the ancestry of 'Ubiquitous' H. pylori from the same region in most of the genome but has nearly fixed single-nucleotide polymorphism differences in 100 genes, many of which encode outer membrane proteins and host interaction factors. Most Hardy strains have a second urease, which uses iron as a cofactor rather than nickel3, and two additional copies of the vacuolating cytotoxin VacA. Hardy strains currently have a limited distribution, including in Indigenous populations in Siberia and the Americas and in lineages that have jumped from humans to other mammals. Analysis of polymorphism data implies that Hardy and Ubiquitous coexisted in the stomachs of modern humans since before we left Africa and that both were dispersed around the world by our migrations. Our results also show that highly distinct adaptive strategies can arise and be maintained stably within bacterial populations, even in the presence of continuous genetic exchange between strains.

Population Genomic Scans for Natural Selection and Demography

Uncovering the fundamental processes that shape genomic variation in natural populations is a primary objective of population genetics. These processes include demographic effects such as past changes in effective population size or gene flow between structured populations. Furthermore, genomic variation is affected by selection on nonneutral genetic variants, for example, through the adaptation of beneficial alleles or balancing selection that maintains genetic variation. In this article, we discuss the characterization of these processes using population genetic models, and we review methods developed on the basis of these models to unravel the underlying processes from modern population genomic data sets. We briefly discuss the conditions in which these approaches can be used to infer demography or identify specific nonneutral genetic variants and cases in which caution is warranted. Moreover, we summarize the challenges of jointly inferring demography and selective processes that affect neutral variation genome-wide.

Nested admixture during and after the Trans-Atlantic Slave Trade on the island of São Tomé

Human admixture history is rarely a simple process in which distinct populations, previously isolated for a long time, come into contact once to form an admixed population. In this study, we aim to reconstruct the complex admixture histories of the population of São Tomé, an island in the Gulf of Guinea that was the site of the first slave-based plantation economy, and experienced successive waves of forced and deliberate migration from Africa. We examined 2.5 million SNPs newly genotyped in 96 São Toméans and found that geography alone cannot explain the observed patterns of genetic differentiation within the island. We defined five genetic groups in São Tomé based on the hypothesis that individuals sharing the most haplotypes are more likely to share similar genetic histories. Using Identical-by-Descent and different local ancestry inference methods, we inferred shared ancestries between 70 African and European populations and each São Toméan genetic group. We identified admixture events between admixed groups that were previously isolated on the island, showing how recently admixed populations can be themselves the sources of other admixture events. This study demonstrates how complex admixture and isolation histories during and after the Transatlantic Slave-Trade shaped extant individual genetic patterns at a local scale in Africa.

Population connectivity and size reductions in the Anthropocene: the consequence of landscapes and historical bottlenecks in white forsythia fragmented habitats

BACKGROUND

White forsythia (Abeliophyllum distichum) is an endangered Korean Peninsula endemic that has been subjected to recent population genomics studies using SNPs via RAD sequencing. Here, we primarily employed the often underutilized haplotype information from RAD loci to further describe the species' previously uninvestigated haplotype-based genomic variation and structure, and genetic-geographic characteristics and gene flow patterns among its five earlier identified genetic groups. We also inferred the time of past events that may have impacted the effective population size of these groups, as well as the species' potential future distribution amidst the warming climate and anthropogenic threats.

RESULTS

Our findings emphasized the recognition of the species' regional patterns of genetic structure, and the role of topography and its associated gene flow patterns as some of the possible factors that may have influenced the species' present-day fragmented population distribution. The inferred bottleneck events during the Anthropocene, some of which aligned with the time of historical catastrophic events on the Peninsula (e.g., the Korean War), were revealed to have contributed to the generally low effective population size of its five lineages, particularly those with marginal distributional range. Future distribution under both optimistic and pessimistic climatic scenarios suggests unlikely suitable habitats for these populations to expand from their current range limits, at least in the next 80 years.

CONCLUSIONS

The small effective population size and landscape-driven limited gene flow among white forsythia populations will remain a big challenge for the conservation management of the species' already fragmented population distribution. To help mitigate these impacts, the merging of various research approaches and the use of genomic data to their full potential is recommended to provide the optimized knowledge-based tools for the conservation of this endangered species, and other similar plants under pressure.

Genetic diversity of the largest African mole-rat genus, Bathyergus. One, two or four species?

Recent advances in sequencing technology and phylogenetic methods allow us to solve puzzling taxonomic questions using detailed analyses of genetic diversity of populations and gene flow between them. The genus of solitary-living dune mole-rat, Bathyergus, is quite unique among six genera of African mole-rats. The animals are by far the largest and the only scratch digging mole-rat genus possessing a skull less adapted to digging, grooved upper incisors, and more surface locomotor activity. Most authors recognize two species of dune mole-rats, B. suillus and B. janetta, but according to others, the genus is monotypic. In addition, recent molecular studies have revealed cryptic genetic diversity and suggested the existence of up to four species. In our study, we used mitochondrial and genome-wide nuclear data collected throughout the distribution of the genus to investigate the number of species. In agreement with previous studies, we found Bathyergus to be differentiated into several distinct lineages, but we also found evidence for a degree of gene flow between some of them. Furthermore, we confirmed that B. janetta is nested within B. suillus, making the latter paraphyletic and we documented an instance of local mitochondrial introgression between these two nominal species. Phylogeographic structure of the genus was found to be very shallow. Although traditionally dated to the Miocene, we found the first split within the genus to be much younger estimated to 0.82 Ma before present. Genealogical distinctiveness of some lineages was very low, and the coancestry matrix showed extensive sharing of closely related haplotypes throughout the genus. Accordingly, Infomap clustering on the matrix showed all populations to form a single cluster. Overall, our study tends to support the existence of only one species of Bathyergus namely, B. suillus. Environmental niche modelling confirmed its dependence on sandy soils and the preference for soils with relatively high carbon content. Bayesian skyline plots indicate recent population decline in the janetta lineage, probably related to global environmental change.

Signals of selection and ancestry in independently feral Gallus gallus populations

Recent work indicates that feralisation is not a simple reversal of domestication, and therefore raises questions about the predictability of evolution across replicated feral populations. In the present study we compare genes and traits of two independently established feral populations of chickens (Gallus gallus) that inhabit archipelagos within the Pacific and Atlantic regions to test for evolutionary parallelism and/or divergence. We find that feral populations from each region are genetically closer to one another than other domestic breeds, despite their geographical isolation and divergent colonisation histories. Next, we used genome scans to identify genomic regions selected during feralisation (selective sweeps) in two independently feral populations from Bermuda and Hawaii. Three selective sweep regions (each identified by multiple detection methods) were shared between feral populations, and this overlap is inconsistent with a null model in which selection targets are randomly distributed throughout the genome. In the case of the Bermudian population, many of the genes present within the selective sweeps were either not annotated or of unknown function. Of the nine genes that were identifiable, five were related to behaviour, with the remaining genes involved in bone metabolism, eye development and the immune system. Our findings suggest that a subset of feralisation loci (i.e. genomic targets of recent selection in feral populations) are shared across independently established populations, raising the possibility that feralisation involves some degree of parallelism or convergence and the potential for a shared feralisation 'syndrome'.

Graphite: painting genomes using a colored de Bruijn graph

The recent growth of microbial sequence data allows comparisons at unprecedented scales, enabling the tracking of strains, mobile genetic elements, or genes. Querying a genome against a large reference database can easily yield thousands of matches that are tedious to interpret and pose computational challenges. We developed Graphite that uses a colored de Bruijn graph (cDBG) to paint query genomes, selecting the local best matches along the full query length. By focusing on the best genomic match of each query region, Graphite reduces the number of matches while providing the most promising leads for sequence tracking or genomic forensics. When applied to hundreds of Campylobacter genomes we found extensive gene sharing, including a previously undetected C. coli plasmid that matched a C. jejuni chromosome. Together, genome painting using cDBGs as enabled by Graphite, can reveal new biological phenomena by mitigating computational hurdles.