Mapping genomic loci implicates genes and synaptic biology in schizophrenia

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Centenarian Exomes as a Tool for Evaluating the Clinical Relevance of Germline Tumor Suppressor Mutations

Objectives:

The aim of the present study was to evaluate the clinical relevance of mutations in tumor suppressor genes using whole-exome sequencing data from centenarians and young healthy individuals.

Methods:

Two pools, one of centenarians and one of young individuals, were constructed and whole-exome sequencing was performed. We examined the whole-exome sequencing data of Bulgarian individuals for carriership of tumor suppressor gene variants.

Results:

Of all variants annotated in both pools, 5080 (0.06%) are variants in tumor suppressor genes but only 46 show significant difference in allele frequencies between the two studied groups. Four variants (0.004%) are pathogenic/risk factors according to single nucleotide polymorphism database: rs1566734 in PTPRJ, rs861539 in XRCC3, rs203462 in AKAP10, and rs486907 in RNASEL.

Discussion:

Based on their high minor allele frequencies and presence in the centenarian group, we could reclassify them from pathogenic/risk factors to benign. Our study shows that centenarian exomes can be used for re-evaluating the clinically uncertain variants.

Ancient mitochondrial DNA pathogenic variants putatively associated with mitochondrial disease

Mitochondrial DNA variants associated with diseases are widely studied in contemporary populations, but their prevalence has not yet been investigated in ancient populations. The publicly available AmtDB database contains 1443 ancient mtDNA Eurasian genomes from different periods. The objective of this study was to use this data to establish the presence of pathogenic mtDNA variants putatively associated with mitochondrial diseases in ancient populations. The clinical significance, pathogenicity prediction and contemporary frequency of mtDNA variants were determined using online platforms. The analyzed ancient mtDNAs contain six variants designated as being "confirmed pathogenic" in modern patients. The oldest of these, m.7510T>C in the MT-TS1 gene, was found in a sample from the Neolithic period, dated 5800-5400 BCE. All six have well established clinical association, and their pathogenic effect is corroborated by very low population frequencies in contemporary populations. Analysis of the geographic location of the ancient samples, contemporary epidemiological trends and probable haplogroup association indicate diverse spatiotemporal dynamics of these variants. The dynamics in the prevalence and distribution is conceivably result of de novo mutations or human migrations and subsequent evolutionary processes. In addition, ten variants designated as possibly or likely pathogenic were found, but the clinical effect of these is not yet well established and further research is warranted. All detected mutations putatively associated with mitochondrial disease in ancient mtDNA samples are in tRNA coding genes. Most of these mutations are in a mt-tRNA type (Model 2) that is characterized by loss of D-loop/T-loop interaction. Exposing pathogenic variants in ancient human populations expands our understanding of their origin and prevalence dynamics.

Ancient human mitochondrial genomes from Bronze Age Bulgaria: new insights into the genetic history of Thracians

One of the best documented Indo-European civilizations that inhabited Bulgaria is the Thracians, who lasted for more than five millennia and whose origin and relationships with other past and present-day populations are debated among researchers. Here we report 25 new complete mitochondrial genomes of ancient individuals coming from three necropolises located in different regions of Bulgaria – Shekerdja mogila, Gabrova mogila and Bereketska mogila – dated to II-III millennium BC. The identified mtDNA haplogroup composition reflects the mitochondrial variability of Western Eurasia. In particular, within the ancient Eurasian genetic landscape, Thracians locate in an intermediate position between Early Neolithic farmers and Late Neolithic-Bronze Age steppe pastoralists, supporting the scenario that the Balkan region has been a link between Eastern Europe and the Mediterranean since the prehistoric time. Spatial Principal Component Analysis (sPCA) performed on Thracian and modern mtDNA sequences, confirms the pattern highlighted on ancient populations, overall indicating that the maternal gene pool of Thracians reflects their central geographical position at the gateway of Europe.

Analysis of shared heritability in common disorders of the brain

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

Estimation of Genetic Correlation via Linkage Disequilibrium Score Regression and Genomic Restricted Maximum Likelihood

Genetic correlation is a key population parameter that describes the shared genetic architecture of complex traits and diseases. It can be estimated by current state-of-art methods, i.e., linkage disequilibrium score regression (LDSC) and genomic restricted maximum likelihood (GREML). The massively reduced computing burden of LDSC compared to GREML makes it an attractive tool, although the accuracy (i.e., magnitude of standard errors) of LDSC estimates has not been thoroughly studied. In simulation, we show that the accuracy of GREML is generally higher than that of LDSC. When there is genetic heterogeneity between the actual sample and reference data from which LD scores are estimated, the accuracy of LDSC decreases further. In real data analyses estimating the genetic correlation between schizophrenia (SCZ) and body mass index, we show that GREML estimates based on ∼150,000 individuals give a higher accuracy than LDSC estimates based on ∼400,000 individuals (from combined meta-data). A GREML genomic partitioning analysis reveals that the genetic correlation between SCZ and height is significantly negative for regulatory regions, which whole genome or LDSC approach has less power to detect. We conclude that LDSC estimates should be carefully interpreted as there can be uncertainty about homogeneity among combined meta-datasets. We suggest that any interesting findings from massive LDSC analysis for a large number of complex traits should be followed up, where possible, with more detailed analyses with GREML methods, even if sample sizes are lesser.

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.

Mitochondrial DNA Suggests a Western Eurasian Origin for Ancient (Proto-) Bulgarians

Ancient (proto-) Bulgarians have long been thought of as a Turkic population. However, evidence found in the past three decades shows that this is not the case. Until now, this evidence has not included ancient mitochondrial DNA (mtDNA) analysis. To fill this void, we collected human remains from the 8th to the 10th century AD located in three necropolises in Bulgaria: Nojarevo (Silistra region) and Monastery of Mostich (Shumen region), both in northeastern Bulgaria, and Tuhovishte (Satovcha region) in southwestern Bulgaria. The phylogenetic analysis of 13 ancient DNA samples (extracted from teeth) identified 12 independent haplotypes, which we further classified into mtDNA haplogroups found in present-day European and western Eurasian populations. Our results suggest a western Eurasian matrilineal origin for proto-Bulgarians, as well as a genetic similarity between proto- and modern Bulgarians. Our future work will provide additional data that will further clarify proto-Bulgarian origins, thereby adding new clues to the current understanding of European genetic evolution.

Epigenetic alterations in patients with type 2 diabetes mellitus

Epigenetic changes, in particular DNA methylation processes, play a role in the pathogenesis and progression of type 2 diabetes mellitus (T2DM) linking genetic and environmental factors. To clarify this role, we have analyzed in patients with different duration of T2DM: (i) expression levels of methyl-CpG-binding domain protein 2 (MBD2) as marker of DNA methylation, and ii) methylation changes in 22 genes connected to cellular stress and toxicity. We have analyzed MBD2 mRNA expression levels in16 patients and 12 controls and the methylation status of stress and toxicity genes in four DNA pools: (i) controls; (ii) newly-diagnosed T2DM patients; (iii) patients with T2DM duration of <5 years and (iv) of >5 years. The MBD2 expression levels were 10.4-times increased on average in T2DM patients compared to controls. Consistent increase in DNA methylation fraction with the increase in T2DM duration was observed in Prdx2 and SCARA3 genes, connected to oxidative stress protection and in BRCA1 and Tp53 tumor-suppressor genes. In conclusion, increased MBD2 expression in patients indicated general dysregulation of DNA methylation in T2DM. The elevated methylation of Prdx2 and SCARA3 genes suggests disturbance in oxidative stress protection in T2DM. The increased methylation of BRCA1 and Tp53 genes unraveled an epigenetic cause for T2DM related increase in cancer risk.

Modeling Linkage Disequilibrium Increases Accuracy of Polygenic Risk Scores

Polygenic risk scores have shown great promise in predicting complex disease risk and will become more accurate as training sample sizes increase. The standard approach for calculating risk scores involves linkage disequilibrium (LD)-based marker pruning and applying a p value threshold to association statistics, but this discards information and can reduce predictive accuracy. We introduce LDpred, a method that infers the posterior mean effect size of each marker by using a prior on effect sizes and LD information from an external reference panel. Theory and simulations show that LDpred outperforms the approach of pruning followed by thresholding, particularly at large sample sizes. Accordingly, predicted R(2) increased from 20.1% to 25.3% in a large schizophrenia dataset and from 9.8% to 12.0% in a large multiple sclerosis dataset. A similar relative improvement in accuracy was observed for three additional large disease datasets and for non-European schizophrenia samples. The advantage of LDpred over existing methods will grow as sample sizes increase.

The Role of Recent Admixture in Forming the Contemporary West Eurasian Genomic Landscape

Over the past few years, studies of DNA isolated from human fossils and archaeological remains have generated considerable novel insight into the history of our species. Several landmark papers have described the genomes of ancient humans across West Eurasia, demonstrating the presence of large-scale, dynamic population movements over the last 10,000 years, such that ancestry across present-day populations is likely to be a mixture of several ancient groups [1, 2, 3, 4, 5, 6, 7]. While these efforts are bringing the details of West Eurasian prehistory into increasing focus, studies aimed at understanding the processes behind the generation of the current West Eurasian genetic landscape have been limited by the number of populations sampled or have been either too regional or global in their outlook [8, 9, 10, 11]. Here, using recently described haplotype-based techniques [11], we present the results of a systematic survey of recent admixture history across Western Eurasia and show that admixture is a universal property across almost all groups. Admixture in all regions except North Western Europe involved the influx of genetic material from outside of West Eurasia, which we date to specific time periods. Within Northern, Western, and Central Europe, admixture tended to occur between local groups during the period 300 to 1200 CE. Comparisons of the genetic profiles of West Eurasians before and after admixture show that population movements within the last 1,500 years are likely to have maintained differentiation among groups. Our analysis provides a timeline of the gene flow events that have generated the contemporary genetic landscape of West Eurasia.

•

Recent admixture events involved outside groups at the edges of West Eurasia

•

Admixture within Europe tended to fall within the European Migration Period

•

West Eurasian genetic structure today is likely to have been maintained by admixture

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

In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide-variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases

Regulatory and coding variants are known to be enriched with associations identified by genome-wide association studies (GWASs) of complex disease, but their contributions to trait heritability are currently unknown. We applied variance-component methods to imputed genotype data for 11 common diseases to partition the heritability explained by genotyped SNPs (hg(2)) across functional categories (while accounting for shared variance due to linkage disequilibrium). Extensive simulations showed that in contrast to current estimates from GWAS summary statistics, the variance-component approach partitions heritability accurately under a wide range of complex-disease architectures. Across the 11 diseases DNaseI hypersensitivity sites (DHSs) from 217 cell types spanned 16% of imputed SNPs (and 24% of genotyped SNPs) but explained an average of 79% (SE = 8%) of hg(2) from imputed SNPs (5.1× enrichment; p = 3.7 × 10(-17)) and 38% (SE = 4%) of hg(2) from genotyped SNPs (1.6× enrichment, p = 1.0 × 10(-4)). Further enrichment was observed at enhancer DHSs and cell-type-specific DHSs. In contrast, coding variants, which span 1% of the genome, explained <10% of hg(2) despite having the highest enrichment. We replicated these findings but found no significant contribution from rare coding variants in independent schizophrenia cohorts genotyped on GWAS and exome chips. Our results highlight the value of analyzing components of heritability to unravel the functional architecture of common disease.

Ancient human genomes suggest three ancestral populations for present-day Europeans

We sequenced the genomes of a ~7,000 year old farmer from Germany and eight ~8,000 year old hunter-gatherers from Luxembourg and Sweden. We analyzed these and other ancient genomes^1–4 with 2,345 contemporary humans to show that most present Europeans derive from at least three highly differentiated populations: West European Hunter-Gatherers (WHG), who contributed ancestry to all Europeans but not to Near Easterners; Ancient North Eurasians (ANE) related to Upper Paleolithic Siberians³, who contributed to both Europeans and Near Easterners; and Early European Farmers (EEF), who were mainly of Near Eastern origin but also harbored WHG-related ancestry. We model these populations’ deep relationships and show that EEF had ~44% ancestry from a “Basal Eurasian” population that split prior to the diversification of other non-African lineages.