26
|
Wijsman CA, van Opstal AM, Kan HE, Maier AB, Westendorp RGJ, Slagboom PE, Webb AG, Mooijaart SP, van Heemst D. Proton magnetic resonance spectroscopy shows lower intramyocellular lipid accumulation in middle-aged subjects predisposed to familial longevity. Am J Physiol Endocrinol Metab 2012; 302:E344-8. [PMID: 22094471 DOI: 10.1152/ajpendo.00455.2011] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Families predisposed to longevity show enhanced glucose tolerance and skeletal muscle insulin sensitivity compared with controls, independent of body composition and physical activity. Intramyocellular lipid (IMCL) accumulation in skeletal muscle has been associated with insulin resistance. Here, we assessed whether subjects enriched for familial longevity have lower IMCL levels. We determined IMCL levels in 48 subjects from the Leiden Longevity Study, comprising 24 offspring of nonagenarian siblings and 24 partners thereof as control subjects. IMCL levels were assessed noninvasively using short echo time proton magnetic resonance spectroscopy ((1)H-MRS) of the tibialis anterior muscle with a 7 Tesla human MR scanner. IMCL levels were calculated relative to the total creatine (tCr) CH3 signal. Physical activity was assessed using the International Physical Activity Questionnaire (IPAQ). After correction for age, sex, BMI, and physical activity, offspring of long-lived nonagenarian siblings tended to show lower IMCL levels compared with controls (IMCL/tCr: 3.1 ± 0.5 vs. 4.5 ± 0.5, respectively, P = 0.051). In a pairwise comparison, this difference reached statistical significance (P = 0.038). We conclude that offspring of nonagenarian siblings predisposed to longevity show lower IMCL levels compared with environmentally matched control subjects. Future research should focus on assessing what mechanisms may explain the lower IMCL levels in familial longevity.
Collapse
|
27
|
Talens RP, Jukema JW, Trompet S, Kremer D, Westendorp RGJ, Lumey LH, Sattar N, Putter H, Slagboom PE, Heijmans BT. Hypermethylation at loci sensitive to the prenatal environment is associated with increased incidence of myocardial infarction. Int J Epidemiol 2011; 41:106-15. [PMID: 22101166 DOI: 10.1093/ije/dyr153] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
BACKGROUND Human epidemiological studies suggest that small size at birth and food deprivation during gestation confer an excess risk of coronary heart diseases (CHD) in adulthood, frequently in a sex-specific manner. Prior epigenetic studies indicate that such prenatal conditions are marked by persistent and sometimes sex-specific changes in DNA methylation. Here, we have investigated the association between DNA methylation and myocardial infarction (MI) at six loci sensitive to prenatal nutrition, anticipating potential sex-specificity. Method Within the placebo group of the PROSPER trial on pravastatin and the risk of CHD, we compared all individuals who were event free at baseline and developed MI during 3 years' follow-up (n = 122) with a similar-sized control group. Methylation at IL10, LEP, ABCA1, IGF2, INS and GNASAS was measured in DNA extracted from leucocytes using mass spectrometry. RESULTS DNA methylation at GNASAS was modestly higher in MI cases compared with controls (P = 0.030). A significant sex interaction was observed for INS (P = 0.014) and GNASAS (P = 0.031). Higher DNA methylation at these loci was associated with MI among women (INS: +2.5%, P = 0.002; GNASAS: +4.2%, P = 0.001). Hypermethylation at one locus and at both loci was associated with odds ratios (ORs) of 2.8 and 8.6, respectively (P(trend) = 3.0 × 10(-4)). No association was observed among men. CONCLUSIONS The risk of MI in women is associated with DNA methylation marks at specific loci previously shown to be sensitive to prenatal conditions. This observation may reflect a developmental component of MI.
Collapse
|
28
|
Martella F, Vermunt JK, Beekman M, Westendorp RGJ, Slagboom PE, Houwing-Duistermaat JJ. A mixture model with random-effects components for classifying sibling pairs. Stat Med 2011; 30:3252-64. [PMID: 21905068 DOI: 10.1002/sim.4365] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2010] [Accepted: 07/21/2011] [Indexed: 11/09/2022]
Abstract
In healthy aging research, typically multiple health outcomes are measured, representing health status. The aim of this paper was to develop a model-based clustering approach to identify homogeneous sibling pairs according to their health status. Model-based clustering approaches will be considered on the basis of linear mixed effect model for the mixture components. Class memberships of siblings within pairs are allowed to be correlated, and within a class the correlation between siblings is modeled using random sibling pair effects. We propose an expectation-maximization algorithm for maximum likelihood estimation. Model performance is evaluated via simulations in terms of estimating the correct parameters, degree of agreement, and the ability to detect the correct number of clusters. The performance of our model is compared with the performance of standard model-based clustering approaches. The methods are used to classify sibling pairs from the Leiden Longevity Study according to their health status. Our results suggest that homogeneous healthy sibling pairs are associated with a longer life span. Software is available for fitting the new models.
Collapse
|
29
|
Yusuf E, Bijsterbosch J, Slagboom PE, Rosendaal FR, Huizinga TWJ, Kloppenburg M. Body mass index and alignment and their interaction as risk factors for progression of knees with radiographic signs of osteoarthritis. Osteoarthritis Cartilage 2011; 19:1117-22. [PMID: 21722745 DOI: 10.1016/j.joca.2011.06.001] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/26/2010] [Revised: 05/18/2011] [Accepted: 06/05/2011] [Indexed: 02/02/2023]
Abstract
OBJECTIVE To investigate in which way body mass index (BMI) and alignment affect the risk for knee osteoarthritis (OA) progression. METHODS Radiographs of 181 knees from 155 patients (85% female, mean age 60 years) with radiographic signs of OA were analyzed at baseline and after 6 years. Progression was defined as 1-point increase in joint space narrowing score in the medial or lateral tibiofemoral (TF) compartment or having knee prosthesis during the follow-up for knees with a Kellgren and Lawrence score ≥ 1 at baseline. BMI at baseline was classified as normal (<25 kg/m(2)), overweight (25-30) and obese (>30). Knee alignment on baseline radiographs was categorized as normal (TF angle between 182° and 184°) and malalignment (<182° or >184°). We estimated the risk ratio (RR) with 95% confidence interval for knee OA progression for overweight and obese patients and for malaligned knees relative to normal using generalized estimating equations (GEE). Additionally, we estimated the added effect when BMI and malalignment were present together on progression of knee OA. Adjustments were made for age and sex. RESULTS Seventy-six knees (42%) showed progression: 27 in lateral and 66 in medial compartment. Knees from overweight and obese patients had an increased risk for progression (RR 2.4 (1.-3.6) and 2.9 (1.7-4.1), respectively). RRs of progression for malaligned, varus and valgus knee were 2.0 (1.3-2.8), 2.3 (1.4-3.1), and 1.7 (0.97-2.6), respectively. When BMI and malalignment were included in one model, the effect of overweight, obesity and malalignment did not change. The added effect when overweight and malalignment were present was 17%. CONCLUSION Overweight is associated with progression of knee OA and shows a small interaction with alignment. Losing weight might be helpful in preventing the progression of knee OA.
Collapse
|
30
|
Waarsing JH, Kloppenburg M, Slagboom PE, Kroon HM, Houwing-Duistermaat JJ, Weinans H, Meulenbelt I. Osteoarthritis susceptibility genes influence the association between hip morphology and osteoarthritis. ARTHRITIS AND RHEUMATISM 2011; 63:1349-54. [PMID: 21400473 DOI: 10.1002/art.30288] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
OBJECTIVE The identified osteoarthritis (OA) susceptibility genes are mainly active in skeletal development and could thus affect joint geometry. Because nonoptimal joint geometry is a risk factor for the development of OA, we investigated if and how the path that leads from nonoptimal joint geometry to OA of the hip is influenced by these genes. METHODS The shape of the hips of subjects in the Genetics, Osteoarthritis and Progression Study, consisting of sibling pairs with symptomatic OA at multiple joint locations, was quantified by applying a statistical shape model to radiographs. Shape aspects (modes) were correlated to OA characteristics. We then tested for the association of shape modes with OA susceptibility single-nucleotide polymorphisms (SNPs) of GDF5, FRZB, and DIO2. RESULTS Four of 23 shape modes (mode 1, mode 17, mode 18, and mode 21) were strongly associated with OA characteristics. We observed a significant interaction between carrier status of DIO2 rs12885300 and hip OA characteristics for mode 1 (P = 0.005). This indicates that this specific aspect of hip shape correlates with OA characteristics only in carriers of the susceptibility allele. CONCLUSION Our results suggest that it is more likely that the rs12885300 SNP of DIO2 increases the vulnerability of cartilage to nonoptimal bone shapes rather than directly influencing the formation of these shapes.
Collapse
|
31
|
Slagboom PE, Beekman M, Passtoors WM, Deelen J, Vaarhorst AAM, Boer JM, van den Akker EB, van Heemst D, de Craen AJM, Maier AB, Rozing M, Mooijaart SP, Heijmans BT, Westendorp RGJ. Genomics of human longevity. Philos Trans R Soc Lond B Biol Sci 2011; 366:35-42. [PMID: 21115528 PMCID: PMC3001312 DOI: 10.1098/rstb.2010.0284] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
In animal models, single-gene mutations in genes involved in insulin/IGF and target of rapamycin signalling pathways extend lifespan to a considerable extent. The genetic, genomic and epigenetic influences on human longevity are expected to be much more complex. Strikingly however, beneficial metabolic and cellular features of long-lived families resemble those in animals for whom the lifespan is extended by applying genetic manipulation and, especially, dietary restriction. Candidate gene studies in humans support the notion that human orthologues from longevity genes identified in lower species do contribute to longevity but that the influence of the genetic variants involved is small. Here we discuss how an integration of novel study designs, labour-intensive biobanking, deep phenotyping and genomic research may provide insights into the mechanisms that drive human longevity and healthy ageing, beyond the associations usually provided by molecular and genetic epidemiology. Although prospective studies of humans from the cradle to the grave have never been performed, it is feasible to extract life histories from different cohorts jointly covering the molecular changes that occur with age from early development all the way up to the age at death. By the integration of research in different study cohorts, and with research in animal models, biological research into human longevity is thus making considerable progress.
Collapse
|
32
|
Panoutsopoulou K, Southam L, Elliott KS, Wrayner N, Zhai G, Beazley C, Thorleifsson G, Arden NK, Carr A, Chapman K, Deloukas P, Doherty M, McCaskie A, Ollier WER, Ralston SH, Spector TD, Valdes AM, Wallis GA, Wilkinson JM, Arden E, Battley K, Blackburn H, Blanco FJ, Bumpstead S, Cupples LA, Day-Williams AG, Dixon K, Doherty SA, Esko T, Evangelou E, Felson D, Gomez-Reino JJ, Gonzalez A, Gordon A, Gwilliam R, Halldorsson BV, Hauksson VB, Hofman A, Hunt SE, Ioannidis JPA, Ingvarsson T, Jonsdottir I, Jonsson H, Keen R, Kerkhof HJM, Kloppenburg MG, Koller N, Lakenberg N, Lane NE, Lee AT, Metspalu A, Meulenbelt I, Nevitt MC, O'Neill F, Parimi N, Potter SC, Rego-Perez I, Riancho JA, Sherburn K, Slagboom PE, Stefansson K, Styrkarsdottir U, Sumillera M, Swift D, Thorsteinsdottir U, Tsezou A, Uitterlinden AG, van Meurs JBJ, Watkins B, Wheeler M, Mitchell S, Zhu Y, Zmuda JM, Zeggini E, Loughlin J. Insights into the genetic architecture of osteoarthritis from stage 1 of the arcOGEN study. Ann Rheum Dis 2010; 70:864-7. [PMID: 21177295 PMCID: PMC3070286 DOI: 10.1136/ard.2010.141473] [Citation(s) in RCA: 82] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Objectives The genetic aetiology of osteoarthritis has not yet been elucidated. To enable a well-powered genome-wide association study (GWAS) for osteoarthritis, the authors have formed the arcOGEN Consortium, a UK-wide collaborative effort aiming to scan genome-wide over 7500 osteoarthritis cases in a two-stage genome-wide association scan. Here the authors report the findings of the stage 1 interim analysis. Methods The authors have performed a genome-wide association scan for knee and hip osteoarthritis in 3177 cases and 4894 population-based controls from the UK. Replication of promising signals was carried out in silico in five further scans (44 449 individuals), and de novo in 14 534 independent samples, all of European descent. Results None of the association signals the authors identified reach genome-wide levels of statistical significance, therefore stressing the need for corroboration in sample sets of a larger size. Application of analytical approaches to examine the allelic architecture of disease to the stage 1 genome-wide association scan data suggests that osteoarthritis is a highly polygenic disease with multiple risk variants conferring small effects. Conclusions Identifying loci conferring susceptibility to osteoarthritis will require large-scale sample sizes and well-defined phenotypes to minimise heterogeneity.
Collapse
|
33
|
Rozing MP, Houwing-Duistermaat JJ, Slagboom PE, Beekman M, Frölich M, de Craen AJM, Westendorp RGJ, van Heemst D. Familial longevity is associated with decreased thyroid function. J Clin Endocrinol Metab 2010; 95:4979-84. [PMID: 20739380 DOI: 10.1210/jc.2010-0875] [Citation(s) in RCA: 93] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
CONTEXT A relation between low thyroid activity and prolonged life span in humans has been observed. Several studies have demonstrated hereditary and genetic influences on thyroid function. OBJECTIVE The objective of the study was to test whether low thyroid activity associated with extreme longevity constitutes a heritable phenotype, which could contribute to the familial longevity observed in the Leiden Longevity Study. DESIGN This was a cross-sectional study. SETTING The study was conducted at a university hospital in the city of Leiden, The Netherlands. PARTICIPANTS Eight hundred fifty-nine nonagenarian siblings (median age 92.9 yr) from 421 long-lived families participated in the study. Families were recruited from the entire Dutch population if at least two long-lived siblings were alive and fulfilled the age criterion of age of 89 yr or older for males and 91 yr or older for females. There were no selection criteria on health or demographic characteristics. INTERVENTION Blood samples were taken for determination of serum parameters of thyroid function. MAIN OUTCOME MEASURE We calculated the family mortality history score of the parents of the nonagenarian siblings and related this to thyroid function parameters in the nonagenarian siblings. RESULTS We found that a lower family mortality history score (less mortality) of the parents of nonagenarian siblings was associated with higher serum TSH levels (P = 0.005) and lower free T(4) levels (P = 0.002) as well as lower free T(3) levels (P = 0.034) in the nonagenarian siblings. CONCLUSIONS Our findings support the previous observation that low thyroid activity in humans constitutes a heritable phenotype that contributes to exceptional familial longevity observed in the Leiden Longevity Study.
Collapse
|
34
|
Callegaro A, Meulenbelt I, Kloppenburg M, Slagboom PE, Houwing-Duistermaat JJ. Allele-sharing statistics using information on family history. Ann Hum Genet 2010; 74:547-54. [PMID: 20707788 DOI: 10.1111/j.1469-1809.2010.00602.x] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
When conducting genetic studies for complex traits, large samples are commonly required to detect new genetic factors. A possible strategy to decrease the sample size is to reduce heterogeneity using available information. In this paper we propose a new class of model-free linkage analysis statistics which takes into account the information given by the ungenotyped affected relatives (positive family history). This information is included into the scoring function of classical allele-sharing statistics. We studied pedigrees of affected sibling pairs with one ungenotyped affected relative. We show that, for rare allele common complex diseases, the proposed method increases the expected power to detect linkage. Allele-sharing methods were applied to the symptomatic osteoarthritis GARP study where taking into account the family-history increased considerably the evidence of linkage in the region of the DIO2 susceptibility locus.
Collapse
|
35
|
Middeldorp CM, Slof-Op 't Landt MCT, Medland SE, van Beijsterveldt CEM, Bartels M, Willemsen G, Hottenga JJ, de Geus EJC, Suchiman HED, Dolan CV, Neale MC, Slagboom PE, Boomsma DI. Anxiety and depression in children and adults: influence of serotonergic and neurotrophic genes? GENES BRAIN AND BEHAVIOR 2010; 9:808-16. [PMID: 20633049 PMCID: PMC3151552 DOI: 10.1111/j.1601-183x.2010.00619.x] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
There are two major hypotheses regarding the etiology of anxiety and depression: the mono-amine hypothesis and the hypothesis of an abnormal stress response acting partly via reduced neurogenesis. Association studies have focused on genes involved in these processes, but with inconclusive results. This study investigated the effect of 45 single nucleotide polymorphisms (SNPs) in genes encoding for serotonin receptors 1A, 1D, 2A, catechol-O-methyltransferase (COMT), tryptophane hydroxylase type 2 (TPH2), brain derived neurotrophic factor (BDNF), PlexinA2 and regulators of G-protein-coupled signaling (RGS) 2, 4, 16. Anxious depression (A/D) symptoms were assessed five times in 11 years in over 11 000 adults with 1504 subjects genotyped and at age 7, 10, 12 and during adolescence in over 20 000 twins with 1078 subjects genotyped. In both cohorts, a longitudinal model with one latent factor loading on all A/D measures over time was analysed. The genetic association effect modeled at the level of this latent factor was 60% and 70% heritable in the children and adults, respectively, and explained around 50% of the total phenotypic variance. Power analyses showed that the samples contained 80% power to detect an effect explaining between 1.4% and 3.6% of the variance. However, no SNP showed a consistent effect on A/D. To conclude, this longitudinal study in children and adults found no association of SNPs in the serotonergic system or core regulators of neurogenesis with A/D. Overall, there has been no convincing evidence, so far, for a role of genetic variation in these pathways in the development of anxiety and depression.
Collapse
|
36
|
Slof-Op 't Landt MCT, Dolan CV, Rebollo-Mesa I, Bartels M, van Furth EF, van Beijsterveldt CEM, Meulenbelt I, Slagboom PE, Boomsma DI. Sex differences in sum scores may be hard to interpret: the importance of measurement invariance. Assessment 2009; 16:415-23. [PMID: 19762518 DOI: 10.1177/1073191109344827] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
In most assessment instruments, distinct items are designed to measure a trait, and the sum score of these items serves as an approximation of an individual's trait score. In interpreting group differences with respect to sum scores, the instrument should measure the same underlying trait across groups (e.g., male/female, young/old). Differences with respect to the sum score should accurately reflect differences in the latent trait of interest. A necessary condition for this is that the instrument is measurement invariant. In the current study, the authors illustrate a stepwise approach for testing measurement invariance with respect to sex in a four-item instrument designed to assess disordered eating behavior in a large epidemiological sample (1,195 men and 1,507 women). This approach can be applied to other phenotypes for which group differences are expected. Any analysis of such variables may be subject to measurement bias if a lack of measurement invariance between grouping variables goes undetected.
Collapse
|
37
|
Kettunen J, Perola M, Martin NG, Cornes BK, Wilson SG, Montgomery GW, Benyamin B, Harris JR, Boomsma D, Willemsen G, Hottenga JJ, Slagboom PE, Christensen K, Kyvik KO, Sørensen TIA, Pedersen NL, Magnusson PKE, Andrew T, Spector TD, Widen E, Silventoinen K, Kaprio J, Palotie A, Peltonen L. Multicenter dizygotic twin cohort study confirms two linkage susceptibility loci for body mass index at 3q29 and 7q36 and identifies three further potential novel loci. Int J Obes (Lond) 2009; 33:1235-42. [PMID: 19721450 DOI: 10.1038/ijo.2009.168] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
OBJECTIVE To identify common loci and potential genetic variants affecting body mass index (BMI, kg m(-2)) in study populations originating from Europe. DESIGN We combined genome-wide linkage scans of six cohorts from Australia, Denmark, Finland, the Netherlands, Sweden and the United Kingdom with an approximately 10-cM microsatellite marker map. Variance components linkage analysis was carried out with age, sex and country of origin as covariates. SUBJECTS The GenomEUtwin consortium consists of twin cohorts from eight countries (Australia, Denmark, the Netherlands, Finland, Italy, Norway, Sweden and the United Kingdom) with a total data collection of more than 500,000 monozygotic and dizygotic (DZ) twin pairs. Variance due to early-life events and the environment is reduced within twin pairs, which makes DZ pairs highly valuable for linkage studies of complex traits. This study totaled 4401 European-originated twin families (10,535 individuals) from six countries (Australia, Denmark, the Netherlands, Finland, Sweden and the United Kingdom). RESULTS We found suggestive evidence for a quantitative trait locus on 3q29 and 7q36 in the combined sample of DZ twins (multipoint logarithm of odds score (MLOD) 2.6 and 2.4, respectively). Two individual cohorts showed strong evidence independently for three additional loci: 16q23 (MLOD=3.7) and 2p24 (MLOD=3.4) in the Dutch cohort and 20q13 (MLOD=3.2) in the Finnish cohort. CONCLUSION Linkage analysis of the combined data in this large twin cohort study provided evidence for suggestive linkage to BMI. In addition, two cohorts independently provided significant evidence of linkage to three new loci. The results of our study suggest a smaller environmental variance between DZ twins than full siblings, with a corresponding increase in heritability for BMI as well as an increase in linkage signal in well-replicated regions. The results are consistent with the possibility of locus heterogeneity for some genomic regions, and indicate a lack of major common quantitative trait locus variants affecting BMI in European populations.
Collapse
|
38
|
Bos SD, Kloppenburg M, Suchiman E, van Beelen E, Slagboom PE, Meulenbelt I. The role of plasma cytokine levels, CRP and Selenoprotein S gene variation in OA. Osteoarthritis Cartilage 2009; 17:621-6. [PMID: 19036615 DOI: 10.1016/j.joca.2008.10.007] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/28/2008] [Accepted: 10/18/2008] [Indexed: 02/02/2023]
Abstract
OBJECTIVE Investigating the association between plasma levels of cytokines and chemokines, Selenoprotein S (SELS) gene variation and osteoarthritis (OA) subtypes. METHODS The genetics of osteoarthritis and progression (GARP) study consists of 191 sibling pairs with symptomatic OA at multiple joint sites. We have measured plasma levels of 17 cytokines and chemokines and genetic variation at the SELS gene. RESULTS Nine out of 17 serum markers could be assessed quantitatively, whereas eight markers were assessed qualitatively. Principal component analysis (PCA) on the quantitatively assessed markers and serum high sensitive C-reactive protein (S-HsCRP) revealed that three components underlie 61% of the total plasma variation. Three single nucleotide polymorphisms (SNPs) in the SELS gene revealed four common haplotypes, one of which, GAG (frequency 3.5%) showed significant association to an anti-inflammatory (P=0.019) and acute phase related (P=0.036) component. OA subtype analysis showed that one component (mainly representing chemokine variation) was significantly associated to hand OA and disc degeneration (P=0.029 and P=0.010 respectively) as well as a physical component score (PCS) (P=0.042). The CRP related component also showed a strong association to the PCS (P=0.007). SELS haplotypes showed no association to OA subtypes in the GARP study. CONCLUSION Genetic variation in the SELS gene associates to components representing inflammatory signaling. Another component, representing chemokine variation, showed association to hand OA and disc degeneration in the GARP study indicating chemokines may contribute to OA pathogenesis.
Collapse
|
39
|
Albers HM, Kurreeman FAS, Stoeken-Rijsbergen G, Brinkman DMC, Kamphuis SSM, van Rossum MAJ, Girschick HJ, Wouters C, Saurenmann RK, Hoppenreijs E, Slagboom P, Houwing-Duistermaat JJ, Verduijn W, Huizinga TWJ, Ten Cate R, Toes REM, Schilham MW. Association of the autoimmunity locus 4q27 with juvenile idiopathic arthritis. ACTA ACUST UNITED AC 2009; 60:901-4. [PMID: 19248117 DOI: 10.1002/art.24296] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
OBJECTIVE Juvenile idiopathic arthritis (JIA) is characterized by chronic arthritis and an autoimmune etiology. In several autoimmune diseases, including rheumatoid arthritis (RA), an association with the 4q27 locus has been reported. We undertook this study to investigate the possible role of the 4q27 locus in JIA. METHODS A case-control association study was conducted, with a total of 655 Caucasian JIA patients and 791 healthy controls divided into 2 independent sample sets. The rs6822844 marker in the 4q27 locus was genotyped. RESULTS In the first and larger sample set, a 5% decrease in T allele frequency was observed in patients compared with controls (allelic odds ratio [OR] 0.72 [95% confidence interval 0.55-0.95], P = 0.019), and in the second set, a 3% decrease was observed (allelic OR 0.81 [95% confidence interval 0.61-1.09], P = 0.169). The combined data set generated an OR of 0.76 (95% confidence interval 0.62-0.93, P = 7.08 x 10(-3)). When the different JIA subtypes were analyzed individually, significant decreases were seen in the subtypes with a polyarticular course of disease (extended oligoarthritis [P = 0.019] and rheumatoid factor-negative polyarthritis [P = 0.038]). CONCLUSION Our findings suggest that the 4q27 locus, previously reported to be associated with RA, type 1 diabetes mellitus, celiac disease, and psoriatic arthritis, is also associated with susceptibility to JIA.
Collapse
|
40
|
Houwing-Duistermaat JJ, Callegaro A, Beekman M, Westendorp RG, Slagboom PE, van Houwelingen JC. Weighted statistics for aggregation and linkage analysis of human longevity in selected families: the Leiden Longevity Study. Stat Med 2009; 28:140-51. [PMID: 18759372 DOI: 10.1002/sim.3421] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Typically long-lived sibling pairs have been collected for linkage analysis of human longevity and information on life span of first-degree relatives is available to assess familial aggregation of life span. We propose a new weighted statistic for aggregation analysis, which tests for a relationship between a family history of excessive survival of the sibships of the long-lived pairs and the survival of their parents and their offspring. For linkage analysis, we derive a new weighted score statistic from a simple gamma frailty model, which assigns more weight to excessive long-lived pairs. We apply the methods to data from the Leiden Longevity Study, which consists of sibling pairs of age 90 years or above and their first-degree relatives. The pairs have been genotyped for microsatellite markers in a candidate region. Association was present between survival within the sibships and survival of the offspring, but not with the parental generation. For linkage analysis, weighting increased the value of the test statistic, but the result was not statistically significant. About the methods we conclude that the statistic for aggregation provides insight into clustering of life span and the statistic for linkage provides a new tool to include demographic information into the analysis.
Collapse
|
41
|
Rodriguez-Lopez J, Pombo-Suarez M, Loughlin J, Tsezou A, Blanco FJ, Meulenbelt I, Slagboom PE, Valdes AM, Spector TD, Gomez-Reino JJ, Gonzalez A. Association of a nsSNP in ADAMTS14 to some osteoarthritis phenotypes. Osteoarthritis Cartilage 2009; 17:321-7. [PMID: 18790654 DOI: 10.1016/j.joca.2008.07.012] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/08/2008] [Accepted: 07/29/2008] [Indexed: 02/02/2023]
Abstract
OBJECTIVE To investigate the effect in OA (Osteoarthritis) susceptibility of putative damaging changes in ADAM (A Disintegrin And Metalloprotease) and ADAMTS (ADAM with ThromboSpondin motif) proteases. METHODS Non-synonymous single nucleotide polymorphisms (nsSNP) in 18 ADAMTS and 31 ADAM genes were analyzed with two software applications for prediction of functional damage. Four putative damaging nsSNP were found in ADAMTS2, ADAMTS14, ADAMTS16 and ADAM12, respectively. These nsSNPs were analyzed in case-control sample collections with a variety of phenotypes totalling 3217 OA patients and 2214 healthy controls, all of them Caucasians. RESULTS No statistically significant differences were found in ADAMTS2, ADAMTS16 and ADAM12 nsSNPs. Conversely, the rare allele of the rs4747096 nsSNP in ADAMTS14 was overrepresented in women requiring joint replacement because of knee OA (O.R.(M-H) (odds ratio. Mantel-Haenszel)=1.41, 95% C.I.=1.1-1.8; P=0.002) and in patients with symptomatic hand OA (O.R.=1.37, 95% C.I.=1.0-1.9; P=0.047). A non significant increase in the frequency of the same allele was also found in patients with hip OA requiring prosthesis (O.R.(M-H)=1.14, 95% C.I.=1.0-1.3; P=0.08). No association was found with other OA phenotypes. CONCLUSION Our findings implicate ADAMTS14 in OA, specifically in knee OA requiring joint replacement in women and, possibly, in hand OA. Independent association of ADAMTS14 genetic variation to knee OA in women has been communicated. ADAMTS14 involvement, if confirmed, will open a new area of interest in OA pathogenesis because of its role in the maturation of collagen fibers.
Collapse
|
42
|
Slof-Op 't Landt MCT, Bartels M, van Furth EF, van Beijsterveldt CEM, Meulenbelt I, Slagboom PE, Boomsma DI. Genetic influences on disordered eating behaviour are largely independent of body mass index. Acta Psychiatr Scand 2008; 117:348-56. [PMID: 18081919 DOI: 10.1111/j.1600-0447.2007.01132.x] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
OBJECTIVE Prior studies suggest eating disorders and related characteristics are moderately to substantially heritable. We are interested in identifying the genes underlying disordered eating behaviour (DEB), and want to know how much of the genetic influence underlying DEB is attributable to genetic influences on body mass index (BMI). METHOD Bivariate analyses were performed, in adolescent twins and siblings, to estimate the genetic and environmental contributions for DEB, BMI, and their overlap. RESULTS Shared genetic risk factors explained the overlap between BMI and DEB (genetic correlation was 0.43 in women, 0.51 in men). DEB was highly heritable in women (a(2) = 0.65; a(2) independent of BMI = 0.53) and moderately heritable in men (a(2) = 0.39; a(2) independent of BMI = 0.29). BMI was highly heritable in both men (a(2) = 0.76) and women (a(2) = 0.80). CONCLUSION The entire correlation between DEB and BMI was explained by shared genetic risk, but the majority of genetic influences on DEB were due to genetic effects independent of BMI.
Collapse
|
43
|
Trompet S, de Craen AJM, Slagboom P, Shepherd J, Blauw GJ, Murphy MB, Bollen ELEM, Buckley BM, Ford I, Gaw A, Macfarlane PW, Packard CJ, Stott DJ, Westendorp RGJ, Jukema JW. Lymphotoxin-alpha C804A polymorphism is a risk factor for stroke. The PROSPER study. Exp Gerontol 2008; 43:801-5. [PMID: 18504081 DOI: 10.1016/j.exger.2008.04.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2008] [Revised: 04/08/2008] [Accepted: 04/15/2008] [Indexed: 11/30/2022]
Abstract
Inflammation plays a prominent role in the development of atherosclerosis, which is the most important risk factor for vascular events. Lymphotoxin-alpha (LTA) is a pro-inflammatory cytokine and is found to be expressed in atherosclerotic lesions. We investigated the association between the C804A polymorphism within the LTA gene and coronary and cerebrovascular events in 5804 participants of the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER). The primary endpoint was the combined endpoint of death from coronary heart disease, non-fatal myocardial infarction, and clinical stroke. Secondary endpoints were the coronary and cerebrovascular components separately. All associations were assessed with a Cox-proportional hazards model adjusted for sex, age, pravastatin use, and country. Our overall analysis showed a significant association between the C804A polymorphism and the primary endpoint (p = 0.03). After stratification for gender, this association was found only in males. Furthermore, we found that the association between the C804A polymorphism and the primary endpoint was mainly attributable to clinical strokes (p = 0.02). The C804A polymorphism in the LTA gene associates with clinical stroke, especially in men. But further research is warranted to confirm our results.
Collapse
|
44
|
Shifman S, Bhomra A, Smiley S, Wray NR, James MR, Martin NG, Hettema JM, An SS, Neale MC, van den Oord EJCG, Kendler KS, Chen X, Boomsma DI, Middeldorp CM, Hottenga JJ, Slagboom PE, Flint J. A whole genome association study of neuroticism using DNA pooling. Mol Psychiatry 2008; 13:302-12. [PMID: 17667963 PMCID: PMC4004964 DOI: 10.1038/sj.mp.4002048] [Citation(s) in RCA: 128] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
We describe a multistage approach to identify single nucleotide polymorphisms (SNPs) associated with neuroticism, a personality trait that shares genetic determinants with major depression and anxiety disorders. Whole genome association with 452 574 SNPs was performed on DNA pools from approximately 2000 individuals selected on extremes of neuroticism scores from a cohort of 88 142 people from southwest England. The most significant SNPs were then genotyped on independent samples to replicate findings. We were able to replicate association of one SNP within the PDE4D gene in a second sample collected by our laboratory and in a family-based test in an independent sample; however, the SNP was not significantly associated with neuroticism in two other independent samples. We also observed an enrichment of low P-values in known regions of copy number variations. Simulation indicates that our study had approximately 80% power to identify neuroticism loci in the genome with odds ratio (OR)>2, and approximately 50% power to identify small effects (OR=1.5). Since we failed to find any loci accounting for more than 1% of the variance, the heritability of neuroticism probably arises from many loci each explaining much less than 1%. Our findings argue the need for much larger samples than anticipated in genetic association studies and that the biological basis of emotional disorders is extremely complex.
Collapse
|
45
|
Trompet S, de Craen AJM, Slagboom P, Shepherd J, Blauw GJ, Murphy MB, Bollen ELEM, Buckley BM, Ford I, Gaw A, Macfarlane PW, Packard CJ, Stott DJ, Jukema JW, Westendorp RGJ. Genetic variation in the interleukin-1 beta-converting enzyme associates with cognitive function. The PROSPER study. Brain 2008; 131:1069-77. [PMID: 18304957 DOI: 10.1093/brain/awn023] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Inflammation is thought to play an important role in the development of cognitive decline and dementia in old age. The interleukin-1 signalling pathway may play a prominent role in this process. The gene encoding for interleukin-1 beta-converting enzyme (ICE) is likely to influence IL-1 beta levels. Inhibition of ICE decreases the age-related increase in IL-1 beta levels and may therefore improve memory function. We assessed whether genetic variation in the ICE gene associates with cognitive function in an elderly population. All 5804 participants of the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER) were genotyped for the 10643GC, 9323GA, 8996AG and 5352GA polymorphisms in the ICE gene. Cross-sectional associations between the polymorphisms and cognitive function were assessed with linear regression. Longitudinal associations between polymorphisms, haplotypes and cognitive function were assessed with linear mixed models. All associations were adjusted for sex, age, education, country, treatment with pravastatin and version of test where appropriate. Subjects carrying the variants 10643C and 5352A allele had significantly lower IL-1 beta production levels (P < 0.01). Furthermore, we demonstrated that homozygous carriers of the 10643C and the 5352A allele performed better on all executive function tests at baseline and during follow-up compared to homozygous carriers of the wild-type allele (all P < 0.02). The haplotype with two variants present (10643C and 5352A) was associated with better executive function (all P < 0.02) compared to the reference haplotype without variants. For memory function the same trend was observed, although not significant. Genetic variation in the ICE gene is associated with better performance on cognitive function and lower IL-1 beta production levels. This suggests that low levels of IL-1 beta are protective for memory and learning deficits. Inhibition of ICE may therefore be an important therapeutic target for maintaining cognitive function in old age.
Collapse
|
46
|
Passtoors WM, Beekman M, Gunn D, Boer JM, Heijmans BT, Westendorp RGJ, Zwaan BJ, Slagboom PE. Genomic studies in ageing research: the need to integrate genetic and gene expression approaches. J Intern Med 2008; 263:153-66. [PMID: 18226093 DOI: 10.1111/j.1365-2796.2007.01904.x] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Genome-wide and hypothesis-based approaches to the study of ageing and longevity have been dominated by genetic investigations. To identify essential mechanisms of a complex trait such as ageing in higher species, a holistic understanding of interacting pathways is required. More information on such interactions is expected to be obtained from global gene expression analysis if combined with genetic studies. Genetic sequence variation often provides a functional gene marker for the trait, whereas a gene expression profile may provide a quantitative biomarker representing complex cellular pathway interactions contributing to the trait. Thus far, gene expression studies have associated multiple pathways to ageing including mitochondrial electron transport and the oxidative stress response. However, most of the studies are underpowered to detect small age-changes. A systematic survey of gene expression changes as a function of age in human individuals and animal models is lacking. Well designed gene expression studies, especially at the level of biological processes, will provide hypotheses on gene-environmental interactions determining biological ageing rate. Cross-sectional studies monitoring the profile as a chronological marker of ageing must be integrated with prospective studies indicating which profiles represent biomarkers preceding and predicting physiological decline and mortality. New study designs such as the Leiden Longevity Study, including two generations of subjects from longevity families, aim to achieve these combined approaches.
Collapse
|
47
|
Middeldorp CM, Hottenga JJ, Slagboom PE, Sullivan PF, de Geus EJC, Posthuma D, Willemsen G, Boomsma DI. Linkage on chromosome 14 in a genome-wide linkage study of a broad anxiety phenotype. Mol Psychiatry 2008; 13:84-9. [PMID: 17700576 PMCID: PMC4205275 DOI: 10.1038/sj.mp.4002061] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Several linkage studies on anxiety have been carried out in samples ascertained through probands with panic disorder. The results indicated that using a broad anxiety phenotype instead of a DSM-IV anxiety disorder diagnosis might enhance the chance of finding a linkage signal. In the current study, a genome-wide linkage analysis was performed on anxiety measured with a self-report questionnaire whose scores are highly correlated with DSM-IV anxiety disorders. The self-report questionnaire was included in five surveys of a longitudinal study of the Netherlands Twin Register. Genotype and phenotype data were available for 1602 twins and siblings. To estimate identity by descent , additional genotype data for 564 parents and 22 siblings were used. Linkage analyses were carried out using MERLIN-regress on the average anxiety scores across time. A linkage signal (logarithm of odds score 3.4, empirical P-value 0.07) was obtained at chromosome 14 for marker D14S65 at 105 cM (90% confidence interval, 99-115 cM bounded by markers D14S1434 and D14S985). This finding replicates a linkage finding for a broad anxiety phenotype in a clinically based sample, indicating that the region might harbor a quantitative trait locus associated with the whole spectrum of general anxiety, that is from the normal to the clinical range. Moreover, genome-wide linkage and association studies on emotionality in mice obtained significant results in a syntenic region on mouse chromosome 12. Two homolog genes lie in this region -Dlk1 (delta-like 1 homolog, Drosophila) and Rtl1 (retrotransposon-like 1). Future association studies of these genes are warranted.
Collapse
|
48
|
Colas F, Meulenbelt I, Houwing-Duistermaat JJ, Kloppenburg M, Watt I, van Rooden SM, Visser M, Marinus J, van Hilten JJ, Slagboom P, Kok JN. Reliability of cluster results for different types of time adjustments in complex disease research. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2008; 2008:4601-4604. [PMID: 19163741 DOI: 10.1109/iembs.2008.4650238] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
We aim to identify subtypes of diseases like Osteoarthritis (OA) and Parkinson's Disease (PD) that present clinical heterogeneity. We do so by searching for homogeneous clusters in values of markers that reflect the severity of the disease. In the current paper we consider two important items for a cluster analysis. First, as time can contribute largely to the measured variability in the data, we search for the most appropriate way to adjust for it. Second, as we aim for reliable cluster analyses, cluster results should exhibit robustness to little change in the data. To investigate these issues, we transform the data by adding noise of different levels before cluster modeling and we rely on a chi(2)-based measure of association to compare cluster results for different types of time adjustment. The results of our experiments suggest to adjust data for a logarithmic age effect for OA and a square root effect of the disease duration for PD because these adjustments lead more reliable cluster results.
Collapse
|
49
|
Bos SD, Suchiman HED, Kloppenburg M, Houwing-Duistermaat JJ, le Graverand MPH, Seymour AB, Kroon HM, Slagboom PE, Meulenbelt I. Allelic variation at the C-reactive protein gene associates to both hand osteoarthritis severity and serum high sensitive C-reactive protein levels in the GARP study. Ann Rheum Dis 2007; 67:877-9. [PMID: 18055473 DOI: 10.1136/ard.2007.079228] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
OBJECTIVE To gain more insight into the role of genetic variation of the C-reactive protein (CRP) gene in serum CRP levels and osteoarthritis (OA). METHODS Serum high sensitive CRP (S-HsCRP) levels were measured in the Genetics of osteoARthritis and Progression (GARP) study. Furthermore, to assess genetic variation of the CRP gene, genotypes of five tagging single nucleotide polymorphisms were assessed in the GARP study and a random control sample. RESULTS A significant and consistent relation between S-HsCRP levels and observed haplotypes was identified. Additionally, a CRP haplotype, which also associated to a significantly higher expected phenotypic mean S-HsCRP level, was associated to severe hand OA. This haplotype was tagged by a single nucleotide polymorphism (rs3091244). Carriers of this allele have an increased risk for the presence of severe hand OA with an OR of 2.3 (95% confidence interval 1.2 to 4.3, p = 0.009). CONCLUSIONS A haplotype of the CRP gene, associated to high basal S-HsCRP level, is also associated to severity of hand OA, indicating that innate high basal S-HsCRP levels may influence OA onset.
Collapse
|
50
|
Trompet S, Pons D, DE Craen AJM, Slagboom P, Shepherd J, Blauw GJ, Murphy MB, Cobbe SM, Bollen ELEM, Buckley BM, Ford I, Hyland M, Gaw A, Macfarlane PW, Packard CJ, Norrie J, Perry IJ, Stott DJ, Sweeney BJ, Twomey C, Westendorp RGJ, Jukema JW. Genetic Variation in the Interleukin-10 Gene Promoter and Risk of Coronary and Cerebrovascular Events: The PROSPER Study. Ann N Y Acad Sci 2007; 1100:189-98. [PMID: 17460178 DOI: 10.1196/annals.1395.018] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Proinflammatory cytokines, like interleukin-6 (IL-6) and tumor necrosis factor-alpha (TNF-alpha), are implicated in the development of atherosclerosis. The role of anti-inflammatory cytokines, like IL-10, is largely unknown. We investigated the association of four single nucleotide polymorphisms (SNPs) in the promoter region of the IL-10 gene (4259AG, -1082GA, -592CA, and -2849GA), with coronary and cerebrovascular disease in participants of the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER) trial. All associations were assessed with Cox proportional hazards models adjusted for sex, age, pravastatin use, and country. Haplotype analysis of the four SNPs showed a significant association between haplotype 4 (containing the -592A variant allele) and risk of coronary events (P = 0.019). Moreover, analysis of separate SNPs found a significant association between -2849AA carriers with incident stroke (HR (95%CI) 1.50 (1.04-2.17), P value = 0.02). Our study suggests that not only proinflammatory processes contribute to atherosclerosis, but that also anti-inflammatory cytokines may play an important role.
Collapse
|