Punctuated bursts in human male demography inferred from 1,244 worldwide Y-chromosome sequences

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.

Exploring Y-chromosomal STRs and SNPs for forensic and genetic insights in the Jiangsu Han population

Y-chromosome short tandem repeats (Y-STRs) and single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used for individual identification, forensic applications, and the study of paternal lineage history. This study analyzed the genetic polymorphism and paternal genetic structure of the Han population in Jiangsu Province by examining 374 unrelated male individuals using 29 Y-STRs and 183 Y-SNPs. Forensic parameters were calculated, and the discriminatory power of five Y-STR systems (MHT, EXT, PPY12, Yfiler, and Y29) was compared. Genetic structure was assessed in the context of the Jiangsu Han and other Chinese populations. Results showed that the Y29 system had the highest discriminatory capacity, identifying 374 unique haplotypes with HD and DC values of 1. Seven major haplogroups (C, D, J, K, O, Q, R) and 83 terminal haplogroups were identified, with haplogroup O being the most predominant (approximately 85%). Subdivision of haplogroup O revealed that the Jiangsu Han population exhibits genetic characteristics of both Southern and Northern Han groups. Population genetic analyses further confirmed that the Jiangsu Han clustered closely with Southern Han populations. Genetic admixture results revealed that the Jiangsu Han population derives 89% of their ancestry from Southern Han populations, while retaining 11% Northern Han contributions. The study provides valuable insights into the genetic structure of the Jiangsu Han population, with significant implications for forensic genetics, anthropological research, and broader population genetics studies.

A generalized structured coalescent for purifying selection without recombination

Purifying selection is a critical factor in shaping genetic diversity. Current theoretical models mostly address scenarios of either very weak or strong selection, leaving a significant gap in our knowledge. The effects of purifying selection on patterns of genomic diversity remain poorly understood when selection against deleterious mutations is weak to moderate, particularly when recombination is limited or absent. In this study, we extend an existing approach, the fitness-class coalescent, to incorporate arbitrary levels of purifying selection in haploid populations. This model offers a comprehensive framework for exploring the influence of purifying selection in a wide range of demographic scenarios. Moreover, our research reveals potential sources of qualitative and quantitative biases in demographic inference, highlighting the significant risk of attributing genetic patterns to past demographic events rather than purifying selection. This work expands our understanding of the complex interplay between selection, drift, and population dynamics, and how purifying selection distorts demographic inference.

Multiple origins and phenotypic implications of an extended human pseudoautosomal region shown by analysis of the UK Biobank

The 2.7-Mb major pseudoautosomal region (PAR1) on the short arms of the human X and Y chromosomes plays a critical role in meiotic sex chromosome segregation and male fertility and has been regarded as evolutionarily stable. However, some European Y chromosomes belonging to Y haplogroups (Y-Hgs) R1b and I2a carry an ∼115-kb extension (ePAR [extended PAR]) arising from X-Y non-allelic homologous recombination (NAHR). To investigate the diversity, history, and dynamics of ePAR formation, we screened for its presence, and that of the predicted reciprocal X chromosome deletion, among ∼218,300 46,XY males of the UK Biobank (UKB), a cohort associated with longitudinal clinical data. The UKB incidence of ePAR is ∼0.77%, and that of the deletion is ∼0.02%. We found that Y-Hg I2a sub-lineages accounted for nearly 90% of ePAR cases but, by Y haplotyping and breakpoint sequencing, determined that, in total, there have been at least 18 independent ePAR origins, associated with nine different Y-Hgs. We found examples of ePAR linked to Y-Hg K among men of self-declared Pakistani ancestry and Y-Hg E1, typical of men with African ancestry, showing that ePAR is not restricted to Europeans. ePAR formation is likely random, with high frequencies in some Y-Hgs arising through drift and male-mediated expansions. Sequencing recombination junction fragments identified likely reciprocal events, and the heterogeneity of ePAR and X-deletion junctions highlighted the recurrent nature of the NAHR events. A phenome-wide association study revealed an association between ePAR and elevated levels of circulating IGF-1 as well as musculoskeletal phenotypes.

Paleolithic divergence and multiple Neolithic expansions of ancestral nomadic emperor-related paternal lineages

The reconstruction of demographic history using ancient and modern genomic resources reveals extensive interactions and admixture between ancient nomadic pastoralists and the social organizations of the Chinese Central Plain. However, the extent to which Y-chromosome genetic legacies from nomadic emperor-related ancestral lineages influence the Chinese paternal gene pool remains unclear. Here, we genotype 2717 ethnolinguistically diverse samples belonging to C2a lineages, perform whole-genome sequencing on 997 representative samples, and integrate these data with ancient genomic sequences. We reconstruct the evolutionary histories of Northern Zhou-, Qing emperor-, and pastoralist-related lineages to assess their genetic impact on modern Chinese populations. This reassembled fine-scale Y-chromosome phylogeny identifies deep divergence and five Neolithic expansion events contributing differently to the formation of northern Chinese populations. Phylogeographic modeling indicates that the nomadic empires of the Northern Zhou and Qing dynasties genetically originated from the Mongolian Plateau. Phylogenetic topology and shared haplotype patterns show that three upstream ancestors of Northern Zhou (C2a1a1b1a2a1b-FGC28857), Donghu tribe (C2a1a1b1-F1756), and Qing (C2a1a3a2-F10283) emperor-related lineages expanded during the middle Neolithic, contributing significantly to genetic flow between ancient northeastern Asians and modern East Asians. Notably, this study reveals limited direct contributions of Emperor Wu of Northern Zhou's lineages to modern East Asians.

Sex-Biased Admixture Followed by Isolation and Adaptive Evolution Shaped the Genomic and Blood Pressure Diversity of the LopNur People

The LopNur people are an ethnic group living on the edge of the Taklamakan Desert, and they are believed to demonstrate a unique genetic makeup due to their isolation and limited contact with neighboring populations. However, a lack of genetic studies on the LopNur people has resulted in limited knowledge about their ancestral origins and demographic history. Here, we conducted the first whole-genome sequencing study of 164 LopNur individuals (LOP) to gain insight into their genetic history and adaptive evolution in an isolated desert area. Our analysis revealed that the present-day LOP have experienced a complex history of admixture followed by long-term isolation, with their ancestry derived from East Asia (∼41.46%), West Eurasia (∼26.43%), Siberia (∼24.27%), and South Asia (∼7.82%). Notably, a remarkable sex-biased admixture occurred between Western males and Eastern females. In addition to complex admixture followed by long-term geographic isolation and further recent migrations, adaptive evolution jointly formed the gene pool and phenotypic diversity of the present-day LOP. Intriguingly, our analysis suggests that the USP35-GAB2 region may be correlated with blood pressure in LOP, based on a joint analysis of genomics and blood pressure data. Moreover, we identified two variants, rs7387065, and rs2229437, located on CSMD1 and PRCP, respectively. These variants exhibited frequency differences between Asian and European populations and were reported to be associated with antihypertensive drug absorption. Our results provide new insight into the complex history of the LOP, an admixed and isolated ethnic group residing at the crossroads of East and West, a case with ancient admixture, long-term isolation, adaptive evolution, and sex-biased gene flow.

Evaluating genome-wide and targeted forensic sequencing approaches to kinship determination

Kinship determination is a valuable tool in forensic genetics, with applications including familial searching, disaster victim identification, and investigative genetic genealogy. Conventional typing of small numbers of autosomal short tandem repeats (STRs) confidently identifies only first-degree relatives. Massively parallel sequencing (MPS) can access more STRs and resolve alleles identical by length but differing in sequence (isoalleles), which may increase the power of kinship estimation, particularly when combined with additional sequenced single nucleotide polymorphism (SNP) loci, as in the ForenSeq DNA Signature Prep kit. MPS sequencing of ∼10,000 SNPs is available in the ForenSeq Kintelligence kit, promising detection of more distant kin, while SNP chips carrying hundreds of thousands of markers increase resolution still further. Here we evaluate these different resolutions in a set of pedigrees, and via simulations. As expected, the key factor influencing the precision of kinship estimation is the number of markers analysed and MPS-based analysis of STRs increases resolution, with the full set of ForenSeq DNA Signature Prep kit markers allowing detection of third-degree relatives. Since SNP chips include non-autosomal (X- and Y-chromosomal, and mitochondrial [mtDNA]) markers, we ask how these perform within the pedigrees, cross-referencing to Y-STR sequence data. We highlight the importance of understanding haplogroup resolutions in the increasingly complex Y and mtDNA phylogenies, to avoid false exclusions. Incorporation of X-SNPs allows tracing of X-chromosome segments within families. These different approaches can add value to kinship estimation, but some require simpler bioinformatic interfaces to make them more widely accessible in practice, and also access to appropriate allele frequency data to avoid problems associated with ancestry mis-specification.

Perspectives and opportunities in forensic human, animal, and plant integrative genomics in the Pangenome era

The Human Pangenome Reference Consortium, the Chinese Pangenome Consortium, and other plant and animal pangenome projects have announced the completion of pilot work aimed at constructing high-quality, haplotype-resolved reference graph genomes representative of global ethno-linguistically different populations or different plant and animal species. These graph-based, gapless pangenome references, which are enriched in terms of genomic diversity, completeness, and contiguity, have the potential for enhancing long-read sequencing (LRS)-based genomic research, as well as improving mappability and variant genotyping on traditional short-read sequencing platforms. We comprehensively discuss the advancements in pangenome-based genomic integrative genomic discoveries across forensic-related species (humans, animals, and plants) and summarize their applications in variant identification and forensic genomics, epigenetics, transcriptomics, and microbiome research. Recent developments in multiplexed array sequencing have introduced a highly efficient and programmable technique to overcome the limitations of short forensic marker lengths in LRS platforms. This technique enables the concatenation of short RNA transcripts and DNA fragments into LRS-optimal molecules for sequencing, assembly, and genotyping. The integration of new pangenome reference coordinates and corresponding computational algorithms will benefit forensic integrative genomics by facilitating new marker identification, accurate genotyping, high-resolution panel development, and the updating of statistical algorithms. This review highlights the necessity of integrating LRS-based platforms, pangenome-based study designs, and graph-based pangenome references in short-read mapping and LRS-based innovations to achieve precision forensic science.

Estimation of population-specific values of theta for PowerPlex Y23 profiles

We examine 31,011 PPY23 profiles at the population, metapopulation and world levels. Most haplotypes appear only once but a few have higher counts, including a set of 23 matching profiles in Delhi, India and a set of 16 matching profiles in Burkina Faso with one additional matching American African profile. We estimate FSTvalues to be used as "theta" (θ) in match probability calculations, following the method we used in our earlier survey of autosomal STR data. Match probability estimates using FˆST or the κ method of Brenner for a previously unseen profile are similar but differ for any profile previously seen.

North Pontic crossroads: Mobility in Ukraine from the Bronze Age to the early modern period

The North Pontic region, which encompasses present-day Ukraine, was a crossroads of migration, connecting the vast Eurasian Steppe with Central Europe. We generated shotgun-sequenced genomic data for 91 individuals dating from around 7000 BCE to 1800 CE to study migration and mobility history in the region, with a particular focus on historically attested migrating groups during the Iron Age and the medieval period. We infer a high degree of temporal heterogeneity in ancestry, with fluctuating genetic affinities to different present-day Eurasian groups. We also infer high heterogeneity in ancestry within geographically, culturally, and socially defined groups. Despite this, we find that ancestry components which are widespread in Eastern and Central Europe have been present in the Ukraine region since the Bronze Age. In short, our study reveals a diverse range of ancestries in the Ukraine region through time as a result of frequent movements, assimilation, and contacts.

Characterization of Y chromosome diversity in newfoundland and labrador: evidence for a structured founding population

The population of Newfoundland and Labrador (NL) is largely derived from settlers who migrated primarily from England and Ireland in the 1700s-1800s. Previously described as an isolated founder population, based on historical and demographic studies, data on the genetic ancestry of this population remains fragmentary. Here we describe the largest investigation of patrilineal ancestry in NL. To determine the paternal genetic structure of the population, 1,110 Y chromosomes from an NL-based cohort were analyzed using 5,761 Y-specific SNPs. We identified 160 distinct terminal haplogroups, the majority of which (71.4%) belong to the R1b haplogroup. When compared with global reference populations, the NL population haplogroup composition and frequencies primarily resemble those observed in English and Irish ancestral source populations. There is also evidence of genetic contributions from Basque, French, Portuguese, and Spanish fishermen and early settlers who frequented NL. Interestingly, the observed population structure shows geographical and religious clustering that can be associated with the settlement of the ancestral source populations from predominantly Protestant, England, and Catholic, Ireland respectively. For example, the R1b-M222 haplogroup, seen in people of Irish descent, is found clustered in the Irish-settled Southeast region of NL. The clustering and expansion of Y haplogroups in conjunction with the geographical and religious clusters illustrate that limited subsequent in-migration, geographic isolation, and societal factors have contributed to the genetic substructure of the NL population and its designation as a founder population.

Genomic insights into Neolithic founding paternal lineages around the Qinghai-Xizang Plateau using integrated YanHuang resource

Indigenous populations of the Qinghai-Xizang Plateau exhibit unique high-altitude adaptations, especially within Tibeto-Burman (TB) groups. However, the paternal genetic heritage of eastern Plateau regions remains less explored. We present one integrative Y chromosome dataset of 9,901 modern and ancient individuals, including whole Y chromosome sequences from 1,297 individuals and extensive Y-SNP/STR genotype data. We reveal the Paleolithic common origin and following divergence of Qinghai-Xizang Plateau ancestors from East Asian lowlands, marked by subsequent isolation and Holocene expansion involving local hunter-gatherers and millet-farming communities. We identified two key TB-related founding lineages, D-Z31591 and O-CTS4658, which underwent significant expansions around 5,000 years ago on the Qinghai-Xizang Plateau and its eastern Tibetan-Yi Corridor. The genetic legacy of these TB lineages highlights crucial migration pathways linking the Plateau and lowland southwestern China. Our findings align paternal genetic structures with East Asian geography and linguistic groups, underscoring the utility of Y chromosome analyses in unraveling complex paternal histories.

Nanopore sequencing with unique molecular identifiers enables accurate mutation analysis and haplotyping in the complex lipoprotein(a) KIV-2 VNTR

BACKGROUND

Repetitive genome regions, such as variable number of tandem repeats (VNTR) or short tandem repeats (STR), are major constituents of the uncharted dark genome and evade conventional sequencing approaches. The protein-coding LPA kringle IV type-2 (KIV-2) VNTR (5.6 kb per unit, 1-40 units per allele) is a medically highly relevant example with a particularly intricate structure, multiple haplotypes, intragenic homologies, and an intra-VNTR STR. It is the primary regulator of plasma lipoprotein(a) [Lp(a)] concentrations, an important cardiovascular risk factor. Lp(a) concentrations vary widely between individuals and ancestries. Multiple variants and functional haplotypes in the LPA gene and especially in the KIV-2 VNTR strongly contribute to this variance.

METHODS

We evaluated the performance of amplicon-based nanopore sequencing with unique molecular identifiers (UMI-ONT-Seq) for SNP detection, haplotype mapping, VNTR unit consensus sequence generation, and copy number estimation via coverage-corrected haplotypes quantification in the KIV-2 VNTR. We used 15 human samples and low-level mixtures (0.5 to 5%) of KIV-2 plasmids as a validation set. We then applied UMI-ONT-Seq to extract KIV-2 VNTR haplotypes in 48 multi-ancestry 1000 Genome samples and analyzed at scale a poorly characterized STR within the KIV-2 VNTR.

RESULTS

UMI-ONT-Seq detected KIV-2 SNPs down to 1% variant level with high sensitivity, specificity, and precision (0.977 ± 0.018; 1.000 ± 0.0005; 0.993 ± 0.02) and accurately retrieved the full-length haplotype of each VNTR unit. Human variant levels were highly correlated with next-generation sequencing (R2 = 0.983) without bias across the whole variant level range. Six reads per UMI produced sequences of each KIV-2 unit with Q40 quality. The KIV-2 repeat number determined by coverage-corrected unique haplotype counting was in close agreement with droplet digital PCR (ddPCR), with 70% of the samples falling even within the narrow confidence interval of ddPCR. We then analyzed 62,679 intra-KIV-2 STR sequences and explored KIV-2 SNP haplotype patterns across five ancestries.

CONCLUSIONS

UMI-ONT-Seq accurately retrieves the SNP haplotype and precisely quantifies the VNTR copy number of each repeat unit of the complex KIV-2 VNTR region across multiple ancestries. This study utilizes the KIV-2 VNTR, presenting a novel and potent tool for comprehensive characterization of medically relevant complex genome regions at scale.

Low Genetic Impact of the Roman Occupation of Britain in Rural Communities

The Roman period saw the empire expand across Europe and the Mediterranean, including much of what is today Great Britain. While there is written evidence of high mobility into and out of Britain for administrators, traders, and the military, the impact of imperialism on local, rural population structure, kinship, and mobility is invisible in the textual record. The extent of genetic change that occurred in Britain during the Roman military occupation remains underexplored. Here, using genome-wide data from 52 ancient individuals from eight sites in Cambridgeshire covering the period of Roman occupation, we show low levels of genetic ancestry differentiation between Romano-British sites and indications of larger populations than in the Bronze Age and Neolithic. We find no evidence of long-distance migration from elsewhere in the Empire, though we do find one case of possible temporary mobility within a family unit during the Late Romano-British period. We also show that the present-day patterns of genetic ancestry composition in Britain emerged after the Roman period.

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.

Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

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.

Inferences about the population history of Rangifer tarandus from Y chromosome and mtDNA phylogenies

Reindeer, called caribou in North America, has a circumpolar distribution and all extant populations belong to the same species (Rangifer tarandus). It has survived the Holocene thanks to its immense adaptability and successful coexistence with humans in different forms of hunting and herding cultures. Here, we examine the paternal and maternal history of Rangifer based on robust Y-chromosomal and mitochondrial DNA (mtDNA) trees representing Eurasian tundra reindeer, Finnish forest reindeer, Svalbard reindeer, Alaska tundra caribou, and woodland caribou. We first assembled Y-chromosomal contigs, representing 1.3 Mb of single-copy Y regions. Based on 545 Y-chromosomal and 458 mtDNA SNPs defined in 55 males, maximum parsimony trees were created. We observed two well separated clades in both phylogenies: the "EuroBeringian clade" formed by animals from Arctic Islands, Eurasia, and a few from North America and the "North American clade" formed only by caribou from North America. The time calibrated Y tree revealed an expansion and dispersal of lineages across continents after the Last Glacial Maximum. We show for the first time unique paternal lineages in Svalbard reindeer and Finnish forest reindeer and reveal a circumscribed Y haplogroup in Fennoscandian tundra reindeer. The Y chromosome in domesticated reindeer is markedly diverse indicating that several male lineages have undergone domestication and less intensive selection on males. This study places R. tarandus onto the list of species with resolved Y and mtDNA phylogenies and builds the basis for studies of the distribution and origin of paternal and maternal lineages in the future.

Evolutionary profiles and complex admixture landscape in East Asia: New insights from modern and ancient Y chromosome variation perspectives

Human Y-chromosomes are characterized by nonrecombination and uniparental inheritance, carrying traces of human history evolution and admixture. Large-scale population-specific genomic sources based on advanced sequencing technologies have revolutionized our understanding of human Y chromosome diversity and its anthropological and forensic applications. Here, we reviewed and meta-analyzed the Y chromosome genetic diversity of modern and ancient people from China and summarized the patterns of founding lineages of spatiotemporally different populations associated with their origin, expansion, and admixture. We emphasized the strong association between our identified founding lineages and language-related human dispersal events correlated with the Sino-Tibetan, Altaic, and southern Chinese multiple-language families related to the Hmong-Mien, Tai-Kadai, Austronesian, and Austro-Asiatic languages. We subsequently summarize the recent advances in translational applications in forensic and anthropological science, including paternal biogeographical ancestry inference (PBGAI), surname investigation, and paternal history reconstruction. Whole-Y sequencing or high-resolution panels with high coverage of terminal Y chromosome lineages are essential for capturing the genomic diversity of ethnolinguistically diverse East Asians. Generally, we emphasized the importance of including more ethnolinguistically diverse, underrepresented modern and spatiotemporally different ancient East Asians in human genetic research for a comprehensive understanding of the paternal genetic landscape of East Asians with a detailed time series and for the reconstruction of a reference database in the PBGAI, even including new technology innovations of Telomere-to-Telomere (T2T) for new genetic variation discovery.

Patrilineal segmentary systems provide a peaceful explanation for the post-Neolithic Y-chromosome bottleneck

Studies have found a pronounced decline in male effective population sizes worldwide around 3000-5000 years ago. This bottleneck was not observed for female effective population sizes, which continued to increase over time. Until now, this remarkable genetic pattern was interpreted as the result of an ancient structuring of human populations into patrilineal groups (gathering closely related males) violently competing with each other. In this scenario, violence is responsible for the repeated extinctions of patrilineal groups, leading to a significant reduction in male effective population size. Here, we propose an alternative hypothesis by modelling a segmentary patrilineal system based on anthropological literature. We show that variance in reproductive success between patrilineal groups, combined with lineal fission (i.e., the splitting of a group into two new groups of patrilineally related individuals), can lead to a substantial reduction in the male effective population size without resorting to the violence hypothesis. Thus, a peaceful explanation involving ancient changes in social structures, linked to global changes in subsistence systems, may be sufficient to explain the reported decline in Y-chromosome diversity.

The maternal U1 haplogroup in the Koraga tribe as a correlate of their North Dravidian linguistic affinity

Introduction: The Koraga tribe are an isolated endogamous tribal group found in the southwest coastal region of India. The Koraga language shares inherited grammatical features with North Dravidian languages. To seek a possible genetic basis for this exceptionality and understand the maternal lineage pattern, we have aimed to reconstruct the inter-population and intra-population relationships of the Koraga tribal population by using mtDNA markers for the hypervariable regions along with a partial coding region sequence analysis. Methods and Results: Amongst the 96 individuals studied, we observe 11 haplogroups, of which a few are shared and others are unique to the clans Soppu, Oṇṭi and Kuṇṭu. In addition to several deep rooted Indian-specific lineages of macrohaplogroups M and U, we observe a high frequency of the U1 lineage (∼38%), unique to the Koraga. A Bayesian analysis of the U1 clade shows that the Koraga tribe share their maternal lineage with ancestral populations of the Caucasus at the cusp of the Last Glacial Maximum. Discussion: Our study suggests that the U1 lineage found in the Indian subcontinent represents a remnant of a post-glacial dispersal. The presence of West Asian U1 when viewed along with historical linguistics leads us to hypothesise that Koraga represents a mother tongue retained by a vanquished population group that fled southward at the demise of the Indus civilisation as opposed to a father tongue, associated with a particular paternal lineage.

Differentiated genomic footprints suggest isolation and long-distance migration of Hmong-Mien populations

BACKGROUND

The underrepresentation of Hmong-Mien (HM) people in Asian genomic studies has hindered our comprehensive understanding of the full landscape of their evolutionary history and complex trait architecture. South China is a multi-ethnic region and indigenously settled by ethnolinguistically diverse HM, Austroasiatic (AA), Tai-Kadai (TK), Austronesian (AN), and Sino-Tibetan (ST) people, which is regarded as East Asia's initial cradle of biodiversity. However, previous fragmented genetic studies have only presented a fraction of the landscape of genetic diversity in this region, especially the lack of haplotype-based genomic resources. The deep characterization of demographic history and natural-selection-relevant genetic architecture of HM people was necessary.

RESULTS

We reported one HM-specific genomic resource and comprehensively explored the fine-scale genetic structure and adaptative features inferred from the genome-wide SNP data of 440 HM individuals from 33 ethnolinguistic populations, including previously unreported She. We identified solid genetic differentiation between HM people and Han Chinese at 7.64‒15.86 years ago (kya) and split events between southern Chinese inland (Miao/Yao) and coastal (She) HM people in the middle Bronze Age period and the latter obtained more gene flow from Ancient Northern East Asians. Multiple admixture models further confirmed that extensive gene flow from surrounding ST, TK, and AN people entangled in forming the gene pool of Chinese coastal HM people. Genetic findings of isolated shared unique ancestral components based on the sharing alleles and haplotypes deconstructed that HM people from the Yungui Plateau carried the breadth of previously unknown genomic diversity. We identified a direct and recent genetic connection between Chinese inland and Southeast Asian HM people as they shared the most extended identity-by-descent fragments, supporting the long-distance migration hypothesis. Uniparental phylogenetic topology and network-based phylogenetic relationship reconstruction found ancient uniparental founding lineages in southwestern HM people. Finally, the population-specific biological adaptation study identified the shared and differentiated natural selection signatures among inland and coastal HM people associated with physical features and immune functions. The allele frequency spectrum of cancer susceptibility alleles and pharmacogenomic genes showed significant differences between HM and northern Chinese people.

CONCLUSIONS

Our extensive genetic evidence combined with the historical documents supported the view that ancient HM people originated from the Yungui regions associated with ancient "Three-Miao tribes" descended from the ancient Daxi-Qujialing-Shijiahe people. Then, some have recently migrated rapidly to Southeast Asia, and some have migrated eastward and mixed respectively with Southeast Asian indigenes, Liangzhu-related coastal ancient populations, and incoming southward ST people. Generally, complex population migration, admixture, and adaptation history contributed to the complicated patterns of population structure of geographically diverse HM people.

Genetic history of Cambridgeshire before and after the Black Death

The extent of the devastation of the Black Death pandemic (1346-1353) on European populations is known from documentary sources and its bacterial source illuminated by studies of ancient pathogen DNA. What has remained less understood is the effect of the pandemic on human mobility and genetic diversity at the local scale. Here, we report 275 ancient genomes, including 109 with coverage >0.1×, from later medieval and postmedieval Cambridgeshire of individuals buried before and after the Black Death. Consistent with the function of the institutions, we found a lack of close relatives among the friars and the inmates of the hospital in contrast to their abundance in general urban and rural parish communities. While we detect long-term shifts in local genetic ancestry in Cambridgeshire, we find no evidence of major changes in genetic ancestry nor higher differentiation of immune loci between cohorts living before and after the Black Death.

Temporal Variation in Introgressed Segments' Length Statistics Computed from a Limited Number of Ancient Genomes Sheds Light on Past Admixture Pulses

Hybridization is recognized as an important evolutionary force, but identifying and timing admixture events between divergent lineages remain a major aim of evolutionary biology. While this has traditionally been done using inferential tools on contemporary genomes, the latest advances in paleogenomics have provided a growing wealth of temporally distributed genomic data. Here, we used individual-based simulations to generate chromosome-level genomic data for a 2-population system and described temporal neutral introgression patterns under a single- and 2-pulse admixture model. We computed 6 summary statistics aiming to inform the timing and number of admixture pulses between interbreeding entities: lengths of introgressed sequences and their variance within genomes, as well as genome-wide introgression proportions and related measures. The first 2 statistics could confidently be used to infer interlineage hybridization history, peaking at the beginning and shortly after an admixture pulse. Temporal variation in introgression proportions and related statistics provided more limited insights, particularly when considering their application to ancient genomes still scant in number. Lastly, we computed these statistics on Homo sapiens paleogenomes and successfully inferred the hybridization pulse from Neanderthal that occurred approximately 40 to 60 kya. The scarce number of genomes dating from this period prevented more precise inferences, but the accumulation of paleogenomic data opens promising perspectives as our approach only requires a limited number of ancient genomes.

Genetic sexing of subadult skeletal remains

When subadult skeletons need to be identified, biological sex diagnosis is one of the first steps in the identification process. Sex assessment of subadults using morphological features is unreliable, and molecular genetic methods were applied in this study. Eighty-three ancient skeletons were used as models for poorly preserved DNA. Three sex-informative markers on the Y and X chromosome were used for sex identification: a qPCR test using the PowerQuant Y target included in PowerQuant System (Promega), the amelogenin test included in ESI 17 Fast STR kit (Promega), and a Y-STR amplification test using the PowerPlex Y-23 kit (Promega). Sex was successfully determined in all but five skeletons. Successful PowerQuant Y-target, Y-amelogenin, and Y-chromosomal STR amplifications proved the presence of male DNA in 35 skeletons, and in 43 subadults female sex was established. No match was found between the genetic profiles of subadult skeletons, and the elimination database and negative control samples produced no profiles, indicating no contamination issue. Our study shows that genetic sex identification is a very successful approach for biological sexing of subadult skeletons whose sex cannot be assessed by anthropological methods. The results of this study are applicable for badly preserved subadult skeletons from routine forensic casework.

Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history

The Two Layer hypothesis is fast becoming the favoured narrative describing East Asian population history. Under this model, hunter-gatherer groups who initially peopled East Asia via a route south of the Himalayas were assimilated by agriculturalist migrants who arrived via a northern route across Eurasia. A lack of ancient samples from tropical East Asia limits the resolution of this model. We consider insight afforded by patterns of variation within the human pathogen Mycobacterium tuberculosis (Mtb) by analysing its phylogeographic signatures jointly with the human Y-chromosome. We demonstrate the Y-chromosome lineages enriched in the traditionally hunter-gatherer groups associated with East Asia's first layer of peopling to display deep roots, low long-term effective population size, and diversity patterns consistent with a southern entry route. These characteristics mirror those of the evolutionarily ancient Mtb lineage 1. The remaining East Asian Y-chromosome lineage is almost entirely absent from traditionally hunter-gatherer groups and displays spatial and temporal characteristics which are incompatible with a southern entry route, and which link it to the development of agriculture in modern-day China. These characteristics mirror those of the evolutionarily modern Mtb lineage 2. This model paves the way for novel host-pathogen coevolutionary research hypotheses in East Asia.

Forensic biogeographical ancestry inference: recent insights and current trends

BACKGROUND

As a powerful complement to the paradigmatic DNA profiling strategy, biogeographical ancestry inference (BGAI) plays a significant part in human forensic investigation especially when a database hit or eyewitness testimony are not available. It indicates one's biogeographical profile based on known population-specific genetic variations, and thus is crucial for guiding authority investigations to find unknown individuals. Forensic biogeographical ancestry testing exploits much of the recent advances in the understanding of human genomic variation and improving of molecular biology.

OBJECTIVE

In this review, recent development of prospective ancestry informative markers (AIMs) and the statistical approaches of inferring biogeographic ancestry from AIMs are elucidated and discussed.

METHODS

We highlight the research progress of three potential AIMs (i.e., single nucleotide polymorphisms, microhaplotypes, and Y or mtDNA uniparental markers) and discuss the prospects and challenges of two methods that are commonly used in BGAI.

CONCLUSION

While BGAI for forensic purposes has been thriving in recent years, important challenges, such as ethics and responsibilities, data completeness, and ununified standards for evaluation, remain for the use of biogeographical ancestry information in human forensic investigations. To address these issues and fully realize the value of BGAI in forensic investigation, efforts should be made not only by labs/institutions around the world independently, but also by inter-lab/institution collaborations.

Benchmarking of human Y-chromosomal haplogroup classifiers with whole-genome and whole-exome sequence data

In anthropological, medical, and forensic studies, the nonrecombinant region of the human Y chromosome (NRY) enables accurate reconstruction of pedigree relationships and retrieval of ancestral information. Using high-throughput sequencing (HTS) data, we present a benchmarking analysis of command-line tools for NRY haplogroup classification. The evaluation was performed using paired Illumina data from whole-genome sequencing (WGS) and whole-exome sequencing (WES) experiments from 50 unrelated donors. Additionally, as a validation, we also used paired WGS/WES datasets of 54 individuals from the 1000 Genomes Project. Finally, we evaluated the tools on data from third-generation HTS obtained from a subset of donors and one reference sample. Our results show that WES, despite typically offering less genealogical resolution than WGS, is an effective method for determining the NRY haplogroup. Y-LineageTracker and Yleaf showed the highest accuracy for WGS data, classifying precisely 98% and 96% of the samples, respectively. Yleaf outperforms all benchmarked tools in the WES data, classifying approximately 90% of the samples. Yleaf, Y-LineageTracker, and pathPhynder can correctly classify most samples (88%) sequenced with third-generation HTS. As a result, Yleaf provides the best performance for applications that use WGS and WES. Overall, our study offers researchers with a guide that allows them to select the most appropriate tool to analyze the NRY region using both second- and third-generation HTS data.

South Asian maternal and paternal lineages in southern Thailand and the role of sex-biased admixture

Previous genome-wide studies have reported South Asian (SA) ancestry in several Mainland Southeast Asian (MSEA) populations; however, additional details concerning population history, in particular the role of sex-specific aspects of the SA admixture in MSEA populations can be addressed with uniparental markers. Here, we generated ∼2.3 mB sequences of the male-specific portions of the Y chromosome (MSY) of a Tai-Kadai (TK)-speaking Southern Thai group (SouthernThai_TK), and complete mitochondrial (mtDNA) genomes of the SouthernThai_TK and an Austronesian (AN)-speaking Southern Thai (SouthernThai_AN) group. We identified new mtDNA haplogroups, e.g. Q3, E1a1a1, B4a1a and M7c1c3 that have not previously reported in Thai populations, but are frequent in Island Southeast Asia and Oceania, suggesting interactions between MSEA and these regions. SA prevalent mtDNA haplogroups were observed at frequencies of ~35-45% in the Southern Thai groups; both of them showed more genetic relatedness to Austroasiatic (AA) speaking Mon than to any other group. For MSY, SouthernThai_TK had ~35% SA prevalent haplogroups and exhibited closer genetic affinity to Central Thais. We also analyzed published data from other MSEA populations and observed SA ancestry in some additional MSEA populations that also reflects sex-biased admixture; in general, most AA- and AN-speaking groups in MSEA were closer to SA than to TK groups based on mtDNA, but the opposite pattern was observed for the MSY. Overall, our results of new genetic lineages and sex-biased admixture from SA to MSEA groups attest to the additional value that uniparental markers can add to studies of genome-wide variation.

Assembly of 43 human Y chromosomes reveals extensive complexity and variation

The prevalence of highly repetitive sequences within the human Y chromosome has prevented its complete assembly to date1 and led to its systematic omission from genomic analyses. Here we present de novo assemblies of 43 Y chromosomes spanning 182,900 years of human evolution and report considerable diversity in size and structure. Half of the male-specific euchromatic region is subject to large inversions with a greater than twofold higher recurrence rate compared with all other chromosomes2. Ampliconic sequences associated with these inversions show differing mutation rates that are sequence context dependent, and some ampliconic genes exhibit evidence for concerted evolution with the acquisition and purging of lineage-specific pseudogenes. The largest heterochromatic region in the human genome, Yq12, is composed of alternating repeat arrays that show extensive variation in the number, size and distribution, but retain a 1:1 copy-number ratio. Finally, our data suggest that the boundary between the recombining pseudoautosomal region 1 and the non-recombining portions of the X and Y chromosomes lies 500 kb away from the currently established1 boundary. The availability of fully sequence-resolved Y chromosomes from multiple individuals provides a unique opportunity for identifying new associations of traits with specific Y-chromosomal variants and garnering insights into the evolution and function of complex regions of the human genome.

The complete sequence of a human Y chromosome

The human Y chromosome has been notoriously difficult to sequence and assemble because of its complex repeat structure that includes long palindromes, tandem repeats and segmental duplications1-3. As a result, more than half of the Y chromosome is missing from the GRCh38 reference sequence and it remains the last human chromosome to be finished4,5. Here, the Telomere-to-Telomere (T2T) consortium presents the complete 62,460,029-base-pair sequence of a human Y chromosome from the HG002 genome (T2T-Y) that corrects multiple errors in GRCh38-Y and adds over 30 million base pairs of sequence to the reference, showing the complete ampliconic structures of gene families TSPY, DAZ and RBMY; 41 additional protein-coding genes, mostly from the TSPY family; and an alternating pattern of human satellite 1 and 3 blocks in the heterochromatic Yq12 region. We have combined T2T-Y with a previous assembly of the CHM13 genome4 and mapped available population variation, clinical variants and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes.

Comprehensive insights into the genetic background of Chinese populations using Y chromosome markers

China is located in East Asia. With a high genetic and cultural diversity, human migration in China has always been a hot topic of genetics research. To explore the origins and migration routes of Chinese males, 3333 Chinese individuals (Han, Hui, Mongolia, Yi and Kyrgyz) with 27 Y-STRs and 143 Y-SNPs from published literature were analysed. Our data showed that there are five dominant haplogroups (O2-M122, O1-F265, C-M130, N-M231, R-M207) in China. Combining analysis of haplogroup frequencies, geographical positions and time with the most recent common ancestor (TMRCA), we found that haplogroups C-M130, N-M231 and R1-M173 and O1a-M175 probably migrated into China via the northern route. Interestingly, we found that haplogroup C*-M130 in China may originate in South Asia, whereas the major subbranches C2a-L1373 and C2b-F1067 migrated from northern China. The results of BATWING showed that the common ancestry of Y haplogroup in China can be traced back to 17 000 years ago, which was concurrent with global temperature increases after the Last Glacial Maximum.

Genomic history of coastal societies from eastern South America

Sambaqui (shellmound) societies are among the most intriguing archaeological phenomena in pre-colonial South America, extending from approximately 8,000 to 1,000 years before present (yr BP) across 3,000 km on the Atlantic coast. However, little is known about their connection to early Holocene hunter-gatherers, how this may have contributed to different historical pathways and the processes through which late Holocene ceramists came to rule the coast shortly before European contact. To contribute to our understanding of the population history of indigenous societies on the eastern coast of South America, we produced genome-wide data from 34 ancient individuals as early as 10,000 yr BP from four different regions in Brazil. Early Holocene hunter-gatherers were found to lack shared genetic drift among themselves and with later populations from eastern South America, suggesting that they derived from a common radiation and did not contribute substantially to later coastal groups. Our analyses show genetic heterogeneity among contemporaneous Sambaqui groups from the southeastern and southern Brazilian coast, contrary to the similarity expressed in the archaeological record. The complex history of intercultural contact between inland horticulturists and coastal populations becomes genetically evident during the final horizon of Sambaqui societies, from around 2,200 yr BP, corroborating evidence of cultural change.

Y chromosome sequence and epigenomic reconstruction across human populations

Recent advances in long-read sequencing technologies have allowed the generation and curation of more complete genome assemblies, enabling the analysis of traditionally neglected chromosomes, such as the human Y chromosome (chrY). Native DNA was sequenced on a MinION Oxford Nanopore Technologies sequencing device to generate genome assemblies for seven major chrY human haplogroups. We analyzed and compared the chrY enrichment of sequencing data obtained using two different selective sequencing approaches: adaptive sampling and flow cytometry chromosome sorting. We show that adaptive sampling can produce data to create assemblies comparable to chromosome sorting while being a less expensive and time-consuming technique. We also assessed haplogroup-specific structural variants, which would be otherwise difficult to study using short-read sequencing data only. Finally, we took advantage of this technology to detect and profile epigenetic modifications among the considered haplogroups. Altogether, we provide a framework to study complex genomic regions with a simple, fast, and affordable methodology that could be applied to larger population genomics datasets.

Demographic history and genetic structure in pre-Hispanic Central Mexico

Aridoamerica and Mesoamerica are two distinct cultural areas in northern and central Mexico, respectively, that hosted numerous pre-Hispanic civilizations between 2500 BCE and 1521 CE. The division between these regions shifted southward because of severe droughts ~1100 years ago, which allegedly drove a population replacement in central Mexico by Aridoamerican peoples. In this study, we present shotgun genome-wide data from 12 individuals and 27 mitochondrial genomes from eight pre-Hispanic archaeological sites across Mexico, including two at the shifting border of Aridoamerica and Mesoamerica. We find population continuity that spans the climate change episode and a broad preservation of the genetic structure across present-day Mexico for the past 2300 years. Lastly, we identify a contribution to pre-Hispanic populations of northern and central Mexico from two ancient unsampled "ghost" populations.

Multiple founding paternal lineages inferred from the newly-developed 639-plex Y-SNP panel suggested the complex admixture and migration history of Chinese people

BACKGROUND

Non-recombining regions of the Y-chromosome recorded the evolutionary traces of male human populations and are inherited haplotype-dependently and male-specifically. Recent whole Y-chromosome sequencing studies have identified previously unrecognized population divergence, expansion and admixture processes, which promotes a better understanding and application of the observed patterns of Y-chromosome genetic diversity.

RESULTS

Here, we developed one highest-resolution Y-chromosome single nucleotide polymorphism (Y-SNP) panel targeted for uniparental genealogy reconstruction and paternal biogeographical ancestry inference, which included 639 phylogenetically informative SNPs. We genotyped these loci in 1033 Chinese male individuals from 33 ethnolinguistically diverse populations and identified 256 terminal Y-chromosomal lineages with frequency ranging from 0.0010 (singleton) to 0.0687. We identified six dominant common founding lineages associated with different ethnolinguistic backgrounds, which included O2a2b1a1a1a1a1a1a1-M6539, O2a1b1a1a1a1a1a1-F17, O2a2b1a1a1a1a1b1a1b-MF15397, O2a2b2a1b1-A16609, O1b1a1a1a1b2a1a1-F2517, and O2a2b1a1a1a1a1a1-F155. The AMOVA and nucleotide diversity estimates revealed considerable differences and high genetic diversity among ethnolinguistically different populations. We constructed one representative phylogenetic tree among 33 studied populations based on the haplogroup frequency spectrum and sequence variations. Clustering patterns in principal component analysis and multidimensional scaling results showed a genetic differentiation between Tai-Kadai-speaking Li, Mongolic-speaking Mongolian, and other Sinitic-speaking Han Chinese populations. Phylogenetic topology inferred from the BEAST and Network relationships reconstructed from the popART further showed the founding lineages from culturally/linguistically diverse populations, such as C2a/C2b was dominant in Mongolian people and O1a/O1b was dominant in island Li people. We also identified many lineages shared by more than two ethnolinguistically different populations with a high proportion, suggesting their extensive admixture and migration history.

CONCLUSIONS

Our findings indicated that our developed high-resolution Y-SNP panel included major dominant Y-lineages of Chinese populations from different ethnic groups and geographical regions, which can be used as the primary and powerful tool for forensic practice. We should emphasize the necessity and importance of whole sequencing of more ethnolinguistically different populations, which can help identify more unrecognized population-specific variations for the promotion of Y-chromosome-based forensic applications.

MPYS Modulates Fatty Acid Metabolism and Immune Tolerance at Homeostasis Independent of Type I IFNs

MPYS/STING (stimulator of IFN genes) senses cyclic dinucleotides (CDNs), generates type I IFNs, and plays a critical role in infection, inflammation, and cancer. In this study, analyzing genotype and haplotype data from the 1000 Genomes Project, we found that the R71H-G230A-R293Q (HAQ) MPYS allele frequency increased 57-fold in East Asians compared with sub-Saharan Africans. Meanwhile, the G230A-R293Q (AQ) allele frequency decreased by 98% in East Asians compared with sub-Saharan Africans. We propose that the HAQ and AQ alleles underwent a natural selection during the out-of-Africa migration. We used mouse models of HAQ and AQ to investigate the underlying mechanism. We found that the mice carrying the AQ allele, which disappeared in East Asians, had normal CDN-type I IFN responses. Adult AQ mice, however, had less fat mass than did HAQ or wild-type mice on a chow diet. AQ epididymal adipose tissue had increased regulatory T cells and M2 macrophages with protein expression associated with enhanced fatty acid oxidation. Conditional knockout mice and adoptive cell transfer indicate a macrophage and regulatory T cell-intrinsic role of MPYS in fatty acid metabolism. Mechanistically, AQ/IFNAR1-/- mice had a similar lean phenotype as for the AQ mice. MPYS intrinsic tryptophan fluorescence revealed that the R71H change increased MPYS hydrophilicity. Lastly, we found that the second transmembrane (TM) and the TM2-TM3 linker region of MPYS interact with activated fatty acid, fatty acyl-CoA. In summary, studying the evolution of the human MPYS gene revealed an MPYS function in modulating fatty acid metabolism that may be critical during the out-of-Africa migration.

Genetic diversity and forensic application of Y-filer STRs in four major ethnic groups of Pakistan

17 Y-chromosomal STRs which are part of the Yfiler Amplification Kit were investigated in 493 unrelated Pakistani individuals belonging to the Punjabi, Sindhi, Baloch, and Pathan ethnic groups. We have assessed the forensic parameters and population genetic structure for each group. Among the 493 unrelated individuals from four ethnic groups (128 Baloch, 122 Pathan, 108 Punjabi, and 135 Sindhi), 82 haplotypes were observed with haplotype diversity (HD) of 0.9906 in Baloch, 102 haplotypes with HD value of 0.9957 in Pathans, 80 haplotypes with HD value of 0.9924 in Punjabi, and 105 haplotypes with HD value of 0.9945 in the Sindhi population. The overall gene diversity for Baloch, Pathan, Punjabi, and Sindhi populations was 0.6367, 0.6479, 0.6657, and 0.6112, respectively. The results had shown us that Pakistani populations do not have a unique set of genes but share the genetic affinity with regional (Central Asia and Northern India) populations. The observed low gene diversity (heterozygosity) values may be because of endogamy trends and this observation is equally supported by the results of forensic parameters which are mostly static across 4 combinations (minimal STRs, extended 11 Y-STRs, Powerplex 12 Y System, and Yfiler 17 Y-STRs) of STRs in these four populations.

The population genomic legacy of the second plague pandemic

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.¹ It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).² Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17^th-19^th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large F_ST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.

Tracing Bai-Yue Ancestry in Aboriginal Li People on Hainan Island

As the most prevalent aboriginal group on Hainan Island located between South China and the mainland of Southeast Asia, the Li people are believed to preserve some unique genetic information due to their isolated circumstances, although this has been largely uninvestigated. We performed the first whole-genome sequencing of 55 Hainan Li (HNL) individuals with high coverage (∼30-50×) to gain insight into their genetic history and potential adaptations. We identified the ancestry enriched in HNL (∼85%) is well preserved in present-day Tai-Kadai speakers residing in South China and North Vietnam, that is, Bai-Yue populations. A lack of admixture signature due to the geographical restriction exacerbated the bottleneck in the present-day HNL. The genetic divergence among Bai-Yue populations began ∼4,000-3,000 years ago when the proto-HNL underwent migration and the settling of Hainan Island. Finally, we identified signatures of positive selection in the HNL, some outstanding examples included FADS1 and FADS2 related to a diet rich in polyunsaturated fatty acids. In addition, we observed that malaria-driven selection had occurred in the HNL, with population-specific variants of malaria-related genes (e.g., CR1) present. Interestingly, HNL harbors a high prevalence of malaria leveraged gene variants related to hematopoietic function (e.g., CD3G) that may explain the high incidence of blood disorders such as B-cell lymphomas in the present-day HNL. The results have advanced our understanding of the genetic history of the Bai-Yue populations and have provided new insights into the adaptive scenarios of the Li people.

Limited Effect of Y Chromosome Variation on Coronary Artery Disease and Mortality in UK Biobank-Brief Report

BACKGROUND

The effect of genetic variation in the male-specific region of the Y chromosome (MSY) on coronary artery disease and cardiovascular risk factors has been disputed. In this study, we systematically assessed the association of MSY genetic variation on these traits using a kin-cohort analysis of family disease history in the largest sample to date.

METHODS

We tested 90 MSY haplogroups against coronary artery disease, hypertension, blood pressure, classical lipid levels, and all-cause mortality in up to 152 186 unrelated, genomically British individuals from UK Biobank. Unlike previous studies, we did not adjust for heritable lifestyle factors (to avoid collider bias) and instead adjusted for geographic variables and socioeconomic deprivation, given the link between MSY haplogroups and geography. For family history traits, subject MSY haplogroups were tested against father and mother disease as validation and negative control, respectively.

RESULTS

Our models find little evidence for an effect of any MSY haplogroup on cardiovascular risk in participants. Parental models confirm these findings.

CONCLUSIONS

Kin-cohort analysis of the Y chromosome uniquely allows for discoveries in subjects to be validated using family history data. Despite our large sample size, improved models, and parental validation, there is little evidence to suggest cardiovascular risk in UK Biobank is influenced by genetic variation in MSY.

Human Y chromosome sequences from Q Haplogroup reveal a South American settlement pre-18,000 years ago and a profound genomic impact during the Younger Dryas

The settlement of the Americas has been the focus of incessant debate for more than 100 years, and open questions regarding the timing and spatial patterns of colonization still remain today. Phylogenetic studies with complete human Y chromosome sequences are used as a highly informative tool to investigate the history of human populations in a given time frame. To study the phylogenetic relationships of Native American lineages and infer the settlement history of the Americas, we analyzed Y chromosome Q Haplogroup, which is a Pan-American haplogroup and represents practically all Native American lineages in Mesoamerica and South America. We built a phylogenetic tree for Q Haplogroup based on 102 whole Y chromosome sequences, of which 13 new Argentine sequences were provided by our group. Moreover, 1,072 new single nucleotide polymorphisms (SNPs) that contribute to its resolution and diversity were identified. Q-M848 is known to be the most frequent autochthonous sub-haplogroup of the Americas. The present is the first genomic study of Q Haplogroup in which current knowledge on Q-M848 sub-lineages is contrasted with the historical, archaeological and linguistic data available. The divergence times, spatial structure and the SNPs found here as novel for Q-Z780, a less frequent sub-haplogroup autochthonous of the Americas, provide genetic support for a South American settlement before 18,000 years ago. We analyzed how environmental events that occurred during the Younger Dryas period may have affected Native American lineages, and found that this event may have caused a substantial loss of lineages. This could explain the current low frequency of Q-Z780 (also perhaps of Q-F4674, a third possible sub-haplogroup autochthonous of the Americas). These environmental events could have acted as a driving force for expansion and diversification of the Q-M848 sub-lineages, which show a spatial structure that developed during the Younger Dryas period.

Assessing the utility of quantitative and qualitative metrics in the DNA quantification process of skeletal remains for autosomal and Y-chromosome STR amplification purposes

In historical cases, ancient DNA investigations and missing persons identification, teeth or bone samples are often the only and almost always the best biological material available for DNA typing. On the other hand, DNA obtained from bone material may be characterized by a high degradation index (DI) or its low content, or DNA tests cannot be repeated due to bone piece size limitation. That is often the effect of the environment in which the material was placed and the time during which exposure to unfavorable environmental factors took place. Therefore, it is very important to use appropriate procedures related to STR analysis. For our study, we selected 80 challenging bone samples. The amount of DNA was compared in qPCR using Quantifiler™ Trio DNA Quantification Kit and Investigator® Quantiplex® Pro RGQ. All qPCR results were confirmed by PCR-CE. The results of DNA concentrations and the assigned degradation index (DI) differed significantly within analyzed samples (~10%). Additionally, the Y-chromosome DI also differed from the autosomal DI in the samples. The difference in degradation indexes could explain the lower Y-chromosome amplification success rate compared to autosomal e.g. during human identification process. The results indicate that performing two DNA quantifications with the use of two different kits (primers sets) allows for a much more precise evaluation of the DNA quality and quantity in the isolate. We suggest that at least one of two suggested DNA concentration measurements should be based on an additional determination of the Y chromosome degradation index. Altogether, it allows for rational isolate management, especially when the volume is limited and the sample is unique.

Human Occupation of the North American Colorado Plateau ∼37,000 Years Ago

Calibrating human population dispersals across Earth’s surface is fundamental to assessing rates and timing of anthropogenic impacts and distinguishing ecological phenomena influenced by humans from those that were not. Here, we describe the Hartley mammoth locality, which dates to 38,900–36,250 cal BP by AMS 14C analysis of hydroxyproline from bone collagen. We accept the standard view that elaborate stone technology of the Eurasian Upper Paleolithic was introduced into the Americas by arrival of the Native American clade ∼16,000 cal BP. It follows that if older cultural sites exist in the Americas, they might only be diagnosed using nuanced taphonomic approaches. We employed computed tomography (CT and μCT) and other state-of-the-art methods that had not previously been applied to investigating ancient American sites. This revealed multiple lines of taphonomic evidence suggesting that two mammoths were butchered using expedient lithic and bone technology, along with evidence diagnostic of controlled (domestic) fire. That this may be an ancient cultural site is corroborated by independent genetic evidence of two founding populations for humans in the Americas, which has already raised the possibility of a dispersal into the Americas by people of East Asian ancestry that preceded the Native American clade by millennia. The Hartley mammoth locality thus provides a new deep point of chronologic reference for occupation of the Americas and the attainment by humans of a near-global distribution.

Substantial but Misunderstood Human Sexual Dimorphism Results Mainly From Sexual Selection on Males and Natural Selection on Females

Human sexual dimorphism has been widely misunderstood. A large literature has underestimated the effect of differences in body composition and the role of male contest competition for mates. It is often assumed that sexually dimorphic traits reflect a history of sexual selection, but natural selection frequently builds different phenotypes in males and females. The relatively small sex difference in stature (∼7%) and its decrease during human evolution have been widely presumed to indicate decreased male contest competition for mates. However, females likely increased in stature relative to males in order to successfully deliver large-brained neonates through a bipedally-adapted pelvis. Despite the relatively small differences in stature and body mass (∼16%), there are marked sex differences in body composition. Across multiple samples from groups with different nutrition, males typically have 36% more lean body mass, 65% more muscle mass, and 72% more arm muscle than women, yielding parallel sex differences in strength. These sex differences in muscle and strength are comparable to those seen in primates where sexual selection, arising from aggressive male mating competition, has produced high levels of dimorphism. Body fat percentage shows a reverse pattern, with females having ∼1.6 times more than males and depositing that fat in different body regions than males. We argue that these sex differences in adipose arise mainly from natural selection on women to accumulate neurodevelopmental resources.

High-Selectivity Single-Nucleotide Variant Capture Technology Based on the DNA Reaction Network

The extremely low abundance of circulating tumor DNA in blood samples has limited the development of liquid biopsy techniques for the early diagnosis of major diseases. In this study, we demonstrate a DRN-based screening technique, SCREEN, which achieves the specific capture and enrichment of low abundance SNV nucleic acid samples without selective amplification. The SCREEN technique achieved a 108-fold increase in the abundance of single-nucleotide variant (SNV) nucleic acids from highly homologous mixtures (from 0.01% to 1.08%) and has been shown to significantly increase the abundance of SNV nucleic acids from 0.1% to 51% further through two rounds of capture. As a highly effective pre-enrichment technique, SCREEN has demonstrated the ability to enhance NGS in detecting an ultralow abundance SNV nucleic acid powerfully and has high compatibility with existing molecular diagnostic methods.