Protein glycosylation is a ubiquitous modification that affects the structure and function of proteins. Our recent genome wide association study identified transcription factor HNF1A as an important regulator of plasma protein glycosylation. To evaluate the potential impact of epigenetic regulation of HNF1A on protein glycosylation we analyzed CpG methylation in 810 individuals. The association between methylation of four CpG sites and the composition of plasma and IgG glycomes was analyzed. Several statistically significant associations were observed between HNF1A methylation and plasma glycans, while there were no significant associations with IgG glycans. The most consistent association with HNF1A methylation was observed with the increase in the proportion of highly branched glycans in the plasma N-glycome. The hypothesis that inactivation of HNF1A promotes branching of glycans was supported by the analysis of plasma N-glycomes in 61 patients with inactivating mutations in HNF1A, where the increase in plasma glycan branching was also observed. This study represents the first demonstration of epigenetic regulation of plasma glycome composition, suggesting a potential mechanism by which epigenetic deregulation of the glycome may contribute to disease development.

Glycosylation of immunoglobulin G (IgG) influences IgG effector function by modulating binding to Fc receptors. To identify genetic loci associated with IgG glycosylation, we quantitated N-linked IgG glycans using two approaches. After isolating IgG from human plasma, we performed 77 quantitative measurements of N-glycosylation using ultra-performance liquid chromatography (UPLC) in 2,247 individuals from four European discovery populations. In parallel, we measured IgG N-glycans using MALDI-TOF mass spectrometry (MS) in a replication cohort of 1,848 Europeans. Meta-analysis of genome-wide association study (GWAS) results identified 9 genome-wide significant loci (P<2.27×10−9) in the discovery analysis and two of the same loci (B4GALT1 and MGAT3) in the replication cohort. Four loci contained genes encoding glycosyltransferases (ST6GAL1, B4GALT1, FUT8, and MGAT3), while the remaining 5 contained genes that have not been previously implicated in protein glycosylation (IKZF1, IL6ST-ANKRD55, ABCF2-SMARCD3, SUV420H1, and SMARCB1-DERL3). However, most of them have been strongly associated with autoimmune and inflammatory conditions (e.g., systemic lupus erythematosus, rheumatoid arthritis, ulcerative colitis, Crohn's disease, diabetes type 1, multiple sclerosis, Graves' disease, celiac disease, nodular sclerosis) and/or haematological cancers (acute lymphoblastic leukaemia, Hodgkin lymphoma, and multiple myeloma). Follow-up functional experiments in haplodeficient Ikzf1 knock-out mice showed the same general pattern of changes in IgG glycosylation as identified in the meta-analysis. As IKZF1 was associated with multiple IgG N-glycan traits, we explored biomarker potential of affected N-glycans in 101 cases with SLE and 183 matched controls and demonstrated substantial discriminative power in a ROC-curve analysis (area under the curve = 0.842). Our study shows that it is possible to identify new loci that control glycosylation of a single plasma protein using GWAS. The results may also provide an explanation for the reported pleiotropy and antagonistic effects of loci involved in autoimmune diseases and haematological cancer.

Author Summary

After analysing glycans attached to human immunoglobulin G in 4,095 individuals, we performed the first genome-wide association study (GWAS) of the glycome of an individual protein. Nine genetic loci were found to associate with glycans with genome-wide significance. Of these, four were enzymes that directly participate in IgG glycosylation, thus the observed associations were biologically founded. The remaining five genetic loci were not previously implicated in protein glycosylation, but the most of them have been reported to be relevant for autoimmune and inflammatory conditions and/or haematological cancers. A particularly interesting gene, IKZF1 was found to be associated with multiple IgG N-glycans. This gene has been implicated in numerous diseases, including systemic lupus erythematosus (SLE). We analysed N-glycans in 101 cases with SLE and 183 matched controls and demonstrated their substantial biomarker potential. Our study shows that it is possible to identify new loci that control glycosylation of a single plasma protein using GWAS. Our results may also provide an explanation for opposite effects of some genes in autoimmune diseases and haematological cancer.

We have used targeted genomic sequencing of high-complexity DNA pools based on long-range PCR and deep DNA sequencing by the SOLiD technology. The method was used for sequencing of 286 kb from four chromosomal regions with quantitative trait loci (QTL) influencing blood plasma lipid and uric acid levels in DNA pools of 500 individuals from each of five European populations. The method shows very good precision in estimating allele frequencies as compared with individual genotyping of SNPs (r2=0.95, P<10−16). Validation shows that the method is able to identify novel SNPs and estimate their frequency in high-complexity DNA pools. In our five populations, 17% of all SNPs and 61% of structural variants are not available in the public databases. A large fraction of the novel variants show a limited geographic distribution, with 62% of the novel SNPs and 59% of novel structural variants being detected in only one of the populations. The large number of population-specific novel SNPs underscores the need for comprehensive sequencing of local populations in order to identify the causal variants of human traits.

The limited proportion of complex trait variance identified in genome-wide association studies may reflect the limited power of single SNP analyses to detect either rare causative alleles or those of small effect. Motivated by studies that demonstrate that loci contributing to trait variation may contain a number of different alleles, we have developed an analytical approach termed Regional Genomic Relationship Mapping that, like linkage-based family methods, integrates variance contributed by founder gametes within a pedigree. This approach takes advantage of very distant (and unrecorded) relationships, and this greatly increases the power of the method, compared with traditional pedigree-based linkage analyses. By integrating variance contributed by founder gametes in the population, our approach provides an estimate of the Regional Heritability attributable to a small genomic region (e.g. 100 SNP window covering ca. 1 Mb of DNA in a 300000 SNP GWAS) and has the power to detect regions containing multiple alleles that individually contribute too little variance to be detectable by GWAS as well as regions with single common GWAS-detectable SNPs. We use genome-wide SNP array data to obtain both a genome-wide relationship matrix and regional relationship (“identity by state" or IBS) matrices for sequential regions across the genome. We then estimate a heritability for each region sequentially in our genome-wide scan. We demonstrate by simulation and with real data that, when compared to traditional (“individual SNP") GWAS, our method uncovers new loci that explain additional trait variation. We analysed data from three Southern European populations and from Orkney for exemplar traits – serum uric acid concentration and height. We show that regional heritability estimates are correlated with results from genome-wide association analysis but can capture more of the genetic variance segregating in the population and identify additional trait loci.

Genome-wide association studies (GWAS) aim to detect single nucleotide polymorphisms (SNP) associated with trait variation. However, due to the large number of tests, standard analysis techniques impose highly stringent significance thresholds, leaving potentially associated SNPs undetected, and much of the trait genetic variation unexplained. Pathway- and network-based methodologies applied to GWAS aim to detect associations missed by standard single-marker approaches. The complex and non-random architecture of the genome makes it a challenge to derive an appropriate testing framework for such methodologies. We developed a rapid and simple permutation approach that uses GWAS SNP association results to establish the significance of pathway associations while accounting for the linkage disequilibrium structure of SNPs and the clustering of functionally related elements in the genome. All SNPs used in the GWAS are placed in a “circular genome” according to their location. Then the complete set of SNP association P values are permuted by rotation with respect to the genomic locations of the SNPs. Once these “simulated” P values are assigned, the joint gene P values are calculated using Fisher’s combination test, and the association of pathways is tested using the hypergeometric test. The circular genomic permutation approach was applied to a human genome-wide association dataset. The data consists of 719 individuals from the ORCADES study genotyped for ∼300,000 SNPs and measured for 51 traits ranging from physical to biochemical measurements. KEGG pathways (n = 225) were used as the sets of pathways to be tested. Our results demonstrate that the circular genomic permutations provide robust association P values. The non-permuted hypergeometric analysis generates ∼1400 pathway-trait combination results with an association P value more significant than P ≤ 0.05, whereas applying circular genomic permutation reduces the number of significant results to a more credible 40% of that value. The circular permutation software (“genomicper”) is available as an R package at http://cran.r-project.org/.

We surveyed gene–gene interactions (epistasis) in human body mass index (BMI) in four European populations (n<1200) via exhaustive pair-wise genome scans where interactions were computed as F ratios by testing a linear regression model fitting two single-nucleotide polymorphisms (SNPs) with interactions against the one without. Before the association tests, BMI was corrected for sex and age, normalised and adjusted for relatedness. Neither single SNPs nor SNP interactions were genome-wide significant in either cohort based on the consensus threshold (P=5.0E−08) and a Bonferroni corrected threshold (P=1.1E−12), respectively. Next we compared sub genome-wide significant SNP interactions (P<5.0E−08) across cohorts to identify common epistatic signals, where SNPs were annotated to genes to test for gene ontology (GO) enrichment. Among the epistatic genes contributing to the commonly enriched GO terms, 19 were shared across study cohorts of which 15 are previously published genome-wide association loci, including CDH13 (cadherin 13) associated with height and SORCS2 (sortilin-related VPS10 domain containing receptor 2) associated with circulating insulin-like growth factor 1 and binding protein 3. Interactions between the 19 shared epistatic genes and those involving BMI candidate loci (P<5.0E−08) were tested across cohorts and found eight replicated at the SNP level (P<0.05) in at least one cohort, which were further tested and showed limited replication in a separate European population (n>5000). We conclude that genome-wide analysis of epistasis in multiple populations is an effective approach to provide new insights into the genetic regulation of BMI but requires additional efforts to confirm the findings.

Stature is a classical and highly heritable complex trait, with 80%–90% of variation explained by genetic factors. In recent years, genome-wide association studies (GWAS) have successfully identified many common additive variants influencing human height; however, little attention has been given to the potential role of recessive genetic effects. Here, we investigated genome-wide recessive effects by an analysis of inbreeding depression on adult height in over 35,000 people from 21 different population samples. We found a highly significant inverse association between height and genome-wide homozygosity, equivalent to a height reduction of up to 3 cm in the offspring of first cousins compared with the offspring of unrelated individuals, an effect which remained after controlling for the effects of socio-economic status, an important confounder (χ2 = 83.89, df = 1; p = 5.2×10−20). There was, however, a high degree of heterogeneity among populations: whereas the direction of the effect was consistent across most population samples, the effect size differed significantly among populations. It is likely that this reflects true biological heterogeneity: whether or not an effect can be observed will depend on both the variance in homozygosity in the population and the chance inheritance of individual recessive genotypes. These results predict that multiple, rare, recessive variants influence human height. Although this exploratory work focuses on height alone, the methodology developed is generally applicable to heritable quantitative traits (QT), paving the way for an investigation into inbreeding effects, and therefore genetic architecture, on a range of QT of biomedical importance.

Author Summary

Studies investigating the extent to which genetics influences human characteristics such as height have concentrated mainly on common variants of genes, where having one or two copies of a given variant influences the trait or risk of disease. This study explores whether a different type of genetic variant might also be important. We investigate the role of recessive genetic variants, where two identical copies of a variant are required to have an effect. By measuring genome-wide homozygosity—the phenomenon of inheriting two identical copies at a given point of the genome—in 35,000 individuals from 21 European populations, and by comparing this to individual height, we found that the more homozygous the genome, the shorter the individual. The offspring of first cousins (who have increased homozygosity) were predicted to be up to 3 cm shorter on average than the offspring of unrelated parents. Height is influenced by the combined effect of many recessive variants dispersed across the genome. This may also be true for other human characteristics and diseases, opening up a new way to understand how genetic variation influences our health.

Serum concentrations of low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), triglycerides (TGs) and total cholesterol (TC) are important heritable risk factors for cardiovascular disease. Although genome-wide association studies (GWASs) of circulating lipid levels have identified numerous loci, a substantial portion of the heritability of these traits remains unexplained. Evidence of unexplained genetic variance can be detected by combining multiple independent markers into additive genetic risk scores. Such polygenic scores, constructed using results from the ENGAGE Consortium GWAS on serum lipids, were applied to predict lipid levels in an independent population-based study, the Rotterdam Study-II (RS-II). We additionally tested for evidence of a shared genetic basis for different lipid phenotypes. Finally, the polygenic score approach was used to identify an alternative genome-wide significance threshold before pathway analysis and those results were compared with those based on the classical genome-wide significance threshold. Our study provides evidence suggesting that many loci influencing circulating lipid levels remain undiscovered. Cross-prediction models suggested a small overlap between the polygenic backgrounds involved in determining LDL-C, HDL-C and TG levels. Pathway analysis utilizing the best polygenic score for TC uncovered extra information compared with using only genome-wide significant loci. These results suggest that the genetic architecture of circulating lipids involves a number of undiscovered variants with very small effects, and that increasing GWAS sample sizes will enable the identification of novel variants that regulate lipid levels.

Variation in the apolipoprotein E gene (APOE) has been reported to be associated with longevity in humans. The authors assessed the allelic distribution of APOE isoforms ε2, ε3, and ε4 among 10,623 participants from 15 case-control and cohort studies of age-related macular degeneration (AMD) in populations of European ancestry (study dates ranged from 1990 to 2009). The authors included only the 10,623 control subjects from these studies who were classified as having no evidence of AMD, since variation within the APOE gene has previously been associated with AMD. In an analysis stratified by study center, gender, and smoking status, there was a decreasing frequency of the APOE ε4 isoform with increasing age (χ2 for trend = 14.9 (1 df); P = 0.0001), with a concomitant increase in the ε3 isoform (χ2 for trend = 11.3 (1 df); P = 0.001). The association with age was strongest in ε4 homozygotes; the frequency of ε4 homozygosity decreased from 2.7% for participants aged 60 years or less to 0.8% for those over age 85 years, while the proportion of participants with the ε3/ε4 genotype decreased from 26.8% to 17.5% across the same age range. Gender had no significant effect on the isoform frequencies. This study provides strong support for an association of the APOE gene with human longevity.

Myopia is a complex genetic disorder and a common cause of visual impairment among working age adults. Genome-wide association studies have identified susceptibility loci on chromosomes 15q14 and 15q25 in Caucasian populations of European ancestry. Here, we present a confirmation and meta-analysis study in which we assessed whether these two loci are also associated with myopia in other populations. The study population comprised 31 cohorts from the Consortium of Refractive Error and Myopia (CREAM) representing 4 different continents with 55,177 individuals; 42,845 Caucasians and 12,332 Asians. We performed a meta-analysis of 14 single nucleotide polymorphisms (SNPs) on 15q14 and 5 SNPs on 15q25 using linear regression analysis with spherical equivalent as a quantitative outcome, adjusted for age and sex. We calculated the odds ratio (OR) of myopia versus hyperopia for carriers of the top-SNP alleles using a fixed effects meta-analysis. At locus 15q14, all SNPs were significantly replicated, with the lowest P value 3.87 × 10−12 for SNP rs634990 in Caucasians, and 9.65 × 10−4 for rs8032019 in Asians. The overall meta-analysis provided P value 9.20 × 10−23 for the top SNP rs634990. The risk of myopia versus hyperopia was OR 1.88 (95 % CI 1.64, 2.16, P < 0.001) for homozygous carriers of the risk allele at the top SNP rs634990, and OR 1.33 (95 % CI 1.19, 1.49, P < 0.001) for heterozygous carriers. SNPs at locus 15q25 did not replicate significantly (P value 5.81 × 10−2 for top SNP rs939661). We conclude that common variants at chromosome 15q14 influence susceptibility for myopia in Caucasian and Asian populations world-wide.

Electronic supplementary material

The online version of this article (doi:10.1007/s00439-012-1176-0) contains supplementary material, which is available to authorized users.

Complement factor D catalyzes a critical step in the alternative complement activation pathway. The authors report a significant elevation in plasma CFD concentrations in age-related macular degeneration (AMD) patients compared with controls and a weak genetic association between CFD gene variants and AMD.

Purpose.

To examine the role of complement factor D (CFD) in age-related macular degeneration (AMD) by analysis of genetic association, copy number variation, and plasma CFD concentrations.

Methods.

Single nucleotide polymorphisms (SNPs) in the CFD gene were genotyped and the results analyzed by binary logistic regression. CFD gene copy number was analyzed by gene copy number assay. Plasma CFD was measured by an enzyme-linked immunosorbent assay.

Results.

Genetic association was found between CFD gene SNP rs3826945 and AMD (odds ratio 1.44; P = 0.028) in a small discovery case-control series (462 cases and 325 controls) and replicated in a combined cohorts meta-analysis of 4765 cases and 2693 controls, with an odds ratio of 1.11 (P = 0.032), with the association almost confined to females. Copy number variation in the CFD gene was identified in 13 out of 640 samples examined but there was no difference in frequency between AMD cases (1.3%) and controls (2.7%). Plasma CFD concentration was measured in 751 AMD cases and 474 controls and found to be elevated in AMD cases (P = 0.00025). The odds ratio for those in the highest versus lowest quartile for plasma CFD was 1.81. The difference in plasma CFD was again almost confined to females.

Conclusions.

CFD regulates activation of the alternative complement pathway, which is implicated in AMD pathogenesis. The authors found evidence for genetic association between a CFD gene SNP and AMD and a significant increase in plasma CFD concentration in AMD cases compared with controls, consistent with a role for CFD in AMD pathogenesis.

QRS interval on the electrocardiogram reflects ventricular depolarization and conduction time, and is a risk factor for mortality, sudden death, and heart failure. We performed a genome-wide association meta-analysis in 40,407 European-descent individuals from 14 studies, with further genotyping in 7170 additional Europeans, and identified 22 loci associated with QRS duration (P < 5 × 10−8). These loci map in or near genes in pathways with established roles in ventricular conduction such as sodium channels, transcription factors, and calcium-handling proteins, but also point to novel biologic processes, such as kinase inhibitors and genes related to tumorigenesis. We demonstrate that SCN10A, a gene at our most significant locus, is expressed in the mouse ventricular conduction system, and treatment with a selective SCN10A blocker prolongs QRS duration. These findings extend our current knowledge of ventricular depolarization and conduction.

Brachial circumference (BC), also known as upper arm or mid arm circumference, can be used as an indicator of muscle mass and fat tissue, which are distributed differently in men and women. Analysis of anthropometric measures of peripheral fat distribution such as BC could help in understanding the complex pathophysiology behind overweight and obesity. The purpose of this study is to identify genetic variants associated with BC through a large-scale genome-wide association scan (GWAS) meta-analysis. We used fixed-effects meta-analysis to synthesise summary results across 14 GWAS discovery and 4 replication cohorts comprising overall 22,376 individuals (12,031 women and 10,345 men) of European ancestry. Individual analyses were carried out for men, women, and combined across sexes using linear regression and an additive genetic model: adjusted for age and adjusted for age and BMI. We prioritised signals for follow-up in two-stages. We did not detect any signals reaching genome-wide significance. The FTO rs9939609 SNP showed nominal evidence for association (p<0.05) in the age-adjusted strata for men and across both sexes. In this first GWAS meta-analysis for BC to date, we have not identified any genome-wide significant signals and do not observe robust association of previously established obesity loci with BC. Large-scale collaborations will be necessary to achieve higher power to detect loci underlying BC.

The level of body iron storage and the erythropoietic need for iron are indicated by the serum levels of ferritin and soluble transferrin receptor (sTfR), respectively. A meta-analysis of five genome-wide association studies on sTfR and ferritin revealed novel association to the PCSK7 and TMPRSS6 loci for sTfR and the HFE locus for both parameters. The PCSK7 association was the most significant (rs236918, P = 1.1 × 10E−27) suggesting that proprotein convertase 7, the gene product of PCSK7, may be involved in sTfR generation and/or iron homeostasis. Conditioning the sTfR analyses on transferrin saturation abolished the HFE signal and substantially diminished the TMPRSS6 signal while the PCSK7 association was unaffected, suggesting that the former may be mediated by transferrin saturation whereas the PCSK7-associated effect on sTfR generation appears to be more direct.

Genome-wide association studies have identified 32 loci associated with body mass index (BMI), a measure that does not allow distinguishing lean from fat mass. To identify adiposity loci, we meta-analyzed associations between ~2.5 million SNPs and body fat percentage from 36,626 individuals, and followed up the 14 most significant (P<10−6) independent loci in 39,576 individuals. We confirmed the previously established adiposity locus in FTO (P=3×10−26), and identified two new loci associated with body fat percentage, one near IRS1 (P=4×10−11) and one near SPRY2 (P=3×10−8). Both loci harbour genes with a potential link to adipocyte physiology, of which the locus near IRS1 shows an intriguing association pattern. The body-fat-decreasing allele associates with decreased IRS1 expression and with an impaired metabolic profile, including decreased subcutaneous-to-visceral fat ratio, increased insulin resistance, dyslipidemia, risk of diabetes and coronary artery disease, and decreased adiponectin levels. Our findings provide new insights into adiposity and insulin resistance.

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P=5.4×10−60) and 9q31.2 (P=2.2×10−33), we identified 30 novel menarche loci (all P<5×10−8) and found suggestive evidence for a further 10 loci (P<1.9×10−6). New loci included four previously associated with BMI (in/near FTO, SEC16B, TRA2B and TMEM18), three in/near other genes implicated in energy homeostasis (BSX, CRTC1, and MCHR2), and three in/near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and MAGENTA pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing.

The genetic structure of human populations is important in population genetics, forensics and medicine. Using genome-wide scans and individuals with all four grandparents born in the same settlement, we here demonstrate remarkable geographical structure across 8–30 km in three different parts of rural Europe. After excluding close kin and inbreeding, village of origin could still be predicted correctly on the basis of genetic data for 89–100% of individuals.

Genome-wide association (GWA) studies have identified a number of loci underlying variation in human serum uric acid (SUA) levels with the SLC2A9 gene having the largest effect identified so far. Gene-gene interactions (epistasis) are largely unexplored in these GWA studies. We performed a full pair-wise genome scan in the Italian MICROS population (n = 1201) to characterise epistasis signals in SUA levels. In the resultant epistasis profile, no SNP pairs reached the Bonferroni adjusted threshold for the pair-wise genome-wide significance. However, SLC2A9 was found interacting with multiple loci across the genome, with NFIA - SLC2A9 and SLC2A9 - ESRRAP2 being significant based on a threshold derived for interactions between GWA significant SNPs and the genome and jointly explaining 8.0% of the phenotypic variance in SUA levels (3.4% by interaction components). Epistasis signal replication in a CROATIAN population (n = 1772) was limited at the SNP level but improved dramatically at the gene ontology level. In addition, gene ontology terms enriched by the epistasis signals in each population support links between SUA levels and neurological disorders. We conclude that GWA epistasis analysis is useful despite relatively low power in small isolated populations.

Genome analysis provides a powerful approach to test for evidence of genetic variation within and between geographical regions and local populations. Copy number variants which comprise insertions, deletions and duplications of genomic sequence provide one such convenient and informative source. Here, we investigate copy number variants from genome wide scans of single nucleotide polymorphisms in three European population isolates, the island of Vis in Croatia, the islands of Orkney in Scotland and the South Tyrol in Italy. We show that whereas the overall copy number variant frequencies are similar between populations, their distribution is highly specific to the population of origin, a finding which is supported by evidence for increased kinship correlation for specific copy number variants within populations.

Background

Plasma levels of coagulation factors VII (FVII), VIII (FVIII), and von Willebrand factor (vWF) influence risk of hemorrhage and thrombosis. We conducted genome-wide association studies to identify new loci associated with plasma levels.

Methods and Results

Setting includes 5 community-based studies for discovery comprising 23,608 European-ancestry participants: ARIC, CHS, B58C, FHS, and RS. All had genome-wide single nucleotide polymorphism (SNP) scans and at least 1 phenotype measured: FVII activity/antigen, FVIII activity, and vWF antigen. Each study used its genotype data to impute to HapMap SNPs and independently conducted association analyses of hemostasis measures using an additive genetic model. Study findings were combined by meta-analysis. Replication was conducted in 7,604 participants not in the discovery cohort. For FVII, 305 SNPs exceeded the genome-wide significance threshold of 5.0×10-8 and comprised 5 loci on 5 chromosomes: 2p23 (smallest p-value 6.2×10-24), 4q25 (3.6×10-12), 11q12 (2.0×10-10), 13q34 (9.0×10-259), and 20q11.2 (5.7×10-37). Loci were within or near genes, including 4 new candidate genes and F7 (13q34). For vWF, 400 SNPs exceeded the threshold and marked 8 loci on 6 chromosomes: 6q24 (1.2×10-22), 8p21 (1.3×10-16), 9q34 (<5.0×10-324), 12p13 (1.7×10-32), 12q23 (7.3×10-10), 12q24.3 (3.8×10-11), 14q32 (2.3×10-10) and 19p13.2 (1.3×10-9). All loci were within genes, including 6 new candidate genes, as well as ABO (9q34) and VWF (12p13). For FVIII, 5 loci were identified and overlapped vWF findings. Nine of the 10 new findings replicated.

Conclusions

New genetic associations were discovered outside previously known biologic pathways and may point to novel prevention and treatment targets of hemostasis disorders.

Over half of all proteins are glycosylated, and alterations in glycosylation have been observed in numerous physiological and pathological processes. Attached glycans significantly affect protein function; but, contrary to polypeptides, they are not directly encoded by genes, and the complex processes that regulate their assembly are poorly understood. A novel approach combining genome-wide association and high-throughput glycomics analysis of 2,705 individuals in three population cohorts showed that common variants in the Hepatocyte Nuclear Factor 1α (HNF1α) and fucosyltransferase genes FUT6 and FUT8 influence N-glycan levels in human plasma. We show that HNF1α and its downstream target HNF4α regulate the expression of key fucosyltransferase and fucose biosynthesis genes. Moreover, we show that HNF1α is both necessary and sufficient to drive the expression of these genes in hepatic cells. These results reveal a new role for HNF1α as a master transcriptional regulator of multiple stages in the fucosylation process. This mechanism has implications for the regulation of immunity, embryonic development, and protein folding, as well as for our understanding of the molecular mechanisms underlying cancer, coronary heart disease, and metabolic and inflammatory disorders.

Author Summary

By combining recently developed high-throughput glycan analysis with genome-wide association study, we performed the first comprehensive analysis of common genetic polymorphisms that affect protein glycosylation. Over half of all proteins are glycosylated; but, due to difficulties in glycan analysis and the absence of a genetic template for their synthesis, knowledge about the complex processes that regulate glycan assembly is still limited. We demonstrated that HNF1α regulates the expression of key fucosyltransferase and fucose biosynthesis genes and acts as a master regulator of plasma protein fucosylation. Proper protein fucosylation is essential in numerous processes including inflammation, cancer, and coronary heart disease, thus the identification of a master regulator of plasma protein fucosylation has important implications for understanding both normal biological functions and disease processes.

Chronic kidney disease (CKD) is a significant public health problem, and recent genetic studies have identified common CKD susceptibility variants. The CKDGen consortium performed a meta-analysis of genome-wide association data in 67,093 Caucasian individuals from 20 population-based studies to identify new susceptibility loci for reduced renal function, estimated by serum creatinine (eGFRcrea), cystatin C (eGFRcys), and CKD (eGFRcrea <60 ml/min/1.73m2; n = 5,807 CKD cases). Follow-up of the 23 genome-wide significant loci (p<5×10−8) in 22,982 replication samples identified 13 novel loci for renal function and CKD (in or near LASS2, GCKR, ALMS1, TFDP2, DAB2, SLC34A1, VEGFA, PRKAG2, PIP5K1B, ATXN2, DACH1, UBE2Q2, and SLC7A9) and 7 creatinine production and secretion loci (CPS1, SLC22A2, TMEM60, WDR37, SLC6A13, WDR72, BCAS3). These results further our understanding of biologic mechanisms of kidney function by identifying loci potentially influencing nephrogenesis, podocyte function, angiogenesis, solute transport, and metabolic functions of the kidney.

Background

Plasma levels of coagulation factors VII (FVII), VIII (FVIII), and von Willebrand factor (vWF) influence risk of hemorrhage and thrombosis. We conducted genome-wide association studies to identify new loci associated with plasma levels.

Methods and Results

Setting includes 5 community-based studies for discovery comprising 23,608 European-ancestry participants: ARIC, CHS, B58C, FHS, and RS. All had genome-wide single nucleotide polymorphism (SNP) scans and at least 1 phenotype measured: FVII activity/antigen, FVIII activity, and vWF antigen. Each study used its genotype data to impute to HapMap SNPs and independently conducted association analyses of hemostasis measures using an additive genetic model. Study findings were combined by meta-analysis. Replication was conducted in 7,604 participants not in the discovery cohort. For FVII, 305 SNPs exceeded the genome-wide significance threshold of 5.0×10−8 and comprised 5 loci on 5 chromosomes: 2p23 (smallest p-value 6.2×10−24), 4q25 (3.6×10−12), 11q12 (2.0×10−10), 13q34 (9.0×10−259), and 20q11.2 (5.7×10−37). Loci were within or near genes, including 4 new candidate genes and F7 (13q34). For vWF, 400 SNPs exceeded the threshold and marked 8 loci on 6 chromosomes: 6q24 (1.2×10−22), 8p21 (1.3×10−16), 9q34 (<5.0×10−324), 12p13 (1.7×10−32), 12q23 (7.3×10−10), 12q24.3 (3.8×10−11), 14q32 (2.3×10−10) and 19p13.2 (1.3×10−9). All loci were within genes, including 6 new candidate genes, as well as ABO (9q34) and VWF (12p13). For FVIII, 5 loci were identified and overlapped vWF findings. Nine of the 10 new findings replicated.

Conclusions

New genetic associations were discovered outside previously known biologic pathways and may point to novel prevention and treatment targets of hemostasis disorders.