The rate and nature of mitochondrial DNA mutations in human pedigrees.

We examined the rate and nature of mitochondrial DNA (mtDNA) mutations in humans using sequence data from 64,806 contemporary Icelanders from 2,548 matrilines. Based on 116,663 mother-child transmissions, 8,199 mutations were detected, providing robust rate estimates by nucleotide type, functional impact, position, and different alleles at the same position. We thoroughly document the true extent of hypermutability in mtDNA, mainly affecting the control region but also some coding-region variants. The results reveal the impact of negative selection on viable deleterious mutations, including rapidly mutating disease-associated 3243A>G and 1555A>G and pre-natal selection that most likely occurs during the development of oocytes. Finally, we show that the fate of new mutations is determined by a drastic germline bottleneck, amounting to an average of 3 mtDNA units effectively transmitted from mother to child.

Complex effects of sequence variants on lipid levels and coronary artery disease.

Many sequence variants have additive effects on blood lipid levels and, through that, on the risk of coronary artery disease (CAD). We show that variants also have non-additive effects and interact to affect lipid levels as well as affecting variance and correlations. Variance and correlation effects are often signatures of epistasis or gene-environmental interactions. These complex effects can translate into CAD risk. For example, Trp154Ter in FUT2 protects against CAD among subjects with the A1 blood group, whereas it associates with greater risk of CAD in others. His48Arg in ADH1B interacts with alcohol consumption to affect lipid levels and CAD. The effect of variants in TM6SF2 on blood lipids is greatest among those who never eat oily fish but absent from those who often do. This work demonstrates that variants that affect variance of quantitative traits can allow for the discovery of epistasis and interactions of variants with the environment.

The genetic history of Scandinavia from the Roman Iron Age to the present.

We investigate a 2,000-year genetic transect through Scandinavia spanning the Iron Age to the present, based on 48 new and 249 published ancient genomes and genotypes from 16,638 modern individuals. We find regional variation in the timing and magnitude of gene flow from three sources: the eastern Baltic, the British-Irish Isles, and southern Europe. British-Irish ancestry was widespread in Scandinavia from the Viking period, whereas eastern Baltic ancestry is more localized to Gotland and central Sweden. In some regions, a drop in current levels of external ancestry suggests that ancient immigrants contributed proportionately less to the modern Scandinavian gene pool than indicated by the ancestry of genomes from the Viking and Medieval periods. Finally, we show that a north-south genetic cline that characterizes modern Scandinavians is mainly due to the differential levels of Uralic ancestry and that this cline existed in the Viking Age and possibly earlier.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations.

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Genetic links between ovarian ageing, cancer risk and de novo mutation rates.

Human genetic studies of common variants have provided substantial insight into the biological mechanisms that govern ovarian ageing1. Here we report analyses of rare protein-coding variants in 106,973 women from the UK Biobank study, implicating genes with effects around five times larger than previously found for common variants (ETAA1, ZNF518A, PNPLA8, PALB2 and SAMHD1). The SAMHD1 association reinforces the link between ovarian ageing and cancer susceptibility1, with damaging germline variants being associated with extended reproductive lifespan and increased all-cause cancer risk in both men and women. Protein-truncating variants in ZNF518A are associated with shorter reproductive lifespan-that is, earlier age at menopause (by 5.61 years) and later age at menarche (by 0.56 years). Finally, using 8,089 sequenced trios from the 100,000 Genomes Project (100kGP), we observe that common genetic variants associated with earlier ovarian ageing associate with an increased rate of maternally derived de novo mutations. Although we were unable to replicate the finding in independent samples from the deCODE study, it is consistent with the expected role of DNA damage response genes in maintaining the genetic integrity of germ cells. This study provides evidence of genetic links between age of menopause and cancer risk.

Large-scale plasma proteomics comparisons through genetics and disease associations.

High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people2, for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity.

A polymorphism in IRF4 affects human pigmentation through a tyrosinase-dependent MITF/TFAP2A pathway.

Sequence polymorphisms linked to human diseases and phenotypes in genome-wide association studies often affect noncoding regions. A SNP within an intron of the gene encoding Interferon Regulatory Factor 4 (IRF4), a transcription factor with no known role in melanocyte biology, is strongly associated with sensitivity of skin to sun exposure, freckles, blue eyes, and brown hair color. Here, we demonstrate that this SNP lies within an enhancer of IRF4 transcription in melanocytes. The allele associated with this pigmentation phenotype impairs binding of the TFAP2A transcription factor that, together with the melanocyte master regulator MITF, regulates activity of the enhancer. Assays in zebrafish and mice reveal that IRF4 cooperates with MITF to activate expression of Tyrosinase (TYR), an essential enzyme in melanin synthesis. Our findings provide a clear example of a noncoding polymorphism that affects a phenotype by modulating a developmental gene regulatory network.

PSMD11 loss-of-function variants correlate with a neurobehavioral phenotype, obesity, and increased interferon response.

Primary proteasomopathies have recently emerged as a new class of rare early-onset neurodevelopmental disorders (NDDs) caused by pathogenic variants in the PSMB1, PSMC1, PSMC3, or PSMD12 proteasome genes. Proteasomes are large multi-subunit protein complexes that maintain cellular protein homeostasis by clearing ubiquitin-tagged damaged, misfolded, or unnecessary proteins. In this study, we have identified PSMD11 as an additional proteasome gene in which pathogenic variation is associated with an NDD-causing proteasomopathy. PSMD11 loss-of-function variants caused early-onset syndromic intellectual disability and neurodevelopmental delay with recurrent obesity in 10 unrelated children. Our findings demonstrate that the cognitive impairment observed in these individuals could be recapitulated in Drosophila melanogaster with depletion of the PMSD11 ortholog Rpn6, which compromised reversal learning. Our investigations in subject samples further revealed that PSMD11 loss of function resulted in impaired 26S proteasome assembly and the acquisition of a persistent type I interferon (IFN) gene signature, mediated by the integrated stress response (ISR) protein kinase R (PKR). In summary, these data identify PSMD11 as an additional member of the growing family of genes associated with neurodevelopmental proteasomopathies and provide insights into proteasomal biology in human health.

Actionable Genotypes and Their Association with Life Span in Iceland.

BACKGROUND: In 2021, the American College of Medical Genetics and Genomics (ACMG) recommended reporting actionable genotypes in 73 genes associated with diseases for which preventive or therapeutic measures are available. Evaluations of the association of actionable genotypes in these genes with life span are currently lacking. METHODS: We assessed the prevalence of coding and splice variants in genes on the ACMG Secondary Findings, version 3.0 (ACMG SF v3.0), list in the genomes of 57,933 Icelanders. We assigned pathogenicity to all reviewed variants using reported evidence in the ClinVar database, the frequency of variants, and their associations with disease to create a manually curated set of actionable genotypes (variants). We assessed the relationship between these genotypes and life span and further examined the specific causes of death among carriers. RESULTS: Through manual curation of 4405 sequence variants in the ACMG SF v3.0 genes, we identified 235 actionable genotypes in 53 genes. Of the 57,933 participants, 2306 (4.0%) carried at least one actionable genotype. We found shorter median survival among persons carrying actionable genotypes than among noncarriers. Specifically, we found that carrying an actionable genotype in a cancer gene was associated with survival that was 3 years shorter than that among noncarriers, with causes of death among carriers attributed primarily to cancer-related conditions. Furthermore, we found evidence of association between carrying an actionable genotype in certain genes in the cardiovascular disease group and a reduced life span. CONCLUSIONS: On the basis of the ACMG SF v3.0 guidelines, we found that approximately 1 in 25 Icelanders carried an actionable genotype and that carrying such a genotype was associated with a reduced life span. (Funded by deCODE Genetics-Amgen.).

The nature of Neanderthal introgression revealed by 27,566 Icelandic genomes.

Human evolutionary history is rich with the interbreeding of divergent populations. Most humans outside of Africa trace about 2% of their genomes to admixture from Neanderthals, which occurred 50-60 thousand years ago1. Here we examine the effect of this event using 14.4 million putative archaic chromosome fragments that were detected in fully phased whole-genome sequences from 27,566 Icelanders, corresponding to a range of 56,388-112,709 unique archaic fragments that cover 38.0-48.2% of the callable genome. On the basis of the similarity with known archaic genomes, we assign 84.5% of fragments to an Altai or Vindija Neanderthal origin and 3.3% to Denisovan origin; 12.2% of fragments are of unknown origin. We find that Icelanders have more Denisovan-like fragments than expected through incomplete lineage sorting. This is best explained by Denisovan gene flow, either into ancestors of the introgressing Neanderthals or directly into humans. A within-individual, paired comparison of archaic fragments with syntenic non-archaic fragments revealed that, although the overall rate of mutation was similar in humans and Neanderthals during the 500 thousand years that their lineages were separate, there were differences in the relative frequencies of mutation types-perhaps due to different generation intervals for males and females. Finally, we assessed 271 phenotypes, report 5 associations driven by variants in archaic fragments and show that the majority of previously reported associations are better explained by non-archaic variants.

FLT3 stop mutation increases FLT3 ligand level and risk of autoimmune thyroid disease.

Autoimmune thyroid disease is the most common autoimmune disease and is highly heritable1. Here, by using a genome-wide association study of 30,234 cases and 725,172 controls from Iceland and the UK Biobank, we find 99 sequence variants at 93 loci, of which 84 variants are previously unreported2-7. A low-frequency (1.36%) intronic variant in FLT3 (rs76428106-C) has the largest effect on risk of autoimmune thyroid disease (odds ratio (OR) = 1.46, P = 2.37 × 10-24). rs76428106-C is also associated with systemic lupus erythematosus (OR = 1.90, P = 6.46 × 10-4), rheumatoid factor and/or anti-CCP-positive rheumatoid arthritis (OR = 1.41, P = 4.31 × 10-4) and coeliac disease (OR = 1.62, P = 1.20 × 10-4). FLT3 encodes fms-related tyrosine kinase 3, a receptor that regulates haematopoietic progenitor and dendritic cells. RNA sequencing revealed that rs76428106-C generates a cryptic splice site, which introduces a stop codon in 30% of transcripts that are predicted to encode a truncated protein, which lacks its tyrosine kinase domains. Each copy of rs76428106-C doubles the plasma levels of the FTL3 ligand. Activating somatic mutations in FLT3 are associated with acute myeloid leukaemia8 with a poor prognosis and rs76428106-C also predisposes individuals to acute myeloid leukaemia (OR = 1.90, P = 5.40 × 10-3). Thus, a predicted loss-of-function germline mutation in FLT3 causes a reduction in full-length FLT3, with a compensatory increase in the levels of its ligand and an increased disease risk, similar to that of a gain-of-function mutation.

Monoallelic de novo variants in DDX17 cause a neurodevelopmental disorder.

DDX17 is an RNA helicase shown to be involved in critical processes during the early phases of neuronal differentiation. Globally, we compiled a case-series of 11 patients with neurodevelopmental phenotypes harbouring de novo monoallelic variants in DDX17. All 11 patients in our case series had a neurodevelopmental phenotype, whereby intellectual disability, delayed speech and language, and motor delay predominated. We performed in utero cortical electroporation in the brain of developing mice, assessing axon complexity and outgrowth of electroporated neurons, comparing wild-type and Ddx17 knockdown. We then undertook ex vivo cortical electroporation on neuronal progenitors to quantitatively assess axonal development at a single cell resolution. Mosaic ddx17 crispants and heterozygous knockouts in Xenopus tropicalis were generated for assessment of morphology, behavioural assays, and neuronal outgrowth measurements. We further undertook transcriptomic analysis of neuroblastoma SH-SY5Y cells, to identify differentially expressed genes in DDX17-KD cells compared to controls. Knockdown of Ddx17 in electroporated mouse neurons in vivo showed delayed neuronal migration as well as decreased cortical axon complexity. Mouse primary cortical neurons revealed reduced axon outgrowth upon knockdown of Ddx17 in vitro. The axon outgrowth phenotype was replicated in crispant ddx17 tadpoles and in heterozygotes. Heterozygous tadpoles had clear neurodevelopmental defects and showed an impaired neurobehavioral phenotype. Transcriptomic analysis identified a statistically significant number of differentially expressed genes involved in neurodevelopmental processes in DDX17-KD cells compared to control cells. We have identified potential neurodevelopment disease-causing variants in a gene not previously associated with genetic disease, DDX17. We provide evidence for the role of the gene in neurodevelopment in both mammalian and non-mammalian species and in controlling the expression of key neurodevelopment genes.

Genetic predisposition to mosaic Y chromosome loss in blood.

Mosaic loss of chromosome Y (LOY) in circulating white blood cells is the most common form of clonal mosaicism1-5, yet our knowledge of the causes and consequences of this is limited. Here, using a computational approach, we estimate that 20% of the male population represented in the UK Biobank study (n = 205,011) has detectable LOY. We identify 156 autosomal genetic determinants of LOY, which we replicate in 757,114 men of European and Japanese ancestry. These loci highlight genes that are involved in cell-cycle regulation and cancer susceptibility, as well as somatic drivers of tumour growth and targets of cancer therapy. We demonstrate that genetic susceptibility to LOY is associated with non-haematological effects on health in both men and women, which supports the hypothesis that clonal haematopoiesis is a biomarker of genomic instability in other tissues. Single-cell RNA sequencing identifies dysregulated expression of autosomal genes in leukocytes with LOY and provides insights into why clonal expansion of these cells may occur. Collectively, these data highlight the value of studying clonal mosaicism to uncover fundamental mechanisms that underlie cancer and other ageing-related diseases.

A partial loss-of-function variant in STAT6 protects against T2 asthma.

BACKGROUND: Signal Transducer and Activator of Transcription 6 (STAT6) is central to Type 2 (T2) inflammation and common non-coding variants at the STAT6 locus associate with various T2 inflammatory traits, including diseases, and its pathway is widely targeted in asthma treatment. OBJECTIVE: To test the association of a rare missense variant in STAT6, p.L406P, with T2 inflammatory traits, including the risk of asthma and allergic diseases, and to characterize its functional consequences in cell culture. METHODS: We tested association of p.L406P with plasma protein levels, white blood cell counts and the risk of asthma and allergic phenotypes. We tested significant associations in other cohorts using a burden test. The effects of p.L406P on STAT6 protein function were examined in cell lines and by comparing CD4+ T-cell responses from carriers and non-carriers of the variant. RESULTS: p.L406P associated with reduced plasma levels of STAT6 and IgE as well as with lower eosinophil and basophil counts in blood. It also protected against asthma, mostly driven by severe T2 high asthma. We showed that p.L406P led to lower IL-4-induced activation in luciferase reporter assays and lower levels of STAT6 in CD4+ T cells. We identified multiple genes with expression that was affected by the p.L406P genotype upon IL-4 treatment of CD4+ T cells; the effect was consistent with a weaker IL-4 response in carriers than non-carriers of p.L406P. CONCLUSIONS: We report a partial loss-of-function variant in STAT6, resulting in dampened IL-4 responses and protection from T2 high asthma, implicating STAT6 as an attractive therapeutic target.

Longitudinal metabolite and protein trajectories prior to diabetes mellitus diagnosis in Danish blood donors: a nested case-control study.

AIMS/HYPOTHESIS: Metabolic risk factors and plasma biomarkers for diabetes have previously been shown to change prior to a clinical diabetes diagnosis. However, these markers only cover a small subset of molecular biomarkers linked to the disease. In this study, we aimed to profile a more comprehensive set of molecular biomarkers and explore their temporal association with incident diabetes. METHODS: We performed a targeted analysis of 54 proteins and 171 metabolites and lipoprotein particles measured in three sequential samples spanning up to 11 years of follow-up in 324 individuals with incident diabetes and 359 individuals without diabetes in the Danish Blood Donor Study (DBDS) matched for sex and birth year distribution. We used linear mixed-effects models to identify temporal changes before a diabetes diagnosis, either for any incident diabetes diagnosis or for type 1 and type 2 diabetes mellitus diagnoses specifically. We further performed linear and non-linear feature selection, adding 28 polygenic risk scores to the biomarker pool. We tested the time-to-event prediction gain of the biomarkers with the highest variable importance, compared with selected clinical covariates and plasma glucose. RESULTS: We identified two proteins and 16 metabolites and lipoprotein particles whose levels changed temporally before diabetes diagnosis and for which the estimated marginal means were significant after FDR adjustment. Sixteen of these have not previously been described. Additionally, 75 biomarkers were consistently higher or lower in the years before a diabetes diagnosis. We identified a single temporal biomarker for type 1 diabetes, IL-17A/F, a cytokine that is associated with multiple other autoimmune diseases. Inclusion of 12 biomarkers improved the 10-year prediction of a diabetes diagnosis (i.e. the area under the receiver operating curve increased from 0.79 to 0.84), compared with clinical information and plasma glucose alone. CONCLUSIONS/INTERPRETATION: Systemic molecular changes manifest in plasma several years before a diabetes diagnosis. A particular subset of biomarkers shows distinct, time-dependent patterns, offering potential as predictive markers for diabetes onset. Notably, these biomarkers show shared and distinct patterns between type 1 diabetes and type 2 diabetes. After independent replication, our findings may be used to develop new clinical prediction models.

Exploring the Genetic Risk of Childhood Daytime Urinary Incontinence: A Genome-Wide Association Study.

PURPOSE: Childhood incontinence is stigmatized and underprioritized, and a basic understanding of its pathogenesis is missing. Our goal was to identify risk-conferring genetic variants in daytime urinary incontinence (DUI). MATERIALS AND METHODS: We conducted a genome-wide association study in the Danish iPSYCH2015 cohort. Cases (3024) were identified through DUI diagnosis codes and redeemed prescriptions for DUI medication in patients aged 5 to 20 years. Controls (30,240), selected from the same sample, were matched to cases on sex and psychiatric diagnoses, if any, and down-sampled to a 1:10 case:control ratio. Replication was performed in the Icelandic deCODE cohort (5475 cases/287,773 controls). Single-nucleotide polymorphism heritability was calculated using the genome-based restricted maximum likelihood method. Cross-trait genetic correlation was estimated using linkage disequilibrium score regression. Polygenic risk scores generated with LDpred2-auto and BOLT-LMM were assessed for association. RESULTS: Variants on chromosome 6 (rs12210989, odds ratio [OR] 1.24, 95% CI 1.17-1.32, P = 3.21 × 10-12) and 20 (rs4809801, OR 1.18, 95% CI 1.11-1.25, P = 3.66 × 10-8) reached genome-wide significance and implicated the PRDM13 and RIPOR3 genes. Chromosome 6 findings were replicated (P = .024, OR 1.09, 95% CI 1.01-1.16). Liability scale heritability ranged from 10.20% (95% CI 6.40%-14.00%) to 15.30% (95% CI 9.66%-20.94%). DUI and nocturnal enuresis showed positive genetic correlation (rg = 1.28 ± 0.38, P = .0007). DUI was associated with attention-deficit/hyperactivity disorder (OR 1.098, 95% CI 1.046-1.152, P < .0001) and BMI (OR 1.129, 95% CI 1.081-1.178, P < .0001) polygenic risk. CONCLUSIONS: Common genetic variants contribute to the risk of childhood DUI, and genes important in neuronal development and detrusor smooth muscle activity were implicated. These findings may help guide identification of new treatment targets.

Genome-wide association study on 13 167 individuals identifies regulators of blood CD34+cell levels.

Stem cell transplantation is a cornerstone in the treatment of blood malignancies. The most common method to harvest stem cells for transplantation is by leukapheresis, requiring mobilization of CD34+ hematopoietic stem and progenitor cells (HSPCs) from the bone marrow into the blood. Identifying the genetic factors that control blood CD34+ cell levels could reveal new drug targets for HSPC mobilization. Here we report the first large-scale, genome-wide association study on blood CD34+ cell levels. Across 13 167 individuals, we identify 9 significant and 2 suggestive associations, accounted for by 8 loci (PPM1H, CXCR4, ENO1-RERE, ITGA9, ARHGAP45, CEBPA, TERT, and MYC). Notably, 4 of the identified associations map to CXCR4, showing that bona fide regulators of blood CD34+ cell levels can be identified through genetic variation. Further, the most significant association maps to PPM1H, encoding a serine/threonine phosphatase never previously implicated in HSPC biology. PPM1H is expressed in HSPCs, and the allele that confers higher blood CD34+ cell levels downregulates PPM1H. Through functional fine-mapping, we find that this downregulation is caused by the variant rs772557-A, which abrogates an MYB transcription factor-binding site in PPM1H intron 1 that is active in specific HSPC subpopulations, including hematopoietic stem cells, and interacts with the promoter by chromatin looping. Furthermore, PPM1H knockdown increases the proportion of CD34+ and CD34+90+ cells in cord blood assays. Our results provide the first large-scale analysis of the genetic architecture of blood CD34+ cell levels and warrant further investigation of PPM1H as a potential inhibition target for stem cell mobilization.

Genetic variants associated with syncope implicate neural and autonomic processes.

AIMS: Syncope is a common and clinically challenging condition. In this study, the genetics of syncope were investigated to seek knowledge about its pathophysiology and prognostic implications. METHODS AND RESULTS: This genome-wide association meta-analysis included 56 071 syncope cases and 890 790 controls from deCODE genetics (Iceland), UK Biobank (United Kingdom), and Copenhagen Hospital Biobank Cardiovascular Study/Danish Blood Donor Study (Denmark), with a follow-up assessment of variants in 22 412 cases and 286 003 controls from Intermountain (Utah, USA) and FinnGen (Finland). The study yielded 18 independent syncope variants, 17 of which were novel. One of the variants, p.Ser140Thr in PTPRN2, affected syncope only when maternally inherited. Another variant associated with a vasovagal reaction during blood donation and five others with heart rate and/or blood pressure regulation, with variable directions of effects. None of the 18 associations could be attributed to cardiovascular or other disorders. Annotation with regard to regulatory elements indicated that the syncope variants were preferentially located in neural-specific regulatory regions. Mendelian randomization analysis supported a causal effect of coronary artery disease on syncope. A polygenic score (PGS) for syncope captured genetic correlation with cardiovascular disorders, diabetes, depression, and shortened lifespan. However, a score based solely on the 18 syncope variants performed similarly to the PGS in detecting syncope risk but did not associate with other disorders. CONCLUSION: The results demonstrate that syncope has a distinct genetic architecture that implicates neural regulatory processes and a complex relationship with heart rate and blood pressure regulation. A shared genetic background with poor cardiovascular health was observed, supporting the importance of a thorough assessment of individuals presenting with syncope.