Genetic determinants of heel bone properties: Genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n 5 14 260), velocity of sound (VOS; n 5 15 514) and BMD (n 5 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n 5 11 452) and new genotyping in 15 cohorts (de novo n 5 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 3 108) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 3 1014). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 3 106 also had the expected direction of association with any fracture (P < 0.05), including threeSNPswithP < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, thisGWAstudy reveals the effect of several genescommon to central DXA-derivedBMDand heel ultrasound/DXAmeasures and points to anewgenetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Adverse childhood experiences and the health of university students in eight provinces of Vietnam

Recent systematic reviews have emphasized the need for more research into the health and social impacts of adverse childhood experiences (ACEs) in the Asia-Pacific region. This cross-sectional study was conducted with 2099 young adult students in 8 medical universities throughout Vietnam. An anonymous, self-report questionnaire included the World Health Organization ACE-International Questionnaire and standardized measures of mental and physical health. Three quarters (76%) of the students reported at least one exposure to ACEs; 21% had 4 or more ACEs. The most commonly reported adversities were emotional abuse, physical abuse, and witnessing a household member being treated violently (42.3%, 39.9%, and 34.6%, respectively). Co-occurrence of ACEs had dose–response relationships with poor mental health, suicidal ideation, and low physical health–related quality of life. This first multisite study of ACEs among Vietnamese university students provided evidence that childhood adversity is common and is significantly linked with impaired health and well-being into the early adult years

ANKH and susceptibility to and severity of ankylosing spondylitis

Objective. Unconfirmed reports describe association of ankylosing spondylitis (AS) with several candidate genes including ANKH. Cellular export of inorganic pyrophosphate is regulated by the ANK protein, and mutant mice (ank/ank), which have a premature stop codon in the 3′ end of the ank gene, develop severe ankylosis. We tested the association between single-nucleotide polymorphisms (SNP) in these genes and susceptibility to AS in a population of patients with AS. We investigated the role of these genes in terms of functional (BASFI) and metrological (BASMI) measures, and the association with radiological severity (mSASSS). Methods. Our study was conducted on 355 patients with AS and 95 ethnically matched healthy controls. AS was defined according to the modified New York criteria. Four SNP in ANKH (rs27356, rs26307, rs25957, and rs28006) were genotyped. Association analysis was performed using Cochrane-Armitage and linear regression tests for dichotomous and quantitative variables. Analyses of Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), BASFI, and mSASSS were controlled for sex and disease duration. Results. None of the 4 markers showed significant single-locus disease associations (p > 0.05), suggesting that ANKH was not a major determinant of AS susceptibility in our population. No association was observed between these SNP and age at symptom onset, BASDAI, BASFI, BASMI, or mSASSS. Conclusion. These results confirm data in white Europeans that ANKH is probably not a major determinant of susceptibility to AS. ANKH polymorphisms do not markedly influence AS disease severity, as measured by BASMI and mSASSS. The Journal of Rheumatology

'Sink or swim': An evaluation of the clinical characteristics of individuals with high bone mass

Summary High bone mineral density on routine dual energy X-ray absorptiometry (DXA) may indicate an underlying skeletal dysplasia. Two hundred fifty-eight individuals with unexplained high bone mass (HBM), 236 relatives (41% with HBM) and 58 spouses were studied. Cases could not float, had mandible enlargement, extra bone, broad frames, larger shoe sizes and increased body mass index (BMI). HBM cases may harbour an underlying genetic disorder. Introduction High bone mineral density is a sporadic incidental finding on routine DXA scanning of apparently asymptomatic individuals. Such individuals may have an underlying skeletal dysplasia, as seen in LRP5 mutations. We aimed to characterize unexplained HBM and determine the potential for an underlying skeletal dysplasia. Methods Two hundred fifty-eight individuals with unexplained HBM (defined as L1 Z-score ≥ +3.2 plus total hip Z-score ≥ +1.2, or total hip Z-score ≥ +3.2) were recruited from 15 UK centres, by screening 335,115 DXA scans. Unexplained HBM affected 0.181% of DXA scans. Next 236 relatives were recruited of whom 94 (41%) had HBM (defined as L1 Z-score + total hip Z-score ≥ +3.2). Fifty-eight spouses were also recruited together with the unaffected relatives as controls. Phenotypes of cases and controls, obtained from clinical assessment, were compared using random-effects linear and logistic regression models, clustered by family, adjusted for confounders, including age and sex. Results Individuals with unexplained HBM had an excess of sinking when swimming (7.11 [3.65, 13.84], p < 0.001; adjusted odds ratio with 95% confidence interval shown), mandible enlargement (4.16 [2.34, 7.39], p < 0.001), extra bone at tendon/ligament insertions (2.07 [1.13, 3.78], p = 0.018) and broad frame (3.55 [2.12, 5.95], p < 0.001). HBM cases also had a larger shoe size (mean difference 0.4 [0.1, 0.7] UK sizes, p = 0.009) and increased BMI (mean difference 2.2 [1.3, 3.1] kg/m 2, p < 0.001). Conclusion Individuals with unexplained HBM have an excess of clinical characteristics associated with skeletal dysplasia and their relatives are commonly affected, suggesting many may harbour an underlying genetic disorder affecting bone mass.

Risedronate in adults with osteogenesis imperfecta type I: Increased bone mineral density and decreased bone turnover, but high fracture rate persists

Summary Bisphosphonates can increase bone mineral density (BMD) in children with osteogenesis imperfecta (OI). In this study of adults with OI type I, risedronate increased BMD at lumbar spine (but not total hip) and decreased bone turnover. However, the fracture rate in these patients remained high. Introduction Intravenous bisphosphonates given to children with OI can increase BMD and reduce fracture incidence. Oral and/or intravenous bisphosphonates may have similar effects in adults with OI. We completed an observational study of the effect of risedronate in adults with OI type I. Methods Thirty-two adults (mean age, 39 years) with OI type I were treated with risedronate (total dose, 35 mg weekly) for 24 months. Primary outcome measures were BMD changes at lumbar spine (LS) and total hip (TH). Secondary outcome measures were fracture incidence, bone pain, and change in bone turnover markers (serum procollagen type I aminopropeptide (P1NP) and bone ALP). A meta-analysis of published studies of oral bisphosphonates in adults and children with OI was performed. Results Twenty-seven participants (ten males and seventeen females) completed the study. BMD increased at LS by 3.9% (0.815 vs. 0.846 g/cm 2, p=0.007; mean Z-score, -1.93 vs. -1.58, p=0.002), with no significant change at TH. P1NP fell by 37% (p=0.00041), with no significant change in bone ALP (p=0.15). Bone pain did not change significantly (p=0.6). Fracture incidence remained high, with 25 clinical fractures and 10 major fractures in fourteen participants (0.18 major fractures per person per year), with historical data of 0.12 fractures per person per year. The meta-analysis did not demonstrate a significant difference in fracture incidence in patients with OI treated with oral bisphosphonates. Conclusions Risedronate in adults with OI type I results in modest but significant increases in BMD at LS, and decreased bone turnover. However, this may be insufficient to make a clinically significant difference to fracture incidence.

Identification of IL6R and chromosome 11q13.5 as risk loci for asthma

We aimed to identify novel genetic variants affecting asthma risk, since these might provide novel insights into molecular mechanisms underlying the disease. We did a genome-wide association study (GWAS) in 2669 physician-diagnosed asthmatics and 4528 controls from Australia. Seven loci were prioritised for replication after combining our results with those from the GABRIEL consortium (n=26 475), and these were tested in an additional 25 358 independent samples from four in-silico cohorts. Quantitative multi-marker scores of genetic load were constructed on the basis of results from the GABRIEL study and tested for association with asthma in our Australian GWAS dataset. Two loci were confirmed to associate with asthma risk in the replication cohorts and reached genome-wide significance in the combined analysis of all available studies (n=57 800): rs4129267 (OR 1·09, combined p= 2·4×10-8) in the interleukin-6 receptor (IL6R) gene and rs7130588 (OR 1·09, p=1·8×10-8) on chromosome 11q13.5 near the leucine-rich repeat containing 32 gene (LRRC32, also known as GARP). The 11q13.5 locus was significantly associated with atopic status among asthmatics (OR 1·33, p=7×10-4), suggesting that it is a risk factor for allergic but not non-allergic asthma. Multi-marker association results are consistent with a highly polygenic contribution to asthma risk, including loci with weak effects that might be shared with other immune-related diseases, such as NDFIP1, HLA-B, LPP, and BACH2. The IL6R association further supports the hypothesis that cytokine signalling dysregulation affects asthma risk, and raises the possibility that an IL6R antagonist (tocilizumab) may be effective to treat the disease, perhaps in a genotype-dependent manner. Results for the 11q13.5 locus suggest that it directly increases the risk of allergic sensitisation which, in turn, increases the risk of subsequent development of asthma. Larger or more functionally focused studies are needed to characterise the many loci with modest effects that remain to be identified for asthma. National Health and Medical Research Council of Australia. A full list of funding sources is provided in the webappendix. © 2011 Elsevier Ltd.

A genome-wide association study of sleep habits and insomnia

Several aspects of sleep behavior such as timing, duration and quality have been demonstrated to be heritable. To identify common variants that influence sleep traits in the population, we conducted a genome-wide association study of six sleep phenotypes assessed by questionnaire in a sample of 2,323 individuals from the Australian Twin Registry. Genotyping was performed on the Illumina 317, 370, and 610K arrays and the SNPs in common between platforms were used to impute non-genotyped SNPs. We tested for association with more than 2,000,000 common polymorphisms across the genome. While no SNPs reached the genome-wide significance threshold, we identified a number of associations in plausible candidate genes. Most notably, a group of SNPs in the third intron of the CACNA1C gene ranked as most significant in the analysis of sleep latency (P = 1.3 x 10(-)(6)). We attempted to replicate this association in an independent sample from the Chronogen Consortium (n = 2,034), but found no evidence of association (P = 0.73). We have identified several other suggestive associations that await replication in an independent sample. We did not replicate the results from previous genome-wide analyses of self-reported sleep phenotypes after correction for multiple testing.

Genome-wide association study of major depressive disorder: New results, meta-analysis, and lessons learned

Major depressive disorder (MDD) is a common complex disorder with a partly genetic etiology. We conducted a genome-wide association study of the MDD2000+ sample (2431 cases, 3673 screened controls and >1 M imputed single-nucleotide polymorphisms (SNPs)). No SNPs achieved genome-wide significance either in the MDD2000+ study, or in meta-analysis with two other studies totaling 5763 cases and 6901 controls. These results imply that common variants of intermediate or large effect do not have main effects in the genetic architecture of MDD. Suggestive but notable results were: (a) gene-based tests suggesting roles for adenylate cyclase 3 (ADCY3, 2p23.3) and galanin (GAL, 11q13.3); published functional evidence relates both of these to MDD and serotonergic signaling; (b) support for the bipolar disorder risk variant SNP rs1006737 in CACNA1C (P=0.020, odds ratio=1.10), and; (c) lack of support for rs2251219, a SNP identified in a meta-analysis of affective disorder studies (P=0.51). We estimate that sample sizes 1.8- to 2.4-fold greater are needed for association studies of MDD compared with those for schizophrenia to detect variants that explain the same proportion of total variance in liability. Larger study cohorts characterized for genetic and environmental risk factors accumulated prospectively are likely to be needed to dissect more fully the etiology of MDD.

Genomic inflation factors under polygenic inheritance

Population structure, including population stratification and cryptic relatedness, can cause spurious associations in genome-wide association studies (GWAS). Usually, the scaled median or mean test statistic for association calculated from multiple single-nucleotide-polymorphisms across the genome is used to assess such effects, and 'genomic control' can be applied subsequently to adjust test statistics at individual loci by a genomic inflation factor. Published GWAS have clearly shown that there are many loci underlying genetic variation for a wide range of complex diseases and traits, implying that a substantial proportion of the genome should show inflation of the test statistic. Here, we show by theory, simulation and analysis of data that in the absence of population structure and other technical artefacts, but in the presence of polygenic inheritance, substantial genomic inflation is expected. Its magnitude depends on sample size, heritability, linkage disequilibrium structure and the number of causal variants. Our predictions are consistent with empirical observations on height in independent samples of ~4000 and ~133,000 individuals.

GWAS of butyrylcholinesterase activity identifies four novel loci, independent effects within BCHE and secondary associations with metabolic risk factors

Serum butyrylcholinesterase (BCHE) activity is associated with obesity, blood pressure and biomarkers of cardiovascular and diabetes risk. We have conducted a genome-wide association scan to discover genetic variants affecting BCHE activity, and to clarify whether the associations between BCHE activity and cardiometabolic risk factors are caused by variation in BCHE or whether BCHE variation is secondary to the metabolic abnormalities. We measured serum BCHE in adolescents and adults from three cohorts of Australian twin and family studies. The genotypes from approximately 2.4 million single-nucleotide polymorphisms (SNPs) were available in 8791 participants with BCHE measurements. We detected significant associations with BCHE activity at three independent groups of SNPs at the BCHE locus (P = 5.8 x 10(-262), 7.8 x 10(-47), 2.9 x 10(-12)) and at four other loci: RNPEP (P = 9.4 x 10(-16)), RAPH1-ABI2 (P = 4.1 x 10(-18)), UGT1A1 (P = 4.0 x 10(-8)) and an intergenic region on chromosome 8 (P = 1.4 x 10(-8)). These loci affecting BCHE activity were not associated with metabolic risk factors. On the other hand, SNPs in genes previously associated with metabolic risk had effects on BCHE activity more often than can be explained by chance. In particular, SNPs within FTO and GCKR were associated with BCHE activity, but their effects were partly mediated by body mass index and triglycerides, respectively. We conclude that variation in BCHE activity is due to multiple variants across the spectrum from uncommon/large effect to common/small effect, and partly results from (rather than causes) metabolic abnormalities.

Common SNPs explain a large proportion of the heritability for human height

SNPs discovered by genome-wide association studies (GWASs) account for only a small fraction of the genetic variation of complex traits in human populations. Where is the remaining heritability? We estimated the proportion of variance for human height explained by 294,831 SNPs genotyped on 3,925 unrelated individuals using a linear model analysis, and validated the estimation method with simulations based on the observed genotype data. We show that 45% of variance can be explained by considering all SNPs simultaneously. Thus, most of the heritability is not missing but has not previously been detected because the individual effects are too small to pass stringent significance tests. We provide evidence that the remaining heritability is due to incomplete linkage disequilibrium between causal variants and genotyped SNPs, exacerbated by causal variants having lower minor allele frequency than the SNPs explored to date.

The shared genetics of migraine and anxious depression

OBJECTIVES To investigate: - (1) whether shared genetic factors influence migraine and anxious depression; - (2) whether the genetic architecture of migraine depends on anxious depression; - (3) whether the association between migraine and anxious depression is causal. BACKGROUND Migraine and anxious depression frequently occur together, but little is known about the mechanisms causing this association. METHODS A twin study was conducted to model the genetic architecture of migraine and anxious depression and the covariance between them. Anxious depression was also added to the model as a moderator variable to examine whether anxious depression affects the genetic architecture of migraine. Causal models were explored with the co-twin control method. RESULTS Modest but significant phenotypic (rP=0.28), genetic (rG=0.30), and nonshared environmental (rE=0.26) correlations were found between the 2 traits. Interestingly, the heritability of migraine depended on the level of anxious depression: the higher the anxious depression score, the lower the relative contribution of genetic factors to the individual differences in migraine susceptibility. The observed risk patterns in discordant twins are most consistent with a bidirectional causal relationship. CONCLUSIONS These findings confirm the genetic association between migraine and anxious depression and are consistent with a syndromic association between the 2 traits. This highlights the importance of taking comorbidity into account in genetic studies of migraine, especially in the context of selection for large-scale genotyping efforts. Genetic studies may be most effective when migraine with and without comorbid anxious depression are treated as separate phenotypes.

Common variants in the trichohyalin gene are associated with straight hair in Europeans

Hair morphology is highly differentiated between populations and among people of European ancestry. Whereas hair morphology in East Asian populations has been studied extensively, relatively little is known about the genetics of this trait in Europeans. We performed a genome-wide association scan for hair morphology (straight, wavy, curly) in three Australian samples of European descent. All three samples showed evidence of association implicating the Trichohyalin gene (TCHH), which is expressed in the developing inner root sheath of the hair follicle, and explaining approximately 6% of variance (p=1.5x10(-31)). These variants are at their highest frequency in Northern Europeans, paralleling the distribution of the straight-hair EDAR variant in Asian populations.

A high-density association screen of 155 ion transport genes for involvement with common migraine

The clinical overlap between monogenic Familial Hemiplegic Migraine (FHM) and common migraine subtypes, and the fact that all three FHM genes are involved in the transport of ions, suggest that ion transport genes may underlie susceptibility to common forms of migraine. To test this leading hypothesis, we examined common variation in 155 ion transport genes using 5257 single nucleotide polymorphisms (SNPs) in a Finnish sample of 841 unrelated migraine with aura cases and 884 unrelated non-migraine controls. The top signals were then tested for replication in four independent migraine case-control samples from the Netherlands, Germany and Australia, totalling 2835 unrelated migraine cases and 2740 unrelated controls. SNPs within 12 genes (KCNB2, KCNQ3, CLIC5, ATP2C2, CACNA1E, CACNB2, KCNE2, KCNK12, KCNK2, KCNS3, SCN5A and SCN9A) with promising nominal association (0.00041 < P < 0.005) in the Finnish sample were selected for replication. Although no variant remained significant after adjusting for multiple testing nor produced consistent evidence for association across all cohorts, a significant epistatic interaction between KCNB2 SNP rs1431656 (chromosome 8q13.3) and CACNB2 SNP rs7076100 (chromosome 10p12.33) (pointwise P = 0.00002; global P = 0.02) was observed in the Finnish case-control sample. We conclude that common variants of moderate effect size in ion transport genes do not play a major role in susceptibility to common migraine within these European populations, although there is some evidence for epistatic interaction between potassium and calcium channel genes, KCNB2 and CACNB2. Multiple rare variants or trans-regulatory elements of these genes are not ruled out.

Zygosity diagnosis in the absence of genotypic data: An approach using latent class analysis

For zygosity diagnosis in the absence of genotypic data, or in the recruitment phase of a twin study where only single twins from same-sex pairs are being screened, or to provide a test for sample duplication leading to the false identification of a dizygotic pair as monozygotic, the appropriate analysis of respondents' answers to questions about zygosity is critical. Using data from a young adult Australian twin cohort (N = 2094 complete pairs and 519 singleton twins from same-sex pairs with complete responses to all zygosity items), we show that application of latent class analysis (LCA), fitting a 2-class model, yields results that show good concordance with traditional methods of zygosity diagnosis, but with certain important advantages. These include the ability, in many cases, to assign zygosity with specified probability on the basis of responses of a single informant (advantageous when one zygosity type is being oversampled); and the ability to quantify the probability of misassignment of zygosity, allowing prioritization of cases for genotyping as well as identification of cases of probable laboratory error. Out of 242 twins (from 121 like-sex pairs) where genotypic data were available for zygosity confirmation, only a single case was identified of incorrect zygosity assignment by the latent class algorithm. Zygosity assignment for that single case was identified by the LCA as uncertain (probability of being a monozygotic twin only 76%), and the co-twin's responses clearly identified the pair as dizygotic (probability of being dizygotic 100%). In the absence of genotypic data, or as a safeguard against sample duplication, application of LCA for zygosity assignment or confirmation is strongly recommended.