Defining the relationship between lipids and immune traits in TwinsUK

Introduction

Lipid metabolism plays a crucial role in both innate and adaptive immune responses. Blood lipids are well known biomarkers of heart disease and atherosclerosis. Given that immune responses are involved in atherosclerosis development and progression, it is of utmost importance to explore the mutual relationship between circulating lipid measures and immunology markers, to better manage cardiovascular disease (CVD).

Purpose

The aim of this work was to evaluate relationships of distinct immune cell subtypes and circulating lipids in a general population setting.

Methods

We investigated the associations between 3,734 immune traits and circulating lipids [LDL-C, HDL-C, total cholesterol (TC), triglycerides (TG), non-HDL-C) in 2,056 individuals from the TwinsUK cohort. Immune traits were defined by high-dimensional flow-cytometry which enabled in-depth assessment of all major lymphocyte and myeloid subsets in peripheral blood. Information on immune cell subset frequency (FREQ) and mean fluorescence intensity (MFI) was available. A linear mixed model regression analysis adjusted for family structure, gender, age and BMI was used to examine the associations between immune traits and lipids. We then undertook a sub-analysis in 608 individuals and 131 lipid metabolites quantified using nuclear magnetic resonance spectroscopy to investigate the associations of specific immune traits and lipoprotein particles.

Results

We observed 18 significant associations (Bonferroni P-value <7.3x10–6) between 9 MFI immune traits and lipids. These involved associations between circulating lipids and 7 different myeloid or lymphoid subsets of immune cells including monocytes and natural killer cells. The top associations were found between CD14+ classical monocytes and non-HDL-C, LDL-C, HDL-C and TG. Increased levels of non-HDL-C (P=4.7x10–14), HDL-C (P=1.9x10–13) and TG (P=8.5x10–10) were associated with lower expression of CD14+ monocytes whereas increased LDL-C was associated with elevated expression of both CD14+ monocytes (P=1.0x10–13) and CD16+CD14+ monocytes (P=2.7x10–8). No significant association was detected between FREQ traits and lipids. We also observed 39 associations (Bonferroni P-value <3.1x10–4) between 19 lipid metabolites and CD14+ classical monocytes. The top associations were linked to different HDL particles (Free cholesterol in L-HDL: P=3.5x10–6; Total cholesterol in L-HDL: P=3.5x10–6; Cholesterol esters in L-HDL: P=4.1x10–6). Different subclasses of VLDL particles were also associated with CD14+ monocytes.

Conclusion

The present study provides the first comprehensive picture of lipid-immune cell associations in the circulatory system. We observed a link between expression levels of CD14+ classical monocytes and circulating lipid metabolites. The same lipid metabolites were previously associated with risk of developing myocardial infarction and stroke speculating a link between CD14+ classical monocytes and CVD risk.

Funding Acknowledgement

Type of funding sources: Public grant(s) – EU funding. Main funding source(s): SYSCID - A Systems medicine approach to chronic inflammatory diseases

The gut microbiome: what the oncologist ought to know

The gut microbiome (GM) has been implicated in a vast number of human pathologies and has become a focus of oncology research over the past 5 years. The normal gut microbiota imparts specific function in host nutrient metabolism, xenobiotic and drug metabolism, maintenance of structural integrity of the gut mucosal barrier, immunomodulation and protection against pathogens. Strong evidence is emerging to support the effects of the GM on the development of some malignancies but also on responses to cancer therapies, most notably, immune checkpoint inhibition. Tools for manipulating the GM including dietary modification, probiotics and faecal microbiota transfer (FMT) are in development. Current understandings of the many complex interrelationships between the GM, cancer, the immune system, nutrition and medication are ultimately based on a combination of short‐term clinical trials and observational studies, paired with an ever-evolving understanding of cancer biology. The next generation of personalised cancer therapies focusses on molecular and phenotypic heterogeneity, tumour evolution and immune status; it is distinctly possible that the GM will become an increasingly central focus amongst them. The aim of this review is to provide clinicians with an overview of microbiome science and our current understanding of the role the GM plays in cancer.

Diagnostic value of cutaneous manifestation of SARS-CoV-2 infection

Background

One of the challenging aspects of SARS‐CoV‐2 infection is its diverse multisystemic disease presentation.

Objectives

To evaluate the diagnostic value of cutaneous manifestations of SARS‐CoV‐2 infection and investigate their duration and timing in relation to other COVID‐19 symptoms.

Methods

We used data from 336 847 UK users of the COVID Symptom Study app to assess the diagnostic value of body rash or an acral rash in SARS‐CoV‐2 infection, and data from an independent online survey of 11 544 respondents to investigate skin‐specific symptoms and collect their photographs.

Results

Using data from the app, we show significant association between skin rashes and a positive swab test result (odds ratio 1·67, 95% confidence interval 1·42–1·97). Strikingly, among the respondents of the independent online survey, we found that 17% of SARS‐CoV‐2‐positive cases reported skin rashes as the first presentation, and 21% as the only clinical sign of COVID‐19. Together with the British Association of Dermatologists, we have compiled a catalogue of images of the most common skin manifestations of COVID‐19 from 400 individuals (https://covidskinsigns.com), which we have made publicly available to assist clinicians in recognition of this early clinical feature of COVID‐19.

Conclusions

Skin rashes cluster with other COVID‐19 symptoms, are predictive of a positive swab test, and occur in a significant number of cases, either alone or before other classical symptoms. Recognizing rashes is important in identifying new and earlier cases of COVID‐19.

What is already known about this topic?

Several studies conducted in hospital settings reported that patients with COVID‐19 presented with unusual skin rashes, including urticarial rashes, vesicular lesions and, less frequently, chilblains in fingers or toes.

What does this study add?

We confirmed, in a community‐based setting that also includes milder forms of the disease, that the presence of a skin rash is predictive of SARS‐CoV‐2 infection.

We provide a website with photos of skin manifestations to help healthcare professionals in diagnosing COVID‐19.

Skin rashes should be taken into account to provide a quick COVID‐19 diagnosis to curb the spread of the disease.

Linked Comment: Naldi. Br J Dermatol 2021; 184:793–794.

Optimal symptom combinations to aid COVID-19 case identification: Analysis from a community-based, prospective, observational cohort

Objectives

Diagnostic work-up following any COVID-19 associated symptom will lead to extensive testing, potentially overwhelming laboratory capacity whilst primarily yielding negative results. We aimed to identify optimal symptom combinations to capture most cases using fewer tests with implications for COVID-19 vaccine developers across different resource settings and public health.

Methods

UK and US users of the COVID-19 Symptom Study app who reported new-onset symptoms and an RT-PCR test within seven days of symptom onset were included. Sensitivity, specificity, and number of RT-PCR tests needed to identify one case (test per case [TPC]) were calculated for different symptom combinations. A multi-objective evolutionary algorithm was applied to generate combinations with optimal trade-offs between sensitivity and specificity.

Findings

UK and US cohorts included 122,305 (1,202 positives) and 3,162 (79 positive) individuals. Within three days of symptom onset, the COVID-19 specific symptom combination (cough, dyspnoea, fever, anosmia/ageusia) identified 69% of cases requiring 47 TPC. The combination with highest sensitivity (fatigue, anosmia/ageusia, cough, diarrhoea, headache, sore throat) identified 96% cases requiring 96 TPC.

Interpretation

We confirmed the significance of COVID-19 specific symptoms for triggering RT-PCR and identified additional symptom combinations with optimal trade-offs between sensitivity and specificity that maximize case capture given different resource settings.

Optimal symptom combinations to aid COVID-19 case identification: analysis from a community-based, prospective, observational cohort

OBJECTIVES

Diagnostic work-up following any COVID-19 associated symptom will lead to extensive testing, potentially overwhelming laboratory capacity whilst primarily yielding negative results. We aimed to identify optimal symptom combinations to capture most cases using fewer tests with implications for COVID-19 vaccine developers across different resource settings and public health.

METHODS

UK and US users of the COVID-19 Symptom Study app who reported new-onset symptoms and an RT-PCR test within seven days of symptom onset were included. Sensitivity, specificity, and number of RT-PCR tests needed to identify one case (test per case [TPC]) were calculated for different symptom combinations. A multi-objective evolutionary algorithm was applied to generate combinations with optimal trade-offs between sensitivity and specificity.

FINDINGS

UK and US cohorts included 122,305 (1,202 positives) and 3,162 (79 positive) individuals. Within three days of symptom onset, the COVID-19 specific symptom combination (cough, dyspnoea, fever, anosmia/ageusia) identified 69% of cases requiring 47 TPC. The combination with highest sensitivity (fatigue, anosmia/ageusia, cough, diarrhoea, headache, sore throat) identified 96% cases requiring 96 TPC.

INTERPRETATION

We confirmed the significance of COVID-19 specific symptoms for triggering RT-PCR and identified additional symptom combinations with optimal trade-offs between sensitivity and specificity that maximize case capture given different resource settings.

Pain in knee osteoarthritis is associated with variation in the neurokinin 1/substance P receptor (TACR1) gene

BACKGROUND

Substance P (SP) is a pain- and inflammation-related neuropeptide which preferentially binds to the neurokinin receptor 1 (NK₁ ). SP and NK₁ receptors have been implicated in joint pain, inflammation and damage in animal models and human studies of osteoarthritis (OA). The aim of this study was to test if genetic variation at the neurokinin 1 receptor gene (TACR1) is associated with pain in individuals with radiographic knee OA.

METHODS

Participants from the Genetics of OA and Lifestyle study were used for the discovery group (n = 1615). Genotype data for six SNPs selected to cover most variation in the TACR1 gene were used to test for an association with symptomatic OA. Replication analysis was performed using data from the Chingford 1000 Women Study, Hertfordshire Cohort Study, Tasmanian Older Adult Cohort Study and the Clearwater OA Study. In total, n = 1715 symptomatic OA and n = 735 asymptomatic OA individuals were analysed.

RESULTS

Out of six SNPs tested in the TACR1 gene, one (rs11688000) showed a nominally significant association with a decreased risk of symptomatic OA in the discovery cohort. This was then replicated in four additional cohorts. After adjusting for age, gender, body mass index and radiographic severity, the G (minor) allele at rs11688000 was associated with a decreased risk of symptomatic OA compared to asymptomatic OA cases (p = 9.90 × 10^-4 , OR = 0.79 95% 0.68-0.90 after meta-analysis).

CONCLUSIONS

This study supports a contribution from the TACR1 gene in human OA pain, supporting further investigation of this gene's function in OA.

SIGNIFICANCE

This study contributes to the knowledge of the genetics of painful osteoarthritis, a condition which affects millions of individuals worldwide. Specifically, a contribution from the TACR1 gene to modulating pain sensitivity in osteoarthritis is suggested.

Untangling the relationship between diet and visceral fat mass through blood metabolomics and gut microbiome profiling

BACKGROUND/OBJECTIVES

Higher visceral fat mass (VFM) is associated with an increased risk for developing cardio-metabolic diseases. The mechanisms by which an unhealthy diet pattern may influence visceral fat (VF) development has yet to be examined through cutting-edge multi-omic methods. Therefore, our objective was to examine the dietary influences on VFM and identify gut microbiome and metabolite profiles that link food intakes to VFM.

SUBJECTS/METHODS

In 2218 twins with VFM, food intake and metabolomics data available we identified food intakes most strongly associated with VFM in 50% of the sample, then constructed and tested the 'VFM diet score' in the remainder of the sample. Using linear regression (adjusted for covariates, including body mass index and total fat mass), we investigated associations between the VFM diet score, the blood metabolomics profile and the fecal microbiome (n=889), and confirmed these associations with VFM. We replicated top findings in monozygotic (MZ) twins discordant (⩾1 s.d. apart) for VFM, matched for age, sex and the baseline genetic sequence.

RESULTS

Four metabolites were associated with the VFM diet score and VFM: hippurate, alpha-hydroxyisovalerate, bilirubin (Z,Z) and butyrylcarnitine. We replicated associations between VFM and the diet score (beta (s.e.): 0.281 (0.091); P=0.002), butyrylcarnitine (0.199 (0.087); P=0.023) and hippurate (-0.297 (0.095); P=0.002) in VFM-discordant MZ twins. We identified a single species, Eubacterium dolichum to be associated with the VFM diet score (0.042 (0.011), P=8.47 × 10^-5), VFM (0.057 (0.019), P=2.73 × 10^-3) and hippurate (-0.075 (0.032), P=0.021). Moreover, higher blood hippurate was associated with elevated adipose tissue expression neuroglobin, with roles in cellular oxygen homeostasis (0.016 (0.004), P=9.82x10^-6).

CONCLUSIONS

We linked a dietary VFM score and VFM to E. dolichum and four metabolites in the blood. In particular, the relationship between hippurate, a metabolite derived from microbial metabolism of dietary polyphenols, and reduced VFM, the microbiome and increased adipose tissue expression of neuroglobin provides potential mechanistic insight into the influence of diet on VFM.

Gut microbiome diversity and high-fibre intake are related to lower long-term weight gain

Background:

Cross-sectional studies suggest that the microbes in the human gut have a role in obesity by influencing the human body’s ability to extract and store calories. The aim of this study was to assess if there is a correlation between change in body weight over time and gut microbiome composition.

Methods:

We analysed 16S ribosomal RNA gene sequence data derived from the faecal samples of 1632 healthy females from TwinsUK to investigate the association between gut microbiome measured cross-sectionally and longitudinal weight gain (adjusted for caloric intake and baseline body mass index). Dietary fibre intake was investigated as a possible modifier.

Results:

Less than half of the variation in long-term weight change was found to be heritable (h²=0.41 (0.31, 0.47)). Gut microbiota diversity was negatively associated with long-term weight gain, whereas it was positively correlated with fibre intake. Nine bacterial operational taxonomic units (OTUs) were significantly associated with weight gain after adjusting for covariates, family relatedness and multiple testing (false discovery rate <0.05). OTUs associated with lower long-term weight gain included those assigned to Ruminococcaceae (associated in mice with improved energy metabolism) and Lachnospiraceae. A Bacterioides species OTU was associated with increased risk of weight gain but this appears to be driven by its correlation with lower levels of diversity.

Conclusions:

High gut microbiome diversity, high-fibre intake and OTUs implicated in animal models of improved energy metabolism are all correlated with lower term weight gain in humans independently of calorie intake and other confounders.

Purine Arabinosides as Inhibitors of Human Haemopoietic Progenitor Cells

Four purine arabinosides that inhibit varicella-zoster virus (VZV) replication in vitro were tested as inhibitors of colony formation by progenitor cells from normal human bone marrow. In general, erythroid burst forming cells (BFU-E) were more sensitive to inhibition by these compounds than were either erythroid colony forming cells (CFU-E) or granulocyte/macrophage colony forming cells (CFU-GM). A 50% reduction in colony formation (IC50) was observed for BFU-E in the presence of 8 μM 6-methoxypurine arabinoside. Adenine arabinoside and hypoxanthine arabinoside had IC50 values of 1 μM and 4 μM respectively, whereas 6-ethoxypurine arabinoside was not inhibitory (IC50 > 50 μM). Enzyme studies showed that both 6-methoxypurine arabinoside and adenine arabinoside were converted to hypoxanthine arabinoside by adenosine deaminase. 6-Ethoxypurine arabinoside was a much less efficient substrate. When the BFU-E assays were performed in the presence of an inhibitor of adenosine deaminase, 6-methoxypurine arabinoside became non-inhibitory. In contrast, adenine arabinoside became much more inhibitory (IC50 = 0.03 μM). The potency of hypoxanthine arabinoside was unaffected. Thus, incubation of 6-methoxypurine arabinoside and adenine arabinoside under conditions appropriate for the BFU-E assay resulted in the in situ conversion of these compounds to hypoxanthine arabinoside. Biotransformation of compounds must be considered in the assessment of toxicity in vitro.

Quantitative genetics of circulating Dickkopf-related protein 1 (DKK1) in community-based sample of UK twins

Dickkopf-related protein 1 (DKK1) is a major inhibitor of Wnt signalling pathway but also plays an important role in bone formation. Its circulating levels appear to correlate significantly with plasma levels of inflammatory factors, fractalkine and IL-6. This study, using a large sample of UK twins, showed that the variation of each of these factors and correlation between them was explained by the genetic factors, and indicated possible association with DKK1 gene variants.

INTRODUCTION

DKK1 is involved in the development of several inflammatory conditions related to bone and joint degradation. Our objectives were to explore the genetic contribution (heritability) to circulating DKK1 variation and its correlation with other inflammatory cytokines, interleukin 6 (IL-6) and fractalkine, and to test whether the DKK1 heritability could be attributable to single nucleotide polymorphisms (SNPs) mapped to DKK1, IL-6 and FRCT genes.

METHODS

The study included a large community-based sample of 4939 women drawn from the general UK population. Plasma samples were analysed for circulating levels of DKK1, IL-6 and fractalkine (FRCT); 65 SNPs of DKK1, IL-6 and FRCT candidate genes, with MAF >0.1, were examined. We applied variance component analysis to evaluate contribution of putative genetic (including above SNPs) and environmental factors to variation of DKK1, and its correlation with IL-6 and FRCT.

RESULTS

Putative genetic factors explained 42.2 ± 2 % of the total variation of circulating DKK1 levels, and were also significant for fractalkine and IL-6 variations. Most importantly, we report significant phenotypic (0.208 ± 0.006-0.459 ± 0.007) and genetic (0.338 ± 0.069-0.617 ± 0.033) correlations between these molecules. We found evidence suggestive of association between the DKK1 and its structural genes variants.

CONCLUSIONS

Circulating DKK1 levels correlated significantly with levels of IL-6 and FRCT, known risk factors for several inflammatory processes suggesting a potential role of DKK1 in inflammation and tissue injury. Our results suggest the contribution of genetic factors in inter-individual variation of DKK1 levels in human population. However, further studies are required to determine genetic polymorphisms affecting DKK1 variation and its correlation with IL-6 and FRCT.

Genome-wide association study of lifetime cannabis use based on a large meta-analytic sample of 32 330 subjects from the International Cannabis Consortium

Cannabis is the most widely produced and consumed illicit psychoactive substance worldwide. Occasional cannabis use can progress to frequent use, abuse and dependence with all known adverse physical, psychological and social consequences. Individual differences in cannabis initiation are heritable (40-48%). The International Cannabis Consortium was established with the aim to identify genetic risk variants of cannabis use. We conducted a meta-analysis of genome-wide association data of 13 cohorts (N=32 330) and four replication samples (N=5627). In addition, we performed a gene-based test of association, estimated single-nucleotide polymorphism (SNP)-based heritability and explored the genetic correlation between lifetime cannabis use and cigarette use using LD score regression. No individual SNPs reached genome-wide significance. Nonetheless, gene-based tests identified four genes significantly associated with lifetime cannabis use: NCAM1, CADM2, SCOC and KCNT2. Previous studies reported associations of NCAM1 with cigarette smoking and other substance use, and those of CADM2 with body mass index, processing speed and autism disorders, which are phenotypes previously reported to be associated with cannabis use. Furthermore, we showed that, combined across the genome, all common SNPs explained 13-20% (P<0.001) of the liability of lifetime cannabis use. Finally, there was a strong genetic correlation (rg=0.83; P=1.85 × 10(-8)) between lifetime cannabis use and lifetime cigarette smoking implying that the SNP effect sizes of the two traits are highly correlated. This is the largest meta-analysis of cannabis GWA studies to date, revealing important new insights into the genetic pathways of lifetime cannabis use. Future functional studies should explore the impact of the identified genes on the biological mechanisms of cannabis use.

Plasma protein biomarkers of Alzheimer's disease endophenotypes in asymptomatic older twins: early cognitive decline and regional brain volumes

There is great interest in blood-based markers of Alzheimer's disease (AD), especially in its pre-symptomatic stages. Therefore, we aimed to identify plasma proteins whose levels associate with potential markers of pre-symptomatic AD. We also aimed to characterise confounding by genetics and the effect of genetics on blood proteins in general. Panel-based proteomics was performed using SOMAscan on plasma samples from TwinsUK subjects who are asymptomatic for AD, measuring the level of 1129 proteins. Protein levels were compared with 10-year change in CANTAB-paired associates learning (PAL; n = 195), and regional brain volumes (n = 34). Replication of proteins associated with regional brain volumes was performed in 254 individuals from the AddNeuroMed cohort. Across all the proteins measured, genetic factors were found to explain ~26% of the variability in blood protein levels on average. The plasma level of the mitogen-activated protein kinase (MAPK) MAPKAPK5 protein was found to positively associate with the 10-year change in CANTAB-PAL in both the individual and twin difference context. The plasma level of protein MAP2K4 was found to suggestively associate negatively (Q < 0.1) with the volume of the left entorhinal cortex. Future studies will be needed to assess the specificity of MAPKAPK5 and MAP2K4 to eventual conversion to AD.

Individuals with high bone mass have an increased prevalence of radiographic knee osteoarthritis

We previously reported an association between high bone mass (HBM) and a bone-forming phenotype of radiographic hip osteoarthritis (OA). As knee and hip OA have distinct risk factors, in this study we aimed to determine (i) whether HBM is also associated with knee OA, and (ii) whether the HBM knee OA phenotype demonstrates a similar pattern of radiographic features to that observed at the hip. HBM cases (defined by DXA BMD Z-scores) from the UK-based HBM study were compared with unaffected family controls and general population controls from the Chingford and Hertfordshire cohort studies. A single blinded observer graded AP weight-bearing knee radiographs for features of OA (Kellgren-Lawrence score, osteophytes, joint space narrowing (JSN), sclerosis) using an atlas. Analyses used logistic regression, adjusting a priori for age and gender, and additionally for BMI as a potential mediator of the HBM-OA association, using Stata v12. 609 HBM knees in 311 cases (mean age 60.8years, 74% female) and 1937 control knees in 991 controls (63.4years, 81% female) were analysed. The prevalence of radiographic knee OA, defined as Kellgren-Lawrence grade≥2, was increased in cases (31.5% vs. 20.9%), with age and gender adjusted OR [95% CI] 2.38 [1.81, 3.14], p<0.001. The association between HBM and osteophytosis was stronger than that for JSN, both before and after adjustment for BMI which attenuated the ORs for knee OA and osteophytes in cases vs. controls by approximately 50%. Our findings support a positive association between HBM and knee OA. This association was strongest for osteophytes, suggesting HBM confers a general predisposition to a subtype of OA characterised by increased bone formation.

Subclinical deformities of the hip are significant predictors of radiographic osteoarthritis and joint replacement in women. A 20 year longitudinal cohort study

OBJECTIVE

Femoroacetabular Impingement (FAI) and Acetabular Dysplasia are common deformities, which have been implicated as a major cause of hip osteoarthritis (OA). We examined whether these subtle deformities of the hip are associated with the development of radiographic OA and total hip replacement (THR) in women.

DESIGN

A population-based, longitudinal cohort of 1003 women underwent pelvis radiographs at years 2 and 20. Alpha Angle, Triangular Index Height, Lateral Centre Edge (LCE) angle and Extrusion Index were measured. An alpha angle of greater than 65° was defined as Cam-type FAI. Radiographic OA and the presence of a THR were then determined at 20 years.

RESULTS

Cam-type FAI was significantly associated with the development of radiographic OA. Each degree increase in alpha angle above 65° was associated with an increase in risk of 5% (Odds Ratio (OR) 1.05 [95% confidence interval (CI) 1.01-1.09]) for radiographic OA and 4% (OR 1.04 [95% CI 1.00-1.08]) for THR. For Acetabular Dysplasia, each degree reduction in LCE angle below 28° was associated with an increase in risk of 13.0% (OR 0.87 [95% CI 0.78-0.96]) for radiographic OA and 18% (OR 0.82 [95% CI 0.75-0.89]) for THR.

CONCLUSIONS

This study demonstrates that Cam-type FAI and mild Acetabular Dysplasia are predictive of subsequent OA and THR in a large female population cohort. These are independent of age, BMI and joint space and significantly improve current predictive models of hip OA development.

Large scale meta-analysis of urinary C-terminal telopeptide, serum cartilage oligomeric protein and matrix metalloprotease degraded type II collagen and their role in prevalence, incidence and progression of osteoarthritis

OBJECTIVE

To evaluate the role of three cartilage-derived biomarkers on osteoarthritis (OA): urinary C-terminal telopeptide (uCTX-II), serum cartilage oligomeric protein (sCOMP), and serum MMP degraded type II collagen (sC2M).

SUBJECTS AND METHODS

Samples from 3582 individuals from the Rotterdam Study, the Genetics osteoArthritis and Progression (GARP), the Chingford Study and the TwinsUK cohort were assayed using enzyme-linked immune sorbent assays. Log10 of concentration levels were correlated with risk of hip, hand and knee OA, hip and knee OA severity and incidence, and progression of knee OA, adjusting for age, gender and body mass index (BMI). Results were meta-analysed to assess overall significance.

RESULTS

After adjusting for covariates, sCOMP was associated with knee OA and hip and knee OA incidence. Furthermore, sC2M was associated with knee OA incidence and progression. After adjustment for multiple tests (Bonferroni P < 0.002) only the association between sCOMP and knee OA remained significant (odds ratio (OR) = 3.26 (95%CI 1.63-10.1) P = 0.0008 for each standard deviation (SD) increase in biomarker levels). Levels of uCTX-II were significantly associated with risk of hand, hip and knee OA, progression and incidence of knee OA. A receiver operating characteristics (ROC) analysis showed a consistent improvement in prediction of knee OA progression from an average area under the curve (AUC) is 0.646 for age, sex and BMI alone to an AUC = 0.668 including uCTX-II for prediction.

CONCLUSIONS

uCTX-II is the most informative biochemical marker for prediction of OA. Both sCOMP and C2M showed some association with OA, thus indicating that they are descriptive of disease activity.

The coupling of bone and cartilage turnover in osteoarthritis: opportunities for bone antiresorptives and anabolics as potential treatments?

Osteoarthritis (OA) is the most common form of arthritic disease, and a major cause of disability and impaired quality of life in the elderly. OA is a complex disease of the entire joint, affecting bone, cartilage and synovium that thereby presents multiple targets for treatment. This manuscript will summarise emerging observations from cell biology, preclinical and preliminary clinical trials that elucidate interactions between the bone and cartilage components in particular. Bone and cartilage health are tightly associated. Ample evidence has been found for bone changes during progression of OA including, but not limited to, increased turnover in the subchondral bone, undermineralisation of the trabecular structure, osteophyte formation, bone marrow lesions and sclerosis of the subchondral plate. Meanwhile, a range of investigations has shown positive effects on cartilage health when bone resorption is suppressed, or deterioration of the cartilage when resorption is increased. Known bone therapies, namely oestrogens, selective oestrogen receptor modifiers (SERMs), bisphosphonates, strontium ranelate, calcitonin and parathyroid hormone, might prove useful for treating two critical tissue components of the OA joint, the bone and the cartilage. An optimal treatment for OA likely targets at least these two tissue components. The patient subgroups for whom these therapies are most appropriate have yet to be fully defined but would likely include, at a minimum, those with high bone turnover.

Gradient Boosting as a SNP Filter: an Evaluation Using Simulated and Hair Morphology Data

Typically, genome-wide association studies consist of regressing the phenotype on each SNP separately using an additive genetic model. Although statistical models for recessive, dominant, SNP-SNP, or SNP-environment interactions exist, the testing burden makes an evaluation of all possible effects impractical for genome-wide data. We advocate a two-step approach where the first step consists of a filter that is sensitive to different types of SNP main and interactions effects. The aim is to substantially reduce the number of SNPs such that more specific modeling becomes feasible in a second step. We provide an evaluation of a statistical learning method called "gradient boosting machine" (GBM) that can be used as a filter. GBM does not require an a priori specification of a genetic model, and permits inclusion of large numbers of covariates. GBM can therefore be used to explore multiple GxE interactions, which would not be feasible within the parametric framework used in GWAS. We show in a simulation that GBM performs well even under conditions favorable to the standard additive regression model commonly used in GWAS, and is sensitive to the detection of interaction effects even if one of the interacting variables has a zero main effect. The latter would not be detected in GWAS. Our evaluation is accompanied by an analysis of empirical data concerning hair morphology. We estimate the phenotypic variance explained by increasing numbers of highest ranked SNPs, and show that it is sufficient to select 10K-20K SNPs in the first step of a two-step approach.

A study of pathophysiological factors associated with gastro-esophageal reflux disease in twins discordant for gastro-esophageal reflux symptoms

BACKGROUND

Differences in lower esophageal sphincter (LES) and peristaltic function and in transient LES relaxations (TLESR) have been described in patients with gastro-esophageal reflux disease (GERD). However, some of these differences may be the result of chronic GERD rather than being an underlying contributory factor.

METHODS

Twins discordant for GERD symptoms, i.e., only one twin had GERD symptoms, underwent standard LES and esophageal body manometry, and then using a sleeve sensor prolonged LES and pH monitoring, 30 min before and 60 min after a 250 mL 1200 kcal lipid meal.

KEY RESULTS

Eight monozygotic and 24 dizygotic female twins were studied. Although there was no difference in preprandial LES pressure (symptomatic 13.2 ± 7.1 mmHg vs asymptomatic 15.1 ± 6.2 mmHg, P = 0.4), LES pressure fell further postprandially in symptomatic twins (LES pressure area under the curve 465 ± 126 vs 331 ± 141 mmHg h, P < 0.01). 12/37 (32%) of acid reflux episodes in symptomatic twins occurred due to low LES pressure or deep inspiration/strain and 0/17 in asymptomatic twins (P = 0.01). There was no difference between symptomatic and asymptomatic twins in: peristaltic amplitude, ineffective esophageal body motility, hiatus hernia prevalence, or LES length. There was also no difference in TLESR frequency preprandially (symptomatic median 1(range 0-2) vs asymptomatic 0(0-2), P = 0.08) or postprandially (2.5(1-8) vs 3(1-6), P = 0.81).

CONCLUSIONS & INFERENCES

Twins with GERD symptoms had lower postprandial LES pressure and given the close genetic link between the twins, it is possible that such differences are caused by GERD. Acid reflux episodes associated with a hypotensive LES were seen in symptomatic, but not in asymptomatic twins.

Association between cartilage and bone biomarkers and incidence of radiographic knee osteoarthritis (RKOA) in UK females: a prospective study

OBJECTIVE

There is a need to find biochemical markers that would identify people with increased risk of developing radiographic knee osteoarthritis (RKOA). The aim of this study was to evaluate the ability of cartilage and bone biomarkers (cartilage oligomeric matrix protein (COMP), aggrecan, cellular inhibitor of apoptosis protein (cIAP), N-telopeptide-to-helix (NTx)) to predict RKOA incidence in a 10-year follow-up of UK females from the Chingford community study.

METHOD

Joint space narrowing (JSN), osteophytes (OSP) and Kellgren-Lawrence (K/L) grades were scored from radiographs of both knees at study baseline and 10 years later in 1,003 women aged 45-64. Circulating levels of biomarkers and demographic variables were measured at baseline. Statistical association analysis was conducted between the potential predictor factors measured at baseline and documentation of RKOA at 10-year follow-up.

RESULTS

Age and body mass index (BMI), were significant predictors of incidence of RKOA as assessed by K/L and OSP. Considering biomarkers, independent significant association was found between COMP circulating levels and K/L scores (Odd Ratio (OR) = 2.87, 95% Confidence Interval (CI) = 1.19-6.89, P = 0.018). Significant negative association was detected between aggrecan plasma concentrations and JSN, with OR = 0.37 (95% CI 0.15-0.89), P = 0.026.

CONCLUSIONS

Aggrecan and COMP circulating levels contribute to identification of phenotype-specific RKOA incidence. These data suggest potentially protective role of aggrecan in cartilage loss, as measured by JSN. High COMP levels are risk factors for development of RKOA, as assessed by K/L scores.

Neuropathic features of joint pain: a community-based study

OBJECTIVE

Quantitative sensory testing (QST) and questionnaire-based assessments have been used to demonstrate features of neuropathic pain in subjects with musculoskeletal pain. However, their direct relationship has not been investigated in the community. The purpose of this study was to conduct an observational study to describe the characteristics of joint pain and to examine the relationship between QST measures and the PainDETECT Questionnaire (PD-Q).

METHODS

Warm detection, heat pain, and mechanical pain thresholds as well as mechanical pain sensitivity over the sternum were determined and the PD-Q scores were calculated in a cross-sectional study of 462 participants in the Chingford Study. Comparisons were made between subjects with and those without joint pain. Logistic regression modeling was used to describe the association between neuropathic pain features, as determined by the PD-Q score, and each of the QST measures individually, adjusting for age, body mass index, and use of pain-modifying medications.

RESULTS

A total of 66.2% of the subjects reported recent joint pain, with a median average pain severity of 5 of 10. There was increased sensitivity to painful stimuli in the group with pain as compared to the pain-free group, and this persisted after stratification by pain-modifying medication use. While only 6.7% of subjects had possible neuropathic pain features and 1.9% likely neuropathic pain features according to the standard PD-Q thresholds, features of neuropathic pain were common and were present in >50% of those reporting pain of at least moderate severity. Heat pain thresholds and mechanical pain sensitivity were significantly associated with features of neuropathic pain identified using the PD-Q, with an odds ratio (OR) of 0.88 (95% confidence interval [95% CI] 0.79-0.97; P = 0.011) and an OR of 1.24 (95% CI 1.04-1.48; P = 0.018), respectively.

CONCLUSION

QST measures and the PD-Q identified features of neuropathic pain in subjects in this community-based study, with significant overlap between the findings of the two techniques.

Alcohol-induced metabolomic differences in humans

Alcohol consumption is one of the world's major risk factors for disease development. But underlying mechanisms by which moderate-to-heavy alcohol intake causes damage are poorly understood and biomarkers are sub-optimal. Here, we investigated metabolite concentration differences in relation to alcohol intake in 2090 individuals of the KORA F4 and replicated results in 261 KORA F3 and up to 629 females of the TwinsUK adult bioresource. Using logistic regression analysis adjusted for age, body mass index, smoking, high-density lipoproteins and triglycerides, we identified 40/18 significant metabolites in males/females with P-values <3.8E-04 (Bonferroni corrected) that differed in concentrations between moderate-to-heavy drinkers (MHD) and light drinkers (LD) in the KORA F4 study. We further identified specific profiles of the 10/5 metabolites in males/females that clearly separated LD from MHD in the KORA F4 cohort. For those metabolites, the respective area under the receiver operating characteristic curves were 0.812/0.679, respectively, thus providing moderate-to-high sensitivity and specificity for the discrimination of LD to MHD. A number of alcohol-related metabolites could be replicated in the KORA F3 and TwinsUK studies. Our data suggests that metabolomic profiles based on diacylphosphatidylcholines, lysophosphatidylcholines, ether lipids and sphingolipids form a new class of biomarkers for excess alcohol intake and have potential for future epidemiological and clinical studies.

A higher alkaline dietary load is associated with greater indexes of skeletal muscle mass in women

Conservation of muscle mass is important for fall and fracture prevention but further understanding of the causes of age-related muscle loss is required. This study found a more alkaline diet was positively associated with muscle mass in women suggesting a role for dietary acid-base load in muscle loss.

INTRODUCTION

Conservation of skeletal muscle is important for preventing falls and fractures but age-related loss of muscle mass occurs even in healthy individuals. However, the mild metabolic acidosis associated with an acidogenic dietary acid-base load could influence loss of muscle mass.

METHODS

We investigated the association between fat-free mass (FFM), percentage FFM (FFM%) and fat-free mass index (FFMI, weight/height²), measured using dual-energy X-ray absorptiometry in 2,689 women aged 18-79 years from the TwinsUK Study, and dietary acid-base load. Body composition was calculated according to quartile of potential renal acid load and adjusted for age, physical activity, misreporting and smoking habit (FFM, FFMI also for fat mass) and additionally with percentage protein.

RESULTS

Fat-free mass was positively associated with a more alkalinogenic dietary load (comparing quartile 1 vs 4: FFM 0.79 kg P < 0.001, FFM% 1.06 % <0.001, FFMI 0.24 kg/m² P = 0.002), and with the ratio of fruits and vegetables to potential acidogenic foods.

CONCLUSIONS

We observed a small but significant positive association between a more alkaline diet and muscle mass indexes in healthy women that was independent of age, physical activity and protein intake equating to a scale of effect between a fifth and one half of the observed relationship with 10 years of age. Although protein is important for maintenance of muscle mass, eating fruits and vegetables that supply adequate amounts of potassium and magnesium are also relevant. The results suggest a potential role for diet in the prevention of muscle loss.

GDF5 single-nucleotide polymorphism rs143383 is associated with lumbar disc degeneration in Northern European women

Objective

Lumbar disc degeneration (LDD) is a serious social and medical problem which has been shown to be highly heritable. It has similarities with peripheral joint osteoarthritis (OA) in terms of both epidemiology and pathologic processes. A few known genetic variants have been identified using a candidate gene approach, but many more are thought to exist. GDF5 is a gene whose variants have been shown to play a role in skeletal height as well as predisposing to peripheral joint OA. In vitro, the gene product growth differentiation factor 5 has been shown to promote growth and repair of animal disc. This study was undertaken to investigate whether the GDF5 gene plays a role in LDD.

Methods

We investigated whether the 5′ upstream single-nucleotide polymorphism (SNP) variant rs143383 was associated with LDD, using plain radiography and magnetic resonance imaging to identify disc space narrowing and osteophytes, in 5 population cohorts from Northern Europe.

Results

An association between LDD and the SNP rs143383 was identified in women, with the same risk allele as in knee and hip OA (odds ratio 1.72 [95% confidence interval 1.15–2.57], P = 0.008).

Conclusion

Our findings in 5 population cohorts from Northern Europe indicate that a variant in the GDF5 gene is a risk factor for LDD in women. Many more such variants are predicted to exist, but this result highlights the growth and differentiation cellular pathway as a possible route to a better understanding of the process behind lumbar disc degeneration.

Abstract OT3-3-01: Eniluracil + 5-fluorouracil + leucovorin (EFL) vs. capecitabine phase 2 trial for metastatic breast cancer

Based on a modified dosing protocol designed to optimize efficacy, an open-label EFL vs. capecitabine (4:3 randomization) Phase 2 trial for metastatic breast cancer is in progress. Eniluracil inactivates dihydropyrimidine dehydrogenase, thereby preventing the formation of α-fluoro-β-alanine, and conferring 100% oral bioavailability and a 5 hr half-life on 5-fluorouracil (5-FU). Study drugs are administered orally for 1st- or 2nd-line treatment for metastatic disease in patients previously treated with an anthracycline and a taxane. Arm 1: eniluracil (40 mg) taken 11–16 hr before 5-FU (30 mg/m2); leucovorin (30 mg) taken with 5-FU and the next day. The regimen is administered once/week for 3 weeks/4 weeks. Arm 2: capecitabine (1000 mg/m2) taken bid for 14 days/21days. Arm 2 patients with disease progression could crossover to take EFL in Arm X. Two sites in the USA and 19 in Russia are enrolling. Currently, 115 patients (21% are 1st-line, 70% had previous 5-FU treatment) are enrolled and 83 have had tumor assessments. EFL was well tolerated with no unexpected toxicities. As of May 2012, there were 11, 7, & 1 partial responses in Arms 1, 2, & X, respectively. The primary endpoint, progression-free survival, will be determined approximately 7.5 months after the trial is enrolled with 140 evaluable patients.

Citation Information: Cancer Res 2012;72(24 Suppl):Abstract nr OT3-3-01.

Common variants in genes encoding adiponectin (ADIPOQ) and its receptors (ADIPOR1/2), adiponectin concentrations, and diabetes incidence in the Diabetes Prevention Program

AIMS

Baseline adiponectin concentrations predict incident Type 2 diabetes mellitus in the Diabetes Prevention Program. We tested the hypothesis that common variants in the genes encoding adiponectin (ADIPOQ) and its receptors (ADIPOR1, ADIPOR2) would associate with circulating adiponectin concentrations and/or with diabetes incidence in the Diabetes Prevention Program population.

METHODS

Seventy-seven tagging single-nucleotide polymorphisms (SNPs) in ADIPOQ (24), ADIPOR1 (22) and ADIPOR2 (31) were genotyped. Associations of SNPs with baseline adiponectin concentrations were evaluated using linear modelling. Associations of SNPs with diabetes incidence were evaluated using Cox proportional hazards modelling.

RESULTS

Thirteen of 24 ADIPOQ SNPs were significantly associated with baseline adiponectin concentrations. Multivariable analysis including these 13 SNPs revealed strong independent contributions of rs17366568, rs1648707, rs17373414 and rs1403696 with adiponectin concentrations. However, no ADIPOQ SNPs were directly associated with diabetes incidence. Two ADIPOR1 SNPs (rs1342387 and rs12733285) were associated with ∼18% increased diabetes incidence for carriers of the minor allele without differences across treatment groups, and without any relationship with adiponectin concentrations.

CONCLUSIONS

  ADIPOQ SNPs are significantly associated with adiponectin concentrations in the Diabetes Prevention Program cohort. This observation extends prior observations from unselected populations of European descent into a broader multi-ethnic population, and confirms the relevance of these variants in an obese/dysglycaemic population. Despite the robust relationship between adiponectin concentrations and diabetes risk in this cohort, variants in ADIPOQ that relate to adiponectin concentrations do not relate to diabetes risk in this population. ADIPOR1 variants exerted significant effects on diabetes risk distinct from any effect of adiponectin concentrations.

Genome-wide association analysis of coffee drinking suggests association with CYP1A1/CYP1A2 and NRCAM

Coffee consumption is a model for addictive behavior. We performed a meta-analysis of genome-wide association studies (GWASs) on coffee intake from 8 Caucasian cohorts (N=18 176) and sought replication of our top findings in a further 7929 individuals. We also performed a gene expression analysis treating different cell lines with caffeine. Genome-wide significant association was observed for two single-nucleotide polymorphisms (SNPs) in the 15q24 region. The two SNPs rs2470893 and rs2472297 (P-values=1.6 × 10(-11) and 2.7 × 10(-11)), which were also in strong linkage disequilibrium (r(2)=0.7) with each other, lie in the 23-kb long commonly shared 5' flanking region between CYP1A1 and CYP1A2 genes. CYP1A1 was found to be downregulated in lymphoblastoid cell lines treated with caffeine. CYP1A1 is known to metabolize polycyclic aromatic hydrocarbons, which are important constituents of coffee, whereas CYP1A2 is involved in the primary metabolism of caffeine. Significant evidence of association was also detected at rs382140 (P-value=3.9 × 10(-09)) near NRCAM-a gene implicated in vulnerability to addiction, and at another independent hit rs6495122 (P-value=7.1 × 10(-09))-an SNP associated with blood pressure-in the 15q24 region near the gene ULK3, in the meta-analysis of discovery and replication cohorts. Our results from GWASs and expression analysis also strongly implicate CAB39L in coffee drinking. Pathway analysis of differentially expressed genes revealed significantly enriched ubiquitin proteasome (P-value=2.2 × 10(-05)) and Parkinson's disease pathways (P-value=3.6 × 10(-05)).

The natural history of radiographic knee osteoarthritis: a fourteen-year population-based cohort study

OBJECTIVE

To establish the natural history of radiographic knee osteoarthritis (OA) over 14 years in a community-based cohort.

METHODS

We examined women from the Chingford Women's Study, a community-based cohort followed up for more than 14 years. We selected women for whom bilateral radiographs of the knees (with the legs in full extension) were obtained at approximately 5-year intervals. Radiographs were scored for OA in a blinded manner, using Kellgren/Lawrence (K/L) grades. Descriptive statistics and odds ratios (ORs) were used to compare the incidence, worsening, and progression of radiographic knee OA.

RESULTS

A complete radiography series was available for 561 of the original 1,003 subjects enrolled in the study. The median age of these subjects at baseline was 53 years (interquartile range 48-58 years). At baseline, 13.7% of the subjects had radiographic knee OA (K/L grade≥2) in at least one knee, and the prevalence increased to 47.8% by year 15. The annual cumulative incidence of radiographic knee OA was 2.3% between baseline and year 15. The annual rates of disease progression and worsening between baseline and year 15 were 2.8% and 3.0%, respectively. Subjects with a K/L grade of 1 at baseline were more likely to experience worsening by year 15 compared with subjects with a baseline grade of 0 (OR 4.5, 95% confidence interval 2.7-7.4).

CONCLUSION

This is the longest natural history study of radiographic knee OA to date. The results showed relatively low rates for the incidence and progression of radiographic knee OA; more than half of all subjects had no radiographic evidence of knee OA over a 15-year period of time. Subjects with a baseline K/L grade of 1 were more likely than subjects with other baseline K/L grades to experience worsening of knee OA.

A twin study of body dysmorphic concerns

BACKGROUND

Dysmorphic concern refers to an excessive preoccupation with a perceived or slight defect in physical appearance. It lies on a continuum of severity from no or minimal concerns to severe concerns over one's appearance. The present study examined the heritability of dysmorphic concerns in a large sample of twins.

METHOD

Twins from the St Thomas UK twin registry completed a valid and reliable self-report measure of dysmorphic concerns, which also includes questions about perceived body odour and malfunction. Twin modelling methods (female twins only, n=3544) were employed to decompose the variance in the liability to dysmorphic concerns into additive genetic, shared and non-shared environmental factors.

RESULTS

Model-fitting analyses showed that genetic factors accounted for approximately 44% [95% confidence intervals (CI) 36-50%] of the variance in dysmorphic concerns, with non-shared environmental factors and measurement error accounting for the remaining variance (56%; 95% CI 50-63%). Shared environmental factors were negligible. The results remained unchanged when excluding individuals reporting an objective medical condition/injury accounting for their concern in physical appearance.

CONCLUSIONS

Over-concern with a perceived or slight defect in physical appearance is a heritable trait, with non-shared environmental factors also playing an important role in its causation. The results are relevant for various psychiatric disorders characterized by excessive concerns in body appearance, odour or function, including but not limited to body dysmorphic disorder.

Prevalence of reported knee pain over twelve years in a community-based cohort

OBJECTIVE

To describe the temporal patterns of knee pain in a community-based cohort over 12 years.

METHODS

Data on self-reported knee pain at 4 time points over 12 years were analyzed in participants from the Chingford Women's Study of osteoarthritis (OA) and osteoporosis. Pain status was defined as any pain in the preceding month and pain on most days in the preceding month. This status was used to classify participants according to pain patterns of asymptomatic, persistent, incident, or intermittent pain. Multinomial logistic regression was used to identify baseline predictors for each pain pattern.

RESULTS

Among the 489 women with complete followup data, the median age at baseline was 52 years (interquartile range [IQR] 48-58 years), the median body mass index (BMI) was 24.39 kg/m(2) (IQR 22.46-27.20), and 11.7% of the women had a Kellgren/Lawrence radiographic OA severity grade of ≥2 in at least one knee. Among subjects reporting any pain in the preceding month versus those reporting pain on most days in the preceding month, 9% versus 2% had persistent pain, 24% versus 16% had incident pain, and 29% versus 18% had intermittent pain. A higher BMI was predictive of persistent pain (odds ratio [OR] 1.14, 95% confidence interval [95% CI] 1.04-1.25) and incident pain (OR 1.10, 95% CI 1.02-1.18). The presence of radiographic knee OA was predictive of persistent pain (OR 3.70, 95% CI 1.34-10.28; P = 0.012), and reported knee injury was predictive of both persistent pain (OR 4.13, 95% CI 1.34-12.66; P = 0.013) and intermittent pain (OR 4.25, 95% CI 1.81-9.98; P = 0.001).

CONCLUSION

Significant variability in the temporal fluctuation of self-reported knee pain was seen in this community-based prospective study over a period of 12 years, with few women consistently reporting knee pain at each time point. Distinct baseline predictors for each pain pattern were identified and may explain the observed heterogeneity of self-reported knee pain when pain status is measured at only one time point.

Genetic and environmental risk factors for sexual distress and its association with female sexual dysfunction

BACKGROUND

The DSM-V Working Group is currently re-evaluating distress as a primary diagnostic criterion for female sexual dysfunction (FSD). Here, for the first time, we explored the epidemiology of sexual distress and its putative aetiological relationship to FSD by estimating the influence of genetic and environmental risk factors.

METHOD

Questionnaire data on a representative sample of 930 British female twins using validated scales of FSD and sexual distress were subject to variance components analyses to quantify latent genetic and environmental factors influencing phenotypic variation and covariation. Multiple regression analyses were used to identify other potential risk factors of sexual distress.

RESULTS

Of 319 women with any sexual problems, only 36.5% reported distress. Of women classified as functional, 16.5% felt sexual distress. Sexual distress had a heritability of 44% [95% confidence interval (CI) 0.33-0.54]. Bivariate analysis suggested that the majority (91% CI 86-99%) of the covariance between sexual distress and FSD was due to unique environmental effects common to both traits. Associations were found between sexual distress and other risk variables, including relationship dissatisfaction [odds ratio (OR) 1.6, p<0.001], anxiety sensitivity and obsessive-compulsive symptomatology (OR 1.2, p<0.01, for both).

CONCLUSIONS

There seems to be a weak phenotypic and genetic basis for including sexual distress as a diagnostic indicator of FSD. Instead, the data indicate that unrelated psychological factors play an important role in sexual distress and tentatively suggest that sexual distress is less a consequence of FSD and more related to general anxiety among women.

Does obesity predict knee pain over fourteen years in women, independently of radiographic changes?

OBJECTIVE

To examine longitudinal patterns in body mass index (BMI) over 14 years and its association with knee pain in the Chingford Study.

METHODS

We studied a total of 594 women with BMI data from clinic visits at years (Y) 1, 5, 10, and 15. Knee pain at Y15 was assessed by questionnaire. Associations between BMI over 14 years and knee pain at Y15 were examined using logistic regression.

RESULTS

BMI significantly increased from Y1 to Y15 (P < 0.0005) with medians (interquartile ranges) of 24.5 kg/m(2) (22.5-27.2 kg/m(2) ) and 26.5 kg/m(2) (23.9-30.1 kg/m(2) ), respectively. At Y15, 45.1% of subjects had knee pain. A greater BMI at Y1 (odds ratio [OR] 1.34, 95% confidence interval [95% CI] 1.05-1.69), at Y15 (OR 1.34, 95% CI 1.10-1.61), and change in BMI over 15 years (OR 1.40, 95% CI 1.00-1.93) were significant predictors of knee pain at Y15 (P < 0.05). BMI change was associated with bilateral (OR 1.61, 95% CI 1.05-1.76, P = 0.024) but not unilateral knee pain (OR 1.22, 95% CI 0.73-1.76, P = 0.298). The association between BMI change and knee pain was independent of radiographic knee osteoarthritis (OA). The strength of association between BMI and knee pain at Y15 was similar during followup measurements.

CONCLUSION

Over 14 years, a higher BMI predicts knee pain at Y15 in women, independently of radiographic knee OA. When adjusted, the association was significant in bilateral, not unilateral, knee pain, suggesting alternative pathologic mechanisms may exist. The longitudinal effect of BMI on knee pain at Y15 is equally important at any time point, which may assist reducing the population burden of knee pain.

Homozygous deletion of the UGT2B17 gene is not associated with osteoporosis risk in elderly Caucasian women

Previously, homozygous deletion of the UGT2B17 gene has shown association with hip fracture. Using a high-throughput qRT-PCR assay, we genotyped UGT2B17 copy number variation (CNV) in 1,347 elderly Caucasian women and examined for effects on bone phenotypes. We found no evidence of association between UGT2B17 CNV and osteoporosis risk in this population.

INTRODUCTION

Genetic studies of osteoporosis commonly examine SNPs in candidate genes or whole genome analyses, but insertions and deletions of DNA, collectively called CNV, also comprise a large amount of the genetic variability between individuals. Previously, homozygous deletion of the UGT2B17 gene in CNV 4q13.2, which encodes an enzyme that mediates the glucuronidation of steroid hormones, has shown association with the risk of hip fracture.

METHODS

We used a quantitative real-time PCR assay for genotyping the UGT2B17 CNV in a well-characterized population study of 1,347 Caucasian women aged 75.2 ± 2.7 years (mean ± SD), to assess the effect of the CNV on bone mass density (BMD) at the total hip site and osteoporosis risk.

RESULTS

The UGT2B17 CNV distribution was consistent with the expected Hardy-Weinberg distribution and not different from frequencies previously reported in a Caucasian population. Data from ANCOVA of age- and weight-adjusted BMD for UGT2B17 CNV genotype showed no significant difference between genotype groups. Individuals with homozygous or heterozygous deletion of the UGT2B17 gene showed no increased risk of incident fragility fracture.

CONCLUSIONS

These data suggest that quantitative real-time PCR is a rapid and efficient technique for determination of candidate CNVs, including the UGT2B17 CNV; however, we found no evidence of an effect of UGT2B17 CNV on osteoporosis risk in elderly Caucasian women.

Large-scale meta-analysis of interleukin-1 beta and interleukin-1 receptor antagonist polymorphisms on risk of radiographic hip and knee osteoarthritis and severity of knee osteoarthritis

OBJECTIVE

To clarify the role of common genetic variation in the Interleukin-1β (IL1B) and Interleukin-1R antagonist (IL1RN) genes on risk of knee and hip osteoarthritis (OA) and severity of knee OA by means of large-scale meta-analyses.

METHODS

We searched PubMed for articles assessing the role of IL1B and IL1RN polymorphisms/haplotypes on the risk of hip and/or knee OA. Novel data were included from eight unpublished studies. Meta-analyses were performed using fixed- and random-effects models with a total of 3595 hip OA and 5013 knee OA cases, and 6559 and 9132 controls respectively. The role of ILRN haplotypes on radiographic severity of knee OA was tested in 1918 cases with Kellgren-Lawrence (K/L) 1 or 2 compared to 199 cases with K/L 3 or 4.

RESULTS

The meta-analysis of six published studies retrieved from the literature search and eight unpublished studies showed no evidence of association between common genetic variation in the IL1B or IL1RN genes and risk of hip OA or knee OA (P>0.05 for rs16944, rs1143634, rs419598 and haplotype C-G-C (rs1143634, rs16944 and rs419598) previously implicated in risk of hip OA). The C-T-A haplotype formed by rs419598, rs315952 and rs9005, previously implicated in radiographic severity of knee OA, was associated with reduced severity of knee OA (odds ratio (OR)=0.71 95%CI 0.56-0.91; P=0.006, I(2)=74%), and achieved borderline statistical significance in a random-effects model (OR=0.61 95%CI 0.35-1.06 P=0.08).

CONCLUSION

Common genetic variation in the Interleukin-1 region is not associated with prevalence of hip or knee OA but our data suggest that IL1RN might have a role in severity of knee OA.

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.

Significant differences in UK and US female bone density reference ranges

In the United Kingdom (UK), T- and Z-scores are usually calculated using reference ranges derived from United States (US) populations. In the UK arm of a recent randomised trial (International Breast Cancer Intervention Study II (IBIS-II)), substantially, fewer women than expected were recruited into the osteopenic (-2.5<T-score<-1.0) and osteoporotic (T-score<-2.5) arms of the study. The comparison with data from two independent studies showed that UK women aged >45 years with a typical body mass index of 28 kg m(-2) have spine and hip bone mineral density (BMD) 0.6 standard deviation higher than their US counterparts.

INTRODUCTION

Dual energy X-ray absorptiometry (DXA) is widely used for the diagnosis of osteoporosis and to investigate the effect of pharmacological treatments on BMD. In both routine and research settings, it is important that DXA results are correctly interpreted.

METHODS

T- and Z-scores for the first 650 UK Caucasian women enrolled in the IBIS-II study were compared with data from two independent studies of unrelated, unselected UK Caucasian women: (1) 2,382 women aged 18 to 79 recruited to the Twins UK Adult Twin Registry; (2) 431 women aged 21 to 84 with no risk factors for osteoporosis recruited at Guy's Hospital. All DXA measurements were performed on Hologic densitometers. Subjects were divided into six age bands, and T- and Z-scores were calculated using the manufacturer's US reference range for the spine and the National Health and Nutrition Examination Survey III reference range for the femoral neck and total hip.

RESULTS

The overall mean Z-scores for the IBIS-II, Twin, and Guy's groups were: spine: +0.61, +0.29, +0.33; femoral neck: +0.42, +0.36, +0.45; total hip: +0.65, +0.38, +0.39 (all p<0.001 compared with the expected value of 0). The mean body weight of subjects in the three studies was 74.4, 65.5, and 65.4 kg, respectively. Analysis revealed a highly significant relationship between Z-score and weight at each BMD site with a slope of 0.03 kg(-1).

CONCLUSIONS

In general, US spine and hip reference ranges are not suitable for the calculation of Z-scores in UK women. For some research study designs, the differences may significantly influence the pattern of subject recruitment.