The correlation between CpG methylation and gene expression is driven by sequence variants.

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

Predicted loss and gain of function mutations in ACO1 are associated with erythropoiesis.

Hemoglobin is the essential oxygen-carrying molecule in humans and is regulated by cellular iron and oxygen sensing mechanisms. To search for novel variants associated with hemoglobin concentration, we performed genome-wide association studies of hemoglobin concentration using a combined set of 684,122 individuals from Iceland and the UK. Notably, we found seven novel variants, six rare coding and one common, at the ACO1 locus associating with either decreased or increased hemoglobin concentration. Of these variants, the missense Cys506Ser and the stop-gained Lys334Ter mutations are specific to eight and ten generation pedigrees, respectively, and have the two largest effects in the study (EffectCys506Ser = -1.61 SD, CI95 = [-1.98, -1.35]; EffectLys334Ter = 0.63 SD, CI95 = [0.36, 0.91]). We also find Cys506Ser to associate with increased risk of persistent anemia (OR = 17.1, P = 2 × 10-14). The strong bidirectional effects seen in this study implicate ACO1, a known iron sensing molecule, as a major homeostatic regulator of hemoglobin concentration.

Meta-analysis of Icelandic and UK data sets identifies missense variants in SMO, IL11, COL11A1 and 13 more new loci associated with osteoarthritis.

Osteoarthritis has a highly negative impact on quality of life because of the associated pain and loss of joint function. Here we describe the largest meta-analysis so far of osteoarthritis of the hip and the knee in samples from Iceland and the UK Biobank (including 17,151 hip osteoarthritis patients, 23,877 knee osteoarthritis patients, and more than 562,000 controls). We found 23 independent associations at 22 loci in the additive meta-analyses, of which 16 of the loci were novel: 12 for hip and 4 for knee osteoarthritis. Two associations are between rare or low-frequency missense variants and hip osteoarthritis, affecting the genes SMO (rs143083812, frequency 0.11%, odds ratio (OR) = 2.8, P = 7.9 × 10-12, p.Arg173Cys) and IL11 (rs4252548, frequency 2.08%, OR = 1.30, P = 2.1 × 10-11, p.Arg112His). A common missense variant in the COL11A1 gene also associates with hip osteoarthritis (rs3753841, frequency 61%, P = 5.2 × 10-10, OR = 1.08, p.Pro1284Leu). In addition, using a recessive model, we confirm an association between hip osteoarthritis and a variant of CHADL1 (rs117018441, P = 1.8 × 10-25, OR = 5.9). Furthermore, we observe a complex relationship between height and risk of osteoarthritis.