Gene expression profiling reveals a downregulation in immune-associated genes in patients with AS

Objective: To identify differentially expressed genes in peripheral blood mononuclear cells (PBMCs) from patients with ankylosing spondylitis (AS) compared with healthy individuals. Methods: RNA was extracted from PBMCs collected from 18 patients with active disease and 18 gender-matched and age-matched controls. Expression profiles of these cells were determined using microarray. Candidate genes with differential expressions were confirmed in the same samples using quantitative reverse transcription-PCR (qRT-PCR). These genes were then validated in a different sample cohort of 35 patients with AS and 18 controls by qRT-PCR. Results: Microarray analysis identified 452 genes detected with 485 probes which were differentially expressed between patients with AS and controls. Underexpression of NR4A2, tumour necrosis factor AIP3 (TNFAIP3) and CD69 was confirmed. These genes were further validated in a different sample group in which the patients with AS had a wider range of disease activity. Predictive algorithms were also developed from the expression data using receiver-operating characteristic curves, which demonstrated that the three candidate genes have ∼80% power to predict AS according to their expression levels. Conclusions: The findings show differences in global gene expression patterns between patients with AS and controls, suggesting an immunosuppressive phenotype in the patients. Furthermore, downregulated expression of three immune-related genes was confirmed. These candidate genes were also shown to be strong predictive markers for AS.

Ankylosing spondylitis is associated with the anthrax toxin receptor 2 gene (ANTXR2)

Objectives ANTXR2 variants have been associated with ankylosing spondylitis (AS) in two previous genome-wide association studies (GWAS) (p∼9×10-8). However, a genome-wide significant association (p<5×10-8) was not observed. We conducted a more comprehensive analysis of ANTXR2 in an independent UK sample to confirm and refine this association. Methods A replication study was carried out with 2978 cases and 8365 controls. Then, these were combined with non-overlapping samples from the two previous GWAS in a meta-analysis. Human leukocyte antigen (HLA)-B27 stratification was also performed to test for ANTXR2-HLA-B27 interaction. Results Out of nine single nucleotide polymorphisms (SNP) in the study, five SNPs were nominally associated (p<0.05) with AS in the replication dataset. In the meta-analysis, eight SNPs showed evidence of association, the strongest being with rs12504282 (OR=0.88, p=6.7×10-9). Seven of these SNPs showed evidence for association in the HLA-B27-positive subgroup, but none was associated with HLA-B27-negative AS. However, no statistically significant interaction was detected between HLA-B27 and ANTXR2 variants. Conclusions ANTXR2 variants are clearly associated with AS. The top SNPs from two previous GWAS (rs4333130 and rs4389526) and this study (rs12504282) are in strong linkage disequilibrium (r2≥0.76). All are located near a putative regulatory region. Further studies are required to clarify the role played by these ANTXR2 variants in AS.

The chromosome 16q region associated with ankylosing spondylitis includes the candidate gene tumour necrosis factor receptor type 1-associated death domain (TRADD)

Objective: To replicate and refine the reported association of ankylosing spondylitis (AS) with two nonsynonymous single nucleotide polymorphisms (nsSNPs) on chromosome 16q22.1. Methods: Firstly, 730 independent UK patients with AS were genotyped for rs9939768 and rs6979 and allele frequencies were compared with 2879 previously typed historic disease controls. Secondly, the two data sets were combined in meta-analyses. Finally, 5 tagging SNPs, located between rs9939768 and rs6979, were analysed in 1604 cases and 1020 controls. Results: The association of rs6979 with AS was replicated, p=0.03, OR=1.14 (95% CI 1.01 to 1.28), and a trend for association with rs9939768 detected, p=0.06, OR=1.25 (95% CI 0.99 to 1.57). Meta-analyses revealed association of both SNPs with AS, p=0.0008, OR=1.31 (95% CI 1.12 to 1.54) and p=0.0009, OR=1.15 (95% CI 1.06 to 1.23) for rs9939768 and rs6979, respectively. New associations with rs9033 and rs868213 (p=0.00002, OR=1.23 (95% CI 1.12 to 1.36) and p=0.00002 OR=1.45 (95% CI 1.22 to 1.72), respectively, were identified. Conclusions: The region on chromosome 16 that has been replicated in the present work is interesting as the highly plausible candidate gene, tumour necrosis factor receptor type 1 (TNFR1)-associated death domain (TRADD), is located between rs9033 and rs868213. It will require additional work to identify the primary genetic association(s) with AS.

Prospective meta-analysis of interleukin 1 gene complex polymorphisms confirms associations with ankylosing spondylitis

Objectives: The aim of the current study was to determine the contribution of interleukin (IL) 1 gene cluster polymorphisms previously implicated in susceptibility for ankylosing spondylitis (AS) to AS susceptibility in different populations worldwide. Methods: Nine polymorphisms in the IL1 gene cluster members IL1A (rs2856836, rs17561 and rs1894399), IL1B (rs16944), IL1F10 (rs3811058) and IL1RN (rs419598, the IL1RA VNTR, rs315952 and rs315951) were genotyped in 2675 AS cases and 2592 healthy controls recruited in 12 different centres in 10 countries. Association of variants with AS was tested by Mantel-Haenszel random effects analysis. Results: Strong association was observed with three single nucleotide polymorphisms (SNPs) in the IL1A gene (rs2856836, rs17561, rs1894399, p = 0.0036, 0.000019 and 0.0003, respectively). There was no evidence of significant heterogeneity of effects between centres, and no evidence of non-combinability of findings. The population attributable risk fraction of these variants in Caucasians is estimated at 4-6%. Conclusions: This study confirms that IL1A is associated with susceptibility to AS. Association of the other IL1 gene complex members could not be excluded in specific populations. Prospective meta-analysis is a useful tool in confirmation studies of genes associated with complex genetic disorders such as AS, providing sufficiently large sample sizes to produce robust findings often not achieved in smaller individual cohorts.

VH gene usage in immunoglobulin E responses of seasonal rhinitis patients allergic to grass pollen is oligoclonal and antigen driven

Background: IgE is the pivotal-specific effector molecule of allergic reactions yet it remains unclear whether the elevated production of IgE in atopic individuals is due to superantigen activation of B cell populations, increased antibody class switching to IgE or oligoclonal allergen-driven IgE responses. Objectives: To increase our understanding of the mechanisms driving IgE responses in allergic disease we examined immunoglobulin variable regions of IgE heavy chain transcripts from three patients with seasonal rhinitis due to grass pollen allergy. Methods: Variable domain of heavy chain-epsilon constant domain 1 cDNAs were amplified from peripheral blood using a two-step semi-nested PCR, cloned and sequenced. Results: The VH gene family usage in subject A was broadly based, but there were two clusters of sequences using genes VH 3-9 and 3-11 with unusually low levels of somatic mutations, 0-3%. Subject B repeatedly used VH 1-69 and subject C repeatedly used VH 1-02, 1-46 and 5a genes. Most clones were highly mutated being only 86-95% homologous to their germline VH gene counterparts and somatic mutations were more abundant at the complementarity determining rather than framework regions. Multiple sequence alignment revealed both repeated use of particular VH genes as well as clonal relatedness among clusters of IgE transcripts. Conclusion: In contrast to previous studies we observed no preferred VH gene common to IgE transcripts of the three subjects allergic to grass pollen. Moreover, most of the VH gene characteristics of the IgE transcripts were consistent with oligoclonal antigen-driven IgE responses.

A Chlamydia trachomatis strain with a chemically generated amino acid substitution (P370L) in the cthtrA gene shows reduced elementary body production

Background Chlamydia (C.) trachomatis is the most prevalent bacterial sexually transmitted infection worldwide and the leading cause of preventable blindness. Genetic approaches to investigate C. trachomatis have been only recently developed due to the organism’s intracellular developmental cycle. HtrA is a critical stress response serine protease and chaperone for many bacteria and in C. trachomatis has been previously shown to be important for heat stress and the replicative phase of development using a chemical inhibitor of the CtHtrA activity. In this study, chemically-induced SNVs in the cthtrA gene that resulted in amino acid substitutions (A240V, G475E, and P370L) were identified and characterized. Methods SNVs were initially biochemically characterized in vitro using recombinant protein techniques to confirm a functional impact on proteolysis. The C. trachomatis strains containing the SNVs with marked reductions in proteolysis were investigated in cell culture to identify phenotypes that could be linked to CtHtrA function. Results The strain harboring the SNV with the most marked impact on proteolysis (cthtrAP370L) was detected to have a significant reduction in the production of infectious elementary bodies. Conclusions This provides genetic evidence that CtHtrA is critical for the C. trachomatis developmental cycle.

Genetic association and gene expression studies suggest that genetic variants in the SYNE1 and TNF genes are related to menstrual migraine

BACKGROUND: Menstrual migraine (MM) encompasses pure menstrual migraine (PMM) and menstrually-related migraine (MRM). This study was aimed at investigating genetic variants that are potentially related to MM, specifically undertaking genotyping and mRNA expression analysis of the ESR1, PGR, SYNE1 and TNF genes in MM cases and non-migraine controls. METHODS: A total of 37 variants distributed across 14 genes were genotyped in 437 DNA samples (282 cases and 155 controls). In addition levels of gene expression were determined in 74 cDNA samples (41 cases and 33 controls). Association and correlation analysis were performed using Plink and RStudio. RESULTS: SNPs rs3093664 and rs9371601 in TNF and SYNE1 genes respectively, were significantly associated with migraine in the MM population (p = 0.008; p = 0.009 respectively). Analysis of qPCR results found no significant difference in levels of gene expression between cases and controls. However, we found a significant correlation between the expression of ESR1 and SYNE1, ESR1 and PGR and TNF and SYNE1 in samples taken during the follicular phase of the menstrual cycle. CONCLUSIONS: Our results show that SNPs rs9371601 and rs3093664 in the SYNE1 and TNF genes respectively, are associated with MM. The present study also provides strong evidence to support the correlation of ESR1, PGR, SYNE1 and TNF gene expression in MM.

Mobile gene silencing in Arabidopsis is regulated by hydrogen peroxide

In plants and nematodes, RNAi can spread from cells from which it is initiated to other cells in the organism. The underlying mechanism controlling the mobility of RNAi signals is not known, especially in the case of plants. A genetic screen designed to recover plants impaired in the movement but not the production or effectiveness of the RNAi signal identified RCI3, which encodes a hydrogen peroxide (H2O2)-producing type III peroxidase, as a key regulator of silencing mobility in Arabidopsis thaliana. Silencing initiated in the roots of rci3 plants failed to spread into leaf tissue or floral tissue. Application of exogenous H2O2 reinstated the spread in rci3 plants and accelerated it in wild-type plants. The addition of catalase or MnO2, which breaks down H2O2, slowed the spread of silencing in wild-type plants. We propose that endogenous H2O2, under the control of peroxidases, regulates the spread of gene silencing by altering plasmodesmata permeability through remodelling of local cell wall structure, and may play a role in regulating systemic viral defence.

Impulsivity-related cognition in alcohol dependence: Is it moderated by DRD2/ANKK1 gene status and executive dysfunction?

Perceived impaired control over alcohol use is a key cognitive construct in alcohol dependence that has been related prospectively to treatment outcome and may mediate the risk for problem drinking conveyed by impulsivity in non-dependent drinkers. The aim of the current study was to investigate whether perceived impaired control may mediate the association between impulsivity-related measures (derived from the Short-form Eysenck Personality Questionnaire-Revised) and alcohol-dependence severity in alcohol-dependent drinkers. Furthermore, the extent to which this hypothesized relationship was moderated by genetic risk (Taq1A polymorphism in the DRD2/ANKK1 gene cluster) and verbal fluency as an indicator of executive cognitive ability (Controlled Oral Word Association Test) was also examined. A sample of 143 alcohol-dependent inpatients provided an extensive clinical history of their alcohol use, gave 10ml of blood for DNA analysis, and completed self-report measures relating to impulsivity, impaired control and severity of dependence. As hypothesized, perceived impaired control (partially) mediated the association between impulsivity-related measures and alcohol-dependence severity. This relationship was not moderated by the DRD2/ANKK1 polymorphism or verbal fluency. These results suggest that, in alcohol dependence, perceived impaired control is a cognitive mediator of impulsivity-related constructs that may be unaffected by DRD2/ANKK1 and neurocognitive processes underlying the retrieval of verbal information

Association of a disintegrin and metalloprotease 33 (ADAM33) gene polymorphisms with the risk of COPD: An updated meta-analysis of 2,644 cases and 4,804 controls

A series of observational studies have been made to investigate the association of the ADAM33 gene polymorphisms with the risk of COPD, but their results were conflicting. Therefore, we performed an updated meta-analysis to quantitatively summarize the associations of ADAM33 gene polymorphisms with the risk of COPD. Thirteen case–control studies referring to nine SNPs were identified: V4 (rs2787094), T+1 (rs2280089), T2 (rs2280090), T1 (rs2280091), S2 (rs528557), S1 (rs3918396), Q−1 (rs612709), F+1 (rs511898) and ST+5 (rs597980). A dominant model (AA+Aa vs. aa), recessive model (AA vs. Aa+aa), additive model (AA vs. aa) and allelic model (A vs. a) were used to evaluate the association of ADAM33 polymorphism with the risk of COPD. The results indicated that significant associations were found for ADAM33 T1, T2, S1, Q−1, F+1 and ST+5 polymorphisms associated with the risk of COPD in different populations. However, no significant associations were found for V4, T+1 and S2 polymorphisms with the risk of COPD in all genetic models, even in the subgroup analysis by ethnicity. This meta-analysis provided evidence that the ADAM33 T1, T2, S1, Q−1, F+1 and ST+5 six locus polymorphisms association with the risk of COPD. Furthermore, T2, Q−1 and ST+5 indicated an association with the risk of COPD in the European populations, whereas T1, T2, S1, F+1 and Q−1 indicated an association with the risk of COPD in the Asian populations.

Genome-wide association study for sight-threatening diabetic retinopathy reveals association with genetic variation near the GRB2 gene

Aims/hypothesis Diabetic retinopathy is a serious complication of diabetes mellitus and can lead to blindness. A genetic component, in addition to traditional risk factors, has been well described although strong genetic factors have not yet been identified. Here, we aimed to identify novel genetic risk factors for sight-threatening diabetic retinopathy using a genome-wide association study. Methods Retinopathy was assessed in white Australians with type 2 diabetes mellitus. Genome-wide association analysis was conducted for comparison of cases of sight-threatening diabetic retinopathy (n = 336) with diabetic controls with no retinopathy (n = 508). Top ranking single nucleotide polymorphisms were typed in a type 2 diabetes replication cohort, a type 1 diabetes cohort and an Indian type 2 cohort. A mouse model of proliferative retinopathy was used to assess differential expression of the nearby candidate gene GRB2 by immunohistochemistry and quantitative western blot. Results The top ranked variant was rs3805931 with p = 2.66 × 10−7, but no association was found in the replication cohort. Only rs9896052 (p = 6.55 × 10−5) was associated with sight-threatening diabetic retinopathy in both the type 2 (p = 0.035) and the type 1 (p = 0.041) replication cohorts, as well as in the Indian cohort (p = 0.016). The study-wide meta-analysis reached genome-wide significance (p = 4.15 × 10−8). The GRB2 gene is located downstream of this variant and a mouse model of retinopathy showed increased GRB2 expression in the retina. Conclusions/interpretation Genetic variation near GRB2 on chromosome 17q25.1 is associated with sight-threatening diabetic retinopathy. Several genes in this region are promising candidates and in particular GRB2 is upregulated during retinal stress and neovascularisation.

Founder p.Arg 446* mutation in the PDHX gene explains over half of cases with congenital lactic acidosis in Roma children

Investigation of 31 of Roma patients with congenital lactic acidosis (CLA) from Bulgaria identified homozygosity for the R446* mutation in the PDHX gene as the most common cause of the disorder in this ethnic group. It accounted for around 60% of patients in the study and over 25% of all CLA cases referred to the National Genetic Laboratory in Bulgaria. The detection of a homozygous patient from Hungary and carriers among population controls from Romania and Slovakia suggests a wide spread of the mutation in the European Roma population. The clinical phenotype of the twenty R446* homozygotes was relatively homogeneous, with lactic acidosis crisis in the first days or months of life as the most common initial presentation (15/20 patients) and delayed psychomotor development and/or seizures in infancy as the leading manifestations in a smaller group (5/20 patients). The subsequent clinical picture was dominated by impaired physical growth and a very consistent pattern of static cerebral palsy-like encephalopathy with spasticity and severe to profound mental retardation seen in over 80% of cases. Most patients had a positive family history. We propose testing for the R446* mutation in PDHX as a rapid first screening in Roma infants with metabolic acidosis. It will facilitate and accelerate diagnosis in a large proportion of cases, allow early rehabilitation to alleviate the chronic clinical course, and prevent further affected births in high-risk families.

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.

Elucidating the chromosome 9 association with AS; CARD9 is a candidate gene

Ankylosing spondylitis (AS) is polygenic with contributions from the immunologically relevant genes HLA-B27, ERAP1 and IL23R. A recent genome-wide association screen (GWAS) identified associations (P0.005) with the non-synonymous single-nucleotide polymorphisms (nsSNPs), rs4077515 and rs3812571, in caspase recruitment domain-containing protein 9 (CARD9) and small nuclear RNA-activating complex polypeptide 4 (SNAPC4) on chromosome 9q that had previously been linked to AS. We replicated these associations in a study of 730 AS patients compared with 2879 historic disease controls (rs4077515 P0.0004, odds ratio (OR)1.2, 95% confidence interval (CI)1.1-1.4; rs3812571 P0.0003, OR1.2, 95% CI1.1-1.4). Meta-analysis revealed strong associations of both SNPs with AS, rs4077515 P0.000005, OR1.2, 95% CI1.1-1.3 and rs3812571 P0.000006, OR1.2, 95% CI1.1-1.3. We then typed 1604 AS cases and 1020 controls for 13 tagging SNPs; 6 showed at least nominal association, 5 of which were in CARD9. We imputed genotypes for 13 additional SNPs but none was more strongly associated with AS than the tagging SNPs. Finally, interrogation of an mRNA expression database revealed that the SNPs most strongly associated with AS (or in strong linkage disequilibrium) were those most associated with CARD9 expression. CARD9 is a plausible candidate for AS given its central role in the innate immune response.