MiMIR: R-shiny application to infer risk factors and endpoints from Nightingale Health's 1H-NMR Metabolomics data

MOTIVATION

1H-NMR metabolomics is rapidly becoming a standard resource in large epidemiological studies to acquire metabolic profiles in large numbers of samples in a relatively low-priced and standardized manner. Concomitantly, metabolomics-based models are increasingly developed that capture disease risk or clinical risk factors. These developments raise the need for user-friendly toolbox to inspect new 1H-NMR metabolomics data and project a wide array of previously established risk models.

RESULTS

We present MiMIR (Metabolomics-based Models for Imputing Risk), a graphical user interface that provides an intuitive framework for ad-hoc statistical analysis of Nightingale Health's 1H-NMR metabolomics data and allows for the projection and calibration of 24 pre-trained metabolomics-based models, without any pre-required programming knowledge.

AVAILABILITY

The R-shiny package is available in CRAN or downloadable at https://github.com/DanieleBizzarri/MiMIR, together with an extensive user manual (also available as Supplementary Documents to the paper).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Variants of FOXO3 and RPA3 genes affecting IGF-1 levels alter the risk of development of primary osteoarthritis

INTRODUCTION

Pathologically high growth hormone (GH) and insulin-like growth factor-1 (IGF-1) levels in patients with acromegaly are associated with arthropathy. Several studies highlight the potential role of the GH/IGF-1 axis in primary osteoarthritis (OA). We aimed to disentangle the role of IGF-1 levels in primary OA pathogenesis.

METHODS

Patients from the Genetics osteoARthritis and Progression (GARP) Study with familial, generalized, symptomatic OA (n = 337, mean age: 59.8 ± 7.4 years, 82% female) were compared to Leiden Longevity Study (LLS) controls (n = 456, mean age: 59.8 ± 6.8 years, 51% female). Subjects were clinically and radiographically assessed, serum IGF-1 levels were measured, and 10 quantitative trait loci (QTL) in the FOXO3, IGFBP3/TNS3, RPA3, SPOCK2 genes, previously related to serum IGF-1 levels, were genotyped. Linear or binary logistic generalized estimating equation models were performed.

RESULTS

Serum IGF-1 levels were increased in OA patients, with male patients exhibiting the strongest effect (males OR = 1.10 (1.04-1.17), P=0.002 vs females OR = 1.04 (1.01-1.07), P = 0.02). Independent of the increased IGF-1 levels, male carriers of the minor allele of FOXO3 QTL rs4946936 had a lower risk to develop hip OA (OR = 0.41 (0.18-0.90), P = 0.026). Additionally, independent of IGF-1 levels, female carriers of the minor alleles of RPA3 QTL rs11769597 had a higher risk to develop knee OA (OR = 1.90 (1.20-2.99), P = 0.006).

CONCLUSION

Patients with primary OA had significantly higher IGF-1 levels compared to controls. Moreover, SNPs in the FOXO3 and RPA3 genes were associated with an altered risk of OA. Therefore, altered IGF-1 levels affect the development of OA, and are potentially the result of the pathophysiological OA process.

Basal cell carcinoma genetic susceptibility increases the rate of skin ageing: a Mendelian randomization study

BACKGROUND

Onset of basal cell carcinoma (BCC) is connected to skin ageing, but it is unclear whether higher BCC genetic susceptibility drives skin ageing.

OBJECTIVES

To investigate whether loci increasing genetic susceptibility to BCC also drive multiple features of skin ageing, independently of confounding factors, using Mendelian randomization.

METHODS

A Mendelian randomization study was conducted in older adults from the Leiden Longevity Study (N = 604). A total of 25 BCC loci, selected based on a published genome-wide association study on BCC (P-value < 5 × 10^-8 ), were used as genetic instruments for the calculation of a standardized (mean = 0, SD = 1) weighted BCC genetic risk score. Based on facial photographs, we determined perceived age, and skin wrinkling and pigmented spot grading.

RESULTS

A higher BCC genetic risk score was associated with a higher perceived age (adjusted for chronological age and sex) of 0.88 years (95% CI: 0.44, 1.31; P-value = 7.1e^-5 ), greater wrinkling by 0.14 grades (95% CI: 0.05, 0.23; P-value = 2.3e^-3 ), and greater pigmented spots by 0.17 grades (95% CI: 0.08, 0.25; P-value = 1.1e^-4 ). These findings were weakened but still present after exclusion of gene variants in MC1R and IRF4 which have potential pleiotropic effects.

CONCLUSIONS

Mechanisms influenced by genetic loci increasing susceptibility to BCC also drive skin ageing suggesting shared biology and shared targets for interventions.

Men Have a Stronger Monocyte-Derived Cytokine Production Response upon Stimulation with the Gram-Negative Stimulus Lipopolysaccharide than Women: A Pooled Analysis Including 15 Study Populations

The incidence of bacterial infections and sepsis, as well as the mortality risk from sepsis, is sex specific. These clinical findings have been attributed to sex differences in immune responsiveness. The aim of the present study was to investigate sex differences in monocyte-derived cytokine production response upon stimulation with the gram-negative stimulus lipopolysaccharide (LPS) using cytokine data from 15 study populations. Individual data on ex vivo cytokine production response upon stimulation with LPS in whole blood were available for 4,020 subjects originating from these 15 study populations, either from the general population or from patient populations with specific diseases. Men had a stronger cytokine production response than women to LPS for tumour necrosis factor-α, interleukin (IL)-6, IL-12, IL-1β, IL-1RA, and IL-10, but not for interferon-γ. The granulocyte-macrophage colony-stimulating factor production response was lower in men than in women. These sex differences were independent of chronological age. As men had higher monocyte concentrations, we normalized the cytokine production responses for monocyte concentration. After normalization, the sex differences in cytokine production response to LPS disappeared, except for IL-10, for which the production response was lower in men than in women. A sex-based approach to interpreting immune responsiveness is crucial.

Systemic Age-Associated DNA Hypermethylation of ELOVL2 Gene: In Vivo and In Vitro Evidences of a Cell Replication Process

Epigenetic remodeling is one of the major features of the aging process. We recently demonstrated that DNA methylation of ELOVL2 and FHL2 CpG islands is highly correlated with age in whole blood. Here we investigated several aspects of age-associated hypermethylation of ELOVL2 and FHL2. We showed that ELOVL2 methylation is significantly different in primary dermal fibroblast cultures from donors of different ages. Using epigenomic data from public resources, we demonstrated that most of the tissues show ELOVL2 and FHL2 hypermethylation with age. Interestingly, ELOVL2 hypermethylation was not found in tissues with very low replication rate. We demonstrated that ELOVL2 hypermethylation is associated with in vitro cell replication rather than with senescence. We confirmed intra-individual hypermethylation of ELOVL2 and FHL2 in longitudinally assessed participants from the Doetinchem Cohort Study. Finally we showed that, although the methylation of the two loci is not associated with longevity/mortality in the Leiden Longevity Study, ELOVL2 methylation is associated with cytomegalovirus status in nonagenarians, which could be informative of a higher number of replication events in a fraction of whole-blood cells. Collectively, these results indicate that ELOVL2 methylation is a marker of cell divisions occurring during human aging.

Neo-cartilage engineered from primary chondrocytes is epigenetically similar to autologous cartilage, in contrast to using mesenchymal stem cells

OBJECTIVES

To compare the epigenetic landscape of 3D cell models of human primary articular chondrocytes (hPACs) and human bone-marrow derived mesenchymal stem cells (hBMSCs) and their respective autologous articular cartilage.

DESIGN

Using Illumina Infinium HumanMethylation450 BeadChip arrays, the DNA methylation landscape of the different cell sources and autologous cartilage was determined. Pathway enrichment was analyzed using DAVID.

RESULTS

Principal Component Analysis (PCA) of methylation data revealed separate clustering of hBMSC samples. Between hBMSCs and autologous cartilage 86,881 cytosine-phosphate-guanine dinucleotides (CpGs) (20.2%), comprising 3,034 differentially methylated regions (DMRs; Δβ > 0.1; with the same direction of effect), were significantly differentially methylated. In contrast, between hPACs and autologous cartilage only 5,706 CpGs (1.33%) were differentially methylated. Of interest was the finding of the transcriptionally active, hyper-methylation of a Cartilage Intermediate Layer Protein (CILP) annotated DMR (Δβ = 0.16) in PAC-cartilage, corresponding to a profound decrease in CILP expression after in vitro culturing of hPACs as compared to autologous cartilage.

CONCLUSIONS

In vitro engineered neo-cartilage tissue from primary chondrocytes, hPACs, exhibits a DNA methylation landscape that is almost identical (99% similarity) to autologous cartilage, in contrast to neo-cartilage engineered from bone marrow-derived mesenchymal stem cells (MSCs). Although hBMSCs are widely used for cartilage engineering purposes the effects of these vast differences on cartilage regeneration and long term consequences of implantation, are not known. The use of hBMSCs or hPACs for future cartilage tissue regeneration purposes should therefore be investigated in more depth in future endeavors to better understand the consequences of the differential methylome on neo-cartilage.

Both low circulating insulin-like growth factor-1 and high-density lipoprotein cholesterol are associated with hair loss in middle-aged women

BACKGROUND

Multiple biomarkers have been associated with hair loss in women, but studies have shown inconsistent results.

OBJECTIVES

We investigated the associations between markers of cardiovascular disease risk (e.g. serum lipid levels and hypertension) and ageing [e.g. 25-hydroxyvitamin D and insulin-like growth factor (IGF)] with hair loss in a population of middle-aged women.

METHODS

In a random subgroup of 323 middle-aged women (mean age 61·5 years) from the Leiden Longevity Study, hair loss was graded by three assessors using the Sinclair scale; women with a mean score > 1·5 were classified as cases with hair loss.

RESULTS

Every 1 SD increase in high-density lipoprotein (HDL) cholesterol was associated with a 0·65-times lower risk [95% confidence interval (CI) 0·46-0·91] of hair loss. For IGF-1 the risk was 0·68 times lower (95% CI 0·48-0·97) per 1 SD increase, independently of the other studied variables. Women with both IGF-1 and HDL cholesterol levels below the medians of the study population had a 3·47-times higher risk (95% CI 1·30-9·25) of having hair loss.

CONCLUSIONS

Low HDL cholesterol and IGF-1 were associated with a higher risk of hair loss in women. However, further studies are required to infer causal relationships.

Getting personal: Endogenous adenosine receptor signaling in lymphoblastoid cell lines

Genetic differences between individuals that affect drug action form a challenge in drug therapy. Many drugs target G protein-coupled receptors (GPCRs), and a number of receptor variants have been noted to impact drug efficacy. This, however, has never been addressed in a systematic way, and, hence, we studied real-life genetic variation of receptor function in personalized cell lines. As a showcase we studied adenosine A2A receptor (A2AR) signaling in lymphoblastoid cell lines (LCLs) derived from a family of four from the Netherlands Twin Register (NTR), using a non-invasive label-free cellular assay. The potency of a partial agonist differed significantly for one individual. Genotype comparison revealed differences in two intron SNPs including rs2236624, which has been associated with caffeine-induced sleep disorders. While further validation is needed to confirm genotype-specific effects, this set-up clearly demonstrated that LCLs are a suitable model system to study genetic influences on A2AR response in particular and GPCR responses in general.

An Internet-Based Physical Activity Intervention to Improve Quality of Life of Inactive Older Adults: A Randomized Controlled Trial

Background

Increasing physical activity is a viable strategy for improving both the health and quality of life of older adults.

Objective

The aim of this study was to assess if an Internet-based intervention aimed to increase physical activity was effective in improving quality of life of inactive older adults. In addition, we analyzed the effect of the intervention on quality of life among those participants who successfully reached their individually targeted increase in daily physical activity as indicated by the intervention program, as well as the dose-response effect of increasing physical activity on quality of life.

Methods

The intervention was tested in a randomized controlled trial and was comprised of an Internet program—DirectLife (Philips)—aimed at increasing physical activity using monitoring and feedback by accelerometry and feedback by digital coaching (n=119). The control group received no intervention (n=116). Participants were inactive 60-70-year-olds and were recruited from the general population. Quality of life and physical activity were measured at baseline and after 3 months using the Research ANd Development 36-item health survey (RAND-36) and wrist-worn triaxial accelerometer, respectively.

Results

After 3 months, a significant improvement in quality of life was seen in the intervention group compared to the control group for RAND-36 subscales on emotional and mental health (2.52 vs -0.72, respectively; P=.03) and health change (8.99 vs 2.03, respectively; P=.01). A total of 50 of the 119 participants (42.0%) in the intervention group successfully reached their physical activity target and showed a significant improvement in quality of life compared to the control group for subscales on emotional and mental health (4.31 vs -0.72, respectively; P=.009) and health change (11.06 vs 2.03, respectively; P=.004). The dose-response analysis showed that there was a significant association between increase in minutes spent in moderate-to-vigorous physical activity (MVPA) and increase in quality of life.

Conclusions

Our study shows that an Internet-based physical activity program was effective in improving quality of life in 60-70-year-olds after 3 months, particularly in participants that reached their individually targeted increase in daily physical activity.

Trial Registration

Nederlands Trial Register: NTR 3045; http://www.trialregister.nl/trialreg/admin/rctview.asp?TC=3045 (Archived by WebCite at http://www.webcitation.org/6fobg2sjJ)

Non-response to (statin) therapy: the importance of distinguishing non-responders from non-adherers in pharmacogenetic studies

PURPOSE

In pharmacogenetic research, genetic variation in non-responders and high responders is compared with the aim to identify the genetic loci responsible for this variation in response. However, an important question is whether the non-responders are truly biologically non-responsive or actually non-adherent? Therefore, the aim of this study was to describe, within the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER), characteristics of both non-responders and high responders of statin treatment in order to possibly discriminate non-responders from non-adherers.

METHODS

Baseline characteristics of non-responders to statin therapy (≤10 % LDL-C reduction) were compared with those of high responders (>40 % LDL-C reduction) through a linear regression analysis. In addition, pharmacogenetic candidate gene analysis was performed to show the effect of excluding non-responders from the analysis.

RESULTS

Non-responders to statin therapy were younger (p = 0.001), more often smoked (p < 0.001), had a higher alcohol consumption (p < 0.001), had lower LDL cholesterol levels (p < 0.001), had a lower prevalence of hypertension (p < 0.001), and had lower cognitive function (p = 0.035) compared to subjects who highly responded to pravastatin treatment. Moreover, excluding non-responders from pharmacogenetic studies yielded more robust results, as standard errors decreased.

CONCLUSION

Our results suggest that non-responders to statin therapy are more likely to actually be non-adherers, since they have more characteristics that are viewed as indicators of high self-perceived health and low disease awareness, possibly making the subjects less adherent to study medication. We suggest that in pharmacogenetic research, extreme non-responders should be excluded to overcome the problem that non-adherence is investigated instead of non-responsiveness.

Effect of calendar age on physical performance: A comparison of standard clinical measures with instrumented measures in middle-aged to older adults

BACKGROUND

Decline in physical performance is highly prevalent during aging. Identification of sensitive markers of age-related changes in physical performance is important for early detection, development of therapeutic strategies and insight into underlying mechanisms. We studied the association of calendar age and familial longevity with standard clinical and instrumented measures of physical performance in a cohort of healthy middle-aged to older adults.

METHODS

Cross-sectional analysis within the Leiden Longevity Study consisting of offspring of nonagenarian siblings and their partners (n=300, mean age (SD) 65.3 (6.7) years). Standard clinical measures were 25-meter walking speed and total duration of the chair stand test (CST). Instrumented measures were determined using a body fixed sensor. Dependence of physical performance on calendar age and familial longevity (offspring versus partner status) was analyzed using linear and logistic regression, respectively, adjusted for gender and height.

RESULTS

Higher calendar age was associated with slower walking speed and longer duration of the CST (standardized β (95% CI) -.024 (-.042; -.006) and .035 (.014;.056), respectively). Instrumented measures showed similar effect sizes with strongest associations for gait stability and symmetry in mediolateral direction and for the extension and flexion phase of sit-to-stand and stand-to-sit transfers, respectively. No differences were observed between offspring of nonagenarian siblings and their partners.

CONCLUSIONS

Standard clinical and instrumented measures of physical performance are associated with similar effect size to age-related changes in physical performance observable from middle age. The potential added value of instrumented measures for understanding underlying mechanisms requires further attention.

Copy number variation associates with mortality in long-lived individuals: a genome-wide assessment

Copy number variants (CNVs) represent a significant source of genetic variation in the human genome and have been implicated in numerous diseases and complex traits. To date, only a few studies have investigated the role of CNVs in human lifespan. To investigate the impact of CNVs on prospective mortality at the extreme end of life, where the genetic component of lifespan appears most profound, we analyzed genomewide CNV data in 603 Danish nonagenarians and centenarians (mean age 96.9 years, range 90.0–102.5 years). Replication was performed in 500 long‐lived individuals from the Leiden Longevity Study (mean age 93.2 years, range 88.9–103.4 years). First, we assessed the association between the CNV burden of each individual (the number of CNVs, the average CNV length, and the total CNV length) and mortality and found a significant increase in mortality per 10 kb increase in the average CNV length, both for all CNVs (hazard ratio (HR) = 1.024, P = 0.002) and for duplications (HR = 1.011, P = 0.005), as well as per 100 kb increase in the total length of deletions (HR = 1.009, P = 0.0005). Next, we assessed the relation between specific deletions and duplications and mortality. Although no genome–wide significant associations were discovered, we identified six deletions and one duplication that showed consistent association with mortality in both or either of the sexes across both study populations. These results indicate that the genome–wide CNV burden, specifically the average CNV length and the total CNV length, associates with higher mortality in long‐lived individuals.

GWAS for executive function and processing speed suggests involvement of the CADM2 gene

To identify common variants contributing to normal variation in two specific domains of cognitive functioning, we conducted a genome-wide association study (GWAS) of executive functioning and information processing speed in non-demented older adults from the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) consortium. Neuropsychological testing was available for 5429-32,070 subjects of European ancestry aged 45 years or older, free of dementia and clinical stroke at the time of cognitive testing from 20 cohorts in the discovery phase. We analyzed performance on the Trail Making Test parts A and B, the Letter Digit Substitution Test (LDST), the Digit Symbol Substitution Task (DSST), semantic and phonemic fluency tests, and the Stroop Color and Word Test. Replication was sought in 1311-21860 subjects from 20 independent cohorts. A significant association was observed in the discovery cohorts for the single-nucleotide polymorphism (SNP) rs17518584 (discovery P-value=3.12 × 10(-8)) and in the joint discovery and replication meta-analysis (P-value=3.28 × 10(-9) after adjustment for age, gender and education) in an intron of the gene cell adhesion molecule 2 (CADM2) for performance on the LDST/DSST. Rs17518584 is located about 170 kb upstream of the transcription start site of the major transcript for the CADM2 gene, but is within an intron of a variant transcript that includes an alternative first exon. The variant is associated with expression of CADM2 in the cingulate cortex (P-value=4 × 10(-4)). The protein encoded by CADM2 is involved in glutamate signaling (P-value=7.22 × 10(-15)), gamma-aminobutyric acid (GABA) transport (P-value=1.36 × 10(-11)) and neuron cell-cell adhesion (P-value=1.48 × 10(-13)). Our findings suggest that genetic variation in the CADM2 gene is associated with individual differences in information processing speed.

Detecting dispersed duplications in high-throughput sequencing data using a database-free approach

MOTIVATION

Dispersed duplications (DDs) such as transposon element insertions and copy number variations are ubiquitous in the human genome. They have attracted the interest of biologists as well as medical researchers due to their role in both evolution and disease. The efforts of discovering DDs in high-throughput sequencing data are currently dominated by database-oriented approaches that require pre-existing knowledge of the DD elements to be detected.

RESULTS

We present DD_DETECTION, a database-free approach to finding DD events in high-throughput sequencing data. DD_DETECTION is able to detect DDs purely from paired-end read alignments. We show in a comparative study that this method is able to compete with database-oriented approaches in recovering validated transposon insertion events. We also experimentally validate the predictions of DD_DETECTION on a human DNA sample, showing that it can find not only duplicated elements present in common databases but also DDs of novel type.

AVAILABILITY AND IMPLEMENTATION

The software presented in this article is open source and available from https://bitbucket.org/mkroon/dd_detection.

Human longevity is characterised by high thyroid stimulating hormone secretion without altered energy metabolism

Few studies have included subjects with the propensity to reach old age in good health, with the aim to disentangle mechanisms contributing to staying healthier for longer. The hypothalamic-pituitary-thyroid (HPT) axis maintains circulating levels of thyroid stimulating hormone (TSH) and thyroid hormone (TH) in an inverse relationship. Greater longevity has been associated with higher TSH and lower TH levels, but mechanisms underlying TSH/TH differences and longevity remain unknown. The HPT axis plays a pivotal role in growth, development and energy metabolism. We report that offspring of nonagenarians with at least one nonagenarian sibling have increased TSH secretion but similar bioactivity of TSH and similar TH levels compared to controls. Healthy offspring and spousal controls had similar resting metabolic rate and core body temperature. We propose that pleiotropic effects of the HPT axis may favour longevity without altering energy metabolism.

Longitudinal weight differences, gene expression and blood biomarkers in BMI-discordant identical twins

Background

BMI discordant monozygotic (MZ) twins allows an examination of the causes and consequences of adiposity in a genetically controlled design. Few studies have examined longitudinal BMI discordance in MZ pairs.

Objectives

To study the development over time of BMI discordance in adolescent and adult MZ twin pairs, and to examine lifestyle, metabolic, inflammatory, and gene expression differences associated with concurrent and long-term BMI discordance in MZ pairs.

Subjects/Methods

BMI data from 2775 MZ twin pairs, collected in eight longitudinal surveys and a biobank project between 1991 and 2011, were analyzed to characterize longitudinal discordance. Lifestyle characteristics were compared within discordant pairs (ΔBMI ≥ 3 kg/m²) and biomarkers (lipids, glucose, insulin, CRP, fibrinogen, IL-6, TNF-α and sIL-6R and liver enzymes AST, ALT and GGT) and gene expression were compared in peripheral blood from discordant pairs who participated in the NTR biobank project.

Results

The prevalence of discordance ranged from 3.2% in 1991 (mean age=17, SD=2.4) to 17.4% (N=202 pairs) in 2009 (mean age=35, SD=15), and was 16.5% (N=174) among pairs participating in the biobank project (mean age=35, SD=12). Of 699 MZ with BMI data from 3-5 time points, 17 pairs (2.4%) were long-term discordant (at all available time points; mean follow-up range=6.4 years). Concurrently discordant pairs showed significant differences in self-ratings of which twin eats most (p=2.3×10⁻¹³), but not in leisure time exercise activity (p=0.28) and smoking (p>0.05). Ten out of 14 biomarkers showed significantly more unfavorable levels in the heavier of twin of the discordant pairs (p-values < 0.001); most of these biomarker differences were largest in longitudinally discordant pairs. No significant gene expression differences were identified, although high ranking genes were enriched for Gene Ontology (GO) terms highlighting metabolic gene regulation and inflammation pathways.

Conclusions

BMI discordance is uncommon in adolescent identical pairs but increases with higher pair-mean of BMI at older ages, although long-term BMI discordance is rare. In discordant pairs, the heavier twin had a more unfavorable blood biomarker profile than the genetically matched leaner twin, in support of causal effects of obesity.

Developmental determinants in non-communicable chronic diseases and ageing

Prenatal and peri-natal events play a fundamental role in health, development of diseases and ageing (Developmental Origins of Health and Disease (DOHaD)). Research on the determinants of active and healthy ageing is a priority to: (i) inform strategies for reducing societal and individual costs of an ageing population and (ii) develop effective novel prevention strategies. It is important to compare the trajectories of respiratory diseases with those of other chronic diseases.

Variants of the IL-10 gene associate with muscle strength in elderly from rural Africa: a candidate gene study

Recently, it has been shown that the capacity of the innate immune system to produce cytokines relates to skeletal muscle mass and strength in older persons. The interleukin-10 (IL-10) gene regulates the production capacities of IL-10 and tumour necrosis factor-α (TNF-α). In rural Ghana, IL-10 gene variants associated with different production capacities of IL-10 and TNF-α are enriched compared with Caucasian populations. In this setting, we explored the association between these gene variants and muscle strength. Among 554 Ghanaians aged 50 years and older, we determined 20 single nucleotide polymorphisms in the IL-10 gene, production capacities of IL-10 and TNF-α in whole blood upon stimulation with lipopolysaccharide (LPS) and handgrip strength as a proxy for skeletal muscle strength. We distinguished pro-inflammatory haplotypes associated with low IL-10 production capacity and anti-inflammatory haplotypes with high IL-10 production capacity. We found that distinct haplotypes of the IL-10 gene associated with handgrip strength. A pro-inflammatory haplotype with a population frequency of 43.2% was associated with higher handgrip strength (P = 0.015). An anti-inflammatory haplotype with a population frequency of 7.9% was associated with lower handgrip strength (P = 0.006). In conclusion, variants of the IL-10 gene contributing to a pro-inflammatory cytokine response associate with higher muscle strength, whereas those with anti-inflammatory response associate with lower muscle strength. Future research needs to elucidate whether these effects of variation in the IL-10 gene are exerted directly through its role in the repair of muscle tissue or indirectly through its role in the defence against infectious diseases.

A genome-wide association study of anorexia nervosa

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

Sex-specific effects of naturally occurring variants in the dopamine receptor D2 locus on insulin secretion and type 2 diabetes susceptibility

AIMS

Modulation of dopamine receptor D2 (DRD2) activity affects insulin secretion in both rodents and isolated pancreatic β-cells. We hypothesized that single nucleotide polymorphisms in the DRD2/ANKK1 locus may affect susceptibility to type 2 diabetes in humans.

METHODS

Four potentially functional variants in the coding region of the DRD2/ANKK1 locus (rs1079597, rs6275, rs6277, rs1800497) were genotyped and analysed for type 2 diabetes susceptibility in up to 25 000 people (8148 with type 2 diabetes and 17687 control subjects) from two large independent Dutch cohorts and one Danish cohort. In addition, 340 Dutch subjects underwent a 2-h hyperglycaemic clamp to investigate insulin secretion. Since sexual dimorphic associations related to DRD2 polymorphisms have been previously reported, we also performed a gender-stratified analysis.

RESULTS

rs1800497 at the DRD2/ANKK1 locus was associated with a significantly increased risk for type 2 diabetes in women (odds ratio 1.14 (1.06-1.23); P = 4.1*10⁴) but not in men (odds ratio 1.00 (95% CI 0.93-1.07); P = 0.92) or the combined group. Although rs1800497 was not associated with insulin secretion, we did find another single nucleotide polymorphism in this locus, rs6275, to be associated with increased first-phase glucose-stimulated insulin secretion in women (P = 5.5*10⁴) but again not in men (P = 0.34).

CONCLUSION

The present data identify DRD2/ANKK1 as a potential sex-specific type 2 diabetes susceptibility gene.

Meta-analysis on blood transcriptomic studies identifies consistently coexpressed protein-protein interaction modules as robust markers of human aging

The bodily decline that occurs with advancing age strongly impacts on the prospects for future health and life expectancy. Despite the profound role of age in disease etiology, knowledge about the molecular mechanisms driving the process of aging in humans is limited. Here, we used an integrative network-based approach for combining multiple large-scale expression studies in blood (2539 individuals) with protein–protein Interaction (PPI) data for the detection of consistently coexpressed PPI modules that may reflect key processes that change throughout the course of normative aging. Module detection followed by a meta-analysis on chronological age identified fifteen consistently coexpressed PPI modules associated with chronological age, including a highly significant module (P = 3.5 × 10⁻³⁸) enriched for ‘T-cell activation’ marking age-associated shifts in lymphocyte blood cell counts (R² = 0.603; P = 1.9 × 10⁻¹⁰). Adjusting the analysis in the compendium for the ‘T-cell activation’ module showed five consistently coexpressed PPI modules that robustly associated with chronological age and included modules enriched for ‘Translational elongation’, ‘Cytolysis’ and ‘DNA metabolic process’. In an independent study of 3535 individuals, four of five modules consistently associated with chronological age, underpinning the robustness of the approach. We found three of five modules to be significantly enriched with aging-related genes, as defined by the GenAge database, and association with prospective survival at high ages for one of the modules including ASF1A. The hereby-detected age-associated and consistently coexpressed PPI modules therefore may provide a molecular basis for future research into mechanisms underlying human aging.

Drug-gene interactions and the search for missing heritability: a cross-sectional pharmacogenomics study of the QT interval

Variability in response to drug use is common and heritable, suggesting that genome-wide pharmacogenomics studies may help explain the 'missing heritability' of complex traits. Here, we describe four independent analyses in 33 781 participants of European ancestry from 10 cohorts that were designed to identify genetic variants modifying the effects of drugs on QT interval duration (QT). Each analysis cross-sectionally examined four therapeutic classes: thiazide diuretics (prevalence of use=13.0%), tri/tetracyclic antidepressants (2.6%), sulfonylurea hypoglycemic agents (2.9%) and QT-prolonging drugs as classified by the University of Arizona Center for Education and Research on Therapeutics (4.4%). Drug-gene interactions were estimated using covariable-adjusted linear regression and results were combined with fixed-effects meta-analysis. Although drug-single-nucleotide polymorphism (SNP) interactions were biologically plausible and variables were well-measured, findings from the four cross-sectional meta-analyses were null (Pinteraction>5.0 × 10(-8)). Simulations suggested that additional efforts, including longitudinal modeling to increase statistical power, are likely needed to identify potentially important pharmacogenomic effects.

Genome-wide linkage scan in affected sibling pairs identifies novel susceptibility region for venous thromboembolism: Genetics In Familial Thrombosis study

BACKGROUND

Venous thromboembolism (VTE) is a multicausal disorder involving environmental and genetic risk factors. In many thrombophilic families the clustering of thrombotic events cannot be explained by known genetic risk factors, indicating that some remain to be discovered.

OBJECTIVES

We aimed to identify novel thrombosis susceptibility alleles in a large panel of small thrombophilic families: the Genetics In Familial Thrombosis (GIFT) study.

PATIENTS/METHODS

In the GIFT study, 201 families were recruited consisting of 438 siblings with an objectively confirmed VTE at a young age. Multipoint linkage analysis (402 SSR markers) and fine mapping were performed, followed by genotyping of tagging SNPs in positional candidate genes.

RESULTS

Established genetic risk factors such as factor V Leiden, ABO blood group non-O, prothrombin 20210A, fibrinogen gamma 10034T and deficiencies of antithrombin, protein C and protein S were more frequent in GIFT patients than in unselected VTE patients. Linkage supported the presence of novel thrombosis susceptibility loci on 7p21.3-22.2 (LOD score = 3.23) and Xq24-27.3 (LOD score = 1.95). Simulation analysis showed that the chr7 signal was genome-wide statistically significant (P = 0.022). Tagging SNPs (n = 157) in eight positional candidate genes (LOD drop 1.5 regions) were genotyped in GIFT patients and 332 healthy controls. Five chr7 SNPs associated with VTE. SNP THSD7A rs2074597 was responsible for part of the chr7 signal.

CONCLUSIONS

The GIFT panel is rich in established genetic risk factors for VTE, but genetic factors remain unidentified in many families. Genome-wide linkage failed to identify the previously established genetic risk factors for VTE, but identified a novel VTE susceptibility locus on chr7.

Pro-inflammatory capacity of classically activated monocytes relates positively to muscle mass and strength

In mice, monocytes that exhibit a pro-inflammatory profile enter muscle tissue after muscle injury and are crucial for clearance of necrotic tissue and stimulation of muscle progenitor cell proliferation and differentiation. The aim of this study was to test if pro-inflammatory capacity of classically activated (M1) monocytes relates to muscle mass and strength in humans. This study included 191 male and 195 female subjects (mean age 64.2 years (SD 6.4) and 61.9 ± 6.4, respectively) of the Leiden Longevity Study. Pro-inflammatory capacity of M1 monocytes was assessed by ex vivo stimulation of whole blood with Toll-like receptor (TLR) 4 agonist lipopolysaccharide (LPS) and TLR-2/1 agonist tripalmitoyl-S-glycerylcysteine (Pam₃Cys-SK₄), both M1 phenotype activators. Cytokines that stimulate M1 monocyte response (IFN-γ and GM-CSF) as well as cytokines that are secreted by M1 monocytes (IL-6, TNF-α, IL-12, and IL-1β) were measured. Analyses were adjusted for age, height, and body fat mass. Upon stimulation with LPS, the cytokine production capacity of INF-γ, GM-CSF, and TNF-α was significantly positively associated with lean body mass, appendicular lean mass and handgrip strength in men, but not in women. Upon stimulation with Pam₃Cys-SK₄, IL-6; TNF-α; and Il-1β were significantly positively associated with lean body mass and appendicular lean in women, but not in men. Taken together, this study shows that higher pro-inflammatory capacity of M1 monocytes upon stimulation is associated with muscle characteristics and sex dependent.

Association study of candidate genes for the progression of hand osteoarthritis

OBJECTIVE

Although a few consistent osteoarthritis (OA) susceptibility genes have been identified, little is known on OA progression. Since OA progression is clinically the most relevant phenotype, we investigate the association between asporin (ASPN), bone morphogenetic protein 5 (BMP5) and growth differentiation factor 5 (GDF5) polymorphisms and progression of hand OA.

METHODS

Single-nucleotide polymorphisms (SNPs) ASPN rs13301537, BMP5 rs373444 and GDF5 rs143383 were genotyped in 251 hand OA patients from the Genetics osteoARthritis and Progression (GARP) study and 725 controls. In a case-control comparison we assessed the association between these SNPs and radiographic progression of hand OA over 6 years, which was based on change in osteophytes or joint space narrowing (JSN), above the smallest detectable change. SNPs with suggestive evidence for association were further analysed for their effect on progression over 2 years, and for the mean change in osteophytes and JSN.

RESULTS

The minor allele of ASPN SNP rs13301537 was associated with hand OA progression over 6 years (odds ratio (OR) (95% CI) 1.49 (1.06-2.07); P = 0.020). The mean change in osteophytes and JSN was higher in carriers of the minor allele compared to homozygous carriers of the common allele with mean difference of 0.73 (95% CI - 0.07-1.56; P = 0.073) and 0.82 (95% CI 0.12-1.52; P = 0.022), respectively. An association with similar effect size was found between ASPN SNP rs13301537 and 2-year progression, and the mean change in osteophytes and JSN was significantly higher in homozygotes.

CONCLUSION

ASPN is associated with hand OA progression. This gives insight in the pathogenesis of hand OA progression and identified a potential target for therapeutic approaches.

Gene expression analysis of mTOR pathway: association with human longevity

mTOR signalling is implicated in the development of disease and in lifespan extension in model organisms. This pathway has been associated with human diseases such as diabetes and cancer, but has not been investigated for its impact on longevity per se. Here, we investigated whether transcriptional variation within the mTOR pathway is associated with human longevity using whole-blood samples from the Leiden Longevity Study. This is a unique cohort of Dutch families with extended survival across generations, decreased morbidity and beneficial metabolic profiles in middle-age. By comparing mRNA levels of nonagenarians and middle-aged controls, the mTOR signalling gene set was found to associate with old age (P = 4.6 × 10(-7)). Single gene analysis showed that seven of 40 mTOR pathway genes had a significant differential expression of at least 5%. Of these, the RPTOR (Raptor) gene was found to be differentially expressed also when the offspring of nonagenarians was compared with their spouses, indicating association with familial longevity in middle-age. This association was not explained by variation between the groups in the prevalence of type 2 diabetes and cancer or glucose levels. Thus, the mTOR pathway not only plays a role in the regulation of disease and aging in animal models, but also in human health and longevity.

Serum insulin-like growth factor 1 and facial ageing: high levels associate with reduced skin wrinkling in a cross-sectional study

BACKGROUND

Insulin-like growth factor (IGF)-1 is a growth factor that can influence fibroblast functioning, with effects including the inhibition of collagenases and the induction of collagen expression.

OBJECTIVES

To assess whether serum IGF-1, IGF-binding protein (IGFBP)3 and the ratio between IGF-1 and IGFBP3, as a measure of IGF-1 bioavailability, are associated with facial ageing and skin wrinkling.

METHODS

From a random sample comprising 617 subjects from the Leiden Longevity Study, perceived age and skin wrinkling were assessed from facial photographs, and IGF-1 and IGFBP3 were measured in serum. The associations were assessed using linear regression models, adjusted for chronological age, sex, body mass index, smoking and sun exposure.

RESULTS

Across tertiles of the ratio of IGF-1 to IGFBP3, and after adjusting for all potential confounding factors, the mean perceived age decreased from 60·6 years in the lowest tertile to 59·5 years in the highest (P = 0·045). Similarly, the mean skin wrinkling grade decreased from 4·8 in the lowest tertile to 4·5 in the highest (P = 0·011). Adding skin wrinkling as a covariate in the analysis between IGF-1 and perceived age diminished this association.

CONCLUSIONS

This study demonstrates that a higher ratio of IGF-1 to IGFBP3 associates with a lower perceived age, via its association with reduced skin wrinkling. Whether high IGF-1 levels actually delay the accumulation of skin wrinkling now needs investigating.

Brain tissue volumes in familial longevity: the Leiden Longevity Study

Atrophy is one of the major age-related changes in the brain. The absence of brain atrophy in elderly individuals reflects deceleration in the process of biological aging. Moreover, results from human twin studies suggest a large genetic influence on the variance of human brain tissue volumes. To investigate the association of brain volumes with exceptional longevity, we tested whether middle-aged to elderly offspring of nonagenarian siblings have larger brain volumes than their spouses using magnetic resonance imaging. No differences in whole brain, gray matter and white matter volume were found. These brain volumes were associated with chronological age in offspring and control subjects (all P < 0.001). Left amygdalar volume of the offspring was larger (P = 0.03) compared with control subjects [mean volume offspring (cm3) (95% confidence interval, CI) = 1.39 (1.36-1.42), mean volume control subjects (cm3) (95% CI) = 1.32 (1.29-1.35)]. Association of left amygdalar volume with familial longevity was particularly pronounced when offspring with the oldest long-lived parent were compared with control subjects (P = 0.01). Amygdalar volumes were not associated with chronological age in both groups. Our findings suggest that the observed association of a larger left amygdalar volume with familial longevity is not caused by a relative preservation of the left amygdala during the course of aging but most likely a result of early development caused by a genetic familial trait.

The number of p16INK4a positive cells in human skin reflects biological age

Cellular senescence is a defense mechanism in response to molecular damage which accumulates with aging. Correspondingly, the number of senescent cells has been reported to be greater in older than in younger subjects and furthermore associates with age-related pathologies. Inter-individual differences exist in the rate at which a person ages (biological age). Here, we studied whether younger biological age is related to fewer senescent cells in middle-aged individuals with the propensity for longevity, using p16INK4a as a marker for cellular senescence. We observed that a younger biological age associates with lower levels of p16INK4a positive cells in human skin.

Glucocorticoid receptor gene polymorphisms are associated with reduced first-phase glucose-stimulated insulin secretion and disposition index in women, but not in men

AIM

Glucocorticoids are efficacious anti-inflammatory agents, but, in susceptible individuals, these drugs may induce glucose intolerance and diabetes by affecting β-cell function and insulin sensitivity. We assessed whether polymorphisms in the glucocorticoid receptor gene NR3C1 associate with measures of β-cell function and insulin sensitivity derived from hyperglycaemic clamps in subjects with normal or impaired glucose tolerance.

METHODS

A cross-sectional cohort study was conducted in four academic medical centres in the Netherlands and Germany. Four hundred and forty-nine volunteers (188 men; 261 women) were recruited with normal glucose tolerance (n=261) and impaired glucose tolerance (n=188). From 2-h hyperglycaemic clamps, first- and second-phase glucose-stimulated insulin secretion, as well as insulin sensitivity index and disposition index, were calculated. All participants were genotyped for the functional NR3C1 polymorphisms N363S (rs6195), BclI (rs41423247), ER22/23EK (rs6189/6190), 9β A/G (rs6198) and ThtIIII (rs10052957). Associations between these polymorphisms and β-cell function parameters were assessed.

RESULTS

In women, but not in men, the N363S polymorphism was associated with reduced disposition index (P=1.06 10(-4) ). Also only in women, the ER22/23EK polymorphism was associated with reduced first-phase glucose-stimulated insulin secretion (P=0.011) and disposition index (P=0.003). The other single-nucleotide polymorphisms were not associated with β-cell function. Finally, none of the polymorphisms was related to insulin sensitivity.

CONCLUSION

The N363S and ER22/23EK polymorphisms of the NR3C1 gene are negatively associated with parameters of β-cell function in women, but not in men.

Zygosity diagnosis in young twins by parental report

This study reports on zygosity determination in twins of childhood age. Parents responded to questionnaire items dealing with twin similarity in physical characteristics and frequency of mistaking one twin for another by parents, relatives and strangers. The accuracy of zygosity diagnosis was evaluated across twins aged 6, 8, and 10 and across parents. In addition, it was examined whether the use of multiple raters and the use of longitudinal data lead to an improvement of zygosity assignment. Complete data on zygosity questions and on genetic markers or blood profiles were available for 618 twin pairs at the age of 6 years. The method used was predictive discriminant analyses. Agreement between zygosity assigned by the replies to the questions and zygosity determined by DNA markers/blood typing was around 93%. The accuracy of assignment remained constant across age and parents. Analyses of data provided by both parents and collected over multiple ages did not result in better prediction of zygosity. Details on the discriminant function are provided. Twin Research (2000) 3, 134–141.

Proton magnetic resonance spectroscopy shows lower intramyocellular lipid accumulation in middle-aged subjects predisposed to familial longevity

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.

Hypermethylation at loci sensitive to the prenatal environment is associated with increased incidence of myocardial infarction

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.

A mixture model with random-effects components for classifying sibling pairs

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.

Body mass index and alignment and their interaction as risk factors for progression of knees with radiographic signs of osteoarthritis

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.

Accuracy of direct segmental multi-frequency bioimpedance analysis in the assessment of total body and segmental body composition in middle-aged adult population

BACKGROUND & AIMS

Body composition measurement is a valuable tool for assessing nutritional status and physical fitness in a variety of clinical settings. Although bioimpedance analysis (BIA) can easily assess body composition, its accuracy remains unclear. We examined the accuracy of direct segmental multi-frequency BIA technique (DSM-BIA) in assessing different body composition parameters, using dual energy X-ray absorptiometry (DEXA) as a reference standard.

METHODS

A total of 484 middle-aged participants from the Leiden Longevity Study were recruited. Agreements between DSM-BIA and DEXA for total and segmental body composition quantification were assessed using intraclass correlation coefficients and Bland-Altman plots.

RESULTS

Excellent agreements were observed between both techniques in whole body lean mass (ICC female = 0.95, ICC men = 0.96), fat mass (ICC female = 0.97, ICC male = 0.93) and percentage body fat (ICC female = 0.93, ICC male = 0.88) measurements. Similarly, Bland-Altman plots revealed narrow limits of agreements with small biases noted for the whole body lean mass quantification but relatively wider limits for fat mass and percentage body fat quantifications. In segmental lean muscle mass quantification, excellent agreements between methods were demonstrated for the upper limbs (ICC female≥0.91, ICC men≥0.87) and lower limbs (ICC female≥0.83, ICC male≥0.85), with good agreements shown for the trunk measurements (ICC female = 0.73, ICC male = 0.70).

CONCLUSIONS

DSM-BIA is a valid tool for the assessments of total body and segmental body composition in the general middle-aged population, particularly for the quantification of body lean mass.

Osteoarthritis susceptibility genes influence the association between hip morphology and osteoarthritis

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.

Genomics of human longevity

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.

Insights into the genetic architecture of osteoarthritis from stage 1 of the arcOGEN study

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.

Familial longevity is associated with decreased thyroid function

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.

Allele-sharing statistics using information on family history

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.

Anxiety and depression in children and adults: influence of serotonergic and neurotrophic genes?

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.

Sex differences in sum scores may be hard to interpret: the importance of measurement invariance

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.

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

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.

The role of plasma cytokine levels, CRP and Selenoprotein S gene variation in OA

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.

Weighted statistics for aggregation and linkage analysis of human longevity in selected families: the Leiden Longevity Study

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.

Association of a nsSNP in ADAMTS14 to some osteoarthritis phenotypes

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.

Genetic influences on disordered eating behaviour are largely independent of body mass index

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.