Signatures of founder effects, admixture, and selection in the Ashkenazi Jewish population

Mitochondrial DNA Haplogroup K Is Protective Against Autism Spectrum Disorder Risk in Populations of European Ancestry

OBJECTIVE

Accumulative evidence indicates a critical role of mitochondrial function in autism spectrum disorders (ASD), implying that ASD risk may be linked to mitochondrial dysfunction due to DNA (mtDNA) variations. Although a few studies have explored the association between mtDNA variations and ASD, the role of mtDNA in ASD is still unclear. Here, we aimed to investigate whether mitochondrial DNA haplogroups are associated with the risk of ASD.

METHOD

Two European cohorts and an Ashkenazi Jewish (AJ) cohort were analyzed, including 2,062 ASD patients in comparison with 4,632 healthy controls. DNA samples were genotyped using Illumina HumanHap550/610 and Illumina 1M arrays, inclusive of mitochondrial markers. Mitochondrial DNA (mtDNA) haplogroups were identified from genotyping data using HaploGrep2. A mitochondrial genome imputation pipeline was established to detect mtDNA variants. We conducted a case-control study to investigate potential associations of mtDNA haplogroups and variants with the susceptibility to ASD.

RESULTS

We observed that the ancient adaptive mtDNA haplogroup K was significantly associated with decreased risk of ASD by the investigation of 2 European cohorts including a total of 2,006 cases and 4,435 controls (odds ratio = 0.64, P=1.79 × 10-5), and we replicated this association in an Ashkenazi Jewish (AJ) cohort including 56 cases and 197 controls (odds ratio = 0.35, P = 9.46 × 10-3). Moreover, we demonstrate that the mtDNA variants rs28358571, rs28358584, and rs28358280 are significantly associated with ASD risk. Further expression quantitative trait loci (eQTLs) analysis indicated that the rs28358584 and rs28358280 genotypes are associated with expression levels of nearby genes in brain tissues, suggesting those mtDNA variants may confer risk for ASD via regulation of expression levels of genes encoded by the mitochondrial genome.

CONCLUSION

This study helps to shed light on the contribution of mitochondria in ASD and provides new insights into the genetic mechanism underlying ASD, suggesting the potential involvement of mtDNA-encoded proteins in the development of ASD.

PLAIN LANGUAGE SUMMARY

Increasing evidence indicates that mitochondrial dysfunction may be linked to autism spectrum disorder (ASD). This study investigated potential associations of mitochondrial DNA (mtDNA) variants in 2 European and Ashkenazi Jewish cohorts including 2,062 individuals with ASD and 4,632 healthy controls. Researchers found that the ancient mtDNA haplogroup K was linked to a reduced risk of ASD in both European and Ashkenazi Jewish populations. Additionally, specific mtDNA variants were associated with ASD risk and were shown to influence the expression of nearby genes in the brain. These findings highlight the potential involvement of mtDNA in ASD development, offering new insights into the genetic mechanisms underlying the disorder.

GenoTools: An Open-Source Python Package for Efficient Genotype Data Quality Control and Analysis

GenoTools, a Python package, streamlines population genetics research by integrating ancestry estimation, quality control (QC), and genome-wide association studies (GWAS) capabilities into efficient pipelines. By tracking samples, variants, and quality-specific measures throughout fully customizable pipelines, users can easily manage genetics data for large and small studies. GenoTools' "Ancestry" module renders highly accurate predictions, allowing for high-quality ancestry-specific studies, and enables custom ancestry model training and serialization, specified to the user's genotyping or sequencing platform. As the genotype processing engine that powers several large initiatives, including the NIH's Center for Alzheimer's and Related Dementias (CARD) and the Global Parkinson's Genetics Program (GP2). GenoTools was used to process and analyze the UK Biobank and major Alzheimer's Disease (AD) and Parkinson's Disease (PD) datasets with over 400,000 genotypes from arrays and 5000 sequences and has led to novel discoveries in diverse populations. It has provided replicable ancestry predictions, implemented rigorous QC, and conducted genetic ancestry-specific GWAS to identify systematic errors or biases through a single command. GenoTools is a customizable tool that enables users to efficiently analyze and scale genotype data with reproducible and scalable ancestry, QC, and GWAS pipelines.

Insights into Ancestral Diversity in Parkinsons Disease Risk: A Comparative Assessment of Polygenic Risk Scores

Objectives To evaluate and compare different polygenic risk score (PRS) models in predicting Parkinsons disease (PD) across diverse ancestries, focusing on identifying the most suitable approach for each population and potentially contributing to equitable advancements in precision medicine. Methods We constructed a total of 105 PRS across individual level data from seven diverse ancestries. First, a cross-ancestry conventional PRS comparison was implemented by utilizing the 90 known European risk loci with weighted effects from four independent summary statistics including European, East Asian, Latino/Admixed American, and African/Admixed. These models were adjusted by sex, age, and principal components (28 PRS) and by sex, age, and percentage of admixture (28 PRS) for comparison. Secondly, a novel and refined multi-ancestry best-fit PRS approach was then applied across the seven ancestries by leveraging multi-ancestry meta-analyzed summary statistics and using a p-value thresholding approach (49 PRS) to enhance prediction applicability in a global setting. Results European-based PRS models predicted disease status across all ancestries to differing degrees of accuracy. Ashkenazi Jewish had the highest Odds Ratio (OR): 1.96 (95% CI: 1.69-2.25, p < 0.0001) with an AUC (Area Under the Curve) of 68%. Conversely, the East Asian population, despite having fewer predictive variants (84 out of 90), had an OR of 1.37 (95% CI: 1.32-1.42) and an AUC of 62%, illustrating the cross-ancestry transferability of this model. Lower OR alongside broader confidence intervals were observed in other populations, including Africans (OR =1.38, 95% CI: 1.12-1.63, p=0.001). Adjustment by percentage of admixture did not outperform principal components. Multi-ancestry best-fit PRS models improved risk prediction in European, Ashkenazi Jewish, and African ancestries, yet didn't surpass conventional PRS in admixed populations such as Latino/American admixed and African admixed populations. Interpretation The present study represents a novel and comprehensive assessment of PRS performance across seven ancestries in PD, highlighting the inadequacy of a 'one size fits all' approach in genetic risk prediction. We demonstrated that European based PD PRS models are partially transferable to other ancestries and could be improved by a novel best-fit multi-ancestry PRS, especially in non-admixed populations.

Federated learning for multi-omics: A performance evaluation in Parkinson's disease

While machine learning (ML) research has recently grown more in popularity, its application in the omics domain is constrained by access to sufficiently large, high-quality datasets needed to train ML models. Federated learning (FL) represents an opportunity to enable collaborative curation of such datasets among participating institutions. We compare the simulated performance of several models trained using FL against classically trained ML models on the task of multi-omics Parkinson's disease prediction. We find that FL model performance tracks centrally trained ML models, where the most performant FL model achieves an AUC-PR of 0.876 ± 0.009, 0.014 ± 0.003 less than its centrally trained variation. We also determine that the dispersion of samples within a federation plays a meaningful role in model performance. Our study implements several open-source FL frameworks and aims to highlight some of the challenges and opportunities when applying these collaborative methods in multi-omics studies.

Federated Learning for multi-omics: a performance evaluation in Parkinson's disease

While machine learning (ML) research has recently grown more in popularity, its application in the omics domain is constrained by access to sufficiently large, high-quality datasets needed to train ML models. Federated Learning (FL) represents an opportunity to enable collaborative curation of such datasets among participating institutions. We compare the simulated performance of several models trained using FL against classically trained ML models on the task of multi-omics Parkinson's Disease prediction. We find that FL model performance tracks centrally trained ML models, where the most performant FL model achieves an AUC-PR of 0.876 ± 0.009, 0.014 ± 0.003 less than its centrally trained variation. We also determine that the dispersion of samples within a federation plays a meaningful role in model performance. Our study implements several open source FL frameworks and aims to highlight some of the challenges and opportunities when applying these collaborative methods in multi-omics studies.

Genetic risk factor clustering within and across neurodegenerative diseases

Overlapping symptoms and co-pathologies are common in closely related neurodegenerative diseases (NDDs). Investigating genetic risk variants across these NDDs can give further insight into disease manifestations. In this study we have leveraged genome-wide single nucleotide polymorphisms and genome-wide association study summary statistics to cluster patients based on their genetic status across identified risk variants for five NDDs (Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, Lewy body dementia and frontotemporal dementia). The multi-disease and disease-specific clustering results presented here provide evidence that NDDs have more overlapping genetic aetiology than previously expected and how neurodegeneration should be viewed as a spectrum of symptomology. These clustering analyses also show potential subsets of patients with these diseases that are significantly depleted for any known common genetic risk factors suggesting environmental or other factors at work. Establishing that NDDs with overlapping pathologies share genetic risk loci, future research into how these variants might have different effects on downstream protein expression, pathology and NDD manifestation in general is important for refining and treating NDDs.

Chromosome-level genome assembly and whole-genome resequencing of topmouth culter (Culter alburnus) provide insights into the intraspecific variation of its semi-buoyant and adhesive eggs

Topmouth culter (Culter alburnus) is an ecologically and economically important species belonging to the subfamily Culterinae that is native to and widespread in East Asia. Intraspecific variation of semi-buoyant and adhesive eggs in topmouth culter provides an ideal opportunity to investigate the genetic mechanisms of spawning habits underlying the adaptive radiation of cyprinids in East Asia. In this study, we present a chromosome-level genome assembly of topmouth culter and re-sequenced 158 individuals from six locations in China covering three geographical groups and two egg type variations. The topmouth culter genome size was 1.05 Gb, with a contig N50 length of 17.8 Mb and anchored onto 24 chromosomes. Phylogenetic analysis showed that the divergence time of the Culterinae was coinciding with the time of initiation of the Asian monsoon intensification. Gene family evolutionary analysis indicated that the expanded gene families in topmouth culter were associated with dietary adaptation. Population-level genetic analysis indicated clear differentiation among the six populations, which were clustered into three distinct clusters, consistent with their geographical divergence. The historical effective population size of topmouth culter correlated with the Tibetan Plateau uplifting according to the demographic history reconstruction. A selective sweep analysis between adhesive and semi-buoyant egg populations revealed the genes associated with the hydration and adhesiveness of eggs, indicating divergent selection towards different hydrological environments. This study offers a high-resolution genetic resource for further studies on evolutionary adaptation, genetic breeding and conservation of topmouth culter, providing insights into the molecular mechanisms for egg type variation of East Asian cyprinids.

Genome-wide Association Identifies Novel Etiological Insights Associated with Parkinson's Disease in African and African Admixed Populations

Background

Understanding the genetic mechanisms underlying diseases in ancestrally diverse populations is a critical step towards the realization of the global application of precision medicine. The African and African admixed populations enable mapping of complex traits given their greater levels of genetic diversity, extensive population substructure, and distinct linkage disequilibrium patterns.

Methods

Here we perform a comprehensive genome-wide assessment of Parkinson's disease (PD) in 197,918 individuals (1,488 cases; 196,430 controls) of African and African admixed ancestry, characterizing population-specific risk, differential haplotype structure and admixture, coding and structural genetic variation and polygenic risk profiling.

Findings

We identified a novel common risk factor for PD and age at onset at the GBA1 locus (risk, rs3115534-G; OR=1.58, 95% CI = 1.37 - 1.80, P=2.397E-14; age at onset, BETA =-2.004, SE =0.57, P = 0.0005), that was found to be rare in non-African/African admixed populations. Downstream short- and long-read whole genome sequencing analyses did not reveal any coding or structural variant underlying the GWAS signal. However, we identified that this signal mediates PD risk via expression quantitative trait locus (eQTL) mechanisms. While previously identified GBA1 associated disease risk variants are coding mutations, here we suggest a novel functional mechanism consistent with a trend in decreasing glucocerebrosidase activity levels. Given the high population frequency of the underlying signal and the phenotypic characteristics of the homozygous carriers, we hypothesize that this variant may not cause Gaucher disease. Additionally, the prevalence of Gaucher's disease in Africa is low.

Interpretation

The present study identifies a novel African-ancestry genetic risk factor in GBA1 as a major mechanistic basis of PD in the African and African admixed populations. This striking result contrasts to previous work in Northern European populations, both in terms of mechanism and attributable risk. This finding highlights the importance of understanding population-specific genetic risk in complex diseases, a particularly crucial point as the field moves toward precision medicine in PD clinical trials and while recognizing the need for equitable inclusion of ancestrally diverse groups in such trials. Given the distinctive genetics of these underrepresented populations, their inclusion represents a valuable step towards insights into novel genetic determinants underlying PD etiology. This opens new avenues towards RNA-based and other therapeutic strategies aimed at reducing lifetime risk.

Research in Context

Evidence Before this Study

Our current understanding of Parkinson's disease (PD) is disproportionately based on studying populations of European ancestry, leading to a significant gap in our knowledge about the genetics, clinical characteristics, and pathophysiology in underrepresented populations. This is particularly notable in individuals of African and African admixed ancestries. Over the last two decades, we have witnessed a revolution in the research area of complex genetic diseases. In the PD field, large-scale genome-wide association studies in the European, Asian, and Latin American populations have identified multiple risk loci associated with disease. These include 78 loci and 90 independent signals associated with PD risk in the European population, nine replicated loci and two novel population-specific signals in the Asian population, and a total of 11 novel loci recently nominated through multi-ancestry GWAS efforts.Nevertheless, the African and African admixed populations remain completely unexplored in the context of PD genetics.

Added Value of this Study

To address the lack of diversity in our research field, this study aimed to conduct the first genome-wide assessment of PD genetics in the African and African admixed populations. Here, we identified a genetic risk factor linked to PD etiology, dissected African-specific differences in risk and age at onset, characterized known genetic risk factors, and highlighted the utility of the African and African admixed risk haplotype substructure for future fine-mapping efforts. We identified a novel disease mechanism via expression changes consistent with decreased GBA1 activity levels. Future large scale single cell expression studies should investigate the neuronal populations in which expression differences are most prominent. This novel mechanism may hold promise for future efficient RNA-based therapeutic strategies such as antisense oligonucleotides or short interfering RNAs aimed at preventing and decreasing disease risk. We envisage that these data generated under the umbrella of the Global Parkinson's Genetics Program (GP2) will shed light on the molecular mechanisms involved in the disease process and might pave the way for future clinical trials and therapeutic interventions. This work represents a valuable resource in an underserved population, supporting pioneering research within GP2 and beyond. Deciphering causal and genetic risk factors in all these ancestries will help determine whether interventions, potential targets for disease modifying treatment, and prevention strategies that are being studied in the European populations are relevant to the African and African admixed populations.

Implications of all the Available Evidence

We nominate a novel signal impacting GBA1 as the major genetic risk factor for PD in the African and African admixed populations. The present study could inform future GBA1 clinical trials, improving patient stratification. In this regard, genetic testing can help to design trials likely to provide meaningful and actionable answers. It is our hope that these findings may ultimately have clinical utility for this underrepresented population.

Effects of local domestication warrant attention in honey bee population genetics

Honey bees, Apis mellifera, have for millennia been managed and exploited by humans and introduced into most suitable regions worldwide. However, given the lack of records for many introduction events, treating A. mellifera populations as native would predictably bias genetic studies regarding origin and evolution. Here, we used the Dongbei bee, a well-documented population, introduced beyond the natural distribution range approximately 100 years ago, to elucidate the effects of local domestication on animal population genetic analyses. Strong domestication pressure was detected in this population, and the genetic divergence between Dongbei bee and its ancestral subspecies was found to have occurred at the lineage level. Results of phylogenetic and time divergence analyses could consequently be misinterpreted. Proposing new subspecies or lineages and performing analyses of origin should thus strive to eliminate anthropogenic effects. We highlight the need for definitions of landrace and breed in honey bee sciences and make preliminary suggestions.

Polygenic Parkinson's Disease Genetic Risk Score as Risk Modifier of Parkinsonism in Gaucher Disease

BACKGROUND

Biallelic pathogenic variants in GBA1 are the cause of Gaucher disease (GD) type 1 (GD1), a lysosomal storage disorder resulting from deficient glucocerebrosidase. Heterozygous GBA1 variants are also a common genetic risk factor for Parkinson's disease (PD). GD manifests with considerable clinical heterogeneity and is also associated with an increased risk for PD.

OBJECTIVE

The objective of this study was to investigate the contribution of PD risk variants to risk for PD in patients with GD1.

METHODS

We studied 225 patients with GD1, including 199 without PD and 26 with PD. All cases were genotyped, and the genetic data were imputed using common pipelines.

RESULTS

On average, patients with GD1 with PD have a significantly higher PD genetic risk score than those without PD (P = 0.021).

CONCLUSIONS

Our results indicate that variants included in the PD genetic risk score were more frequent in patients with GD1 who developed PD, suggesting that common risk variants may affect underlying biological pathways. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Gaucher disease protects against tuberculosis

Biallelic mutations in the glucocerebrosidase (GBA1) gene cause Gaucher disease, characterized by lysosomal accumulation of glucosylceramide and glucosylsphingosine in macrophages. Gaucher and other lysosomal diseases occur with high frequency in Ashkenazi Jews. It has been proposed that the underlying mutations confer a selective advantage, in particular conferring protection against tuberculosis. Here, using a zebrafish Gaucher disease model, we find that the mutation GBA1 N370S, predominant among Ashkenazi Jews, increases resistance to tuberculosis through the microbicidal activity of glucosylsphingosine in macrophage lysosomes. Consistent with lysosomal accumulation occurring only in homozygotes, heterozygotes remain susceptible to tuberculosis. Thus, our findings reveal a mechanistic basis for protection against tuberculosis by GBA1 N370S and provide biological plausibility for its selection if the relatively mild deleterious effects in homozygotes were offset by significant protection against tuberculosis, a rampant killer of the young in Europe through the Middle Ages into the 19th century.

Clustering of Juvenile Canavan disease in an Indian community due to population bottleneck and isolation: genomic signatures of a founder event

Mild/juvenile Canavan disease (M/JCD) is less frequently reported in the literature and little is known about its pathogenetic mechanisms. We report a comprehensive investigation into the pathogenetic mechanism of a novel NM_000049.4(ASPA):c.526G>A variant in two families. The families belong to Telugu Devanga Chettiar community (TDC) from southern India. TDC has a complex history of migration from their historical origin centuries ago with high endogamy. TDC probably has the highest clustering M/JCD recorded historically (around 24 cases). The pathogenic variant was shown to cause non-classical splicing defect resulting in two different transcripts. The splicing aberration, a loss of function mechanism coupled with a milder missense effect can explain the milder phenotype compared to the infantile-onset CD. The high clustering of an extremely rare form of neurodegenerative disorder with reduced fitness, led us to speculate the possibility of a founder event. Genotyping array of TDC and multiple distinct populations of Indian origin for several population genetic parameters was performed. It yielded robust signatures of a founder event in TDC, such as a high fixation index, increased runs of homozygosity and identity-by-descent in the absence of consanguinity; a large haplotype with high linkage disequilibrium among markers comprising the pathogenic variant; a robust population structure; mutation dating, estimating the age of the potential founder of TDC at around 375 years; possibly a high carrier rate in TDC. This study has not only focused its attention on natural history and pathogenetics but also paves way for carrier screening programs in TDC and future therapeutic studies.

Genome-wide data from medieval German Jews show that the Ashkenazi founder event pre-dated the 14th century

We report genome-wide data from 33 Ashkenazi Jews (AJ), dated to the 14^th century, obtained following a salvage excavation at the medieval Jewish cemetery of Erfurt, Germany. The Erfurt individuals are genetically similar to modern AJ, but they show more variability in Eastern European-related ancestry than modern AJ. A third of the Erfurt individuals carried a mitochondrial lineage common in modern AJ and eight carried pathogenic variants known to affect AJ today. These observations, together with high levels of runs of homozygosity, suggest that the Erfurt community had already experienced the major reduction in size that affected modern AJ. The Erfurt bottleneck was more severe, implying substructure in medieval AJ. Overall, our results suggest that the AJ founder event and the acquisition of the main sources of ancestry pre-dated the 14^th century and highlight late medieval genetic heterogeneity no longer present in modern AJ.

Genomes from a medieval mass burial show Ashkenazi-associated hereditary diseases pre-date the 12th century

We report genome sequence data from six individuals excavated from the base of a medieval well at a site in Norwich, UK. A revised radiocarbon analysis of the assemblage is consistent with these individuals being part of a historically attested episode of antisemitic violence on 6 February 1190 CE. We find that four of these individuals were closely related and all six have strong genetic affinities with modern Ashkenazi Jews. We identify four alleles associated with genetic disease in Ashkenazi Jewish populations and infer variation in pigmentation traits, including the presence of red hair. Simulations indicate that Ashkenazi-associated genetic disease alleles were already at appreciable frequencies, centuries earlier than previously hypothesized. These findings provide new insights into a significant historical crime, into Ashkenazi population history, and into the origins of genetic diseases associated with modern Jewish populations.

Principal Component Analyses (PCA)-based findings in population genetic studies are highly biased and must be reevaluated

Principal Component Analysis (PCA) is a multivariate analysis that reduces the complexity of datasets while preserving data covariance. The outcome can be visualized on colorful scatterplots, ideally with only a minimal loss of information. PCA applications, implemented in well-cited packages like EIGENSOFT and PLINK, are extensively used as the foremost analyses in population genetics and related fields (e.g., animal and plant or medical genetics). PCA outcomes are used to shape study design, identify, and characterize individuals and populations, and draw historical and ethnobiological conclusions on origins, evolution, dispersion, and relatedness. The replicability crisis in science has prompted us to evaluate whether PCA results are reliable, robust, and replicable. We analyzed twelve common test cases using an intuitive color-based model alongside human population data. We demonstrate that PCA results can be artifacts of the data and can be easily manipulated to generate desired outcomes. PCA adjustment also yielded unfavorable outcomes in association studies. PCA results may not be reliable, robust, or replicable as the field assumes. Our findings raise concerns about the validity of results reported in the population genetics literature and related fields that place a disproportionate reliance upon PCA outcomes and the insights derived from them. We conclude that PCA may have a biasing role in genetic investigations and that 32,000-216,000 genetic studies should be reevaluated. An alternative mixed-admixture population genetic model is discussed.

Reconstructing the history of founder events using genome-wide patterns of allele sharing across individuals

Founder events play a critical role in shaping genetic diversity, fitness and disease risk in a population. Yet our understanding of the prevalence and distribution of founder events in humans and other species remains incomplete, as most existing methods require large sample sizes or phased genomes. Thus, we developed ASCEND that measures the correlation in allele sharing between pairs of individuals across the genome to infer the age and strength of founder events. We show that ASCEND can reliably estimate the parameters of founder events under a range of demographic scenarios. We then apply ASCEND to two species with contrasting evolutionary histories: ~460 worldwide human populations and ~40 modern dog breeds. In humans, we find that over half of the analyzed populations have evidence for recent founder events, associated with geographic isolation, modes of sustenance, or cultural practices such as endogamy. Notably, island populations have lower population sizes than continental groups and most hunter-gatherer, nomadic and indigenous groups have evidence of recent founder events. Many present-day groups––including Native Americans, Oceanians and South Asians––have experienced more extreme founder events than Ashkenazi Jews who have high rates of recessive diseases due their known history of founder events. Using ancient genomes, we show that the strength of founder events differs markedly across geographic regions and time––with three major founder events related to the peopling of Americas and a trend in decreasing strength of founder events in Europe following the Neolithic transition and steppe migrations. In dogs, we estimate extreme founder events in most breeds that occurred in the last 25 generations, concordant with the establishment of many dog breeds during the Victorian times. Our analysis highlights a widespread history of founder events in humans and dogs and elucidates some of the demographic and cultural practices related to these events.

A founder event occurs when small numbers of ancestral individuals give rise to a large fraction of the population. Founder events reduce genetic variation and increase the risk of recessive diseases. Despite their importance in evolutionary and disease studies, we still only have a limited comprehension of their prevalence and properties in humans and other species, as most existing methods require large sample sizes or phased genomes. Here, we present a flexible method, ASCEND, to infer the timing and the strength of founder events that is suitable for sparse datasets with few samples or limited coverage. ASCEND provides reliable estimates across a wide range of demographic scenarios. By applying it to data from two species (humans and dogs), we document a widespread history of recent founder events in both species and provide insights about the demographic processes related to these events. Our analysis helps to identify groups with strong founder events that should be prioritized for future studies as they offer a unique opportunity for biological discovery and reducing disease burden through mapping of recessive disease-causing genes and pathways, as previously shown in studies of Ashkenazi Jews and Finns.

Multi-modality machine learning predicting Parkinson's disease

Personalized medicine promises individualized disease prediction and treatment. The convergence of machine learning (ML) and available multimodal data is key moving forward. We build upon previous work to deliver multimodal predictions of Parkinson's disease (PD) risk and systematically develop a model using GenoML, an automated ML package, to make improved multi-omic predictions of PD, validated in an external cohort. We investigated top features, constructed hypothesis-free disease-relevant networks, and investigated drug-gene interactions. We performed automated ML on multimodal data from the Parkinson's progression marker initiative (PPMI). After selecting the best performing algorithm, all PPMI data was used to tune the selected model. The model was validated in the Parkinson's Disease Biomarker Program (PDBP) dataset. Our initial model showed an area under the curve (AUC) of 89.72% for the diagnosis of PD. The tuned model was then tested for validation on external data (PDBP, AUC 85.03%). Optimizing thresholds for classification increased the diagnosis prediction accuracy and other metrics. Finally, networks were built to identify gene communities specific to PD. Combining data modalities outperforms the single biomarker paradigm. UPSIT and PRS contributed most to the predictive power of the model, but the accuracy of these are supplemented by many smaller effect transcripts and risk SNPs. Our model is best suited to identifying large groups of individuals to monitor within a health registry or biobank to prioritize for further testing. This approach allows complex predictive models to be reproducible and accessible to the community, with the package, code, and results publicly available.

Parallel subgenome structure and divergent expression evolution of allo-tetraploid common carp and goldfish

How two subgenomes in allo-tetraploids adapt to coexistence and coordinate through structure and expression evolution requires extensive studies. In the present study, we report an improved genome assembly of allo-tetraploid common carp, an updated genome annotation of allo-tetraploid goldfish and the chromosome-scale assemblies of a progenitor-like diploid Puntius tetrazona and an outgroup diploid Paracanthobrama guichenoti. Parallel subgenome structure evolution in the allo-tetraploids was featured with equivalent chromosome components, higher protein identities, similar transposon divergence and contents, homoeologous exchanges, better synteny level, strong sequence compensation and symmetric purifying selection. Furthermore, we observed subgenome expression divergence processes in the allo-tetraploids, including inter-/intrasubgenome trans-splicing events, expression dominance, decreased expression levels, dosage compensation, stronger expression correlation, dynamic functionalization and balancing of differential expression. The potential disorders introduced by different progenitors in the allo-tetraploids were hypothesized to be alleviated by increasing structural homogeneity and performing versatile expression processes. Resequencing three common carp strains revealed two major ecotypes and uncovered candidate genes relevant to growth and survival rate.

Accelerating Medicines Partnership: Parkinson's Disease. Genetic Resource

BACKGROUND

Whole-genome sequencing data are available from several large studies across a variety of diseases and traits. However, massive storage and computation resources are required to use these data, and to achieve sufficient power for discoveries, harmonization of multiple cohorts is critical.

OBJECTIVES

The Accelerating Medicines Partnership Parkinson's Disease program has developed a research platform for Parkinson's disease (PD) that integrates the storage and analysis of whole-genome sequencing data, RNA expression data, and clinical data, harmonized across multiple cohort studies.

METHODS

The version 1 release contains whole-genome sequencing data derived from 3941 participants from 4 cohorts. Samples underwent joint genotyping by the TOPMed Freeze 9 Variant Calling Pipeline. We performed descriptive analyses of these whole-genome sequencing data using the Accelerating Medicines Partnership Parkinson's Disease platform.

RESULTS

The clinical diagnosis of participants in version 1 release includes 2005 idiopathic PD patients, 963 healthy controls, 64 prodromal subjects, 62 clinically diagnosed PD subjects without evidence of dopamine deficit, and 705 participants of genetically enriched cohorts carrying PD risk-associated GBA variants or LRRK2 variants, of whom 304 were affected. We did not observe significant enrichment of pathogenic variants in the idiopathic PD group, but the polygenic risk score was higher in PD both in nongenetically enriched cohorts and genetically enriched cohorts. The population analysis showed a correlation between genetically enriched cohorts and Ashkenazi Jewish ancestry.

CONCLUSIONS

We describe the genetic component of the Accelerating Medicines Partnership Parkinson's Disease platform, a solution to democratize data access and analysis for the PD research community. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article is a U.S. Government work and is in the public domain in the USA.

The Genetic Analyses of French Canadians of Quebec Facilitate the Characterization of New Cancer Predisposing Genes Implicated in Hereditary Breast and/or Ovarian Cancer Syndrome Families

The French Canadian population of the province of Quebec has been recognized for its contribution to research in medical genetics, especially in defining the role of heritable pathogenic variants in cancer predisposing genes. Multiple carriers of a limited number of pathogenic variants in BRCA1 and BRCA2, the major risk genes for hereditary breast and/or ovarian cancer syndrome families, have been identified in French Canadians, which is in stark contrast to the array of over 2000 different pathogenic variants reported in each of these genes in other populations. As not all such cancer syndrome families are explained by BRCA1 and BRCA2, newly proposed gene candidates identified in other populations have been investigated for their role in conferring risk in French Canadian cancer families. For example, multiple carriers of distinct variants were identified in PALB2 and RAD51D. The unique genetic architecture of French Canadians has been attributed to shared ancestry due to common ancestors of early settlers of this population with origins mainly from France. In this review, we discuss the merits of genetically characterizing cancer predisposing genes in French Canadians of Quebec. We focused on genes that have been implicated in hereditary breast and/or ovarian cancer syndrome families as they have been the most thoroughly characterized cancer syndromes in this population. We describe how genetic analyses of French Canadians have facilitated: (i) the classification of variants in BRCA1 and BRCA2; (ii) the identification and classification of variants in newly proposed breast and/or ovarian cancer predisposing genes; and (iii) the identification of a new breast cancer predisposing gene candidate, RECQL. The genetic architecture of French Canadians provides a unique opportunity to evaluate new candidate cancer predisposing genes regardless of the population in which they were identified.

Genome-wide characterization of human minisatellite VNTRs: population-specific alleles and gene expression differences

Variable Number Tandem Repeats (VNTRs) are tandem repeat (TR) loci that vary in copy number across a population. Using our program, VNTRseek, we analyzed human whole genome sequencing datasets from 2770 individuals in order to detect minisatellite VNTRs, i.e., those with pattern sizes ≥7 bp. We detected 35 638 VNTR loci and classified 5676 as commonly polymorphic (i.e. with non-reference alleles occurring in >5% of the population). Commonly polymorphic VNTR loci were found to be enriched in genomic regions with regulatory function, i.e. transcription start sites and enhancers. Investigation of the commonly polymorphic VNTRs in the context of population ancestry revealed that 1096 loci contained population-specific alleles and that those could be used to classify individuals into super-populations with near-perfect accuracy. Search for quantitative trait loci (eQTLs), among the VNTRs proximal to genes, indicated that in 187 genes expression differences correlated with VNTR genotype. We validated our predictions in several ways, including experimentally, through the identification of predicted alleles in long reads, and by comparisons showing consistency between sequencing platforms. This study is the most comprehensive analysis of minisatellite VNTRs in the human population to date.

Drivers of linkage disequilibrium across a species' geographic range

While linkage disequilibrium (LD) is an important parameter in genetics and evolutionary biology, the drivers of LD remain elusive. Using whole-genome sequences from across a species’ range, we assessed the impact of demographic history and mating system on LD. Both range expansion and a shift from outcrossing to selfing in North American Arabidopsis lyrata were associated with increased average genome-wide LD. Our results indicate that range expansion increases short-distance LD at the farthest range edges by about the same amount as a shift to selfing. However, the extent over which LD in genic regions unfolds was shorter for range expansion compared to selfing. Linkage among putatively neutral variants and between neutral and deleterious variants increased to a similar degree with range expansion, providing support that genome-wide LD was positively associated with mutational load. As a consequence, LD combined with mutational load may decelerate range expansions and set range limits. Finally, a small number of genes were identified as LD outliers, suggesting that they experience selection by either of the two demographic processes. These included genes involved in flowering and photoperiod for range expansion, and the self-incompatibility locus for mating system.

Nearby genomic variants are often co-inherited because of limited recombination. The extent of non-random association of alleles at different loci is called linkage disequilibrium (LD) and is commonly used in genomic analyses, for example to detect regions under selection or to determine effective population size. Here we reversed testing and addressed how demographic history may affect LD within a species. Using genomic data from more than a thousand individuals of North American Arabidopsis lyrata from across the entire species’ range, we quantified the effect of postglacial range expansion and a shift in mating system from outcrossing to selfing on LD. We show that both factors lead to increased LD, and that the maximal effect of range expansion is comparable with a shift in mating system to selfing. Heightened LD involves deleterious mutations, and therefore, LD can also serve as an indicator of mutation accumulation. Furthermore, we provide evidence that some genes experienced stronger increases in LD possibly due to selection associated with the two demographic changes. Our results provide a novel and broad view on the evolutionary factors shaping LD that may also apply to the very many species that underwent postglacial range expansion.

Biallelic loss-of-function variants in PLD1 cause congenital right-sided cardiac valve defects and neonatal cardiomyopathy

Congenital heart disease is the most common type of birth defect, accounting for one-third of all congenital anomalies. Using whole-exome sequencing of 2718 patients with congenital heart disease and a search in GeneMatcher, we identified 30 patients from 21 unrelated families of different ancestries with biallelic phospholipase D1 (PLD1) variants who presented predominantly with congenital cardiac valve defects. We also associated recessive PLD1 variants with isolated neonatal cardiomyopathy. Furthermore, we established that p.I668F is a founder variant among Ashkenazi Jews (allele frequency of ~2%) and describe the phenotypic spectrum of PLD1-associated congenital heart defects. PLD1 missense variants were overrepresented in regions of the protein critical for catalytic activity, and, correspondingly, we observed a strong reduction in enzymatic activity for most of the mutant proteins in an enzymatic assay. Finally, we demonstrate that PLD1 inhibition decreased endothelial-mesenchymal transition, an established pivotal early step in valvulogenesis. In conclusion, our study provides a more detailed understanding of disease mechanisms and phenotypic expression associated with PLD1 loss of function.

Pancreatic cancer pathology viewed in the light of evolution

One way to understand ductal adenocarcinoma of the pancreas (pancreatic cancer) is to view it as unimaginably large numbers of evolving living organisms interacting with their environment. This “evolutionary view” creates both expected and surprising perspectives in all stages of neoplastic progression. Advances in the field will require greater attention to this critical evolutionary prospective.

Population Histories of the United States Revealed through Fine-Scale Migration and Haplotype Analysis

The population of the United States is shaped by centuries of migration, isolation, growth, and admixture between ancestors of global origins. Here, we assemble a comprehensive view of recent population history by studying the ancestry and population structure of more than 32,000 individuals in the US using genetic, ancestral birth origin, and geographic data from the National Geographic Genographic Project. We identify migration routes and barriers that reflect historical demographic events. We also uncover the spatial patterns of relatedness in subpopulations through the combination of haplotype clustering, ancestral birth origin analysis, and local ancestry inference. Examples of these patterns include substantial substructure and heterogeneity in Hispanics/Latinos, isolation-by-distance in African Americans, elevated levels of relatedness and homozygosity in Asian immigrants, and fine-scale structure in European descents. Taken together, our results provide detailed insights into the genetic structure and demographic history of the diverse US population.

High-resolution inference of genetic relationships among Jewish populations

Recent studies have used genome-wide single-nucleotide polymorphisms (SNPs) to investigate relationships among various Jewish populations and their non-Jewish historical neighbors, often focusing on small subsets of populations from a limited geographic range or relatively small samples within populations. Here, building on the significant progress that has emerged from genomic SNP studies in the placement of Jewish populations in relation to non-Jewish populations, we focus on population structure among Jewish populations. In particular, we examine Jewish population-genetic structure in samples that span much of the historical range of Jewish populations in Europe, the Middle East, North Africa, and South Asia. Combining 429 newly genotyped samples from 29 Jewish and 3 non-Jewish populations with previously reported genotypes on Jewish and non-Jewish populations, we investigate variation in 2789 individuals from 114 populations at 486,592 genome-wide autosomal SNPs. Using multidimensional scaling analysis, unsupervised model-based clustering, and population trees, we find that, genetically, most Jewish samples fall into four major clusters that largely represent four culturally defined groupings, namely the Ashkenazi, Mizrahi, North African, and Sephardi subdivisions of the Jewish population. We detect high-resolution population structure, including separation of the Ashkenazi and Sephardi groups and distinctions among populations within the Mizrahi and North African groups. Our results refine knowledge of Jewish population-genetic structure and contribute to a growing understanding of the distinctive genetic ancestry evident in closely related but historically separate Jewish communities.

Disparities in discovery of pathogenic variants for autosomal recessive non-syndromic hearing impairment by ancestry

Hearing impairment (HI) is characterized by extensive genetic heterogeneity. To determine the population-specific contribution of known autosomal recessive nonsyndromic (ARNS)HI genes and variants to HI etiology; pathogenic and likely pathogenic (PLP) ARNSHI variants were selected from ClinVar and the Deafness Variation Database and their frequencies were obtained from gnomAD for seven populations. ARNSHI prevalence due to PLP variants varies greatly by population ranging from 96.9 affected per 100,000 individuals for Ashkenazi Jews to 5.2 affected per 100,000 individuals for Africans/African Americans. For Europeans, Finns have the lowest prevalence due to ARNSHI PLP variants with 9.5 affected per 100,000 individuals. For East Asians, Latinos, non-Finish Europeans, and South Asians, ARNSHI prevalence due to PLP variants ranges from 17.1 to 33.7 affected per 100,000 individuals. ARNSHI variants that were previously reported in a single ancestry or family were observed in additional populations, e.g., USH1C p.(Q723*) reported in a Chinese family was the most prevalent pathogenic variant observed in gnomAD for African/African Americans. Variability between populations is due to how extensively ARNSHI has been studied, ARNSHI prevalence and ancestry specific ARNSHI variant architecture which is impacted by population history. Our study demonstrates that additional gene and variant discovery studies are necessary for all populations and particularly for individuals of African ancestry.

Substructured Population Growth in the Ashkenazi Jews Inferred with Approximate Bayesian Computation

The Ashkenazi Jews (AJ) are a population isolate sharing ancestry with both European and Middle Eastern populations that has likely resided in Central Europe since at least the tenth century. Between the 11th and 16th centuries, the AJ population expanded eastward leading to two culturally distinct communities in Western/Central and Eastern Europe. Our aim was to determine whether the western and eastern groups are genetically distinct, and if so, what demographic processes contributed to population differentiation. We used Approximate Bayesian Computation to choose among models of AJ history and to infer demographic parameter values, including divergence times, effective population sizes, and levels of gene flow. For the ABC analysis, we used allele frequency spectrum and identical by descent-based statistics to capture information on a wide timescale. We also mitigated the effects of ascertainment bias when performing ABC on SNP array data by jointly modeling and inferring SNP discovery. We found that the most likely model was population differentiation between Eastern and Western AJ ∼400 years ago. The differentiation between the Eastern and Western AJ could be attributed to more extreme population growth in the Eastern AJ (0.250 per generation) than the Western AJ (0.069 per generation).

A Transient Pulse of Genetic Admixture from the Crusaders in the Near East Identified from Ancient Genome Sequences

During the medieval period, hundreds of thousands of Europeans migrated to the Near East to take part in the Crusades, and many of them settled in the newly established Christian states along the Eastern Mediterranean coast. Here, we present a genetic snapshot of these events and their aftermath by sequencing the whole genomes of 13 individuals who lived in what is today known as Lebanon between the 3^rd and 13^th centuries CE. These include nine individuals from the “Crusaders’ pit” in Sidon, a mass burial in South Lebanon identified from the archaeology as the grave of Crusaders killed during a battle in the 13^th century CE. We show that all of the Crusaders’ pit individuals were males; some were Western Europeans from diverse origins, some were locals (genetically indistinguishable from present-day Lebanese), and two individuals were a mixture of European and Near Eastern ancestries, providing direct evidence that the Crusaders admixed with the local population. However, these mixtures appear to have had limited genetic consequences since signals of admixture with Europeans are not significant in any Lebanese group today—in particular, Lebanese Christians are today genetically similar to local people who lived during the Roman period which preceded the Crusades by more than four centuries.

Family history of cancer, Ashkenazi Jewish ancestry, and pancreatic cancer risk

Background

Individuals with a family history of cancer may be at increased risk of pancreatic cancer. Ashkenazi Jewish (AJ) individuals carry increased risk for pancreatic cancer and other cancer types.

Methods

We examined the association between family history of cancer, AJ heritage, and incident pancreatic cancer in 49 410 male participants of the prospective Health Professionals Follow-up Study. Hazard ratios (HRs) were estimated using multivariable-adjusted Cox proportional hazards models.

Results

During 1.1 million person-years (1986–2016), 452 participants developed pancreatic cancer. Increased risk of pancreatic cancer was observed in individuals with a family history of pancreatic (HR, 2.79; 95% confidence interval [CI], 1.28–6.07) or breast cancer (HR, 1.40; 95% CI, 1.01–1.94). There was a trend towards higher risk of pancreatic cancer in relation to a family history of colorectal cancer (HR, 1.21; 95% CI, 0.95–1.55) or AJ heritage (HR, 1.29; 95% CI, 0.94–1.77). The risk was highly elevated among AJ men with a family history of breast or colorectal cancer (HR, 2.61 [95% CI, 1.41–4.82] and 1.92 [95% CI, 1.05–3.49], respectively).

Conclusion

Family history of pancreatic cancer was associated with increased risk of this malignancy. Family history of breast or colorectal cancer was associated with the increased risk among AJ men.

Assessment of coding region variants in Kuwaiti population: implications for medical genetics and population genomics

Consanguineous populations of the Arabian Peninsula have been underrepresented in global efforts that catalogue human exome variability. We sequenced 291 whole exomes of unrelated, healthy native Arab individuals from Kuwait to a median coverage of 45X and characterised 170,508 single-nucleotide variants (SNVs), of which 21.7% were ‘personal’. Up to 12% of the SNVs were novel and 36% were population-specific. Half of the SNVs were rare and 54% were missense variants. The study complemented the Greater Middle East Variome by way of reporting many additional Arabian exome variants. The study corroborated Kuwaiti population genetic substructures previously derived using genome-wide genotype data and illustrated the genetic relatedness among Kuwaiti population subgroups, Middle Eastern, European and Ashkenazi Jewish populations. The study mapped 112 rare and frequent functional variants relating to pharmacogenomics and disorders (recessive and common) to the phenotypic characteristics of Arab population. Comparative allele frequency data and carrier distributions of known Arab mutations for 23 disorders seen among Arabs, of putative OMIM-listed causal mutations for 12 disorders observed among Arabs but not yet characterized for genetic basis in Arabs, and of 17 additional putative mutations for disorders characterized for genetic basis in Arab populations are presented for testing in future Arab studies.

A Genomic Variation Map Provides Insights into the Genetic Basis of Spring Chinese Cabbage (Brassica rapa ssp. pekinensis) Selection

Chinese cabbage is the most consumed leafy crop in East Asian countries. However, premature bolting induced by continuous low temperatures severely decreases the yield and quality of the Chinese cabbage, and therefore restricts its planting season and geographic distribution. In the past 40 years, spring Chinese cabbage with strong winterness has been selected to meet the market demand. Here, we report a genome variation map of Chinese cabbage generated from the resequencing data of 194 geographically diverse accessions of three ecotypes. In-depth analyses of the selection sweeps and genome-wide patterns revealed that spring Chinese cabbage was selected from a specific population of autumn Chinese cabbage around the area of Shandong peninsula in northern China. We identified 23 genomic loci that underwent intensive selection, and further demonstrated by gene expression and haplotype analyses that the incorporation of elite alleles of VERNALISATION INSENTIVE 3.1 (BrVIN3.1) and FLOWER LOCUS C 1 (BrFLC1) is a determinant genetic source of variation during selection. Moreover, we showed that the quantitative response of BrVIN3.1 to cold due to the sequence variations in the cis elements of the BrVIN3.1 promoter significantly contributes to bolting-time variation in Chinese cabbage. Collectively, our study provides valuable insights into the genetic basis of spring Chinese cabbage selection and will facilitate the breeding of bolting-resistant varieties by molecular-marker-assisted selection, transgenic or gene editing approaches.

Contribution of rare inherited and de novo variants in 2,871 congenital heart disease probands

Congenital heart disease (CHD) is the leading cause of mortality from birth defects. Here, exome sequencing of a single cohort of 2,871 CHD probands, including 2,645 parent-offspring trios, implicated rare inherited mutations in 1.8%, including a recessive founder mutation in GDF1 accounting for ∼5% of severe CHD in Ashkenazim, recessive genotypes in MYH6 accounting for ∼11% of Shone complex, and dominant FLT4 mutations accounting for 2.3% of Tetralogy of Fallot. De novo mutations (DNMs) accounted for 8% of cases, including ∼3% of isolated CHD patients and ∼28% with both neurodevelopmental and extra-cardiac congenital anomalies. Seven genes surpassed thresholds for genome-wide significance, and 12 genes not previously implicated in CHD had >70% probability of being disease related. DNMs in ∼440 genes were inferred to contribute to CHD. Striking overlap between genes with damaging DNMs in probands with CHD and autism was also found.

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.

The Genetic Legacy of Zoroastrianism in Iran and India: Insights into Population Structure, Gene Flow, and Selection

Zoroastrianism is one of the oldest extant religions in the world, originating in Persia (present-day Iran) during the second millennium BCE. Historical records indicate that migrants from Persia brought Zoroastrianism to India, but there is debate over the timing of these migrations. Here we present genome-wide autosomal, Y chromosome, and mitochondrial DNA data from Iranian and Indian Zoroastrians and neighboring modern-day Indian and Iranian populations and conduct a comprehensive genome-wide genetic analysis in these groups. Using powerful haplotype-based techniques, we find that Zoroastrians in Iran and India have increased genetic homogeneity relative to other sampled groups in their respective countries, consistent with their current practices of endogamy. Despite this, we infer that Indian Zoroastrians (Parsis) intermixed with local groups sometime after their arrival in India, dating this mixture to 690–1390 CE and providing strong evidence that Iranian Zoroastrian ancestry was maintained primarily through the male line. By making use of the rich information in DNA from ancient human remains, we also highlight admixture in the ancestors of Iranian Zoroastrians dated to 570 BCE–746 CE, older than admixture seen in any other sampled Iranian group, consistent with a long-standing isolation of Zoroastrians from outside groups. Finally, we report results, and challenges, from a genome-wide scan to identify genomic regions showing signatures of positive selection in present-day Zoroastrians that might correlate to the prevalence of particular diseases among these communities.

Differential analysis of mutations in the Jewish population and their implications for diseases

Sequencing large cohorts of ethnically homogeneous individuals yields genetic insights with implications for the entire population rather than a single individual. In order to evaluate the genetic basis of certain diseases encountered at high frequency in the Ashkenazi Jewish population (AJP), as well as to improve variant annotation among the AJP, we examined the entire exome, focusing on specific genes with known clinical implications in 128 Ashkenazi Jews and compared these data to other non-Jewish populations (European, African, South Asian and East Asian). We targeted American College of Medical Genetics incidental finding recommended genes and the Catalogue of Somatic Mutations in Cancer (COSMIC) germline cancer-related genes. We identified previously known disease-causing variants and discovered potentially deleterious variants in known disease-causing genes that are population specific or substantially more prevalent in the AJP, such as in the ATP and HGFAC genes associated with colorectal cancer and pancreatic cancer, respectively. Additionally, we tested the advantage of utilizing the database of the AJP when assigning pathogenicity to rare variants of independent whole-exome sequencing data of 49 Ashkenazi Jew early-onset breast cancer (BC) patients. Importantly, population-based filtering using our AJP database enabled a reduction in the number of potential causal variants in the BC cohort by 36%. Taken together, population-specific sequencing of the AJP offers valuable, clinically applicable information and improves AJP filter annotation.

The time and place of European admixture in Ashkenazi Jewish history

The Ashkenazi Jewish (AJ) population is important in genetics due to its high rate of Mendelian disorders. AJ appeared in Europe in the 10^th century, and their ancestry is thought to comprise European (EU) and Middle-Eastern (ME) components. However, both the time and place of admixture are subject to debate. Here, we attempt to characterize the AJ admixture history using a careful application of new and existing methods on a large AJ sample. Our main approach was based on local ancestry inference, in which we first classified each AJ genomic segment as EU or ME, and then compared allele frequencies along the EU segments to those of different EU populations. The contribution of each EU source was also estimated using GLOBETROTTER and haplotype sharing. The time of admixture was inferred based on multiple statistics, including ME segment lengths, the total EU ancestry per chromosome, and the correlation of ancestries along the chromosome. The major source of EU ancestry in AJ was found to be Southern Europe (≈60–80% of EU ancestry), with the rest being likely Eastern European. The inferred admixture time was ≈30 generations ago, but multiple lines of evidence suggest that it represents an average over two or more events, pre- and post-dating the founder event experienced by AJ in late medieval times. The time of the pre-bottleneck admixture event, which was likely Southern European, was estimated to ≈25–50 generations ago.

The Ashkenazi Jewish population has resided in Europe for much of its 1000-year existence. However, its ethnic and geographic origins are controversial, due to the scarcity of reliable historical records. Previous genetic studies have found links to Middle-Eastern and European ancestries, but the admixture history has not been studied in detail yet, partly due to technical difficulties in disentangling signals from multiple admixture events. Here, we present an in-depth analysis of the sources of European gene flow and the time of admixture events by using multiple new and existing methods and extensive simulations. Our results suggest a model of at least two events of European admixture. One event slightly pre-dated a late medieval founder event and was likely from a Southern European source. Another event post-dated the founder event and likely occurred in Eastern Europe. These results, as well as the methods introduced, will be highly valuable for geneticists and other researchers interested in Ashkenazi Jewish origins.

Genetic structure in the Sherpa and neighboring Nepalese populations

Background

We set out to describe the fine-scale population structure across the Eastern region of Nepal. To date there is relatively little known about the genetic structure of the Sherpa residing in Nepal and their genetic relationship with the Nepalese. We assembled dense genotype data from a total of 1245 individuals representing Nepal and a variety of different populations resident across the greater Himalayan region including Tibet, China, India, Pakistan, Kazakhstan, Uzbekistan, Tajikistan and Kirghizstan. We performed analysis of principal components, admixture and homozygosity.

Results

We identified clear substructure across populations resident in the Himalayan arc, with genetic structure broadly mirroring geographical features of the region. Ethnic subgroups within Nepal show distinct genetic structure, on both admixture and principal component analysis. We detected differential proportions of ancestry from northern Himalayan populations across Nepalese subgroups, with the Nepalese Rai, Magar and Tamang carrying the greatest proportions of Tibetan ancestry.

Conclusions

We show that populations dwelling on the Himalayan plateau have had a clear impact on the Northern Indian gene pool. We illustrate how the Sherpa are a remarkably isolated population, with little gene flow from surrounding Nepalese populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3469-5) contains supplementary material, which is available to authorized users.

Factor structure of cognition and functional capacity in two studies of schizophrenia and bipolar disorder: Implications for genomic studies

OBJECTIVE

Impairments in cognition and everyday functioning are common in schizophrenia and bipolar disorder (BPD). In this article, we present factor analyses of cognitive and functional capacity (FC) measures based on 2 studies of schizophrenia (SCZ) and bipolar I disorder (BPI) using similar methods. The overall goal of these analyses was to determine whether performance-based assessments should be examined individually, or aggregated on the basis of the correlational structure of the tests, as well as to evaluate the similarity of factor structures of SCZ and BPI.

METHOD

Veterans Affairs Cooperative Studies Program Study #572 (Harvey et al., 2014) evaluated cognitive and FC measures among 5,414 BPI and 3,942 SCZ patients. A 2nd study evaluated similar neuropsychological (NP) and FC measures among 368 BPI and 436 SCZ patients. Principal components analysis, as well as exploratory and CFAs, were used to examine the data.

RESULTS

Analyses in both datasets suggested that NP and FC measures were explained by a single underlying factor in BPI and SCZ patients, both when analyzed separately or as in a combined sample. The factor structure in both studies was similar, with or without inclusion of FC measures; homogeneous loadings were observed for that single factor across cognitive and FC domains across the samples.

CONCLUSION

The empirically derived factor model suggests that NP performance and FC are best explained as a single latent trait applicable to people with SCZ and BPD. This single measure may enhance the robustness of the analyses relating genomic data to performance-based phenotypes.

A Genome-Wide Search for Greek and Jewish Admixture in the Kashmiri Population

The Kashmiri population is an ethno-linguistic group that resides in the Kashmir Valley in northern India. A longstanding hypothesis is that this population derives ancestry from Jewish and/or Greek sources. There is historical and archaeological evidence of ancient Greek presence in India and Kashmir. Further, some historical accounts suggest ancient Hebrew ancestry as well. To date, it has not been determined whether signatures of Greek or Jewish admixture can be detected in the Kashmiri population. Using genome-wide genotyping and admixture detection methods, we determined there are no significant or substantial signs of Greek or Jewish admixture in modern-day Kashmiris. The ancestry of Kashmiri Tibetans was also determined, which showed signs of admixture with populations from northern India and west Eurasia. These results contribute to our understanding of the existing population structure in northern India and its surrounding geographical areas.

In Search of the jüdische Typus: A Proposed Benchmark to Test the Genetic Basis of Jewishness Challenges Notions of "Jewish Biomarkers"

The debate as to whether Jewishness is a biological trait inherent from an "authentic" "Jewish type" (jüdische Typus) ancestor or a system of beliefs has been raging for over two centuries. While the accumulated biological and anthropological evidence support the latter argument, recent genetic findings, bolstered by the direct-to-consumer genetic industry, purport to identify Jews or quantify one's Jewishness from genomic data. To test the merit of claims that Jews and non-Jews are genetically distinguishable, we propose a benchmark where genomic data of Jews and non-Jews are hybridized over two generations and the observed and predicted Jewishness of the terminal offspring according to either the Orthodox religious law (Halacha) or the Israeli Law of Return are compared. Members of academia, the public, and 23andMe were invited to use the benchmark to test claims that Jews are genetically distinct from non-Jews. Here, we report the findings from these trials. We also compare the genomic similarity of ∼300 individuals from nearly thirty Afro-Eurasian Jewish communities to a simulated jüdische Typus population. The results are discussed in light of modern trends in the genetics of Jews and related fields and provide a tentative answer to the ageless question "who is a Jew?"

The genetic history of Cochin Jews from India

Cochin Jews form a small and unique community on the Malabar coast in southwest India. While the arrival time of any putative Jewish ancestors of the community has been speculated to have taken place as far back as biblical times (King Solomon’s era), a Jewish community in the Malabar coast has been documented only since the 9th century CE. Here, we explore the genetic history of Cochin Jews by collecting and genotyping 21 community members and combining the data with that of 707 individuals from 72 other Indian, Jewish, and Pakistani populations, together with additional individuals from worldwide populations. We applied comprehensive genome-wide analyses based on principal component analysis, F_ST, ADMIXTURE, identity-by-descent sharing, admixture linkage disequilibrium decay, haplotype sharing, allele sharing autocorrelation decay and contrasting the X chromosome with the autosomes. We find that, as reported by several previous studies, the genetics of Cochin Jews resembles that of local Indian populations. However, we also identify considerable Jewish genetic ancestry that is not present in any other Indian or Pakistani populations (with the exception of the Jewish Bene Israel, which we characterized previously). Combined, Cochin Jews have both Jewish and Indian ancestry. Specifically, we detect a significant recent Jewish gene flow into this community 13–22 generations (~470–730 years) ago, with contributions from Yemenite, Sephardi, and Middle-Eastern Jews, in accordance with historical records. Genetic analyses also point to high endogamy and a recent population bottleneck in this population, which might explain the increased prevalence of some recessive diseases in Cochin Jews.

Localizing Ashkenazic Jews to Primeval Villages in the Ancient Iranian Lands of Ashkenaz

The Yiddish language is over 1,000 years old and incorporates German, Slavic, and Hebrew elements. The prevalent view claims Yiddish has a German origin, whereas the opposing view posits a Slavic origin with strong Iranian and weak Turkic substrata. One of the major difficulties in deciding between these hypotheses is the unknown geographical origin of Yiddish speaking Ashkenazic Jews (AJs). An analysis of 393 Ashkenazic, Iranian, and mountain Jews and over 600 non-Jewish genomes demonstrated that Greeks, Romans, Iranians, and Turks exhibit the highest genetic similarity with AJs. The Geographic Population Structure analysis localized most AJs along major primeval trade routes in northeastern Turkey adjacent to primeval villages with names that may be derived from "Ashkenaz." Iranian and mountain Jews were localized along trade routes on the Turkey's eastern border. Loss of maternal haplogroups was evident in non-Yiddish speaking AJs. Our results suggest that AJs originated from a Slavo-Iranian confederation, which the Jews call "Ashkenazic" (i.e., "Scythian"), though these Jews probably spoke Persian and/or Ossete. This is compatible with linguistic evidence suggesting that Yiddish is a Slavic language created by Irano-Turko-Slavic Jewish merchants along the Silk Roads as a cryptic trade language, spoken only by its originators to gain an advantage in trade. Later, in the 9th century, Yiddish underwent relexification by adopting a new vocabulary that consists of a minority of German and Hebrew and a majority of newly coined Germanoid and Hebroid elements that replaced most of the original Eastern Slavic and Sorbian vocabularies, while keeping the original grammars intact.

Genome-wide association study of schizophrenia in Ashkenazi Jews

Schizophrenia is a common, clinically heterogeneous disorder associated with lifelong morbidity and early mortality. Several genetic variants associated with schizophrenia have been identified, but the majority of the heritability remains unknown. In this study, we report on a case-control sample of Ashkenazi Jews (AJ), a founder population that may provide additional insights into genetic etiology of schizophrenia. We performed a genome-wide association analysis (GWAS) of 592 cases and 505 controls of AJ ancestry ascertained in the US. Subsequently, we performed a meta-analysis with an Israeli AJ sample of 913 cases and 1640 controls, followed by a meta-analysis and polygenic risk scoring using summary results from Psychiatric GWAS Consortium 2 schizophrenia study. The U.S. AJ sample showed strong evidence of polygenic inheritance (pseudo-R(2) ∼9.7%) and a SNP-heritability estimate of 0.39 (P = 0.00046). We found no genome-wide significant associations in the U.S. sample or in the combined US/Israeli AJ meta-analysis of 1505 cases and 2145 controls. The strongest AJ specific associations (P-values in 10(-6) -10(-7) range) were in the 22q 11.2 deletion region and included the genes TBX1, GLN1, and COMT. Supportive evidence (meta P < 1 × 10(-4) ) was also found for several previously identified genome-wide significant findings, including the HLA region, CNTN4, IMMP2L, and GRIN2A. The meta-analysis of the U.S. sample with the PGC2 results provided initial genome-wide significant evidence for six new loci. Among the novel potential susceptibility genes is PEPD, a gene involved in proline metabolism, which is associated with a Mendelian disorder characterized by developmental delay and cognitive deficits.

Microarray Analysis of Copy Number Variants on the Human Y Chromosome Reveals Novel and Frequent Duplications Overrepresented in Specific Haplogroups

BACKGROUND

The human Y chromosome is almost always excluded from genome-wide investigations of copy number variants (CNVs) due to its highly repetitive structure. This chromosome should not be forgotten, not only for its well-known relevance in male fertility, but also for its involvement in clinical phenotypes such as cancers, heart failure and sex specific effects on brain and behaviour.

RESULTS

We analysed Y chromosome data from Affymetrix 6.0 SNP arrays and found that the signal intensities for most of 8179 SNP/CN probes in the male specific region (MSY) discriminated between a male, background signals in a female and an isodicentric male containing a large deletion of the q-arm and a duplication of the p-arm of the Y chromosome. Therefore, this SNP/CN platform is suitable for identification of gain and loss of Y chromosome sequences. In a set of 1718 males, we found 25 different CNV patterns, many of which are novel. We confirmed some of these variants by PCR or qPCR. The total frequency of individuals with CNVs was 14.7%, including 9.5% with duplications, 4.5% with deletions and 0.7% exhibiting both. Hence, a novel observation is that the frequency of duplications was more than twice the frequency of deletions. Another striking result was that 10 of the 25 detected variants were significantly overrepresented in one or more haplogroups, demonstrating the importance to control for haplogroups in genome-wide investigations to avoid stratification. NO-M214(xM175) individuals presented the highest percentage (95%) of CNVs. If they were not counted, 12.4% of the rest included CNVs, and the difference between duplications (8.9%) and deletions (2.8%) was even larger.

CONCLUSIONS

Our results demonstrate that currently available genome-wide SNP platforms can be used to identify duplications and deletions in the human Y chromosome. Future association studies of the full spectrum of Y chromosome variants will demonstrate the potential involvement of gain or loss of Y chromosome sequence in different human phenotypes.

The peopling of the African continent and the diaspora into the new world

Africa is the birthplace of anatomically modern humans, and is the geographic origin of human migration across the globe within the last 100,000 years. The history of African populations has consisted of a number of demographic events that have influenced patterns of genetic and phenotypic variation across the continent. With the increasing amount of genomic data and corresponding developments in computational methods, researchers are able to explore long-standing evolutionary questions, expanding our understanding of human history within and outside of Africa. This review will summarize some of the recent findings regarding African demographic history, including the African Diaspora, and will briefly explore their implications for disease susceptibility in populations of African descent.