Olga L. Posukh

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

Here we report the Simons Genome Diversity Project data set: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioural modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that of other non-Africans.

Genomic evidence for the Pleistocene and recent population history of Native Americans

Genetic history of Native Americans Several theories have been put forth as to the origin and timing of when Native American ancestors entered the Americas. To clarify this controversy, Raghavan et al. examined the genomic variation among ancient and modern individuals from Asia and the Americas. There is no evidence for multiple waves of entry or recurrent gene flow with Asians in northern populations. The earliest migrations occurred no earlier than 23,000 years ago from Siberian ancestors. Amerindians and Athabascans originated from a single population, splitting approximately 13,000 years ago. Science, this issue 10.1126/science.aab3884 Genetic variation within ancient and extant Native American populations informs on their migration into the Americas. INTRODUCTION The consensus view on the peopling of the Americas is that ancestors of modern Native Americans entered the Americas from Siberia via the Bering Land Bridge and that this occurred at least ~14.6 thousand years ago (ka). However, the number and timing of migrations into the Americas remain controversial, with conflicting interpretations based on anatomical and genetic evidence. RATIONALE In this study, we address four major unresolved issues regarding the Pleistocene and recent population history of Native Americans: (i) the timing of their divergence from their ancestral group, (ii) the number of migrations into the Americas, (iii) whether there was ~15,000 years of isolation of ancestral Native Americans in Beringia (Beringian Incubation Model), and (iv) whether there was post-Pleistocene survival of relict populations in the Americas related to Australo-Melanesians, as suggested by apparent differences in cranial morphologies between some early (“Paleoamerican”) remains and those of more recent Native Americans. We generated 31 high-coverage modern genomes from the Americas, Siberia, and Oceania; 23 ancient genomic sequences from the Americas dating between ~0.2 and 6 ka; and SNP chip genotype data from 79 present-day individuals belonging to 28 populations from the Americas and Siberia. The above data sets were analyzed together with published modern and ancient genomic data from worldwide populations, after masking some present-day Native Americans for recent European admixture. RESULTS Using three different methods, we determined the divergence time for all Native Americans (Athabascans and Amerindians) from their Siberian ancestors to be ~20 ka, and no earlier than ~23 ka. Furthermore, we dated the divergence between Athabascans (northern Native American branch, together with northern North American Amerindians) and southern North Americans and South and Central Americans (southern Native American branch) to be ~13 ka. Similar divergence times from East Asian populations and a divergence time between the two branches that is close in age to the earliest well-established archaeological sites in the Americas suggest that the split between the branches occurred within the Americas. We additionally found that several sequenced Holocene individuals from the Americas are related to present-day populations from the same geographical regions, implying genetic continuity of ancient and modern populations in some parts of the Americas over at least the past 8500 years. Moreover, our results suggest that there has been gene flow between some Native Americans from both North and South America and groups related to East Asians and Australo-Melanesians, the latter possibly through an East Asian route that might have included ancestors of modern Aleutian Islanders. Last, using both genomic and morphometric analyses, we found that historical Native American groups such as the Pericúes and Fuego-Patagonians were not “relicts” of Paleoamericans, and hence, our results do not support an early migration of populations directly related to Australo-Melanesians into the Americas. CONCLUSION Our results provide an upper bound of ~23 ka on the initial divergence of ancestral Native Americans from their East Asian ancestors, followed by a short isolation period of no more than ~8000 years, and subsequent entrance and spread across the Americas. The data presented are consistent with a single-migration model for all Native Americans, with later gene flow from sources related to East Asians and, indirectly, Australo-Melanesians. The single wave diversified ~13 ka, likely within the Americas, giving rise to the northern and southern branches of present-day Native Americans. Population history of present-day Native Americans. The ancestors of all Native Americans entered the Americas as a single migration wave from Siberia (purple) no earlier than ~23 ka, separate from the Inuit (green), and diversified into “northern” and “southern” Native American branches ~13 ka. There is evidence of post-divergence gene flow between some Native Americans and groups related to East Asians/Inuit and Australo-Melanesians (yellow). How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative “Paleoamerican” relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

Global diversity, population stratification, and selection of human copy-number variation.

Duplications and deletions in the human genome Duplications and deletions can lead to variation in copy number for genes and genomic loci among humans. Such variants can reveal evolutionary patterns and have implications for human health. Sudmant et al. examined copy-number variation across 236 individual genomes from 125 human populations. Deletions were under more selection, whereas duplications showed more population-specific structure. Interestingly, Oceanic populations retain large duplications postulated to have originated in an ancient Denisovan lineage. Science, this issue 10.1126/science.aab3761 Copy-number variation reveals how selection affects the human genome across the globe. INTRODUCTION Most studies of human genetic variation have focused on single-nucleotide variants (SNVs). However, copy-number variants (CNVs) affect more base pairs of DNA among humans, and yet our understanding of CNV diversity among human populations is limited. RATIONALE We aimed to understand the pattern, selection, and diversity of copy-number variation by analyzing deeply sequenced genomes representing the diversity of all humans. We compared the selective constraints of deletions versus duplications to understand population stratification in the context of the ancestral human genome and to assess differences in CNV load between African and non-African populations. RESULTS We sequenced 236 individual genomes from 125 distinct human populations and identified 14,467 autosomal CNVs and 545 X-linked CNVs with a sequence read-depth approach. Deletions exhibit stronger selective pressure and are better phylogenetic markers of population relationships than duplication polymorphisms. We identified 1036 population-stratified copy-number–variable regions, 295 of which intersect coding regions and 199 of which exhibit extreme signatures of differentiation. Duplicated loci were 1.8-fold more likely to be stratified than deletions but were poorly correlated with flanking genetic diversity. Among these, we highlight a duplication polymorphism restricted to modern Oceanic populations yet also present in the genome of the archaic Denisova hominin. This 225–kilo–base pair (kbp) duplication includes two microRNA genes and is almost fixed among human Papuan-Bougainville genomes. The data allowed us to reconstruct the ancestral human genome and create a more accurate evolutionary framework for the gain and loss of sequences during human evolution. We identified 571 loci that segregate in the human population and another 2026 loci of fixed-copy 2 in all human genomes but absent from the reference genome. The total deletion and duplication load between African and non-African population groups showed no difference after we account for ancestral sequences missing from the human reference. However, we did observe that the relative number of base pairs affected by CNVs compared to single-nucleotide polymorphisms is higher among non-Africans than Africans. CONCLUSION Deletions, duplications, and CNVs have shaped, to different extents, the genetic diversity of human populations by the combined forces of mutation, selection, and demography. Figure Global human CNV diversity and archaic introgression of a chromosome 16 duplication. (Left) The geographic coordinates of populations sampled are indicated on a world map (colored dots). The pie charts show the continental population allele frequency of a single ~225-kbp duplication polymorphism found exclusively among Oceanic populations and an archaic Denisova. (Right) The ancestral structure of this duplication locus (1) and the Denisova duplication structure (2) are shown in relation to their position on chromosome 16. We estimate that the duplication emerged ~440 thousand years ago (ka) in the Denisova and then introgressed into ancestral Papuan populations ~40 ka. In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide–variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

High Levels of Y-Chromosome Differentiation among Native Siberian Populations and the Genetic Signature of a Boreal Hunter-Gatherer Way of Life

We examined genetic variation on the nonrecombining portion of the Y chromosome (NRY) to investigate the paternal population structure of indigenous Siberian groups and to reconstruct the historical events leading to the peopling of Siberia. A set of 62 biallelic markers on the NRY were genotyped in 1432 males representing 18 Siberian populations, as well as nine populations from Central and East Asia and one from European Russia. A subset of these markers defines the 18 major NRY haplogroups (A-R) recently described by the Y Chromosome Consortium (YCC 2002). While only four of these 18 major NRY haplogroups accounted for ~ 95% of Siberian Y-chromosome variation, native Siberian populations differed greatly in their haplogroup composition and exhibited the highest F ST value for any region of the world. When we divided our Siberian sample into four geographic regions versus five major linguistic groupings, analyses of molecular variance (AMOVA) indicated higher F ST and F CT values for linguistic groups than for geographic groups. Mantel tests also supported the existence of NRY genetic patterns that were correlated with language, indicating that language affiliation might be a better predictor of the genetic affinity among Siberians than their present geographic position. The combined results, including those from a nested cladistic analysis, underscored the important role of directed dispersals, range expansions, and long-distance colonizations bound by common ethnic and linguistic affiliation in shaping the genetic landscape of Siberia. The Siberian pattern of reduced haplogroup diversity within populations combined with high levels of differentiation among populations may be a general feature characteristic of indigenous groups that have small effective population sizes and that have been isolated for long periods of time.

Y Chromosome Markers and Trans-Bering Strait Dispersals

Five polymorphisms involving two paternally inherited loci were surveyed in 38 world populations (n = 1,631) to investigate the origins of Native Americans. One of the six Y chromosome combination haplotypes (1T) was found at relatively high frequencies (17.8-75.0%) in nine Native American populations (n = 206) representing the three major linguistic divisions in the New World. Overall, these data do not support the Greenberg et al. (1986) tripartite model for the early peopling of the Americas. The 1T haplotype was also discovered at a low frequency in Siberian Eskimos (3/22), Chukchi (1/6), and Evens (1/65) but was absent from 17 other Asian populations (n = 987). The perplexing presence of the 1T haplotype in northeastern Siberia may be due to back-migration from the New World to Asia.

First molecular screening of deafness in the Altai Republic population

BackgroundWe studied the molecular basis of NSHL in Republic of Altai (South Siberia, Russia). The Altaians are the indigenous Asian population of the Altai Mountain region considered as a melting-pot and a dispersion center for world-wide human expansions in the past.MethodsA total of 76 patients of Altaian, Russian or mixed ethnicity and 130 Altaian controls were analyzed by PCR-DHPLC and sequencing in the GJB2 gene. The GJB6 deletion and the common non-syndromic deafness-causing mitochondrial mutations were also tested when appropriate.Results8.3% of the Altaian chromosomes were carrying GJB2 mutations versus 46.9% of the Russian chromosomes. The 235delC mutation was predominant among Altaians, whereas the 35delG mutation was most prevalent among Russian patients.ConclusionWe found an Asian-specific GJB2 diversity among Altaians, and different GJB2 contribution for deafness in the Altaian and Russian patients. The high carrier frequency of 235delC in Altaians (4.6%) is probably defined by gene drift/founder effect in a particular group. The question whether the Altai region could be one of founder sources for the 235delC mutation widespread in Asia is open.

Autosomal recessive deafness 1A (DFNB1A) in Yakut population isolate in Eastern Siberia: extensive accumulation of the splice site mutation IVS1+1G>A in GJB2 gene as a result of founder effect.

Hereditary forms of hearing impairment (HI) caused by GJB2 (Cx26) mutations are the frequent sensory disorders registered among newborns in various human populations. In this study, we present data on the molecular, audiological and population features of autosomal recessive deafness 1A (DFNB1A) associated with the donor splicing site IVS1+1G>A mutation of GJB2 gene in Yakut population isolate of the Sakha Republic (Yakutia) located in Eastern Siberia (Russian Federation). The Yakut population exhibits high frequency of some Mendelian disorders, which are rare in other populations worldwide. Mutational analysis of GJB2 gene in 86 unrelated Yakut patients with congenital HI without other clinical features has been performed. In this study, we registered a large cohort of Yakut patients homozygous for the IVS1+1G>A mutation (70 unrelated deaf subjects in total). Detailed audiological analysis of 40 deaf subjects with genotype IVS1+1G>A/IVS1+1G>A revealed significant association of this genotype with mostly symmetrical bilateral severe to profound HI (85% severe-to-profound HI versus 15% mild-to-moderate HI, P<0.05). The highest among six investigated Eastern Siberian populations carrier frequency of the IVS1+1G>A mutation (11.7%) has been found in Yakut population. Reconstruction of 140 haplotypes with IVS1+1G>A mutation demonstrates the common origin of all mutant chromosomes found in Yakuts. The age of mutation was estimated to be approximately 800 years. These findings characterize Eastern Siberia as the region with the most extensive accumulation of the IVS1+1G>A mutation in the world as a result of founder effect.

Carrier frequency of GJB2 gene mutations c.35delG, c.235delC and c.167delT among the populations of Eurasia.

Hearing impairment is one of the most common disorders of sensorineural function and the incidence of profound prelingual deafness is about 1 per 1000 at birth. GJB2 gene mutations make the largest contribution to hereditary hearing impairment. The spectrum and prevalence of some GJB2 mutations are known to be dependent on the ethnic origin of the population. This study presents data on the carrier frequencies of major GJB2 mutations, c.35delG, c.167delT and c.235delC, among 2308 healthy persons from 18 various populations of Eurasia: Russians, Bashkirs, Tatars, Chuvashes, Udmurts, Komi-Permyaks and Mordvins (Volga-Ural region of Russia); Belarusians and Ukrainians (East Europe); Abkhazians, Avars, Cherkessians and Ingushes (Caucasus); Kazakhs, Uighurs and Uzbeks (Central Asia); and Yakuts and Altaians (Siberia). The data on c.35delG and c.235delC mutation prevalence in the studied ethnic groups can be used to investigate the prospective founder effect in the origin and prevalence of these mutations in Eurasia and consequently in populations around the world.

Russian Old Believers: Genetic Consequences of Their Persecution and Exile, as Shown by Mitochondrial DNA Evidence

Abstract In 1653, the Patriarch Nikon modified liturgical practices to bring the Russian Orthodox Church in line with those of the Eastern (Greek) Orthodox Church, from which it had split 200 years earlier. The Old Believers (staroveri) rejected these changes and continued to worship using the earlier practices. These actions resulted in their persecution by the Russian Orthodox Church, which forced them into exile across Siberia. Given their history, we investigate whether populations of Old Believers have diverged genetically from other Slavic populations as a result of their isolation. We also examine whether the three Old Believer populations analyzed in this study are part of a single gene pool (founder population) or are instead derived from heterogeneous sources. As part of this analysis, we survey the mitochondrial DNAs (mtDNAs) of 189 Russian Old Believer individuals from three populations in Siberia and 201 ethnic Russians from different parts of Siberia for phylogenetically informative mutations in the coding and noncoding regions. Our results indicate that the Old Believers have not significantly diverged genetically from other Slavic populations over the 200–300 years of their isolation in Siberia. However, they do show some unique patterns of mtDNA variation relative to other Slavic groups, such as a high frequency of subhaplogroup U4, a surprisingly low frequency of haplogroup H, and low frequencies of the rare East Eurasian subhaplogroup D5.

Spectrum and Frequency of the GJB2 Gene Pathogenic Variants in a Large Cohort of Patients with Hearing Impairment Living in a Subarctic Region of Russia (the Sakha Republic)

Pathogenic variants in the GJB2 gene, encoding connexin 26, are known to be a major cause of hearing impairment (HI). More than 300 allelic variants have been identified in the GJB2 gene. Spectrum and allelic frequencies of the GJB2 gene vary significantly among different ethnic groups worldwide. Until now, the spectrum and frequency of the pathogenic variants in exon 1, exon 2 and the flanking intronic regions of the GJB2 gene have not been described thoroughly in the Sakha Republic (Yakutia), which is located in a subarctic region in Russia. The complete sequencing of the non-coding and coding regions of the GJB2 gene was performed in 393 patients with HI (Yakuts—296, Russians—51, mixed and other ethnicities—46) and in 187 normal hearing individuals of Yakut (n = 107) and Russian (n = 80) populations. In the total sample (n = 580), we revealed 12 allelic variants of the GJB2 gene, 8 of which were recessive pathogenic variants. Ten genotypes with biallelic recessive pathogenic variants in the GJB2 gene (in a homozygous or a compound heterozygous state) were found in 192 out of 393 patients (48.85%). We found that the most frequent GJB2 pathogenic variant in the Yakut patients was c.-23+1G>A (51.82%) and that the second most frequent was c.109G>A (2.37%), followed by c.35delG (1.64%). Pathogenic variants с.35delG (22.34%), c.-23+1G>A (5.31%), and c.313_326del14 (2.12%) were found to be the most frequent among the Russian patients. The carrier frequencies of the c.-23+1G>A and с.109G>A pathogenic variants in the Yakut control group were 10.20% and 2.80%, respectively. The carrier frequencies of с.35delG and c.101T>C were identical (2.5%) in the Russian control group. We found that the contribution of the GJB2 gene pathogenic variants in HI in the population of the Sakha Republic (48.85%) was the highest among all of the previously studied regions of Asia. We suggest that extensive accumulation of the c.-23+1G>A pathogenic variant in the indigenous Yakut population (92.20% of all mutant chromosomes in patients) and an extremely high (10.20%) carrier frequency in the control group may indicate a possible selective advantage for the c.-23+1G>A carriers living in subarctic climate.