Olof Sigurdardottir

Large-scale whole-genome sequencing of the Icelandic population

Here we describe the insights gained from sequencing the whole genomes of 2,636 Icelanders to a median depth of 20×. We found 20 million SNPs and 1.5 million insertions-deletions (indels). We describe the density and frequency spectra of sequence variants in relation to their functional annotation, gene position, pathway and conservation score. We demonstrate an excess of homozygosity and rare protein-coding variants in Iceland. We imputed these variants into 104,220 individuals down to a minor allele frequency of 0.1% and found a recessive frameshift mutation in MYL4 that causes early-onset atrial fibrillation, several mutations in ABCB4 that increase risk of liver diseases and an intronic variant in GNAS associating with increased thyroid-stimulating hormone levels when maternally inherited. These data provide a study design that can be used to determine how variation in the sequence of the human genome gives rise to human diversity.

Variants with large effects on blood lipids and the role of cholesterol and triglycerides in coronary disease

Sequence variants affecting blood lipids and coronary artery disease (CAD) may enhance understanding of the atherogenicity of lipid fractions. Using a large resource of whole-genome sequence data, we examined rare and low-frequency variants for association with non-HDL cholesterol, HDL cholesterol, LDL cholesterol, and triglycerides in up to 119,146 Icelanders. We discovered 13 variants with large effects (within ANGPTL3, APOB, ABCA1, NR1H3, APOA1, LIPC, CETP, LDLR, and APOC1) and replicated 14 variants. Five variants within PCSK9, APOA1, ANGPTL4, and LDLR associate with CAD (33,090 cases and 236,254 controls). We used genetic risk scores for the lipid fractions to examine their causal relationship with CAD. The non-HDL cholesterol genetic risk score associates most strongly with CAD (P = 2.7 × 10−28), and no other genetic risk score associates with CAD after accounting for non-HDL cholesterol. The genetic risk score for non-HDL cholesterol confers CAD risk beyond that of LDL cholesterol (P = 5.5 × 10−8), suggesting that targeting atherogenic remnant cholesterol may reduce cardiovascular risk.

Variant ASGR1 Associated with a Reduced Risk of Coronary Artery Disease

BACKGROUND Several sequence variants are known to have effects on serum levels of non-high-density lipoprotein (HDL) cholesterol that alter the risk of coronary artery disease. METHODS We sequenced the genomes of 2636 Icelanders and found variants that we then imputed into the genomes of approximately 398,000 Icelanders. We tested for association between these imputed variants and non-HDL cholesterol levels in 119,146 samples. We then performed replication testing in two populations of European descent. We assessed the effects of an implicated loss-of-function variant on the risk of coronary artery disease in 42,524 case patients and 249,414 controls from five European ancestry populations. An augmented set of genomes was screened for additional loss-of-function variants in a target gene. We evaluated the effect of an implicated variant on protein stability. RESULTS We found a rare noncoding 12-base-pair (bp) deletion (del12) in intron 4 of ASGR1, which encodes a subunit of the asialoglycoprotein receptor, a lectin that plays a role in the homeostasis of circulating glycoproteins. The del12 mutation activates a cryptic splice site, leading to a frameshift mutation and a premature stop codon that renders a truncated protein prone to degradation. Heterozygous carriers of the mutation (1 in 120 persons in our study population) had a lower level of non-HDL cholesterol than noncarriers, a difference of 15.3 mg per deciliter (0.40 mmol per liter) (P=1.0×10(-16)), and a lower risk of coronary artery disease (by 34%; 95% confidence interval, 21 to 45; P=4.0×10(-6)). In a larger set of sequenced samples from Icelanders, we found another loss-of-function ASGR1 variant (p.W158X, carried by 1 in 1850 persons) that was also associated with lower levels of non-HDL cholesterol (P=1.8×10(-3)). CONCLUSIONS ASGR1 haploinsufficiency was associated with reduced levels of non-HDL cholesterol and a reduced risk of coronary artery disease. (Funded by the National Institutes of Health and others.).

Skin-infiltrating T cells and cytokine expression in Icelandic horses affected with insect bite hypersensitivity: A possible role for regulatory T cells

Equine insect bite hypersensitivity (IBH) is a seasonally recurrent, pruritic skin disorder caused by an IgE-mediated reaction to salivary proteins of biting flies, predominantly of the genus Culicoides. The aim of this study was to define T cell subsets and cytokine profile in the skin of IBH-affected Icelandic horses with particular focus on the balance between T helper (Th) 1, Th2 and T regulatory (Treg) cells. Distribution and number of CD4+, CD8+ and Forkhead box P3 (FoxP3)+ T cells were characterized by immunohistochemical staining in lesional and non-lesional skin of moderately and severely IBH-affected horses (n=14) and in the skin of healthy control horses (n=10). Using real-time quantitative reverse transcription-polymerase chain reaction, mRNA expression levels of Th2 cytokines (Interleukin (IL)-4, IL-5, IL-13), Th1 cytokines (Interferon-γ), regulatory cytokines (Transforming Growth Factor β1, IL-10) and the Treg transcription factor FoxP3 were measured in skin and blood samples. Furthermore, Culicoides nubeculosus specific serum IgE levels were assessed. Lesions of IBH-affected horses contained significantly higher numbers of CD4+ cells than skin of healthy control horses. Furthermore, the total number of T cells (CD4+ and CD8+) was significantly increased in lesional compared to non-lesional skin and there was a tendency (p=0.07) for higher numbers of CD4+ cells in lesional compared to non-lesional skin. While the number of FoxP3+ T cells did not differ significantly between the groups, the ratio of Foxp3 to CD4+ cells was significantly lower in lesions of severely IBH-affected horses than in moderately affected or control horses. Interestingly, differences in FoxP3 expression were more striking at the mRNA level. FoxP3 mRNA levels were significantly reduced in lesional skin, compared both to non-lesional and to healthy skin and were also significantly lower in non-lesional compared to healthy skin. Expression levels of IL-13, but not IL-4 or IL-5, were significantly elevated in lesional and non-lesional skin of IBH-affected horses. IL-10 levels were lower in lesional compared to non-lesional skin (p=0.06) and also lower (p=0.06) in the blood of IBH-affected than of healthy horses. No significant changes were observed regarding blood expression levels of Th1 and Th2 cytokines or FoxP3. Finally, IBH-affected horses had significantly higher Culicoides nubeculosus specific serum IgE levels than control horses. The presented data suggest that an imbalance between Th2 and Treg cells is a characteristic feature in IBH. Treatment strategies for IBH should thus aim at restoring the balance between Th2 and Treg cells.

Common and rare variants associated with kidney stones and biochemical traits

Kidney stone disease is a complex disorder with a strong genetic component. We conducted a genome-wide association study of 28.3 million sequence variants detected through whole-genome sequencing of 2,636 Icelanders that were imputed into 5,419 kidney stone cases, including 2,172 cases with a history of recurrent kidney stones, and 279,870 controls. We identify sequence variants associating with kidney stones at ALPL (rs1256328[T], odds ratio (OR)=1.21, P=5.8 × 10−10) and a suggestive association at CASR (rs7627468[A], OR=1.16, P=2.0 × 10−8). Focusing our analysis on coding sequence variants in 63 genes with preferential kidney expression we identify two rare missense variants SLC34A1 p.Tyr489Cys (OR=2.38, P=2.8 × 10−5) and TRPV5 p.Leu530Arg (OR=3.62, P=4.1 × 10−5) associating with recurrent kidney stones. We also observe associations of the identified kidney stone variants with biochemical traits in a large population set, indicating potential biological mechanism.

A genome-wide association study yields five novel thyroid cancer risk loci

The great majority of thyroid cancers are of the non-medullary type. Here we report findings from a genome-wide association study of non-medullary thyroid cancer, including in total 3,001 patients and 287,550 controls from five study groups of European descent. Our results yield five novel loci (all with Pcombined<3 × 10−8): 1q42.2 (rs12129938 in PCNXL2), 3q26.2 (rs6793295 a missense mutation in LRCC34 near TERC), 5q22.1 (rs73227498 between NREP and EPB41L4A), 10q24.33 (rs7902587 near OBFC1), and two independently associated variants at 15q22.33 (rs2289261 and rs56062135; both in SMAD3). We also confirm recently published association results from a Chinese study of a variant on 5p15.33 (rs2736100 near the TERT gene) and present a stronger association result for a moderately correlated variant (rs10069690; OR=1.20, P=3.2 × 10−7) based on our study of individuals of European ancestry. In combination, these results raise several opportunities for future studies of the pathogenesis of thyroid cancer.

A rare IL33 loss-of-function mutation reduces blood eosinophil counts and protects from asthma

IL-33 is a tissue-derived cytokine that induces and amplifies eosinophilic inflammation and has emerged as a promising new drug target for asthma and allergic disease. Common variants at IL33 and IL1RL1, encoding the IL-33 receptor ST2, associate with eosinophil counts and asthma. Through whole-genome sequencing and imputation into the Icelandic population, we found a rare variant in IL33 (NM_001199640:exon7:c.487-1G>C (rs146597587-C), allele frequency = 0.65%) that disrupts a canonical splice acceptor site before the last coding exon. It is also found at low frequency in European populations. rs146597587-C associates with lower eosinophil counts (β = -0.21 SD, P = 2.5×10–16, N = 103,104), and reduced risk of asthma in Europeans (OR = 0.47; 95%CI: 0.32, 0.70, P = 1.8×10–4, N cases = 6,465, N controls = 302,977). Heterozygotes have about 40% lower total IL33 mRNA expression than non-carriers and allele-specific analysis based on RNA sequencing and phased genotypes shows that only 20% of the total expression is from the mutated chromosome. In half of those transcripts the mutation causes retention of the last intron, predicted to result in a premature stop codon that leads to truncation of 66 amino acids. The truncated IL-33 has normal intracellular localization but neither binds IL-33R/ST2 nor activates ST2-expressing cells. Together these data demonstrate that rs146597587-C is a loss of function mutation and support the hypothesis that IL-33 haploinsufficiency protects against asthma.

Fatal infection in two Icelandic stallions caused by Halicephalobus gingivalis (Nematoda: Rhabditida)

Opportunistic infections with the free living nematode Halicephalobus gingivalis are infrequently reported in horses but the cases are widespread geographically. The nematodes are believed to penetrate wounds and subsequently reproduce within the host tissues. This paper reports two cases of a fatal disease in stallions of the Icelandic breed in Iceland. Case 1: a stallion, which sustained injuries to the mouth after an accident, developed severe neurological signs and had to be euthanatized. Histological examination revealed mild inflammation and malacia in the cerebellum associated with the presence of numerous H. gingivalis nematodes. Case 2: a stallion that started swerving to one side and lost balance was euthanatized due to lack of response to therapy and rapid deterioration. Histological examination revealed numerous H. gingivalis nematodes in the cerebellum, brain stem, cervical spinal cord and in the meninges, with minimal reactive changes. In case 1 the infection presumably was acquired by nematodes from soil penetrating through wounds in the mouth. The mode of the H. gingivalis infection in case 2 is uncertain. These are the first cases of H. gingivalis infection reported from Iceland and the second report from the Nordic countries.

Common and rare variants associating with serum levels of creatine kinase and lactate dehydrogenase

Creatine kinase (CK) and lactate dehydrogenase (LDH) are widely used markers of tissue damage. To search for sequence variants influencing serum levels of CK and LDH, 28.3 million sequence variants identified through whole-genome sequencing of 2,636 Icelanders were imputed into 63,159 and 98,585 people with CK and LDH measurements, respectively. Here we describe 13 variants associating with serum CK and 16 with LDH levels, including four that associate with both. Among those, 15 are non-synonymous variants and 12 have a minor allele frequency below 5%. We report sequence variants in genes encoding the enzymes being measured (CKM and LDHA), as well as in genes linked to muscular (ANO5) and immune/inflammatory function (CD163/CD163L1, CSF1, CFH, HLA-DQB1, LILRB5, NINJ1 and STAB1). A number of the genes are linked to the mononuclear/phagocyte system and clearance of enzymes from the serum. This highlights the variety in the sources of normal diversity in serum levels of enzymes.

A Splice Region Variant in LDLR Lowers Non-high Density Lipoprotein Cholesterol and Protects against Coronary Artery Disease

Through high coverage whole-genome sequencing and imputation of the identified variants into a large fraction of the Icelandic population, we found four independent signals in the low density lipoprotein receptor gene (LDLR) that associate with levels of non-high density lipoprotein cholesterol (non-HDL-C) and coronary artery disease (CAD). Two signals are novel with respect to association with non-HDL-C and are represented by non-coding low frequency variants (between 2–4% frequency), the splice region variant rs72658867-A in intron 14 and rs17248748-T in intron one. These two novel associations were replicated in three additional populations. Both variants lower non-HDL-C levels (rs72658867-A, non-HDL-C effect = -0.44 mmol/l, P adj = 1.1 × 10−80 and rs17248748-T, non-HDL-C effect = -0.13 mmol/l, P adj = 1.3 × 10−12) and confer protection against CAD (rs72658867-A, OR = 0.76 and P adj = 2.7 × 10−8 and rs17248748-T, OR = 0.92 and P adj = 0.022). The LDLR splice region variant, rs72658867-A, located at position +5 in intron 14 (NM_000527:c.2140+5G>A), causes retention of intron 14 during transcription and is expected to produce a truncated LDL receptor lacking domains essential for function of the receptor. About half of the transcripts generated from chromosomes carrying rs72658867-A are characterized by this retention of the intron. The same variant also increases LDLR mRNA expression, however, the wild type transcripts do not exceed levels in non-carriers. This demonstrates that sequence variants that disrupt the LDL receptor can lower non-HDL-C and protect against CAD.