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.

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.).

Humoral Immune Response to SARS-CoV-2 in Iceland.

BACKGROUND: Little is known about the nature and durability of the humoral immune response to infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). METHODS: We measured antibodies in serum samples from 30,576 persons in Iceland, using six assays (including two pan-immunoglobulin [pan-Ig] assays), and we determined that the appropriate measure of seropositivity was a positive result with both pan-Ig assays. We tested 2102 samples collected from 1237 persons up to 4 months after diagnosis by a quantitative polymerase-chain-reaction (qPCR) assay. We measured antibodies in 4222 quarantined persons who had been exposed to SARS-CoV-2 and in 23,452 persons not known to have been exposed. RESULTS: Of the 1797 persons who had recovered from SARS-CoV-2 infection, 1107 of the 1215 who were tested (91.1%) were seropositive; antiviral antibody titers assayed by two pan-Ig assays increased during 2 months after diagnosis by qPCR and remained on a plateau for the remainder of the study. Of quarantined persons, 2.3% were seropositive; of those with unknown exposure, 0.3% were positive. We estimate that 0.9% of Icelanders were infected with SARS-CoV-2 and that the infection was fatal in 0.3%. We also estimate that 56% of all SARS-CoV-2 infections in Iceland had been diagnosed with qPCR, 14% had occurred in quarantined persons who had not been tested with qPCR (or who had not received a positive result, if tested), and 30% had occurred in persons outside quarantine and not tested with qPCR. CONCLUSIONS: Our results indicate that antiviral antibodies against SARS-CoV-2 did not decline within 4 months after diagnosis. We estimate that the risk of death from infection was 0.3% and that 44% of persons infected with SARS-CoV-2 in Iceland were not diagnosed by qPCR.

Parental influence on human germline de novo mutations in 1,548 trios from Iceland.

The characterization of mutational processes that generate sequence diversity in the human genome is of paramount importance both to medical genetics and to evolutionary studies. To understand how the age and sex of transmitting parents affect de novo mutations, here we sequence 1,548 Icelanders, their parents, and, for a subset of 225, at least one child, to 35× genome-wide coverage. We find 108,778 de novo mutations, both single nucleotide polymorphisms and indels, and determine the parent of origin of 42,961. The number of de novo mutations from mothers increases by 0.37 per year of age (95% CI 0.32-0.43), a quarter of the 1.51 per year from fathers (95% CI 1.45-1.57). The number of clustered mutations increases faster with the mother's age than with the father's, and the genomic span of maternal de novo mutation clusters is greater than that of paternal ones. The types of de novo mutation from mothers change substantially with age, with a 0.26% (95% CI 0.19-0.33%) decrease in cytosine-phosphate-guanine to thymine-phosphate-guanine (CpG>TpG) de novo mutations and a 0.33% (95% CI 0.28-0.38%) increase in C>G de novo mutations per year, respectively. Remarkably, these age-related changes are not distributed uniformly across the genome. A striking example is a 20 megabase region on chromosome 8p, with a maternal C>G mutation rate that is up to 50-fold greater than the rest of the genome. The age-related accumulation of maternal non-crossover gene conversions also mostly occurs within these regions. Increased sequence diversity and linkage disequilibrium of C>G variants within regions affected by excess maternal mutations indicate that the underlying mutational process has persisted in humans for thousands of years. Moreover, the regional excess of C>G variation in humans is largely shared by chimpanzees, less by gorillas, and is almost absent from orangutans. This demonstrates that sequence diversity in humans results from evolving interactions between age, sex, mutation type, and genomic location.

Sequence variants associating with urinary biomarkers.

Urine dipstick tests are widely used in routine medical care to diagnose kidney and urinary tract and metabolic diseases. Several environmental factors are known to affect the test results, whereas the effects of genetic diversity are largely unknown. We tested 32.5 million sequence variants for association with urinary biomarkers in a set of 150 274 Icelanders with urine dipstick measurements. We detected 20 association signals, of which 14 are novel, associating with at least one of five clinical entities defined by the urine dipstick: glucosuria, ketonuria, proteinuria, hematuria and urine pH. These include three independent glucosuria variants at SLC5A2, the gene encoding the sodium-dependent glucose transporter (SGLT2), a protein targeted pharmacologically to increase urinary glucose excretion in the treatment of diabetes. Two variants associating with proteinuria are in LRP2 and CUBN, encoding the co-transporters megalin and cubilin, respectively, that mediate proximal tubule protein uptake. One of the hematuria-associated variants is a rare, previously unreported 2.5 kb exonic deletion in COL4A3. Of the four signals associated with urine pH, we note that the pH-increasing alleles of two variants (POU2AF1, WDR72) associate significantly with increased risk of kidney stones. Our results reveal that genetic factors affect variability in urinary biomarkers, in both a disease dependent and independent context.

popSTR2 enables clinical and population-scale genotyping of microsatellites.

SUMMARY: popSTR2 is an update and augmentation of our previous work 'popSTR: a population-based microsatellite genotyper'. To make genotyping sensitive to inter-sample differences, we supply a kernel to estimate sample-specific slippage rates. For clinical sequencing purposes, a panel of known pathogenic repeat expansions is provided along with a script that scans and flags for manual inspection markers indicative of a pathogenic expansion. Like its predecessor, popSTR2 allows for joint genotyping of samples at a population scale. We now provide a binning method that makes the microsatellite genotypes more amenable to analysis within standard association pipelines and can increase association power. AVAILABILITY AND IMPLEMENTATION: https://github.com/DecodeGenetics/popSTR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Rare SCARB1 mutations associate with high-density lipoprotein cholesterol but not with coronary artery disease.

Aims: Scavenger receptor Class B Type 1 (SR-BI) is a major receptor for high-density lipoprotein (HDL) that promotes hepatic uptake of cholesterol from HDL. A rare mutation p.P376L, in the gene encoding SR-BI, SCARB1, was recently reported to associate with elevated HDL cholesterol (HDL-C) and increased risk of coronary artery disease (CAD), suggesting that increased HDL-C caused by SR-BI impairment might be an independent marker of cardiovascular risk. We tested the hypothesis that alleles in or close to SCARB1 that associate with elevated levels of HDL-C also associate with increased risk of CAD in the relatively homogeneous population of Iceland. Methods and results: Using a large resource of whole-genome sequenced Icelanders, we identified thirteen SCARB1 coding mutations that we examined for association with HDL-C (n = 136 672). Three rare SCARB1 mutations, encoding p.G319V, p.V111M, and p.V32M (combined allelic frequency = 0.2%) associate with elevated levels of HDL-C (p.G319V: ß = 11.1 mg/dL, P = 8.0 × 10-7; p.V111M: ß = 8.3 mg/dL, P = 1.1 × 10-6; p.V32M: ß = 10.2 mg/dL, P = 8.1 × 10-4). These mutations do not associate with CAD (36 886 cases/306 268 controls) (odds ratio = 0.90, 95% confidence interval 0.67-1.22, P = 0.49), despite effects on HDL-C comparable to that reported for p.P376L, both in terms of direction and magnitude. Furthermore, HDL-C raising alleles of three common SCARB1 non-coding variants, including one previously unreported (rs61941676-C: ß = 1.25 mg/dL, P = 1.7 × 10-18), and of one low frequency coding variant (p.V135I) that independently associate with higher HDL-C, do not confer increased risk of CAD. Conclusion: Elevated HDL-C due to genetically compromised SR-BI function is not a marker of CAD risk.

Identification of Lynch syndrome risk variants in the Romanian population.

Two familial forms of colorectal cancer (CRC), Lynch syndrome (LS) and familial adenomatous polyposis (FAP), are caused by rare mutations in DNA mismatch repair genes (MLH1, MSH2, MSH6, PMS2) and the genes APC and MUTYH, respectively. No information is available on the presence of high-risk CRC mutations in the Romanian population. We performed whole-genome sequencing of 61 Romanian CRC cases with a family history of cancer and/or early onset of disease, focusing the analysis on candidate variants in the LS and FAP genes. The frequencies of all candidate variants were assessed in a cohort of 688 CRC cases and 4567 controls. Immunohistochemical (IHC) staining for MLH1, MSH2, MSH6, and PMS2 was performed on tumour tissue. We identified 11 candidate variants in 11 cases; six variants in MLH1, one in MSH6, one in PMS2, and three in APC. Combining information on the predicted impact of the variants on the proteins, IHC results and previous reports, we found three novel pathogenic variants (MLH1:p.Lys84ThrfsTer4, MLH1:p.Ala586CysfsTer7, PMS2:p.Arg211ThrfsTer38), and two novel variants that are unlikely to be pathogenic. Also, we confirmed three previously published pathogenic LS variants and suggest to reclassify a previously reported variant of uncertain significance to pathogenic (MLH1:c.1559-1G>C).

COPA syndrome in an Icelandic family caused by a recurrent missense mutation in COPA.

BACKGROUND: Rare missense mutations in the gene encoding coatomer subunit alpha (COPA) have recently been shown to cause autoimmune interstitial lung, joint and kidney disease, also known as COPA syndrome, under a dominant mode of inheritance. CASE PRESENTATION: Here we describe an Icelandic family with three affected individuals over two generations with a rare clinical presentation of lung and joint disease and a histological diagnosis of follicular bronchiolitis. We performed whole-genome sequencing (WGS) of the three affected as well as three unaffected members of the family, and searched for rare genotypes associated with disease using 30,067 sequenced Icelanders as a reference population. We assessed all coding and splicing variants, prioritizing variants in genes known to cause interstitial lung disease. We detected a heterozygous missense mutation, p.Glu241Lys, in the COPA gene, private to the affected family members. The mutation occurred de novo in the paternal germline of the index case and was absent from 30,067 Icelandic genomes and 141,353 individuals from the genome Aggregation Database (gnomAD). The mutation occurs within the conserved and functionally important WD40 domain of the COPA protein. CONCLUSIONS: This is the second report of the p.Glu241Lys mutation in COPA, indicating the recurrent nature of the mutation. The mutation was reported to co-segregate with COPA syndrome in a large family from the USA with five affected members, and classified as pathogenic. The two separate occurrences of the p.Glu241Lys mutation in cases and its absence from a large number of sequenced genomes confirms its role in the pathogenesis of the COPA syndrome.

Compound heterozygous mutations in UBA5 causing early-onset epileptic encephalopathy in two sisters.

BACKGROUND: Epileptic encephalopathies are a group of childhood epilepsies that display high phenotypic and genetic heterogeneity. The recent, extensive use of next-generation sequencing has identified a large number of genes in epileptic encephalopathies, including UBA5 in which biallelic mutations were first described as pathogenic in 2016 (Colin E et al., Am J Hum Genet 99(3):695-703, 2016. Muona M et al., Am J Hum Genet 99(3):683-694, 2016). UBA5 encodes an activating enzyme for a post-translational modification mechanism known as ufmylation, and is the first gene from the ufmylation pathway that is linked to disease. CASE PRESENTATION: We sequenced the genomes of two sisters with early-onset epileptic encephalopathy along with their unaffected parents in an attempt to find a genetic cause for their condition. The sisters, born in 2004 and 2006, presented with infantile spasms at six months of age, which later progressed to recurrent, treatment-resistant seizures. We detected a compound heterozygous genotype in UBA5 in the sisters, a genotype not seen elsewhere in an Icelandic reference set of 30,067 individuals nor in public databases. One of the mutations, c.684G > A, is a paternally inherited exonic splicing mutation, occuring at the last nucleotide of exon 7 of UBA5. The mutation is predicted to disrupt the splice site, resulting in loss-of-function of one allele of UBA5. The second mutation is a maternally inherited missense mutation, p.Ala371Thr, previously reported as pathogenic when in compound heterozygosity with a loss-of-function mutation in UBA5 and is believed to produce a hypomorphic allele. Supportive of this, we have identified three adult Icelanders homozygous for the p.Ala371Thr mutation who show no signs of neurological disease. CONCLUSIONS: We describe compound heterozygous mutations in the UBA5 gene in two sisters with early-onset epileptic encephalopathy. To our knowledge, this is the first description of mutations in UBA5 since the initial discovery that pathogenic biallelic variants in the gene cause early-onset epileptic encephalopathy. We further provide confirmatory evidence that p.Ala371Thr is a hypomorphic mutation, by presenting three adult homozygotes who show no signs of neurological disease.

A population-based survey of FBN1 variants in Iceland reveals underdiagnosis of Marfan syndrome.

Marfan syndrome (MFS) is an autosomal dominant condition characterized by aortic aneurysm, skeletal abnormalities, and lens dislocation, and is caused by variants in the FBN1 gene. To explore causes of MFS and the prevalence of the disease in Iceland we collected information from all living individuals with a clinical diagnosis of MFS in Iceland (n = 32) and performed whole-genome sequencing of those who did not have a confirmed genetic diagnosis (27/32). Moreover, to assess a potential underdiagnosis of MFS in Iceland we attempted a genotype-based approach to identify individuals with MFS. We interrogated deCODE genetics' database of 35,712 whole-genome sequenced individuals to search for rare sequence variants in FBN1. Overall, we identified 15 pathogenic or likely pathogenic variants in FBN1 in 44 individuals, only 22 of whom were previously diagnosed with MFS. The most common of these variants, NM_000138.4:c.8038 C > T p.(Arg2680Cys), is present in a multi-generational pedigree, and was found to stem from a single forefather born around 1840. The p.(Arg2680Cys) variant associates with a form of MFS that seems to have an enrichment of abdominal aortic aneurysm, suggesting that this may be a particularly common feature of p.(Arg2680Cys)-associated MFS. Based on these combined genetic and clinical data, we show that MFS prevalence in Iceland could be as high as 1/6,600 in Iceland, compared to 1/10,000 based on clinical diagnosis alone, which indicates underdiagnosis of this actionable genetic disorder.

Variant in the synaptonemal complex protein SYCE2 associates with pregnancy loss through effect on recombination.

Two-thirds of all human conceptions are lost, in most cases before clinical detection. The lack of detailed understanding of the causes of pregnancy losses constrains focused counseling for future pregnancies. We have previously shown that a missense variant in synaptonemal complex central element protein 2 (SYCE2), in a key residue for the assembly of the synaptonemal complex backbone, associates with recombination traits. Here we show that it also increases risk of pregnancy loss in a genome-wide association analysis on 114,761 women with reported pregnancy loss. We further show that the variant associates with more random placement of crossovers and lower recombination rate in longer chromosomes but higher in the shorter ones. These results support the hypothesis that some pregnancy losses are due to failures in recombination. They further demonstrate that variants with a substantial effect on the quality of recombination can be maintained in the population.

Screening for Rare Coding Variants That Associate With the QTc Interval in Iceland.

Background Long-QT syndrome (LQTS) is a cardiac repolarization abnormality that can lead to sudden cardiac death. The most common causes are rare coding variants in the genes KCNQ1, KCNH2, and SCN5A. The data on LQTS epidemiology are limited, and information on expressivity and penetrance of pathogenic variants is sparse. Methods and Results We screened for rare coding variants associated with the corrected QT (QTc) interval in Iceland. We explored the frequency of the identified variants, their penetrance, and their association with severe events. Twelve variants were associated with the QTc interval. Five in KCNQ1, 3 in KCNH2, 2 in cardiomyopathy genes MYBPC3 and PKP2, and 2 in genes where coding variants have not been associated with the QTc interval, ISOC1 and MYOM2. The combined carrier frequency of the 8 variants in the previously known LQTS genes was 530 per 100 000 individuals (1:190). p.Tyr315Cys and p.Leu273Phe in KCNQ1 were associated with having a mean QTc interval longer than 500 ms (P=4.2×10-7; odds ratio [OR], 38.6; P=8.4×10-10, OR, 26.5; respectively), and p.Leu273Phe was associated with sudden cardiac death (P=0.0034; OR, 2.99). p.Val215Met in KCNQ1 was carried by 1 in 280 Icelanders, had a smaller effect on the QTc interval (P=1.8×10-44; effect, 22.8 ms), and did not associate with severe clinical events. Conclusions The carrier frequency of associating variants in LQTS genes was higher than previous estimates of the prevalence of LQTS. The variants have variable effects on the QTc interval, and carriers of p.Tyr315Cys and p.Leu273Phe have a more severe disease than carriers of p.Val215Met. These data could lead to improved identification, risk stratification, and a more precise clinical approach to those with QTc prolongation.

Sequence variants affecting the genome-wide rate of germline microsatellite mutations.

Microsatellites are polymorphic tracts of short tandem repeats with one to six base-pair (bp) motifs and are some of the most polymorphic variants in the genome. Using 6084 Icelandic parent-offspring trios we estimate 63.7 (95% CI: 61.9-65.4) microsatellite de novo mutations (mDNMs) per offspring per generation, excluding one bp repeats motifs (homopolymers) the estimate is 48.2 mDNMs (95% CI: 46.7-49.6). Paternal mDNMs occur at longer repeats than maternal ones, which are in turn larger with a mean size of 3.4 bp vs 3.1 bp for paternal ones. mDNMs increase by 0.97 (95% CI: 0.90-1.04) and 0.31 (95% CI: 0.25-0.37) per year of father's and mother's age at conception, respectively. Here, we find two independent coding variants that associate with the number of mDNMs transmitted to offspring; The minor allele of a missense variant (allele frequency (AF) = 1.9%) in MSH2, a mismatch repair gene, increases transmitted mDNMs from both parents (effect: 13.1 paternal and 7.8 maternal mDNMs). A synonymous variant (AF = 20.3%) in NEIL2, a DNA damage repair gene, increases paternally transmitted mDNMs (effect: 4.4 mDNMs). Thus, the microsatellite mutation rate in humans is in part under genetic control.

Sequence variant affects GCSAML splicing, mast cell specific proteins, and risk of urticaria.

Urticaria is a skin disorder characterized by outbreaks of raised pruritic wheals. In order to identify sequence variants associated with urticaria, we performed a meta-analysis of genome-wide association studies for urticaria with a total of 40,694 cases and 1,230,001 controls from Iceland, the UK, Finland, and Japan. We also performed transcriptome- and proteome-wide analyses in Iceland and the UK. We found nine sequence variants at nine loci associating with urticaria. The variants are at genes participating in type 2 immune responses and/or mast cell biology (CBLB, FCER1A, GCSAML, STAT6, TPSD1, ZFPM1), the innate immunity (C4), and NF-κB signaling. The most significant association was observed for the splice-donor variant rs56043070[A] (hg38: chr1:247556467) in GCSAML (MAF = 6.6%, OR = 1.24 (95%CI: 1.20-1.28), P-value = 3.6 × 10-44). We assessed the effects of the variants on transcripts, and levels of proteins relevant to urticaria pathophysiology. Our results emphasize the role of type 2 immune response and mast cell activation in the pathogenesis of urticaria. Our findings may point to an IgE-independent urticaria pathway that could help address unmet clinical need.

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality.

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

Genetic architecture of band neutrophil fraction in Iceland.

The characteristic lobulated nuclear morphology of granulocytes is partially determined by composition of nuclear envelope proteins. Abnormal nuclear morphology is primarily observed as an increased number of hypolobulated immature neutrophils, called band cells, during infection or in rare envelopathies like Pelger-Huët anomaly. To search for sequence variants affecting nuclear morphology of granulocytes, we performed a genome-wide association study using band neutrophil fraction from 88,101 Icelanders. We describe 13 sequence variants affecting band neutrophil fraction at nine loci. Five of the variants are at the Lamin B receptor (LBR) locus, encoding an inner nuclear membrane protein. Mutations in LBR are linked to Pelger-Huët anomaly. In addition, we identify cosegregation of a rare stop-gain sequence variant in LBR and Pelger Huët anomaly in an Icelandic eight generation pedigree, initially reported in 1963. Two of the other loci include genes which, like LBR, play a role in the nuclear membrane function and integrity. These GWAS results highlight the role proteins of the inner nuclear membrane have as important for neutrophil nuclear morphology.

Population-level deficit of homozygosity unveils CPSF3 as an intellectual disability syndrome gene.

Predicting the pathogenicity of biallelic missense variants can be challenging. Here, we use a deficit of observed homozygous carriers of missense variants, versus an expected number in a set of 153,054 chip-genotyped Icelanders, to identify potentially pathogenic genotypes. We follow three missense variants with a complete deficit of homozygosity and find that their pathogenic effect in homozygous state ranges from severe childhood disease to early embryonic lethality. One of these variants is in CPSF3, a gene not previously linked to disease. From a set of clinically sequenced Icelanders, and by sequencing archival samples targeted through the Icelandic genealogy, we find four homozygous carriers. Additionally, we find two homozygous carriers of Mexican descent of another missense variant in CPSF3. All six homozygous carriers of missense variants in CPSF3 show severe intellectual disability, seizures, microcephaly, and abnormal muscle tone. Here, we show how the absence of certain homozygous genotypes from a large population set can elucidate causes of previously unexplained recessive diseases and early miscarriage.

Allele frequency of variants reported to cause adenine phosphoribosyltransferase deficiency.

Adenine phosphoribosyltransferase deficiency is a rare, autosomal recessive disorder of purine metabolism that causes nephrolithiasis and progressive chronic kidney disease. The small number of reported cases indicates an extremely low prevalence, although it has been suggested that missed diagnoses may play a role. We assessed the prevalence of APRT deficiency based on the frequency of causally-related APRT sequence variants in a diverse set of large genomic databases. A thorough search was carried out for all APRT variants that have been confirmed as pathogenic under recessive mode of inheritance, and the frequency of the identified variants examined in six population genomic databases: the deCODE genetics database, the UK Biobank, the 100,000 Genomes Project, the Genome Aggregation Database, the Human Genetic Variation Database and the Korean Variant Archive. The estimated frequency of homozygous genotypes was calculated using the Hardy-Weinberg equation. Sixty-two pathogenic APRT variants were identified, including six novel variants. Most common were the missense variants c.407T>C (p.(Met136Thr)) in Japan and c.194A>T (p.(Asp65Val)) in Iceland, as well as the splice-site variant c.400 + 2dup (p.(Ala108Glufs*3)) in the European population. Twenty-nine variants were detected in at least one of the six genomic databases. The highest cumulative minor allele frequency (cMAF) of pathogenic variants outside of Japan and Iceland was observed in the Irish population (0.2%), though no APRT deficiency cases have been reported in Ireland. The large number of cases in Japan and Iceland is consistent with a founder effect in these populations. There is no evidence for widespread underdiagnosis based on the current analysis.