Matters Arising: Zhu et al, "A novel gene variation of TNFα associated with ankylosing spondylitis: A reconfirmed study"

In a recently published study1 involving 79 ankylosing spondylitis (AS) patients and 132 unrelated healthy blood donors, Zhu et al report association of a tumour necrosis factor α (TNFα) single nucleotide polymorphism (SNP), rs1799724, with AS.1 I have many concerns with this paper and its conclusions...

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.

Whole blood transcriptional profiling in ankylosing spondylitis identifies novel candidate genes that might contribute to the inflammatory and tissue-destructive disease aspects

Introduction: A number of genetic-association studies have identified genes contributing to ankylosing spondylitis (AS) susceptibility but such approaches provide little information as to the gene activity changes occurring during the disease process. Transcriptional profiling generates a 'snapshot' of the sampled cells' activity and thus can provide insights into the molecular processes driving the disease process. We undertook a whole-genome microarray approach to identify candidate genes associated with AS and validated these gene-expression changes in a larger sample cohort. Methods: A total of 18 active AS patients, classified according to the New York criteria, and 18 gender- and age-matched controls were profiled using Illumina HT-12 whole-genome expression BeadChips which carry cDNAs for 48,000 genes and transcripts. Class comparison analysis identified a number of differentially expressed candidate genes. These candidate genes were then validated in a larger cohort using qPCR-based TaqMan low density arrays (TLDAs). Results: A total of 239 probes corresponding to 221 genes were identified as being significantly different between patients and controls with a P-value <0.0005 (80% confidence level of false discovery rate). Forty-seven genes were then selected for validation studies, using the TLDAs. Thirteen of these genes were validated in the second patient cohort with 12 downregulated 1.3- to 2-fold and only 1 upregulated (1.6-fold). Among a number of identified genes with well-documented inflammatory roles we also validated genes that might be of great interest to the understanding of AS progression such as SPOCK2 (osteonectin) and EP300, which modulate cartilage and bone metabolism. Conclusions: We have validated a gene expression signature for AS from whole blood and identified strong candidate genes that may play roles in both the inflammatory and joint destruction aspects of the disease.

Fall rates in hospital rehabilitation units after individualised patient and staff education programmes: A pragmatic, stepped-wedge, cluster-randomised controlled trial

- Background Falls are the most frequent adverse events that are reported in hospitals. We examined the effectiveness of individualised falls-prevention education for patients, supported by training and feedback for staff, delivered as a ward-level programme. - Methods Eight rehabilitation units in general hospitals in Australia participated in this stepped-wedge, cluster-randomised study, undertaken during a 50 week period. Units were randomly assigned to intervention or control groups by use of computer-generated, random allocation sequences. We included patients admitted to the unit during the study with a Mini-Mental State Examination (MMSE) score of more than 23/30 to receive individualised education that was based on principles of changes in health behaviour from a trained health professional, in addition to usual care. We provided information about patients' goals, feedback about the ward environment, and perceived barriers to engagement in falls-prevention strategies to staff who were trained to support the uptake of strategies by patients. The coprimary outcome measures were patient rate of falls per 1000 patient-days and the proportion of patients who were fallers. All analyses were by intention to treat. This trial is registered with the Australian New Zealand Clinical Trials registry, number ACTRN12612000877886). - Findings Between Jan 13, and Dec 27, 2013, 3606 patients were admitted to the eight units (n=1983 control period; n=1623 intervention period). There were fewer falls (n=196, 7·80/1000 patient-days vs n=380, 13·78/1000 patient-days, adjusted rate ratio 0·60 [robust 95% CI 0·42–0·94], p=0·003), injurious falls (n=66, 2·63/1000 patient-days vs 131, 4·75/1000 patient-days, 0·65 [robust 95% CI 0·42–0·88], p=0·006), and fallers (n=136 [8·38%] vs n=248 [12·51%] adjusted odds ratio 0·55 [robust 95% CI 0·38 to 0·81], p=0·003) in the intervention compared with the control group. There was no significant difference in length of stay (intervention median 11 days [IQR 7–19], control 10 days [6–18]). - Interpretation Individualised patient education programmes combined with training and feedback to staff added to usual care reduces the rates of falls and injurious falls in older patients in rehabilitation hospital-units.

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.

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.

A novel near-infrared spectroscopy and chemometrics method for rapid analysis of several chemical components and antioxidant activity of mint (mentha haplocalyx briq.) samples

A novel near-infrared spectroscopy (NIRS) method has been researched and developed for the simultaneous analyses of the chemical components and associated properties of mint (Mentha haplocalyx Briq.) tea samples. The common analytes were: total polysaccharide content, total flavonoid content, total phenolic content, and total antioxidant activity. To resolve the NIRS data matrix for such analyses, least squares support vector machines was found to be the best chemometrics method for prediction, although it was closely followed by the radial basis function/partial least squares model. Interestingly, the commonly used partial least squares was unsatisfactory in this case. Additionally, principal component analysis and hierarchical cluster analysis were able to distinguish the mint samples according to their four geographical provinces of origin, and this was further facilitated with the use of the chemometrics classification methods-K-nearest neighbors, linear discriminant analysis, and partial least squares discriminant analysis. In general, given the potential savings with sampling and analysis time as well as with the costs of special analytical reagents required for the standard individual methods, NIRS offered a very attractive alternative for the simultaneous analysis of mint samples.

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.

Runs of homozygosity and a cluster of vulvar cancer in young Australian Aboriginal women

Objective A cluster of vulvar cancer exists in young Aboriginal women living in remote communities in Arnhem Land, Australia. A genetic case–control study was undertaken involving 30 cases of invasive vulvar cancer and its precursor lesion, high-grade vulvar intraepithelial neoplasia (VIN), and 61 controls, matched for age and community of residence. It was hypothesized that this small, isolated population may exhibit increased autozygosity, implicating recessive effects as a possible mechanism for increased susceptibility to vulvar cancer. Methods Genotyping data from saliva samples were used to identify runs of homozygosity (ROH) in order to calculate estimates of genome-wide homozygosity. Results No evidence of an effect of genome-wide homozygosity on vulvar cancer and VIN in East Arnhem women was found, nor was any individual ROH found to be significantly associated with case status. This study found further evidence supporting an association between previous diagnosis of CIN and diagnosis of vulvar cancer or VIN, but found no association with any other medical history variable. Conclusions These findings do not eliminate the possibility of genetic risk factors being involved in this cancer cluster, but rather suggest that alternative analytical strategies and genetic models should be explored.

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.

Mosaic PPM1D mutations are associated with predisposition to breast and ovarian cancer

Improved sequencing technologies offer unprecedented opportunities for investigating the role of rare genetic variation in common disease. However, there are considerable challenges with respect to study design, data analysis and replication. Using pooled next-generation sequencing of 507 genes implicated in the repair of DNA in 1,150 samples, an analytical strategy focused on protein-truncating variants (PTVs) and a large-scale sequencing case-control replication experiment in 13,642 individuals, here we show that rare PTVs in the p53-inducible protein phosphatase PPM1D are associated with predisposition to breast cancer and ovarian cancer. PPM1D PTV mutations were present in 25 out of 7,781 cases versus 1 out of 5,861 controls (P = 1.12 × 10-5), including 18 mutations in 6,912 individuals with breast cancer (P = 2.42 × 10-4) and 12 mutations in 1,121 individuals with ovarian cancer (P = 3.10 × 10-9). Notably, all of the identified PPM1D PTVs were mosaic in lymphocyte DNA and clustered within a 370-base-pair region in the final exon of the gene, carboxy-terminal to the phosphatase catalytic domain. Functional studies demonstrate that the mutations result in enhanced suppression of p53 in response to ionizing radiation exposure, suggesting that the mutant alleles encode hyperactive PPM1D isoforms. Thus, although the mutations cause premature protein truncation, they do not result in the simple loss-of-function effect typically associated with this class of variant, but instead probably have a gain-of-function effect. Our results have implications for the detection and management of breast and ovarian cancer risk. More generally, these data provide new insights into the role of rare and of mosaic genetic variants in common conditions, and the use of sequencing in their identification.

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.