Identification of IDUA and WNT16 phosphorylation-related non-synonymous polymorphisms for bone mineral density in meta-analyses of genome-wide association studies

Protein phosphorylation regulates a wide variety of cellular processes. Thus, we hypothesize that single-nucleotide polymorphisms (SNPs) that may modulate protein phosphorylation could affect osteoporosis risk. Based on a previous conventional genome-wide association (GWA) study, we conducted a three-stage meta-analysis targeting phosphorylation-related SNPs (phosSNPs) for femoral neck (FN)-bone mineral density (BMD), total hip (HIP)-BMD, and lumbar spine (LS)-BMD phenotypes. In stage 1, 9593 phosSNPs were meta-analyzed in 11,140 individuals of various ancestries. Genome-wide significance (GWS) and suggestive significance were defined by α = 5.21 × 10–6 (0.05/9593) and 1.00 × 10–4, respectively. In stage 2, nine stage 1–discovered phosSNPs (based on α = 1.00 × 10–4) were in silico meta-analyzed in Dutch, Korean, and Australian cohorts. In stage 3, four phosSNPs that replicated in stage 2 (based on α = 5.56 × 10–3, 0.05/9) were de novo genotyped in two independent cohorts. IDUA rs3755955 and rs6831280, and WNT16 rs2707466 were associated with BMD phenotypes in each respective stage, and in three stages combined, achieving GWS for both FN-BMD (p = 8.36 × 10–10, p = 5.26 × 10–10, and p = 3.01 × 10–10, respectively) and HIP-BMD (p = 3.26 × 10–6, p = 1.97 × 10–6, and p = 1.63 × 10–12, respectively). Although in vitro studies demonstrated no differences in expressions of wild-type and mutant forms of IDUA and WNT16B proteins, in silico analyses predicts that WNT16 rs2707466 directly abolishes a phosphorylation site, which could cause a deleterious effect on WNT16 protein, and that IDUA phosSNPs rs3755955 and rs6831280 could exert indirect effects on nearby phosphorylation sites. Further studies will be required to determine the detailed and specific molecular effects of these BMD-associated non-synonymous variants. © 2015 American Society for Bone and Mineral Research.

Disease-associated polymorphisms in ERAP1 do not alter endoplasmic reticulum stress in patients with ankylosing spondylitis

The mechanism by which human leukocyte antigen B27 (HLA-B27) contributes to ankylosing spondylitis (AS) remains unclear. Genetic studies demonstrate that association with and interaction between polymorphisms of endoplasmic reticulum aminopeptidase 1 (ERAP1) and HLA-B27 influence the risk of AS. It has been hypothesised that ERAP1-mediated HLA-B27 misfolding increases endoplasmic reticulum (ER) stress, driving an interleukin (IL) 23-dependent, pro-inflammatory immune response. We tested the hypothesis that AS-risk ERAP1 variants increase ER-stress and concomitant pro-inflammatory cytokine production in HLA-B27 + but not HLA-B27-AS patients or controls. Forty-nine AS cases and 22 healthy controls were grouped according to HLA-B27 status and AS-associated ERAP1 rs30187 genotypes: HLA-B27 + ERAP1 risk, HLA-B27 + ERAP1 protective, HLA-B27-ERAP1 risk and HLA-B27-ERAP1 protective. Expression levels of ER-stress markers GRP78 (8 kDa glucose-regulated protein), CHOP (C/EBP-homologous protein) and inflammatory cytokines were determined in peripheral blood mononuclear cell and ileal biopsies. We found no differences in ER-stress gene expression between HLA-B27 + and HLA-B27-cases or healthy controls, or between cases or controls stratified by carriage of ERAP1 risk or protective alleles in the presence or absence of HLA-B27. No differences were observed between expression of IL17A or TNF (tumour necrosis factor) in HLA-B27 + ERAP1 risk, HLA-B27 + ERAP1 protective and HLA-B27-ERAP1 protective cases. These data demonstrate that aberrant ERAP1 activity and HLA-B27 carriage does not alter ER-stress levels in AS, suggesting that ERAP1 and HLA-B27 may influence disease susceptibility through other mechanisms. © 2015 Macmillan Publishers Limited.

Endoplasmic reticulum aminopeptidases in the pathogenesis of ankylosing spondylitis

There has been significant progress in our understanding of the pathogenesis of AS. The advent of genome-wide association studies has increased the known loci associated with AS to more than 40. The endoplasmic reticulum resident aminopeptidases (ERAP) 1 and 2 were identified in this manner and are of particular interest. There appears to be a genetic as well as a functional interaction of ERAP1 and 2 with HLA-B27 based on the known functions of these molecules. Recent studies on the structure, immunological effects and the peptide-trimming properties of ERAP 1 and 2 have helped to provide insight into their pathogenic potential in AS. In this review, we explore the role of ERAP 1 and 2 in the pathogenesis of AS. © The Author 2015.

ERAP2 functional knockout in humans does not alter surface heavy chains or HLA-B27, inflammatory cytokines or endoplasmic reticulum stress markers

Introduction Single nucleotide polymorphisms in ERAP2 are strongly associated with ankylosing spondylitis (AS). One AS-associated single nucleotide polymorphism, rs2248374, causes a truncated ERAP2 protein that is degraded by nonsense-mediated decay. Approximately 25% of the populations of European ancestry are therefore natural ERAP2 knockouts. We investigated the effect of this associated variant on HLA class I allele presentation, surface heavy chains, endoplasmic reticulum (ER) stress markers and cytokine gene transcription in AS. Methods Patients with AS and healthy controls with either AA or GG homozygous status for rs2248374 were studied. Antibodies to CD14, CD19-ECD, HLA-A-B-C, Valpha7.2, CD161, anti-HC10 and anti-HLA-B27 were used to analyse peripheral blood mononuclear cells. Expression levels of ER stress markers (GRP78 and CHOP) and proinflammatory genes (tumour necrosis factor (TNF), IL6, IL17 and IL22) were assessed by qPCR. Results There was no significant difference in HLAclass I allele presentation or major histocompatibility class I heavy chains or ER stress markers GRP78 and CHOP or proinflammatory gene expression between genotypes for rs2248374 either between cases, between cases and controls, and between controls. Discussion Large differences were not seen in HLAB27 expression or cytokine levels between subjects with and without ERAP2 in AS cases and controls. This suggests that ERAP2 is more likely to influence AS risk through other mechanisms.

Association of FOXE1 Polyalanine Repeat Region with Papillary Thyroid Cancer

CONTEXT: Polyalanine tract variations in transcription factors have been identified for a wide spectrum of developmental disorders. The thyroid transcription factor forkhead factor E1 (FOXE1) contains a polymorphic polyalanine tract with 12-22 alanines. Single-nucleotide polymorphisms (SNP) close to this locus are associated with papillary thyroid cancer (PTC), and a strong linkage disequilibrium block extends across this region. OBJECTIVE: The objective of the study was to assess whether the FOXE1 polyalanine repeat region was associated with PTC and to assess the effect of polyalanine repeat region variants on protein expression, DNA binding, and transcriptional function on FOXE1-responsive promoters. DESIGN: This was a case-control study. SETTING: The study was conducted at a tertiary referral hospital. PATIENTS AND METHODS: The FOXE1 polyalanine repeat region and tag SNP were genotyped in 70 PTC, with a replication in a further 92 PTC, and compared with genotypes in 5767 healthy controls (including 5667 samples from the Wellcome Trust Case Control Consortium). In vitro studies were performed to examine the protein expression, DNA binding, and transcriptional function for FOXE1 variants of different polyalanine tract lengths. RESULTS: All the genotyped SNP were in tight linkage disequilibrium, including the FOXE1 polyalanine repeat region. We confirmed the strong association of rs1867277 with PTC (overall P = 1 × 10(-7), odds ratio 1.84, confidence interval 1.31-2.57). rs1867277 was in tight linkage disequilibrium with the FOXE1 polyalanine repeat region (r(2) = 0.95). FOXE1(16Ala) was associated with PTC with an odds ratio of 2.23 (confidence interval 1.42-3.50; P = 0.0005). Functional studies in vitro showed that FOXE1(16Ala) was transcriptionally impaired compared with FOXE1(14Ala), which was not due to differences in protein expression or DNA binding. CONCLUSIONS: We have confirmed the previous association of FOXE1 with PTC. Our data suggest that the coding polyalanine expansion in FOXE1 may be responsible for the observed association between FOXE1 and PTC.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10-8) level: 14q24.2 (rs227425, P-value 3.98 × 10-13, SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10-9, CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n 5 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis. © The Author 2013. Published by Oxford University Press.

Novel TBK1 truncating mutation in a familial amyotrophic lateral sclerosis patient of Chinese origin

Missense and frameshift mutations in TRAF family member-associated NF-kappa-B activator (TANK)-binding kinase 1 (TBK1) have been reported in European sporadic and familial amyotrophic lateral sclerosis (ALS) cohorts. To assess the role of TBK1 in ALS patient cohorts of wider ancestry, we have analyzed whole-exome sequence data from an Australian cohort of familial ALS (FALS) patients and controls. We identified a novel TBK1 deletion (c.1197delC) in a FALS patient of Chinese origin. This frameshift mutation (p.L399fs) likely results in a truncated protein that lacks functional domains required for adapter protein binding, as well as protein activation and structural integrity. No novel or reported TBK1 mutations were identified in FALS patients of European ancestry. This is the first report of a TBK1 mutation in an ALS patient of Asian origin and indicates that sequence variations in TBK1 are a rare cause of FALS in Australia. © 2015 Elsevier Inc.

Whole-exome re-sequencing in a family quartet identifies POP1 mutations as the cause of a novel skeletal dysplasia

Recent advances in DNA sequencing have enabled mapping of genes for monogenic traits in families with small pedigrees and even in unrelated cases. We report the identification of disease-causing mutations in a rare, severe, skeletal dysplasia, studying a family of two healthy unrelated parents and two affected children using whole-exome sequencing. The two affected daughters have clinical and radiographic features suggestive of anauxetic dysplasia (OMIM 607095), a rare form of dwarfism caused by mutations of RMRP. However, mutations of RMRP were excluded in this family by direct sequencing. Our studies identified two novel compound heterozygous loss-of-function mutations in POP1, which encodes a core component of the RNase mitochondrial RNA processing (RNase MRP) complex that directly interacts with the RMRP RNA domains that are affected in anauxetic dysplasia. We demonstrate that these mutations impair the integrity and activity of this complex and that they impair cell proliferation, providing likely molecular and cellular mechanisms by which POP1 mutations cause this severe skeletal dysplasia. © 2011 Glazov et al.

Mutational analysis of the insulin‑like growth factor 1 receptor tyrosine kinase domain in non-small cell lung cancer patients

The insulin‑like growth factor 1 receptor (IGF1R) pathway plays an important role in the pathogenesis of non‑small cell lung cancer (NSCLC) and also provides a mechanism of resistance to targeted therapies. IGF1R is therefore an ideal therapeutic target and several inhibitors have entered clinical trials. However, thus far the response to these inhibitors has been poor, highlighting the importance of predictive biomarkers to identify patient cohorts who will benefit from these targeted agents. It is well‑documented that mutations and/or deletions in the epidermal growth factor receptor (EGFR) tyrosine kinase (TK) domain predict sensitivity of NSCLC patients to EGFR TK inhibitors. Single‑nucleotide polymorphisms (SNPs) in the IGF pathway have been associated with disease, including breast and prostate cancer. The aim of the present study was to elucidate whether the IGF1R TK domain harbours SNPs, somatic mutations or deletions in NSCLC patients and correlates the mutation status to patient clinicopathological data and prognosis. Initially 100 NSCLC patients were screened for mutations/deletions in the IGF1R TK domain (exons 16‑21) by sequencing analysis. Following the identification of SNP rs2229765, a further 98 NSCLC patients and 866 healthy disease‑free control patients were genotyped using an SNP assay. The synonymous SNP (rs2229765) was the only aberrant base change identified in the IGF1R TK domain of 100 NSCLC patients initially analysed. SNP rs2229765 was detected in exon 16 and was found to have no significant association between IGF1R expression and survival. The GA genotype was identified in 53.5 and 49.4% of NSCLC patients and control individuals, respectively. No significant difference was found in the genotype (P=0.5487) or allele (P=0.9082) frequencies between the case and control group. The present findings indicate that in contrast to the EGFR TK domain, the IGF1R TK domain is not frequently mutated in NSCLC patients. The synonymous SNP (rs2229765) had no significant association between IGF1R expression and survival in the cohort of NSCLC patients.

Association of a Notch 3 gene polymorphism with migraine susceptibility

Introduction Cerebral autosomal dominant arteriopathy with subcortical infarcts and leucoencephalopathy (CADASIL) shares common symptoms with migraine. Most CADASIL causative mutations occur in exons 3 and 4 of the Notch 3 gene. This study investigated the role of C381T (rs 3815188) and G684A (rs 1043994) single nucleotide polymorphisms (SNP) in exons 3 and 4, respectively, of the Notch 3 gene in migraine. Results The first part of the study, in a population of 275 migraineurs and 275 control individuals, found a significant association between the C381T variant and migraine, specifically in migraine without aura (MO) sufferers. The G684A variant was also found to be significantly associated with migraine, specifically in migraine with aura (MA) sufferers. A follow-up study in 300 migraineurs and 300 control individuals did not show replicated association of the C381T variant with migraineurs. However, the G684A variant was again shown to be significantly associated with migraine, specifically with MA. Conclusion Further investigation of the G684A variant and the Notch 3 gene is warranted to understand their role in migraine.

Campylobacter jejuni and campylobacter coli genotyping by high-resolution melting analysis of a flaA fragment

The highly variable flagellin-encoding flaA gene has long been used for genotyping Campylobacter jejuni and Campylobacter coli. High-resolution melting (HRM) analysis is emerging as an efficient and robust method for discriminating DNA sequence variants. The objective of this study was to apply HRM analysis to flaA-based genotyping. The initial aim was to identify a suitable flaA fragment. It was found that the PCR primers commonly used to amplify the flaA short variable repeat (SVR) yielded a mixed PCR product unsuitable for HRM analysis. However, a PCR primer set composed of the upstream primer used to amplify the fragment used for flaA restriction fragment length polymorphism (RFLP) analysis and the downstream primer used for flaA SVR amplification generated a very pure PCR product, and this primer set was used for the remainder of the study. Eighty-seven C. jejuni and 15 C. coli isolates were analyzed by flaA HRM and also partial flaA sequencing. There were 47 flaA sequence variants, and all were resolved by HRM analysis. The isolates used had previously also been genotyped using single-nucleotide polymorphisms (SNPs), binary markers, CRISPR HRM, and flaA RFLP. flaAHRManalysis provided resolving power multiplicative to the SNPs, binary markers, and CRISPR HRM and largely concordant with the flaA RFLP. It was concluded that HRM analysis is a promising approach to genotyping based on highly variable genes.

HapMap tag-SNP analysis confirms a role for COMT in schizophrenia risk and reveals a novel association

Catechol-O-methyl transferase (COMT) encodes an enzyme involved in the metabolism of dopamine and maps to a commonly deleted region that increases schizophrenia risk. A non-synonymous polymorphism (rs4680) in COMT has been previously found to be associated with schizophrenia and results in altered activity levels of COMT. Using a haplotype block-based gene-tagging approach we conducted an association study of seven COMT single nucleotide polymorphisms (SNPs) in 160 patients with a DSM-IV diagnosis of schizophrenia and 250 controls in an Australian population. Two polymorphisms including rs4680 and rs165774 were found to be significantly associated with schizophrenia. The rs4680 results in a Val/Met substitution but the strongest association was shown by the novel SNP, rs165774, which may still be functional even though it is located in intron five. Individuals with schizophrenia were more than twice as likely to carry the GG genotype compared to the AA genotype for both the rs165774 and rs4680 SNPs. This association was slightly improved when males were analysed separately possibly indicating a degree of sexual dimorphism. Our results confirm that COMT is a good candidate for schizophrenia risk, by replicating the association with rs4680 and identifying a novel SNP association.

Association between polymorphisms in the progesterone receptor gene and endometriosis

The progesterone receptor (PR) is a candidate gene for the development of endometriosis, a complex disease with strong hormonal features, common in women of reproductive age. We typed the 306 base pair Alu insertion (AluIns) polymorphism in intron G of PR in 101 individuals, estimated linkage disequilibrium (LD) between five single-nucleotide polymorphisms (SNPs) across the PR locus in 980 Australian triads (endometriosis case and two parents) and used transmission disequilibrium testing (TDT) for association with endometriosis. The five SNPs showed strong pairwise LD, and the AluIns was highly correlated with proximal SNPs rs1042839 (Δ2 = 0.877, D9 = 1.00, P < 0.0001) and rs500760 (Δ2 = 0.438, D9 = 0.942, P < 0.0001). TDT showed weak evidence of allelic association between endometriosis and rs500760 (P = 0.027) but not in the expected direction. We identified a common susceptibility haplotype GGGCA across the five SNPs (P = 0.0167) in the whole sample, but likelihood ratio testing of haplotype transmission and non-transmission of the AluIns and flanking SNPs showed no significant pattern. Further, analysis of our results pooled with those from two previous studies suggested that neither the T2 allele of the AluIns nor the T1/T2 genotype was associated with endometriosis.

A three-nucleotide mutation altering the Maize streak virus Rep pRBR-interaction motif reduces symptom severity in maize and partially reverts at high frequency without restoring pRBR-Rep binding

Geminivirus infectivity is thought to depend on interactions between the virus replication-associated proteins Rep or RepA and host retinoblastoma-related proteins (pRBR), which control cell-cycle progression. It was determined that the substitution of two amino acids in the Maize streak virus (MSV) RepA pRBR-interaction motif (LLCNE to LLCLK) abolished detectable RepA-pRBR interaction in yeast without abolishing infectivity in maize. Although the mutant virus was infectious in maize, it induced less severe symptoms than the wild-type virus. Sequence analysis of progeny viral DNA isolated from infected maize enabled detection of a high-frequency single-nucleotide reversion of C(601)A in the 3 nt mutated sequence of the Rep gene. Although it did not restore RepA-pRBR interaction in yeast, sequence-specific PCR showed that, in five out of eight plants, the C(601)A reversion appeared by day 10 post-inoculation. In all plants, the C(601)A revertant eventually completely replaced the original mutant population, indicating a high selection pressure for the single-nucleotide reversion. Apart from potentially revealing an alternative or possibly additional function for the stretch of DNA that encodes the apparently non-essential pRBR-interaction motif of MSV Rep, the consistent emergence and eventual dominance of the C(601)A revertant population might provide a useful tool for investigating aspects of MSV biology, such as replication, mutation and evolution rates, and complex population phenomena, such as competition between quasispecies and population turnover. © 2005 SGM.