Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies. © 2011 Duncan et al.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ~2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2×10-9). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3×10-12, and -0.16 SD per G allele, P = 1.2×10-15, respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3×10-9), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9×10-6 and rs2707466: OR = 1.22, P = 7.2×10-6). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16-/- mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5×10-13<P<5.9×10-4) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture. © 2012 Zheng et al.

Mimicry of 21. 5 KDa myelin basic protein peptide by a Maedi Visna virus polymerase peptide does not contribute to the pathogenesis of encephalitis in sheep

Epitope mimicry is the theory that an infectious agent such as a virus causes pathological effects via mimicry of host proteins and thus elicits a cross-reactive immune response to host tissues. Weise and Carnegie (1988) found a region of sequence similarity between the pol gene of the Maedi Visna virus (MVV), which induces demyelinating encephalitis in sheep, and myelin basic protein (MBP), which is known to induce experimental allergic encephalitis (EAE) in laboratory animals. In this study, cross-reactions between sera raised in sheep against synthetic peptides of MVV (TGKIPWILLPGR) and 21.5 kDa MBP (SGKVPWLKRPGR) were demonstrated using enzyme-linked immunosorbant assay (ELISA) and thin layer chromatography (TLC) immunoprobing. The antibody responses of MVV-infected sheep were investigated using ELISA against the peptides, and MBP protein, immunoprobing of the peptides on TPC plates and Western blotting against MBP. Slight significant reactions to the 21.5 kDa MBP peptide (P < 0.001) and to a lesser extent sheep MBP (P < 0.004) were detected in ELISA. The MBP peptide evoked stronger responses from more sera than the MVV peptide on immunoprobed TLC plates. On the Western blots, eight of the 23 sheep with Visna had serum reactivity to MBP. This slight reaction to MBP in MVV-infected sheep is of interest because of the immune responses to MBP evident in multiple sclerosis and EAE, but its relevance in Visna is limited since no correlation with disease severity was observed. The cell-mediated immune responses of MVV-infected sheep against similar peptides was assessed. The peptides did not stimulate proliferation of peripheral blood lymphocytes of MVV-infected sheep. Since the MVV peptide was not recognised by antibodies or T lymphocytes from MVV-infected and encephalic sheep, it was concluded that epitope mimicry of this 21.5 kDa MBP peptide by the similar MVV pol peptide was not contributing to the immunopathogenesis of Visna. The slight antibody response to MBP and the MBP peptide can be attributed to by-stander effects of the immunopathology of MVV-induced encephalitis.

Mimicry between HTLV-I and myelin basic protein: No response in HTLV-I-associated myelopathy patients

The reactivity to a peptide from the HTLV-I polyprotein (FKLPGLNSR) and a similar sequence from myelin basic protein (MBP) (FKLGGRDSR) was examined in relation to the proposal that mimicry of MBP by HTLV-I could be involved in autoimmune responses in HTLV-I-associated myelopathy (HAM). It was found that rabbit antibodies raised against the HTLV-I peptide recognised both peptides, with a titre of 1/10240 to the HTLV-I peptide and 1/5220 to the MBP peptide. Human sera from HAM patients and a HTLV-I carrier without HAM showed slightly higher responses to the HTLV-I peptide compared to the responses from uninfected human sera. HAM patients had greater responses to the HTLV-I peptide than to the similar MBP peptide and an unrelated bovine MBP peptide. There was no recognition of the peptides by peripheral blood lymphocytes from HAM patients or a HTLV-I carrier without HAM. It was concluded that although cross-reactivity was demonstrated in rabbits and the HTLV-I peptide was recognised by sera from HAM patients, the epitope does not appear to evoke a mimicking response to the similar region in MBP. Hence it is not likely to be involved in the pathogenesis of HAM through molecular mimicry.

The effect of LRP5 polymorphisms on bone mineral density is apparent in childhood

Bone mass acquired during childhood is the primary determinant of adult bone mineral density (BMD) and osteoporosis risk. Bone accrual is subject to genetic influences. Activating and inactivating LRP5 gene mutations elicit extreme bone phenotypes, while more common LRP5 polymorphisms are associated with normal variation of BMD. Our aim was to test the hypothesis that LRP5 gene polymorphisms influence bone mass acquisition during childhood. The association between LRP5 gene polymorphisms and bone size and mineralization was examined in 819 unrelated British Caucasian children (n = 429 boys) aged 9 years. Height, weight, pubertal status (where available), total-body and spinal bone area, bone mineral content (BMC), BMD, and area-adjusted BMC (aBMC) were assessed. Dual-energy X-ray absorptiometry (DXA)-gene associations were assessed by linear regression, with adjustment for age, gender, pubertal status, and body size parameters. There were 140, 79, 12, and 2 girls who achieved Tanner stages I-IV, respectively, and 179 and 32 boys who achieved Tanner stages I and II, respectively. The rs2306862 (N740N) coding polymorphism in exon 10 of the LRP5 gene was associated with spinal BMD and aBMC (each P = 0.01) and total-body BMD and aBMC (P = 0.04 and 0.03, respectively). Adjusting for pubertal stage strengthened associations between this polymorphism and spinal BMD and aBMC (P = 0.01 and 0.002, respectively). Individuals homozygous for the T allele had greater spinal BMD and aBMC scores than those homozygous for the C allele. A dose effect was apparent as the mean spinal BMD and aBMC of heterozygous TC individuals were intermediate between those of their TT and CC counterparts. The N740N polymorphism in exon 10 of LRP5 was associated with spinal BMD and aBMC in pre- and early pubertal children. These results indicate that LRP5 influences volumetric bone density in childhood, possibly through effects on trabecular bone formation.

Suggestive linkage of 2p22-25 and 11q12-13 with low bone mineral density at the lumbar spine in the Irish population

Osteoporosis is a disease characterized by low bone mineral density (BMD) and poor bone quality. Peak bone density is achieved by the third decade of life, after which bone is maintained by a balanced cycle of bone resorption and synthesis. Age-related bone loss occurs as the bone resorption phase outweighs the bone synthesis phase of bone metabolism. Heritability accounts for up to 90% of the variability in BMD. Chromosomal loci including 1p36, 2p22-25, 11q12-13, parathyroid hormone receptor type 1 (PTHR1), interleukin-6 (IL-6), interleukin 1 alpha (IL-1α) and type II collagen A1/vitamin D receptor (COL11A1/VDR) have been linked or shown suggestive linkage with BMD in other populations. To determine whether these loci predispose to low BMD in the Irish population, we investigated 24 microsatellite markers at 7 chromosomal loci by linkage studies in 175 Irish families of probands with primary low BMD (T-score ≤ -1.5). Nonparametric analysis was performed using the maximum likelihood variance estimation and traditional Haseman-Elston tests on the Mapmaker/Sibs program. Suggestive evidence of linkage was observed with lumbar spine BMD at 2p22-25 (maximum LOD score 2.76) and 11q12-13 (MLS 2.55). One region, 1p36, approached suggestive linkage with femoral neck BMD (MLS 2.17). In addition, seven markers achieved LOD scores > 1.0, D2S149, D11S1313, D11S987, D11S1314 including those encompassing the PTHR1 (D3S3559, D3S1289) for lumbar spine BMD and D2S149 for femoral neck BMD. Our data suggest that genes within a these chromosomal regions are contributing to a predisposition to low BMD in the Irish population.

Dependence of light attenuation and backscattering on collagen concentration and chondrocyte density in agarose scaffolds

Optical coherence tomography (OCT) has been applied for high resolution imaging of articular cartilage. However, the contribution of individual structural elements of cartilage on OCT signal has not been thoroughly studied. We hypothesize that both collagen and chondrocytes, essential structural components of cartilage, act as important light scatterers and that variation in their concentrations can be detected by OCT through changes in backscattering and attenuation. To evaluate this hypothesis, we established a controlled model system using agarose scaffolds embedded with variable collagen concentrations and chondrocyte densities. Using OCT, we measured the backscattering coefficient (µb) and total attenuation coefficient (µt) in these scaffolds. Along our hypothesis, light backscattering and attenuation in agarose were dependent on collagen concentration and chondrocyte density. Significant correlations were found between µt and chondrocyte density (ρ = 0.853, p < 0.001) and between µt and collagen concentration (ρ = 0.694, p < 0.001). µb correlated significantly with chondrocyte density (ρ = 0.504, p < 0.001) but not with collagen concentration (ρ = 0.103, p = 0.422) of the scaffold. Thus, quantitation of light backscattering and, especially, attenuation could be valuable when evaluating the integrity of soft tissues, such as articular cartilage with OCT.

A suitable method for journal bearing wear measurement

Purpose – Ideally, there is no wear in hydrodynamic lubrication regime. A small amount of wear occurs during start and stop of the machines and the amount of wear is so small that it is difficult to measure with accuracy. Various wear measuring techniques have been used where out-of-roundness was found to be the most reliable method of measuring small wear quantities in journal bearings. This technique was further developed to achieve higher accuracy in measuring small wear quantities. The method proved to be reliable as well as inexpensive. The paper aims to discuss these issues. Design/methodology/approach – In an experimental study, the effect of antiwear additives was studied on journal bearings lubricated with oil containing solid contaminants. The test duration was too long and the wear quantities achieved were too small. To minimise the test duration, short tests of about 90 min duration were conducted and wear was measured recording changes in variety of parameters related to weight, geometry and wear debris. The out-of-roundness was found to be the most effective method. This method was further refined by enlarging the out-of-roundness traces on a photocopier. The method was proved to be reliable and inexpensive. Findings – Study revealed that the most commonly used wear measurement techniques such as weight loss, roughness changes and change in particle count were not adequate for measuring small wear quantities in journal bearings. Out-of-roundness method with some refinements was found to be one of the most reliable methods for measuring small wear quantities in journal bearings working in hydrodynamic lubrication regime. By enlarging the out-of-roundness traces and determining the worn area of the bearing cross-section, weight loss in bearings was calculated, which was repeatable and reliable. Research limitations/implications – This research is a basic in nature where a rudimentary solution has been developed for measuring small wear quantities in rotary devices such as journal bearings. The method requires enlarging traces on a photocopier and determining the shape of the worn area on an out-of-roundness trace on a transparency, which is a simple but a crude method. This may require an automated procedure to determine the weight loss from the out-of-roundness traces directly. This method can be very useful in reducing test duration and measuring wear quantities with higher precision in situations where wear quantities are very small. Practical implications – This research provides a reliable method of measuring wear of circular geometry. The Talyrond equipment used for measuring the change in out-of-roundness due to wear of bearings indicates that this equipment has high potential to be used as a wear measuring device also. Measurement of weight loss from the traces is an enhanced capability of this equipment and this research may lead to the development of a modified version of Talyrond type of equipment for wear measurements in circular machine components. Originality/value – Wear measurement in hydrodynamic bearings requires long duration tests to achieve adequate wear quantities. Out-of-roundness is one of the geometrical parameters that changes with progression of wear in a circular shape components. Thus, out-of-roundness is found to be an effective wear measuring parameter that relates to change in geometry. Method of increasing the sensitivity and enlargement of out-of-roundness traces is original work through which area of worn cross-section can be determined and weight loss can be derived for materials of known density with higher precision.

Effects of Tamoxifen and oestrogen on histology and radiographic density in high and low mammographic density human breast tissues maintained in murine tissue engineering chambers

Mammographic density (MD) is a strong risk factor for breast cancer. It is altered by exogenous endocrine treatments, including hormone replacement therapy and Tamoxifen. Such agents also modify breast cancer (BC) risk. However, the biomolecular basis of how systemic endocrine therapy modifies MD and MD-associated BC risk is poorly understood. This study aims to determine whether our xenograft biochamber model can be used to study the effectiveness of therapies aimed at modulating MD, by examine the effects of Tamoxifen and oestrogen on histologic and radiographic changes in high and low MD tissues maintained within the biochamber model. High and low MD human tissues were precisely sampled under radiographic guidance from prophylactic mastectomy fresh specimens of high-risk women, then inserted into separate vascularized murine biochambers. The murine hosts were concurrently implanted with Tamoxifen, oestrogen or placebo pellets, and the high and low MD biochamber tissues maintained in the murine host environment for 3 months, before the high and low MD biochamber tissues were harvested for histologic and radiographic analyses. The radiographic density of high MD tissue maintained in murine biochambers was decreased in Tamoxifen-treated mice compared to oestrogen-treated mice (p = 0.02). Tamoxifen treatment of high MD tissue in SCID mice led to a decrease in stromal (p = 0.009), and an increase in adipose (p = 0.023) percent areas, compared to placebo-treated mice. No histologic or radiographic differences were observed in low MD biochamber tissue with any treatment. High MD biochamber tissues maintained in mice implanted with Tamoxifen, oestrogen or placebo pellets had dynamic and measurable histologic compositional and radiographic changes. This further validates the dynamic nature of the MD xenograft model, and suggests the biochamber model may be useful for assessing the underlying molecular pathways of Tamoxifen-reduced MD, and in testing of other pharmacologic interventions in a preclinical model of high MD.

Alteration of the C-terminal ligand specificity of the Erbin PDZ domain by allosteric mutational effects

Modulation of protein binding specificity is important for basic biology and for applied science. Here we explore how binding specificity is conveyed in PDZ (postsynaptic density protein-95/discs large/zonula occludens-1) domains, small interaction modules that recognize various proteins by binding to an extended C terminus. Our goal was to engineer variants of the Erbin PDZ domain with altered specificity for the most C-terminal position (position 0) where a Val is strongly preferred by the wild-type domain. We constructed a library of PDZ domains by randomizing residues in direct contact with position 0 and in a loop that is close to but does not contact position 0. We used phage display to select for PDZ variants that bind to 19 peptide ligands differing only at position 0. To verify that each obtained PDZ domain exhibited the correct binding specificity, we selected peptide ligands for each domain. Despite intensive efforts, we were only able to evolve Erbin PDZ domain variants with selectivity for the aliphatic C-terminal side chains Val, Ile and Leu. Interestingly, many PDZ domains with these three distinct specificities contained identical amino acids at positions that directly contact position 0 but differed in the loop that does not contact position 0. Computational modeling of the selected PDZ domains shows how slight conformational changes in the loop region propagate to the binding site and result in different binding specificities. Our results demonstrate that second-sphere residues could be crucial in determining protein binding specificity.

The dependence of computed tomography number to relative electron density conversion on phantom geometry and its impact on planned dose

A computed tomography number to relative electron density (CT-RED) calibration is performed when commissioning a radiotherapy CT scanner by imaging a calibration phantom with inserts of specified RED and recording the CT number displayed. In this work, CT-RED calibrations were generated using several commercially available phantoms to observe the effect of phantom geometry on conversion to electron density and, ultimately, the dose calculation in a treatment planning system. Using an anthropomorphic phantom as a gold standard, the CT number of a material was found to depend strongly on the amount and type of scattering material surrounding the volume of interest, with the largest variation observed for the highest density material tested, cortical bone. Cortical bone gave a maximum CT number difference of 1,110 when a cylindrical insert of diameter 28 mm scanned free in air was compared to that in the form of a 30 × 30 cm2 slab. The effect of using each CT-RED calibration on planned dose to a patient was quantified using a commercially available treatment planning system. When all calibrations were compared to the anthropomorphic calibration, the largest percentage dose difference was 4.2 % which occurred when the CT-RED calibration curve was acquired with heterogeneity inserts removed from the phantom and scanned free in air. The maximum dose difference observed between two dedicated CT-RED phantoms was ±2.1 %. A phantom that is to be used for CT-RED calibrations must have sufficient water equivalent scattering material surrounding the heterogeneous objects that are to be used for calibration.

Experimental and theoretical characterization of 5,10- diminoporphodimethenes: Dearomatized porphyrinoids from palladium-catalyzed hydrazinations of 5,10-diarylporphyrins

A new dearomatized porphyrinoid, 5,10-diiminoporphodimethene (5,10-DIPD), has been prepared by palladium-catalyzed hydrazination of 5,10-dibromo-15,20-bis(3,5-di-tert-butylphenyl)porphyrin and its nickel(II) complex, by using ethyl and 4-methoxybenzyl carbazates. The oxidative dearomatization of the porphyrin ring occurs in high yield. Further oxidation with 2,3-dichloro-5,6-dicyanobenzoquinone forms the corresponding 5,10-bis(azocarboxylates), thereby restoring the porphyrin aromaticity. The UV/visible spectra of the NiII DIPDs exhibit remarkable redshifts of the lowest-energy bands to 780 nm, and differential pulse voltammetry reveals a contracted electrochemical HOMO–LUMO gap of 1.44 V. Density functional theory (DFT) was used to calculate the optimized geometries and frontier molecular orbitals of model 5,10-DIPD Ni7c and 5,10-bis(azocarboxylate) Ni8c. The conformations of the carbamate groups and the configurations of the CNZ unit were considered in conjunction with the NOESY spectra, to generate the global minimum geometry and two other structures with slightly higher energies. In the absence of solution data regarding conformations, ten possible local minimum conformations were considered for Ni8c. Partition of the porphyrin macrocycle into tri- and monopyrrole fragments in Ni7c and the inclusion of terminal conjugating functional groups generate unique frontier molecular orbital distributions and a HOMO–LUMO transition with a strong element of charge transfer from the monopyrrole ring. Time-dependent DFT calculations were performed for the three lowest-energy structures of Ni7c and Ni8c, and weighting according to their energies allowed the prediction of the electronic spectra. The calculations reproduce the lower-energy regions of the spectra and the overall forms of the spectra with high accuracy, but agreement is not as good in the Soret region below 450 nm.

Cardiometabolic effects of genetic upregulation of the interleukin 1 receptor antagonist: A Mendelian randomisation analysis

Background To investigate potential cardiovascular and other effects of long-term pharmacological interleukin 1 (IL-1) inhibition, we studied genetic variants that produce inhibition of IL-1, a master regulator of inflammation. Methods We created a genetic score combining the effects of alleles of two common variants (rs6743376 and rs1542176) that are located upstream of IL1RN, the gene encoding the IL-1 receptor antagonist (IL-1Ra; an endogenous inhibitor of both IL-1α and IL-1β); both alleles increase soluble IL-1Ra protein concentration. We compared effects on inflammation biomarkers of this genetic score with those of anakinra, the recombinant form of IL-1Ra, which has previously been studied in randomised trials of rheumatoid arthritis and other inflammatory disorders. In primary analyses, we investigated the score in relation to rheumatoid arthritis and four cardiometabolic diseases (type 2 diabetes, coronary heart disease, ischaemic stroke, and abdominal aortic aneurysm; 453 411 total participants). In exploratory analyses, we studied the relation of the score to many disease traits and to 24 other disorders of proposed relevance to IL-1 signalling (746 171 total participants). Findings For each IL1RN minor allele inherited, serum concentrations of IL-1Ra increased by 0·22 SD (95% CI 0·18–0·25; 12·5%; p=9·3 × 10−33), concentrations of interleukin 6 decreased by 0·02 SD (−0·04 to −0·01; −1·7%; p=3·5 × 10−3), and concentrations of C-reactive protein decreased by 0·03 SD (−0·04 to −0·02; −3·4%; p=7·7 × 10−14). We noted the effects of the genetic score on these inflammation biomarkers to be directionally concordant with those of anakinra. The allele count of the genetic score had roughly log-linear, dose-dependent associations with both IL-1Ra concentration and risk of coronary heart disease. For people who carried four IL-1Ra-raising alleles, the odds ratio for coronary heart disease was 1·15 (1·08–1·22; p=1·8 × 10−6) compared with people who carried no IL-1Ra-raising alleles; the per-allele odds ratio for coronary heart disease was 1·03 (1·02–1·04; p=3·9 × 10−10). Per-allele odds ratios were 0·97 (0·95–0·99; p=9·9 × 10−4) for rheumatoid arthritis, 0·99 (0·97–1·01; p=0·47) for type 2 diabetes, 1·00 (0·98–1·02; p=0·92) for ischaemic stroke, and 1·08 (1·04–1·12; p=1·8 × 10−5) for abdominal aortic aneurysm. In exploratory analyses, we observed per-allele increases in concentrations of proatherogenic lipids, including LDL-cholesterol, but no clear evidence of association for blood pressure, glycaemic traits, or any of the 24 other disorders studied. Modelling suggested that the observed increase in LDL-cholesterol could account for about a third of the association observed between the genetic score and increased coronary risk. Interpretation Human genetic data suggest that long-term dual IL-1α/β inhibition could increase cardiovascular risk and, conversely, reduce the risk of development of rheumatoid arthritis. The cardiovascular risk might, in part, be mediated through an increase in proatherogenic lipid concentrations. Funding UK Medical Research Council, British Heart Foundation, UK National Institute for Health Research, National Institute for Health Research Cambridge Biomedical Research Centre, European Research Council, and European Commission Framework Programme 7.

Meta-analysis of genome-wide studies identifies MEF2C SNPs associated with bone mineral density at forearm

Background Forearm fractures affect 1.7 million individuals worldwide each year and most occur earlier in life than hip fractures. While the heritability of forearm bone mineral density (BMD) and fracture is high, their genetic determinants are largely unknown. Aim To identify genetic variants associated with forearm BMD and forearm fractures. Methods BMD at distal radius, measured by dualenergy x-ray absorptiometry, was tested for association with common genetic variants. We conducted a metaanalysis of genome-wide association studies for BMD in 5866 subjects of European descent and then selected the variants for replication in 715 Mexican American samples. Gene-based association was carried out to supplement the single-nucleotide polymorphism (SNP) association test. We then tested the BMD-associated SNPs for association with forearm fracture in 2023 cases and 3740 controls. Results We found that five SNPs in the introns of MEF2C were associated with forearm BMD at a genome-wide significance level (p<5×10-8) in meta-analysis (lead SNP, rs11951031[T] -0.20 SDs per allele, p=9.01×10-9). The gene-based association test suggested an association between MEF2C and forearm BMD ( p=0.003). The association between MEF2C variants and risk of fracture did not achieve statistical significance (SNP rs12521522[A]: OR=1.14 (95% CI 0.92 to 1.35), p=0.14). Meta-analysis also revealed two genome-wide suggestive loci at CTNNA2 and 6q23.2. Conclusions These findings demonstrate that variants at MEF2C were associated with forearm BMD, implicating this gene in the determination of BMD at forearm.

Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci

Ankylosing spondylitis is a common, highly heritable inflammatory arthritis affecting primarily the spine and pelvis. In addition to HLA-B*27 alleles, 12 loci have previously been identified that are associated with ankylosing spondylitis in populations of European ancestry, and 2 associated loci have been identified in Asians. In this study, we used the Illumina Immunochip microarray to perform a case-control association study involving 10,619 individuals with ankylosing spondylitis (cases) and 15,145 controls. We identified 13 new risk loci and 12 additional ankylosing spondylitis-associated haplotypes at 11 loci. Two ankylosing spondylitis-associated regions have now been identified encoding four aminopeptidases that are involved in peptide processing before major histocompatibility complex (MHC) class I presentation. Protective variants at two of these loci are associated both with reduced aminopeptidase function and with MHC class I cell surface expression.