Hafdis T. Helgadottir

Sequence variants affecting eosinophil numbers associate with asthma and myocardial infarction

Eosinophils are pleiotropic multifunctional leukocytes involved in initiation and propagation of inflammatory responses and thus have important roles in the pathogenesis of inflammatory diseases. Here we describe a genome-wide association scan for sequence variants affecting eosinophil counts in blood of 9,392 Icelanders. The most significant SNPs were studied further in 12,118 Europeans and 5,212 East Asians. SNPs at 2q12 (rs1420101), 2q13 (rs12619285), 3q21 (rs4857855), 5q31 (rs4143832) and 12q24 (rs3184504) reached genome-wide significance (P = 5.3 × 10−14, 5.4 × 10−10, 8.6 × 10−17, 1.2 × 10−10 and 6.5 × 10−19, respectively). A SNP at IL1RL1 associated with asthma (P = 5.5 × 10−12) in a collection of ten different populations (7,996 cases and 44,890 controls). SNPs at WDR36, IL33 and MYB that showed suggestive association with eosinophil counts were also associated with atopic asthma (P = 4.2 × 10−6, 2.2 × 10−5 and 2.4 × 10−4, respectively). We also found that a nonsynonymous SNP at 12q24, in SH2B3, associated significantly (P = 8.6 × 10−8) with myocardial infarction in six different populations (6,650 cases and 40,621 controls).

Sequence variants at the TERT-CLPTM1L locus associate with many cancer types

The common sequence variants that have recently been associated with cancer risk are particular to a single cancer type or at most two. Following up on our genome-wide scan of basal cell carcinoma, we found that rs401681[C] on chromosome 5p15.33 satisfied our threshold for genome-wide significance (OR = 1.25, P = 3.7 × 10−12). We tested rs401681 for association with 16 additional cancer types in over 30,000 cancer cases and 45,000 controls and found association with lung cancer (OR = 1.15, P = 7.2 × 10−8) and urinary bladder, prostate and cervix cancer (ORs = 1.07−1.31, all P < 4 × 10−4). However, rs401681[C] seems to confer protection against cutaneous melanoma (OR = 0.88, P = 8.0 × 10−4). Notably, most of these cancer types have a strong environmental component to their risk. Investigation of the region led us to rs2736098[A], which showed stronger association with some cancer types. However, neither variant could fully account for the association of the other. rs2736098 corresponds to A305A in the telomerase reverse transcriptase (TERT) protein and rs401681 is in an intron of the CLPTM1L gene.

Common variants on 9q22.33 and 14q13.3 predispose to thyroid cancer in European populations

In order to search for sequence variants conferring risk of thyroid cancer we conducted a genome-wide association study in 192 and 37,196 Icelandic cases and controls, respectively, followed by a replication study in individuals of European descent. Here we show that two common variants, located on 9q22.33 and 14q13.3, are associated with the disease. Overall, the strongest association signals were observed for rs965513 on 9q22.33 (OR = 1.75; P = 1.7 × 10−27) and rs944289 on 14q13.3 (OR = 1.37; P = 2.0 × 10−9). The gene nearest to the 9q22.33 locus is FOXE1 (TTF2) and NKX2-1 (TTF1) is among the genes located at the 14q13.3 locus. Both variants contribute to an increased risk of both papillary and follicular thyroid cancer. Approximately 3.7% of individuals are homozygous for both variants, and their estimated risk of thyroid cancer is 5.7-fold greater than that of noncarriers. In a study on a large sample set from the general population, both risk alleles are associated with low concentrations of thyroid stimulating hormone (TSH), and the 9q22.33 allele is associated with low concentration of thyroxin (T4) and high concentration of triiodothyronine (T3).

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.

Mutations in BRIP1 confer high risk of ovarian cancer

Ovarian cancer causes more deaths than any other gynecologic malignancy in developed countries. Sixteen million sequence variants, identified through whole-genome sequencing of 457 Icelanders, were imputed to 41,675 Icelanders genotyped using SNP chips, as well as to their relatives. Sequence variants were tested for association with ovarian cancer (N of affected individuals = 656). We discovered a rare (0.41% allelic frequency) frameshift mutation, c.2040_2041insTT, in the BRIP1 (FANCJ) gene that confers an increase in ovarian cancer risk (odds ratio (OR) = 8.13, P = 2.8 × 10−14). The mutation was also associated with increased risk of cancer in general and reduced lifespan by 3.6 years. In a Spanish population, another frameshift mutation in BRIP1, c.1702_1703del, was seen in 2 out of 144 subjects with ovarian cancer and 1 out of 1,780 control subjects (P = 0.016). This allele was also associated with breast cancer (seen in 6/927 cases; P = 0.0079). Ovarian tumors from heterozygous carriers of the Icelandic mutation show loss of the wild-type allele, indicating that BRIP1 behaves like a classical tumor suppressor gene in ovarian cancer.

A rare variant in MYH6 is associated with high risk of sick sinus syndrome

Through complementary application of SNP genotyping, whole-genome sequencing and imputation in 38,384 Icelanders, we have discovered a previously unidentified sick sinus syndrome susceptibility gene, MYH6, encoding the alpha heavy chain subunit of cardiac myosin. A missense variant in this gene, c.2161C>T, results in the conceptual amino acid substitution p.Arg721Trp, has an allelic frequency of 0.38% in Icelanders and associates with sick sinus syndrome with an odds ratio = 12.53 and P = 1.5 × 10−29. We show that the lifetime risk of being diagnosed with sick sinus syndrome is around 6% for non-carriers of c.2161C>T but is approximately 50% for carriers of the c.2161C>T variant.

Identification of low-frequency and rare sequence variants associated with elevated or reduced risk of type 2 diabetes.

Through whole-genome sequencing of 2,630 Icelanders and imputation into 11,114 Icelandic cases and 267,140 controls followed by testing in Danish and Iranian samples, we discovered 4 previously unreported variants affecting risk of type 2 diabetes (T2D). A low-frequency (1.47%) variant in intron 1 of CCND2, rs76895963[G], reduces risk of T2D by half (odds ratio (OR) = 0.53, P = 5.0 × 10−21) and is correlated with increased CCND2 expression. Notably, this variant is also associated with both greater height and higher body mass index (1.17 cm per allele, P = 5.5 × 10−12 and 0.56 kg/m2 per allele, P = 6.5 × 10−7, respectively). In addition, two missense variants in PAM, encoding p.Asp563Gly (frequency of 4.98%) and p.Ser539Trp (frequency of 0.65%), confer moderately higher risk of T2D (OR = 1.23, P = 3.9 × 10−10 and OR = 1.47, P = 1.7 × 10−5, respectively), and a rare (0.20%) frameshift variant in PDX1, encoding p.Gly218Alafs*12, associates with high risk of T2D (OR = 2.27, P = 7.3 × 10−7).

Nonsense mutation in the LGR4 gene is associated with several human diseases and other traits

Low bone mineral density (BMD) is used as a parameter of osteoporosis. Genome-wide association studies of BMD have hitherto focused on BMD as a quantitative trait, yielding common variants of small effects that contribute to the population diversity in BMD. Here we use BMD as a dichotomous trait, searching for variants that may have a direct effect on the risk of pathologically low BMD rather than on the regulation of BMD in the healthy population. Through whole-genome sequencing of Icelandic individuals, we found a rare nonsense mutation within the leucine-rich-repeat-containing G-protein-coupled receptor 4 (LGR4) gene (c.376C>T) that is strongly associated with low BMD, and with osteoporotic fractures. This mutation leads to termination of LGR4 at position 126 and fully disrupts its function. The c.376C>T mutation is also associated with electrolyte imbalance, late onset of menarche and reduced testosterone levels, as well as an increased risk of squamous cell carcinoma of the skin and biliary tract cancer. Interestingly, the phenotype of carriers of the c.376C>T mutation overlaps that of Lgr4 mutant mice.

Discovery of common variants associated with low TSH levels and thyroid cancer risk

To search for sequence variants conferring risk of nonmedullary thyroid cancer, we focused our analysis on 22 SNPs with a P < 5 × 10−8 in a genome-wide association study on levels of thyroid stimulating hormone (TSH) in 27,758 Icelanders. Of those, rs965513 has previously been shown to associate with thyroid cancer. The remaining 21 SNPs were genotyped in 561 Icelandic individuals with thyroid cancer (cases) and up to 40,013 controls. Variants suggestively associated with thyroid cancer (P < 0.05) were genotyped in an additional 595 non-Icelandic cases and 2,604 controls. After combining the results, three variants were shown to associate with thyroid cancer: rs966423 on 2q35 (OR = 1.34; Pcombined = 1.3 × 10−9), rs2439302 on 8p12 (OR = 1.36; Pcombined = 2.0 × 10−9) and rs116909374 on 14q13.3 (OR = 2.09; Pcombined = 4.6 × 10−11), a region previously reported to contain an uncorrelated variant conferring risk of thyroid cancer. A strong association (P = 9.1 × 10−91) was observed between rs2439302 on 8p12 and expression of NRG1, which encodes the signaling protein neuregulin 1, in blood.

A study based on whole-genome sequencing yields a rare variant at 8q24 associated with prostate cancer

In Western countries, prostate cancer is the most prevalent cancer of men and one of the leading causes of cancer-related death in men. Several genome-wide association studies have yielded numerous common variants conferring risk of prostate cancer. Here, we analyzed 32.5 million variants discovered by whole-genome sequencing 1,795 Icelanders. We identified a new low-frequency variant at 8q24 associated with prostate cancer in European populations, rs188140481[A] (odds ratio (OR) = 2.90; Pcombined = 6.2 × 10−34), with an average risk allele frequency in controls of 0.54%. This variant is only very weakly correlated (r2 ≤ 0.06) with previously reported risk variants at 8q24, and its association remains significant after adjustment for all known risk-associated variants. Carriers of rs188140481[A] were diagnosed with prostate cancer 1.26 years younger than non-carriers (P = 0.0059). We also report results for a previously described HOXB13 variant (rs138213197[T]), confirming it as a prostate cancer risk variant in populations from across Europe.