ABSTRACT
The extent to which prehistoric migrations of farmers influenced the genetic pool of western North Africans remains unclear. Archaeological evidence suggests that the Neolithization process may have happened through the adoption of innovations by local Epipaleolithic communities or by demic diffusion from the Eastern Mediterranean shores or Iberia. Here, we present an analysis of individuals' genome sequences from Early and Late Neolithic sites in Morocco and from Early Neolithic individuals from southern Iberia. We show that Early Neolithic Moroccans (â¼5,000 BCE) are similar to Later Stone Age individuals from the same region and possess an endemic element retained in present-day Maghrebi populations, confirming a long-term genetic continuity in the region. This scenario is consistent with Early Neolithic traditions in North Africa deriving from Epipaleolithic communities that adopted certain agricultural techniques from neighboring populations. Among Eurasian ancient populations, Early Neolithic Moroccans are distantly related to Levantine Natufian hunter-gatherers (â¼9,000 BCE) and Pre-Pottery Neolithic farmers (â¼6,500 BCE). Late Neolithic (â¼3,000 BCE) Moroccans, in contrast, share an Iberian component, supporting theories of trans-Gibraltar gene flow and indicating that Neolithization of North Africa involved both the movement of ideas and people. Lastly, the southern Iberian Early Neolithic samples share the same genetic composition as the Cardial Mediterranean Neolithic culture that reached Iberia â¼5,500 BCE. The cultural and genetic similarities between Iberian and North African Neolithic traditions further reinforce the model of an Iberian migration into the Maghreb.
Subject(s)
Ethnicity/genetics , Genome, Human , Human Migration/history , Africa, Northern , Agriculture/history , Chromosomes, Human, Y/genetics , DNA, Mitochondrial/genetics , Ethnicity/history , Europe , Gene Flow , Gene Library , Genetics, Population , History, Ancient , Humans , Middle East , Morocco , Sequence Analysis, DNA , Spain/ethnologyABSTRACT
Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation â¼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.
Subject(s)
Genome/genetics , Hylobates/classification , Hylobates/genetics , Karyotype , Phylogeny , Animals , Evolution, Molecular , Hominidae/classification , Hominidae/genetics , Humans , Molecular Sequence Data , Retroelements/genetics , Selection, Genetic , Transcription Termination, GeneticABSTRACT
Sequencing the genomes of extinct hominids has reshaped our understanding of modern human origins. Here, we analyze â¼120 kb of exome-captured Y-chromosome DNA from a Neandertal individual from El Sidrón, Spain. We investigate its divergence from orthologous chimpanzee and modern human sequences and find strong support for a model that places the Neandertal lineage as an outgroup to modern human Y chromosomes-including A00, the highly divergent basal haplogroup. We estimate that the time to the most recent common ancestor (TMRCA) of Neandertal and modern human Y chromosomes is â¼588 thousand years ago (kya) (95% confidence interval [CI]: 447-806 kya). This is â¼2.1 (95% CI: 1.7-2.9) times longer than the TMRCA of A00 and other extant modern human Y-chromosome lineages. This estimate suggests that the Y-chromosome divergence mirrors the population divergence of Neandertals and modern human ancestors, and it refutes alternative scenarios of a relatively recent or super-archaic origin of Neandertal Y chromosomes. The fact that the Neandertal Y we describe has never been observed in modern humans suggests that the lineage is most likely extinct. We identify protein-coding differences between Neandertal and modern human Y chromosomes, including potentially damaging changes to PCDH11Y, TMSB4Y, USP9Y, and KDM5D. Three of these changes are missense mutations in genes that produce male-specific minor histocompatibility (H-Y) antigens. Antigens derived from KDM5D, for example, are thought to elicit a maternal immune response during gestation. It is possible that incompatibilities at one or more of these genes played a role in the reproductive isolation of the two groups.
Subject(s)
Chromosomes, Human, Y/genetics , Neanderthals/genetics , Animals , Cadherins/genetics , DNA/genetics , Genetic Variation , Histone Demethylases/genetics , Humans , Male , Minor Histocompatibility Antigens/genetics , Mutation, Missense , Pan troglodytes/genetics , Protocadherins , Sequence Analysis, DNA , Spain , Ubiquitin Thiolesterase/geneticsABSTRACT
Difference in male and female effective population sizes has, at times, been attributed to both sexes having unequal variance in their number of offspring. Such difference is paralleled by the relative effective sizes of autosomes, sex chromosomes, and mitochondrial DNA. I develop a simple framework to calculate the inbreeding effective population sizes for loci with different modes of inheritance. In this framework, I separate the effects due to mating strategy and those due to genetic transmission. I then show that, in addition to differences in the variance in offspring number, skew in the male/female effective sizes can also be caused by family composition. This approach can be used to illustrate the effect of induced behaviors on the relative male and female effective population sizes. In particular, I show the impact of the one-child policy formerly implemented in the People's Republic of China on the relative male and female effective population sizes. Furthermore, I argue that, under some strong constraints on family structure, the concepts of male and female effective population sizes are invalid.
Subject(s)
Consanguinity , Genetics, Population , Population Density , China , Female , Humans , Male , Reproduction , Sex RatioABSTRACT
We report the discovery of an African American Y chromosome that carries the ancestral state of all SNPs that defined the basal portion of the Y chromosome phylogenetic tree. We sequenced â¼240 kb of this chromosome to identify private, derived mutations on this lineage, which we named A00. We then estimated the time to the most recent common ancestor (TMRCA) for the Y tree as 338 thousand years ago (kya) (95% confidence interval = 237-581 kya). Remarkably, this exceeds current estimates of the mtDNA TMRCA, as well as those of the age of the oldest anatomically modern human fossils. The extremely ancient age combined with the rarity of the A00 lineage, which we also find at very low frequency in central Africa, point to the importance of considering more complex models for the origin of Y chromosome diversity. These models include ancient population structure and the possibility of archaic introgression of Y chromosomes into anatomically modern humans. The A00 lineage was discovered in a large database of consumer samples of African Americans and has not been identified in traditional hunter-gatherer populations from sub-Saharan Africa. This underscores how the stochastic nature of the genealogical process can affect inference from a single locus and warrants caution during the interpretation of the geographic location of divergent branches of the Y chromosome phylogenetic tree for the elucidation of human origins.
Subject(s)
Black or African American/genetics , Chromosomes, Human, Y , DNA, Mitochondrial/genetics , Africa South of the Sahara , Africa, Central , Cell Lineage , Female , Fossils , Genetics, Population/methods , Humans , Male , Phylogeny , Polymorphism, Single NucleotideABSTRACT
Signals of archaic admixture have been identified through comparisons of the draft Neanderthal and Denisova genomes with those of living humans. Studies of individual loci contributing to these genome-wide average signals are required for characterization of the introgression process and investigation of whether archaic variants conferred an adaptive advantage to the ancestors of contemporary human populations. However, no definitive case of adaptive introgression has yet been described. Here we provide a DNA sequence analysis of the innate immune gene STAT2 and show that a haplotype carried by many Eurasians (but not sub-Saharan Africans) has a sequence that closely matches that of the Neanderthal STAT2. This haplotype, referred to as N, was discovered through a resequencing survey of the entire coding region of STAT2 in a global sample of 90 individuals. Analyses of publicly available complete genome sequence data show that haplotype N shares a recent common ancestor with the Neanderthal sequence (~80 thousand years ago) and is found throughout Eurasia at an average frequency of ~5%. Interestingly, N is found in Melanesian populations at ~10-fold higher frequency (~54%) than in Eurasian populations. A neutrality test that controls for demography rejects the hypothesis that a variant of N rose to high frequency in Melanesia by genetic drift alone. Although we are not able to pinpoint the precise target of positive selection, we identify nonsynonymous mutations in ERBB3, ESYT1, and STAT2-all of which are part of the same 250 kb introgressive haplotype-as good candidates.
Subject(s)
Genetics, Population , Haplotypes/genetics , Neanderthals/genetics , Phylogeny , STAT2 Transcription Factor/genetics , Selection, Genetic , Adaptation, Biological/genetics , Animals , Base Sequence , Evolution, Molecular , Gene Components , Genotype , Humans , Linkage Disequilibrium , Molecular Sequence Data , Mutation/genetics , Papua New Guinea , Receptor, ErbB-3/genetics , Sequence Analysis, DNAABSTRACT
Analyses of ancient DNA from extinct humans reveal signals of at least two independent hybridization events in the history of non-African populations. To date, there are very few examples of specific genetic variants that have been rigorously identified as introgressive. Here, we survey DNA sequence variation in the OAS gene cluster on chromosome 12 and provide strong evidence that a haplotype extending for ~185 kb introgressed from Neandertals. This haplotype is nearly restricted to Eurasians and is estimated to have diverged from the Neandertal sequence ~125 kya. Despite the potential for novel functional variation, the observed frequency of this haplotype is consistent with neutral introgression. This is the second locus in the human genome, after STAT2, carrying distinct haplotypes that appear to have introgressed separately from both Neandertals and Denisova.
Subject(s)
2',5'-Oligoadenylate Synthetase/genetics , Genetic Loci/immunology , Neanderthals/genetics , Polymorphism, Single Nucleotide , Animals , Evolution, Molecular , Genetic Speciation , Gorilla gorilla/genetics , Haplotypes , Humans , Linkage Disequilibrium , Models, Genetic , Multigene Family , Pan troglodytes/genetics , PhylogeographyABSTRACT
A long-debated question concerns the fate of archaic forms of the genus Homo: did they go extinct without interbreeding with anatomically modern humans, or are their genes present in contemporary populations? This question is typically focused on the genetic contribution of archaic forms outside of Africa. Here we use DNA sequence data gathered from 61 noncoding autosomal regions in a sample of three sub-Saharan African populations (Mandenka, Biaka, and San) to test models of African archaic admixture. We use two complementary approximate-likelihood approaches and a model of human evolution that involves recent population structure, with and without gene flow from an archaic population. Extensive simulation results reject the null model of no admixture and allow us to infer that contemporary African populations contain a small proportion of genetic material (≈ 2%) that introgressed ≈ 35 kya from an archaic population that split from the ancestors of anatomically modern humans ≈ 700 kya. Three candidate regions showing deep haplotype divergence, unusual patterns of linkage disequilibrium, and small basal clade size are identified and the distributions of introgressive haplotypes surveyed in a sample of populations from across sub-Saharan Africa. One candidate locus with an unusual segment of DNA that extends for >31 kb on chromosome 4 seems to have introgressed into modern Africans from a now-extinct taxon that may have lived in central Africa. Taken together our results suggest that polymorphisms present in extant populations introgressed via relatively recent interbreeding with hominin forms that diverged from the ancestors of modern humans in the Lower-Middle Pleistocene.
Subject(s)
Black People/genetics , Gene Flow/genetics , Hominidae/genetics , Africa , Animals , Base Sequence , Consanguinity , Genetic Loci/genetics , Genetics, Population , Geography , Humans , Likelihood Functions , Models, GeneticABSTRACT
While whole genome sequencing (WGS) of cell-free DNA (cfDNA) holds enormous promise for detection of molecular residual disease (MRD), its performance is limited by WGS error rate. Here we introduce AccuScan, an efficient cfDNA WGS technology that enables genome-wide error correction at single read-level, achieving an error rate of 4.2 × 10-7, which is about two orders of magnitude lower than a read-centric de-noising method. The application of AccuScan to MRD demonstrated analytical sensitivity down to 10-6 circulating variant allele frequency at 99% sample-level specificity. AccuScan showed 90% landmark sensitivity (within 6 weeks after surgery) and 100% specificity for predicting relapse in colorectal cancer. It also showed 67% sensitivity and 100% specificity in esophageal cancer using samples collected within one week after surgery. When AccuScan was applied to monitor immunotherapy in melanoma patients, the circulating tumor DNA (ctDNA) levels and dynamic profiles were consistent with clinical outcomes. Overall, AccuScan provides a highly accurate WGS solution for MRD detection, empowering ctDNA detection at parts per million range without requiring high sample input or personalized reagents.
ABSTRACT
Recent analysis of DNA extracted from two Eurasian forms of archaic human shows that more genetic variants are shared with humans currently living in Eurasia than with anatomically modern humans in sub-Saharan Africa. Although these genome-wide average measures of genetic similarity are consistent with the hypothesis of archaic admixture in Eurasia, analyses of individual loci exhibiting the signal of archaic introgression are needed to test alternative hypotheses and investigate the admixture process. Here, we provide a detailed sequence analysis of the innate immune gene OAS1, a locus with a divergent Melanesian haplotype that is very similar to the Denisova sequence from the Altai region of Siberia. We resequenced a 7-kb region encompassing the OAS1 gene in 88 individuals from six Old World populations (San, Biaka, Mandenka, French Basque, Han Chinese, and Papua New Guineans) and discovered previously unknown and ancient genetic variation. The 5' region of this gene has unusual patterns of diversity, including 1) higher levels of nucleotide diversity in Papuans than in sub-Saharan Africans, 2) very deep ancestry with an estimated time to the most recent common ancestor of >3 myr, and 3) a basal branching pattern with Papuan individuals on either side of the rooted network. A global geographic survey of >1,500 individuals showed that the divergent Papuan haplotype is nearly restricted to populations from eastern Indonesia and Melanesia. Polymorphic sites within this haplotype are shared with the draft Denisova genome over a span of â¼90 kb and are associated with an extended block of linkage disequilibrium, supporting the hypothesis that this haplotype introgressed from an archaic source that likely lived in Eurasia.
Subject(s)
2',5'-Oligoadenylate Synthetase/genetics , Black People/genetics , Hybridization, Genetic , Polymorphism, Genetic , 5' Untranslated Regions/genetics , Animals , Asian People/genetics , Chromosomes, Human, Pair 12/genetics , Genetic Speciation , Genotyping Techniques , Haplotypes , Hominidae/genetics , Humans , Linkage Disequilibrium , Melanesia , Models, Genetic , Native Hawaiian or Other Pacific Islander/genetics , Sequence Analysis, DNAABSTRACT
Sub-Saharan Africa has consistently been shown to be the most genetically diverse region in the world. Despite the fact that a substantial portion of this variation is partitioned between groups practicing a variety of subsistence strategies and speaking diverse languages, there is currently no consensus on the genetic relationships of sub-Saharan African populations. San (a subgroup of KhoeSan) and many Pygmy groups maintain hunter-gatherer lifestyles and cluster together in autosomal-based analysis, whereas non-Pygmy Niger-Kordofanian speakers (non-Pygmy NKs) predominantly practice agriculture and show substantial genetic homogeneity despite their wide geographic range throughout sub-Saharan Africa. However, KhoeSan, who speak a set of relatively unique click-based languages, have long been thought to be an early branch of anatomically modern humans based on phylogenetic analysis. To formally test models of divergence among the ancestors of modern African populations, we resequenced a sample of San, Eastern, and Western Pygmies and non-Pygmy NKs individuals at 40 nongenic (â¼2 kb) regions and then analyzed these data within an Approximate Bayesian Computation (ABC) framework. We find substantial support for a model of an early divergence of KhoeSan ancestors from a proto-Pygmy-non-Pygmy NKs group â¼110 thousand years ago over a model incorporating a proto-KhoeSan-Pygmy hunter-gatherer divergence from the ancestors of non-Pygmy NKs. The results of our analyses are consistent with previously identified signals of a strong bottleneck in Mbuti Pygmies and a relatively recent expansion of non-Pygmy NKs. We also develop a number of methodologies that utilize "pseudo-observed" data sets to optimize our ABC-based inference. This approach is likely to prove to be an invaluable tool for demographic inference using genome-wide resequencing data.
Subject(s)
Biological Evolution , Black People/genetics , Genetic Variation , Africa South of the Sahara , DNA, Mitochondrial/genetics , Genetics, Population , Humans , Phylogeny , Sequence Analysis, DNAABSTRACT
Increasing phylogenetic resolution of the Y chromosome haplogroup tree has led to finer temporal and spatial resolution for studies of human migration. Haplogroup T, initially known as K2 and defined by mutation M70, is found at variable frequencies across West Asia, Africa, and Europe. While several SNPs were recently discovered that extended the length of the branch leading to haplogroup T, only two SNPs are known to mark internal branches of haplogroup T. This low level of phylogenetic resolution has hindered studies of the origin and dispersal of this interesting haplogroup, which is found in Near Eastern non-Jewish populations, Jewish populations from several communities, and in the patrilineage of President Thomas Jefferson. Here we map 10 new SNPs that, together with the previously known SNPs, mark 11 lineages and two large subclades (T1a and T1b) of haplogroup T. We also report a new SNP that links haplogroups T and L within the major framework of Y chromosome evolution. Estimates of the timing of the branching events within haplogroup T, along with a comprehensive geographic survey of the major T subclades, suggest that this haplogroup began to diversify in the Near East -25 kya. Our survey also points to a complex history of dispersal of this rare and informative haplogroup within the Near East and from the Near East to Europe and sub-Saharan Africa. The presence of T1a2 chromosomes in Near Eastern Jewish and non-Jewish populations may reflect early exiles between the ancient lands of Israel and Babylon. The presence of different subclades of T chromosomes in Europe may be explained by both the spread of Neolithic farmers and the later dispersal of Jews from the Near East. Finally, the moderately high frequency (-18%) of T1b* chromosomes in the Lemba of southern Africa supports the hypothesis of a Near Eastern, but not necessarily a Jewish, origin for their paternal line.
Subject(s)
Chromosomes, Human, Y/genetics , Emigration and Immigration/statistics & numerical data , Genetics, Population , Haplotypes/genetics , Microsatellite Repeats/genetics , Africa , Bayes Theorem , Europe , Genotype , History, Ancient , Humans , Israel , Jews/genetics , Jews/history , Male , Middle East , Polymorphism, Single Nucleotide , TimeABSTRACT
Comparisons of levels of variability on the autosomes and X chromosome can be used to test hypotheses about factors influencing patterns of genomic variation. While a tremendous amount of nucleotide sequence data from across the genome is now available for multiple human populations, there has been no systematic effort to examine relative levels of neutral polymorphism on the X chromosome versus autosomes. We analyzed approximately 210 kb of DNA sequencing data representing 40 independent noncoding regions on the autosomes and X chromosome from each of 90 humans from six geographically diverse populations. We correct for differences in mutation rates between males and females by considering the ratio of within-human diversity to human-orangutan divergence. We find that relative levels of genetic variation are higher than expected on the X chromosome in all six human populations. We test a number of alternative hypotheses to explain the excess polymorphism on the X chromosome, including models of background selection, changes in population size, and sex-specific migration in a structured population. While each of these processes may have a small effect on the relative ratio of X-linked to autosomal diversity, our results point to a systematic difference between the sexes in the variance in reproductive success; namely, the widespread effects of polygyny in human populations. We conclude that factors leading to a lower male versus female effective population size must be considered as important demographic variables in efforts to construct models of human demographic history and for understanding the forces shaping patterns of human genomic variability.
Subject(s)
Evolution, Molecular , Genetic Variation , Chromosomes, Human/genetics , Chromosomes, Human, X/genetics , Female , Genetics, Population , Humans , Male , Models, Genetic , Mutation , Population Density , Selection, Genetic , Sequence Analysis, DNA , Sex CharacteristicsABSTRACT
A history of Pleistocene population expansion has been inferred from the frequency spectrum of polymorphism in the mitochondrial DNA (mtDNA) of many human populations. Similar patterns are not typically observed for autosomal and X-linked loci. One explanation for this discrepancy is a recent population bottleneck, with different rates of recovery for haploid and autosomal loci as a result of their different effective population sizes. This hypothesis predicts that mitochondrial and Y chromosomal DNA will show a similar skew in the frequency spectrum in populations that have experienced a recent increase in effective population size. We test this hypothesis by resequencing 6.6 kb of noncoding Y chromosomal DNA and 780 basepairs of the mtDNA cytochrome c oxidase subunit III (COIII) gene in 172 males from 5 African populations. Four tests of population expansion are employed for each locus in each population: Fu's Fs statistic, the R(2) statistic, coalescent simulations, and the mismatch distribution. Consistent with previous results, patterns of mtDNA polymorphism better fit a model of constant population size for food-gathering populations and a model of population expansion for food-producing populations. In contrast, none of the tests reveal evidence of Y chromosome growth for either food-gatherers or food-producers. The distinct mtDNA and Y chromosome polymorphism patterns most likely reflect sex-biased demographic processes in the recent history of African populations. We hypothesize that males experienced smaller effective population sizes and/or lower rates of migration during the Bantu expansion, which occurred over the last 5,000 years.
Subject(s)
Chromosomes, Human, Y/genetics , DNA, Mitochondrial/genetics , Electron Transport Complex IV/genetics , Polymorphism, Genetic , Population Density , Africa , Genetics, Population , Humans , MaleABSTRACT
It has been known for over a decade that a majority of men who self report as members of the Jewish priesthood (Cohanim) carry a characteristic Y chromosome haplotype termed the Cohen Modal Haplotype (CMH). The CMH has since been used to trace putative Jewish ancestral origins of various populations. However, the limited number of binary and STR Y chromosome markers used previously did not provide the phylogenetic resolution needed to infer the number of independent paternal lineages that are encompassed within the Cohanim or their coalescence times. Accordingly, we have genotyped 75 binary markers and 12 Y-STRs in a sample of 215 Cohanim from diverse Jewish communities, 1,575 Jewish men from across the range of the Jewish Diaspora, and 2,099 non-Jewish men from the Near East, Europe, Central Asia, and India. While Cohanim from diverse backgrounds carry a total of 21 Y chromosome haplogroups, 5 haplogroups account for 79.5% of Cohanim Y chromosomes. The most frequent Cohanim lineage (46.1%) is marked by the recently reported P58 T->C mutation, which is prevalent in the Near East. Based on genotypes at 12 Y-STRs, we identify an extended CMH on the J-P58* background that predominates in both Ashkenazi and non-Ashkenazi Cohanim and is remarkably absent in non-Jews. The estimated divergence time of this lineage based on 17 STRs is 3,190 +/- 1,090 years. Notably, the second most frequent Cohanim lineage (J-M410*, 14.4%) contains an extended modal haplotype that is also limited to Ashkenazi and non-Ashkenazi Cohanim and is estimated to be 4.2 +/- 1.3 ky old. These results support the hypothesis of a common origin of the CMH in the Near East well before the dispersion of the Jewish people into separate communities, and indicate that the majority of contemporary Jewish priests descend from a limited number of paternal lineages.
Subject(s)
Chromosomes, Human, Y/genetics , Jews/genetics , Judaism , Ethnicity/genetics , Family Characteristics , Gene Frequency , Genetic Markers , Geography , Haplotypes , Humans , Male , Microsatellite Repeats/genetics , Models, GeneticABSTRACT
A 2.4-kb stretch within the RRM2P4 region of the X chromosome, previously sequenced in a sample of 41 globally distributed humans, displayed both an ancient time to the most recent common ancestor (e.g., a TMRCA of approximately 2 million years) and a basal clade composed entirely of Asian sequences. This pattern was interpreted to reflect a history of introgressive hybridization from archaic hominins (most likely Asian Homo erectus) into the anatomically modern human genome. Here, we address this hypothesis by resequencing the 2.4-kb RRM2P4 region in 131 African and 122 non-African individuals and by extending the length of sequence in a window of 16.5 kb encompassing the RRM2P4 pseudogene in a subset of 90 individuals. We find that both the ancient TMRCA and the skew in non-African representation in one of the basal clades are essentially limited to the central 2.4-kb region. We define a new summary statistic called the minimum clade proportion (pmc), which quantifies the proportion of individuals from a specified geographic region in each of the two basal clades of a binary gene tree, and then employ coalescent simulations to assess the likelihood of the observed central RRM2P4 genealogy under two alternative views of human evolutionary history: recent African replacement (RAR) and archaic admixture (AA). A molecular-clock-based TMRCA estimate of 2.33 million years is a statistical outlier under the RAR model; however, the large variance associated with this estimate makes it difficult to distinguish the predictions of the human origins models tested here. The pmc summary statistic, which has improved power with larger samples of chromosomes, yields values that are significantly unlikely under the RAR model and fit expectations better under a range of archaic admixture scenarios.
Subject(s)
Chromosomes, Human, X/genetics , Genealogy and Heraldry , Models, Genetic , DNA, Intergenic/genetics , Demography , Genetic Variation , Humans , Likelihood Functions , Phylogeny , Sequence Analysis, DNAABSTRACT
We report the sequences of 1,244 human Y chromosomes randomly ascertained from 26 worldwide populations by the 1000 Genomes Project. We discovered more than 65,000 variants, including single-nucleotide variants, multiple-nucleotide variants, insertions and deletions, short tandem repeats, and copy number variants. Of these, copy number variants contribute the greatest predicted functional impact. We constructed a calibrated phylogenetic tree on the basis of binary single-nucleotide variants and projected the more complex variants onto it, estimating the number of mutations for each class. Our phylogeny shows bursts of extreme expansion in male numbers that have occurred independently among each of the five continental superpopulations examined, at times of known migrations and technological innovations.
Subject(s)
Chromosomes, Human, Y , Demography , Haplotypes , Humans , Male , Mutation , Phylogeny , Polymorphism, Single NucleotideABSTRACT
The highly structured distribution of Y-chromosome haplogroups suggests that current patterns of variation may be informative of past population processes. However, limited phylogenetic resolution, particularly of subclades within haplogroup K, has obscured the relationships of lineages that are common across Eurasia. Here we genotype 13 new highly informative single-nucleotide polymorphisms in a worldwide sample of 4413 males that carry the derived allele at M526, and reconstruct an NRY haplogroup tree with significantly higher resolution for the major clade within haplogroup K, K-M526. Although K-M526 was previously characterized by a single polytomy of eight major branches, the phylogenetic structure of haplogroup K-M526 is now resolved into four major subclades (K2a-d). The largest of these subclades, K2b, is divided into two clusters: K2b1 and K2b2. K2b1 combines the previously known haplogroups M, S, K-P60 and K-P79, whereas K2b2 comprises haplogroups P and its subhaplogroups Q and R. Interestingly, the monophyletic group formed by haplogroups R and Q, which make up the majority of paternal lineages in Europe, Central Asia and the Americas, represents the only subclade with K2b that is not geographically restricted to Southeast Asia and Oceania. Estimates of the interval times for the branching events between M9 and P295 point to an initial rapid diversification process of K-M526 that likely occurred in Southeast Asia, with subsequent westward expansions of the ancestors of haplogroups R and Q.
Subject(s)
Chromosomes, Human, Y , Genetic Variation , Genetics, Population , Haplotypes , Phylogeny , Asia, Southeastern , Ethnicity/genetics , Gene Frequency , Humans , Male , Mutation , Phylogeography , Polymorphism, GeneticABSTRACT
The ratio of X-linked to autosomal diversity was estimated from an analysis of six human genome sequences and found to deviate from the expected value of 0.75. However, the direction of this deviation depends on whether a particular sequence is close to or far from the nearest gene. This pattern may be explained by stronger locally acting selection on X-linked genes compared with autosomal genes, combined with larger effective population sizes for females than for males.