Sigurjon A. Gudjonsson

Rate of de novo mutations and the importance of father/'s age to disease risk

Mutations generate sequence diversity and provide a substrate for selection. The rate of de novo mutations is therefore of major importance to evolution. We conducted a study of genomewide mutation rate by sequencing the entire genomes of 78 Icelandic parent-offspring trios at high coverage. Here we show that in our samples, with an average father’s age of 29.7, the average de novo mutation rate is 1.20×10−8 per nucleotide per generation. Most strikingly, the diversity in mutation rate of single-nucleotide polymorphism (SNP) is dominated by the age of the father at conception of the child. The effect is an increase of about 2 mutations per year. After accounting for random Poisson variation, father’s age is estimated to explain nearly all of the remaining variation in the de novo mutation counts. These observations shed light on the importance of the father’s age on the risk of diseases such as schizophrenia and autism.

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

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.

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

Identification of a large set of rare complete human knockouts

Loss-of-function mutations cause many mendelian diseases. Here we aimed to create a catalog of autosomal genes that are completely knocked out in humans by rare loss-of-function mutations. We sequenced the whole genomes of 2,636 Icelanders and imputed the sequence variants identified in this set into 101,584 additional chip-genotyped and phased Icelanders. We found a total of 6,795 autosomal loss-of-function SNPs and indels in 4,924 genes. Of the genotyped Icelanders, 7.7% are homozygotes or compound heterozygotes for loss-of-function mutations with a minor allele frequency (MAF) below 2% in 1,171 genes (complete knockouts). Genes that are highly expressed in the brain are less often completely knocked out than other genes. Homozygous loss-of-function offspring of two heterozygous parents occurred less frequently than expected (deficit of 136 per 10,000 transmissions for variants with MAF <2%, 95% confidence interval (CI) = 10–261).

A sequence variant on 17q21 is associated with age at onset and severity of asthma

A sequence variant (rs7216389-T) near the ORMDL3 gene on chromosome 17q21 was recently found to be associated with childhood asthma. We sought to evaluate the effect of rs7216389-T on asthma subphenotypes and its correlation with expression levels of neighboring genes. The association of rs7216389-T with asthma was replicated in six European and one Asian study cohort (N=4917 cases N=34u2009589 controls). In addition, we found that the association of rs7216389-T was confined to cases with early onset of asthma, particularly in early childhood (age: 0–5 years OR=1.51, P=6.89·10−9) and adolescence (age: 14–17 years OR=1.71, P=5.47·10−9). A weaker association was observed for onset between 6 and 13 years of age (OR=1.17, P=0.035), but none for adult-onset asthma (OR=1.07, P=0.12). Cases were further stratified by sex, asthma severity and atopy status. An association with greater asthma severity was observed among early-onset asthma cases (P=0.0012), but no association with sex or atopy status was observed among the asthma cases. An association between sequence variants and the expression of genes in the 17q21 region was assessed in white blood cell RNA samples collected from Icelandic individuals (n=743). rs7216389 associated with the expression of GSDMB and ORMDL3 genes. However, other sequence variants showing a weaker association with asthma compared with that of rs7216389 were more strongly associated with the expression of both genes. Thus, the contribution of rs7216389-T to the development of asthma is unlikely to operate only through an impact on the expression of ORMDL3 or GSDMB genes.

Weighting sequence variants based on their annotation increases power of whole-genome association studies

The consensus approach to genome-wide association studies (GWAS) has been to assign equal prior probability of association to all sequence variants tested. However, some sequence variants, such as loss-of-function and missense variants, are more likely than others to affect protein function and are therefore more likely to be causative. Using data from whole-genome sequencing of 2,636 Icelanders and the association results for 96 quantitative and 123 binary phenotypes, we estimated the enrichment of association signals by sequence annotation. We propose a weighted Bonferroni adjustment that controls for the family-wise error rate (FWER), using as weights the enrichment of sequence annotations among association signals. We show that this weighted adjustment increases the power to detect association over the standard Bonferroni correction. We use the enrichment of associations by sequence annotation we have estimated in Iceland to derive significance thresholds for other populations with different numbers and combinations of sequence variants.

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.

HLA class II sequence variants influence tuberculosis risk in populations of European ancestry

Mycobacterium tuberculosis infections cause 9 million new tuberculosis cases and 1.5 million deaths annually. To identify variants conferring risk of tuberculosis, we tested 28.3 million variants identified through whole-genome sequencing of 2,636 Icelanders for association with tuberculosis (8,162 cases and 277,643 controls), pulmonary tuberculosis (PTB) and M. tuberculosis infection. We found association of three variants in the region harboring genes encoding the class II human leukocyte antigens (HLAs): rs557011[T] (minor allele frequency (MAF) = 40.2%), associated with M. tuberculosis infection (odds ratio (OR) = 1.14, P = 3.1 × 10−13) and PTB (OR = 1.25, P = 5.8 × 10−12), and rs9271378[G] (MAF = 32.5%), associated with PTB (OR = 0.78, P = 2.5 × 10−12)—both located between HLA-DQA1 and HLA-DRB1—and a missense variant encoding p.Ala210Thr in HLA-DQA1 (MAF = 19.1%, rs9272785), associated with M. tuberculosis infection (P = 9.3 × 10−9, OR = 1.14). We replicated association of these variants with PTB in samples of European ancestry from Russia and Croatia (P < 5.9 × 10−4). These findings show that the HLA class II region contributes to genetic risk of tuberculosis, possibly through reduced presentation of protective M. tuberculosis antigens to T cells.

Sequence variants in the PTCH1 gene associate with spine bone mineral density and osteoporotic fractures

Bone mineral density (BMD) is a measure of osteoporosis and is useful in evaluating the risk of fracture. In a genome-wide association study of BMD among 20,100 Icelanders, with follow-up in 10,091 subjects of European and East-Asian descent, we found a new BMD locus that harbours the PTCH1 gene, represented by rs28377268 (freq. 11.4–22.6%) that associates with reduced spine BMD (P=1.0 × 10−11, β=−0.09). We also identified a new spine BMD signal in RSPO3, rs577721086 (freq. 6.8%), that associates with increased spine BMD (P=6.6 × 10−10, β=0.14). Importantly, both variants associate with osteoporotic fractures and affect expression of the PTCH1 and RSPO3 genes that is in line with their influence on BMD and known biological function of these genes. Additional new BMD signals were also found at the AXIN1 and SOST loci and a new lead SNP at the EN1 locus.