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.

Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1

We report a genome-wide association study for open-angle glaucoma (OAG) blindness using a discovery cohort of 590 individuals with severe visual field loss (cases) and 3,956 controls. We identified associated loci at TMCO1 (rs4656461[G] odds ratio (OR) = 1.68, P = 6.1 × 10-10) and CDKN2B-AS1 (rs4977756[A] OR = 1.50, P = 4.7 × 10-9). We replicated these associations in an independent cohort of cases with advanced OAG (rs4656461 P = 0.010; rs4977756 P = 0.042) and two additional cohorts of less severe OAG (rs4656461 combined discovery and replication P = 6.00 × 10-14, OR = 1.51, 95% CI 1.35-1.68; rs4977756 combined P = 1.35 × 10-14, OR = 1.39, 95% CI 1.28-1.51). We show retinal expression of genes at both loci in human ocular tissues. We also show that CDKN2A and CDKN2B are upregulated in the retina of a rat model of glaucoma. © 2011 Nature America, Inc. All rights reserved.

Genomics of ankylosing spondylitis

Ankylosing spondylitis (AS) is the prototypic and most prevalent and debilitating spondyloarthropathy, a group of arthritides where the spine and pelvis are specifically targeted. Unlike many other forms of arthritis in which joint damage is mediated through tissue destruction, in AS uncontrolled bone formation occurs, frequently resulting in joint fusion and consequently significant disability. It is estimated that there are 2.4 million spondyloarthritis sufferers in the U.S., twice as many as rheumatoid arthritis. The pathogenesis of AS is very poorly understood and both genetics and gene expression profiling approaches have been utilized to elucidate the underlying mechanisms and pathways that drive the disease. Using powerful genome-wide association study approaches a number of candidate genes have been found to be associated with AS. However, although such approaches can identify genes that can contribute to the disease process, they do not inform us of the actual changes in gene/cell activity at any point in the disease process. Expression profiling allows us to take a "snapshot" of cellular activity and what gene activity changes are underlying those changes. A number of expression profiling studies have been undertaken in AS, looking at both circulating cells and tissues from affected joints. The results to date have been somewhat disappointing with little consensus on gene activity changes due to the low power of the studies undertaken. Some more recent better powered studies have identified diagnostic expression profiles that do point to a possible role for expression profiling in early AS diagnosis. Future studies will require collaborative approaches to target specific disease stages and sites with larger numbers of samples.

Epidemiology of ankylosing spondylitis: IGAS 2009

The International Genetics of Ankylosing Spondylitis (IGAS) meeting was held in Houston, Texas, July 25, 2009. Sixteen investigators from Asia, Australia, Europe, and North and South America presented the status of their respective cohorts of patients with ankylosing spondylitis (AS). They also reviewed a proposal to examine their patients by single-nucleotide polymorphism (SNP) genotyping on an Illumina Infinium microarray SNP genotyping chip in a case-control cohort exceeding 12,000 samples. This chip will type 200,000 SNP selected from the most strongly associated variants identified in genome-wide association studies of inflammatory diseases, including inflammatory bowel disease, psoriasis, and ankylosing spondylitis.

The multiple sclerosis whole blood mRNA transcriptome and genetic associations indicate dysregulation of specific T cell pathways in pathogenesis

Multiple sclerosis (MS) is an autoimmune disease with a genetic component, caused at least in part by aberrant lymphocyte activity. The whole blood mRNA transcriptome was measured for 99 untreated MS patients: 43 primary progressive MS, 20 secondary progressive MS, 36 relapsing remitting MS and 45 age-matched healthy controls. The ANZgene Multiple Sclerosis Genetics Consortium genotyped more than 300 000 SNPs for 115 of these samples. Transcription from genes on translational regulation, oxidative phosphorylation, immune synapse and antigen presentation pathways was markedly increased in all forms of MS. Expression of genes tagging T cells was also upregulated (P < 10-12) in MS. A T cell gene signature predicts disease state with a concordance index of 0.79 with age and gender as co-variables, but the signature is not associated with clinical course or disability. The ANZgene genome wide association screen identified two novel regions with genome wide significance: one encoding the T cell co-stimulatory molecule, CD40; the other a region on chromosome 12q13-14. The CD40 haplotype associated with increased MS susceptibility has decreased gene expression in MS (P < 0.0007). The second MS susceptibility region includes 17 genes on 12q13-14 in tight linkage disequilibrium. Of these, only 13 are expressed in leukocytes, and of these the expression of one, FAM119B, is much lower in the susceptibility haplotype (P tdthomlt; 10-14). Overall, these data indicate dysregulation of T cells can be detected in the whole blood of untreated MS patients, and supports targeting of activated T cells in therapy for all forms of MS.

Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10−8). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10−4, Bonferroni corrected), of which six reached P < 5 × 10−8, including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.

Common variants in the region around Osterix are associated with bone mineral density and growth in childhood

Peak bone mass achieved in adolescence is a determinant of bone mass in later life. In order to identify genetic variants affecting bone mineral density (BMD), we performed a genome-wide association study of BMD and related traits in 1518 children from the Avon Longitudinal Study of Parents and Children (ALSPAC). We compared results with a scan of 134 adults with high or low hip BMD. We identified associations with BMD in an area of chromosome 12 containing the Osterix (SP7) locus, a transcription factor responsible for regulating osteoblast differentiation (ALSPAC: P = 5.8 × 10-4; Australia: P = 3.7 × 10-4). This region has previously shown evidence of association with adult hip and lumbar spine BMD in an Icelandic population, as well as nominal association in a UK population. A meta-analysis of these existing studies revealed strong association between SNPs in the Osterix region and adult lumbar spine BMD (P = 9.9 × 10-11). In light of these findings, we genotyped a further 3692 individuals from ALSPAC who had whole body BMD and confirmed the association in children as well (P = 5.4 × 10-5). Moreover, all SNPs were related to height in ALSPAC children, but not weight or body mass index, and when height was included as a covariate in the regression equation, the association with total body BMD was attenuated. We conclude that genetic variants in the region of Osterix are associated with BMD in children and adults probably through primary effects on growth.

Genome-wide association study identifies eight loci associated with blood pressure

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5 million genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N 71,225 European ancestry, N 12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N = 29,136). We identified association between systolic or diastolic blood pressure and common variants in eight regions near the CYP17A1 (P = 7 × 10 24), CYP1A2 (P = 1 × 10 23), FGF5 (P = 1 × 10 21), SH2B3 (P = 3 × 10 18), MTHFR (P = 2 × 10 13), c10orf107 (P = 1 × 10 9), ZNF652 (P = 5 × 10 9) and PLCD3 (P = 1 × 10 8) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.

Non-B27 MHC associations of ankylosing spondylitis

Ankylosing spondylitis (AS) has been associated with human leukocyte antigen (HLA)-B27 for over 30 years; however, the mechanism of action has remained elusive. Although many studies have reported associations between AS and other genes in the major histocompatibility complex (MHC) in AS, no conclusive results have emerged. To investigate the contribution of non-B27 MHC genes to AS, a large cohort of AS families and controls were B27 typed and genotyped across the region. Interrogation of the data identified a region of 270kb, lying from 31952649 to 32221738 base pairs from the p-telomere of chromosome 6 and containing 23 genes, which is likely to include genes involved with susceptibility to AS.

Combined analysis of three whole genome linkage scans for Ankylosing Spondylitis

Objective. Ankylosing spondylitis (AS) is a debilitating chronic inflammatory condition with a high degree of familiality (λs=82) and heritability (>90%) that primarily affects spinal and sacroiliac joints. Whole genome scans for linkage to AS phenotypes have been conducted, although results have been inconsistent between studies and all have had modest sample sizes. One potential solution to these issues is to combine data from multiple studies in a retrospective meta-analysis. Methods: The International Genetics of Ankylosing Spondylitis Consortium combined data from three whole genome linkage scans for AS (n=3744 subjects) to determine chromosomal markers that show evidence of linkage with disease. Linkage markers typed in different centres were integrated into a consensus map to facilitate effective data pooling. We performed a weighted meta-analysis to combine the linkage results, and compared them with the three individual scans and a combined pooled scan. Results: In addition to the expected region surrounding the HLA-B27 gene on chromosome 6, we determined that several marker regions showed significant evidence of linkage with disease status. Regions on chromosome 10q and 16q achieved 'suggestive' evidence of linkage, and regions on chromosomes 1q, 3q, 5q, 6q, 9q, 17q and 19q showed at least nominal linkage in two or more scans and in the weighted meta-analysis. Regions previously associated with AS on chromosome 2q (the IL-1 gene cluster) and 22q (CYP2D6) exhibited nominal linkage in the meta-analysis, providing further statistical support for their involvement in susceptibility to AS. Conclusion: These findings provide a useful guide for future studies aiming to identify the genes involved in this highly heritable condition. . Published by on behalf of the British Society for Rheumatology.

The v-MFG test: Investigating maternal, offspring and maternal-fetal genetic incompatibility effects on disease and viability

The MFG test is a family-based association test that detects genetic effects contributing to disease in offspring, including offspring allelic effects, maternal allelic effects and MFG incompatibility effects. Like many other family-based association tests, it assumes that the offspring survival and the offspring-parent genotypes are conditionally independent provided the offspring is affected. However, when the putative disease-increasing locus can affect another competing phenotype, for example, offspring viability, the conditional independence assumption fails and these tests could lead to incorrect conclusions regarding the role of the gene in disease. We propose the v-MFG test to adjust for the genetic effects on one phenotype, e.g., viability, when testing the effects of that locus on another phenotype, e.g., disease. Using genotype data from nuclear families containing parents and at least one affected offspring, the v-MFG test models the distribution of family genotypes conditional on offspring phenotypes. It simultaneously estimates genetic effects on two phenotypes, viability and disease. Simulations show that the v-MFG test produces accurate genetic effect estimates on disease as well as on viability under several different scenarios. It generates accurate type-I error rates and provides adequate power with moderate sample sizes to detect genetic effects on disease risk when viability is reduced. We demonstrate the v-MFG test with HLA-DRB1 data from study participants with rheumatoid arthritis (RA) and their parents, we show that the v-MFG test successfully detects an MFG incompatibility effect on RA while simultaneously adjusting for a possible viability loss.

Non-inherited maternal HLA alleles are associated with rheumatoid arthritis

Background. Rheumatoid arthritis (RA) is strongly associated with a series of HLA-DRB1 alleles that encode a conserved sequence of amino acids (70Q/R K/R R A A74) in the DRβ1 chain, known as the shared epitope (SE). However 30% of patients are negative for DRB1*04 and 15% are SE-negative. Exposure to these alleles as non-inherited maternal antigens (NIMA) might explain this discrepancy. We undertook a family study to investigate the role of NIMA in RA. Methods. One hundred families, including the RA proband and both parents, were recruited. HLA-DRB1 genotyping was performed using an allele-specific polymerase chain reaction by standard methods. The frequencies of NIMA and non-inherited paternal antigens (NIPA) were compared using contingency tables and a two-tailed P test. We then reviewed four previously published studies of NIMA in RA and conducted an analysis of the combined data Results. We identified 36 families in which the proband was DRB1*04-negative and 13 in which the proband lacked the SE. There was an excess of DRB1*04 and SE NIMA (P=0.05) compared with NIPA. Combined analysis with previous studies showed that 53/231 mothers (23%) versus 25/205 fathers (12%) had a non-inherited DRB1*04 (P=0.003) and 30/99 mothers versus 18/101 fathers had a non-inherited SE allele (P=0.03). Conclusion. A role for HLA NIMA in RA is suggested by these results.

Polymorphisms of the CYP2D6 gene increase susceptibility to ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly familial rheumatic disorder. The sibling recurrence risk ratio for the disease is 63 and heritability assessed in twins > 90%. Although MHC genes, including HLA-B27, contribute only 20-50% of the genetic risk for the disease, no non-MHC gene has yet been convincingly demonstrated to influence either susceptibility to the disease or its phenotypic expression. Previous linkage and association studies have suggested the presence of a susceptibility gene for AS close to, or within, the cytochrome P450 2D6 gene (CYP2D6, debrisoquine hydroxylase) located at chromosome 22q13.1. We performed a linkage study of chromosome 22 in 200 families with AS affected sibling-pairs. Association of alleles of the CYP2D6 gene was examined by both case-control and within-family means. For case-control studies, 617 unrelated individuals with AS (361 probands from sibling-pair and parent-case trio families and 256 unrelated non-familial sporadic cases) and 402 healthy ethnically matched controls were employed. For within-family association studies, 361 families including 161 parent-case trios and 200 affected sibling-pair families were employed. Homozygosity for poor metabolizer alleles was found to be associated with AS. Heterozygosity for the most frequent poor metabolizer allele (CYP2D6*4) was not associated with increased susceptibility to AS. Significant within-family association of CYP2D6*4 alleles and AS was demonstrated. Weak linkage was also demonstrated between CYP2D6 and AS. We postulate that altered metabolism of a natural toxin or antigen by the CYP2D6 gene may increase susceptibility to AS.

Secretor status and ankylosing spondylitis

Objective. To assess the genetic association between secretor status and ankylosing spondylitis (AS). Methods. A restriction digest method for determining secretor status was developed; 166 patients with AS and 216 healthy British white controls were typed for secretor status using this method. Results. The frequency of nonsecretors among patients with AS (47/166, 28%) was not significantly different from controls (72/216, 33%). Conclusion. Secretor status does not influence susceptibility to AS.

Risk analysis of prostate cancer in PRACTICAL, a multinational consortium, using 25 known prostate cancer susceptibility loci

Background Genome-wide association studies have identified multiple genetic variants associated with prostate cancer risk which explain a substantial proportion of familial relative risk. These variants can be used to stratify individuals by their risk of prostate cancer. Methods We genotyped 25 prostate cancer susceptibility loci in 40,414 individuals and derived a polygenic risk score (PRS).We estimated empirical odds ratios (OR) for prostate cancer associated with different risk strata defined by PRS and derived agespecific absolute risks of developing prostate cancer by PRS stratum and family history. Results The prostate cancer risk for men in the top 1% of the PRS distribution was 30.6 (95% CI, 16.4-57.3) fold compared with men in the bottom 1%, and 4.2 (95% CI, 3.2-5.5) fold compared with the median risk. The absolute risk of prostate cancer by age of 85 years was 65.8% for a man with family history in the top 1% of the PRS distribution, compared with 3.7% for a man in the bottom 1%. The PRS was only weakly correlated with serum PSA level (correlation = 0.09). Conclusions Risk profiling can identify men at substantially increased or reduced risk of prostate cancer. The effect size, measured by OR per unit PRS, was higher in men at younger ages and in men with family history of prostate cancer. Incorporating additional newly identified loci into a PRS should improve the predictive value of risk profiles. Impact:We demonstrate that the risk profiling based on SNPs can identify men at substantially increased or reduced risk that could have useful implications for targeted prevention and screening programs.