Genome-wide association study reveals a locus in ADARB2 for complete freedom from headache in Danish Blood Donors.

Headache disorders are the most common disorders of the nervous system. The lifetime prevalence of headache disorders show that some individuals never experience headache. The etiology of complete freedom from headache is not known. To assess genetic variants associated with complete freedom from headache, we performed a genome-wide association study of individuals who have never experienced a headache. We included 63,992 individuals (2,998 individuals with complete freedom from headache and 60,994 controls) from the Danish Blood Donor Study Genomic Cohort. Participants were included in two rounds, from 2015 to 2018 and in 2020. We discovered a genome-wide significant association, with the lead variant rs7904615[G] in ADARB2 (EAF = 27%, OR = 1.20 [1.13-1.27], p = 3.92 × 10-9). The genomic locus was replicated in a non-overlapping cohort of 13,032 individuals (539 individuals with complete freedom from headache and 12,493 controls) from the Danish Blood Donor Study Genomic Cohort (p < 0.05, two-sided). Participants for the replication were included from 2015 to 2020. In conclusion, we show that complete freedom from headache has a genetic component, and we suggest that ADARB2 is involved in complete freedom from headache. The genomic locus was specific for complete freedom from headache and was not associated with any primary headache disorders.

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.

Long-read sequencing of 3,622 Icelanders provides insight into the role of structural variants in human diseases and other traits.

Long-read sequencing (LRS) promises to improve the characterization of structural variants (SVs). We generated LRS data from 3,622 Icelanders and identified a median of 22,636 SVs per individual (a median of 13,353 insertions and 9,474 deletions). We discovered a set of 133,886 reliably genotyped SV alleles and imputed them into 166,281 individuals to explore their effects on diseases and other traits. We discovered an association of a rare deletion in PCSK9 with lower low-density lipoprotein (LDL) cholesterol levels, compared to the population average. We also discovered an association of a multiallelic SV in ACAN with height; we found 11 alleles that differed in the number of a 57-bp-motif repeat and observed a linear relationship between the number of repeats carried and height. These results show that SVs can be accurately characterized at the population scale using LRS data in a genome-wide non-targeted approach and demonstrate how SVs impact phenotypes.

Sequence Variants in TAAR5 and Other Loci Affect Human Odor Perception and Naming.

Olfactory receptor (OR) genes in humans form a special class characterized by unusually high DNA sequence diversity, which should give rise to differences in perception and behavior. In the largest genome-wide association study to date based on olfactory testing, we investigated odor perception and naming with smell tasks performed by 9,122 Icelanders, with replication in a separate sample of 2,204 individuals. We discovered an association between a low-frequency missense variant in TAAR5 and reduced intensity rating of fish odor containing trimethylamine (p.Ser95Pro, pcombined = 5.6 × 10-15). We demonstrate that TAAR5 genotype affects aversion to fish odor, reflected by linguistic descriptions of the odor and pleasantness ratings. We also discovered common sequence variants in two canonical olfactory receptor loci that associate with increased intensity and naming of licorice odor (trans-anethole: lead variant p.Lys233Asn in OR6C70, pcombined = 8.8 × 10-16 and pcombined = 1.4 × 10-9) and enhanced naming of cinnamon (trans-cinnamaldehyde; intergenic variant rs317787-T, pcombined = 5.0 × 10-17). Together, our results show that TAAR5 genotype variation influences human odor responses and highlight that sequence diversity in canonical OR genes can lead to enhanced olfactory ability, in contrast to the view that greater tolerance for mutations in the human OR repertoire leads to diminished function.

Predicted loss and gain of function mutations in ACO1 are associated with erythropoiesis.

Hemoglobin is the essential oxygen-carrying molecule in humans and is regulated by cellular iron and oxygen sensing mechanisms. To search for novel variants associated with hemoglobin concentration, we performed genome-wide association studies of hemoglobin concentration using a combined set of 684,122 individuals from Iceland and the UK. Notably, we found seven novel variants, six rare coding and one common, at the ACO1 locus associating with either decreased or increased hemoglobin concentration. Of these variants, the missense Cys506Ser and the stop-gained Lys334Ter mutations are specific to eight and ten generation pedigrees, respectively, and have the two largest effects in the study (EffectCys506Ser = -1.61 SD, CI95 = [-1.98, -1.35]; EffectLys334Ter = 0.63 SD, CI95 = [0.36, 0.91]). We also find Cys506Ser to associate with increased risk of persistent anemia (OR = 17.1, P = 2 × 10-14). The strong bidirectional effects seen in this study implicate ACO1, a known iron sensing molecule, as a major homeostatic regulator of hemoglobin concentration.

Genome-wide association identifies seven loci for pelvic organ prolapse in Iceland and the UK Biobank.

Pelvic organ prolapse (POP) is a downward descent of one or more of the pelvic organs, resulting in a protrusion of the vaginal wall and/or uterus. We performed a genome-wide association study of POP using data from Iceland and the UK Biobank, a total of 15,010 cases with hospital-based diagnosis code and 340,734 female controls, and found eight sequence variants at seven loci associating with POP (P < 5 × 10-8); seven common (minor allele frequency >5%) and one with minor allele frequency of 4.87%. Some of the variants associating with POP also associated with traits of similar pathophysiology. Of these, rs3820282, which may alter the estrogen-based regulation of WNT4, also associates with leiomyoma of uterus, gestational duration and endometriosis. Rs3791675 at EFEMP1, a gene involved in connective tissue homeostasis, also associates with hernias and carpal tunnel syndrome. Our results highlight the role of connective tissue metabolism and estrogen exposure in the etiology of POP.

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 loss-of-function variant in ALOX15 protects against nasal polyps and chronic rhinosinusitis.

Nasal polyps (NP) are lesions on the nasal and paranasal sinus mucosa and are a risk factor for chronic rhinosinusitis (CRS). We performed genome-wide association studies on NP and CRS in Iceland and the UK (using UK Biobank data) with 4,366 NP cases, 5,608 CRS cases, and >700,000 controls. We found 10 markers associated with NP and 2 with CRS. We also tested 210 markers reported to associate with eosinophil count, yielding 17 additional NP associations. Of the 27 NP signals, 7 associate with CRS and 13 with asthma. Most notably, a missense variant in ALOX15 that causes a p.Thr560Met alteration in arachidonate 15-lipoxygenase (15-LO) confers large genome-wide significant protection against NP (P = 8.0 × 10-27, odds ratio = 0.32; 95% confidence interval = 0.26, 0.39) and CRS (P = 1.1 × 10-8, odds ratio = 0.64; 95% confidence interval = 0.55, 0.75). p.Thr560Met, carried by around 1 in 20 Europeans, was previously shown to cause near total loss of 15-LO enzymatic activity. Our findings identify 15-LO as a potential target for therapeutic intervention in NP and CRS.

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.

Genome-wide association meta-analysis yields 20 loci associated with gallstone disease.

Gallstones are responsible for one of the most common diseases in the Western world and are commonly treated with cholecystectomy. We perform a meta-analysis of two genome-wide association studies of gallstone disease in Iceland and the UK, totaling 27,174 cases and 736,838 controls, uncovering 21 novel gallstone-associated variants at 20 loci. Two distinct low frequency missense variants in SLC10A2, encoding the apical sodium-dependent bile acid transporter (ASBT), associate with an increased risk of gallstone disease (Pro290Ser: OR = 1.36 [1.25-1.49], P = 2.1 × 10-12, MAF = 1%; Val98Ile: OR = 1.15 [1.10-1.20], P = 1.8 × 10-10, MAF = 4%). We demonstrate that lower bile acid transport by ASBT is accompanied by greater risk of gallstone disease and highlight the role of the intestinal compartment of the enterohepatic circulation of bile acids in gallstone disease susceptibility. Additionally, two low frequency missense variants in SERPINA1 and HNF4A and 17 common variants represent novel associations with gallstone disease.

A homozygous loss-of-function mutation leading to CYBC1 deficiency causes chronic granulomatous disease.

Mutations in genes encoding subunits of the phagocyte NADPH oxidase complex are recognized to cause chronic granulomatous disease (CGD), a severe primary immunodeficiency. Here we describe how deficiency of CYBC1, a previously uncharacterized protein in humans (C17orf62), leads to reduced expression of NADPH oxidase's main subunit (gp91phox) and results in CGD. Analyzing two brothers diagnosed with CGD we identify a homozygous loss-of-function mutation, p.Tyr2Ter, in CYBC1. Imputation of p.Tyr2Ter into 155K chip-genotyped Icelanders reveals six additional homozygotes, all with signs of CGD, manifesting as colitis, rare infections, or a severely impaired PMA-induced neutrophil oxidative burst. Homozygosity for p.Tyr2Ter consequently associates with inflammatory bowel disease (IBD) in Iceland (P = 8.3 × 10-8; OR = 67.6), as well as reduced height (P = 3.3 × 10-4; -8.5 cm). Overall, we find that CYBC1 deficiency results in CGD characterized by colitis and a distinct profile of infections indicative of macrophage dysfunction.

A rare missense variant in NR1H4 associates with lower cholesterol levels.

Searching for novel sequence variants associated with cholesterol levels is of particular interest due to the causative role of non-HDL cholesterol levels in cardiovascular disease. Through whole-genome sequencing of 15,220 Icelanders and imputation of the variants identified, we discovered a rare missense variant in NR1H4 (R436H) associating with lower levels of total cholesterol (effect = -0.47 standard deviations or -0.55 mmol L-1, p = 4.21 × 10-10, N = 150,211). Importantly, NR1H4 R436H also associates with lower levels of non-HDL cholesterol and, consistent with this, protects against coronary artery disease. NR1H4 encodes FXR that regulates bile acid homeostasis, however, we do not detect a significant association between R436H and biological markers of liver function. Transcriptional profiling of hepatocytes carrying R436H shows that it is not a loss-of-function variant. Rather, we observe changes in gene expression compatible with effects on lipids. These findings highlight the role of FXR in regulation of cholesterol levels in humans.

A truncating mutation in EPOR leads to hypo-responsiveness to erythropoietin with normal haemoglobin.

The cytokine erythropoietin (EPO), signalling through the EPO receptor (EPO-R), is essential for the formation of red blood cells. We performed a genome-wide association study (GWAS) testing 32.5 million sequence variants for association with serum EPO levels in a set of 4187 individuals. We detect an association between a rare and well imputed stop-gained variant rs370865377[A] (p.Gln82Ter) in EPOR, carried by 1 in 550 Icelanders, and increased serum EPO levels (MAF = 0.09%, Effect = 1.47 SD, P = 3.3 × 10-7). We validated these findings by measuring serum EPO levels in 34 additional pairs of carriers and matched controls and found carriers to have 3.23-fold higher EPO levels than controls (P = 1.7 × 10-6; P combined = 1.6 × 10-11). In contrast to previously reported EPOR mutations, p.Gln82Ter does not associate with haemoglobin levels (Effect = -0.045 SD, P = 0.32, N = 273,160), probably due to a compensatory EPO upregulation in response to EPO-R hypo-responsiveness.

Coding variants in RPL3L and MYZAP increase risk of atrial fibrillation.

Most sequence variants identified hitherto in genome-wide association studies (GWAS) of atrial fibrillation are common, non-coding variants associated with risk through unknown mechanisms. We performed a meta-analysis of GWAS of atrial fibrillation among 29,502 cases and 767,760 controls from Iceland and the UK Biobank with follow-up in samples from Norway and the US, focusing on low-frequency coding and splice variants aiming to identify causal genes. We observe associations with one missense (OR = 1.20) and one splice-donor variant (OR = 1.50) in RPL3L, the first ribosomal gene implicated in atrial fibrillation to our knowledge. Analysis of 167 RNA samples from the right atrium reveals that the splice-donor variant in RPL3L results in exon skipping. We also observe an association with a missense variant in MYZAP (OR = 1.38), encoding a component of the intercalated discs of cardiomyocytes. Both discoveries emphasize the close relationship between the mechanical and electrical function of the heart.

Variants associating with uterine leiomyoma highlight genetic background shared by various cancers and hormone-related traits.

Uterine leiomyomas are common benign tumors of the myometrium. We performed a meta-analysis of two genome-wide association studies of leiomyoma in European women (16,595 cases and 523,330 controls), uncovering 21 variants at 16 loci that associate with the disease. Five variants were previously reported to confer risk of various malignant or benign tumors (rs78378222 in TP53, rs10069690 in TERT, rs1800057 and rs1801516 in ATM, and rs7907606 at OBFC1) and four signals are located at established risk loci for hormone-related traits (endometriosis and breast cancer) at 1q36.12 (CDC42/WNT4), 2p25.1 (GREB1), 20p12.3 (MCM8), and 6q26.2 (SYNE1/ESR1). Polygenic score for leiomyoma, computed using UKB data, is significantly correlated with risk of cancer in the Icelandic population. Functional annotation suggests that the non-coding risk variants affect multiple genes, including ESR1. Our results provide insights into the genetic background of leiomyoma that are shared by other benign and malignant tumors and highlight the role of hormones in leiomyoma growth.

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.

Sequence variant at 4q25 near PITX2 associates with appendicitis.

Appendicitis is one of the most common conditions requiring acute surgery and can pose a threat to the lives of affected individuals. We performed a genome-wide association study of appendicitis in 7,276 Icelandic and 1,139 Dutch cases and large groups of controls. In a combined analysis of the Icelandic and Dutch data, we detected a single signal represented by an intergenic variant rs2129979 [G] close to the gene PITX2 associating with increased risk of appendicitis (OR = 1.15, P = 1.8 × 10-11). We only observe the association in patients diagnosed in adulthood. The marker is close to, but distinct from, a set of markers reported to associate with atrial fibrillation, which have been linked to PITX2. PITX2 has been implicated in determination of right-left symmetry during development. Anomalies in organ arrangement have been linked to increased prevalence of gastrointestinal and intra-abdominal complications, which may explain the effect of rs2129979 on appendicitis risk.

A rare splice donor mutation in the haptoglobin gene associates with blood lipid levels and coronary artery disease.

Common sequence variants at the haptoglobin gene (HP) have been associated with blood lipid levels. Through whole-genome sequencing of 8,453 Icelanders, we discovered a splice donor founder mutation in HP (NM_001126102.1:c.190 + 1G > C, minor allele frequency = 0.56%). This mutation occurs on the HP1 allele of the common copy number variant in HP and leads to a loss of function of HP1. It associates with lower levels of haptoglobin (P = 2.1 × 10-54), higher levels of non-high density lipoprotein cholesterol (ß = 0.26 mmol/l, P = 2.6 × 10-9) and greater risk of coronary artery disease (odds ratio = 1.30, 95% confidence interval: 1.10-1.54, P = 0.0024). Through haplotype analysis and with RNA sequencing, we provide evidence of a causal relationship between one of the two haptoglobin isoforms, namely Hp1, and lower levels of non-HDL cholesterol. Furthermore, we show that the HP1 allele associates with various other quantitative biological traits.

Epigenetic and genetic components of height regulation.

Adult height is a highly heritable trait. Here we identified 31.6 million sequence variants by whole-genome sequencing of 8,453 Icelanders and tested them for association with adult height by imputing them into 88,835 Icelanders. Here we discovered 13 novel height associations by testing four different models including parent-of-origin (|ß|=0.4-10.6 cm). The minor alleles of three parent-of-origin signals associate with less height only when inherited from the father and are located within imprinted regions (IGF2-H19 and DLK1-MEG3). We also examined the association of these sequence variants in a set of 12,645 Icelanders with birth length measurements. Two of the novel variants, (IGF2-H19 and TET1), show significant association with both adult height and birth length, indicating a role in early growth regulation. Among the parent-of-origin signals, we observed opposing parental effects raising questions about underlying mechanisms. These findings demonstrate that common variations affect human growth by parental imprinting.