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.

Utility of temporal artery biopsy samples for genome-wide analysis of giant cell arteritis

Giant Cell Arteritis (GCA) is the most common vasculitis affecting the elderly. Archived formalin-fixed paraffin-embedded (FFPE) temporal artery biopsy (TAB) specimens potentially represent a valuable resource for large-scale genetic analysis of this disease. FFPE TAB samples were obtained from 12 patients with GCA. Extracted TAB DNA was assessed by real time PCR before restoration using the Illumina HD FFPE Restore Kit. Paired FFPE-blood samples were genotyped on the Illumina OmniExpress FFPE microarray. The FFPE samples that passed stringent quality control measures had a mean genotyping success of >97%. When compared with their matching peripheral blood DNA, the mean discordant heterozygote and homozygote single nucleotide polymorphisms calls were 0.0028 and 0.0003, respectively, which is within the accepted tolerance of reproducibility. This work demonstrates that it is possible to successfully obtain high-quality microarray-based genotypes FFPE TAB samples and that this data is similar to that obtained from peripheral blood.

Genetic associations and functional characterization of M1 aminopeptidases and immune-mediated diseases

Endosplasmic reticulum aminopeptidase 1 (ERAP1), endoplasmic reticulum aminopeptidase 2 (ERAP2) and puromycin-sensitive aminopeptidase (NPEPPS) are key zinc metallopeptidases that belong to the oxytocinase subfamily of M1 aminopeptidase family. NPEPPS catalyzes the processing of proteosome-derived peptide repertoire followed by trimming of antigenic peptides by ERAP1 and ERAP2 for presentation on major histocompatibility complex (MHC) Class I molecules. A series of genome-wide association studies have demonstrated associations of these aminopeptidases with a range of immune-mediated diseases such as ankylosing spondylitis, psoriasis, Behçet's disease, inflammatory bowel disease and type I diabetes, and significantly, genetic interaction between some aminopeptidases and HLA Class I loci with which these diseases are strongly associated. In this review, we highlight the current state of understanding of the genetic associations of this class of genes, their functional role in disease, and potential as therapeutic targets.

Genetics of rheumatoid arthritis contributes to biology and drug discovery

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA)1. Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ~10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2, 3, 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation5, cis-acting expression quantitative trait loci6 and pathway analyses7, 8, 9—as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes—to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

Serum vitamin D levels are lower in Australian children and adolescents with type 1 diabetes than in children without diabetes

Vitamin D is synthesised in the skin through the action of UVB radiation (sunlight), and 25-hydroxy vitamin D (25OHD) measured in serum as a marker of vitamin D status. Several studies, mostly conducted in high latitudes, have shown an association between type 1 diabetes mellitus (T1DM) and low serum 25OHD. We conducted a case-control study to determine whether, in a sub-tropical environment with abundant sunlight (latitude 27.5°S), children with T1DM have lower serum vitamin D than children without diabetes. Fifty-six children with T1DM (14 newly diagnosed) and 46 unrelated control children participated in the study. Serum 25OHD, 1,25-dihydroxy vitamin D (1,25(OH)2D) and selected biochemical indices were measured. Vitamin D receptor (VDR) polymorphisms Taq1, Fok1, and Apa1 were genotyped. Fitzpatrick skin classification, self-reported daily hours of outdoor exposure, and mean UV index over the 35d prior to blood collection were recorded. Serum 25OHD was lower in children with T1DM (n=56) than in controls (n=46) [mean (95%CI)=78.7 (71.8-85.6) nmol/L vs. 91.4 (83.5-98.7) nmol/L, p=0.02]. T1DM children had lower self-reported outdoor exposure and mean UV exposure, but no significant difference in distribution of VDR polymorphisms. 25OHD remained lower in children with T1DM after covariate adjustment. Children newly diagnosed with T1DM had lower 1,25(OH)2D [median (IQR)=89 (68-122) pmol/L] than controls [121 (108-159) pmol/L, p=0.03], or children with established diabetes [137 (113-153) pmol/L, p=0.01]. Children with T1DM have lower 25OHD than controls, even in an environment of abundant sunlight. Whether low vitamin D is a risk factor or consequence of T1DM is unknown. © 2012 John Wiley & Sons A/S.

Genome-wide association analysis identifies 11 risk variants associated with the asthma with hay fever phenotype

Background To date, no genome-wide association study (GWAS) has considered the combined phenotype of asthma with hay fever. Previous analyses of family data from the Tasmanian Longitudinal Health Study provide evidence that this phenotype has a stronger genetic cause than asthma without hay fever. Objective We sought to perform a GWAS of asthma with hay fever to identify variants associated with having both diseases. Methods We performed a meta-analysis of GWASs comparing persons with both physician-diagnosed asthma and hay fever (n = 6,685) with persons with neither disease (n = 14,091). Results At genome-wide significance, we identified 11 independent variants associated with the risk of having asthma with hay fever, including 2 associations reaching this level of significance with allergic disease for the first time: ZBTB10 (rs7009110; odds ratio [OR], 1.14; P = 4 × 10−9) and CLEC16A (rs62026376; OR, 1.17; P = 1 × 10−8). The rs62026376:C allele associated with increased asthma with hay fever risk has been found to be associated also with decreased expression of the nearby DEXI gene in monocytes. The 11 variants were associated with the risk of asthma and hay fever separately, but the estimated associations with the individual phenotypes were weaker than with the combined asthma with hay fever phenotype. A variant near LRRC32 was a stronger risk factor for hay fever than for asthma, whereas the reverse was observed for variants in/near GSDMA and TSLP. Single nucleotide polymorphisms with suggestive evidence for association with asthma with hay fever risk included rs41295115 near IL2RA (OR, 1.28; P = 5 × 10−7) and rs76043829 in TNS1 (OR, 1.23; P = 2 × 10−6). Conclusion By focusing on the combined phenotype of asthma with hay fever, variants associated with the risk of allergic disease can be identified with greater efficiency.

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.

Schizophrenia Genetic Variants are not Associated with Intelligence

Background Schizophrenia is associated with lower pre-morbid intelligence (IQ) in addition to (pre-morbid) cognitive decline. Both schizophrenia and IQ are highly heritable traits. Therefore, we hypothesized that genetic variants associated with schizophrenia, including copy number variants (CNVs) and a polygenic schizophrenia (risk) score (PSS), may influence intelligence. Method IQ was estimated with the Wechsler Adult Intelligence Scale (WAIS). CNVs were determined from single nucleotide polymorphism (SNP) data using the QuantiSNP and PennCNV algorithms. For the PSS, odds ratios for genome-wide SNP data were calculated in a sample collected by the Psychiatric Genome-Wide Association Study (GWAS) Consortium (8690 schizophrenia patients and 11 831 controls). These were used to calculate individual PSSs in our independent sample of 350 schizophrenia patients and 322 healthy controls. Results Although significantly more genes were disrupted by deletions in schizophrenia patients compared to controls (p = 0.009), there was no effect of CNV measures on IQ. The PSS was associated with disease status (R 2 = 0.055, p = 2.1 × 10 -7) and with IQ in the entire sample (R 2 = 0.018, p = 0.0008) but the effect on IQ disappeared after correction for disease status. Conclusions Our data suggest that rare and common schizophrenia-associated variants do not explain the variation in IQ in healthy subjects or in schizophrenia patients. Thus, reductions in IQ in schizophrenia patients may be secondary to other processes related to schizophrenia risk. © Cambridge University Press 2013.

Genome-wide association study of intraocular pressure identifies the GLCCI1/ICA1 region as a glaucoma susceptibility locus

To discover quantitative trait loci for intraocular pressure, amajor risk factor for glaucoma and the only modifiable one,weperformed agenome-wide association studyonadiscoverycohort of2175individualsfromSydney, Australia. We found a novel association between intraocular pressure and a common variant at 7p21 near to GLCCI1 and ICA1. The findings in this region were confirmed through two UK replication cohorts totalling 4866 individuals (rs59072263, Pcombined = 1.10 × 10-8). A copy of the G allele at this SNP is associated with an increase in mean IOP of 0.45 mmHg (95%CI = 0.30-0.61 mmHg). These results lend support to the implication of vesicle trafficking and glucocorticoid inducibility pathways in the determination of intraocular pressure and in the pathogenesis of primary open-angle glaucoma.

Genome-wide association analysis identifies 13 new risk loci for schizophrenia

Schizophrenia is an idiopathic mental disorder with a heritable component and a substantial public health impact. We conducted a multi-stage genome-wide association study (GWAS) for schizophrenia beginning with a Swedish national sample (5,001 cases and 6,243 controls) followed by meta-Analysis with previous schizophrenia GWAS (8,832 cases and 12,067 controls) and finally by replication of SNPs in 168 genomic regions in independent samples (7,413 cases, 19,762 controls and 581 parent-offspring trios). We identified 22 loci associated at genome-wide significance; 13 of these are new, and 1 was previously implicated in bipolar disorder. Examination of candidate genes at these loci suggests the involvement of neuronal calcium signaling. We estimate that 8,300 independent, mostly common SNPs (95% credible interval of 6,300-10,200 SNPs) contribute to risk for schizophrenia and that these collectively account for at least 32% of the variance in liability. Common genetic variation has an important role in the etiology of schizophrenia, and larger studies will allow more detailed understanding of this disorder.

A new regulatory variant in the interleukin-6 receptor gene associates with asthma risk

The main genetic determinant of soluble interleukin 6 receptor (sIL-6R) levels is the missense variant rs2228145 that maps to the cleavage site of IL-6R. For each Ala allele, sIL-6R serum levels increase by ∼20 ng ml -1 and asthma risk by 1.09-fold. However, this variant does not explain the total heritability for sIL-6R levels. Additional independent variants in IL6R may therefore contribute to variation in sIL-6R levels and influence asthma risk. We imputed 471 variants in IL6R and tested these for association with sIL-6R serum levels in 360 individuals. An intronic variant (rs12083537) was associated with sIL-6R levels independently of rs4129267 (P=0.0005), a proxy single-nucleotide polymorphism for rs2228145. A significant and consistent association for rs12083537 was observed in a replication panel of 354 individuals (P=0.033). Each rs12083537:A allele increased sIL-6R serum levels by 2.4 ng ml -1. Analysis of mRNA levels in two cohorts did not identify significant associations between rs12083537 and IL6R transcription levels. On the other hand, results from 16 705 asthmatics and 30 809 controls showed that the rs12083537:A allele increased asthma risk by 1.04-fold (P=0.0419). Genetic risk scores based on IL6R regulatory variants may prove useful in explaining variation in clinical response to tocilizumab, an anti-IL-6R monoclonal antibody.

Genetic insights into common pathways and complex relationships among iImmune-mediated diseases

Shared aetiopathogenic factors among immune-mediated diseases have long been suggested by their co-familiality and co-occurrence, and molecular support has been provided by analysis of human leukocyte antigen (HLA) haplotypes and genome-wide association studies. The interrelationships can now be better appreciated following the genotyping of large immune disease sample sets on a shared SNP array: the 'Immunochip'. Here, we systematically analyse loci shared among major immune-mediated diseases. This reveals that several diseases share multiple susceptibility loci, but there are many nuances. The most associated variant at a given locus frequently differs and, even when shared, the same allele often has opposite associations. Interestingly, risk alleles conferring the largest effect sizes are usually disease-specific. These factors help to explain why early evidence of extensive 'sharing' is not always reflected in epidemiological overlap. © 2013 Macmillan Publishers Limited. All rights reserved.

Comparison of the performance of two commercial genome-wide association study genotyping platforms in Han Chinese samples

Most genome-wide association studies to date have been performed in populations of European descent, but there is increasing interest in expanding these studies to other populations. The performance of genotyping chips in Asian populations is not well established. Therefore, we sought to test the performance of widely used fixed-marker, genome-wide association studies chips in the Han Chinese population. Non-HapMap Chinese samples (n = 396) were genotyped using the Illumina OmniExpress and Affymetrix 6.0 platforms, whereas a subset also were genotyped using the Immunochip. Genotyped markers from the Affymetrix 6.0 and Illumina OmniExpress were used for full genome imputation based on the HapMap 2 JPT+CHB (Japanese from Tokyo, Japan and Chinese from Beijing, China) reference panel. The concordance between markers genotypes for the three platforms was very high whether directly genotyped or genotyped and imputed single nucleotide polymorphisms (SNPs; .99.8% for directly genotyped and .99.5% for genotyped and imputed SNPs, respectively) were compared. The OmniExpress chip data enabled more SNPs to be imputed, particularly SNPs with minor allele frequency .5%. The OmniExpress chip achieved better coverage of HapMap SNPs than the Affymetrix 6.0 chip (73.6% vs. 65.9%, respectively, for minor allele frequency .5%). The Affymetrix 6.0 and Illumina OmniExpress chip have similar genotyping accuracy and provide similar accuracy of imputed SNPs. The OmniExpress chip however provides better coverage of Asian HapMap SNPs, although its coverage of HapMap SNPs is moderate. © 2013 Jiang et al.

Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to visceral leishmaniasis

To identify susceptibility loci for visceral leishmaniasis, we undertook genome-wide association studies in two populations: 989 cases and 1,089 controls from India and 357 cases in 308 Brazilian families (1,970 individuals). The HLA-DRB1-HLA-DQA1 locus was the only region to show strong evidence of association in both populations. Replication at this region was undertaken in a second Indian population comprising 941 cases and 990 controls, and combined analysis across the three cohorts for rs9271858 at this locus showed P combined = 2.76 × 10 -17 and odds ratio (OR) = 1.41, 95% confidence interval (CI) = 1.30-1.52. A conditional analysis provided evidence for multiple associations within the HLA-DRB1-HLA-DQA1 region, and a model in which risk differed between three groups of haplotypes better explained the signal and was significant in the Indian discovery and replication cohorts. In conclusion, the HLA-DRB1-HLA-DQA1 HLA class II region contributes to visceral leishmaniasis susceptibility in India and Brazil, suggesting shared genetic risk factors for visceral leishmaniasis that cross the epidemiological divides of geography and parasite species. © 2013 Nature America, Inc. All rights reserved.

Genome-wide association studies

Genome-wide association studies (GWAS) are a powerful hypothesis-free tool for the dissection of susceptibility to common heritable human diseases, including osteoporosis. To date, more than 2000 loci for common human diseases have been identified by GWAS. Success using the GWAS model depends on genetic risk being determined by shared stretches of DNA carried with different frequencies in cases and controls, inherited from ancient ancestors, termed the “common disease–common variant” hypothesis. Not all disease risk is caused by common variants, however, and thus GWAS will not detect all variants involved. Successful GWAS performance requires careful quality control, especially as the effect sizes under study are modest, and there are multiple potential sources of error. Conservative interpretation, use of stringent significance thresholds, and replication in independent cohorts are required to ensure results are robust. Despite these challenging parameters, much has been learnt from GWAS and, as the approach matures and is modified to identify a wider range of variants, significantly more will be learnt about the etiopathogenesis of common diseases such as osteoporosis.