Fast and accurate recurrent event analysis for genome-wide association studies.

Many diseases recur after recovery, for example, recurrences in cancer and infections. However, research is often focused on analysing only time-to-first recurrence, thereby ignoring any subsequent recurrences that may occur after the first. Statistical models for the analysis of recurrent events are available, of which the extended Cox proportional hazards frailty model is the current state-of-the-art. However, this model is too statistically complex for computationally efficient application in high-dimensional data sets, including genome-wide association studies (GWAS). Here, we develop an application for fast and accurate recurrent event analysis in GWAS, called SPARE (SaddlePoint Approximation for Recurrent Event analysis). In SPARE, every DNA variant is tested for association with recurrence risk using a modified score statistic. A saddlepoint approximation is implemented to achieve statistical accuracy. SPARE controls the Type I error, and its statistical power is similar to existing recurrent event models, yet SPARE is significantly faster. An application of SPARE in a recurrent event GWAS on bladder cancer for 6.2 million DNA variants in 1,443 individuals required less than 15 min, whereas existing recurrent event methods would require several weeks.

Genomics Research of Lifetime Depression in the Netherlands: The BIObanks Netherlands Internet Collaboration (BIONIC) Project.

In this cohort profile article we describe the lifetime major depressive disorder (MDD) database that has been established as part of the BIObanks Netherlands Internet Collaboration (BIONIC). Across the Netherlands we collected data on Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) lifetime MDD diagnosis in 132,850 Dutch individuals. Currently, N = 66,684 of these also have genomewide single nucleotide polymorphism (SNP) data. We initiated this project because the complex genetic basis of MDD requires large population-wide studies with uniform in-depth phenotyping. For standardized phenotyping we developed the LIDAS (LIfetime Depression Assessment Survey), which then was used to measure MDD in 11 Dutch cohorts. Data from these cohorts were combined with diagnostic interview depression data from 5 clinical cohorts to create a dataset of N = 29,650 lifetime MDD cases (22%) meeting DSM-5 criteria and 94,300 screened controls. In addition, genomewide genotype data from the cohorts were assembled into a genomewide association study (GWAS) dataset of N = 66,684 Dutch individuals (25.3% cases). Phenotype data include DSM-5-based MDD diagnoses, sociodemographic variables, information on lifestyle and BMI, characteristics of depressive symptoms and episodes, and psychiatric diagnosis and treatment history. We describe the establishment and harmonization of the BIONIC phenotype and GWAS datasets and provide an overview of the available information and sample characteristics. Our next step is the GWAS of lifetime MDD in the Netherlands, with future plans including fine-grained genetic analyses of depression characteristics, international collaborations and multi-omics studies.

Genome-wide association study identifies 74 loci associated with educational attainment.

Educational attainment is strongly influenced by social and other environmental factors, but genetic factors are estimated to account for at least 20% of the variation across individuals. Here we report the results of a genome-wide association study (GWAS) for educational attainment that extends our earlier discovery sample of 101,069 individuals to 293,723 individuals, and a replication study in an independent sample of 111,349 individuals from the UK Biobank. We identify 74 genome-wide significant loci associated with the number of years of schooling completed. Single-nucleotide polymorphisms associated with educational attainment are disproportionately found in genomic regions regulating gene expression in the fetal brain. Candidate genes are preferentially expressed in neural tissue, especially during the prenatal period, and enriched for biological pathways involved in neural development. Our findings demonstrate that, even for a behavioural phenotype that is mostly environmentally determined, a well-powered GWAS identifies replicable associated genetic variants that suggest biologically relevant pathways. Because educational attainment is measured in large numbers of individuals, it will continue to be useful as a proxy phenotype in efforts to characterize the genetic influences of related phenotypes, including cognition and neuropsychiatric diseases.

Transferrin Saturation/Hepcidin Ratio Discriminates TMPRSS6-Related Iron Refractory Iron Deficiency Anemia from Patients with Multi-Causal Iron Deficiency Anemia.

Pathogenic TMPRSS6 variants impairing matriptase-2 function result in inappropriately high hepcidin levels relative to body iron status, leading to iron refractory iron deficiency anemia (IRIDA). As diagnosing IRIDA can be challenging due to its genotypical and phenotypical heterogeneity, we assessed the transferrin saturation (TSAT)/hepcidin ratio to distinguish IRIDA from multi-causal iron deficiency anemia (IDA). We included 20 IRIDA patients from a registry for rare inherited iron disorders and then enrolled 39 controls with IDA due to other causes. Plasma hepcidin-25 levels were measured by standardized isotope dilution mass spectrometry. IDA controls had not received iron therapy in the last 3 months and C-reactive protein levels were <10.0 mg/L. IRIDA patients had significantly lower TSAT/hepcidin ratios compared to IDA controls, median 0.6%/nM (interquartile range, IQR, 0.4-1.1%/nM) and 16.7%/nM (IQR, 12.0-24.0%/nM), respectively. The area under the curve for the TSAT/hepcidin ratio was 1.000 with 100% sensitivity and specificity (95% confidence intervals 84-100% and 91-100%, respectively) at an optimal cut-off point of 5.6%/nM. The TSAT/hepcidin ratio shows excellent performance in discriminating IRIDA from TMPRSS6-unrelated IDA early in the diagnostic work-up of IDA provided that recent iron therapy and moderate-to-severe inflammation are absent. These observations warrant further exploration in a broader IDA population.

Identification of ADHD risk genes in extended pedigrees by combining linkage analysis and whole-exome sequencing.

Attention-deficit/hyperactivity disorder (ADHD) is a common neurodevelopmental disorder with a complex genetic background, hampering identification of underlying genetic risk factors. We hypothesized that combining linkage analysis and whole-exome sequencing (WES) in multi-generation pedigrees with multiple affected individuals can point toward novel ADHD genes. Three families with multiple ADHD-affected members (Ntotal = 70) and apparent dominant inheritance pattern were included in this study. Genotyping was performed in 37 family members, and WES was additionally carried out in 10 of those. Linkage analysis was performed using multi-point analysis in Superlink Online SNP 1.1. From prioritized linkage regions with a LOD score ≥ 2, a total of 24 genes harboring rare variants were selected. Those genes were taken forward and were jointly analyzed in gene-set analyses of exome-chip data using the MAGMA software in an independent sample of patients with persistent ADHD and healthy controls (N = 9365). The gene-set including all 24 genes together, and particularly the gene-set from one of the three families (12 genes), were significantly associated with persistent ADHD in this sample. Among the latter, gene-wide analysis for the AAED1 gene reached significance. A rare variant (rs151326868) within AAED1 segregated with ADHD in one of the families. The analytic strategy followed here is an effective approach for identifying novel ADHD risk genes. Additionally, this study suggests that both rare and more frequent variants in multiple genes act together in contributing to ADHD risk, even in individual multi-case families.

Sex-Specific Regulation of Inflammation and Metabolic Syndrome in Obesity.

OBJECTIVE: Metabolic dysregulation and inflammation are important consequences of obesity and impact susceptibility to cardiovascular disease. Anti-inflammatory therapy in cardiovascular disease is being developed under the assumption that inflammatory pathways are identical in women and men, but it is not known if this is indeed the case. In this study, we assessed the sex-specific relation between inflammation and metabolic dysregulation in obesity. Approach and Results: Three hundred two individuals were included, half with a BMI 27 to 30 kg/m2 and half with a BMI>30 kg/m2, 45% were women. The presence of metabolic syndrome was assessed according to the National Cholesterol Education Program-ATPIII criteria, and inflammation was studied using circulating markers of inflammation, cell counts, and ex vivo cytokine production capacity of isolated immune cells. Additionally, lipidomic and metabolomic data were gathered, and subcutaneous fat biopsies were histologically assessed. Metabolic syndrome is associated with an increased inflammatory profile that profoundly differs between women and men: women with metabolic syndrome show a lower concentration of the anti-inflammatory adiponectin, whereas men show increased levels of several pro-inflammatory markers such as IL (interleukin)-6 and leptin. Adipose tissue inflammation showed similar sex-specific associations with these markers. Peripheral blood mononuclear cells isolated from men, but not women, with metabolic syndrome display enhanced cytokine production capacity. CONCLUSIONS: We identified sex-specific pathways that influence inflammation in obesity. Excessive production of proinflammatory cytokines was observed in men with metabolic syndrome. In contrast, women typically showed reduced levels of the anti-inflammatory adipokine adiponectin. These different mechanisms of inflammatory dysregulation between women and men with obesity argue for sex-specific therapeutic strategies.

Integrating Metabolomics, Genomics, and Disease Pathways in Age-Related Macular Degeneration: The EYE-RISK Consortium.

PURPOSE: The current study aimed to identify metabolites associated with age-related macular degeneration (AMD) by performing the largest metabolome association analysis in AMD to date, as well as aiming to determine the effect of AMD-associated genetic variants on metabolite levels and investigate associations between the identified metabolites and activity of the complement system, one of the main AMD-associated disease pathways. DESIGN: Case-control association analysis of metabolomics data. PARTICIPANTS: Five European cohorts consisting of 2267 AMD patients and 4266 control participants. METHODS: Metabolomics was performed using a high-throughput proton nuclear magnetic resonance metabolomics platform, which allows quantification of 146 metabolite measurements and 79 derivative values. Metabolome-AMD associations were studied using univariate logistic regression analyses. The effect of 52 AMD-associated genetic variants on the identified metabolites was investigated using linear regression. In addition, associations between the identified metabolites and activity of the complement pathway (defined by the C3d-to-C3 ratio) were investigated using linear regression. MAIN OUTCOME MEASURES: Metabolites associated with AMD. RESULTS: We identified 60 metabolites that were associated significantly with AMD, including increased levels of large and extra-large high-density lipoprotein (HDL) subclasses and decreased levels of very low-density lipoprotein (VLDL), amino acids, and citrate. Of 52 AMD-associated genetic variants, 7 variants were associated significantly with 34 of the identified metabolites. The strongest associations were identified for genetic variants located in or near genes involved in lipid metabolism (ABCA1, CETP, APOE, and LIPC) with metabolites belonging to the large and extra-large HDL subclasses. Also, 57 of 60 metabolites were associated significantly with complement activation levels, independent of AMD status. Increased large and extra-large HDL levels and decreased VLDL and amino acid levels were associated with increased complement activation. CONCLUSIONS: Lipoprotein levels were associated with AMD-associated genetic variants, whereas decreased essential amino acids may point to nutritional deficiencies in AMD. We observed strong associations between the vast majority of the AMD-associated metabolites and systemic complement activation levels, independent of AMD status. This may indicate biological interactions between the main AMD disease pathways and suggests that multiple pathways may need to be targeted simultaneously for successful treatment of AMD.

Measurement and genetic architecture of lifetime depression in the Netherlands as assessed by LIDAS (Lifetime Depression Assessment Self-report).

BACKGROUND: Major depressive disorder (MDD) is a common mood disorder, with a heritability of around 34%. Molecular genetic studies made significant progress and identified genetic markers associated with the risk of MDD; however, progress is slowed down by substantial heterogeneity as MDD is assessed differently across international cohorts. Here, we used a standardized online approach to measure MDD in multiple cohorts in the Netherlands and evaluated whether this approach can be used in epidemiological and genetic association studies of depression. METHODS: Within the Biobank Netherlands Internet Collaboration (BIONIC) project, we collected MDD data in eight cohorts involving 31 936 participants, using the online Lifetime Depression Assessment Self-report (LIDAS), and estimated the prevalence of current and lifetime MDD in 22 623 unrelated individuals. In a large Netherlands Twin Register (NTR) twin-family dataset (n ≈ 18 000), we estimated the heritability of MDD, and the prediction of MDD in a subset (n = 4782) through Polygenic Risk Score (PRS). RESULTS: Estimates of current and lifetime MDD prevalence were 6.7% and 18.1%, respectively, in line with population estimates based on validated psychiatric interviews. In the NTR heritability estimates were 0.34/0.30 (s.e. = 0.02/0.02) for current/lifetime MDD, respectively, showing that the LIDAS gives similar heritability rates for MDD as reported in the literature. The PRS predicted risk of MDD (OR 1.23, 95% CI 1.15-1.32, R2 = 1.47%). CONCLUSIONS: By assessing MDD status in the Netherlands using the LIDAS instrument, we were able to confirm previously reported MDD prevalence and heritability estimates, which suggests that this instrument can be used in epidemiological and genetic association studies of depression.

Standardized serum hepcidin values in Dutch children: Set point relative to body iron changes during childhood.

BACKGROUND: Use of serum hepcidin measurements in pediatrics would benefit from standardized age- and sex-specific reference ranges in children, in order to enable the establishment of clinical decision limits that are universally applicable. PROCEDURE: We measured serum hepcidin-25 levels in 266 healthy Dutch children aged 0.3-17 years, using an isotope dilution mass spectrometry assay, standardized with our commutable secondary reference material (RM), assigned by a candidate primary RM. RESULTS: We constructed age- and sex-specific values for serum hepcidin and its ratio with ferritin and transferrin saturation (TSAT). Serum hepcidin levels and hepcidin/ferritin and TSAT/hepcidin ratios were similar for both sexes. Serum hepcidin and hepcidin/ferritin ratio substantially declined after the age of 12 years and TSAT/hepcidin ratio gradually increased with increasing age. Serum hepcidin values for Dutch children <12 years (n = 170) and >12 years (n = 96) were 1.9 nmol/L (median); 0.1-13.1 nmol/L (p2.5-p97.5) and 0.9 nmol/L; 0.0-9.1 nmol/L, respectively. Serum ferritin was the most significant correlate of serum hepcidin in our study population, explaining 15.1% and 7.9% of variance in males and females, respectively. Multivariable linear regression analysis including age, blood sampling time, iron parameters, ALT, CRP, and body mass index as independent variables showed a statistically significant negative association between age as a dichotomous variable (≤12 vs >12 years) and log-transformed serum hepcidin levels in both sexes. CONCLUSIONS: We demonstrate that serum hepcidin relative to indicators of body iron is age dependent in children, suggesting that the set point of serum hepcidin relative to stored and circulating iron changes during childhood.

Cigarette Smoking and the Risk of Cutaneous Melanoma: A Case-Control Study.

OBJECTIVE: Previous studies regarding cigarette smoking causing a lower risk of melanoma are inconclusive. Here, we re-examined melanoma risk in relation to cigarette smoking in a large, case-control study. METHODS: In total 1,157 patients with melanoma diagnosed between 2003 and 2011 in the Netherlands and 5,595 controls from the Nijmegen Biomedical Study were included. Information concerning smoking habits and known risk factors for melanoma were obtained through self-administered questionnaires. Logistic regression analyses stratified by gender were performed to study the risk of cigarette smoking on melanoma risk, adjusted for age, marital status, highest level of education, skin type, sun vacation, use of solarium, time spent outdoors, and sun protective measures. RESULTS: Among men, current and former smokers did not have a higher risk of melanoma compared to never smokers: adjusted odds ratio (OR) = 0.56 (95% confidence interval [CI]: 0.40-0.79) and adjusted OR = 0.50 (95% CI: 0.39-0.64), respectively. With an increasing number of years smoked the risk of melanoma decreased: <20 years: OR = 0.61 (95% CI: 0.46-0.80); 21-40 years: OR = 0.50 (95% CI: 0.37-0.68); >40 years: OR = 0.26 (95% CI: 0.15-0.44). No clear trend was found for the number of cigarettes smoked. Results for females were less clear and not statistically significant (current smoker: adjusted OR = 0.96, 95% CI: 0.74-1.26, former smoker: adjusted OR = 0.89, 95% CI: 0.73-1.08). CONCLUSION: This study shows a strong inverse association between cigarette smoking and melanoma risk in men. Fundamental laboratory research is necessary to investigate the biological relation between smoking cigarettes and melanoma.

Variant ASGR1 Associated with a Reduced Risk of Coronary Artery Disease.

BACKGROUND: Several sequence variants are known to have effects on serum levels of non-high-density lipoprotein (HDL) cholesterol that alter the risk of coronary artery disease. METHODS: We sequenced the genomes of 2636 Icelanders and found variants that we then imputed into the genomes of approximately 398,000 Icelanders. We tested for association between these imputed variants and non-HDL cholesterol levels in 119,146 samples. We then performed replication testing in two populations of European descent. We assessed the effects of an implicated loss-of-function variant on the risk of coronary artery disease in 42,524 case patients and 249,414 controls from five European ancestry populations. An augmented set of genomes was screened for additional loss-of-function variants in a target gene. We evaluated the effect of an implicated variant on protein stability. RESULTS: We found a rare noncoding 12-base-pair (bp) deletion (del12) in intron 4 of ASGR1, which encodes a subunit of the asialoglycoprotein receptor, a lectin that plays a role in the homeostasis of circulating glycoproteins. The del12 mutation activates a cryptic splice site, leading to a frameshift mutation and a premature stop codon that renders a truncated protein prone to degradation. Heterozygous carriers of the mutation (1 in 120 persons in our study population) had a lower level of non-HDL cholesterol than noncarriers, a difference of 15.3 mg per deciliter (0.40 mmol per liter) (P=1.0×10(-16)), and a lower risk of coronary artery disease (by 34%; 95% confidence interval, 21 to 45; P=4.0×10(-6)). In a larger set of sequenced samples from Icelanders, we found another loss-of-function ASGR1 variant (p.W158X, carried by 1 in 1850 persons) that was also associated with lower levels of non-HDL cholesterol (P=1.8×10(-3)). CONCLUSIONS: ASGR1 haploinsufficiency was associated with reduced levels of non-HDL cholesterol and a reduced risk of coronary artery disease. (Funded by the National Institutes of Health and others.).

Identification of a novel susceptibility locus at 13q34 and refinement of the 20p12.2 region as a multi-signal locus associated with bladder cancer risk in individuals of European ancestry.

Candidate gene and genome-wide association studies (GWAS) have identified 15 independent genomic regions associated with bladder cancer risk. In search for additional susceptibility variants, we followed up on four promising single-nucleotide polymorphisms (SNPs) that had not achieved genome-wide significance in 6911 cases and 11 814 controls (rs6104690, rs4510656, rs5003154 and rs4907479, P < 1 × 10(-6)), using additional data from existing GWAS datasets and targeted genotyping for studies that did not have GWAS data. In a combined analysis, which included data on up to 15 058 cases and 286 270 controls, two SNPs achieved genome-wide statistical significance: rs6104690 in a gene desert at 20p12.2 (P = 2.19 × 10(-11)) and rs4907479 within the MCF2L gene at 13q34 (P = 3.3 × 10(-10)). Imputation and fine-mapping analyses were performed in these two regions for a subset of 5551 bladder cancer cases and 10 242 controls. Analyses at the 13q34 region suggest a single signal marked by rs4907479. In contrast, we detected two signals in the 20p12.2 region-the first signal is marked by rs6104690, and the second signal is marked by two moderately correlated SNPs (r(2) = 0.53), rs6108803 and the previously reported rs62185668. The second 20p12.2 signal is more strongly associated with the risk of muscle-invasive (T2-T4 stage) compared with non-muscle-invasive (Ta, T1 stage) bladder cancer (case-case P ≤ 0.02 for both rs62185668 and rs6108803). Functional analyses are needed to explore the biological mechanisms underlying these novel genetic associations with risk for bladder cancer.

Identification of novel genetic Loci associated with thyroid peroxidase antibodies and clinical thyroid disease.

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto's thyroiditis), as well as autoimmune hyperthyroidism (Graves' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10(-8)) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68-2.81, P = 8.1×10(-8)), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26-1.82, P = 2.9×10(-6)), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66-0.89, P = 6.5×10(-4)). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves' disease (OR: 1.37, 95% CI 1.22-1.54, P = 1.2×10(-7) and OR: 1.25, 95% CI 1.12-1.39, P = 6.2×10(-5)). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18-2.10, P = 1.9×10(-3)). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.

Genome-wide association study yields variants at 20p12.2 that associate with urinary bladder cancer.

Genome-wide association studies (GWAS) of urinary bladder cancer (UBC) have yielded common variants at 12 loci that associate with risk of the disease. We report here the results of a GWAS of UBC including 1670 UBC cases and 90 180 controls, followed by replication analysis in additional 5266 UBC cases and 10 456 controls. We tested a dataset containing 34.2 million variants, generated by imputation based on whole-genome sequencing of 2230 Icelanders. Several correlated variants at 20p12, represented by rs62185668, show genome-wide significant association with UBC after combining discovery and replication results (OR = 1.19, P = 1.5 × 10(-11) for rs62185668-A, minor allele frequency = 23.6%). The variants are located in a non-coding region approximately 300 kb upstream from the JAG1 gene, an important component of the Notch signaling pathways that may be oncogenic or tumor suppressive in several forms of cancer. Our results add to the growing number of UBC risk variants discovered through GWAS.

Known susceptibility SNPs for sporadic prostate cancer show a similar association with "hereditary" prostate cancer.

BACKGROUND: More than 70 single nucleotide polymorphisms (SNPs) have been reported to be associated with prostate cancer (PC) risk; these were mainly identified in the general population with predominantly sporadic PC (SPC). Previous studies have suggested similar associations between a selection of these SNPs and hereditary PC (HPC). Our aim was to evaluate the effect of all known PC risk SNPs and their discriminative value for SPC and HPC. METHODS: Seventy-four PC susceptibility SNPs (reported in literature up to June 2014) were genotyped in a population-based series of 620 SPC patients, 312 HPC patients from the national Dutch registry and 1819 population-based referents. Association analyses were performed using logistic regression, focusing on directional consistency of the odds ratios (ORs) with those in the original reports, that is, whether the OR was in the same direction as in the original report. Discriminative performance was evaluated by a genetic risk score used in logistic regression and receiver operating characteristic (ROC) curve analyses. RESULTS: Directional consistency was seen for 62 SNPs in SPC and 64 SNPs in HPC, 56 of which overlapped. ORs were mostly higher for HPC with 22 ORs >1.25 versus 5 for SPC. Discriminative performance was better for HPC with an area under the ROC curve of 0.73 versus 0.64 for SPC. CONCLUSIONS: A large overlap was found for the associations between low-penetrance susceptibility SNPs and SPC and HPC, suggesting a similarity in genetic etiology. This warrants a reconsideration of "HPC" and a restrictive policy toward prostate-specific antigen testing in men with a positive family history. Genetic risk scores might be used for PC risk stratification on the population level.

Iron and hepcidin as risk factors in atherosclerosis: what do the genes say?

BACKGROUND: Previous reports suggested a role for iron and hepcidin in atherosclerosis. Here, we evaluated the causality of these associations from a genetic perspective via (i) a Mendelian randomization (MR) approach, (ii) study of association of atherosclerosis-related single nucleotide polymorphisms (SNPs) with iron and hepcidin, and (iii) estimation of genomic correlations between hepcidin, iron and atherosclerosis. RESULTS: Analyses were performed in a general population sample. Iron parameters (serum iron, serum ferritin, total iron-binding capacity and transferrin saturation), serum hepcidin and genome-wide SNP data were available for N = 1,819; non-invasive measurements of atherosclerosis (NIMA), i.e., presence of plaque, intima media thickness and ankle-brachial index (ABI), for N = 549. For the MR, we used 12 iron-related SNPs that were previously identified in a genome-wide association meta-analysis on iron status, and assessed associations of individual SNPs and quartiles of a multi-SNP score with NIMA. Quartile 4 versus quartile 1 of the multi-SNP score showed directionally consistent associations with the hypothesized direction of effect for all NIMA in women, indicating that increased body iron status is a risk factor for atherosclerosis in women. We observed no single SNP associations that fit the hypothesized directions of effect between iron and NIMA, except for rs651007, associated with decreased ferritin concentration and decreased atherosclerosis risk. Two of six NIMA-related SNPs showed association with the ratio hepcidin/ferritin, suggesting that an increased hepcidin/ferritin ratio increases atherosclerosis risk. Genomic correlations were close to zero, except for hepcidin and ferritin with ABI at rest [-0.27 (SE 0.34) and -0.22 (SE 0.35), respectively] and ABI after exercise [-0.29 (SE 0.34) and -0.30 (0.35), respectively]. The negative sign indicates an increased atherosclerosis risk with increased hepcidin and ferritin concentrations. CONCLUSIONS: Our results suggest a potential causal role for hepcidin and ferritin in atherosclerosis, and may indicate that iron status is causally related to atherosclerosis in women.

Serum hepcidin is associated with presence of plaque in postmenopausal women of a general population.

OBJECTIVE: Iron and the iron regulatory hormone hepcidin, major determinant of body iron distribution, are hypothesized to play a role in cardiovascular disease. Here, we assess the associations of hepcidin as well as ferritin, iron, total iron-binding capacity, and transferrin saturation (ie, iron parameters) with noninvasive measurements of atherosclerosis in men and women of a population-based cohort. APPROACH AND RESULTS: We included 766 participants of the Nijmegen Biomedical Study aged 46 to 67 years for whom serum measurements of hepcidin, iron parameters, and noninvasive measurements of atherosclerosis were available. Noninvasive measurements of atherosclerosis were presence of plaque, ankle-brachial index, and intima-media thickness. We performed multivariable logistic and linear regression analyses using quartiles of hepcidin and iron parameters. Analyses were stratified by sex and adjusted for several demographic, clinical, and biochemical determinants, including traditional risk factors of cardiovascular disease based on the Framingham risk score. Hepcidin and the hepcidin/ferritin ratio, reflecting hepcidin expression relative to iron stores, were significantly associated with the presence of plaque in women (adjusted odds ratios for quartile 4 versus quartile 1 [95% confidence intervals] of 3.07 [1.36-6.90] and 2.31 [1.03-5.18], respectively). The hepcidin/ferritin ratio was significantly and negatively associated with ankle-brachial index at rest in men and women (adjusted ß for quartile 4 versus quartile 1 [95% confidence intervals] of -0.03 [-0.07 to 0.00] and -0.04 [-0.06 to -0.01], respectively). CONCLUSIONS: Our results suggest that the body iron distribution as determined by hepcidin affects the development of atherosclerosis in women.

Associations of common variants in HFE and TMPRSS6 with iron parameters are independent of serum hepcidin in a general population: a replication study.

BACKGROUND: Genome-wide association studies have convincingly shown that single nucleotide polymorphisms (SNPs) in HFE and TMPRSS6 are associated with iron parameters. It was commonly thought that these associations could be explained by the intermediate effect on hepcidin concentration. A recent study in an isolated Italian population, however, concluded that these associations were not exclusively dependent on hepcidin values. We report here the second study to investigate the role of hepcidin in the associations between common variants in HFE and TMPRSS6 with iron parameters. METHODS: We extracted 101 SNPs in HFE and TMPRSS6 from genome-wide imputed SNP data of 1832 individuals from the general population (Nijmegen Biomedical Study). Single locus and haplotype associations with serum iron parameters and hepcidin were studied using linear regression analyses. RESULTS: We found that HFE rs1800562 and TMPRSS6 rs855791 are the main determinants of HFE and TMPRSS6 related variation in serum iron, ferritin, transferrin saturation, and total iron binding capacity. These SNPs are associated with the ratios hepcidin/ferritin (p<1×10(-5)) and hepcidin/transferrin saturation (p<1×10(-3)), but not with serum hepcidin (p>0.2). Adjustment for hepcidin or the ratio hepcidin/ferritin did not decrease the strength of the SNP-iron parameter associations. CONCLUSIONS: Our results do not support an intermediate role for hepcidin in the SNP-iron parameter associations, which confirms previous findings, and indicate a pleiotropic SNP effect on the hepcidin ratios and the iron parameters. Taken together, this suggests that there might be other, yet unknown, serum hepcidin independent mechanisms which play a role in the association of HFE and TMPRSS6 variants with serum iron parameters.

Towards personalized TSH reference ranges: A genetic and population-based approach in three independent cohorts.

BACKGROUND: Serum thyroid-stimulating hormone (TSH) measurement is the diagnostic cornerstone for primary thyroid dysfunction. There is high inter-individual, but limited intra-individual variation in TSH concentrations, largely due to genetic factors. The currently used wide population-based reference intervals may lead to inappropriate management decisions. METHODS: A polygenic score (PGS) including 59 genetic variants was used to calculate genetically-determined TSH reference ranges in a thyroid disease-free cohort (N=6,834). Its effect on reclassification of diagnoses was investigated when compared to using population-based reference ranges. Next, results were validated in a second independent population-based thyroid disease-free cohort (N=3,800). Potential clinical implications were assessed in a third independent population-based cohort including individuals without thyroid disease (N=26,321) as well as individuals on levothyroxine (LT4) treatment (N=1,132). RESULTS: PGS was a much stronger predictor of individual TSH concentrations than FT4 (total variance in TSH concentrations explained 9.2-11.1% vs. 2.4-2.7%, respectively) or any other non-genetic factor (total variance in TSH concentrations explained 0.2-1.8%). Genetically-determined TSH reference ranges differed significantly between PGS quartiles in all cohorts, while the differences in FT4 concentrations were absent or only minor. Up to 24.7-30.1% of individuals, previously classified as having subclinical hypo- and hyperthyroidism when using population-based TSH reference ranges, were reclassified as euthyroid when genetically-determined TSH reference ranges were applied. Individuals in the higher PGS quartiles had a higher probability of being prescribed LT4 treatment compared to individuals from the lower PGS quartiles (3.3% in Q1 vs. 5.2% in Q4, Pfor trend =1.7x10-8). CONCLUSIONS: Individual genetic profiles have potential to personalize TSH reference ranges, with large effects on reclassification of diagnosis and LT4 prescriptions. As the currently used PGS can only predict approximately 10% of inter-individual variation in TSH concentrations, it should be further improved when more genetic variants determining TSH concentrations are identified in future studies.

Causal relationships between risk of venous thromboembolism and 18 cancers: a bidirectional Mendelian randomization analysis.

BACKGROUND: People with cancer experience high rates of venous thromboembolism (VTE). Risk of subsequent cancer is also increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer. METHODS: We used data from large genome-wide association study meta-analyses to perform bidirectional Mendelian randomization analyses to estimate causal associations between genetic liability to VTE and risk of 18 different cancers. RESULTS: We found no conclusive evidence that genetic liability to VTE was causally associated with an increased incidence of cancer, or vice versa. We observed an association between liability to VTE and pancreatic cancer risk [odds ratio for pancreatic cancer: 1.23 (95% confidence interval: 1.08-1.40) per log-odds increase in VTE risk, P = 0.002]. However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence to suggest a causal relationship. CONCLUSIONS: These findings do not support the hypothesis that genetic liability to VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesize evidence for these mechanisms.