A mouse with an N-ethyl-N-nitrosourea (ENU) induced Trp589Arg Galnt3 mutation represents a model for hyperphosphataemic familial tumoural calcinosis

Mutations of UDP-N-acetyl-alpha-D-galactosamine polypeptide N-acetyl galactosaminyl transferase 3 (GALNT3) result in familial tumoural calcinosis (FTC) and the hyperostosis-hyperphosphataemia syndrome (HHS), which are autosomal recessive disorders characterised by soft-tissue calcification and hyperphosphataemia. To facilitate in vivo studies of these heritable disorders of phosphate homeostasis, we embarked on establishing a mouse model by assessing progeny of mice treated with the chemical mutagen N-ethyl-N-nitrosourea (ENU), and identified a mutant mouse, TCAL, with autosomal recessive inheritance of ectopic calcification, which involved multiple tissues, and hyperphosphataemia; the phenotype was designated TCAL and the locus, Tcal. TCAL males were infertile with loss of Sertoli cells and spermatozoa, and increased testicular apoptosis. Genetic mapping localized Tcal to chromosome 2 (62.64-71.11 Mb) which contained the Galnt3. DNA sequence analysis identified a Galnt3 missense mutation (Trp589Arg) in TCAL mice. Transient transfection of wild-type and mutant Galnt3-enhanced green fluorescent protein (EGFP) constructs in COS-7 cells revealed endoplasmic reticulum retention of the Trp589Arg mutant and Western blot analysis of kidney homogenates demonstrated defective glycosylation of Galnt3 in Tcal/Tcal mice. Tcal/Tcal mice had normal plasma calcium and parathyroid hormone concentrations; decreased alkaline phosphatase activity and intact Fgf23 concentrations; and elevation of circulating 1,25-dihydroxyvitamin D. Quantitative reverse transcriptase-PCR (qRT-PCR) revealed that Tcal/Tcal mice had increased expression of Galnt3 and Fgf23 in bone, but that renal expression of Klotho, 25-hydroxyvitamin D-1α-hydroxylase (Cyp27b1), and the sodium-phosphate co-transporters type-IIa and -IIc was similar to that in wild-type mice. Thus, TCAL mice have the phenotypic features of FTC and HHS, and provide a model for these disorders of phosphate metabolism. © 2012 Esapa et al.

Association of STAT3 and TNFRSF1A with ankylosing spondylitis in Han Chinese

Objectives: Recent association studies by the Australo-Anglo-American Spondyloarthritis Consortium (TASC) in Caucasian European populations from Australia, North America and the UK have identified a number of genes as being associated with ankylosing spondylitis (AS). A candidate gene study in a Han Chinese population was performed based on these findings to identify associated genes in this population. Methods: A case-control study was performed in a Han Chinese population of patients with AS (n=775) and controls (n=1587) from Shanghai and Nanjing. All patients met the modified New York criteria for AS. The cases and controls were genotyped for 115 single nucleotide polymorphisms (SNPs) tagging IL23R, ERAP1, STAT3, JAK2, TNFRSF1A and TRADD, as well as other confirmation SNPs from the TASC study, using the Sequenom iPlex and the ABI OpenArray platforms. Statistical analysis of genotyped SNPs was performed using the Cochran - Armitage test for trend and meta-analysis was performed using METAL. SNPs in AS-associated genes in this study were then imputed using MaCH, and association with AS tested by logistic regression. Results: SNPs in TNFRSF1A (rs4149577, p=8.2×10-4), STAT3 (rs2293152, p=0.0015; rs1053005, p=0.017) and ERAP1 (rs27038, p=0.0091; rs27037, p=0.0092) were significantly associated with AS in Han Chinese. Association was also observed between AS and the intergenic region 2p15 (rs10865331, p=0.023). The lack of association between AS and IL23R in Han Chinese was confirmed (all SNPs p>0.1). Conclusions: The study results demonstrate for the first time that genetic polymorphisms in STAT3, TNFRSF1A and 2p15 are associated with AS in Han Chinese, suggesting common pathogenic mechanisms for the disease in Chinese and Caucasian European populations. Furthermore, previous findings demonstrating that ERAP1, but not IL23R, is associated with AS in Chinese patients were confirmed.

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.

ERAP2 is associated with ankylosing spondylitis in HLA-B27-positive and HLA-B27-negative patients

The association of endoplasmic reticulum aminopeptidase 2 (ERAP2) with ankylosing spondylitis (AS) was recently described in the large International Genetics of AS Consortium Immunochip study...

Is there a higher genetic load of susceptibility loci in familial ankylosing spondylitis?

Objective Several genetic risk variants for ankylosing spondylitis (AS) have been identified in genome-wide association studies. Our objective was to examine whether familial AS cases have a higher genetic load of these susceptibility variants. Methods Overall, 502 AS patients were examined, consisting of 312 patients who had first-degree relatives (FDRs) with AS (familial) and 190 patients who had no FDRs with AS or spondylarthritis (sporadic). All patients and affected FDRs fulfilled the modified New York criteria for AS. The patients were recruited from 2 US cohorts (the North American Spondylitis Consortium and the Prospective Study of Outcomes in Ankylosing Spondylitis) and from the UK-Oxford cohort. The frequencies of AS susceptibility loci in IL-23R, IL1R2, ANTXR2, ERAP-1, 2 intergenic regions on chromosomes 2p15 and 21q22, and HLA-B27 status as determined by the tag single-nucleotide polymorphism (SNP) rs4349859 were compared between familial and sporadic cases of AS. Association between SNPs and multiplex status was assessed by logistic regression controlling for sibship size. Results HLA-B27 was significantly more prevalent in familial than sporadic cases of AS (odds ratio 4.44 [95% confidence interval 2.06, 9.55], P = 0.0001). Furthermore, the AS risk allele at chromosome 21q22 intergenic region showed a trend toward higher frequency in the multiplex cases (P = 0.08). The frequency of the other AS risk variants did not differ significantly between familial and sporadic cases, either individually or combined. Conclusion HLA-B27 is more prevalent in familial than sporadic cases of AS, demonstrating higher familial aggregation of AS in patients with HLA-B27 positivity. The frequency of the recently described non-major histocompatibility complex susceptibility loci is not markedly different between the sporadic and familial cases of AS.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5×10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6×10-8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98-1.14; p=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p=0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

Association of IL23R and ERAP1 genes with ankylosing spondylitis in a Portuguese population

Objective: Association between ankylosing spondylitis (AS) and two genes, ERAP1 and IL23R, has recently been reported in North American and British populations. The population attributable risk fraction for ERAP1 in this study was 25%, and for IL23R, 9%. Confirmation of these findings to ERAP1 in other ethnic groups has not yet been demonstrated. We sought to test the association between single nucleotide polymorphisms (SNPs) in these genes and susceptibility to AS among a Portuguese population. We also investigated the role of these genes in clinical manifestations of AS, including age of symptom onset, the Bath Ankylosing Spondylitis Disease Activity, Metrology and Functional Indices, and the modified Stoke Ankylosing Spondylitis Spinal Score. Methods: The study was conducted on 358 AS cases and 285 ethnically matched Portuguese healthy controls. AS was defined according to the modified New York Criteria. Genotyping of IL23R and ERAP1 allelic variants was carried out with TaqMan allelic discrimination assays. Association analysis was performed using the Cochrane-Armitage and linear regression tests of genotypes as implemented in PLINK for dichotomous and quantitative variables respectively. A meta-analysis for Portuguese and previously published Spanish IL23R data was performed using the StatsDirect® Statistical tools, by fixed and random effects models. Results: A total of 14 nsSNPs markers (8 for IL23R, 5 for ERAPl, 1 for LN-PEP) were analysed. Three markers (2 for IL23R and 1 for ERAP1) showed significant single-locus disease associations, confirming that the association of these genes with AS in the Portuguese population. The strongest associated SNP in IL23R was rs1004819 (OR=1.4, p=0.0049), and in ERAP1 was rs30187 (OR=1.26, p=0.035). The population attributable risk fractions in the Portuguese population for these SNPs are 11% and 9.7% respectively. No association was seen with any SNP in LN-PEP, which flanks ERAP1 and was associated with AS in the British population. No association was seen with clinical manifestations of AS. Conclusions: These results show that IL23R and ERAP1 genes are also associated with susceptibility to AS in the Portuguese population, and that they contribute a significant proportion of the population risk for this disease.

Atmospheric pressure gas plasma-induced colorectal cancer cell death is mediated by Nox2-ASK1 apoptosis pathways and oxidative stress is mitigated by Srx-Nrf2 anti-oxidant system

Atmospheric pressure gas plasma (AGP) generates reactive oxygen species (ROS) that induce apoptosis in cultured cancer cells. The majority of cancer cells develop a ROS-scavenging anti-oxidant system regulated by Nrf2, which confers resistance to ROS-mediated cancer cell death. Generation of ROS is involved in the AGP-induced cancer cell death of several colorectal cancer cells (Caco2, HCT116 and SW480) by activation of ASK1-mediated apoptosis signaling pathway without affecting control cells (human colonic sub-epithelial myofibroblasts; CO18, human fetal lung fibroblast; MRC5 and fetal human colon; FHC). However, the identity of an oxidase participating in AGP-induced cancer cell death is unknown. Here, we report that AGP up-regulates the expression of Nox2 (NADPH oxidase) to produce ROS. RNA interference designed to target Nox2 effectively inhibits the AGP-induced ROS production and cancer cell death. In some cases both colorectal cancer HT29 and control cells showed resistance to AGP treatment. Compared to AGP-sensitive Caco2 cells, HT29 cells show a higher basal level of the anti-oxidant system transcriptional regulator Nrf2 and its target protein sulfiredoxin (Srx) which are involved in cellular redox homeostasis. Silencing of both Nrf2 and Srx sensitized HT29 cells, leads to ROS overproduction and decreased cell viability. This indicates that in HT29 cells, Nrf2/Srx axis is a protective factor against AGP-induced oxidative stress. The inhibition of Nrf2/Srx signaling should be considered as a central target in drug-resistant colorectal cancer treatments.

Cardiometabolic effects of genetic upregulation of the interleukin 1 receptor antagonist: A Mendelian randomisation analysis

Background To investigate potential cardiovascular and other effects of long-term pharmacological interleukin 1 (IL-1) inhibition, we studied genetic variants that produce inhibition of IL-1, a master regulator of inflammation. Methods We created a genetic score combining the effects of alleles of two common variants (rs6743376 and rs1542176) that are located upstream of IL1RN, the gene encoding the IL-1 receptor antagonist (IL-1Ra; an endogenous inhibitor of both IL-1α and IL-1β); both alleles increase soluble IL-1Ra protein concentration. We compared effects on inflammation biomarkers of this genetic score with those of anakinra, the recombinant form of IL-1Ra, which has previously been studied in randomised trials of rheumatoid arthritis and other inflammatory disorders. In primary analyses, we investigated the score in relation to rheumatoid arthritis and four cardiometabolic diseases (type 2 diabetes, coronary heart disease, ischaemic stroke, and abdominal aortic aneurysm; 453 411 total participants). In exploratory analyses, we studied the relation of the score to many disease traits and to 24 other disorders of proposed relevance to IL-1 signalling (746 171 total participants). Findings For each IL1RN minor allele inherited, serum concentrations of IL-1Ra increased by 0·22 SD (95% CI 0·18–0·25; 12·5%; p=9·3 × 10−33), concentrations of interleukin 6 decreased by 0·02 SD (−0·04 to −0·01; −1·7%; p=3·5 × 10−3), and concentrations of C-reactive protein decreased by 0·03 SD (−0·04 to −0·02; −3·4%; p=7·7 × 10−14). We noted the effects of the genetic score on these inflammation biomarkers to be directionally concordant with those of anakinra. The allele count of the genetic score had roughly log-linear, dose-dependent associations with both IL-1Ra concentration and risk of coronary heart disease. For people who carried four IL-1Ra-raising alleles, the odds ratio for coronary heart disease was 1·15 (1·08–1·22; p=1·8 × 10−6) compared with people who carried no IL-1Ra-raising alleles; the per-allele odds ratio for coronary heart disease was 1·03 (1·02–1·04; p=3·9 × 10−10). Per-allele odds ratios were 0·97 (0·95–0·99; p=9·9 × 10−4) for rheumatoid arthritis, 0·99 (0·97–1·01; p=0·47) for type 2 diabetes, 1·00 (0·98–1·02; p=0·92) for ischaemic stroke, and 1·08 (1·04–1·12; p=1·8 × 10−5) for abdominal aortic aneurysm. In exploratory analyses, we observed per-allele increases in concentrations of proatherogenic lipids, including LDL-cholesterol, but no clear evidence of association for blood pressure, glycaemic traits, or any of the 24 other disorders studied. Modelling suggested that the observed increase in LDL-cholesterol could account for about a third of the association observed between the genetic score and increased coronary risk. Interpretation Human genetic data suggest that long-term dual IL-1α/β inhibition could increase cardiovascular risk and, conversely, reduce the risk of development of rheumatoid arthritis. The cardiovascular risk might, in part, be mediated through an increase in proatherogenic lipid concentrations. Funding UK Medical Research Council, British Heart Foundation, UK National Institute for Health Research, National Institute for Health Research Cambridge Biomedical Research Centre, European Research Council, and European Commission Framework Programme 7.

An N-ethyl-n-nitrosourea induced corticotropin-releasing hormone promoter mutation provides a mouse model for endogenous glucocorticoid excess

Cushing's syndrome, which is characterized by excessive circulating glucocorticoid concentrations, maybe due to ACTH-dependent or -independent causes that include anterior pituitary and adrenal cortical tumors, respectively. ACTH secretion is stimulated by CRH, and we report a mouse model for Cushing's syndrome due to an N-ethyl-N-nitrosourea (ENU) induced Crh mutation at -120 bp of the promoter region, which significantly increased luciferase reporter activity and was thus a gain-of-function mutation. Crh -120/+ mice, when compared with wild-type littermates, had obesity, muscle wasting, thin skin, hair loss, and elevated plasma and urinary concentrations of corticosterone. In addition, Crh-120/+ mice had hyperglycemia, hyperfructosaminemia, hyperinsulinemia, hypercholesterolemia, hypertriglyceridemia, and hyperleptinemia but normal adiponectin. Crh -120/+ mice also had low bone mineral density, hypercalcemia, hypercalciuria, and decreased concentrations of plasma PTH and osteocalcin. Bone histomorphometry revealed Crh-120/+ mice to have significant reductions in mineralizing surface area, mineral apposition, bone formation rates, osteoblast number, and the percentage of corticoendosteal bone covered by osteoblasts, which was accompanied by an increase in adipocytes in the bone marrow. Thus, a mouse model for Cushing's syndrome has been established, and this will help in further elucidating the pathophysiological effects of glucocorticoid excess and in evaluating treatments for corticosteroid-induced osteoporosis.

Major histocompatibility complex associations of ankylosing spondylitis are complex and involve further epistasis with ERAP1

Ankylosing spondylitis (AS) is a common, highly heritable, inflammatory arthritis for which HLA-B*27 is the major genetic risk factor, although its role in the aetiology of AS remains elusive. To better understand the genetic basis of the MHC susceptibility loci, we genotyped 7,264 MHC SNPs in 22,647 AS cases and controls of European descent. We impute SNPs, classical HLA alleles and amino-acid residues within HLA proteins, and tested these for association to AS status. Here we show that in addition to effects due to HLA-B*27 alleles, several other HLA-B alleles also affect susceptibility. After controlling for the associated haplotypes in HLA-B, we observe independent associations with variants in the HLA-A, HLA-DPB1 and HLA-DRB1 loci. We also demonstrate that the ERAP1 SNP rs30187 association is not restricted only to carriers of HLA-B*27 but also found in HLA-B*40:01 carriers independently of HLA-B*27 genotype.

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.

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.

Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci

Ankylosing spondylitis is a common, highly heritable inflammatory arthritis affecting primarily the spine and pelvis. In addition to HLA-B*27 alleles, 12 loci have previously been identified that are associated with ankylosing spondylitis in populations of European ancestry, and 2 associated loci have been identified in Asians. In this study, we used the Illumina Immunochip microarray to perform a case-control association study involving 10,619 individuals with ankylosing spondylitis (cases) and 15,145 controls. We identified 13 new risk loci and 12 additional ankylosing spondylitis-associated haplotypes at 11 loci. Two ankylosing spondylitis-associated regions have now been identified encoding four aminopeptidases that are involved in peptide processing before major histocompatibility complex (MHC) class I presentation. Protective variants at two of these loci are associated both with reduced aminopeptidase function and with MHC class I cell surface expression.

Using genome-Wwide complex trait analysis to quantify 'missing heritability' in Parkinson's disease

Genome-wide association studies (GWASs) have been successful at identifying single-nucleotide polymorphisms (SNPs) highly associated with common traits; however, a great deal of the heritable variation associated with common traits remains unaccounted for within the genome. Genome-wide complex trait analysis (GCTA) is a statistical method that applies a linear mixed model to estimate phenotypic variance of complex traits explained by genome-wide SNPs, including those not associated with the trait in a GWAS. We applied GCTA to 8 cohorts containing 7096 case and 19 455 control individuals of European ancestry in order to examine the missing heritability present in Parkinson's disease (PD). We meta-analyzed our initial results to produce robust heritability estimates for PD types across cohorts. Our results identify 27% (95% CI 17-38, P = 8.08E - 08) phenotypic variance associated with all types of PD, 15% (95% CI -0.2 to 33, P = 0.09) phenotypic variance associated with early-onset PD and 31% (95% CI 17-44, P = 1.34E - 05) phenotypic variance associated with late-onset PD. This is a substantial increase from the genetic variance identified by top GWAS hits alone (between 3 and 5%) and indicates there are substantially more risk loci to be identified. Our results suggest that although GWASs are a useful tool in identifying the most common variants associated with complex disease, a great deal of common variants of small effect remain to be discovered. © Published by Oxford University Press 2012.