Hrefna Johannsdottir

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.

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.

Loss-of-function variants in ABCA7 confer risk of Alzheimer's disease.

We conducted a search for rare, functional variants altering susceptibility to Alzheimers disease that exploited knowledge of common variants associated with the same disease. We found that loss-of-function variants in ABCA7 confer risk of Alzheimers disease in Icelanders (odds ratio (OR) = 2.12, P = 2.2 × 10−13) and discovered that the association replicated in study groups from Europe and the United States (combined OR = 2.03, P = 6.8 × 10−15).

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.

Chromosome 5 imbalance mapping in breast tumors from BRCA1 and BRCA2 mutation carriers and sporadic breast tumors

Comparative genomic hybridization (CGH) analysis has shown that chromosome 5q deletions are the most frequent aberration in breast tumors from BRCA1 mutation carriers. To map the location of putative 5q tumor suppressor gene(s), 26 microsatellite markers covering chromosome 5 were used in loss of heterozygosity (LOH) analysis of breast tumors from BRCA1 (n = 42) and BRCA2 mutation carriers (n = 67), as well as in sporadic cases (n = 65). High‐density array CGH was also used to map chromosome 5 imbalance in 10 BRCA1 tumors. A high LOH frequency was found in BRCA1 tumors (range 19–82%), as compared to BRCA2 and sporadic tumors (ranges 11–44% and 7–43%, respectively). In all, 11 distinct chromosome 5 regions with LOH were observed, the most frequent being 5q35.3 (82%), 5q14.2 (71%) and 5q33.1 (69%) in BRCA1 tumors; 5q35.3 (44%), 5q31.3 (43%) and 5q13.3 (43%) in BRCA2 tumors and 5q31.3 (43%) in sporadic tumors. Array CGH analysis confirmed the very high frequency of 5q deletions, including candidate tumor suppressor genes such as XRCC4, RAD50, RASA1, APC and PPP2R2B. In addition, 2 distinct homozygous deletions were identified, spanning regions of 0.7–1.5 Mbp on 5q12.1 and 5q12.3‐q13.1, respectively. These regions include only a few genes, most notably BRCC3/DEPDC1B (pleckstrin/G protein interacting and RhoGAP domains) and PIK3R1 (PI3 kinase P85 regulatory subunit). Significant association (p ≤ 0.05) was found between LOH at certain 5q regions and factors of poor prognosis, including negative estrogen and progesterone receptor status, high grade, large tumor size and high portion of cells in S‐phase. In conclusion, our results confirm a very high prevalence of chromosome 5q alterations in BRCA1 tumors, pinpointing new regions and genes that should be further investigated.

Loss-of-function variants in ATM confer risk of gastric cancer

Gastric cancer is a serious health problem worldwide, with particularly high prevalence in eastern Asia. Genome-wide association studies (GWAS) in Asian populations have identified several loci that associate with gastric cancer risk. Here we report a GWAS of gastric cancer in a European population, using information on 2,500 population-based gastric cancer cases and 205,652 controls. We found a new gastric cancer association with loss-of-function mutations in ATM (gene test, P = 8.0 × 10−12; odds ratio (OR) = 4.74). The combination of the loss-of-function variants p.Gln852*, p.Ser644* and p.Tyr103* (combined minor allele frequency (MAF) = 0.3%) also associates with pancreatic and prostate cancers (OR = 3.81 and 2.18, respectively) and gives an indication of risk of breast and colorectal cancers (OR = 1.82 and 1.97, respectively). Cancers in those carrying loss-of-function ATM mutations are diagnosed at a significantly earlier age than in non-carriers. Our results confirm an association between gastric cancer in Europeans and three loci previously reported in Asians, MUC1, PRKAA1 and PSCA, refine the association signal at PRKAA1 and support a pathogenic role for the tandem repeat identified in MUC1.

A common variant at 8q24.21 is associated with renal cell cancer

Renal cell carcinoma (RCC) represents between 80 and 90% of kidney cancers. Previous genome-wide association studies of RCC have identified five variants conferring risk of the disease. Here we report the results from a discovery RCC genome-wide association study and replication analysis, including a total of 2,411 patients and 71,497 controls. One variant, rs35252396[CG] located at 8q24.21, is significantly associated with RCC after combining discovery and replication results (OR=1.27, P(combined)=5.4 × 10(-11)) and has an average risk allele frequency in controls of 46%. rs35252396[CG] does not have any strongly correlated variants in the genome and is located within a region predicted to have regulatory functions in several cell lines, including six originating from the kidney. This is the first RCC variant reported at 8q24.21 and it is largely independent (r(2)≤0.02) of the numerous previously reported cancer risk variants at this locus.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus

Background Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. Aim To identify CNVs associated with osteoporotic bone fracture risk. Method We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. Results A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p=8.69×10−5). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p=0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. Conclusions These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.