Thorvaldur Jonsson

Genome-wide association yields new sequence variants at seven loci that associate with measures of obesity

Obesity results from the interaction of genetic and environmental factors. To search for sequence variants that affect variation in two common measures of obesity, weight and body mass index (BMI), both of which are highly heritable, we performed a genome-wide association (GWA) study with 305,846 SNPs typed in 25,344 Icelandic, 2,998 Dutch, 1,890 European Americans and 1,160 African American subjects and combined the results with previously published results from the Diabetes Genetics Initiative (DGI) on 3,024 Scandinavians. We selected 43 variants in 19 regions for follow-up in 5,586 Danish individuals and compared the results to a genome-wide study on obesity-related traits from the GIANT consortium. In total, 29 variants, some correlated, in 11 chromosomal regions reached a genome-wide significance threshold of P < 1.6 × 10−7. This includes previously identified variants close to or in the FTO, MC4R, BDNF and SH2B1 genes, in addition to variants at seven loci not previously connected with obesity.

Common variants on chromosomes 2q35 and 16q12 confer susceptibility to estrogen receptor-positive breast cancer.

Familial clustering studies indicate that breast cancer risk has a substantial genetic component. To identify new breast cancer risk variants, we genotyped approximately 300,000 SNPs in 1,600 Icelandic individuals with breast cancer and 11,563 controls using the Illumina Hap300 platform. We then tested selected SNPs in five replication sample sets. Overall, we studied 4,554 affected individuals and 17,577 controls. Two SNPs consistently associated with breast cancer: ∼25% of individuals of European descent are homozygous for allele A of rs13387042 on chromosome 2q35 and have an estimated 1.44-fold greater risk than noncarriers, and for allele T of rs3803662 on 16q12, about 7% are homozygous and have a 1.64-fold greater risk. Risk from both alleles was confined to estrogen receptor–positive tumors. At present, no genes have been identified in the linkage disequilibrium block containing rs13387042. rs3803662 is near the 5′ end of TNRC9 , a high mobility group chromatin–associated protein whose expression is implicated in breast cancer metastasis to bone.

Many sequence variants affecting diversity of adult human height

Adult human height is one of the classical complex human traits. We searched for sequence variants that affect height by scanning the genomes of 25,174 Icelanders, 2,876 Dutch, 1,770 European Americans and 1,148 African Americans. We then combined these results with previously published results from the Diabetes Genetics Initiative on 3,024 Scandinavians and tested a selected subset of SNPs in 5,517 Danes. We identified 27 regions of the genome with one or more sequence variants showing significant association with height. The estimated effects per allele of these variants ranged between 0.3 and 0.6 cm and, taken together, they explain around 3.7% of the population variation in height. The genes neighboring the identified loci cluster in biological processes related to skeletal development and mitosis. Association to three previously reported loci are replicated in our analyses, and the strongest association was with SNPs in the ZBTB38 gene.

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.

Parental origin of sequence variants associated with complex diseases.

Effects of susceptibility variants may depend on from which parent they are inherited. Although many associations between sequence variants and human traits have been discovered through genome-wide associations, the impact of parental origin has largely been ignored. Here we show that for 38,167 Icelanders genotyped using single nucleotide polymorphism (SNP) chips, the parental origin of most alleles can be determined. For this we used a combination of genealogy and long-range phasing. We then focused on SNPs that associate with diseases and are within 500 kilobases of known imprinted genes. Seven independent SNP associations were examined. Five—one with breast cancer, one with basal-cell carcinoma and three with type 2 diabetes—have parental-origin-specific associations. These variants are located in two genomic regions, 11p15 and 7q32, each harbouring a cluster of imprinted genes. Furthermore, we observed a novel association between the SNP rs2334499 at 11p15 and type 2 diabetes. Here the allele that confers risk when paternally inherited is protective when maternally transmitted. We identified a differentially methylated CTCF-binding site at 11p15 and demonstrated correlation of rs2334499 with decreased methylation of that site.

Common variants on chromosome 5p12 confer susceptibility to estrogen receptor-positive breast cancer.

We carried out a genome-wide association study of breast cancer predisposition with replication and refinement studies involving 6,145 cases and 33,016 controls and identified two SNPs (rs4415084 and rs10941679) on 5p12 that confer risk, preferentially for estrogen receptor (ER)-positive tumors (OR = 1.27, P = 2.5 × 10−12 for rs10941679). The nearest gene, MRPS30, was previously implicated in apoptosis, ER-positive tumors and favorable prognosis. A recently reported signal in FGFR2 was also found to associate specifically with ER-positive breast cancer.

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

Genome-wide association study identifies sequence variants on 6q21 associated with age at menarche

Earlier menarche correlates with shorter adult height and higher childhood body fat. We conducted a genome-wide association study of age at menarche (AAM) on 15,297 Icelandic women. Combined analysis with replication sets from Iceland, Denmark and the Netherlands (N = 10,040) yielded a significant association between rs314280[T] on 6q21, near the LIN28B gene, and AAM (effect = 1.2 months later per allele; P = 1.8 × 10−14). A second SNP within the same linkage disequilibrium (LD) block, rs314277, splits rs314280[T] into two haplotypes with different effects (0.9 months and 1.9 months per allele). These variants have been associated with greater adult height. The association with adult height did not account for the association with AAM or vice versa. Other variants, previously associated with height, did not associate significantly with AAM. Given the link between body fat and AAM, we also assessed 11 variants recently associated with higher body mass index (BMI) and 5 of those associated with earlier AAM.

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.