Search

R1a-M420 is one of the most widely spread Y-chromosome haplogroups; however, its substructure within Europe and Asia has remained poorly characterized. Using a panel of 16 244 male subjects from 126 populations sampled across Eurasia, we identified 2923 R1a-M420 Y-chromosomes and analyzed them to a highly granular phylogeographic resolution. Whole Y-chromosome sequence analysis of eight R1a and five R1b individuals suggests a divergence time of ∼25 000 (95% CI: 21 300-29 000) years ago and a coalescence time within R1a-M417 of ∼5800 (95% CI: 4800-6800) years. The spatial frequency distributions of R1a sub-haplogroups conclusively indicate two major groups, one found primarily in Europe and the other confined to Central and South Asia. Beyond the major European versus Asian dichotomy, we describe several younger sub-haplogroups. Based on spatial distributions and diversity patterns within the R1a-M420 clade, particularly rare basal branches detected primarily within Iran and eastern Turkey, we conclude that the initial episodes of haplogroup R1a diversification likely occurred in the vicinity of present-day Iran. [ABSTRACT FROM AUTHOR]

The Genographic Project is studying the genetic signatures of ancient human migrations and creating an open-source research database. It allows members of the public to participate in a real-time anthropological genetics study by submitting personal samples for analysis and donating the genetic results to the database. We report our experience from the first 18 months of public participation in the Genographic Project, during which we have created the largest standardized human mitochondrial DNA (mtDNA) database ever collected, comprising 78,590 genotypes. Here, we detail our genotyping and quality assurance protocols including direct sequencing of the mtDNA HVS-I, genotyping of 22 coding-region SNPs, and a series of computational quality checks based on phylogenetic principles. This database is very informative with respect to mtDNA phylogeny and mutational dynamics, and its size allows us to develop a nearest neighbor-based methodology for mtDNA haplogroup prediction based on HVS-I motifs that is superior to classic rule-based approaches. We make available to the scientific community and general public two new resources: a periodically updated database comprising all data donated by participants, and the nearest neighbor haplogroup prediction tool. [ABSTRACT FROM AUTHOR]

Background: Recent advances in the understanding of the maternal and paternal heritage of south and southwest Asian populations have highlighted their role in the colonization of Eurasia by anatomically modern humans. Further understanding requires a deeper insight into the topology of the branches of the Indian mtDNA phylogenetic tree, which should be contextualized within the phylogeography of the neighboring regional mtDNA variation. Accordingly, we have analyzed mtDNA control and coding region variation in 796 Indian (including both tribal and caste populations from different parts of India) and 436 Iranian mtDNAs. The results were integrated and analyzed together with published data from South, Southeast Asia and West Eurasia. Results: Four new Indian-specific haplogroup M sub-clades were defined. These, in combination with two previously described haplogroups, encompass approximately one third of the haplogroup M mtDNAs in India. Their phylogeography and spread among different linguistic phyla and social strata was investigated in detail. Furthermore, the analysis of the Iranian mtDNA pool revealed patterns of limited reciprocal gene flow between Iran and the Indian sub-continent and allowed the identification of different assemblies of shared mtDNA sub-clades. Conclusions: Since the initial peopling of South and West Asia by anatomically modern humans, when this region may well have provided the initial settlers who colonized much of the rest of Eurasia, the gene flow in and out of India of the maternally transmitted mtDNA has been surprisingly limited. Specifically, our analysis of the mtDNA haplogroups, which are shared between Indian and Iranian populations and exhibit coalescence ages corresponding to around the early Upper Paleolithic, indicates that they are present in India largely as Indian-specific sub-lineages. In contrast, other ancient Indian-specific variants of M and R are very rare outside the sub-continent. [ABSTRACT FROM AUTHOR]

Key Clinical Message A patient with syndromic Duane retraction syndrome harbors a chromosome 811.1q13.2 inversion and 8p11.1-q12.3 marker chromosome containing subregions with differing mosaicism and allele frequencies. This case highlights the potential requirement for multiple genetic methods to gain insight into genotype–phenotype correlation, and ultimately into molecular mechanisms that underlie human disease.

A correction to the article "Most of the Extant mtDNA Boundaries in South and Southwest Asia Were Likely Shaped During the Initial Settlement of Eurasia by Anatomically Modern Humans," by Mait Metspalu et al is presented.

A response by Anna Olivieri and colleagues to a letter to the editor about their article "The mtDNA Legacy of the Levantine Early Upper Paleolithic in Africa" in the December 15, 2006 issue is presented.

We set out to identify the origins of the Árpád Dynasty based on genome sequencing of DNA derived from the skeletal remains of Hungarian King Béla III (1172-1196) and eight additional individuals (six males, two females) originally interred at the Royal Basilica of Székesfehérvár. Y-chromosome analysis established that two individuals, Béla III and HU52 assign to haplogroups R-Z2125 whose distribution centres near South Central Asia with subsidiary expansions in the regions of modern Iran, the Volga Ural region and the Caucasus. Out of a cohort of 4340 individuals from these geographic areas, we acquired whole-genome data from 208 individuals derived for the R-Z2123 haplogroup. From these data we have established that the closest living kin of the Árpád Dynasty are R-SUR51 derived modern day Bashkirs predominantly from the Burzyansky and Abzelilovsky districts of Bashkortostan in the Russian Federation. Our analysis also reveals the existence of SNPs defining a novel Árpád Dynasty specific haplogroup R-ARP. Framed within the context of a high resolution R-Z2123 phylogeny, the ancestry of the first Hungarian royal dynasty traces to the region centering near Northern Afghanistan about 4500 years ago and identifies the Bashkirs as their closest kin, with a separation date between the two populations at the beginning of the first millennium CE.

Purpose: Increased implementation of complex genetic technologies in clinical practice emphasizes the urgency of genomic literacy and proficiency for medical professionals. We evaluated our genomic education model., Methods: We assessed the 5-day, extended format program, encompassing lectures, videos, interactive tests, practice cases, and clinical exercises. Pre- and post questionnaires assessed knowledge change, using t-tests to compare groups. Satisfaction on program completion and after 3 years were evaluated. Implementation in other centers determined acceptability., Results: During 2012-2018, 774 clinicians from multiple disciplines and career stages attended 35 programs; 334 (43%) attended the 5-day extended format. Evaluations showed significant improvement of genomic literacy (mean 15.05/100 points, p < 0.001). Residents initially had higher scores than specialists (pre: 66.3 ± 17.3 vs. 58.7 ± 16.6, respectively, p = 0.002); both significantly improved, with specialists "catching up" (post: 79.1 ± 17.2 vs. 75.7 ± 15.9, nonsignificant (NS)); there was a similar trend between fellows and subspecialists (pre: 70 ± 18 vs. 59.4 ± 16.4, respectively, p = 0.007; post: 78.6 ± 16.4 vs. 73.2 ± 17.7, respectively, NS). Younger specialists (≤10 years residency) had significantly higher pre- and post scores. Absolute improvement in scores did not depend on medical specialties., Conclusion: Our program is effective in improving genomics literacy for clinicians, irrespective of career length or expertise, and could be a model for improving skills in practical genomics for all medical professionals.

Purpose: Mutations in the gene HSD17B3 encoding the 17-beta hydroxysteroid dehydrogenase 3 enzyme cause testosterone insufficiency leading to XY disorders of sex development. In this study the clinical and molecular characteristics of three patients from consanguineous families are elucidated., Methods: We identified three patients from two unrelated families with XY DSD and a novel homozygous HSD17B3:c. 673G>A mutation. The effect of the mutation on splicing was determined in RNA extracted from the testis of one patient., Results: Three patients presented at ages 0.1, 8 and 0.7 years with ambiguous genitalia and an XY Karyotype. Endocrine workup showed normal cortisol and mineralocorticoid levels with a low testosterone/androstenedione ratio. Whole-exome sequencing, carried out in the first family, revealed a homozygous novel mutation in the HSD17B3 gene: c. 673G>A, p. V225M. The same mutation was found by Sanger sequencing in the third unrelated patient. Haplotype analysis of a 4 Mb region surrounding the HSD17B3 gene on chromosome 9 revealed that the mutation resides on the same allele in all three patients. The mutation, being the first nucleic acid on exon 10, affects splicing and causes exon 10 skipping in one of our patients' testes., Conclusion: The novel homozygous c. 673G>A, p. V225M mutation in the 17HSDB3 gene is likely a founder mutation and causes severe XY-DSD. It changes a conserved amino acid residue, and also alters 17HSDB3 gene transcription by causing skipping of exon 10, thereby contributing to an imbalance in the relevant protein isoforms and consequently, significant decreased 17HDSB3 enzymatic activity.

Purpose: Expanded preconception carrier screening (ECS) identifies at-risk couples (ARCs) for multiple diseases. ECS reports currently include only pathogenic/likely pathogenic variants (P/LPVs). Variants of unknown significance (VUS) are not reported, unlike genomic or chromosomal array test results in other post/prenatal settings. Couples who are P/LP and VUS carriers (P/LP*VUS) may be at risk, particularly in genes with high P/LP carrier rates. We examined the possible contribution of P/LP*cVUS (coding, nonsynonymous VUS) matings to ECS yield in an Ashkenazi Jewish cohort, a population with well-established preconception screening., Methods: We analyzed 672 Ashkenazi Jewish genome sequences (225,456 virtual matings) for variants in three different gene sets and calculated the rates of P/LP*P/LP and P/LP*cVUS matings., Results: Across 180 genes, we identified 4671 variants: 144 (3.1%) P/LP and 1963 (42%) VUS. Across gene sets, the proportion of P/LP*P/LP and P/LP*cVUS ARCs was 2.7-3.8% and 6.8-7.5%, respectively., Conclusion: Disregarding VUS in ECS may miss ARCs. Even if only 10% of couples currently classified as P/LP*cVUS are ultimately reclassified as P/LP*P/LP, ECS yield would increase by ≈20%. While current understanding of VUS precludes VUS reporting in ECS, these findings underscore the importance of VUS reclassification. This will crucially depend on enlarging population frequency databases, especially of affected individuals.

Objectives: The gene encoding glucose transporter 3 (GLUT3, SLC2A3 ) is present in the human population at variable copy number. An overt disease phenotype of SLC2A3 copy number variants has not been reported; however, deletion of SLC2A3 has been previously reported to protect carriers from rheumatoid arthritis, implicating GLUT3 as a therapeutic target in rheumatoid arthritis. Here we aim to perform functional analysis of GLUT3 copy number variants in immune cells, and test the reported protective association of the GLUT3 copy number variants for rheumatoid arthritis in a genetic replication study., Methods: Cells from genotyped healthy controls were analyzed for SLC2A3 /GLUT3 expression and glycolysis capacity. We genotyped the SLC2A3 copy number variant in four independent cohorts of rheumatoid arthritis and controls and one cohort of multiple sclerosis and controls., Results: Heterozygous deletion of SLC2A3 correlates directly with expression levels of GLUT3 and influences glycolysis rates in the human immune system. The frequency of the SLC2A3 copy number variant is not different between rheumatoid arthritis, multiple sclerosis and control groups., Conclusions: Despite a robust SLC2A3 gene copy number dependent phenotype, our study of large groups of rheumatoid arthritis cases and controls provides no evidence for rheumatoid arthritis disease protection in deletion carriers. These data emphasize the importance of well powered replication studies to confirm or refute genetic associations, particularly for relatively rare variants.

Surrounded by speakers of Indo-European, Dravidian and Tibeto-Burman languages, around 11 million Munda (a branch of Austroasiatic language family) speakers live in the densely populated and genetically diverse South Asia. Their genetic makeup holds components characteristic of South Asians as well as Southeast Asians. The admixture time between these components has been previously estimated on the basis of archaeology, linguistics and uniparental markers. Using genome-wide genotype data of 102 Munda speakers and contextual data from South and Southeast Asia, we retrieved admixture dates between 2000-3800 years ago for different populations of Munda. The best modern proxies for the source populations for the admixture with proportions 0.29/0.71 are Lao people from Laos and Dravidian speakers from Kerala in India. The South Asian population(s), with whom the incoming Southeast Asians intermixed, had a smaller proportion of West Eurasian genetic component than contemporary proxies. Somewhat surprisingly Malaysian Peninsular tribes rather than the geographically closer Austroasiatic languages speakers like Vietnamese and Cambodians show highest sharing of IBD segments with the Munda. In addition, we affirmed that the grouping of the Munda speakers into North and South Munda based on linguistics is in concordance with genome-wide data.

The rugged topography of the Himalayan region has hindered large-scale human migrations, population admixture and assimilation. Such complexity in geographical structure might have facilitated the existence of several small isolated communities in this region. We have genotyped about 850,000 autosomal markers among 35 individuals belonging to the four major populations inhabiting the Himalaya and adjoining regions. In addition, we have genotyped 794 individuals belonging to 16 ethnic groups from the same region, for uniparental (mitochondrial and Y chromosomal DNA) markers. Our results in the light of various statistical analyses suggest a closer link of the Himalayan and adjoining populations to East Asia than their immediate geographical neighbours in South Asia. Allele frequency-based analyses likely support the existence of a specific ancestry component in the Himalayan and adjoining populations. The admixture time estimate suggests a recent westward migration of populations living to the East of the Himalaya. Furthermore, the uniparental marker analysis among the Himalayan and adjoining populations reveal the presence of East, Southeast and South Asian genetic signatures. Interestingly, we observed an antagonistic association of Y chromosomal haplogroups O3 and D clines with the longitudinal distance. Thus, we summarise that studying the Himalayan and adjoining populations is essential for a comprehensive reconstruction of the human evolutionary and ethnolinguistic history of eastern Eurasia.

High-coverage whole-genome sequence studies have so far focused on a limited number of geographically restricted populations, or been targeted at specific diseases, such as cancer. Nevertheless, the availability of high-resolution genomic data has led to the development of new methodologies for inferring population history and refuelled the debate on the mutation rate in humans. Here we present the Estonian Biocentre Human Genome Diversity Panel (EGDP), a dataset of 483 high-coverage human genomes from 148 populations worldwide, including 379 new genomes from 125 populations, which we group into diversity and selection sets. We analyse this dataset to refine estimates of continent-wide patterns of heterozygosity, long- and short-distance gene flow, archaic admixture, and changes in effective population size through time as well as for signals of positive or balancing selection. We find a genetic signature in present-day Papuans that suggests that at least 2% of their genome originates from an early and largely extinct expansion of anatomically modern humans (AMHs) out of Africa. Together with evidence from the western Asian fossil record, and admixture between AMHs and Neanderthals predating the main Eurasian expansion, our results contribute to the mounting evidence for the presence of AMHs out of Africa earlier than 75,000 years ago., Competing Interests: The authors declare no competing financial interests.

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.

The paternal haplogroup (hg) N is distributed from southeast Asia to eastern Europe. The demographic processes that have shaped the vast extent of this major Y chromosome lineage across numerous linguistically and autosomally divergent populations have previously been unresolved. On the basis of 94 high-coverage re-sequenced Y chromosomes, we establish and date a detailed hg N phylogeny. We evaluate geographic structure by using 16 distinguishing binary markers in 1,631 hg N Y chromosomes from a collection of 6,521 samples from 56 populations. The more southerly distributed sub-clade N4 emerged before N2a1 and N3, found mostly in the north, but the latter two display more elaborate branching patterns, indicative of regional contrasts in recent expansions. In particular, a number of prominent and well-defined clades with common N3a3'6 ancestry occur in regionally dissimilar northern Eurasian populations, indicating almost simultaneous regional diversification and expansion within the last 5,000 years. This patrilineal genetic affinity is decoupled from the associated higher degree of language diversity., (Copyright © 2016 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Inherited optic neuropathies are a heterogeneous group of disorders characterized by mild to severe visual loss, colour vision deficit, central or paracentral visual field defects and optic disc pallor. Optic atrophies can be classified into isolated or non-syndromic and syndromic forms. While multiple modes of inheritance have been reported, autosomal dominant optic atrophy and mitochondrial inherited Leber's hereditary optic neuropathy are the most common forms. Optic atrophy type 1, caused by mutations in the OPA1 gene is believed to be the most common hereditary optic neuropathy, and most patients inherit a mutation from an affected parent. In this study we used whole-exome sequencing to investigate the genetic aetiology in a patient affected with isolated optic atrophy. Since the proband was the only affected individual in his extended family, and was a product of consanguineous marriage, homozygosity mapping followed by whole-exome sequencing were pursued. Exome results identified a novel de novo OPA1 mutation in the proband. We conclude, that though de novo OPA1 mutations are uncommon, testing of common optic atrophy-associated genes such as mitochondrial mutations and OPA1 gene sequencing should be performed first in single individuals presenting with optic neuropathy, even when dominant inheritance is not apparent.

Objectives: Recent studies have highlighted the potential of analyses of genomic sharing to produce insight into the demographic processes affecting human populations. We study runs of homozygosity (ROH) in 18 Jewish populations, examining these groups in relation to 123 non-Jewish populations sampled worldwide., Methods: By sorting ROH into 3 length classes (short, intermediate, and long), we evaluate the impact of demographic processes on genomic patterns in Jewish populations., Results: We find that the portion of the genome appearing in long ROH - the length class most directly related to recent consanguinity - closely accords with data gathered from interviews during the 1950s on frequencies of consanguineous unions in various Jewish groups., Conclusion: The high correlation between 1950s consanguinity levels and coverage by long ROH explains differences across populations in ROH patterns. The dissection of ROH into length classes and the comparison to consanguinity data assist in understanding a number of additional phenomena, including similarities of Jewish populations to Middle Eastern, European, and Central and South Asian non-Jewish populations in short ROH patterns, relative lengths of identity-by-descent tracts in different Jewish groups, and the "population isolate" status of the Ashkenazi Jews., (© 2017 S. Karger AG, Basel.)

Genetic syndromes involving both brain and eye abnormalities are numerous and include syndromes such as Warburg micro syndrome, Kaufman oculocerebrofacial syndrome, Cerebro-oculo-facio-skeletal syndrome, Kahrizi syndrome and others. Using exome sequencing, we have been able to identify homozygous mutation p.(Tyr39Cys) in MED25 as the cause of a syndrome characterized by eye, brain, cardiac and palatal abnormalities as well as growth retardation, microcephaly and severe intellectual disability in seven patients from four unrelated families, all originating from the same village. The protein encoded by MED25 belongs to Mediator complex or MED complex, which is an evolutionary conserved multi-subunit RNA polymerase II transcriptional regulator complex. The MED25 point mutation is located in the von Willebrand factor type A (MED25 VWA) domain which is responsible for MED25 recruitment into the Mediator complex; co-immunoprecipitation experiment demonstrated that this mutation dramatically impairs MED25 interaction with the Mediator complex in mammalian cells.

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males., (© 2015 Karmin et al.; Published by Cold Spring Harbor Laboratory Press.)

Lissencephaly comprises a heterogeneous group of developmental brain disorders of varying severity, involving abnormal cortical gyration. We studied a highly consanguineous Israeli Moslem family with a lethal form of autosomal recessive lissencephaly with cerebellar hypoplasia (LCH). Using microarray-based homozygosity mapping in the reported family, combined with whole exome sequencing in one affected infant, we identified a homozygous splice site mutation g.IVS8+1G>A in cyclin-dependent kinase 5 (CDK5), causing complete skipping of exon 8, and leading to a frame shift and premature stop codon (p.V162SfsX19). The mutation co-segregated with the disease phenotype in all 29 study participants (4 patients and 25 healthy relatives), and was not identified in 200 ethnically matched control chromosomes. The p.V162SfsX19 mutation causes lack of endogenous CDK5 expression in affected dermal fibroblasts and brain tissue at the mRNA and protein levels, consistent with nonsense-mediated mRNA decay. Functional analysis of the p.V162SfsX19 mutation, using a yeast complementation assay, showed loss-of-function of the mutant CDK5 gene product, thereby implicating its role in the pathogenesis of autosomal recessive LCH in the studied family.

Congenital pancreatic lipase (PNLIP) deficiency is a rare monoenzymatic form of exocrine pancreatic failure characterized by decreased absorption of dietary fat and greasy voluminous stools, but apparent normal development and an overall good state of health. While considered to be an autosomal recessive state affecting a few dozens of individuals world-wide and involving the PNLIP gene, no causative mutations for this phenotype were so far reported. Here, we report the identification of the homozygote missense mutation, Thr221Met [c.662C>T], in two brothers from a consanguineous family of Arab ancestry. The observed genotypes among the family members were concordant with an autosomal recessive mode of inheritance but moreover a clear segregation between the genotype state and the serum PNLIP activity was evident. Based on biophysical computational tools, we suggest the mutation disrupts the protein's stability and impairs its normal function. Although the role of PNLIP is well established, our observations provide genetic evidence that PNLIP mutations are causative for this phenotype.