A recent bottleneck of Y chromosome diversity coincides with a global change in culture

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.

Ancient genome of the Chinese Emperor Wu of Northern Zhou

Emperor Wu (, Wudi) of the Xianbei-led Northern Zhou dynasty, named Yuwen Yong (, 543-578 CE), was a highly influential emperor who reformed the system of regional troops, pacified the Turks, and unified the northern part of the country. His genetic profile and physical characteristics, including his appearance and potential diseases, have garnered significant interest from the academic community and the public. In this study, we have successfully generated a 0.343×-coverage genome of Wudi with 1,011,419 single-nucleotide polymorphisms (SNPs) on the 1240k panel. By analyzing pigmentation-relevant SNPs and conducting cranial CT-based facial reconstruction, we have determined that Wudi possessed a typical East or Northeast Asian appearance. Furthermore, pathogenic SNPs suggest Wudi faced an increased susceptibility to certain diseases, such as stroke. Wudi shared the closest genetic relationship with ancient Khitan and Heishui Mohe samples and modern Daur and Mongolian populations but also showed additional affinity with Yellow River (YR) farmers. We estimated that Wudi derived 61% of his ancestry from ancient Northeast Asians (ANAs) and nearly one-third from YR farmer-related groups. This can likely be attributed to continuous intermarriage between Xianbei royal families, and local Han aristocrats.1,2 Furthermore, our study has revealed genetic diversities among available ancient Xianbei individuals from different regions, suggesting that the formation of the Xianbei was a dynamic process influenced by admixture with surrounding populations.

Population genetics reveals new introgression in the nucleus herd of min pigs

OBJECTIVE

Min pigs are a unique genetic resource among local pig breeds in China. They have more excellent characteristics in cold and stress resistance, good meat quality, and a high reproductive rate. However, the genetic structure and driving factors remain unclear in the nucleus herd. In this study, the genetic diversity of Min pigs was studied to reveal the formation mechanism of its unique genetic structure. We hope to protect and develop the genetic resources of Min pigs.

METHODS

We analyzed different types of genes to identify the genetic structure and gene introgression pattern of Min pigs. The nuclear DNA dataset includes information on 21 microsatellite loci and 6 Y-chromosome genes, and the mitochondrial D-loop gene is selected to represent maternal lineages. The above genes are all from the nucleus herd of Min pigs.

RESULTS

The results of genetic structure identification and analysis of potential exogenous gene introgression patterns indicate that the nucleus herd of Min pigs maintains a high level of genetic diversity (polymorphism information content = 0.713, expected heterozygosity = 0.662, observed heterozygosity = 0.612). Compared with other Asian pig breeds, the formation of Min pig breeds is more special. Gene introgression from European pig breeds to Min pigs has occurred, which is characterized by complete introgression of paternal genes and incomplete introgression of maternal genes.

CONCLUSION

Gene introgression caused by cross-breeding is not the main factor leading to the formation of the current genetic structure of Min pigs, but this process has increased the level of genetic diversity in the nucleus herd. Compared with the influence of gene introgression, our research suggest that artificial selection and environmental adaptive evolution make Min pigs form unique genetic characteristics.

Reconstructing the ancestral gene pool to uncover the origins and genetic links of Hmong-Mien speakers

BACKGROUND

Hmong-Mien (HM) speakers are linguistically related and live primarily in China, but little is known about their ancestral origins or the evolutionary mechanism shaping their genomic diversity. In particular, the lack of whole-genome sequencing data on the Yao population has prevented a full investigation of the origins and evolutionary history of HM speakers. As such, their origins are debatable.

RESULTS

Here, we made a deep sequencing effort of 80 Yao genomes, and our analysis together with 28 East Asian populations and 968 ancient Asian genomes suggested that there is a strong genetic basis for the formation of the HM language family. We estimated that the most recent common ancestor dates to 5800 years ago, while the genetic divergence between the HM and Tai-Kadai speakers was estimated to be 8200 years ago. We proposed that HM speakers originated from the Yangtze River Basin and spread with agricultural civilization. We identified highly differentiated variants between HM and Han Chinese, in particular, a deafness-related missense variant (rs72474224) in the GJB2 gene is in a higher frequency in HM speakers than in others.

CONCLUSIONS

Our results indicated complex gene flow and medically relevant variants involved in the HM speakers' evolution history.

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.

Genetic history of the Koryaks and Evens of the Magadan region based on Y chromosome polymorphism data

In order to clarify the history of gene pool formation of the indigenous populations of the Northern Priokhotye (the northern coast of the Sea of Okhotsk), Y-chromosome polymorphisms were studied in the Koryaks and Evens living in the Magadan region. The results of the study showed that the male gene pool of the Koryaks is represented by haplogroups C-B90-B91, N-B202, and Q-B143, which are also widespread in other peoples of Northeastern Siberia, mainly of Paleo-Asiatic origin. High frequency of haplogroup C-B80, typical of other Tungus-Manchurian peoples, is characteristic of the Evens of the Magadan region. The shared components of the gene pools of the Koryaks and Evens are haplogroups R-M17 and I-P37.2 inherited as a result of admixture with Eastern Europeans (mainly Russians). The high frequency of such Y chromosome haplogroups in the Koryaks (16.7 %) and Evens (37.8 %) is indicative of close interethnic contacts during the last centuries, and most probably especially during the Soviet period. The genetic contribution of the European males' Y chromosome significantly prevails over that of maternally inherited mitochondrial DNA. The study of the Y chromosome haplogroup diversity has shown that only relatively young phylogenetic branches have been preserved in the Koryak gene pool. The age of the oldest component of the Koryak gene pool (haplogroup C-B90-B91) is estimated to be about 3.8 thousand years, the age of the younger haplogroups Q-B143 and N-B202 is about 2.8 and 2.4 thousand years, respectively. Haplogroups C-B90-B91 and N-B202 are Siberian in origin, and haplogroup Q-B143 was apparently inherited by the ancestors of the Koryaks and other Paleo-Asiatic peoples from the Paleo-Eskimos as a result of their migrations to Northeast Asia from the Americas. The analysis of microsatellite loci for haplogroup Q-B143 in the Eskimos of Greenland, Canada and Alaska as well as in the indigenous peoples of Northeastern Siberia showed a decrease in genetic diversity from east to west, pointing to the direction of distribution of the Paleo-Eskimo genetic component in the circumpolar region of America and Asia. At the same time, the Evens appeared in the Northern Priokhotye much later (in the XVII century) as a result of the expansion of the Tungusic tribes, which is confirmed by the results of the analysis of haplogroup C-B80 polymorphisms.

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.

Delineating the dispersal of Y-chromosome sub-haplogroup O2a2b-P164 among Austronesian-speaking populations

This article reports on an exploration of the Y-chromosome sub-haplogroup O2a2b-P164 in Austronesian-speaking populations. Moderate to high abundance of the P 164 mutation is seen in the West Pacific including the Amis of Formosa (36%) and the Filipinos of Mindanao (50%) as well as in the Kiritimati of Micronesia (70%), and Tonga and Samoa of West Polynesia (54% and 33%, respectively), and it drops to low frequencies in populations of East Polynesia. The communities of Polynesia and Micronesia exhibit considerable inter- and intra-population haplotype sharing suggesting extensive population affinity. The observed affinities, as well as the ages and diversity values within the P 164 sub-haplogroup among Austronesian-speaking populations signal an ancestral migration route and relationships that link the Amis of Taiwan with distant communities in West and East Polynesia, Micronesia, and the Maori of New Zealand. High resolution sequencing of the Austronesian Y chromosome indicate that the P 164 lineage originated about 19,000 ya and then split into three branches separating the Ami aborigines, Southeast Asian and Polynesian/Micronesian populations about 4700 ya, roughly coinciding with the initiation of the Austronesian diaspora. The Y-chromosomes of all the Polynesian and Micronesian population examined belong to the new FT 257096 haplogroup.

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 polymorphism of Y-chromosome in Kazakh populations from Southern Kazakhstan

BACKGROUND

The Kazakhs are one of the biggest Turkic-speaking ethnic groups, controlling vast swaths of land from the Altai to the Caspian Sea. In terms of area, Kazakhstan is ranked ninth in the world. Northern, Eastern, and Western Kazakhstan have already been studied in relation to genetic polymorphism 27 Y-STR. However, current information on the genetic polymorphism of the Y-chromosome of Southern Kazakhstan is limited only by 17 Y-STR and no geographical study of other regions has been studied at this variation.

RESULTS

The Kazakhstan Y-chromosome Haplotype Reference Database was expanded with 468 Kazakh males from the Zhambyl and Turkestan regions of South Kazakhstan by having their 27 Y-STR loci and 23 Y-SNP markers analyzed. Discrimination capacity (DC = 91.23%), haplotype match probability (HPM = 0.0029) and haplotype diversity (HD = 0.9992) are defined. Most of this Y-chromosome variability is attributed to haplogroups C2a1a1b1-F1756 (2.1%), C2a1a2-M48 (7.3%), C2a1a3-F1918 (33.3%) and C2b1a1a1a-M407 (6%). Median-joining network analysis was applied to understand the relationship between the haplotypes of the three regions. In three genetic layer can be described the position of the populations of the Southern region of Kazakhstan-the geographic Kazakh populations of Kazakhstan, the Kazakh tribal groups, and the people of bordering Asia.

CONCLUSION

The Kazakhstan Y-chromosome Haplotype Reference Database was formed for 27 Y-STR loci with a total sample of 1796 samples of Kazakhs from 16 regions of Kazakhstan. The variability of the Y-chromosome of the Kazakhs in a geographical context can be divided into four main clusters-south, north, east, west. At the same time, in the genetic space of tribal groups, the population of southern Kazakhs clusters with tribes from the same region, and genetic proximity is determined with the populations of the Hazaras of Afghanistan and the Mongols of China.

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.

Fine-scale sampling uncovers the complexity of migrations in 5th-6th century Pannonia

As the collapse of the Western Roman Empire accelerated during the 4th and 5th centuries, arriving "barbarian" groups began to establish new communities in the border provinces of the declining (and eventually former) empire. This was a time of significant cultural and political change throughout not only these border regions but Europe as a whole.1,2 To better understand post-Roman community formation in one of these key frontier zones after the collapse of the Hunnic movement, we generated new paleogenomic data for a set of 38 burials from a time series of three 5th century cemeteries3,4,5 at Lake Balaton, Hungary. We utilized a comprehensive sampling approach to characterize these cemeteries along with data from 38 additional burials from a previously published mid-6th century site6 and analyzed them alongside data from over 550 penecontemporaneous individuals.7,8,9,10,11,12,13,14,15,16,17,18,19 The range of genetic diversity in all four of these local burial communities is extensive and wider ranging than penecontemporaneous Europeans sequenced to date. Despite many commonalities in burial customs and demography, we find that there were substantial differences in genetic ancestry between the sites. We detect evidence of northern European gene flow into the Lake Balaton region. Additionally, we observe a statistically significant association between dress artifacts and genetic ancestry among 5th century genetically female burials. Our analysis shows that the formation of early Medieval communities was a multifarious process even at a local level, consisting of genetically heterogeneous groups.

SNPtotree-Resolving the Phylogeny of SNPs on Non-Recombining DNA

Genetic variants on non-recombining DNA and the hierarchical order in which they accumulate are commonly of interest. This variant hierarchy can be established and combined with information on the population and geographic origin of the individuals carrying the variants to find population structures and infer migration patterns. Further, individuals can be assigned to the characterized populations, which is relevant in forensic genetics, genetic genealogy, and epidemiologic studies. However, there is currently no straightforward method to obtain such a variant hierarchy. Here, we introduce the software SNPtotree v1.0, which uniquely determines the hierarchical order of variants on non-recombining DNA without error-prone manual sorting. The algorithm uses pairwise variant comparisons to infer their relationships and integrates the combined information into a phylogenetic tree. Variants that have contradictory pairwise relationships or ambiguous positions in the tree are removed by the software. When benchmarked using two human Y-chromosomal massively parallel sequencing datasets, SNPtotree outperforms traditional methods in the accuracy of phylogenetic trees for sequencing data with high amounts of missing information. The phylogenetic trees of variants created using SNPtotree can be used to establish and maintain publicly available phylogeny databases to further explore genetic epidemiology and genealogy, as well as population and forensic genetics.

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.

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes

Y chromosome markers can shed light on male-specific population dynamics but for many species no such markers have been discovered and are available yet, despite the potential for recovering Y-linked loci from available genome sequences. Here, we investigated how effective available bioinformatic tools are in recovering informative Y chromosome microsatellites from whole genome sequence data. In order to do so, we initially explored a large dataset of whole genome sequences comprising individuals at various coverages belonging to different species of baboons (genus: Papio) using Y chromosome references belonging to the same genus and more distantly related species (Macaca mulatta). We then further tested this approach by recovering Y-STRs from available Theropithecus gelada genomes using Papio and Macaca Y chromosome as reference sequences. Identified loci were validated in silico by a) comparing within-species relationships of Y chromosome lineages and b) genotyping male individuals in available pedigrees. Each STR was selected not to extend in its variable region beyond 100 base pairs, so that loci can be developed for PCR-based genotyping of non-invasive DNA samples. In addition to assembling a first set of Papio and Theropithecus Y-specific microsatellite markers, we released TYpeSTeR, an easy-to-use script to identify and genotype Y chromosome STRs using population genomic data which can be modulated according to available male reference genomes and genomic data, making it widely applicable across taxa.

Unveiling 2,000 years of differentiation among Tungusic-speaking populations: a revised phylogeny of the paternal founder lineage C2a-M48-SK1061

Previous studies demonstrated Y chromosome haplogroup C2a-M48-SK1061 is the only founding paternal lineage of all Tungusic-speaking populations. To infer the differentiation history of these populations, we studied more sequences and constructed downstream structure of haplogroup C2a-M48-SK1061 with better resolution. In this study, we generated 100 new sequences and co-analyzed 140 sequences of C2a-M48-SK1061 to reconstruct a highly revised phylogenetic tree with age estimates. We also performed the analysis of the geographical distribution and spatial autocorrelation of sub-branches. Dozens of new sub-branches were discovered, many sub-branches were nearly unique for Ewenki, Evens, Oroqen, Xibe, Manchu, Daur, and Mongolian. The topology of these unique sub-branches is the key evidence for understanding the complex evolutionary relationship between different Tungusic-speaking populations. The revised phylogeny provided a clear pattern for the differentiation history of haplogroup C2a-M48-SK1061 in the past 2,000 years. This study showed that the divergence pattern of founder lineage is essential to understanding the differentiation history of populations.

Mutations, sex, and genetic diversity: New arguments for partially directed evolution

A review of the theories of the existence of sexes, genetic diversity, and the distribution of mutations among organisms shows that all these concepts are not a product of random evolution and cannot be explained within the framework of Darwinism. Most mutations are the result of the genome acting on itself. This is an organized process that is implemented very differently in different species, in different places in the genome. Because of the fact that it is not random, this process must be directed and regulated, albeit with complex and not fully understood laws. This means that an additional reason must be included in order to model such mutations during evolution. The assumption of directionality must not only be explicitly included in evolutionary theory but must also occupy a central place in it. In this study an updated model of partially directed evolution is constructed, which is capable of qualitatively explaining the indicated features of evolution. Experiments are described that can confirm or disprove the proposed model.

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.

Human cooperation and evolutionary transitions in individuality

A major evolutionary transition in individuality involves the formation of a cooperative group and the transformation of that group into an evolutionary entity. Human cooperation shares principles with those of multicellular organisms that have undergone transitions in individuality: division of labour, communication, and fitness interdependence. After the split from the last common ancestor of hominoids, early hominins adapted to an increasingly terrestrial niche for several million years. We posit that new challenges in this niche set in motion a positive feedback loop in selection pressure for cooperation that ratcheted coevolutionary changes in sociality, communication, brains, cognition, kin relations and technology, eventually resulting in egalitarian societies with suppressed competition and rapid cumulative culture. The increasing pace of information innovation and transmission became a key aspect of the evolutionary niche that enabled humans to become formidable cooperators with explosive population growth, the ability to cooperate and compete in groups of millions, and emergent social norms, e.g. private property. Despite considerable fitness interdependence, the rise of private property, in concert with population explosion and socioeconomic inequality, subverts potential transition of human groups into evolutionary entities due to resurgence of latent competition and conflict. This article is part of the theme issue 'Human socio-cultural evolution in light of evolutionary transitions'.

The impacts of bronze age in the gene pool of Chinese: Insights from phylogeographics of Y-chromosomal haplogroup N1a2a-F1101

Objectives: Previous studies of archaeology and history suggested that the rise and prosperity of Bronze Age culture in East Asia had made essential contribution to the formation of early state and civilization in this region. However, the impacts in perspective of genetics remain ambiguous. Previous genetic researches indicated the Y-chromosome Q1a1a-M120 and N1a2a-F1101 may be the two most important paternal lineages among the Bronze Age people in ancient northwest China. Here, we investigated the 9,000-years history of haplogroup N1a2a-F1101 with revised phylogenetic tree and spatial autocorrelation analysis.

Materials and Methods: In this study, 229 sequences of N1a2a-F1101 were analyzed. We developed a highly-revised phylogenetic tree with age estimates for N1a2a-F1101. In addition, we also explored the geographical distribution of sub-lineages of N1a2a-F1101, and spatial autocorrelation analysis was conducted for each sub-branch.

Results: The initial differentiation location of N1a2a-F1101 and its most closely related branch, N1a2b-P43, a major lineage of Uralic-speaking populations in northern Eurasia, is likely the west part of northeast China. After ~4 thousand years of bottleneck effect period, haplgroup N1a2a-F1101 experienced continuous expansion during the Chalcolithic age (~ 4.5 kya to 4 kya) and Bronze age (~ 4 kya to 2.5 kya) in northern China. Ancient DNA evidence supported that this haplogroup is the lineage of ruling family of Zhou Dynasty (~ 3 kya-2.2 kya) of ancient China.

Discussion: In general, we proposed that the Bronze Age people in the border area between the eastern Eurasian steppe and northern China not only played a key role in promoting the early state and civilization of China, but also left significant traces in the gene pool of Chinese people.

Palaeogenomics of Upper Palaeolithic to Neolithic European hunter-gatherers

Modern humans have populated Europe for more than 45,000 years1,2. Our knowledge of the genetic relatedness and structure of ancient hunter-gatherers is however limited, owing to the scarceness and poor molecular preservation of human remains from that period3. Here we analyse 356 ancient hunter-gatherer genomes, including new genomic data for 116 individuals from 14 countries in western and central Eurasia, spanning between 35,000 and 5,000 years ago. We identify a genetic ancestry profile in individuals associated with Upper Palaeolithic Gravettian assemblages from western Europe that is distinct from contemporaneous groups related to this archaeological culture in central and southern Europe4, but resembles that of preceding individuals associated with the Aurignacian culture. This ancestry profile survived during the Last Glacial Maximum (25,000 to 19,000 years ago) in human populations from southwestern Europe associated with the Solutrean culture, and with the following Magdalenian culture that re-expanded northeastward after the Last Glacial Maximum. Conversely, we reveal a genetic turnover in southern Europe suggesting a local replacement of human groups around the time of the Last Glacial Maximum, accompanied by a north-to-south dispersal of populations associated with the Epigravettian culture. From at least 14,000 years ago, an ancestry related to this culture spread from the south across the rest of Europe, largely replacing the Magdalenian-associated gene pool. After a period of limited admixture that spanned the beginning of the Mesolithic, we find genetic interactions between western and eastern European hunter-gatherers, who were also characterized by marked differences in phenotypically relevant variants.

An evolutionary case for polygyny to counter demographic collapse

Sex ratio theory suggests why mating practices have become dysfunctional in the West and other regions. Spain, Japan, and over 20 other nations are on course to have their populations halved by 2100, dramatically aging their citizenry. Experts and opinion makers warn that a demographic collapse cannot be absorbed by our current social order; Elon Musk proclaims this to be "the biggest threat to human civilization." Statistics from the Nordic countries-the world's most gender-equal region-indicate that subjective perceptions of the sex ratio in modern environments drive singledom and low reproduction. Scandinavia has the world's highest occurrence of one-person households: 43-46%. In the past decade, the Norwegian fertility rate dropped from 2.0 to 1.5. Sex ratio studies suggest that women's perception of there being few acceptable partners activates a polygynous mindset, which in prosperous, monogamous societies drives promiscuity to the detriment of pair-bonding. More than 6 million years of hominin evolution under promiscuous, polygynous, and monogamous regimes shaped mate preferences that evoke different cultural and behavioral responses as environments change. The Church's imposition of lifelong monogamy contributed to the emergence of the modern world, but if this world's gender-equal societies no longer motivate reproduction, being more open to polygyny could be worth considering as a means for increasing fertility. This article makes this case by exploring hominin mating from our last common ancestor with chimpanzees-through the genus Homo's forager and agricultural periods-to modern Scandinavians. In the past millennium, mating practices have coevolved with the emergence of modernity, necessitating frequent cultural updates. An evolutionary analysis of Nordic works of literature illuminates the ways in which ideological narratives influence reproductive norms. The insights gleaned are considered in the context of people's perceived sex ratio.

Bidirectional dispersals during the peopling of the North American Arctic

It is thought that Paleo-Inuit were the first people that settled the American Arctic about 5000 BP (before the present) from a migration that crossed Beringia from Northeast Asia. It is theorized that this group initially migrated to the North Slopes of Alaska and subsequently expanded eastward, eventually reaching Greenland. A second circumpolar dispersal of Neo-Inuit from the North Slopes associated with the Thule-Inuk culture has been postulated to have extended eastward around 800 BP, totally replacing the original Paleo-Inuit without admixing. Although generally accepted, this migration scenario is incompatible with previously reported indications of east to west gene flow across the American Arctic. Here we report on the Y-chromosome haplogroup and Y-STR diversity of the four circumpolar populations of the Tuva Republic (N = 24), Northeast Siberia (N = 9), Bethel, Alaska (N = 40), and Barrow, Alaska (N = 31). Four haplogroup lineages (Q-NWT01, Q-M3, Q-M346, and Q-M120) were detected, Q-NWT01 and Q-M3 being the most abundant at 11.11 and 66.67% in Northeast Siberia, 32.50 and 65.00% in Bethel, and 67.74 and 32.26% in Barrow, respectively. The same samples genotyped for Y-chromosome SNPs were typed for 17 Y-STYR loci using the AmpFlSTR Yfiler system. Age estimates and diversity values for the Q-NWT01 and Q-M3 mutations suggest extensive movement of male individuals along the entire longitudinal stretch of the American circumpolar region. Throughout the entire region, Q-M3 exhibits a west to east decreasing gradient in age and diversity while Q-NWT01 indicates the opposite with older TMRCA and higher diversity values running from east to west with the most recent estimates in Canada and Alaska. The high age and diversity values in Greenland are congruent with an origin of the Q-NWT01 mutation in the east of the circumpolar range about 2000-3000 ya. This scenario is incompatible with a complete biological replacement starting about 700 BP of Paleo-Inuit like the Dorset by the Thule-Inuit (Neo-Inuit), as is currently thought, and more parsimonious with gene flow carrying the NWT01 mutation from a pre-Thule population to the ancestors of the present-day Inuit.

Understanding the sex inequality in childlessness: an approach using Swedish register data

In most countries, men are more likely to be childless than women. Understanding how this inequality arises is important given the significance of parenthood for individuals' lives. The objective of this study was to explore how three prominent explanations for sex inequalities in childlessness relate to the Sex Gap in Childlessness (SGC) in Sweden. The three explanations examined were sex ratio imbalance (more men than women), mismeasurement of fatherhood (inequalities in registration) and partnership differences (inequality in multi-partner fertility). Administrative register data for cohorts born in 1945-1974 were used. The population was restricted to men and women who were born in Sweden or arrived prior to the age of 15, and all registered childbearing partnerships were examined. To explore the possible significance of the three explanations, counter-factual standardization was used. Of the three explanations examined, the population sex ratio had the largest positive impact on the SGC, while multi-partner fertility had a negative impact. The results show that inequalities in the sex ratio can explain about 20-34% of the SGC depending on cohort. Inequalities in registration of fathers explain about 9-24% of the SGC depending on cohort. Finally, results show that women are slightly more likely to have multiple partners, and that this behaviour has a substantial minimizing effect on the SGC (minimizing it by 6-65%). To the authors' knowledge this was the first paper to estimate the scope of the impacts of these three mechanisms on the SGC. Differences in multi-partner fertility have in many instances been used as an explanation for men's higher childlessness. This study shows that women have slightly more childbearing partners than men, and that this actually leads to a smaller SGC in the studied population.

Analysis of Y chromosome haplogroups in Parkinson's disease

Parkinson's disease is a complex neurodegenerative disorder that is about 1.5 times more prevalent in males than females. Extensive work has been done to identify the genetic risk factors behind Parkinson's disease on autosomes and more recently on Chromosome X, but work remains to be done on the male-specific Y chromosome. In an effort to explore the role of the Y chromosome in Parkinson's disease, we analysed whole-genome sequencing data from the Accelerating Medicines Partnership-Parkinson's disease initiative (1466 cases and 1664 controls), genotype data from NeuroX (3491 cases and 3232 controls) and genotype data from UKBiobank (182 517 controls, 1892 cases and 3783 proxy cases), all consisting of male European ancestry samples. We classified sample Y chromosomes by haplogroup using three different tools for comparison (Snappy, Yhaplo and Y-LineageTracker) and meta-analysed this data to identify haplogroups associated with Parkinson's disease. This was followed up with a Y-chromosome association study to identify specific variants associated with disease. We also analysed blood-based RNASeq data obtained from the Accelerating Medicines Partnership-Parkinson's disease initiative (1020 samples) and RNASeq data obtained from the North American Brain Expression Consortium (171 samples) to identify Y-chromosome genes differentially expressed in cases, controls, specific haplogroups and specific tissues. RNASeq analyses suggest Y-chromosome gene expression differs between brain and blood tissues but does not differ significantly in cases, controls or specific haplogroups. Overall, we did not find any strong associations between Y-chromosome genetics and Parkinson's disease, suggesting the explanation for the increased prevalence in males may lie elsewhere.

Ancient Components and Recent Expansion in the Eurasian Heartland: Insights into the Revised Phylogeny of Y-Chromosomes from Central Asia

In the past two decades, studies of Y chromosomal single nucleotide polymorphisms (Y-SNPs) and short tandem repeats (Y-STRs) have shed light on the demographic history of Central Asia, the heartland of Eurasia. However, complex patterns of migration and admixture have complicated population genetic studies in Central Asia. Here, we sequenced and analyzed the Y-chromosomes of 187 male individuals from Kazakh, Kyrgyz, Uzbek, Karakalpak, Hazara, Karluk, Tajik, Uyghur, Dungan, and Turkmen populations. High diversity and admixture from peripheral areas of Eurasia were observed among the paternal gene pool of these populations. This general pattern can be largely attributed to the activities of ancient people in four periods, including the Neolithic farmers, Indo-Europeans, Turks, and Mongols. Most importantly, we detected the consistent expansion of many minor lineages over the past thousand years, which may correspond directly to the formation of modern populations in these regions. The newly discovered sub-lineages and variants provide a basis for further studies of the contributions of minor lineages to the formation of modern populations in Central Asia.

Genetic insights into the social organization of Neanderthals

Genomic analyses of Neanderthals have previously provided insights into their population history and relationship to modern humans1-8, but the social organization of Neanderthal communities remains poorly understood. Here we present genetic data for 13 Neanderthals from two Middle Palaeolithic sites in the Altai Mountains of southern Siberia: 11 from Chagyrskaya Cave9,10 and 2 from Okladnikov Cave11-making this one of the largest genetic studies of a Neanderthal population to date. We used hybridization capture to obtain genome-wide nuclear data, as well as mitochondrial and Y-chromosome sequences. Some Chagyrskaya individuals were closely related, including a father-daughter pair and a pair of second-degree relatives, indicating that at least some of the individuals lived at the same time. Up to one-third of these individuals' genomes had long segments of homozygosity, suggesting that the Chagyrskaya Neanderthals were part of a small community. In addition, the Y-chromosome diversity is an order of magnitude lower than the mitochondrial diversity, a pattern that we found is best explained by female migration between communities. Thus, the genetic data presented here provide a detailed documentation of the social organization of an isolated Neanderthal community at the easternmost extent of their known range.

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.

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.

Sequence analyses of Malaysian Indigenous communities reveal historical admixture between Hoabinhian hunter-gatherers and Neolithic farmers

Southeast Asia comprises 11 countries that span mainland Asia across to numerous islands that stretch from the Andaman Sea to the South China Sea and Indian Ocean. This region harbors an impressive diversity of history, culture, religion and biology. Indigenous people of Malaysia display substantial phenotypic, linguistic, and anthropological diversity. Despite this remarkable diversity which has been documented for centuries, the genetic history and structure of indigenous Malaysians remain under-studied. To have a better understanding about the genetic history of these people, especially Malaysian Negritos, we sequenced whole genomes of 15 individuals belonging to five indigenous groups from Peninsular Malaysia and one from North Borneo to high coverage (30X). Our results demonstrate that indigenous populations of Malaysia are genetically close to East Asian populations. We show that present-day Malaysian Negritos can be modeled as an admixture of ancient Hoabinhian hunter-gatherers and Neolithic farmers. We observe gene flow from South Asian populations into the Malaysian indigenous groups, but not into Dusun of North Borneo. Our study proposes that Malaysian indigenous people originated from at least three distinct ancestral populations related to the Hoabinhian hunter-gatherers, Neolithic farmers and Austronesian speakers.

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.

Genetic characterization of populations in the Marquesas Archipelago in the context of the Austronesian expansion

Our exploration of the genetic constitution of Nuku Hiva (n = 51), Hiva Oa (n = 28) and Tahuata (n = 8) of the Marquesas Archipelago based on the analyses of genome-wide autosomal markers as well as high-resolution genotyping of paternal and maternal lineages provides us with information on the origins and settlement of these islands at the fringe of the Austronesian expansion. One widespread theme that emerges from this study is the genetic uniformity and relative isolation exhibited by the Marquesas and Society populations. This genetic homogeneity within East Polynesia groups is reflected in their limited average heterozygosity, uniformity of constituents in the Structure analyses, reiteration of complete mtDNA sequences, marked separation from Asian and other Oceanic populations in the PC analyses, limited differentiation in the PCAs and large number of IBD segments in common. Both the f3 and the Outgroup f3 results provide indications of intra-East Polynesian gene flow that may have promoted the observed intra-East Polynesia genetic homogeneity while ALDER analyses indicate that East Polynesia experienced two gene flow episodes, one relatively recent from Europe that coincides roughly with the European incursion into the region and an early one that may represent the original settlement of the islands by Austronesians. Median Network analysis based on high-resolution Y-STR loci under C2a-M208 generates a star-like topology with East Polynesian groups (especially from the Society Archipelago) in central stem positions and individuals from the different populations radiating out one mutational step away while several Samoan and outlier individuals occupy peripheral positions. This arrangement of populations is congruent with dispersals of C2a-M208 Y chromosomes from East Polynesia as a migration hub signaling dispersals in various directions. The equivalent ages of the C2a-M208 lineage of the populations in the Network corroborate an east to west flow of the most abundant Polynesian Y chromosome.

Episodes of diversification and isolation in Island Southeast Asian and Near Oceanian male lineages

Island Southeast Asia (ISEA) and Oceania host one of the world’s richest assemblages of human phenotypic, linguistic, and cultural diversity. Despite this, the region’s male genetic lineages are globally among the last to remain unresolved. We compiled ∼9.7 Mb of Y chromosome (chrY) sequence from a diverse sample of over 380 men from this region, including 152 first reported here. The granularity of this data set allows us to fully resolve and date the regional chrY phylogeny. This new high-resolution tree confirms two main population bursts: multiple rapid diversifications following the region’s initial settlement ∼50 kya, and extensive expansions <6 kya. Notably, ∼40–25 kya the deep rooting local lineages of C-M130, M-P256, and S-B254 show almost no further branching events in ISEA, New Guinea, and Australia, matching a similar pause in diversification seen in maternal mitochondrial DNA lineages. The main local lineages start diversifying ∼25 kya, at the time of the last glacial maximum. This improved chrY topology highlights localized events with important historical implications, including pre-Holocene contact between Mainland and ISEA, potential interactions between Australia and the Papuan world, and a sustained period of diversification following the flooding of the ancient Sunda and Sahul continents as the insular landscape observed today formed. The high-resolution phylogeny of the chrY presented here thus enables a detailed exploration of past isolation, interaction, and change in one of the world’s least understood regions.

An invasive Haemophilus influenzae serotype b infection in an Anglo-Saxon plague victim

Background

The human pathogen Haemophilus influenzae was the main cause of bacterial meningitis in children and a major cause of worldwide infant mortality before the introduction of a vaccine in the 1980s. Although the occurrence of serotype b (Hib), the most virulent type of H. influenzae, has since decreased, reports of infections with other serotypes and non-typeable strains are on the rise. While non-typeable strains have been studied in-depth, very little is known of the pathogen’s evolutionary history, and no genomes dating prior to 1940 were available.

Results

We describe a Hib genome isolated from a 6-year-old Anglo-Saxon plague victim, from approximately 540 to 550 CE, Edix Hill, England, showing signs of invasive infection on its skeleton. We find that the genome clusters in phylogenetic division II with Hib strain NCTC8468, which also caused invasive disease. While the virulence profile of our genome was distinct, its genomic similarity to NCTC8468 points to mostly clonal evolution of the clade since the 6th century. We also reconstruct a partial Yersinia pestis genome, which is likely identical to a published first plague pandemic genome of Edix Hill.

Conclusions

Our study presents the earliest genomic evidence for H. influenzae, points to the potential presence of larger genomic diversity in the phylogenetic division II serotype b clade in the past, and allows the first insights into the evolutionary history of this major human pathogen. The identification of both plague and Hib opens questions on the effect of plague in immunocompromised individuals already affected by infectious diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02580-z.

Placing Ancient DNA Sequences into Reference Phylogenies

Joint phylogenetic analysis of ancient DNA (aDNA) with modern phylogenies is hampered by low sequence coverage and post-mortem deamination, often resulting in overconservative or incorrect assignment. We provide a new efficient likelihood-based workflow, pathPhynder, that takes advantage of all the polymorphic sites in the target sequence. This effectively evaluates the number of ancestral and derived alleles present on each branch and reports the most likely placement of an ancient sample in the phylogeny and a haplogroup assignment, together with alternatives and supporting evidence. To illustrate the application of pathPhynder, we show improved Y chromosome assignments for published aDNA sequences, using a newly compiled Y variation data set (120,908 markers from 2,014 samples) that significantly enhances Y haplogroup assignment for low coverage samples. We apply the method to all published male aDNA samples from Africa, giving new insights into ancient migrations and the relationships between ancient and modern populations. The same software can be used to place samples with large amounts of missing data into other large non-recombining phylogenies such as the mitochondrial tree.

Overview of the Americas’ First Peopling from a Patrilineal Perspective: New Evidence from the Southern Continent

Uniparental genetic systems are unique sex indicators and complement the study of autosomal diversity by providing landmarks of human migrations that repeatedly shaped the structure of extant populations. Our knowledge of the variation of the male-specific region of the Y chromosome in Native Americans is still rather scarce and scattered, but by merging sequence information from modern and ancient individuals, we here provide a comprehensive and updated phylogeny of the distinctive Native American branches of haplogroups C and Q. Our analyses confirm C-MPB373, C-P39, Q-Z780, Q-M848, and Q-Y4276 as the main founding haplogroups and identify traces of unsuccessful (pre-Q-F1096) or extinct (C-L1373*, Q-YP4010*) Y-chromosome lineages, indicating that haplogroup diversity of the founder populations that first entered the Americas was greater than that observed in the Indigenous component of modern populations. In addition, through a diachronic and phylogeographic dissection of newly identified Q-M848 branches, we provide the first Y-chromosome insights into the early peopling of the South American hinterland (Q-BY104773 and Q-BY15730) and on overlying inland migrations (Q-BY139813).

The genetic origin of Daunians and the Pan-Mediterranean southern Italian Iron Age context

The geographical location and shape of Apulia, a narrow land stretching out in the sea at the South of Italy, made this region a Mediterranean crossroads connecting Western Europe and the Balkans. Such movements culminated at the beginning of the Iron Age with the Iapygian civilization which consisted of three cultures: Peucetians, Messapians, and Daunians. Among them, the Daunians left a peculiar cultural heritage, with one-of-a-kind stelae and pottery, but, despite the extensive archaeological literature, their origin has been lost to time. In order to shed light on this and to provide a genetic picture of Iron Age Southern Italy, we collected and sequenced human remains from three archaeological sites geographically located in Northern Apulia (the area historically inhabited by Daunians) and radiocarbon dated between 1157 and 275 calBCE. We find that Iron Age Apulian samples are still distant from the genetic variability of modern-day Apulians, they show a degree of genetic heterogeneity comparable with the cosmopolitan Republican and Imperial Roman civilization, even though a few kilometers and centuries separate them, and they are well inserted into the Iron Age Pan-Mediterranean genetic landscape. Our study provides for the first time a window on the genetic make-up of pre-Roman Apulia, whose increasing connectivity within the Mediterranean landscape, would have contributed to laying the foundation for modern genetic variability. In this light, the genetic profile of Daunians may be compatible with an at least partial autochthonous origin, with plausible contributions from the Balkan peninsula.

Perceived Effectiveness of Flirtation Tactics: The Effects of sex, Mating Context and Individual Differences in US and Norwegian Samples

Flirting involves various signals communicated between individuals. To attract potential mates, men and women exhibit flirtatious behavior to get the attention of, and potentially elicit sexual or romantic interest from, a desired partner. In this first large, preregistered study of judgement of the effectiveness of flirtation tactics based on Sexual Strategies Theory, we considered the effects of flirter's (actor) sex and mating contexts in addition to rater's (participant) sex across two cultures, Norway and the U.S. Culturally relevant covariates such as sociosexuality, extraversion, mate value, age, and religiosity were examined. Participants from Norway (N = 415, 56% women) and the US (N = 577, 69% women) responded to one of four different randomized questionnaires representing a factorial design considering either short-term versus long-term mating context and either female or male sex of actor. We found that sexual availability cues were judged more effective when employed by women in short-term mating contexts. Friendly contact, such as hugs or kissing on the cheek, was not. Cues to generosity and commitment were judged more effective when employed by men in long-term mating contexts. Humor was rated as more effective when used by men and in long-term contexts, and least effective when used by women in short term contexts. However, laughing or giggling at someone's jokes was an effective flirtation tactic for both sexes. Overall, predictions for culturally relevant covariates were not supported, but cultural differences were found in bodily displays, initial contact, and generosity. These findings dovetail neatly with findings from the self-promotion literature, and further support that flirtation is a universal mate signaling strategy.

Is History the Same as Evolution? No. Is it Independent of Evolution? Certainly Not

History is full of violence and oppression within and between groups, and although group conflicts enhance within-group cooperation (mediated by oxytocin, which promotes parochial altruism) the hierarchy within groups ensures that spoils accrue very unevenly. Darwin suggested, and we now know, that sexual selection is as powerful as selection by mortality, and the main purpose of survival is reproduction. Male reproductive skew is greater than that among females in all societies, but the difference became much greater after the hunting-gathering era, and the rise of so-called "civilization" was everywhere a process of predatory expansion, producing kingdoms and empires where top males achieved astounding heights of reproductive success. This was shown by historical and ethnographic data now strongly confirmed by genomic science. Psychological research confirms that group identity, out-group stigmatization, leadership characterized by charisma, the will to power, narcissism, sociopathy, and cruelty, and followership characterized by hypnotic obedience, loss of individuality, and cruelty are integral parts of human nature. We can thank at least ten or twelve millennia of microevolutionary processes such as those described above, all more prominent in males than females. Followers in wars have faced a difficult risk-benefit analysis, but if they survived and won they too could increase their reproductive success through the rape and other sexual exploitation that have accompanied almost all wars. For modern leaders, social monogamy and contraception have separated autocracy from reproductive success, but only partly, and current worldwide autocratic trends still depend on the evolved will to power, obedience, and cruelty.

Genetic connections and convergent evolution of tropical indigenous peoples in Asia

Tropical indigenous peoples in Asia (TIA) attract much attention for their unique appearance, whereas their genetic history and adaptive evolution remain mysteries. We conducted a comprehensive study to characterize the genetic distinction and connection of broad geographical TIAs. Despite the diverse genetic makeup and large interarea genetic differentiation between the TIA groups, we identified a basal Asian ancestry (bASN) specifically shared by these populations. The bASN ancestry was relatively enriched in ancient Asian human genomes dated as early as ∼50,000 years before the present and diminished in more recent history. Notably, the bASN ancestry is unlikely to be derived from archaic hominins. Instead, we suggest it may be better modeled as a survived lineage of the initial peopling of Asia. Shared adaptations inherited from the ancient Asian ancestry were detected among the TIA groups (e.g., LIMS1 for hair morphology, and COL24A1 for bone formation), and they are enriched in neurological functions either at an identical locus (e.g., NKAIN3), or different loci in an identical gene (e.g., TENM4). The bASN ancestry could also have formed the substrate of the genetic architecture of the dark pigmentation observed in the TIA peoples. We hypothesize that phenotypic convergence of the dark pigmentation in TIAs could have resulted from parallel (e.g., DDB1/DAK) or genetic convergence driven by admixture (e.g., MTHFD1 and RAD18), new mutations (e.g., STK11), or notably purifying selection (e.g., MC1R). Our results provide new insights into the initial peopling of Asia and an advanced understanding of the phenotypic convergence of the TIA peoples.

A Genetic Variation in the Y Chromosome Among Modern Japanese Males Related to Several Physiological and Psychological Characteristics

Previous studies in population genetics have proposed that the Y-chromosomal (Y-DNA) haplogroup D ancestor likely originated from Africa. The haplogroup D branch next started Out-of-Africa migration, rapidly expanded across Eurasia, and later diversified in East Asia. Y-DNA haplogroup D-M55, one of the branches of haplogroup D, is only found in modern Japanese males, suggesting that individuals with Y-DNA haplogroup D migrated from the Eurasian continent. Based on previous observations, Y-DNA haplogroup D is expected to be associated with some male characteristics including personality. Therefore, this study investigated whether the Y-DNA haplogroup D-M55 is associated with several physiological and psychological characteristics, including exploratory motivation and human relationship-related perception. We recruited Japanese young adult males and females and investigated the association between Y-DNA haplogroup D-M55, physiological [body mass index (BMI)], and several psychological parameters [perceived number of close friends, behavioral inhibition system/behavioral activation system (BIS/BAS), perceived happiness, and perceived loneliness]. The results indicated that males with haplogroup D-M55 had a higher BMI and more close friends, compared with non-carrier males. Additional multiple regression analyses, which tested the hypothesis that haplogroup D-M55 predicts BMI and perceived number of close friends, confirmed our hypothesis, even after controlling for the potentially confounding variables of age and sex. We also analyzed the gene-gene interaction between haplogroup D-M55 and an autosomal gene polymorphism associated with BMI and human relationships, such as the dopamine D2 receptor gene (DRD2: rs1800497). Results showed gene-gene interactions between haplogroups D-M55 and DRD2 in BMI. Based on these findings, it is demonstrated that Y-DNA haplogroup D is associated with human personality.

New insights into the evolution of human Y chromosome palindromes through mutation and gene conversion

About one-quarter of the euchromatic portion of the male-specific region of the human Y chromosome consists of large duplicated sequences that are organized in eight palindromes (termed P1-P8), which undergo arm-to arm gene conversion, a proposed mechanism for maintaining their sequence integrity. Although the relevance of gene conversion in the evolution of palindromic sequences has been profoundly recognized, the dynamic of this mechanism is still nuanced. To shed light into the evolution of these genomic elements, we performed a high-depth (50×) targeted next-generation sequencing of the palindrome P6 in 157 subjects belonging to the most divergent evolutionary lineages of the Y chromosome. We found 118 new paralogous sequence variants, which were placed into the context of a robust Y chromosome phylogeny based on 7240 SNPs of the X-degenerate region. We mapped along the phylogeny 80 gene conversion events that shaped the diversity of P6 arms during recent human history. In contrast to previous studies, we demonstrated that arm-to-arm gene conversion, which occurs at a rate of 6.01 × 10 -6 conversions/base/year, is not biased toward the retention of the ancestral state of sequences. We also found a significantly lower mutation rate of the arms (6.18 × 10-10 mutations/base/year) compared with the spacer (9.16 × 10-10 mutations/base/year), a finding that may explain the observed higher inter-species conservation of arms, without invoking any bias of conversion. Finally, by formally testing the mutation/conversion balance in P6, we found that the arms of this palindrome reached a steady-state equilibrium between mutation and gene conversion.

Nine Levels of Explanation : A Proposed Expansion of Tinbergen's Four-Level Framework for Understanding the Causes of Behavior

Tinbergen's classic "On Aims and Methods of Ethology" (Zeitschrift für Tierpsychologie, 20, 1963) proposed four levels of explanation of behavior, which he thought would soon apply to humans. This paper discusses the need for multilevel explanation; Huxley and Mayr's prior models, and others that followed; Tinbergen's differences with Lorenz on "the innate"; and Mayr's ultimate/proximate distinction. It synthesizes these approaches with nine levels of explanation in three categories: phylogeny, natural selection, and genomics (ultimate causes); maturation, sensitive period effects, and routine environmental effects (intermediate causes); and hormonal/metabolic processes, neural circuitry, and eliciting stimuli (proximate causes), as a respectful extension of Tinbergen's levels. The proposed classification supports and builds on Tinbergen's multilevel model and Mayr's ultimate/proximate continuum, adding intermediate causes in accord with Tinbergen's emphasis on ontogeny. It requires no modification of Standard Evolutionary Theory or The Modern Synthesis, but shows that much that critics claim was missing was in fact part of Neo-Darwinian theory (so named by J. Mark Baldwin in The American Naturalist in 1896) all along, notably reciprocal causation in ontogeny, niche construction, cultural evolution, and multilevel selection. Updates of classical examples in ethology are offered at each of the nine levels, including the neuroethological and genomic findings Tinbergen foresaw. Finally, human examples are supplied at each level, fulfilling his hope of human applications as part of the biology of behavior. This broad ethological framework empowers us to explain human behavior-eventually completely-and vindicates the idea of human nature, and of humans as a part of nature.

Revisiting the out of Africa event with a deep-learning approach

Anatomically modern humans evolved around 300 thousand years ago in Africa. They started to appear in the fossil record outside of Africa as early as 100 thousand years ago, although other hominins existed throughout Eurasia much earlier. Recently, several studies argued in favor of a single out of Africa event for modern humans on the basis of whole-genome sequence analyses. However, the single out of Africa model is in contrast with some of the findings from fossil records, which support two out of Africa events, and uniparental data, which propose a back to Africa movement. Here, we used a deep-learning approach coupled with approximate Bayesian computation and sequential Monte Carlo to revisit these hypotheses from the whole-genome sequence perspective. Our results support the back to Africa model over other alternatives. We estimated that there are two sequential separations between Africa and out of African populations happening around 60-90 thousand years ago and separated by 13-15 thousand years. One of the populations resulting from the more recent split has replaced the older West African population to a large extent, while the other one has founded the out of Africa populations.