Jouke-Jan Hottenga

Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts

Recent genome-wide association (GWA) studies of lipids have been conducted in samples ascertained for other phenotypes, particularly diabetes. Here we report the first GWA analysis of loci affecting total cholesterol (TC), low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglycerides sampled randomly from 16 population-based cohorts and genotyped using mainly the Illumina HumanHap300-Duo platform. Our study included a total of 17,797–22,562 persons, aged 18–104 years and from geographic regions spanning from the Nordic countries to Southern Europe. We established 22 loci associated with serum lipid levels at a genome-wide significance level (P < 5 × 10−8), including 16 loci that were identified by previous GWA studies. The six newly identified loci in our cohort samples are ABCG5 (TC, P = 1.5 × 10−11; LDL, P = 2.6 × 10−10), TMEM57 (TC, P = 5.4 × 10−10), CTCF-PRMT8 region (HDL, P = 8.3 × 10−16), DNAH11 (LDL, P = 6.1 × 10−9), FADS3-FADS2 (TC, P = 1.5 × 10−10; LDL, P = 4.4 × 10−13) and MADD-FOLH1 region (HDL, P = 6 × 10−11). For three loci, effect sizes differed significantly by sex. Genetic risk scores based on lipid loci explain up to 4.8% of variation in lipids and were also associated with increased intima media thickness (P = 0.001) and coronary heart disease incidence (P = 0.04). The genetic risk score improves the screening of high-risk groups of dyslipidemia over classical risk factors.

Identification of seven loci affecting mean telomere length and their association with disease

Interindividual variation in mean leukocyte telomere length (LTL) is associated with cancer and several age-associated diseases. We report here a genome-wide meta-analysis of 37,684 individuals with replication of selected variants in an additional 10,739 individuals. We identified seven loci, including five new loci, associated with mean LTL (P < 5 × 10−8). Five of the loci contain candidate genes (TERC, TERT, NAF1, OBFC1 and RTEL1) that are known to be involved in telomere biology. Lead SNPs at two loci (TERC and TERT) associate with several cancers and other diseases, including idiopathic pulmonary fibrosis. Moreover, a genetic risk score analysis combining lead variants at all 7 loci in 22,233 coronary artery disease cases and 64,762 controls showed an association of the alleles associated with shorter LTL with increased risk of coronary artery disease (21% (95% confidence interval, 5–35%) per standard deviation in LTL, P = 0.014). Our findings support a causal role of telomere-length variation in some age-related diseases.

Netherlands Twin Register: From Twins to Twin Families

In the late 1980s The Netherlands Twin Register (NTR) was established by recruiting young twins and multiples at birth and by approaching adolescent and young adult twins through city councils. The Adult NTR (ANTR) includes twins, their parents, siblings, spouses and their adult offspring. The number of participants in the ANTR who take part in survey and / or laboratory studies is over 22,000 subjects. A special group of participants consists of sisters who are mothers of twins. In the Young NTR (YNTR), data on more than 50,000 young twins have been collected. Currently we are extending the YNTR by including siblings of twins. Participants in YNTR and ANTR have been phenotyped every 2 to 3 years in longitudinal survey studies, since 1986 and 1991 for the YNTR and ANTR, respectively. The resulting large population-based datasets are used for genetic epidemiological studies and also, for example, to advance phenotyping through the development of new syndrome scales based on existing items from other inventories. New research developments further include brain imaging studies in selected and unselected groups, clinical assessment of psychopathology through interviews, and cross-referencing the NTR database to other national databases. A large biobank enterprise is ongoing in the ANTR in which blood and urine samples are collected for genotyping, expression analysis, and metabolomics studies. In this paper we give an update on the YNTR and ANTR phenotyping and on the ongoing ANTR biobank studies.

Variation in plasminogen-activator-inhibitor-1 gene and risk of meningococcal septic shock

BACKGROUND Some patients infected with Neisseria meningitidis develop septic shock accompanied by fibrin deposition and microthrombus formation in various organs, whereas others develop bacteraemia or meningitis without septic shock. We investigated whether genetic differences in the fibrinolytic system influence the development of meningococcal septic shock. METHODS We investigated 50 patients who survived meningococcal infection, and 131 controls from the same geographical region. Because we had no information on genotypes of patients who died, we also genotyped 183 first-degree relatives of a consecutive series of patients with meningococcal infection for the 4G/5G deletion/insertion polymorphism in the promoter region of the plasminogen-activator-inhibitor-1 gene (PAI-1). The 4G allele is associated with increased gene transcription in cell lines in vitro and with increased PAI-1 concentrations in carriers in vivo. FINDINGS The allele frequencies of 4G and 5G were similar between patients and controls. However, the 4G/4G genotype was present in only 9% of relatives of patients with meningitis compared with 36% of relatives of patients with septic shock. The 5G/5G genotype was more common among relatives of patients with meningitis (31 vs 11%, p=0.001). Patients whose relatives were carriers of the 4G/4G genotype had a six-fold higher risk of developing septic shock than meningitis (odds ratio 5.9 [95% CI 1.9-18.0]) compared with all other genotypes. INTERPRETATION Variation in the PAI-1 gene does not affect the probability of contracting meningococcal infection, but does influence the development of septic shock.

Genome-wide meta-analysis identifies new susceptibility loci for migraine

Migraine is the most common brain disorder, affecting approximately 14% of the adult population, but its molecular mechanisms are poorly understood. We report the results of a meta-analysis across 29 genome-wide association studies, including a total of 23,285 individuals with migraine (cases) and 95,425 population-matched controls. We identified 12 loci associated with migraine susceptibility (P < 5 × 10−8). Five loci are new: near AJAP1 at 1p36, near TSPAN2 at 1p13, within FHL5 at 6q16, within C7orf10 at 7p14 and near MMP16 at 8q21. Three of these loci were identified in disease subgroup analyses. Brain tissue expression quantitative trait locus analysis suggests potential functional candidate genes at four loci: APOA1BP, TBC1D7, FUT9, STAT6 and ATP5B.

Common variants in TMPRSS6 are associated with iron status and erythrocyte volume

We report a genome-wide association study to iron status. We identify an association of SNPs in TPMRSS6 to serum iron (rs855791, combined P = 1.5 × 10−20), transferrin saturation (combined P = 2.2 × 10−23) and erythrocyte mean cell volume (MCV, combined P = 1.1 × 10−10). We also find suggestive evidence of association with blood hemoglobin levels (combined P = 5.3 × 10−7). These findings demonstrate the involvement of TMPRSS6 in control of iron homeostasis and in normal erythropoiesis.

Heritability and genomics of gene expression in peripheral blood

We assessed gene expression profiles in 2,752 twins, using a classic twin design to quantify expression heritability and quantitative trait loci (eQTLs) in peripheral blood. The most highly heritable genes (∼777) were grouped into distinct expression clusters, enriched in gene-poor regions, associated with specific gene function or ontology classes, and strongly associated with disease designation. The design enabled a comparison of twin-based heritability to estimates based on dizygotic identity-by-descent sharing and distant genetic relatedness. Consideration of sampling variation suggests that previous heritability estimates have been upwardly biased. Genotyping of 2,494 twins enabled powerful identification of eQTLs, which we further examined in a replication set of 1,895 unrelated subjects. A large number of non-redundant local eQTLs (6,756) met replication criteria, whereas a relatively small number of distant eQTLs (165) met quality control and replication standards. Our results provide a new resource toward understanding the genetic control of transcription.

Combined Genome Scans for Body Stature in 6,602 European Twins: Evidence for Common Caucasian Loci

Twin cohorts provide a unique advantage for investigations of the role of genetics and environment in the etiology of variation in common complex traits by reducing the variance due to environment, age, and cohort differences. The GenomEUtwin (http://www.genomeutwin.org) consortium consists of eight twin cohorts (Australian, Danish, Dutch, Finnish, Italian, Norwegian, Swedish, and United Kingdom) with the total resource of hundreds of thousands of twin pairs. We performed quantitative trait locus (QTL) analysis of one of the most heritable human complex traits, adult stature (body height) using genome-wide scans performed for 3,817 families (8,450 individuals) derived from twin cohorts from Australia, Denmark, Finland, Netherlands, Sweden, and United Kingdom with an approximate ten-centimorgan microsatellite marker map. The marker maps for different studies differed and they were combined and related to the sequence positions using software developed by us, which is publicly available (https://apps.bioinfo.helsinki.fi/software/cartographer.aspx). Variance component linkage analysis was performed with age, sex, and country of origin as covariates. The covariate adjusted heritability was 81% for stature in the pooled dataset. We found evidence for a major QTL for human stature on 8q21.3 (multipoint logarithm of the odds 3.28), and suggestive evidence for loci on Chromosomes X, 7, and 20. Some evidence of sex heterogeneity was found, however, no obvious female-specific QTLs emerged. Several cohorts contributed to the identified loci, suggesting an evolutionarily old genetic variant having effects on stature in European-based populations. To facilitate the genetic studies of stature we have also set up a website that lists all stature genome scans published and their most significant loci (http://www.genomeutwin.org/stature_gene_map.htm).

Genome-wide Association Study of Smoking Initiation and Current Smoking

For the identification of genes associated with smoking initiation and current smoking, genome-wide association analyses were carried out in 3497 subjects. Significant genes that replicated in three independent samples (n = 405, 5810, and 1648) were visualized into a biologically meaningful network showing cellular location and direct interaction of their proteins. Several interesting groups of proteins stood out, including glutamate receptors (e.g., GRIN2B, GRIN2A, GRIK2, GRM8), proteins involved in tyrosine kinase receptor signaling (e.g., NTRK2, GRB14), transporters (e.g., SLC1A2, SLC9A9) and cell-adhesion molecules (e.g., CDH23). We conclude that a network-based genome-wide association approach can identify genes influencing smoking behavior.

Genome-wide association meta-analysis of human longevity identifies a novel locus conferring survival beyond 90 years of age

The genetic contribution to the variation in human lifespan is ∼25%. Despite the large number of identified disease-susceptibility loci, it is not known which loci influence population mortality. We performed a genome-wide association meta-analysis of 7729 long-lived individuals of European descent (≥85 years) and 16 121 younger controls (<65 years) followed by replication in an additional set of 13 060 long-lived individuals and 61 156 controls. In addition, we performed a subset analysis in cases aged ≥90 years. We observed genome-wide significant association with longevity, as reflected by survival to ages beyond 90 years, at a novel locus, rs2149954, on chromosome 5q33.3 (OR = 1.10, P = 1.74 × 10−8). We also confirmed association of rs4420638 on chromosome 19q13.32 (OR = 0.72, P = 3.40 × 10−36), representing the TOMM40/APOE/APOC1 locus. In a prospective meta-analysis (n = 34 103), the minor allele of rs2149954 (T) on chromosome 5q33.3 associates with increased survival (HR = 0.95, P = 0.003). This allele has previously been reported to associate with low blood pressure in middle age. Interestingly, the minor allele (T) associates with decreased cardiovascular mortality risk, independent of blood pressure. We report on the first GWAS-identified longevity locus on chromosome 5q33.3 influencing survival in the general European population. The minor allele of this locus associates with low blood pressure in middle age, although the contribution of this allele to survival may be less dependent on blood pressure. Hence, the pleiotropic mechanisms by which this intragenic variation contributes to lifespan regulation have to be elucidated.