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.

Performance comparison: exome sequencing as a single test replacing Sanger sequencing.

Next generation sequencing tests are used routinely as first-choice tests in the clinic. However, systematic performance comparing the results of exome sequencing as a single test replacing Sanger sequencing of targeted gene(s) is still lacking. Performance comparison data are critically important for clinical case management. In this study, we compared Sanger-sequencing results of 258 genes to those obtained from next generation sequencing (NGS) using two exome-sequencing enrichment kits: Agilent-SureSelectQXT and Illumina-Nextera. Sequencing was performed on leukocytes and buccal-derived DNA from a single individual, and all 258 genes were sequenced a total of 11 times (using different sequencing methods and DNA sources). Sanger sequencing was completed for all exons, including flanking ± 8 bp regions. For the 258 genes, NGS mean coverage was > 20 × for > 98 and > 91% of the regions targeted by SureSelect and Nextera, respectively. Overall, 449 variants were identified in at least one experiment, and 407/449 (90.6%) were detected by all. Of the 42 discordant variants, 23 were determined as true calls, summing-up to a truth set of 430 variants. Sensitivity of true-variant detection was 99% for Sanger sequencing and 97-100% for the NGS experiments. Mean false-positive rates were 3.7E-6 for Sanger sequencing, 2.5E-6 for SureSelect-NGS and 5.2E-6 for Nextera-NGS. Our findings suggest a high overall concordance between Sanger sequencing and NGS performances. Both methods demonstrated false-positive and false-negative calls. High clinical suspicion for a specific diagnosis should, therefore, override negative results of either Sanger sequencing or NGS.

Ancestral patterns of recessive dystrophic epidermolysis bullosa mutations in Hispanic populations suggest sephardic ancestry.

Recessive dystrophic epidermolysis bullosa (RDEB) is a rare genodermatosis caused by mutations in the gene coding for type VII collagen (COL7A1). More than 800 different pathogenic mutations in COL7A1 have been described to date; however, the ancestral origins of many of these mutations have not been precisely identified. In this study, 32 RDEB patient samples from the Southwestern United States, Mexico, Chile, and Colombia carrying common mutations in the COL7A1 gene were investigated to determine the origins of these mutations and the extent to which shared ancestry contributes to disease prevalence. The results demonstrate both shared European and American origins of RDEB mutations in distinct populations in the Americas and suggest the influence of Sephardic ancestry in at least some RDEB mutations of European origins. Knowledge of ancestry and relatedness among RDEB patient populations will be crucial for the development of future clinical trials and the advancement of novel therapeutics.

Preconception carrier screening yield: effect of variants of unknown significance in partners of carriers with clinically significant variants.

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.

Teaching clinicians practical genomic medicine: 7 years' experience in a tertiary care center.

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.

The spectrum of BRCA1 and BRCA2 pathogenic sequence variants in Middle Eastern, North African, and South European countries.

BRCA1 BRCA2 mutational spectrum in the Middle East, North Africa, and Southern Europe is not well characterized. The unique history and cultural practices characterizing these regions, often involving consanguinity and inbreeding, plausibly led to the accumulation of population-specific founder pathogenic sequence variants (PSVs). To determine recurring BRCA PSVs in these locales, a search in PUBMED, EMBASE, BIC, and CIMBA was carried out combined with outreach to researchers from the relevant countries for unpublished data. We identified 232 PSVs in BRCA1 and 239 in BRCA2 in 25 of 33 countries surveyed. Common PSVs that were detected in four or more countries were c.5266dup (p.Gln1756Profs), c.181T>G (p.Cys61Gly), c.68_69del (p.Glu23Valfs), c.5030_5033del (p.Thr1677Ilefs), c.4327C>T (p.Arg1443Ter), c.5251C>T (p.Arg1751Ter), c.1016dup (p.Val340Glyfs), c.3700_3704del (p.Val1234Glnfs), c.4065_4068del (p.Asn1355Lysfs), c.1504_1508del (p.Leu502Alafs), c.843_846del (p.Ser282Tyrfs), c.798_799del (p.Ser267Lysfs), and c.3607C>T (p.Arg1203Ter) in BRCA1 and c.2808_2811del (p.Ala938Profs), c.5722_5723del (p.Leu1908Argfs), c.9097dup (p.Thr3033Asnfs), c.1310_1313del (p. p.Lys437Ilefs), and c.5946del (p.Ser1982Argfs) for BRCA2. Notably, some mutations (e.g., p.Asn257Lysfs (c.771_775del)) were observed in unrelated populations. Thus, seemingly genotyping recurring BRCA PSVs in specific populations may provide first pass BRCA genotyping platform.

Human Y Chromosome Haplogroup N: A Non-trivial Time-Resolved Phylogeography that Cuts across Language Families.

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.

Reconstructing the demographic history of the Himalayan and adjoining populations.

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.

Consanguinity Rates Predict Long Runs of Homozygosity in Jewish Populations.

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.

Exome sequencing identified mutations in CASK and MYBPC3 as the cause of a complex dilated cardiomyopathy phenotype.

Whole-exome sequencing for clinical applications is now an integral part of medical genetics practice. Though most studies are performed in order to establish diagnoses in individuals with rare and clinically unrecognizable disorders, due to the constantly decreasing costs and commercial availability, whole-exome sequencing has gradually become the initial tool to study patients with clinically recognized disorders when more than one gene is responsible for the phenotype or in complex phenotypes, when variants in more than one gene can be the cause for the disease. Here we report a patient presenting with a complex phenotype consisting of severe, adult-onset, dilated cardiomyopathy, hearing loss and developmental delay, in which exome sequencing revealed two genetic variants that are inherited from a healthy mother: a novel missense variant in the CASK gene, mutations in which cause a spectrum of neurocognitive manifestations, and a second variant, in MYBPC3, that is associated with hereditary cardiomyopathy. We conclude that although the potential for co-occurrence of rare diseases is higher when analyzing undefined phenotypes in consanguineous families, it should also be given consideration in the genetic evaluation of complex phenotypes in non-consanguineous families.

Exome sequencing identified a novel de novo OPA1 mutation in a consanguineous family presenting with optic atrophy.

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.

The genome-wide structure of the Jewish people.

Contemporary Jews comprise an aggregate of ethno-religious communities whose worldwide members identify with each other through various shared religious, historical and cultural traditions. Historical evidence suggests common origins in the Middle East, followed by migrations leading to the establishment of communities of Jews in Europe, Africa and Asia, in what is termed the Jewish Diaspora. This complex demographic history imposes special challenges in attempting to address the genetic structure of the Jewish people. Although many genetic studies have shed light on Jewish origins and on diseases prevalent among Jewish communities, including studies focusing on uniparentally and biparentally inherited markers, genome-wide patterns of variation across the vast geographic span of Jewish Diaspora communities and their respective neighbours have yet to be addressed. Here we use high-density bead arrays to genotype individuals from 14 Jewish Diaspora communities and compare these patterns of genome-wide diversity with those from 69 Old World non-Jewish populations, of which 25 have not previously been reported. These samples were carefully chosen to provide comprehensive comparisons between Jewish and non-Jewish populations in the Diaspora, as well as with non-Jewish populations from the Middle East and north Africa. Principal component and structure-like analyses identify previously unrecognized genetic substructure within the Middle East. Most Jewish samples form a remarkably tight subcluster that overlies Druze and Cypriot samples but not samples from other Levantine populations or paired Diaspora host populations. In contrast, Ethiopian Jews (Beta Israel) and Indian Jews (Bene Israel and Cochini) cluster with neighbouring autochthonous populations in Ethiopia and western India, respectively, despite a clear paternal link between the Bene Israel and the Levant. These results cast light on the variegated genetic architecture of the Middle East, and trace the origins of most Jewish Diaspora communities to the Levant.

Mutations in TAX1BP3 cause dilated cardiomyopathy with septo-optic dysplasia.

We describe a Bedouin family with a novel autosomal recessive syndrome characterized by dilated cardiomyopathy and septo-optic dysplasia. Genetic analysis revealed a homozygous missense mutation in TAX1BP3, which encodes a small PDZ domain containing protein implicated in regulation of the Wnt/ß-catenin signaling pathway, as the causative mutation. The mutation affects a conserved residue located at the core of TAX1BP3 binding pocket and is predicted to impair the nature of a crucial hydrophobic patch, thereby interrupting the structure and stability of the protein, and its ability to interact with other proteins. TAX1BP3 is highly expressed in heart and brain and consistent with the clinical findings observed in our patients; a knockdown of TAX1BP3 causes elongation defects, enlarged pericard, and enlarged head structures in zebrafish embryos. Thus, we describe a new genetic disorder that expands the monogenic cardiomyopathy disease spectrum and suggests that TAX1BP3 is essential for heart and brain development.

Autosomal recessive lissencephaly with cerebellar hypoplasia is associated with a loss-of-function mutation in CDK5.

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.

Homozygous MED25 mutation implicated in eye-intellectual disability syndrome.

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.

A "Copernican" reassessment of the human mitochondrial DNA tree from its root.

Mutational events along the human mtDNA phylogeny are traditionally identified relative to the revised Cambridge Reference Sequence, a contemporary European sequence published in 1981. This historical choice is a continuous source of inconsistencies, misinterpretations, and errors in medical, forensic, and population genetic studies. Here, after having refined the human mtDNA phylogeny to an unprecedented level by adding information from 8,216 modern mitogenomes, we propose switching the reference to a Reconstructed Sapiens Reference Sequence, which was identified by considering all available mitogenomes from Homo neanderthalensis. This "Copernican" reassessment of the human mtDNA tree from its deepest root should resolve previous problems and will have a substantial practical and educational influence on the scientific and public perception of human evolution by clarifying the core principles of common ancestry for extant descendants.

The Basque paradigm: genetic evidence of a maternal continuity in the Franco-Cantabrian region since pre-Neolithic times.

Different lines of evidence point to the resettlement of much of western and central Europe by populations from the Franco-Cantabrian region during the Late Glacial and Postglacial periods. In this context, the study of the genetic diversity of contemporary Basques, a population located at the epicenter of the Franco-Cantabrian region, is particularly useful because they speak a non-Indo-European language that is considered to be a linguistic isolate. In contrast with genome-wide analysis and Y chromosome data, where the problem of poor time estimates remains, a new timescale has been established for the human mtDNA and makes this genome the most informative marker for studying European prehistory. Here, we aim to increase knowledge of the origins of the Basque people and, more generally, of the role of the Franco-Cantabrian refuge in the postglacial repopulation of Europe. We thus characterize the maternal ancestry of 908 Basque and non-Basque individuals from the Basque Country and immediate adjacent regions and, by sequencing 420 complete mtDNA genomes, we focused on haplogroup H. We identified six mtDNA haplogroups, H1j1, H1t1, H2a5a1, H1av1, H3c2a, and H1e1a1, which are autochthonous to the Franco-Cantabrian region and, more specifically, to Basque-speaking populations. We detected signals of the expansion of these haplogroups at ∼4,000 years before present (YBP) and estimated their separation from the pan-European gene pool at ∼8,000 YBP, antedating the Indo-European arrival to the region. Our results clearly support the hypothesis of a partial genetic continuity of contemporary Basques with the preceding Paleolithic/Mesolithic settlers of their homeland.

Identification of a novel mutation in the PNLIP gene in two brothers with congenital pancreatic lipase deficiency.

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.

Human Y chromosome haplogroup L1-M22 traces Neolithic expansion in West Asia and supports the Elamite and Dravidian connection.

West and South Asian populations profoundly influenced Eurasian genetic and cultural diversity. We investigate the genetic history of the Y chromosome haplogroup L1-M22, which, while prevalent in these regions, lacks in-depth study. Robust Bayesian analyses of 165 high-coverage Y chromosomes favor a West Asian origin for L1-M22 ∼20.6 thousand years ago (kya). Moreover, this haplogroup parallels the genome-wide genetic ancestry of hunter-gatherers from the Iranian Plateau and the Caucasus. We characterized two L1-M22 harboring population groups during the Early Holocene. One expanded with the West Asian Neolithic transition. The other moved to South Asia ∼8-6 kya but showed no expansion. This group likely participated in the spread of Dravidian languages. These South Asian L1-M22 lineages expanded ∼4-3 kya, coinciding with the Steppe ancestry introduction. Our findings advance the current understanding of Eurasian historical dynamics, emphasizing L1-M22's West Asian origin, associated population movements, and possible linguistic impacts.

The Caucasus as an asymmetric semipermeable barrier to ancient human migrations.

The Caucasus, inhabited by modern humans since the Early Upper Paleolithic and known for its linguistic diversity, is considered to be important for understanding human dispersals and genetic diversity in Eurasia. We report a synthesis of autosomal, Y chromosome, and mitochondrial DNA (mtDNA) variation in populations from all major subregions and linguistic phyla of the area. Autosomal genome variation in the Caucasus reveals significant genetic uniformity among its ethnically and linguistically diverse populations and is consistent with predominantly Near/Middle Eastern origin of the Caucasians, with minor external impacts. In contrast to autosomal and mtDNA variation, signals of regional Y chromosome founder effects distinguish the eastern from western North Caucasians. Genetic discontinuity between the North Caucasus and the East European Plain contrasts with continuity through Anatolia and the Balkans, suggesting major routes of ancient gene flows and admixture.