Human Disease Variation in the Light of Population Genomics.

Identifying the causes of similarities and differences in genetic disease prevalence among humans is central to understanding disease etiology. While present-day humans are not strongly differentiated, vast amounts of genomic data now make it possible to study subtle patterns of genetic variation. This allows us to trace our genomic history thousands of years into the past and its implications for the distribution of disease-associated variants today. Genomic analyses have shown that demographic processes shaped the distribution and frequency of disease-associated variants over time. Furthermore, local adaptation to new environmental conditions-including pathogens-has generated strong patterns of differentiation at particular loci. Researchers are also beginning to uncover the genetic architecture of complex diseases, affected by many variants of small effect. The field of population genomics thus holds great potential for providing further insights into the evolution of human disease.

The selection landscape and genetic legacy of ancient Eurasians.

The Holocene (beginning around 12,000 years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using a dataset of more than 1,600 imputed ancient genomes1, we modelled the selection landscape during the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify key selection signals related to metabolism, including that selection at the FADS cluster began earlier than previously reported and that selection near the LCT locus predates the emergence of the lactase persistence allele by thousands of years. We also find strong selection in the HLA region, possibly due to increased exposure to pathogens during the Bronze Age. Using ancient individuals to infer local ancestry tracts in over 400,000 samples from the UK Biobank, we identify widespread differences in the distribution of Mesolithic, Neolithic and Bronze Age ancestries across Eurasia. By calculating ancestry-specific polygenic risk scores, we show that height differences between Northern and Southern Europe are associated with differential Steppe ancestry, rather than selection, and that risk alleles for mood-related phenotypes are enriched for Neolithic farmer ancestry, whereas risk alleles for diabetes and Alzheimer's disease are enriched for Western hunter-gatherer ancestry. Our results indicate that ancient selection and migration were large contributors to the distribution of phenotypic diversity in present-day Europeans.

100 ancient genomes show repeated population turnovers in Neolithic Denmark.

Major migration events in Holocene Eurasia have been characterized genetically at broad regional scales1-4. However, insights into the population dynamics in the contact zones are hampered by a lack of ancient genomic data sampled at high spatiotemporal resolution5-7. Here, to address this, we analysed shotgun-sequenced genomes from 100 skeletons spanning 7,300 years of the Mesolithic period, Neolithic period and Early Bronze Age in Denmark and integrated these with proxies for diet (13C and 15N content), mobility (87Sr/86Sr ratio) and vegetation cover (pollen). We observe that Danish Mesolithic individuals of the Maglemose, Kongemose and Ertebølle cultures form a distinct genetic cluster related to other Western European hunter-gatherers. Despite shifts in material culture they displayed genetic homogeneity from around 10,500 to 5,900 calibrated years before present, when Neolithic farmers with Anatolian-derived ancestry arrived. Although the Neolithic transition was delayed by more than a millennium relative to Central Europe, it was very abrupt and resulted in a population turnover with limited genetic contribution from local hunter-gatherers. The succeeding Neolithic population, associated with the Funnel Beaker culture, persisted for only about 1,000 years before immigrants with eastern Steppe-derived ancestry arrived. This second and equally rapid population replacement gave rise to the Single Grave culture with an ancestry profile more similar to present-day Danes. In our multiproxy dataset, these major demographic events are manifested as parallel shifts in genotype, phenotype, diet and land use.

Population genomics of post-glacial western Eurasia.

Western Eurasia witnessed several large-scale human migrations during the Holocene1-5. Here, to investigate the cross-continental effects of these migrations, we shotgun-sequenced 317 genomes-mainly from the Mesolithic and Neolithic periods-from across northern and western Eurasia. These were imputed alongside published data to obtain diploid genotypes from more than 1,600 ancient humans. Our analyses revealed a 'great divide' genomic boundary extending from the Black Sea to the Baltic. Mesolithic hunter-gatherers were highly genetically differentiated east and west of this zone, and the effect of the neolithization was equally disparate. Large-scale ancestry shifts occurred in the west as farming was introduced, including near-total replacement of hunter-gatherers in many areas, whereas no substantial ancestry shifts happened east of the zone during the same period. Similarly, relatedness decreased in the west from the Neolithic transition onwards, whereas, east of the Urals, relatedness remained high until around 4,000 BP, consistent with the persistence of localized groups of hunter-gatherers. The boundary dissolved when Yamnaya-related ancestry spread across western Eurasia around 5,000 BP, resulting in a second major turnover that reached most parts of Europe within a 1,000-year span. The genetic origin and fate of the Yamnaya have remained elusive, but we show that hunter-gatherers from the Middle Don region contributed ancestry to them. Yamnaya groups later admixed with individuals associated with the Globular Amphora culture before expanding into Europe. Similar turnovers occurred in western Siberia, where we report new genomic data from a 'Neolithic steppe' cline spanning the Siberian forest steppe to Lake Baikal. These prehistoric migrations had profound and lasting effects on the genetic diversity of Eurasian populations.

Population genomics of the Viking world.

The maritime expansion of Scandinavian populations during the Viking Age (about AD 750-1050) was a far-flung transformation in world history1,2. Here we sequenced the genomes of 442 humans from archaeological sites across Europe and Greenland (to a median depth of about 1×) to understand the global influence of this expansion. We find the Viking period involved gene flow into Scandinavia from the south and east. We observe genetic structure within Scandinavia, with diversity hotspots in the south and restricted gene flow within Scandinavia. We find evidence for a major influx of Danish ancestry into England; a Swedish influx into the Baltic; and Norwegian influx into Ireland, Iceland and Greenland. Additionally, we see substantial ancestry from elsewhere in Europe entering Scandinavia during the Viking Age. Our ancient DNA analysis also revealed that a Viking expedition included close family members. By comparing with modern populations, we find that pigmentation-associated loci have undergone strong population differentiation during the past millennium, and trace positively selected loci-including the lactase-persistence allele of LCT and alleles of ANKA that are associated with the immune response-in detail. We conclude that the Viking diaspora was characterized by substantial transregional engagement: distinct populations influenced the genomic makeup of different regions of Europe, and Scandinavia experienced increased contact with the rest of the continent.

A comprehensive map of genetic relationships among diagnostic categories based on 48.6 million relative pairs from the Danish genealogy.

For more than half a century, Denmark has maintained population-wide demographic, health care, and socioeconomic registers that provide detailed information on the interaction between all residents and the extensive national social services system. We leverage this resource to reconstruct the genealogy of the entire nation based on all individuals legally residing in Denmark since 1968. We cross-reference 6,691,426 individuals with nationwide health care registers to estimate heritability and genetic correlations of 10 broad diagnostic categories involving all major organs and systems. Heritability estimates for mental disorders were consistently the highest across demographic cohorts (average h2 = 0.406, 95% CI = [0.403, 0.408]), whereas estimates for cancers were the lowest (average h2 = 0.130, 95% CI = [0.125, 0.134]). The average genetic correlation of each of the 10 diagnostic categories with the other nine was highest for gastrointestinal conditions (average rg = 0.567, 95% CI = [0.566, 0.567]) and lowest for urogenital conditions (average rg = 0.386, 95% CI = [0.385, 0.388]). Mental, pulmonary, gastrointestinal, and neurological conditions had similar genetic correlation profiles.

Comparing migraine with and without aura to healthy controls using RNA sequencing.

BACKGROUND: Migraine mechanisms are *These authors contributed equally to this work. only partly known. Some studies have previously described genes differentially expressed between blood from migraineurs and controls. The objective of this study was to describe gene expression in subtypes of migraine outside of attack and in healthy controls. METHODS: We extensively phenotyped 17 migraine without aura and nine migraine with aura female patients, and 20 age-matched female controls. Cubital venous blood was RNA sequenced. Genes differentially expressed between migraineurs (migraine without aura and migraine with aura) and controls, and between migraine without aura and migraine with aura were identified using a case-control design. A co-expression network was constructed to investigate the difference between migraineurs and healthy controls at the network level. RESULTS: We found two differentially expressed genes: NMNAT2 and RETN. Both were differentially expressed between migraine with aura and controls, but they could not be replicated in an independent cohort. Co-expression network analysis resulted in one cluster of highly interconnected genes that was nominally significantly associated with migraine; however, no pathways or gene ontology terms were detected. CONCLUSIONS: We showed no clear distinct difference in gene expression profiles of peripheral blood of migraineurs and controls and were not able to replicate findings from previous studies. A larger sample size may be needed to detect minor differences.

Parental Origin of Interstitial Duplications at 15q11.2-q13.3 in Schizophrenia and Neurodevelopmental Disorders.

Duplications at 15q11.2-q13.3 overlapping the Prader-Willi/Angelman syndrome (PWS/AS) region have been associated with developmental delay (DD), autism spectrum disorder (ASD) and schizophrenia (SZ). Due to presence of imprinted genes within the region, the parental origin of these duplications may be key to the pathogenicity. Duplications of maternal origin are associated with disease, whereas the pathogenicity of paternal ones is unclear. To clarify the role of maternal and paternal duplications, we conducted the largest and most detailed study to date of parental origin of 15q11.2-q13.3 interstitial duplications in DD, ASD and SZ cohorts. We show, for the first time, that paternal duplications lead to an increased risk of developing DD/ASD/multiple congenital anomalies (MCA), but do not appear to increase risk for SZ. The importance of the epigenetic status of 15q11.2-q13.3 duplications was further underlined by analysis of a number of families, in which the duplication was paternally derived in the mother, who was unaffected, whereas her offspring, who inherited a maternally derived duplication, suffered from psychotic illness. Interestingly, the most consistent clinical characteristics of SZ patients with 15q11.2-q13.3 duplications were learning or developmental problems, found in 76% of carriers. Despite their lower pathogenicity, paternal duplications are less frequent in the general population with a general population prevalence of 0.0033% compared to 0.0069% for maternal duplications. This may be due to lower fecundity of male carriers and differential survival of embryos, something echoed in the findings that both types of duplications are de novo in just over 50% of cases. Isodicentric chromosome 15 (idic15) or interstitial triplications were not observed in SZ patients or in controls. Overall, this study refines the distinct roles of maternal and paternal interstitial duplications at 15q11.2-q13.3, underlining the critical importance of maternally expressed imprinted genes in the contribution of Copy Number Variants (CNVs) at this interval to the incidence of psychotic illness. This work will have tangible benefits for patients with 15q11.2-q13.3 duplications by aiding genetic counseling.

Polygenic risk scores for schizophrenia and bipolar disorder associate with addiction.

We use polygenic risk scores (PRSs) for schizophrenia (SCZ) and bipolar disorder (BPD) to predict smoking, and addiction to nicotine, alcohol or drugs in individuals not diagnosed with psychotic disorders. Using PRSs for 144 609 subjects, including 10 036 individuals admitted for in-patient addiction treatment and 35 754 smokers, we find that diagnoses of various substance use disorders and smoking associate strongly with PRSs for SCZ (P = 5.3 × 10-50 -1.4 × 10-6 ) and BPD (P = 1.7 × 10-9 -1.9 × 10-3 ), showing shared genetic etiology between psychosis and addiction. Using standardized scores for SCZ and BPD scaled to a unit increase doubling the risk of the corresponding disorder, the odds ratios for alcohol and substance use disorders range from 1.19 to 1.31 for the SCZ-PRS, and from 1.07 to 1.29 for the BPD-PRS. Furthermore, we show that as regular smoking becomes more stigmatized and less prevalent, these biological risk factors gain importance as determinants of the behavior.

Detection of sharing by descent, long-range phasing and haplotype imputation.

Uncertainty about the phase of strings of SNPs creates complications in genetic analysis, although methods have been developed for phasing population-based samples. However, these methods can only phase a small number of SNPs effectively and become unreliable when applied to SNPs spanning many linkage disequilibrium (LD) blocks. Here we show how to phase more than 1,000 SNPs simultaneously for a large fraction of the 35,528 Icelanders genotyped by Illumina chips. Moreover, haplotypes that are identical by descent (IBD) between close and distant relatives, for example, those separated by ten meioses or more, can often be reliably detected. This method is particularly powerful in studies of the inheritance of recurrent mutations and fine-scale recombinations in large sample sets. A further extension of the method allows us to impute long haplotypes for individuals who are not genotyped.

Identification of rare high-risk copy number variants affecting the dopamine transporter gene in mental disorders.

BACKGROUND: The dopamine transporter, also known as solute carrier 6A3 (SLC6A3), plays an important role in synaptic transmission by regulating the reuptake of dopamine in the synapses. In line with this, variations in the gene encoding this transporter have been linked to both schizophrenia and affective disorders. Recently, copy number variants (CNVs) in SLC6A3 have been identified in healthy subjects but so far, the implication of CNVs affecting this gene in psychiatric diseases has not been addressed. AIMS: In the present study, we aimed to investigate whether CNVs affecting SLC6A3 represent rare high-risk variants of psychiatric disorders. METHODS: We performed a systematic screening for CNVs affecting SLC6A3 in 761 healthy controls, 672 schizophrenia patients, and 194 patients with bipolar disorder in addition to 253 family members from six large pedigrees affected by mental disorders using single nucleotide polymorphism arrays and subsequent verification by real-time polymerase chain reaction. RESULTS: We identified two duplications and one deletion affecting SLC6A3 in the patients, while no such CNVs were identified in any of the controls. The identified CNVs were of different sizes and two affected several genes in addition to SLC6A3. CONCLUSION: Our findings suggest that rare high-risk CNVs affecting the gene encoding the dopamine transporter contribute to the pathogenesis of schizophrenia and affective disorders.

Common variants conferring risk of schizophrenia.

Schizophrenia is a complex disorder, caused by both genetic and environmental factors and their interactions. Research on pathogenesis has traditionally focused on neurotransmitter systems in the brain, particularly those involving dopamine. Schizophrenia has been considered a separate disease for over a century, but in the absence of clear biological markers, diagnosis has historically been based on signs and symptoms. A fundamental message emerging from genome-wide association studies of copy number variations (CNVs) associated with the disease is that its genetic basis does not necessarily conform to classical nosological disease boundaries. Certain CNVs confer not only high relative risk of schizophrenia but also of other psychiatric disorders. The structural variations associated with schizophrenia can involve several genes and the phenotypic syndromes, or the 'genomic disorders', have not yet been characterized. Single nucleotide polymorphism (SNP)-based genome-wide association studies with the potential to implicate individual genes in complex diseases may reveal underlying biological pathways. Here we combined SNP data from several large genome-wide scans and followed up the most significant association signals. We found significant association with several markers spanning the major histocompatibility complex (MHC) region on chromosome 6p21.3-22.1, a marker located upstream of the neurogranin gene (NRGN) on 11q24.2 and a marker in intron four of transcription factor 4 (TCF4) on 18q21.2. Our findings implicating the MHC region are consistent with an immune component to schizophrenia risk, whereas the association with NRGN and TCF4 points to perturbation of pathways involved in brain development, memory and cognition.

A variant associated with nicotine dependence, lung cancer and peripheral arterial disease.

Smoking is a leading cause of preventable death, causing about 5 million premature deaths worldwide each year. Evidence for genetic influence on smoking behaviour and nicotine dependence (ND) has prompted a search for susceptibility genes. Furthermore, assessing the impact of sequence variants on smoking-related diseases is important to public health. Smoking is the major risk factor for lung cancer (LC) and is one of the main risk factors for peripheral arterial disease (PAD). Here we identify a common variant in the nicotinic acetylcholine receptor gene cluster on chromosome 15q24 with an effect on smoking quantity, ND and the risk of two smoking-related diseases in populations of European descent. The variant has an effect on the number of cigarettes smoked per day in our sample of smokers. The same variant was associated with ND in a previous genome-wide association study that used low-quantity smokers as controls, and with a similar approach we observe a highly significant association with ND. A comparison of cases of LC and PAD with population controls each showed that the variant confers risk of LC and PAD. The findings provide a case study of a gene-environment interaction, highlighting the role of nicotine addiction in the pathology of other serious diseases.

Large recurrent microdeletions associated with schizophrenia.

Reduced fecundity, associated with severe mental disorders, places negative selection pressure on risk alleles and may explain, in part, why common variants have not been found that confer risk of disorders such as autism, schizophrenia and mental retardation. Thus, rare variants may account for a larger fraction of the overall genetic risk than previously assumed. In contrast to rare single nucleotide mutations, rare copy number variations (CNVs) can be detected using genome-wide single nucleotide polymorphism arrays. This has led to the identification of CNVs associated with mental retardation and autism. In a genome-wide search for CNVs associating with schizophrenia, we used a population-based sample to identify de novo CNVs by analysing 9,878 transmissions from parents to offspring. The 66 de novo CNVs identified were tested for association in a sample of 1,433 schizophrenia cases and 33,250 controls. Three deletions at 1q21.1, 15q11.2 and 15q13.3 showing nominal association with schizophrenia in the first sample (phase I) were followed up in a second sample of 3,285 cases and 7,951 controls (phase II). All three deletions significantly associate with schizophrenia and related psychoses in the combined sample. The identification of these rare, recurrent risk variants, having occurred independently in multiple founders and being subject to negative selection, is important in itself. CNV analysis may also point the way to the identification of additional and more prevalent risk variants in genes and pathways involved in schizophrenia.

A common inversion under selection in Europeans.

A refined physical map of chromosome 17q21.31 uncovered a 900-kb inversion polymorphism. Chromosomes with the inverted segment in different orientations represent two distinct lineages, H1 and H2, that have diverged for as much as 3 million years and show no evidence of having recombined. The H2 lineage is rare in Africans, almost absent in East Asians but found at a frequency of 20% in Europeans, in whom the haplotype structure is indicative of a history of positive selection. Here we show that the H2 lineage is undergoing positive selection in the Icelandic population, such that carrier females have more children and have higher recombination rates than noncarriers.

Sequence analysis of 17 NRXN1 deletions.

BACKGROUND: Genome instability plays fundamental roles in human evolution and phenotypic variation within our population. This instability leads to genomic rearrangements that are involved in a wide variety of human disorders, including congenital and neurodevelopmental disorders, and cancers. Insight into the molecular mechanisms governing such genomic rearrangements may increase our understanding of disease pathology and evolutionary processes. Here we analyse 17 carriers of non-recurrent deletions in the NRXN1 gene, which have been associated with neurodevelopmental disorders, e.g. schizophrenia, autism and epilepsies. METHODS: 17 non-recurrent NRXN1 deletions identified by GWA were sequenced to map the breakpoints of each. Meme … etc. was used to identify shared patterns between the deletions and compare these were previously studies on non-recurrent deletions. RESULTS: We discovered two novel sequence motifs shared between all 17 NRXN1 deletions and a significantly higher AT nucleotide content at the breakpoints, compared to the overall nucleotide content on chromosome 2. We found different alteration of sequence at the breakpoint; small insertions and duplications giving rise to short microhomology sequences. CONCLUSIONS: No single mechanism seems to be implicated in the deletion events, but the results suggest that NHEJ, FoSTeS or MMBIR is implicated. The two novel sequence motifs together with a high AT content in all in NRXN1 deletions may lead to increased instability leading to a increase susceptibility to a single stranded structures. This favours potentially repaired by NHEJ mechanism of double strand breaks or may leading to replication errors. © 2013 Wiley Periodicals, Inc.

Polygenic profiles define aspects of clinical heterogeneity in attention deficit hyperactivity disorder.

Attention deficit hyperactivity disorder (ADHD) is a complex disorder that manifests variability in long-term outcomes and clinical presentations. The genetic contributions to such heterogeneity are not well understood. Here we show several genetic links to clinical heterogeneity in ADHD in a case-only study of 14,084 diagnosed individuals. First, we identify one genome-wide significant locus by comparing cases with ADHD and autism spectrum disorder (ASD) to cases with ADHD but not ASD. Second, we show that cases with ASD and ADHD, substance use disorder and ADHD, or first diagnosed with ADHD in adulthood have unique polygenic score (PGS) profiles that distinguish them from complementary case subgroups and controls. Finally, a PGS for an ASD diagnosis in ADHD cases predicted cognitive performance in an independent developmental cohort. Our approach uncovered evidence of genetic heterogeneity in ADHD, helping us to understand its etiology and providing a model for studies of other disorders.

Population-Based Risk of Psychiatric Disorders Associated With Recurrent Copy Number Variants.

Importance: Recurrent copy number variants (rCNVs) have been associated with increased risk of psychiatric disorders in case-control studies, but their population-level impact is unknown. Objective: To provide unbiased population-based estimates of prevalence and risk associated with psychiatric disorders for rCNVs and to compare risks across outcomes, rCNV dosage type (deletions or duplications), and locus features. Design, Setting, and Participants: This genetic association study is an analysis of data from the Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH) case-cohort sample of individuals born in Denmark in 1981-2008 and followed up until 2015, including (1) all individuals (n = 92 531) with a hospital discharge diagnosis of attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder, major depressive disorder (MDD), or schizophrenia spectrum disorder (SSD) and (2) a subcohort (n = 50 625) randomly drawn from the source population. Data were analyzed from January 2021 to August 2023. Exposures: Carrier status of deletions and duplications at 27 autosomal rCNV loci was determined from neonatal blood samples genotyped on single-nucleotide variant microarrays. Main Outcomes and Measures: Population-based rCNV prevalence was estimated with a survey model using finite population correction to account for oversampling of cases. Hazard ratio (HR) estimates and 95% CIs for psychiatric disorders were derived using weighted Cox proportional hazard models. Risks were compared across outcomes, dosage type, and locus features using generalized estimating equation models. Results: A total of 3547 rCNVs were identified in 64 735 individuals assigned male at birth (53.8%) and 55 512 individuals assigned female at birth (46.2%) whose age at the end of follow-up ranged from 7.0 to 34.7 years (mean, 21.8 years). Most observed increases in rCNV-associated risk for ADHD, ASD, or SSD were moderate, and risk estimates were highly correlated across these disorders. Notable exceptions included high ASD-associated risk observed for Prader-Willi/Angelman syndrome duplications (HR, 20.8; 95% CI, 7.9-55). No rCNV was associated with increased MDD risk. Also, rCNV-associated risk was positively correlated with locus size and gene constraint but not with dosage type. Comparison with published case-control and community-based studies revealed a higher prevalence of deletions and lower associated increase in risk for several rCNVs in iPSYCH2015. Conclusions and Relevance: This study found that several rCNVs were more prevalent and conferred less risk of psychiatric disorders than estimated previously. Most case-control studies overestimate rCNV-associated risk of psychiatric disorders, likely because of selection bias. In an era where genetics is increasingly being clinically applied, these results highlight the importance of population-based risk estimates for genetics-based predictions.