New genetic loci link adipose and insulin biology to body fat distribution

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10−8). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Genome-wide enrichment analysis between endometriosis and obesity-related traits reveals novel susceptibility loci

Endometriosis is a chronic inflammatory condition in women that results in pelvic pain and subfertility, and has been associated with decreased body mass index (BMI). Genetic variants contributing to the heritable component have started to emerge from genome-wide association studies (GWAS), although the majority remain unknown. Unexpectedly, we observed an intergenic locus on 7p15.2 that was genome-wide significantly associated with both endometriosis and fat distribution (waist-to-hip ratio adjusted for BMI; WHRadjBMI) in an independent meta-GWAS of European ancestry individuals. This led us to investigate the potential overlap in genetic variants underlying the aetiology of endometriosis, WHRadjBMI and BMI using GWAS data. Our analyses demonstrated significant enrichment of common variants between fat distribution and endometriosis (P = 3.7 x 10(-3)), which was stronger when we restricted the investigation to more severe (Stage B) cases (P = 4.5 x 10(-4)). However, no genetic enrichment was observed between endometriosis and BMI (P = 0.79). In addition to 7p15.2, we identify four more variants with statistically significant evidence of involvement in both endometriosis and WHRadjBMI (in/near KIFAP3, CAB39L, WNT4, GRB14); two of these, KIFAP3 and CAB39L, are novel associations for both traits. KIFAP3, WNT4 and 7p15.2 are associated with the WNT signalling pathway; formal pathway analysis confirmed a statistically significant (P = 6.41 x 10(-4)) overrepresentation of shared associations in developmental processes/WNT signalling between the two traits. Our results demonstrate an example of potential biological pleiotropy that was hitherto unknown, and represent an opportunity for functional follow-up of loci and further cross-phenotype comparisons to assess how fat distribution and endometriosis pathogenesis research fields can inform each other.

Anthropometric measures of obesity-their association with type 2 diabetes and hypertension across ethnic groups

Clinical trials have shown that weight reduction with lifestyles can delay or prevent diabetes and reduce blood pressure. An appropriate definition of obesity using anthropometric measures is useful in predicting diabetes and hypertension at the population level. However, there is debate on which of the measures of obesity is best or most strongly associated with diabetes and hypertension and on what are the optimal cut-off values for body mass index (BMI) and waist circumference (WC) in this regard. The aims of the study were 1) to compare the strength of the association for undiagnosed or newly diagnosed diabetes (or hypertension) with anthropometric measures of obesity in people of Asian origin, 2) to detect ethnic differences in the association of undiagnosed diabetes with obesity, 3) to identify ethnic- and sex-specific change point values of BMI and WC for changes in the prevalence of diabetes and 4) to evaluate the ethnic-specific WC cutoff values proposed by the International Diabetes Federation (IDF) in 2005 for central obesity. The study population comprised 28 435 men and 35 198 women, ≥ 25 years of age, from 39 cohorts participating in the DECODA and DECODE studies, including 5 Asian Indian (n = 13 537), 3 Mauritian Indian (n = 4505) and Mauritian Creole (n = 1075), 8 Chinese (n =10 801), 1 Filipino (n = 3841), 7 Japanese (n = 7934), 1 Mongolian (n = 1991), and 14 European (n = 20 979) studies. The prevalence of diabetes, hypertension and central obesity was estimated, using descriptive statistics, and the differences were determined with the χ2 test. The odds ratios (ORs) or  coefficients (from the logistic model) and hazard ratios (HRs, from the Cox model to interval censored data) for BMI, WC, waist-to-hip ratio (WHR), and waist-to-stature ratio (WSR) were estimated for diabetes and hypertension. The differences between BMI and WC, WHR or WSR were compared, applying paired homogeneity tests (Wald statistics with 1 df). Hierarchical three-level Bayesian change point analysis, adjusting for age, was applied to identify the most likely cut-off/change point values for BMI and WC in association with previously undiagnosed diabetes. The ORs for diabetes in men (women) with BMI, WC, WHR and WSR were 1.52 (1.59), 1.54 (1.70), 1.53 (1.50) and 1.62 (1.70), respectively and the corresponding ORs for hypertension were 1.68 (1.55), 1.66 (1.51), 1.45 (1.28) and 1.63 (1.50). For diabetes the OR for BMI did not differ from that for WC or WHR, but was lower than that for WSR (p = 0.001) in men while in women the ORs were higher for WC and WSR than for BMI (both p < 0.05). Hypertension was more strongly associated with BMI than with WHR in men (p < 0.001) and most strongly with BMI than with WHR (p < 0.001), WSR (p < 0.01) and WC (p < 0.05) in women. The HRs for incidence of diabetes and hypertension did not differ between BMI and the other three central obesity measures in Mauritian Indians and Mauritian Creoles during follow-ups of 5, 6 and 11 years. The prevalence of diabetes was highest in Asian Indians, lowest in Europeans and intermediate in others, given the same BMI or WC category. The  coefficients for diabetes in BMI (kg/m2) were (men/women): 0.34/0.28, 0.41/0.43, 0.42/0.61, 0.36/0.59 and 0.33/0.49 for Asian Indian, Chinese, Japanese, Mauritian Indian and European (overall homogeneity test: p > 0.05 in men and p < 0.001 in women). Similar results were obtained in WC (cm). Asian Indian women had lower  coefficients than women of other ethnicities. The change points for BMI were 29.5, 25.6, 24.0, 24.0 and 21.5 in men and 29.4, 25.2, 24.9, 25.3 and 22.5 (kg/m2) in women of European, Chinese, Mauritian Indian, Japanese, and Asian Indian descent. The change points for WC were 100, 85, 79 and 82 cm in men and 91, 82, 82 and 76 cm in women of European, Chinese, Mauritian Indian, and Asian Indian. The prevalence of central obesity using the 2005 IDF definition was higher in Japanese men but lower in Japanese women than in their Asian counterparts. The prevalence of central obesity was 52 times higher in Japanese men but 0.8 times lower in Japanese women compared to the National Cholesterol Education Programme definition. The findings suggest that both BMI and WC predicted diabetes and hypertension equally well in all ethnic groups. At the same BMI or WC level, the prevalence of diabetes was highest in Asian Indians, lowest in Europeans and intermediate in others. Ethnic- and sex-specific change points of BMI and WC should be considered in setting diagnostic criteria for obesity to detect undiagnosed or newly diagnosed diabetes.

Physical activity, fitness, abdominal obesity, and cardiovascular risk factors in Finnish men and women : The National FINRISK 2002 Study

Physical inactivity, low cardiorespiratory fitness, and abdominal obesity are direct and mediating risk factors for cardiovascular disease (CVD). The results of recent studies suggest that individuals with higher levels of physical activity or cardiorespiratory fitness have lower CVD and all-cause mortality than those with lower activity or fitness levels regardless of their level of obesity. The interrelationships of physical activity, fitness, and abdominal obesity with cardiovascular risk factors have not been studied in detail. The aim of this study was to investigate the associations of different types of leisure time physical activity and aerobic fitness with cardiovascular risk factors in a large population of Finnish adults. In addition, a novel aerobic fitness test was implemented and the distribution of aerobic fitness was explored in men and women across age groups. The interrelationships of physical activity, aerobic fitness and abdominal obesity were examined in relation to cardiovascular risk factors. This study was part of the National FINRISK Study 2002, which monitors cardiovascular risk factors in a Finnish adult population. The sample comprised 13 437 men and women aged 25 to 74 years and was drawn from the Population Register as a stratified random sample according to 10-year age groups, gender and area. A separate physical activity study included 9179 subjects, of whom 5 980 participated (65%) in the study. At the study site, weight, height, waist and hip circumferences, and blood pressure were measured, a blood sample was drawn, and an aerobic fitness test was performed. The fitness test estimated maximal oxygen uptake (VO2max) and was based on a non-exercise method by using a heart rate monitor at rest. Waist-to-hip ratio (WHR) was calculated by dividing waist circumference with hip circumference and was used as a measure of abdominal obesity. Participants filled in a questionnaire on health behavior, a history of diseases, and current health status, and a detailed 12-month leisure time physical activity recall. Based on the recall data, relative energy expenditure was calculated using metabolic equivalents, and physical activity was divided into conditioning, non-conditioning, and commuting physical activity. Participants aged 45 to 74 years were later invited to take part in a 2-hour oral glucose tolerance test with fasting insulin and glucose measurements. Based on the oral glucose tolerance test, undiagnosed impaired glucose tolerance and type 2 diabetes were defined. The estimated aerobic fitness was lower among women and decreased with age. A higher estimated aerobic fitness and a lower WHR were independently associated with lower systolic and diastolic blood pressure, lower total cholesterol and triglyceride levels, and with higher high-density lipoprotein (HDL) cholesterol and HDL to total cholesterol ratio. The associations of the estimated aerobic fitness with diastolic blood pressure, triglycerides, and HDL to total cholesterol ratio were stronger in men with a higher WHR. High levels of conditioning and non-conditioning physical activity were associated with lower high-sensitivity C-reactive protein (CRP) levels. High levels of conditioning and overall physical activities were associated with lower insulin and glucose levels. The associations were stronger among women than men. A better self-rated physical fitness was associated with a higher estimated aerobic fitness, lower CRP levels, and lower insulin and glucose levels in men and women. In each WHR third, the risk of impaired glucose tolerance and type 2 diabetes was higher among physically inactive individuals who did not undertake at least 30 minutes of moderate-intensity physical activity on five days per week. These cross-sectional data show that higher levels of estimated aerobic fitness and regular leisure time physical activity are associated with a favorable cardiovascular risk factor profile and that these associations are present at all levels of abdominal obesity. Most of the associations followed a dose-response manner, suggesting that already low levels of physical activity or fitness are beneficial to health and that larger improvements in risk factor levels may be gained from higher activity and fitness levels. The present findings support the recommendation to engage regularly in leisure time physical activity, to pursue a high level of aerobic fitness, and to prevent abdominal obesity.

WinBUGS code to run a network meta-analysis

WinBUGS code and data to reproduce our network meta-analysis from "Control strategies to prevent total hip replacement-related infections: a systematic review and mixed treatment comparison" published in BMJ Open.

Male traits and herd reproductive capability in tropical beef cattle. 1. Experimental design and animal measures

Research into the genetics of whole herd profitability has been a focus of the Beef Cooperative Research Centre for Beef Genetic Technologies over the past decade and it has been identified that measures of male reproduction may offer a potential indirect means of selecting for improved female reproduction. This paper describes the experimental design and provides a descriptive analysis of an array of male traits in Brahman and Tropical Composite genotypes managed under the medium to high stress, semi-extensive to extensive production systems of northern Australia. A total of 1639 Brahman and 2424 Tropical Composite bulls with known pedigrees, bred and raised in northern Australia, were evaluated for a comprehensive range of productive and reproductive traits. These included blood hormonal traits (luteinising hormone, inhibin and insulin-like growth factor-I); growth and carcass traits (liveweight, body condition score, ultrasound scanned 12-13th rib fat, rump P8 fat, eye muscle area and hip height); adaptation traits (flight time and rectal temperature); and a bull breeding soundness evaluation (leg and hoof conformation, sheath score, length of everted prepuce, penile anatomy, scrotal circumference, semen mass activity, sperm motility and sperm morphology). Large phenotypic variation was evident for most traits, with complete overlap between genotypes, indicating that there is likely to be a significant opportunity to improve bull fertility traits through management and bull selection.

Male traits and herd reproductive capability in tropical beef cattle. 2. Genetic parameters of bull traits

A total of 4063 young bulls of two tropical genotypes (1639 Brahman and 2424 Tropical Composite) raised in northern Australia were evaluated for a comprehensive range of production and reproduction traits up to 24 months of age. Prior to weaning, peripheral blood concentrations of luteinising hormone (LH) and inhibin were measured at 4 months of age. At weaning (6 months) blood insulin-like growth factor-1 (IGF-I) and flight time were recorded. Body composition traits of fat depth and eye-muscle area were determined by ultrasonography at 15 months of age when additional measurements of liveweight, hip height and body condition score were recorded. Bull breeding soundness was evaluated at similar to 12, 18 and 24 months of age when measurements of scrotal circumference, sheath score, semen mass activity, progressive motility of individual sperm and percent morphologically normal sperm were recorded. Magnitude of heritability and genetic correlations changed across time for some traits. Heritability of LH, inhibin, IGF-I and of 18-month scrotal circumference, mass activity, progressive motility and percent normal sperm was 0.31, 0.74, 0.44, 0.75, 0.24, 0.15 and 0.25, respectively, for Brahmans and 0.48, 0.72, 0.36, 0.43, 0.13, 0.15 and 0.20, respectively, for Tropical Composites. Inhibin and IGF-I had moderate genetic association with percent normal sperm at 24 months in Brahmans but low to negligible associations in Tropical Composites. Body condition score in Brahmans and sperm motility (mass and individual) traits in both genotypes had moderate to strong genetic correlation with percent normal sperm and may prove useful candidates for indirect selection. There is scope to increase scrotal circumference by selection and this will be associated with favourable correlated responses of improved semen quality in both genotypes. The lack of genetic antagonism among bull traits indicates that selection for improved semen quality will not adversely affect other production traits.

The inferior whorl for detecting diabetic peripheral neuropathy using corneal confocal microscopy

PURPOSE: In vivo corneal confocal microscopy (CCM) is increasingly used as a surrogate endpoint in studies of diabetic polyneuropathy (DPN). However, it is not clear whether imaging the central cornea provides optimal diagnostic utility for DPN. Therefore, we compared nerve morphology in the central cornea and the inferior whorl, a more distal and densely innervated area located inferior and nasal to the central cornea. METHODS: A total of 53 subjects with type 1/type 2 diabetes and 15 age-matched control subjects underwent detailed assessment of neuropathic symptoms (NPS), deficits (neuropathy disability score [NDS]), quantitative sensory testing (vibration perception threshold [VPT], cold and warm threshold [CT/WT], and cold- and heat-induced pain [CIP/HIP]), and electrophysiology (sural and peroneal nerve conduction velocity [SSNCV/PMNCV], and sural and peroneal nerve amplitude [SSNA/PMNA]) to diagnose patients with (DPN+) and without (DPN-) neuropathy. Corneal nerve fiber density (CNFD) and length (CNFL) in the central cornea, and inferior whorl length (IWL) were quantified. RESULTS: Comparing control subjects to DPN- and DPN+ patients, there was a significant increase in NDS (0 vs. 2.6 ± 2.3 vs. 3.3 ± 2.7, P < 0.01), VPT (V; 5.4 ± 3.0 vs. 10.6 ± 10.3 vs. 17.7 ± 11.8, P < 0.01), WT (°C; 37.7 ± 3.5 vs. 39.1 ± 5.1 vs. 41.7 ± 4.7, P < 0.05), and a significant decrease in SSNCV (m/s; 50.2 ± 5.4 vs. 48.4 ± 5.0 vs. 39.5 ± 10.6, P < 0.05), CNFD (fibers/mm2; 37.8 ± 4.9 vs. 29.7 ± 7.7 vs. 27.1 ± 9.9, P < 0.01), CNFL (mm/mm2; 27.5 ± 3.6 vs. 24.4 ± 7.8 vs. 20.7 ± 7.1, P < 0.01), and IWL (mm/mm2; 35.1 ± 6.5 vs. 26.2 ± 10.5 vs. 23.6 ± 11.4, P < 0.05). For the diagnosis of DPN, CNFD, CNFL, and IWL achieved an area under the curve (AUC) of 0.75, 0.74, and 0.70, respectively, and a combination of IWL-CNFD achieved an AUC of 0.76. CONCLUSIONS: The parameters of CNFD, CNFL, and IWL have a comparable ability to diagnose patients with DPN. However, IWL detects an abnormality even in patients without DPN. Combining IWL with CNFD may improve the diagnostic performance of CCM.

On the propagation of a weak shock front: Theory and application

In this paper the kinematics of a weak shock front governed by a hyperbolic system of conservation laws is studied. This is used to develop a method for solving problems, involving the propagation of nonlinear unimodal waves. It consists of first solving the nonlinear wave problem by moving along the bicharacteristics of the system and then fitting the shock into this solution field, so that it satisfies the necessary jump conditions. The kinematics of the shock leads in a natural way to the definition of ldquoshock-raysrdquo, which play the same role as the ldquoraysrdquo in a continuous flow. A special case of a circular cylinder introduced suddenly in a constant streaming flow is studied in detail. The shock fitted in the upstream region propagates with a velocity which is the mean of the velocities of the linear and the nonlinear wave fronts. In the downstream the solution is given by an expansion wave.

On the propagation of a multi-dimensional shock of arbitrary strength

In this paper the kinematics of a curved shock of arbitrary strength has been discussed using the theory of generalised functions. This is the extension of Moslov’s work where he has considered isentropic flow even across the shock. The condition for a nontrivial jump in the flow variables gives the shock manifold equation (sme). An equation for the rate of change of shock strength along the shock rays (defined as the characteristics of the sme) has been obtained. This exact result is then compared with the approximate result of shock dynamics derived by Whitham. The comparison shows that the approximate equations of shock dynamics deviate considerably from the exact equations derived here. In the last section we have derived the conservation form of our shock dynamic equations. These conservation forms would be very useful in numerical computations as it would allow us to derive difference schemes for which it would not be necessary to fit the shock-shock explicitly.

Multi-trait assessment of early-in-life female, male and genomic measures for use in genetic selection to improve female reproductive performance of Brahman cattle

Early-in-life female and male measures with potential to be practical genetic indicators were chosen from earlier analyses and examined together with genomic measures for multi-trait use to improve female reproduction of Brahman cattle. Combinations of measures were evaluated on the genetic gains expected from selection of sires and dams for each of age at puberty (AGECL, i.e. first observation of a corpus luteum), lactation anoestrous interval in 3-year-old cows (LAI), and lifetime annual weaning rate (LAWR, i.e. the weaning rate of cows based on the number of annual matings they experienced over six possible matings). Selection was on an index of comparable records for each combination. Selection intensities were less than theoretically possible but assumed a concerted selection effort was able to be made across the Brahman breed. The results suggested that substantial genetic gains could be possible but need to be confirmed in other data. The estimated increase in LAWR in 10 years, for combinations without or with genomic measures, ranged from 8 to 12 calves weaned per 100 cows from selection of sires, and from 12 to 15 calves weaned per 100 cows from selection of sires and dams. Corresponding reductions in LAI were 60-103 days or 94-136 days, and those for AGECL were 95-125 or 141-176 days, respectively. Coat score (a measure of the sleekness or wooliness of the coat) and hip height in females, and preputial eversion and liveweight in males, were measures that may warrant wider recording for Brahman female reproduction genetic evaluation. Pregnancy-test outcomes from Matings 1 and 2 also should be recorded. Percentage normal sperm may be important to record for reducing LAI and scrotal size and serum insulin-like growth factor-I concentration in heifers at 18 months for reducing AGECL. Use of a genomic estimated breeding value (EBV) in combination with other measures added to genetic gains, especially at genomic EBV accuracies of 40%. Accuracies of genomic EBVs needed to approach 60% for the genomic EBV to be the most important contributor to gains in the combinations of measures studied.

Reliability and validity of objective measures of physical activity in youth with cerebral palsy who are ambulatory

BACKGROUND Physical therapy for youth with cerebral palsy (CP) who are ambulatory includes interventions to increase functional mobility and participation in physical activity (PA). Thus, reliable and valid measures are needed to document PA in youth with CP. OBJECTIVE The purpose of this study was to evaluate the inter-instrument reliability and concurrent validity of 3 accelerometer-based motion sensors with indirect calorimetry as the criterion for measuring PA intensity in youth with CP. METHODS Fifty-seven youth with CP (mean age=12.5 years, SD=3.3; 51% female; 49.1% with spastic hemiplegia) participated. Inclusion criteria were: aged 6 to 20 years, ambulatory, Gross Motor Function Classification System (GMFCS) levels I through III, able to follow directions, and able to complete the full PA protocol. Protocol activities included standardized activity trials with increasing PA intensity (resting, writing, household chores, active video games, and walking at 3 self-selected speeds), as measured by weight-relative oxygen uptake (in mL/kg/min). During each trial, participants wore bilateral accelerometers on the upper arms, waist/hip, and ankle and a portable indirect calorimeter. Intraclass coefficient correlations (ICCs) were calculated to evaluate inter-instrument reliability (left-to-right accelerometer placement). Spearman correlations were used to examine concurrent validity between accelerometer output (activity and step counts) and indirect calorimetry. Friedman analyses of variance with post hoc pair-wise analyses were conducted to examine the validity of accelerometers to discriminate PA intensity across activity trials. RESULTS All accelerometers exhibited excellent inter-instrument reliability (ICC=.94-.99) and good concurrent validity (rho=.70-.85). All accelerometers discriminated PA intensity across most activity trials. LIMITATIONS This PA protocol consisted of controlled activity trials. CONCLUSIONS Accelerometers provide valid and reliable measures of PA intensity among youth with CP.

Portrait of Benjamin Liebermann

This large ambitious painting is one of Liebermann's earliest. His uncle, the Geheime Kommerzienrat (Privy Councillor) Benjamin Liebermann (1847-1935), is elderly and slim, but erect and proudly wearing his honorary business medals, has his left hand on his hip, while his right hand rests on the back of a red velvet chair. He is wearing long sideburns and whiskers. The background is plain. The style is tight and painstaking, but already exhibits Liebermann's preference for visible, impasto brushstrokes. The emphasis is on the face, with no attempt at idealization of the akward, but strong features. The colors are mainly black and brown.

Launch Envelope Optimization of Virtual Sliding Target Guidance Scheme

This paper presents an optimization of the performance of a recently proposed virtual sliding target (VST) guidance scheme in terms of maximization of its launch envelope for three- dimensional (3-D) engagements. The objective is to obtain the launch envelope of the missile using the VST guidance scheme for different lateral launch angles with respect to the line of sight (LOS) and demonstrate its superiority over kinematics-based guidance laws like proportional navigation (PN). The VST scheme uses PN as its basic guidance scheme and exploits the relation between the atmospheric properties, missile aerodynamic characteristics, and the optimal trajectory of the missile. The missile trajectory is shaped by controlling the instantaneous position and the speed of a virtual target which the missile pursues during the midcourse phase. In the proposed method it is shown that an appropriate value of initial position for the virtual target in 3-D, combined with optimized virtual target parameters, can significantly improve the launch envelope performance. The paper presents the formulation of the optimization problem, obtains the approximate models used to make the optimization problem more tractable, and finally presents the optimized performance of the missile in terms of launch envelope and shows significant improvement over kinematic-based guidance laws. The paper also proposes modification to the basic VST scheme. Some simulations using the full-fledged six degrees-of-freedom (6-DOF) models are also presented to validate the models and technique used.

Identification of IDUA and WNT16 phosphorylation-related non-synonymous polymorphisms for bone mineral density in meta-analyses of genome-wide association studies

Protein phosphorylation regulates a wide variety of cellular processes. Thus, we hypothesize that single-nucleotide polymorphisms (SNPs) that may modulate protein phosphorylation could affect osteoporosis risk. Based on a previous conventional genome-wide association (GWA) study, we conducted a three-stage meta-analysis targeting phosphorylation-related SNPs (phosSNPs) for femoral neck (FN)-bone mineral density (BMD), total hip (HIP)-BMD, and lumbar spine (LS)-BMD phenotypes. In stage 1, 9593 phosSNPs were meta-analyzed in 11,140 individuals of various ancestries. Genome-wide significance (GWS) and suggestive significance were defined by α = 5.21 × 10–6 (0.05/9593) and 1.00 × 10–4, respectively. In stage 2, nine stage 1–discovered phosSNPs (based on α = 1.00 × 10–4) were in silico meta-analyzed in Dutch, Korean, and Australian cohorts. In stage 3, four phosSNPs that replicated in stage 2 (based on α = 5.56 × 10–3, 0.05/9) were de novo genotyped in two independent cohorts. IDUA rs3755955 and rs6831280, and WNT16 rs2707466 were associated with BMD phenotypes in each respective stage, and in three stages combined, achieving GWS for both FN-BMD (p = 8.36 × 10–10, p = 5.26 × 10–10, and p = 3.01 × 10–10, respectively) and HIP-BMD (p = 3.26 × 10–6, p = 1.97 × 10–6, and p = 1.63 × 10–12, respectively). Although in vitro studies demonstrated no differences in expressions of wild-type and mutant forms of IDUA and WNT16B proteins, in silico analyses predicts that WNT16 rs2707466 directly abolishes a phosphorylation site, which could cause a deleterious effect on WNT16 protein, and that IDUA phosSNPs rs3755955 and rs6831280 could exert indirect effects on nearby phosphorylation sites. Further studies will be required to determine the detailed and specific molecular effects of these BMD-associated non-synonymous variants. © 2015 American Society for Bone and Mineral Research.