Understanding the complex genetic architecture connecting rheumatoid arthritis, osteoporosis and inflammation: discovering causal pathways.

Rheumatoid arthritis (RA) and osteoporosis (OP) are two comorbid complex inflammatory conditions with evidence of shared genetic background and causal relationships. We aimed to clarify the genetic architecture underlying RA and various OP phenotypes while additionally considering an inflammatory component, C-reactive protein (CRP). Genome-wide association study summary statistics were acquired from the GEnetic Factors for OSteoporosis Consortium, Cohorts for Heart and Aging Research Consortium and UK Biobank. Mendelian randomization (MR) was used to detect the presence of causal relationships. Colocalization analysis was performed to determine shared genetic variants between CRP and OP phenotypes. Analysis of pleiotropy between traits owing to shared causal single nucleotide polymorphisms (SNPs) was performed using PL eiotropic A nalysis under CO mposite null hypothesis (PLACO). MR analysis was suggestive of horizontal pleiotropy between RA and OP traits. RA was a significant causal risk factor for CRP (ß = 0.027, 95% confidence interval = 0.016-0.038). There was no evidence of CRP→OP causal relationship, but horizontal pleiotropy was apparent. Colocalization established shared genomic regions between CRP and OP, including GCKR and SERPINA1 genes. Pleiotropy arising from shared causal SNPs revealed through the colocalization analysis was all confirmed by PLACO. These genes were found to be involved in the same molecular function 'protein binding' (GO:0005515) associated with RA, OP and CRP. We identified three major components explaining the epidemiological relationship among RA, OP and inflammation: (1) Pleiotropy explains a portion of the shared genetic relationship between RA and OP, albeit polygenically; (2) RA contributes to CRP elevation and (3) CRP, which is influenced by RA, demonstrated pleiotropy with OP.

Shared genetic influence on frailty and chronic widespread pain: a study from TwinsUK.

Introduction: frailty is an increased vulnerability to adverse health outcomes, across multiple physiological systems, with both environmental and genetic drivers. The two most commonly used measures are Rockwood's frailty index (FI) and Fried's frailty phenotype (FP). Material and methods: the present study included 3626 individuals from the TwinsUK Adult Twin Registry. We used the classical twin model to determine whether FI and FP share the same latent aetiological factors. We also investigated the relationship between frailty and chronic widespread musculoskeletal pain (CWP), another holistic age-related condition with significant clinical impact. Results: FP and FI shared underlying genetic and environmental aetiology. CWP was associated with both frailty measures, and health deficits appeared to mediate the relationship between phenotypic frailty and pain. Latent genetic factors underpinning CWP were shared with frailty. While frailty was increased in the twins reporting pain, co-twin regression analysis indicated that the relationship between CWP and frailty is reduced after accounting for shared genetic and environmental factors. Conclusions: both measures of frailty tap the same root causes, thus this work helps unify frailty research. We confirmed a strong association between CWP and frailty, and showed a large and significant shared genetic aetiology of both phenomena. Our findings argue against pain being a significant causative factor in the development of frailty, favouring common causation. This study highlights the need to manage CWP in frail individuals and undertake a Comprehensive Geriatric Assessment in individuals presenting with CWP. Finally, the search for genetic factors underpinning CWP and frailty could be aided by integrating measures of pain and frailty.

Environmental rather than genetic factors determine the variation in the age of the infancy to childhood transition: a twins study.

OBJECTIVE: Using a twins study, we sought to assess the contribution of genetic against environmental factor as they affect the age at transition from infancy to childhood (ICT). STUDY DESIGN: The subjects were 56 pairs of monozygotic twins, 106 pairs of dizygotic twins, and 106 pairs of regular siblings (SBs), for a total of 536 children. Their ICT was determined, and a variance component analysis was implemented to estimate components of the familial variance, with simultaneous adjustment for potential covariates. RESULTS: We found substantial contribution of the common environment shared by all types of SBs that explained 27.7% of the total variance in ICT, whereas the common twin environment explained 9.2% of the variance, gestational age 3.5%, and birth weight 1.8%. In addition, 8.7% was attributable to sex difference, but we found no detectable contribution of genetic factors to inter-individual variation in ICT age. CONCLUSIONS: Developmental plasticity impacts much of human growth. Here we show that of the ∼50% of the variance provided to adult height by the ICT, 42.2% is attributable to adaptive cues represented by shared twin and SB environment, with no detectable genetic involvement.

Low back and common widespread pain share common genetic determinants.

Low back (LBP) and chronic widespread musculoskeletal pain (CWP) both have a significant genetic component and are associated with increased body mass index (BMI). We examined whether LBP and CWP share common genetic factors, and to what extent this correlation is modified by the genetic factors influencing BMI. Genetic analysis of binary traits such as pain is not simple, particularly if their risk is associated with age or other quantitative traits. Implementing Falconer's polygenic threshold concept for dichotomous traits inheritance, we developed new software to examine the extent of the genetic influence on LBP and CWP under age and BMI dependence. The analysis was conducted on 3266 and 2256 UK female twins, assessed for LBP and CWP, respectively. Analysis of the liability scores with threshold to LBP and CWP established substantial contribution of genetic factors to their variation (h(2) > 0.60, p<0.004-0.0003) and covariation (p=3.1E-08). Some 39% of the CWP and 70% of the LBP heritability estimates were attributable to genetic effects shared by both phenotypes, and 40% and 67% of the residual variation is caused by environmental factors simultaneously affecting both pain syndromes. However, contribution of BMI to variation/covariation of both pain phenotypes-although statistically highly significant (p∼10-7)-was not determinative.

Association of interleukin-6 gene polymorphisms with hand osteoarthritis and hand osteoporosis.

OBJECTIVE: Several genes, including IL-6 encoding pro-inflammatory cytokines, are involved in development of osteoarthritis and osteoporosis. The association of radiographic hand osteoarthritis (RHOA) and osteoporosis related phenotypes (RHOP) with polymorphisms in IL-6 has been reported inconsistently. The aim of this study was to examine the association, between RHOA and RHOP and IL-6 polymorphisms in two independent samples. METHODS: Two samples: UK females, including 1440 individuals assessed for RHOA and 3470 assessed for RHOP; Chuvash pedigree including 1499 females and males were assessed for RHOP and RHOA. SNPs were genotyped in the IL-6 genomic region, and used in association analysis with RHOA and RHOP phenotypes. RESULTS: RHOP phenotypes showed similar heritability estimates in both samples, ranging from 34.5 ± 5.5% to 61.0 ± 2.4%. RHOA in Chuvash had substantially lower heritability estimates compared to twins (e.g. OSP scores: 11.8 ± 2.3% vs. 39.2 ± 4.1%) with much higher prevalence and considerably stronger correlation with age (r = 0.811 vs. r = 0.505). RHOA in Chuvash sample may be traumatic in nature, caused by heavy and prolonged manual work related to their private farming. There were a number of statistically significant association results with both types of phenotypes. The most consistent result was obtained for JSN in both samples with SNP from the same haploblock. Their combined probability of no association was only p = 0.000003. Additionally, there were SNPs common for both RHOA and RHOP. CONCLUSIONS: We have shown polymorphisms in IL_6 are significantly associated with RHOA and hand RHOP in two samples having different ethnicity and lifestyle. Age × environment × genes interaction appears as an important factor of RHOA manifestation and progression.

Inheritance of quantitative dermatoglyphic traits with asymmetry and diversity in Muzeina Bedouin tribe: a small isolated and consanguineous population from South Sinai.

The genetic factors contribute significantly to the determination of dermatoglyphic traits is well established. However, the controversies in views and findings of this issue are still inconclusive. The present study is an attempt to evaluate the inheritance of quantitative dermatoglyphic traits with asymmetry (DA and FA) and diversity (Div) through sibling correlations. Data include 218 individuals from (88 families) in a small isolate, the nomadic tribe Muzeina with a high degree of consanguinity (0.09) from South Sinai. Statistical analyses include sibling correlations, cross-correlations and genetic correlation (GC)--a ratio of sibling cross-correlation between traits divided on square root of the both traits sibling correlation product. The familial correlation coefficients for quantitative dermatoglyphic traits are perhaps expected lower in such a small isolated and consanguineous population than our previous studied in Indian populations and Chuvashian populations from Russia. These results indicate a simpler genetic basis due to high degree (0.09 inbreeding coefficient) of consanguinity in Muzeina Bedouin tribe. There is no evidence of major gene involvement, although a little genetic effect obtained from familial correlations on asymmetry (DA and FA) and diversity (Div) traits through sibling correlations. The significant interaction between sexes was found, which contradicts with the other populations perhaps due to high level of consanguinity. Lower correlation coefficients than in other non-consanguineous populations for quantitative dermatoglyphic traits indicate a simpler genetic basis due to high degree of inbreeding coefficient (0.09) in Muzeina. Dermatoglyphic asymmetry and diversity traits may be due to environmental factors rather than dominance in Bedouins, although a little genetic effect was found suggests a measure of developmental instability in human (FA).

Inheritance of dermatoglyphic asymmetry and diversity traits in twins based on factor: variance decomposition analysis.

Dermatoglyphic asymmetry and diversity traits from a large number of twins (MZ and DZ) were analyzed based on principal factors to evaluate genetic effects and common familial environmental influences on twin data by the use of maximum likelihood-based Variance decomposition analysis. Sample consists of monozygotic (MZ) twins of two sexes (102 male pairs and 138 female pairs) and 120 pairs of dizygotic (DZ) female twins. All asymmetry (DA and FA) and diversity of dermatoglyphic traits were clearly separated into factors. These are perfectly corroborated with the earlier studies in different ethnic populations, which indicate a common biological validity perhaps exists of the underlying component structures of dermatoglyphic characters. Our heritability result in twins clearly showed that DA_F2 is inherited mostly in dominant type (28.0%) and FA_F1 is additive (60.7%), but no significant difference in sexes was observed for these factors. Inheritance is also very prominent in diversity Factor 1, which is exactly corroborated with our previous findings. The present results are similar with the earlier results of finger ridge count diversity in twin data, which suggested that finger ridge count diversity is under genetic control.

Insights into the pleiotropic relationships between chronic back pain and inflammation-related musculoskeletal conditions: rheumatoid arthritis and osteoporotic abnormalities.

ABSTRACT: The ageing process includes the development of debilitating musculoskeletal (MSK) conditions, including chronic back pain (CBP), rheumatoid arthritis (RA), and osteoporosis (OP). The mechanisms involved in the genetic-epidemiological relationships between these MSK phenotypes are controversial and limited and thus require clarification, in particular, between CBP and the other MSK phenotypes. A cross-sectional statistical analysis was conducted using Europeans from the UK Biobank data collection, including 73,794 CBP, 4883 RA, and 7153 OP cases as well as 242,216 calcaneus bone mineral density scores. C-reactive protein (CRP) was measured for 402,165 subjects in this sample. Genetic correlations were assessed to evaluate shared genetic background between traits. Mendelian randomization was performed to assess a causal relationship between CBP and RA and OP along with other risk factors, such as CRP. Colocalization analysis was conducted to identify shared pleiotropic regions between the examined traits. Bayesian modelling was performed to determine a potential pathway that may explain the interrelationships among these traits. Mendelian randomization analyses revealed that CRP causally predicts CBP only (ß = 0.183, 95% CI = 0.077-0.290, P -value = 0.001). Horizontally pleiotropy appeared to explain the relationship between CBP and RA and OP. Through colocalization analysis, several genomic regions emerged describing common genetic influences between CBP and its proposed risk factors, including HLA-DQA1/HLA-DQB1, APOE , SOX5, and MYH7B as well as Histone 1 genes. We speculate that among other factors, CBP and its MSK comorbidities may arise from common inflammatory mechanisms. Colocalized identified genes may aid in advancing or improving the mode of treatment in patients with CBP.

Association of FTO gene variants with body composition in UK twins.

The association of FTO gene variants with body mass index (BMI) and other obesity characteristics is well established. However, uncertainties remain whether the association is present only in young populations and whether it is attributable to body fat mass specifically. We aimed to clarify these two questions in a large sample (N= 4,523 individuals) of middle-aged and older (range 40-80 years) British female twins. The women were assessed for BMI, waist and hip circumference, total lean (LBM) and fat (FBM) body mass. Since the majority of FTO association signals have been reported in a haploblock bordering 52,355-52,408 kb (on chromosome 16q12.2), we examined five genotyped and 43 imputed SNPs mapped to this block. Canonical correlation and other association analyses showed significant and consistent association between the selected SNP and studied body composition phenotypes, with p-values reaching p= 0.000004. Of particular interest, in addition to the expected significant associations between FTO variants and FBM, we also identified significant associations with LBM. These results suggest that the association between FTO variants and body composition phenotypes is present across a wide range of ages, and that FTO appears primarily to affect the amount of body soft tissue, influencing both fat and lean mass.

Inheritance of dermatoglyphic asymmetry in 500 Indian pedigrees: complex segregation analysis.

The major aim of this study is to determine the mode of inheritance of asymmetry of quantitative dermatoglyphic traits based on principal factors through the application of complex segregation (genetic model fitting) analyses on a large ethnically homogeneous sample of 500 Indian pedigrees (2435 individuals) of two generations. By segregation analysis of the traits- PC1_FA both Mendelian and Environmental models were rejected (< 0.001) with the General model, i.e. that despite presence of significant inheritance (rejection of Environmental model), the nature of inheritance is more complex, than Mendelian one. Although a little genetic effect was observed due to familial correlations on asymmetry traits, no evidence was found of major gene contribution to be involved, but this does not contradict the notion postulated by several earlier authors that asymmetry (fluctuating) provides a measure of developmental instability in human.

Shared Genetic Architecture Between Rheumatoid Arthritis and Varying Osteoporotic Phenotypes.

Rheumatoid arthritis (RA) and low bone mineral density (BMD), an indicator of osteoporosis (OP), appear epidemiologically associated. Shared genetic factors may explain this association. This study aimed to investigate the presence of pleiotropy to clarify the potential genetic association between RA and OP. We examined BMDs at varying skeletal sites reported in UK Biobank as well as OP fracture acquired from the Genetic Factors for Osteoporosis (GEFOS) Consortium and the TwinsUK study. PRSice-2 was used to assess the potential shared genetic overlap between RA and OP. The presence of pleiotropy was examined using colocalization analysis. PRSice-2 revealed that RA was significantly associated with OP fracture (ß = 351.6 ± 83.9, p value = 2.76E-05), total BMD (ß = -1763.5 ± 612.8, p = 4.00E-03), spine BMD (ß = -919.8 ± 264.6, p value = 5.09E-04), and forearm BMD (ß = -66.09 ± 31.40, p value = 3.53E-02). Through colocalization analysis, the same causal genetic variants, associated with both RA and OP, were apparent in 12 genes: PLCL1, BOLL, AC011997.1, TNFAIP3, RP11-158I9.1, CDK6, CHCHD4P2, RP11-505C13.1, PHF19, TRAF1, C5, and C11orf49 with moderate posterior probabilities (>50%). Pleiotropy is involved in the association between RA and OP phenotypes. These findings contribute to the understanding of disease mechanisms and provide insight into possible therapeutic advancements and enhanced screening measures. © 2021 American Society for Bone and Mineral Research (ASBMR).

Parametric model-based statistics for possible genotyping errors and sample stratification in sibling-pair SNP data.

The detection of genotyping errors, based on apparent Mendelian incompatibilities in a sample of sib-pairs, is a complicated problem. In the case of a single marker and unknown parental genotypes, all combinations of sib-pair genotypes are self-consistent. Moreover, the observed deviation from equilibrium genotype frequencies may result from genotyping errors as well as from the sample's stratification. This in turn, may profoundly affect the results of association and linkage analyses, and therefore an estimation of these factors should be done beforehand. Here we present several parametric models, and using likelihood ratio statistics, we suggest a method of combined analysis of genotyping errors and a sample stratification for randomly ascertained sib-pair single nucleotide polymorphism (SNP) data. Specifically, we implemented two models of genotyping errors in either heterozygotes or homozygotes, and two models of sample stratification resulting from either the presence of families of different ethnic origin (e.g., a population admixture) or from a different ethnic origin of the parents in the family (e.g., intermarriage). The power of this method was established by Monte Carlo data simulation. The results clearly suggest that the proposed method is most efficient for detecting genotyping errors in heterozygotes, a common error caused by incorrect SNP data interpretation. We also provide an example of its application to real data.

Lumbar disc degeneration and genetic factors are the main risk factors for low back pain in women: the UK Twin Spine Study.

OBJECTIVE: Low back pain (LBP) is a common musculoskeletal disorder, but it is still unclear which individuals develop it. The authors examined the contribution of genetic factors, lumbar disc degeneration (LDD) and other risk factors in a female sample of the general population. MATERIAL AND METHODS: A cross-sectional study was conducted among 2256 women (371 and 698 monozygotic and dizygotic twin pairs and 29 sibling pairs and 60 singletons) with a mean age of 50 years (18-84). A self-reported validated questionnaire was used to collect back pain data. Risk factors including body weight, smoking, occupation, physical exercise and MRI assessed LDD were measured. Data analysis included logistic regression and variance decomposition. RESULTS: The major factors associated with LBP included genetic background, with OR approximately 6 if the monozygotic co-twin had LBP, or 2.2 if she was a dizygotic co-twin. In addition, LDD and overweight were highly significantly (p<0.001) associated with non-specific LBP. The single most important risk factor was the amount of LDD. After adjustment for other risk factors, the individuals who exhibited advanced LDD (90% vs 10%) had 3.2 higher odds of manifesting LBP. The data also showed a significant (p<0.001) genetic correlation between the LBP and LDD measurements, suggesting that approximately 11-13% of the genetic effects are shared by LDD and LBP. CONCLUSIONS: The main risk factors for reported episodes of severe and disabling LBP in UK women include the degree of LDD as assessed by MRI, being overweight and genetic heritability.

Genetic determinants of circulating levels of tumor necrosis factor receptor II and their association with TNF-RII gene polymorphisms.

BACKGROUND: Tumor necrosis factor alpha (TNFalpha) is a cytokine involved in inflammatory, immune, and metabolic events. TNFalpha signals are mediated through activation of two receptors, one of which is tumor necrosis factor receptor TNF-RII. OBJECTIVE: To examine the effects of genetic and environmental factors on TNF-RII plasma concentration and its association with polymorphisms in the TNF-RII gene locus. METHODS: The levels of sTNF-RII were determined in 897 individuals. The association between sTNF-RII and polymorphisms in its structural gene locus was examined by pedigree-based association analyses (PDT) and transmission disequilibrium tests (TDTs). RESULTS: 49.57% of the adjusted sTNF-RII variability was attributable to genetic factors. sTNF-RII plasma levels were nominally associated with the genomic region spanning TNF-RII promoter and the first intron, represented by rs976881 (p=0.029). Although after correction for multiple testing this PDT signal formally did not reach statistical significance, it was reflected also in series of TDTs and further confirmed by association observed for haplotype of rs976881 with rs590368 (nominal p=0.006) and by ANOVA. CONCLUSIONS: sTNF-RII plasma concentration is determined by both genetic and environmental factors. Our results suggest association between sTNF-RII levels and polymorphisms in vicinity to TNF-RII promoter region. This finding requires further thorough validation in other populations.

Radiographic hand osteoarthritis in two ethnic groups living in the same geographic area.

The aim of the study was to compare the prevalence of radiographic hand osteoarthritis (OA) and its association with age, sex, body mass index in two ethnic groups, Russian and Buryats who reside in the same geographic area. It was a cross-sectional observational study. The study population comprised ethnic Russians (N = 572) and Buryats (N = 327) from the Barguzinsky District of the Buryat Republic, Russian Federation. OA was evaluated in 14 joints of the left hand according to Kellgren and Lawrence's grading system. A diagnosis of OA was determined by the number of affected joints and by the presence of at least one affected joint. Statistical analyses included prevalence estimation, linear and logistic regressions, and chi(2) tests. Our major finding was that individuals of different ethnic groups, residing in the same location, have similar prevalence and severity of radiographic hand OA. Considering the results of our previous study, where we found significant differences in the prevalence and severity of hand OA between the Russian samples (same ethnicity) with diverse places of residence, we conclude that environmental factors play an important role in the development of hand OA.

Family-based study of association between ENPP1 genetic variants and craniofacial morphology.

BACKGROUND: Human craniofacial morphology is characterized by considerable diversity among individuals. The ENPP1 gene is essential for bone physiology. However, the potential effects of its genetic variants on head size phenotypes have not yet been studied. AIM: The aim of this research was to investigate the association of polymorphisms in the ENPP1 locus with normal variability of craniofacial phenotypes. SUBJECTS AND METHODS: Fourteen SNPs and 13 haplotypes in the ENPP1 locus were tested for association with six head size traits in 1042 Western Eurasian individuals. RESULTS: The most significant and consistent association was observed between upper facial height and the polymorphisms located near the promoter region and upstream from ENPP1 gene (p = 0.00009), which remained significant after adjustment for multiple testing. Additionally, association signals were detected between head breadths and lower face height, and markers residing in or close to the promoter and 3' untranslated regions of the ENPP1 gene (p < 0.05). CONCLUSIONS: The findings obtained in this study suggest that the upstream, promoter and 3' untranslated regions in the ENPP1 locus harbor genetic variants affecting different aspects of craniofacial morphology. Further research is required to validate the relevancy of the potentially functional ENPP1 regions to normal and pathologic craniofacial growth.

Self-reported hearing loss questions provide a good measure for genetic studies: a polygenic risk score analysis from UK Biobank.

Age-related hearing impairment (ARHI) is very common in older adults and has major impact on quality of life. The heritability of ARHI has been estimated to be around 50%. The present study aimed to estimate heritability and environmental contributions to liability of ARHI and the extent to which a polygenic risk score (PRS) derived from a recent genome-wide association study of questionnaire items regarding hearing loss using the UK Biobank is predictive of hearing loss in other samples. We examined (1) a sample from TwinsUK who have had hearing ability measured by pure-tone audiogram and the speech-to-noise ratio test as well as questionnaire measures that are comparable with the UK Biobank questionnaire items and (2) European and non-European samples from the UK Biobank which were not part of the original GWAS. Results indicated that the questionnaire items were over 50% heritable in TwinsUK and comparable with the objective hearing measures. In addition, we found very high genetic correlation (0.30-0.84) between the questionnaire responses and objective hearing measures in the TwinsUK sample. Finally, PRS computed from weighted UK Biobank GWAS results were predictive of both questionnaire and objective measures of hearing loss in the TwinsUK sample, as well as questionnaire-measured hearing loss in Europeans but not non-European subpopulations. These results demonstrate the utility of questionnaire-based methods in genetic association studies of hearing loss in adults and highlight the differences in genetic predisposition to ARHI by ethnic background.

The cannabinoid receptor type 2 (CNR2) gene is associated with hand bone strength phenotypes in an ethnically homogeneous family sample.

Genetic variants within the CNR2 gene encoding the cannabinoid receptor CB2 have been shown to be associated with osteoporosis and low bone mineral density (BMD) in case-control studies. We now examined the association of polymorphisms in CNR2 with hand bone strength in an ethnically homogeneous healthy family sample of European origin (Chuvashians) living in Russia. We show that non-synonymous CNR2 SNPs are significantly associated with radiographic hand BMD and breaking bending resistance index (BBRI) by two different transmission disequilibrium tests. For both tests highly significant p values (ranging from 0.007 to 0.008 for hand BMD, and from 0.001 to 0.003 for BBRI) were also obtained with additional SNPs at the CNR2 locus. The associations remained significant after correction for multiple testing. In conclusion, in addition to the association of CNR2 polymorphisms with low BMD at selected clinically relevant skeletal sites, we now report their significant association with hand bone strength phenotypes using a family-based study design implying an even broader impact of genetic variation at the CNR2 locus on bone structure and function.

Anthropometric and bone-related biochemical factors are associated with different haplotypes of ANKH locus.

BACKGROUND: The human homologue of the mouse progressive ankylosis (ANKH) gene is one of the key genetic factors involved in bone mineralization. Previous studies have shown that plasma levels of osteoprotegerin (OPG) and parathyroid hormone (PTH) are associated with the distal region of the ANKH gene, whereas skeletal size measurements are associated with the promoter region. AIM: The present study examines the possible phenotype-haplotype specificity of the associations in these two gene regions. SUBJECTS AND METHODS: The total sample consists of 1249 healthy individuals (mean age = 47.7, SD = 16.8) from 404 nuclear families. Fifteen interrelated anthropometric measurements were transformed into two principal components, reflecting body size and mass. Those, plus circulating levels of PTH and OPG, were subjected to association analysis, using transmission disequilibrium tests (TDTs) with ANKH gene. From 805 to 1150 individuals per SNP were genotyped. RESULTS: In the proximal region (rs3006069-rs835154-rs835141), associations were found between the A-A-C haplotype and the first principal component reflecting body size (p < or = 0.048), whereas another haplotype, G-G-C, was associated with the first principal component, reflecting the body mass (p < or = 0.008). In the distal region of ANKH (rs39968-rs696294-rs875525), the A-A-C haplotype was found to be associated with OPG plasma levels (p < or = 0.001), whereas the G-A-C haplotype was associated with PTH circulating concentrations (p < or = 0.025). CONCLUSION: Taken together, the results show discrimination between the corresponding regions and haplotypes, suggesting trait-specific gene variants that influenced bone-related phenotypic variation in the studied population.

Multi-OMICS analyses of frailty and chronic widespread musculoskeletal pain suggest involvement of shared neurological pathways.

Chronic widespread musculoskeletal pain (CWP) and frailty are prevalent conditions in older people. We have shown previously that interindividual variation in frailty and CWP is genetically determined. We also reported an association of frailty and CWP caused by shared genetic and common environmental factors. The aim of this study was to use omic approaches to identify molecular genetic factors underlying the heritability of frailty and its genetic correlation with CWP. Frailty was quantified through the Rockwood Frailty Index (FI) as a proportion of deficits from 33 binary health deficit questions in 3626 female twins. Common widespread pain was assessed using a screening questionnaire. OMICS analysis included 305 metabolites and whole-genome (>2.5 × 10 SNPs) and epigenome (∼1 × 10 MeDIP-seq regions) assessments performed on fasting blood samples. Using family-based statistical analyses, including path analysis, we examined how FI scores were related to molecular genetic factors and to CWP, taking into account known risk factors such as fat mass and smoking. Frailty Index was significantly correlated with 51 metabolites after correction for multiple testing, with 20 metabolites having P-values between 2.1 × 10 and 4.0 × 10. Three metabolites (uridine, C-glycosyl tryptophan, and N-acetyl glycine) were statistically independent and thought to exert a direct effect on FI. Epiandrosterone sulphate, previously shown to be highly inversely associated with CWP, was found to exert an indirect influence on FI. Bioinformatics analysis of genome-wide association study and EWAS showed that FI and its covariation with CWP was through genomic regions involved in neurological pathways. Neurological pathway involvement accounts for the associated conditions of aging CWP and FI.