Bivariate genome-wide association study suggests that the DARC gene influences lean body mass and age at menarche

Lean body mass (LBM) and age at menarche (AAM) are two important complex traits for human health. The aim of this study was to identify pleiotropic genes for both traits using a powerful bivariate genome-wide association study (GWAS). Two studies, a discovery study and a replication study, were performed. In the discovery study, 909622 single nucleotide polymorphisms (SNPs) were genotyped in 801 unrelated female Han Chinese subjects using the Affymetrix human genome-wide SNP array 6.0 platform. Then, a bivariate GWAS was performed to identify the SNPs that may be important for LBM and AAM. In the replication study, significant findings from the discovery study were validated in 1692 unrelated Caucasian female subjects. One SNP rs3027009 that was bivariately associated with left arm lean mass and AAM in the discovery samples (P=7.26×10^?6) and in the replication samples (P=0.005) was identified. The SNP is located at the upstream of DARC (Duffy antigen receptor for chemokines) gene, suggesting that DARC may play an important role in regulating the metabolisms of both LBM and AAM.     

Bivariate whole-genome linkage scan for bone geometry and total body fat mass

To quantify the genetic correlations between total body fat mass （TBFM） and femoral neck geometric parameters （FNGPs） and, if pos- sible, to detect the specific genomic regions shared by them, bivariate genetic analysis and bivariate whole-genome linkage scan were carried out in a large Caucasian population. All the phenotypes studied were significantly controlled by genetic factors （P 〈 0.001） with the heritabilities ranging from 0.45 to 0.68. Significantly genetic correlations were found between TBFM and CSA （cross-section area）, W （sub-periosteal diameter）, Z （section modulus） and CT （cortical thickness） except between TBFM and BR （buckling ratio）. The peak bivariate LOD scores were 3.23 （20q12）, 2.47 （20p11）, 3.19 （6q27）, 1.68 （20p12）, and 2.47 （7q11） for the five pairs of TBFM and BR, CSA, CT, W, and Z in the entire sample, respectively. Gender-specific bivariate linkage evidences were also found for the five pairs. 6p25 had complete pleiotropic effects on the variations of TBFM ＆ Z in the female sub-population, and 6q27 and 17q11 had coincident link- ages for TBFM ＆ CSA and TBFM ＆ Z in the entire population. We identified moderate genetic correlations and several shared genomic regions between TBFM and FNGPs in a large Caucasian population.     

Suggestion of GLYAT gene underlying variation of bone size and body lean mass as revealed by a bivariate genome-wide association study

Bone and muscle, two major tissue types of musculoskeletal system, have strong genetic determination. Abnormality in bone and/or muscle may cause musculoskeletal diseases such as osteoporosis and sarcopenia. Bone size phenotypes (BSPs), such as hip bone size (HBS), appendicular bone size (ABS), are genetically correlated with body lean mass (mainly muscle mass). However, the specific genes shared by these phenotypes are largely unknown. In this study, we aimed to identify the specific genes with pleiotropic effects on BSPs and appendicular lean mass (ALM). We performed a bivariate genome-wide association study (GWAS) by analyzing ~690,000 SNPs in 1,627 unrelated Han Chinese adults (802 males and 825 females) followed by a replication study in 2,286 unrelated US Caucasians (558 males and 1,728 females). We identified 14 interesting single nucleotide polymorphisms (SNPs) that may contribute to variation of both BSPs and ALM, with p values <10^?6 in discovery stage. Among them, the association of three SNPs (rs2507838, rs7116722, and rs11826261) in/near GLYAT (glycine-N-acyltransferase) gene was replicated in US Caucasians, with p values ranging from 1.89 × 10^?3 to 3.71 × 10^?4 for ALM–ABS, from 5.14 × 10^?3 to 1.11 × 10^?2 for ALM–HBS, respectively. Meta-analyses yielded stronger association signals for rs2507838, rs7116722, and rs11826261, with pooled p values of 1.68 × 10^?8, 7.94 × 10^?8, 6.80 × 10^?8 for ALB–ABS and 1.22 × 10^?4, 9.85 × 10^?5, 3.96 × 10^?4 for ALM–HBS, respectively. Haplotype allele ATA based on these three SNPs was also associated with ALM–HBS and ALM–ABS in both discovery and replication samples. Interestingly, GLYAT was previously found to be essential to glucose metabolism and energy metabolism, suggesting the gene’s dual role in both bone development and muscle growth. Our findings, together with the prior biological evidence, suggest the importance of GLYAT gene in co-regulation of bone phenotypes and body lean mass.     

Development of a surrogate model based on patient weight,bone mass and geometry to predict femoral neck strains and fracture loads

Osteoporosis and related bone fractures are an increasing global burden in our ageing society. Areal bone mineral density assessed through dual energy X-ray absorptiometry (DEXA), the clinically accepted and most used method, is not sufficient to assess fracture risk individually. Finite element (FE) modelling has shown improvements in prediction of fracture risk, better than aBMD from DEXA, but is not practical for widespread clinical use. The aim of this study was to develop an adaptive neural network (ANN)-based surrogate model to predict femoral neck strains and fracture loads obtained from a previously developed population-based FE model. The surrogate model performance was assessed in simulating two loading conditions: the stance phase of gait and a fall.The surrogate model successfully predicted strains estimated by FE (r² = 0.90–0.98 for level gait load case, r² = 0.92–0.96 for the fall load case). Moreover, an ANN model based on three measurements obtainable in clinics (femoral neck length (level gait) or maximum femoral neck diameter (fall), femoral neck bone mass, body weight) was able to give reasonable predictions (r² = 0.84–0.94) for all of the strain metrics and the estimated femoral neck fracture load. Overall, the surrogate model has potential for clinical applications as they are based on simple measures of geometry and bone mass which can be derived from DEXA images, accurately predicting FE model outcomes, with advantages over FE models as they are quicker and easier to perform.     

Fracture mechanics of the femoral neck in a composite bone model: Effects of platen geometry

Load applicator (platen) geometry used for axial load to failure testing of the femoral neck varies between studies and the biomechanical consequences are unknown. The purpose of this study was to determine if load application with a flat versus a conical platen results in differing fracture mechanics. Femurs were aligned in 25° of adduction and an axial compressive force was applied to the femoral heads at a rate of 6 mm/min until failure. Load application with the conical platen resulted in an average ultimate failure load, stiffness, and energy to failure of 9067 N, 4033 N/mm, and 12.12 J, respectively. Load application with the flat platen resulted in a significant (p<0.05) reduction in ultimate failure load (7620 N) and stiffness (2924 N/mm). Energy to failure (12.30 J) was not significantly different (p=0.893). Different fracture patterns were observed for the two platens and the conical platen produced fractures more similar to clinical observations. Use of a flat platen underestimates the strength and stiffness of the femoral neck and inaccurately predicts the associated fracture pattern. These findings must be considered when interpreting the results of prior biomechanical studies on femoral neck fracture and for the development of future femoral neck fracture models.     

A genome-wide association analysis implicates SOX6 as a candidate gene for wrist bone mass

Osteoporosis is a highly heritable common bone disease leading to fractures that severely impair the life quality of patients. Wrist fractures caused by osteoporosis are largely due to the scarcity of wrist bone mass. Here we report the results of a genome-wide association study (GWAS) of wrist bone mineral density (BMD). We examined ～500000 SNP markers in 1000 unrelated homogeneous Caucasian subjects and found a novel allelic association with wrist BMD at rs11023787 in the SOX6 (SRY (sex determining region Y)-box 6) gene (P=9.00×10(-5)). Subjects carrying the C allele of rs11023787 in SOX6 had significantly higher mean wrist BMD values than those with the T allele (0.485:0.462 g cm(-2) for C allele vs. T allele carriers). For validation, we performed SOX6 association for BMD in an independent Chinese sample and found that SNP rs11023787 was significantly associated with wrist BMD in the Chinese sample (P=6.41×10(-3)). Meta-analyses of the GWAS scan and the replication studies yielded P-values of 5.20×10(-6) for rs11023787. Results of this study, together with the functional relevance of SOX6 in cartilage formation, support the SOX6 gene as an important gene for BMD variation.     

BioSMACK: a linux live CD for genome-wide association analyses

Recent advances in high-throughput genotyping technologies have enabled us to conduct a genome-wide association study (GWAS) on a large cohort. However, analyzing millions of single nucleotide polymorphisms (SNPs) is still a difficult task for researchers conducting a GWAS. Several difficulties such as compatibilities and dependencies are often encountered by researchers using analytical tools, during the installation of software. This is a huge obstacle to any research institute without computing facilities and specialists. Therefore, a proper research environment is an urgent need for researchers working on GWAS. We developed BioSMACK to provide a research environment for GWAS that requires no configuration and is easy to use. BioSMACK is based on the Ubuntu Live CD that offers a complete Linux-based operating system environment without installation. Moreover, we provide users with a GWAS manual consisting of a series of guidelines for GWAS and useful examples. BioSMACK is freely available at http://ksnp.cdc. go.kr/biosmack.     

Bioimpedance analysis: a useful technique for assessing appendicular lean soft tissue mass and distribution.

The present study was aimed at evaluating the feasibility and reliability of lower limb skeletal muscle (SM) mass estimates obtained by bioimpedance analysis (BIA). BIA estimates were compared with the estimates obtained by dual-energy X-ray absorptiometry (DXA). Ten normal weight and 10 obese women had BIA and DXA evaluations. Lower limb SM mass was then derived from DXA appendicular lean soft tissue estimates. Lower limb SM mass and SM distribution were also estimated from BIA modeling that fits measured resistance values along the leg. SM mass (mean +/- SD) was 5.8 +/- 1.0 kg by BIA vs. 5.8 +/- 1.1 kg by DXA in normal weight subjects and 7.2 +/- 1.4 kg by BIA vs. 7.2 +/- 1.2 kg by DXA in obese subjects. Mean +/- SD of the absolute value of the relative error was 7.0 +/- 3.4 and 5.9 +/- 3.4% in the two groups, respectively. Similar results were obtained by using five resistance values for the analysis. In conclusion, the proposed BIA model provides an adequate means of evaluating appendicular SM mass.     

Pathway analysis of genome-wide association study for bone mineral density

The aim of this study was to identify the candidate causal single nucleotide polymorphisms (SNPs) and candidate causal mechanisms that contribute to bone mineral density (BMD) and to generate a SNP to gene to pathway hypothesis using an analytical pathway-based approach. We used hip BMD GWAS data of the genotypes of 301,019 SNPs in 5,715 Europeans. ICSNPathway (identify candidate causal SNPs and pathways) analysis was applied to the BMD GWAS dataset. The first stage involved the pre-selection of candidate causal SNPs by linkage disequilibrium analysis and the functional SNP annotation of the most significant SNPs found. The second stage involved the annotation of biological mechanisms for the pre-selected candidate causal SNPs using improved-gene set enrichment analysis. ICSNPathway analysis identified seven candidate SNPs, eight candidate pathways, and seven hypothetical biological mechanisms. Eight pathways are as follows; gamma-hexachlorocyclohexane degradation (nominal p-value < 0.001, false discovery rate (FDR) <0.001), regulation of the smoothened signaling pathway (nominal p-value < 0.001, FDR = 0.016), TACI and BCMA stimulation of B cell immune response (nominal p-value < 0.001, FDR = 0.021), endonuclease activity (nominal p-value = 0.001, FDR = 0,026), regulation of defense response to virus (nominal p-value = 0.001, FDR = 0.028), serine_type_endopeptidase_inhibitor_activity (nominal p-value = 0.001, FDR = 0.044), endoribonuclease activity (nominal p-value = 0.002, FDR = 0.045), and myeloid leukocyte differentiation (nominal p-value = 0.001, FDR = 0.050). The most significant causal pathway was gamma-hexachlorocyclohexane degradation. CYP3A5, PON2, PON3, CMBL, PON1, ALPL, CYP3A43, CYP3A7, ACP6, ACPP, and ALPI (p < 0.05) are involved in the pathway of gamma-hexachlorocyclohexane degradation. Further examination of the gene contents revealed that DBR1, DICER1, EXO1, FEN1, POP1, POP4, RPP30, and RPP38 were involved in 2 of the 8 pathways (p < 0.05). By applying ICSNPathway analysis to BMD GWAS data, we identified seven candidate SNPs and eight pathways involving gamma-hexachlorocyclohexane degradation, which may contribute to low BMD.     

Regulation of lean mass, bone mass, and exercise tolerance by the central melanocortin system

Signaling via the type 4-melanocortin receptor (MC4R) is an important determinant of body weight in mice and humans, where loss of function mutations lead to significant obesity. Humans with mutations in the MC4R experience an increase in lean mass. However, the simultaneous accrual of fat mass in such individuals may contribute to this effect via mechanical loading. We therefore examined the relationship of fat mass and lean mass in mice lacking the type-4 melanocortin receptor (MC4RKO). We demonstrate that MC4RKO mice display increased lean body mass. Further, this is not dependent on changes in adipose mass, as MC4RKO mice possess more lean body mass than diet-induced obese (DIO) wild type mice with equivalent fat mass. To examine potential sources of the increased lean mass in MC4RKO mice, bone mass and strength were examined in MC4RKO mice. Both parameters increase with age in MC4RKO mice, which likely contributes to increases in lean body mass. We functionally characterized the increased lean mass in MC4RKO mice by examining their capacity for treadmill running. MC4R deficiency results in a decrease in exercise performance. No changes in the ratio of oxidative to glycolytic fibers were seen, however MC4RKO mice demonstrate a significantly reduced heart rate, which may underlie their impaired exercise performance. The reduced exercise capacity we report in the MC4RKO mouse has potential clinical ramifications, as efforts to control body weight in humans with melanocortin deficiency may be ineffective due to poor tolerance for physical activity.     

Cortical bone distribution in the femoral neck of strepsirhine primates

The thickness of the inferior and superior cortices of the femoral neck was measured on X-rays of 181 strepsirhine primate femora representing 24 species. Neck length, neck depth and neck-shaft angle were also measured. The strength of the femoral neck in frontal bending was estimated by modeling the neck as a hollow cylinder, with neck depth as the outer diameter and cortical thickness representing the superior and inferior shell dimensions. Results indicate that the inferior cortex is always thicker than the superior cortex. The ratio of superior to inferior cortical thickness is highly variable but distinguishes two of the three locomotor groups in the sample. Vertical clingers and leapers have higher ratios (i.e., a more even distribution of cortical bone) than quadrupeds. The slow climbers tend to have the lowest ratios, although they do not differ significantly from the leapers and quadrupeds. These results do not confirm prior theoretical expectations and reported data for anthropoid primates that link greater asymmetry of the cortical shell to more stereotypical hip excursions. The ratio of superior to inferior cortical thickness is unrelated to body mass, femoral neck length, and neck-shaft angle, calling into question whether the short neck of strepsirhine primates acts as a cantilever beam in bending. On the other hand, the estimated section moduli are highly correlated with body mass and neck length, a correlation that is driven primarily by body mass. In conclusion, we believe that an alternative interpretation to the cantilever beam model is needed to explain the asymmetry in bone distribution in the femoral neck, at least in strepsirhine primates (e.g., a thicker inferior cortex is required to reinforce the strongly curved inferior surface). As in prior studies of cross-sectional geometry of long bones, we found slightly positive allometry of cortical dimensions with body mass.     

Cell adhesion molecules contribute to Alzheimer's disease: multiple pathway analyses of two genome-wide association studies

Alzheimer's disease (AD) is a kind of complex neurological disorder. The complex genetic architecture of AD makes genetic analysis difficult. Fortunately, a pathway-based method to study the existing genome-wide association studies datasets has been applied into AD. However, no shared Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway was reported. In this study, we performed multiple pathway analyses of French AD genome-wide association studies dataset (discovery dataset, n = 7360, 2032 cases and 5328 controls) and Pfizer dataset (validation dataset, n = 2220, 1034 cases and 1186 controls). First, we performed multiple pathway analyses by Hypergeometric test, improved gene set enrichment analysis (IGSEA) and Z-statistic test in KEGG. Using Hypergeometric test, we identified 54 and 25 significant pathways (p < 0.05) in discovery dataset and validation dataset, respectively. Using IGSEA method, we identified three significant pathways in both discovery and validation datasets, respectively. Using Z-statistic test, we identified 19 significant pathways in validation dataset. Among the significant pathways, cell adhesion molecules (CAM) pathway was identified to be the only consistent signal emerging across multiple analyses in KEGG. After permutation and multiple testing corrections, CAM pathway was significant with p = 2.40E-05 (Hypergeometric test) and p = 3.00E-03 (IGSEA) in discovery dataset. In validation dataset, CAM pathway was significant with p = 1.84E-06 (Hypergeometric test), p = 1.00E-02 (IGSEA) and p = 2.81E-03 (Z-statistic test). We replicated the association by multiple pathway analyses in Gene Ontology using Hypergeometric test (WebGestalt), modified Fisher's exact test (DAVID) and Binomial test (PANTHER). Our findings provided further evidence on the association between CAM pathway and AD susceptibility, which would be helpful to study the genetic mechanisms of AD and may significantly assist in the development of therapeutic strategies.     

Geometry and bone mineral density determinants of femoral neck strength changes following exercise

Physical exercise induces spatially heterogeneous adaptation in bone. However, it remains unclear where the changes in BMD and geometry have the greatest impact on femoral neck strength. The aim of this study was to determine the principal BMD-and-geometry changes induced by exercise that have the greatest effect on femoral neck strength. Pre- and post-exercise 3D-DXA images of the proximal femur were collected of male participants from the LIFTMOR-M exercise intervention trial. Meshes with element-by-element correspondence were generated by morphing a template mesh to each bone to calculate changes in BMD and geometry. Finite element (FE) models predicted femoral neck strength changes under single-leg stance and sideways fall load. Partial least squares regression (PLSR) models were developed with BMD-only, geometry-only, and BMD-and-geometry changes to determine the principal modes that explained the greatest variation in neck strength changes. The PLSR models explained over 90% of the strength variation with 3 PLS components using BMD-only (R² > 0.92, RMSE < 0.06 N) and 8 PLS components with geometry-only (R² > 0.93, RMSE < 0.06 N). Changes in the superior neck and distal cortex were most important during single-leg stance while the superior neck, medial head, and lateral trochanter were most important during a sideways fall. Local changes in femoral neck and head geometry could differentiate the exercise groups from the control group. Exercise interventions may target BMD changes in the superior neck, inferior neck, and greater trochanter for improved femoral neck strength in single-leg stance and sideways fall.

     

Alpine skiing is associated with higher femoral neck bone mineral density

Objective:

To evaluate the influence of elite-level alpine skiing on athletes’ skeleton.

Methods:

Thirteen professional alpine skiers (9 males and 4 females with mean age of 22.6 years) and their age- and height matched control subjects were measured with dual energy X-ray absorptiometry (total body, lumbar spine, proximal femur, forearm) and quantitative ultrasound (hand).

Results:

After adjusting for sex, age, weight and height, between-group differences were 15% (p=0.012) for the lumbar spine, 14% (p=0.022) for the femoral neck, 10% (p=0.051) for the total hip, and 11% (p=0.001) for the total body favoring the alpine skiers. However, after controlling for total body lean mass (~muscle mass), the group-differences lost their statistical significance, the borderline 10% difference (p=0.051) in femoral neck BMD excluded.

Conclusion:

Factors contributing to the alpine skiers’ higher BMD may not only include the greater muscle mass (~stronger muscles) of these athletes but also a large number of impacts and possibly other high-frequency features in external loading generated by the high-speed skiing performance.     

The meta-analysis of genome-wide association studies

The pressure to publish novel genetic associations has meant that meta-analysis has been applied to genome-wide association studies without the time for a careful consideration of the methods that are used. This review distinguishes between the use of meta-analysis to validate previously reported genetic associations and its use for gene discovery, and advocates viewing gene discovery as an exploratory screen that requires independent replication instead of treating it as the application of hundreds of thousands of statistical tests. The review considers the use of fixed and random effects meta-analyses, the investigation of between-study heterogeneity, adjustment for confounding, assessing the combined evidence and genomic control, and comments on alternative approaches that have been used in the literature.     

Gene-centric characteristics of genome-wide association studies

Background

The high-throughput genotyping chips have contributed greatly to genome-wide association (GWA) studies to identify novel disease susceptibility single nucleotide polymorphisms (SNPs). The high-density chips are designed using two different SNP selection approaches, the direct gene-centric approach, and the indirect quasi-random SNPs or linkage disequilibrium (LD)-based tagSNPs approaches. Although all these approaches can provide high genome coverage and ascertain variants in genes, it is not clear to which extent these approaches could capture the common genic variants. It is also important to characterize and compare the differences between these approaches.

Methodology/Principal Findings

In our study, by using both the Phase II HapMap data and the disease variants extracted from OMIM, a gene-centric evaluation was first performed to evaluate the ability of the approaches in capturing the disease variants in Caucasian population. Then the distribution patterns of SNPs were also characterized in genic regions, evolutionarily conserved introns and nongenic regions, ontologies and pathways. The results show that, no mater which SNP selection approach is used, the current high-density SNP chips provide very high coverage in genic regions and can capture most of known common disease variants under HapMap frame. The results also show that the differences between the direct and the indirect approaches are relatively small. Both have similar SNP distribution patterns in these gene-centric characteristics.

Conclusions/Significance

This study suggests that the indirect approaches not only have the advantage of high coverage but also are useful for studies focusing on various functional SNPs either in genes or in the conserved regions that the direct approach supports. The study and the annotation of characteristics will be helpful for designing and analyzing GWA studies that aim to identify genetic risk factors involved in common diseases, especially variants in genes and conserved regions.