Variant STAT4 and Response to Ruxolitinib in an Autoinflammatory Syndrome

BACKGROUND

Disabling pansclerotic morphea (DPM) is a rare systemic inflammatory disorder, characterized by poor wound healing, fibrosis, cytopenias, hypogammaglobulinemia, and squamous-cell carcinoma. The cause is unknown, and mortality is high.

METHODS

We evaluated four patients from three unrelated families with an autosomal dominant pattern of inheritance of DPM. Genomic sequencing independently identified three heterozygous variants in a specific region of the gene that encodes signal transducer and activator of transcription 4 (STAT4). Primary skin fibroblast and cell-line assays were used to define the functional nature of the genetic defect. We also assayed gene expression using single-cell RNA sequencing of peripheral-blood mononuclear cells to identify inflammatory pathways that may be affected in DPM and that may respond to therapy.

RESULTS

Genome sequencing revealed three novel heterozygous missense gain-of-function variants in STAT4. In vitro, primary skin fibroblasts showed enhanced interleukin-6 secretion, with impaired wound healing, contraction of the collagen matrix, and matrix secretion. Inhibition of Janus kinase (JAK)-STAT signaling with ruxolitinib led to improvement in the hyperinflammatory fibroblast phenotype in vitro and resolution of inflammatory markers and clinical symptoms in treated patients, without adverse effects. Single-cell RNA sequencing revealed expression patterns consistent with an immunodysregulatory phenotype that were appropriately modified through JAK inhibition.

CONCLUSIONS

Gain-of-function variants in STAT4 caused DPM in the families that we studied. The JAK inhibitor ruxolitinib attenuated the dermatologic and inflammatory phenotype in vitro and in the affected family members. (Funded by the American Academy of Allergy, Asthma, and Immunology Foundation and others.).

Severe delayed hypersensitivity reactions to IL-1 and IL-6 inhibitors link to common HLA-DRB1*15 alleles

OBJECTIVES

Drug reaction with eosinophilia and systemic symptoms (DRESS) is a severe, delayed hypersensitivity reaction (DHR). We observed DRESS to inhibitors of interleukin 1 (IL-1) or IL-6 in a small group of patients with Still's disease with atypical lung disease. We sought to characterise features of patients with Still's disease with DRESS compared with drug-tolerant Still's controls. We analysed human leucocyte antigen (HLA) alleles for association to inhibitor-related DHR, including in a small Kawasaki disease (KD) cohort.

METHODS

In a case/control study, we collected a multicentre series of patients with Still's disease with features of inhibitor-related DRESS (n=66) and drug-tolerant Still's controls (n=65). We retrospectively analysed clinical data from all Still's subjects and typed 94/131 for HLA. European Still's-DRESS cases were ancestry matched to International Childhood Arthritis Genetics Consortium paediatric Still's cases (n=550) and compared for HLA allele frequencies. HLA association also was analysed using Still's-DRESS cases (n=64) compared with drug-tolerant Still's controls (n=30). KD subjects (n=19) were similarly studied.

RESULTS

Still's-DRESS features included eosinophilia (89%), AST-ALT elevation (75%) and non-evanescent rash (95%; 88% involving face). Macrophage activation syndrome during treatment was frequent in Still's-DRESS (64%) versus drug-tolerant Still's (3%; p=1.2×10-14). We found striking enrichment for HLA-DRB1*15 haplotypes in Still's-DRESS cases versus INCHARGE Still's controls (p=7.5×10-13) and versus self-identified, ancestry-matched Still's controls (p=6.3×10-10). In the KD cohort, DRB1*15:01 was present only in those with suspected anakinra reactions.

CONCLUSIONS

DRESS-type reactions occur among patients treated with IL-1/IL-6 inhibitors and strongly associate with common HLA-DRB1*15 haplotypes. Consideration of preprescription HLA typing and vigilance for serious reactions to these drugs are warranted.

Sequence-Based Screening of Patients With Idiopathic Polyarteritis Nodosa, Granulomatosis With Polyangiitis, and Microscopic Polyangiitis for Deleterious Genetic Variants in ADA2

OBJECTIVE

Deficiency of adenosine deaminase 2 (DADA2) is a monogenic form of vasculitis that can resemble polyarteritis nodosa (PAN). This study was undertaken to identify potential disease-causing sequence variants in ADA2 in patients with idiopathic PAN, granulomatosis with polyangiitis (GPA), or microscopic polyangiitis (MPA).

METHODS

Patients with idiopathic PAN (n = 118) and patients with GPA or MPA (n = 1,107) were screened for rare nonsynonymous variants in ADA2 using DNA sequencing methods. ADA-2 enzyme activity was assessed in selected serum samples.

RESULTS

Nine of 118 patients with PAN (7.6%) were identified as having rare nonsynonymous variants in ADA2. Four patients (3.4%) were biallelic for pathogenic or likely pathogenic variants, and 5 patients (4.2%) were monoallelic carriers for 3 variants of uncertain significance and 2 likely pathogenic variants. Serum samples from 2 patients with PAN with biallelic variants were available and showed markedly reduced ADA-2 enzyme activity. ADA-2 enzyme testing of 86 additional patients revealed 1 individual with strongly reduced ADA-2 activity without detectable pathogenic variants. Patients with PAN and biallelic variants in ADA2 were younger at diagnosis than patients with 1 or no variant in ADA2, with no other clinical differences noted. None of the patients with GPA or MPA carried biallelic variants in ADA2.

CONCLUSION

A subset of patients with idiopathic PAN meet genetic criteria for DADA2. Given that tumor necrosis factor inhibition is efficacious in DADA2 but is not conventional therapy for PAN, these findings suggest that ADA-2 testing should strongly be considered in patients with hepatitis B virus-negative idiopathic PAN.

Genetic control of CCL24, POR, and IL23R contributes to the pathogenesis of sarcoidosis

Sarcoidosis is a genetically complex systemic inflammatory disease that affects multiple organs. We present a GWAS of a Japanese cohort (700 sarcoidosis cases and 886 controls) with replication in independent samples from Japan (931 cases and 1,042 controls) and the Czech Republic (265 cases and 264 controls). We identified three loci outside the HLA complex, CCL24, STYXL1-SRRM3, and C1orf141-IL23R, which showed genome-wide significant associations (P < 5.0 × 10⁻⁸) with sarcoidosis; CCL24 and STYXL1-SRRM3 were novel. The disease-risk alleles in CCL24 and IL23R were associated with reduced CCL24 and IL23R expression, respectively. The disease-risk allele in STYXL1-SRRM3 was associated with elevated POR expression. These results suggest that genetic control of CCL24, POR, and IL23R expression contribute to the pathogenesis of sarcoidosis. We speculate that the CCL24 risk allele might be involved in a polarized Th1 response in sarcoidosis, and that POR and IL23R risk alleles may lead to diminished host defense against sarcoidosis pathogens.

Akira Meguro et al. report a genome-wide association study for sarcoidosis—a systemic inflammatory disease—in the Japanese population. They identify 3 non-HLA loci with genome-wide significance, 2 of which have not been previously associated with sarcoidosis in any population.

Common genetic susceptibility loci link PFAPA syndrome, Behçet's disease, and recurrent aphthous stomatitis

Periodic fever, aphthous stomatitis, pharyngitis, and cervical adenitis (PFAPA) syndrome is the most common periodic fever syndrome in children. The disease appears to cluster in families, but the pathogenesis is unknown. We queried two European-American cohorts and one Turkish cohort (total n = 231) of individuals with PFAPA for common variants previously associated with two other oropharyngeal ulcerative disorders, Behçet's disease and recurrent aphthous stomatitis. In a metaanalysis, we found that a variant upstream of IL12A (rs17753641) is strongly associated with PFAPA (OR 2.13, P = 6 × 10-9). We demonstrated that monocytes from individuals who are heterozygous or homozygous for this risk allele produce significantly higher levels of IL-12p70 upon IFN-γ and LPS stimulation than those from individuals without the risk allele. We also found that variants near STAT4, IL10, and CCR1-CCR3 were significant susceptibility loci for PFAPA, suggesting that the pathogenesis of PFAPA involves abnormal antigen-presenting cell function and T cell activity and polarization, thereby implicating both innate and adaptive immune responses at the oropharyngeal mucosa. Our results illustrate genetic similarities among recurrent aphthous stomatitis, PFAPA, and Behçet's disease, placing these disorders on a common spectrum, with recurrent aphthous stomatitis on the mild end, Behçet's disease on the severe end, and PFAPA intermediate. We propose naming these disorders Behçet's spectrum disorders to highlight their relationship. HLA alleles may be factors that influence phenotypes along this spectrum as we found new class I and II HLA associations for PFAPA distinct from Behçet's disease and recurrent aphthous stomatitis.

Pseudodominance of autoinflammatory disease in a single Turkish family explained by co-inheritance of haploinsufficiency of A20 and familial Mediterranean fever

OBJECTIVES

We investigated a Turkish family with multiple patients presenting with familial Mediterranean fever (FMF) and Behçet's disease (BD)-like manifestations. The index case and the two daughters with Behçet-like disease, were previously found to have a TNFAIP3 frameshift mutation. The high number of affected cases in this expanded family could be consistent with a dominantly inherited inflammatory disease, although some individuals had clinical features more consistent with recessively inherited FMF. We sequenced DNA from members of this family to determine whether the TNFAIP3 frameshift mutation and/or MEFV variants could explain this autoinflammatory disease pedigree.

METHODS

Patients were clinically diagnosed to have FMF or BD. Sanger sequence targeting TNFAIP3 exon 5 and MEFV exon 10 was carried out.

RESULTS

The symptomatic mother of the index case and her affected maternal uncle had compound heterozygous FMF-associated MEFV mutations, p.Met680Ile and p.Arg761His. Two affected daughters of the maternal uncle also had compound heterozygous FMF-associated mutations, p.Met680Ile and p.Val726Ala. The index case and her two affected daughters had a TNFAIP3 frameshift mutation (c.799delG; p.Pro268Leufs*19), which is consistent with their HA20 diagnosis, and also carried a heterozygous MEFV p.Arg761His mutation.

CONCLUSIONS

Autoinflammatory disease manifestations in a Turkish family with multiple affected cases could be explained by co-inheritance of pathogenic MEFV variants and a heterozygous HA20-associated mutation. FMF-associated p.Arg761His allele carried with the loss of function TNFAIP3 mutation by all three HA20 patients may contribute to their autoinflammatory phenotype and could also be responsible for their favourable response to colchicine.

Brief Report: Whole-Exome Sequencing to Identify Rare Variants and Gene Networks That Increase Susceptibility to Scleroderma in African Americans

OBJECTIVE

Whole-exome sequencing (WES) studies in systemic sclerosis (SSc) patients of European American (EA) ancestry have identified variants in the ATP8B4 gene and enrichment of variants in genes in the extracellular matrix (ECM)-related pathway that increase SSc susceptibility. This study was undertaken to evaluate the association of the ATP8B4 gene and the ECM-related pathway with SSc in a cohort of African American (AA) patients.

METHODS

SSc patients of AA ancestry were enrolled from 23 academic centers across the US under the Genome Research in African American Scleroderma Patients consortium. Unrelated AA individuals without serologic evidence of autoimmunity who were enrolled in the Howard University Family Study were used as unaffected controls. Functional variants in genes reported in the 2 WES studies in EA patients with SSc were selected for gene association testing using the optimized sequence kernel association test (SKAT-O) and pathway analysis by Ingenuity Pathway Analysis in 379 patients and 411 controls.

RESULTS

Principal components analysis demonstrated that the patients and controls had similar ancestral backgrounds, with roughly equal proportions of mean European admixture. Using SKAT-O, we examined the association of individual genes that were previously reported in EA patients and none remained significant, including ATP8B4 (P = 0.98). However, we confirmed the previously reported association of the ECM-related pathway with enrichment of variants within the COL13A1, COL18A1, COL22A1, COL4A3, COL4A4, COL5A2, PROK1, and SERPINE1 genes (corrected P = 1.95 × 10^-4 ).

CONCLUSION

In the largest genetic study in AA patients with SSc to date, our findings corroborate the role of functional variants that aggregate in a fibrotic pathway and increase SSc susceptibility.

IL1RN Variation Influences Both Disease Susceptibility and Response to Recombinant Human Interleukin-1 Receptor Antagonist Therapy in Systemic Juvenile Idiopathic Arthritis

OBJECTIVE

To determine whether systemic juvenile idiopathic arthritis (JIA) susceptibility loci that were identified by candidate gene studies demonstrate association with systemic JIA in the largest study population assembled to date.

METHODS

Single-nucleotide polymorphisms (SNPs) from 11 previously reported systemic JIA risk loci were examined for association in 9 populations, including 770 patients with systemic JIA and 6,947 controls. The effect of systemic JIA-associated SNPs on gene expression was evaluated in silico in paired whole genome and RNA sequencing data from the lymphoblastoid cell lines (LCLs) of 373 European subjects from the 1000 Genomes Project. Responses of systemic JIA-associated SNPs to anakinra treatment were evaluated in 38 US patients for whom treatment response data were available.

RESULTS

We found no association between the previously reported 26 SNPs and systemic JIA. Expanded analysis of the regions containing the 26 SNPs revealed only 1 significant association: the promoter region of IL1RN (P < 1 × 10-4 ). Systemic JIA-associated SNPs correlated with IL1RN expression in LCLs, with an inverse correlation between systemic JIA risk and IL1RN expression. The presence of homozygous IL1RN high expression alleles correlated strongly with a lack of response to anakinra therapy (odds ratio 28.7 [95% confidence interval 3.2-255.8]).

CONCLUSION

In our study, IL1RN was the only candidate locus associated with systemic JIA. The implicated SNPs are among the strongest known determinants of IL1RN and interleukin-1 receptor antagonist levels, linking low expression with increased systemic JIA risk. Homozygous high expression alleles predicted nonresponsiveness to anakinra therapy, making them ideal candidate biomarkers to guide systemic JIA treatment. This study is an important first step toward the personalized treatment of systemic JIA.

Three variants in the nicotinamide adenine dinucleotide phosphate oxidase complex are associated with HCV-related liver damage

Approximately 71 million people are chronically infected with the hepatitis C virus (HCV), a potentially lethal pathogen. HCV generates oxidative stress correlating with disease severity. HCV proteins increase reactive oxygen species production by stimulating nicotinamide adenine dinucleotide phosphate oxidase (NOX) activity. Reactive oxygen species are necessary for host defense and cell signaling; however, elevated NOX activity contributes to cancer, and NOX overexpression is associated with hepatic fibrosis. Our aim was to investigate whether single nucleotide polymorphisms (SNPs) in NOX family members are associated with HCV-related liver damage. Three hundred and thirty-one individuals of European ancestry and 90 individuals of African ancestry, all diagnosed with HCV, were genotyped for 243 tagSNPs in NOX enzymes and their regulatory factors. Pathology scores were available for 288 Caucasians and 71 Africans, and mortality status was determined for all subjects. SNPs were tested for association with pathology scores and as predictors of mortality. In Africans, homozygosity for the A allele of rs12753665 (neutrophil cytosolic factor 2) and homozygosity for the T allele of rs760519 (neutrophil cytosolic factor 4) were associated with and predictive of higher rates of advanced fibrosis and cirrhosis compared to other genotypes after controlling for age and sex. In Caucasians, homozygosity for the T allele of rs2292464 (dual oxidase 1) was associated with and predictive of decreased periportal inflammation after controlling for age and sex. No SNPs were significant predictors of mortality. Conclusion: In this exploratory study, three NOX-related polymorphisms in two ethnic groups were significantly associated with hepatic inflammation and fibrosis. Future studies investigating these SNPs in larger cohorts of patients with HCV are warranted. (Hepatology Communications 2017;1:973-982).

Aberrant tRNA processing causes an autoinflammatory syndrome responsive to TNF inhibitors

OBJECTIVES

To characterise the clinical features, immune manifestations and molecular mechanisms in a recently described autoinflammatory disease caused by mutations in TRNT1, a tRNA processing enzyme, and to explore the use of cytokine inhibitors in suppressing the inflammatory phenotype.

METHODS

We studied nine patients with biallelic mutations in TRNT1 and the syndrome of congenital sideroblastic anaemia with immunodeficiency, fevers and developmental delay (SIFD). Genetic studies included whole exome sequencing (WES) and candidate gene screening. Patients' primary cells were used for deep RNA and tRNA sequencing, cytokine profiling, immunophenotyping, immunoblotting and electron microscopy (EM).

RESULTS

We identified eight mutations in these nine patients, three of which have not been previously associated with SIFD. Three patients died in early childhood. Inflammatory cytokines, mainly interleukin (IL)-6, interferon gamma (IFN-γ) and IFN-induced cytokines were elevated in the serum, whereas tumour necrosis factor (TNF) and IL-1β were present in tissue biopsies of patients with active inflammatory disease. Deep tRNA sequencing of patients' fibroblasts showed significant deficiency of mature cytosolic tRNAs. EM of bone marrow and skin biopsy samples revealed striking abnormalities across all cell types and a mix of necrotic and normal-appearing cells. By immunoprecipitation, we found evidence for dysregulation in protein clearance pathways. In 4/4 patients, treatment with a TNF inhibitor suppressed inflammation, reduced the need for blood transfusions and improved growth.

CONCLUSIONS

Mutations of TRNT1 lead to a severe and often fatal syndrome, linking protein homeostasis and autoinflammation. Molecular diagnosis in early life will be crucial for initiating anti-TNF therapy, which might prevent some of the severe disease consequences.

Clinical and serological features of systemic sclerosis in a multicenter African American cohort: Analysis of the genome research in African American scleroderma patients clinical database

Racial differences exist in the severity of systemic sclerosis (SSc). To enhance our knowledge about SSc in African Americans, we established a comprehensive clinical database from the largest multicenter cohort of African American SSc patients assembled to date (the Genome Research in African American Scleroderma Patients (GRASP) cohort).African American SSc patients were enrolled retrospectively and prospectively over a 30-year period (1987-2016), from 18 academic centers throughout the United States. The cross-sectional prevalence of sociodemographic, clinical, and serological features was evaluated. Factors associated with clinically significant manifestations of SSc were assessed using multivariate logistic regression analyses.The study population included a total of 1009 African American SSc patients, comprised of 84% women. In total, 945 (94%) patients met the 2013 American College of Rheumatology/European League Against Rheumatism (ACR/EULAR) classification criteria for SSc, with the remaining 64 (6%) meeting the 1980 ACR or CREST (calcinosis, Raynaud's phenomenon, esophageal dysmotility, sclerodactyly, telangiectasia) criteria. While 43% were actively employed, 33% required disability support. The majority (57%) had the more severe diffuse subtype and a young age at symptom onset (39.1 ± 13.7 years), in marked contrast to that reported in cohorts of predominantly European ancestry. Also, 1 in 10 patients had a severe Medsger cardiac score of 4. Pulmonary fibrosis evident on computed tomography (CT) chest was present in 43% of patients and was significantly associated with anti-topoisomerase I positivity. 38% of patients with CT evidence of pulmonary fibrosis had a severe restrictive ventilator defect, forced vital capacity (FVC) ≤50% predicted. A significant association was noted between longer disease duration and higher odds of pulmonary hypertension, telangiectasia, and calcinosis. The prevalence of potentially fatal scleroderma renal crisis was 7%, 3.5 times higher than the 2% prevalence reported in the European League Against Rheumatism Scleroderma Trials and Research (EUSTAR) cohort.Our study emphasizes the unique and severe disease burden of SSc in African Americans compared to those of European ancestry.

Brief Report: Deficiency of Complement 1r Subcomponent in Early-Onset Systemic Lupus Erythematosus: The Role of Disease-Modifying Alleles in a Monogenic Disease

OBJECTIVE

To identify a genetic cause of early-onset systemic lupus erythematosus (SLE) in a large consanguineous family from Turkey and to study the mechanisms of the disease.

METHODS

We performed whole-exome sequencing and single-nucleotide polymorphism array genotyping in family members with and without SLE. Protein and gene expression, cytokine profile, neutrophil extracellular trap (NET) formation, and presence of low-density granulocytes were evaluated in patient primary cells and serum samples.

RESULTS

We identified a novel, homozygous, loss-of-function mutation (p.Pro445Leufs*11) in the C1R gene. Using the Sanger method of DNA sequencing in 14 family members, we confirmed the presence of the mutation in 4 patients with SLE and in an asymptomatic 9-year-old girl. Complement levels were low in sera from patients with truncated C1r protein. Two siblings with SLE who were available for detailed evaluation exhibited strong type I interferon (IFN) inflammatory signatures despite their disease being clinically inactive at the time of sampling. The type I IFN transcriptional signature in the patients' blood correlated with disease expressivity, whereas the neutrophil signature in peripheral blood mononuclear cells was likely associated with disease severity. The female patient with SLE with the most severe phenotype presented with a stronger neutrophil signature, defined by enhanced NET formation and the presence of low-density granulocytes. Analysis of exome data for modifying alleles suggested enrichment of common SLE-associated variants in the more severely affected patients. Lupus-associated HLA alleles or HLA haplotypes were not shared among the 4 affected subjects.

CONCLUSION

Our findings revealed a novel high-penetrance mutation in C1R as the cause of monogenic SLE. Disease expressivity in this family appears to be influenced by additional common and rare genetic variants.

Genetic architecture distinguishes systemic juvenile idiopathic arthritis from other forms of juvenile idiopathic arthritis: clinical and therapeutic implications

OBJECTIVES

Juvenile idiopathic arthritis (JIA) is a heterogeneous group of conditions unified by the presence of chronic childhood arthritis without an identifiable cause. Systemic JIA (sJIA) is a rare form of JIA characterised by systemic inflammation. sJIA is distinguished from other forms of JIA by unique clinical features and treatment responses that are similar to autoinflammatory diseases. However, approximately half of children with sJIA develop destructive, long-standing arthritis that appears similar to other forms of JIA. Using genomic approaches, we sought to gain novel insights into the pathophysiology of sJIA and its relationship with other forms of JIA.

METHODS

We performed a genome-wide association study of 770 children with sJIA collected in nine countries by the International Childhood Arthritis Genetics Consortium. Single nucleotide polymorphisms were tested for association with sJIA. Weighted genetic risk scores were used to compare the genetic architecture of sJIA with other JIA subtypes.

RESULTS

The major histocompatibility complex locus and a locus on chromosome 1 each showed association with sJIA exceeding the threshold for genome-wide significance, while 23 other novel loci were suggestive of association with sJIA. Using a combination of genetic and statistical approaches, we found no evidence of shared genetic architecture between sJIA and other common JIA subtypes.

CONCLUSIONS

The lack of shared genetic risk factors between sJIA and other JIA subtypes supports the hypothesis that sJIA is a unique disease process and argues for a different classification framework. Research to improve sJIA therapy should target its unique genetics and specific pathophysiological pathways.

Dense genotyping of immune-related loci implicates host responses to microbial exposure in Behçet's disease susceptibility

We analyzed 1,900 Turkish Behçet’s disease cases and 1,779 controls genotyped with the Immunochip. The most significantly associated single nucleotide polymorphism (SNP) was rs1050502, a tag SNP for HLA-B*51. In the Turkish discovery set, we identified three novel loci, IL1A-IL1B, IRF8, and CEBPB-PTPN1, with genome-wide significance (P<5×10⁻⁸) by direct genotyping, and ADO-EGR2 by imputation. ADO-EGR2, IRF8, and CEBPB-PTPN1 replicated by genotyping 969 Iranian cases and 826 controls. Imputed data in 608 Japanese cases and 737 controls replicated ADO-EGR2 and IRF8 and meta-analysis additionally identified RIPK2 and LACC1. The disease-associated allele of rs4402765, the lead marker of the IL1A-IL1B locus, was associated with both decreased interleukin-1α and increased interleukin-1β production. ABO non-secretor genotypes of two ancestry-specific FUT2 SNPs showed strong disease association (P=5.89×10⁻¹⁵). Our findings extend shared susceptibility genes with Crohn’s disease and leprosy, and implicate mucosal factors and the innate immune response to microbial exposure in Behçet’s disease susceptibility.

Evaluation of KIR3DL1/KIR3DS1 polymorphism in Behçet's disease

The Behçet's disease (BD)-associated human leukocyte antigen (HLA) allele, HLA-B*51 (B*51), encodes a ligand for a pair of allelic killer immunoglobulin-like receptors (KIR) present on cytotoxic cells-KIR3DL1, which inhibits their cytotoxicity, and KIR3DS1, which activates their cytotoxic activity. We tested whether KIR-regulated mechanisms contribute to BD by testing for association of KIR3DL1/KIR3DS1 genotypes with disease in 1799 BD patients and 1710 healthy controls from Turkey, as well as in different subsets of individuals with HLA-type-defined ligands for the KIR3D receptors. HLA types were imputed from single nucleotide polymorphism genotypes determined with the Immunochip. The presence of inhibitory KIR3DL1 or activating KIR3DS1 alleles did not differ significantly between cases and controls (KIR3DL1: 92.9% vs 93.4%, P_dominant=0.55; KIR3DS1: 42.7% vs 41.0%, P_dominant=0.29). The KIR3DL1/KIR3DS1 alleles were also present at similar frequencies among cases and controls bearing HLA-B with a Bw4 motif; HLA-B with a Bw4 motif with isoleucine at position 80; and HLA-B*51. Our results suggest that pathogenic mechanisms associated with HLA-B*51 do not primarily involve differential interactions with KIR3DL1 and KIR3DS1 receptors. However, due to the complexity of this locus (that is, sequence variation and copy number variation), we cannot exclude a role for other types of KIR variation in the pathogenesis of BD.

Hypophosphatasia and the risk of atypical femur fractures: a case-control study

Background

Case reports have linked adult hypophosphatasia as a possible cause of atypical femur fractures (AFF) associated with bisphosphonate use. Adult hypophosphatasia is an asymptomatic genetic condition which results in low alkaline phosphatase and elevated pyridoxal phosphate. We conducted a case–control study to assess the role of hypophosphatasia and atypical femur fracture.

Methods

We recruited 13 control patients who took long term bisphosphonates without complication and 10 patients who sustained atypical femur fractures (mean bisphosphonate use, 9 years both cohorts). Patients underwent clinical exam and measurement of alkaline phosphatase and pyridoxal phosphate (PLP) levels. In addition, DNA was extracted and the ALPL gene was sequenced in both cohorts.

Results

Low alkaline phosphatase levels (<55 U/L) were seen in 5/10 AFF patients and 5/13 control patients. Two control patients demonstrated low alkaline phosphatase levels and elevated PLP. The alkaline phosphatase (ALPL) gene exons and intron splice sites were sequenced in the atypical femur fracture and control cohorts and no coding mutations were identified in any subjects. Atypical femur fracture patients demonstrated more varus hip alignment (p < 0.048) with no significant difference in mechanical axis.

Conclusions

We found no evidence of hypophosphatasia as a risk factor for atypical femur fractures. Laboratory findings of mildly low alkaline phosphatase activity were equally common in atypical and control cohorts and may be due to long term bisphosphonate use.

Trial registration

Clinicaltrials.gov number NCT01360099. Prospectively registered May 20, 2011. First patient enrolled June 14, 2011.

A single endoplasmic reticulum aminopeptidase-1 protein allotype is a strong risk factor for Behçet's disease in HLA-B*51 carriers

INTRODUCTION

Endoplasmic reticulum aminopeptidase-1 (ERAP1) protein is highly polymorphic with numerous missense amino acid variants. We sought to determine the naturally occurring ERAP1 protein allotypes and their contribution to Behçet's disease.

METHODS

Genotypes of all reported missense ERAP1 gene variants with 1000 Genomes Project EUR superpopulation frequency >1% were determined in 1900 Behçet's disease cases and 1779 controls from Turkey. ERAP1 protein allotypes and their contributions to Behçet's disease risk were determined by haplotype identification and disease association analyses.

RESULTS

One ERAP1 protein allotype with five non-ancestral amino acids was recessively associated with disease (p=3.13×10^-6, OR 2.55, 95% CI 1.70 to 3.82). The ERAP1 association was absent in individuals who lacked HLA-B*51. Individuals who carry HLA-B*51 and who are also homozygous for the haplotype had an increased disease odds compared with those with neither risk factor (p=4.80×10^-20, OR 10.96, 95% CI 5.91 to 20.32).

DISCUSSION

The Behçet's disease-associated ERAP1 protein allotype was previously shown to have poor peptide trimming activity. Combined with its requirement for HLA-B*51, these data suggest that a hypoactive ERAP1 allotype contributes to Behçet's disease risk by altering the peptides available for binding to HLA-B*51.

Additive loss-of-function proteasome subunit mutations in CANDLE/PRAAS patients promote type I IFN production

Autosomal recessive mutations in proteasome subunit β 8 (PSMB8), which encodes the inducible proteasome subunit β5i, cause the immune-dysregulatory disease chronic atypical neutrophilic dermatosis with lipodystrophy and elevated temperature (CANDLE), which is classified as a proteasome-associated autoinflammatory syndrome (PRAAS). Here, we identified 8 mutations in 4 proteasome genes, PSMA3 (encodes α7), PSMB4 (encodes β7), PSMB9 (encodes β1i), and proteasome maturation protein (POMP), that have not been previously associated with disease and 1 mutation in PSMB8 that has not been previously reported. One patient was compound heterozygous for PSMB4 mutations, 6 patients from 4 families were heterozygous for a missense mutation in 1 inducible proteasome subunit and a mutation in a constitutive proteasome subunit, and 1 patient was heterozygous for a POMP mutation, thus establishing a digenic and autosomal dominant inheritance pattern of PRAAS. Function evaluation revealed that these mutations variably affect transcription, protein expression, protein folding, proteasome assembly, and, ultimately, proteasome activity. Moreover, defects in proteasome formation and function were recapitulated by siRNA-mediated knockdown of the respective subunits in primary fibroblasts from healthy individuals. Patient-isolated hematopoietic and nonhematopoietic cells exhibited a strong IFN gene-expression signature, irrespective of genotype. Additionally, chemical proteasome inhibition or progressive depletion of proteasome subunit gene transcription with siRNA induced transcription of type I IFN genes in healthy control cells. Our results provide further insight into CANDLE genetics and link global proteasome dysfunction to increased type I IFN production.

A polymorphism in CCR1/CCR3 is associated with narcolepsy

Etiology of narcolepsy-cataplexy involves multiple genetic and environmental factors. While the human leukocyte antigen (HLA)-DRB1*15:01-DQB1*06:02 haplotype is strongly associated with narcolepsy, it is not sufficient for disease development. To identify additional, non-HLA susceptibility genes, we conducted a genome-wide association study (GWAS) using Japanese samples. An initial sample set comprising 409 cases and 1562 controls was used for the GWAS of 525,196 single nucleotide polymorphisms (SNPs) located outside the HLA region. An independent sample set comprising 240 cases and 869 controls was then genotyped at 37 SNPs identified in the GWAS. We found that narcolepsy was associated with a SNP in the promoter region of chemokine (C-C motif) receptor 1 (CCR1) (rs3181077, P=1.6×10(-5), odds ratio [OR]=1.86). This rs3181077 association was replicated with the independent sample set (P=0.032, OR=1.36). We measured mRNA levels of candidate genes in peripheral blood samples of 38 cases and 37 controls. CCR1 and CCR3 mRNA levels were significantly lower in patients than in healthy controls, and CCR1 mRNA levels were associated with rs3181077 genotypes. In vitro chemotaxis assays were also performed to measure monocyte migration. We observed that monocytes from carriers of the rs3181077 risk allele had lower migration indices with a CCR1 ligand. CCR1 and CCR3 are newly discovered susceptibility genes for narcolepsy. These results highlight the potential role of CCR genes in narcolepsy and support the hypothesis that patients with narcolepsy have impaired immune function.

The immunogenetics of Behçet's disease: A comprehensive review

Behçet's disease is a chronic multisystem inflammatory disorder characterized mainly by recurrent oral ulcers, ocular involvement, genital ulcers, and skin lesions, presenting with remissions and exacerbations. It is thought that both environmental and genetic factors contribute to its onset and development. Although the etiology of Behçet's disease remains unclear, recent immunogenetic findings are providing clues to its pathogenesis. In addition to the positive association of HLA-B*51, which was identified more than four decades ago, and which has since been confirmed in multiple populations, recent studies report additional independent associations in the major histocompatibility complex class I region. HLA-B*15, -B*27, -B*57, and -A*26 are independent risk factors for Behçet's disease, while HLA-B*49 and -A*03 are independent class I alleles that are protective for Behçet's disease. Genome-wide association studies have identified associations with genome-wide significance (P < 5 × 10(-8)) in the IL23R-IL12RB2, IL10, STAT4, CCR1-CCR3, KLRC4, ERAP1, TNFAIP3, and FUT2 loci. In addition, targeted next-generation sequencing has revealed the involvement of rare nonsynonymous variants of IL23R, TLR4, NOD2, and MEFV in Behçet's disease pathogenesis. Significant differences in gene function or mRNA expression associated with the risk alleles of the disease susceptibility loci suggest which genes in a disease-associated locus influence disease pathogenesis. These genes encompass both innate and adaptive immunity and confirm the importance of the predominant polarization towards helper T cell (Th) 1 versus Th2 cells, and the involvement of Th17 cells. In addition, epistasis observed between HLA-B*51 and the risk coding haplotype of the endoplasmic reticulum-associated protease, ERAP1, provides a clue that an HLA class I-peptide presentation-based mechanism contributes to this complex disease.

Brief Report: Cryopyrin-Associated Periodic Syndrome Caused by a Myeloid-Restricted Somatic NLRP3 Mutation

OBJECTIVE

To identify the cause of disease in an adult patient presenting with recent-onset fevers, chills, urticaria, fatigue, and profound myalgia, who was found to be negative for cryopyrin-associated periodic syndrome (CAPS) NLRP3 mutations by conventional Sanger DNA sequencing.

METHODS

We performed whole-exome sequencing and targeted deep sequencing using DNA from the patient's whole blood to identify a possible NLRP3 somatic mutation. We then screened for this mutation in subcloned NLRP3 amplicons from fibroblasts, buccal cells, granulocytes, negatively selected monocytes, and T and B lymphocytes and further confirmed the somatic mutation by targeted sequencing of exon 3.

RESULTS

We identified a previously reported CAPS-associated mutation, p.Tyr570Cys, with a mutant allele frequency of 15% based on exome data. Targeted sequencing and subcloning of NLRP3 amplicons confirmed the presence of the somatic mutation in whole blood at a ratio similar to the exome data. The mutant allele frequency was in the range of 13.3-16.8% in monocytes and 15.2-18% in granulocytes. Notably, this mutation was either absent or present at a very low frequency in B and T lymphocytes, in buccal cells, and in the patient's cultured fibroblasts.

CONCLUSION

Our findings indicate the possibility of myeloid-restricted somatic mosaicism in the pathogenesis of CAPS, underscoring the emerging role of massively parallel sequencing in clinical diagnosis.

Genetic architectures of seropositive and seronegative rheumatic diseases

Rheumatoid arthritis, systemic lupus erythematosus, ankylosing spondylitis and some other rheumatic diseases are genetically complex, with evidence of familial clustering, but not of Mendelian inheritance. These diseases are thought to result from contributions and interactions of multiple genetic and nongenetic risk factors, which have small effects individually. Genome-wide association studies (GWAS) of large collections of data from cases and controls have revealed many genetic factors that contribute to non-Mendelian rheumatic diseases, thus providing insights into associated molecular mechanisms. This Review summarizes methods for the identification of gene variants that influence genetically complex diseases and focuses on what we have learned about the rheumatic diseases for which GWAS have been reported. Our review of the disease-associated loci identified to date reveals greater sharing of risk loci among the groups of seropositive (diseases in which specific autoantibodies are often present) or seronegative diseases than between these two groups. The nature of the shared and discordant loci suggests important similarities and differences among these diseases.

Early-onset stroke and vasculopathy associated with mutations in ADA2

BACKGROUND

We observed a syndrome of intermittent fevers, early-onset lacunar strokes and other neurovascular manifestations, livedoid rash, hepatosplenomegaly, and systemic vasculopathy in three unrelated patients. We suspected a genetic cause because the disorder presented in early childhood.

METHODS

We performed whole-exome sequencing in the initial three patients and their unaffected parents and candidate-gene sequencing in three patients with a similar phenotype, as well as two young siblings with polyarteritis nodosa and one patient with small-vessel vasculitis. Enzyme assays, immunoblotting, immunohistochemical testing, flow cytometry, and cytokine profiling were performed on samples from the patients. To study protein function, we used morpholino-mediated knockdowns in zebrafish and short hairpin RNA knockdowns in U937 cells cultured with human dermal endothelial cells.

RESULTS

All nine patients carried recessively inherited mutations in CECR1 (cat eye syndrome chromosome region, candidate 1), encoding adenosine deaminase 2 (ADA2), that were predicted to be deleterious; these mutations were rare or absent in healthy controls. Six patients were compound heterozygous for eight CECR1 mutations, whereas the three patients with polyarteritis nodosa or small-vessel vasculitis were homozygous for the p.Gly47Arg mutation. Patients had a marked reduction in the levels of ADA2 and ADA2-specific enzyme activity in the blood. Skin, liver, and brain biopsies revealed vasculopathic changes characterized by compromised endothelial integrity, endothelial cellular activation, and inflammation. Knockdown of a zebrafish ADA2 homologue caused intracranial hemorrhages and neutropenia - phenotypes that were prevented by coinjection with nonmutated (but not with mutated) human CECR1. Monocytes from patients induced damage in cocultured endothelial-cell layers.

CONCLUSIONS

Loss-of-function mutations in CECR1 were associated with a spectrum of vascular and inflammatory phenotypes, ranging from early-onset recurrent stroke to systemic vasculopathy or vasculitis. (Funded by the National Institutes of Health Intramural Research Programs and others.).

Cyclic AMP and c-KIT signaling in familial testicular germ cell tumor predisposition

BACKGROUND

Familial testicular germ cell tumors (FTGCTs) are hypothesized to result from the combined interaction of multiple low-penetrance genes. We reported inactivating germline mutations of the cAMP-binding phosphodiesterase 11A (PDE11A) as modifiers of FTGCT risk. Recent genome-wide association studies have identified single-nucleotide polymorphisms in the KITLG gene, the ligand for the cKIT tyrosine kinase receptor, as strong modifiers of susceptibility to both familial and sporadic testicular germ cell tumors.

DESIGN

We studied 94 patients with FTGCTs and 50 at-risk male relatives from 63 unrelated kindreds, in whom the PDE11A gene had been sequenced by investigating the association between KITLG genome-wide association study single-nucleotide polymorphisms rs3782179 and rs4474514 and FTGCT risk in these patients and in 692 controls. We also examined cAMP and c-KIT signaling in testicular tissues and cell lines and extended the studies to 2 sporadic cases, one with a PDE11A defect and one without, as a comparison.

RESULTS

We found a higher frequency of the KITLG risk alleles in FTGCT patients who also had a PDE11A sequence variant, compared with those with a wild-type PDE11A sequence. In NTERA-2 and Tcam-2 cells transfected with the mutated forms of PDE11A (R52T, F258Y, Y727C, R804H, V820M, R867G, and M878V), cAMP levels were significantly higher, and the relative phosphodiesterase activity was lower than in the wild-type cells. KITLG expression was consistently increased in the presence of PDE11A-inactivating defects, both at the RNA and protein levels, in familial testicular germ cell tumors. The 2 sporadic cases that were studied, one with a PDE11A defect and another without, agreed with the data in FTGTCT and in the cell lines.

CONCLUSIONS

Patients with FTGCT and PDE11A defects also carry KITLG risk alleles more frequently. There may be an interaction between cAMP and c-KIT signaling in predisposition to testicular germ cell tumors.

Identification of two new arthritis severity loci that regulate levels of autoantibodies, interleukin-1β, and joint damage in pristane- and collagen-induced arthritis

OBJECTIVE

Cia3 is a locus on rat chromosome 4 that regulates severity and joint damage in collagen- and pristane-induced arthritis (CIA and PIA). This study was undertaken to refine the Cia3 gene-containing interval toward gene identification and obtain insights into its mode of action.

METHODS

Five DA.F344(Cia3) subcongenic rat strains were generated and studied using the PIA and CIA models. Levels of antibodies against type II collagen (both allo- and autoantibodies) were measured. Joints and synovial tissue were collected 32 days after the induction of PIA (chronic stage) for histologic and quantitative polymerase chain reaction analysis of interleukin-1β (IL-1β) and matrix metalloproteinase (MMP) levels.

RESULTS

Three subcongenic strains sharing the centromeric Cia3d interval were protected and 2 subcongenic strains sharing the telomeric Cia3g interval, which did not overlap with Cia3d, were also protected, developing significantly less severe CIA and PIA. Normal joint architecture was preserved in DA.F344(Cia3) and DA.F344(Cia3d) congenic rats with PIA, while DA rats had pronounced synovial hyperplasia, angiogenesis, inflammatory infiltration, and bone or cartilage erosions. The DA.F344(Cia3d) and DA.F344(Cia3g) strains had significantly lower synovial levels of IL-1β (5-fold and nearly 2-fold, respectively [the latter not reaching statistical significance]), MMP-1 (expressed predominantly in DA rats), MMP-3 (79-fold and 8-fold, respectively), and MMP-14 (21-fold and 1.4-fold, respectively) and reduced levels of pathogenic autoantibodies against type II collagen, compared with DA rats.

CONCLUSION

We have identified 2 new arthritis severity and articular damage loci within Cia3. These loci regulate pathogenic processes in 2 different models of rheumatoid arthritis, and the identification of these genes has the potential to generate new targets for therapies aimed at reducing disease severity and articular damage, and may additionally have prognostic value.

Cold urticaria, immunodeficiency, and autoimmunity related to PLCG2 deletions

BACKGROUND

Mendelian analysis of disorders of immune regulation can provide insight into molecular pathways associated with host defense and immune tolerance.

METHODS

We identified three families with a dominantly inherited complex of cold-induced urticaria, antibody deficiency, and susceptibility to infection and autoimmunity. Immunophenotyping methods included flow cytometry, analysis of serum immunoglobulins and autoantibodies, lymphocyte stimulation, and enzymatic assays. Genetic studies included linkage analysis, targeted Sanger sequencing, and next-generation whole-genome sequencing.

RESULTS

Cold urticaria occurred in all affected subjects. Other, variable manifestations included atopy, granulomatous rash, autoimmune thyroiditis, the presence of antinuclear antibodies, sinopulmonary infections, and common variable immunodeficiency. Levels of serum IgM and IgA and circulating natural killer cells and class-switched memory B cells were reduced. Linkage analysis showed a 7-Mb candidate interval on chromosome 16q in one family, overlapping by 3.5 Mb a disease-associated haplotype in a smaller family. This interval includes PLCG2, encoding phospholipase Cγ(2) (PLCγ(2)), a signaling molecule expressed in B cells, natural killer cells, and mast cells. Sequencing of complementary DNA revealed heterozygous transcripts lacking exon 19 in two families and lacking exons 20 through 22 in a third family. Genomic sequencing identified three distinct in-frame deletions that cosegregated with disease. These deletions, located within a region encoding an autoinhibitory domain, result in protein products with constitutive phospholipase activity. PLCG2-expressing cells had diminished cellular signaling at 37°C but enhanced signaling at subphysiologic temperatures.

CONCLUSIONS

Genomic deletions in PLCG2 cause gain of PLCγ(2) function, leading to signaling abnormalities in multiple leukocyte subsets and a phenotype encompassing both excessive and deficient immune function. (Funded by the National Institutes of Health Intramural Research Programs and others.).

A novel unstable duplication upstream of HAS2 predisposes to a breed-defining skin phenotype and a periodic fever syndrome in Chinese Shar-Pei dogs

Hereditary periodic fever syndromes are characterized by recurrent episodes of fever and inflammation with no known pathogenic or autoimmune cause. In humans, several genes have been implicated in this group of diseases, but the majority of cases remain unexplained. A similar periodic fever syndrome is relatively frequent in the Chinese Shar-Pei breed of dogs. In the western world, Shar-Pei have been strongly selected for a distinctive thick and heavily folded skin. In this study, a mutation affecting both these traits was identified. Using genome-wide SNP analysis of Shar-Pei and other breeds, the strongest signal of a breed-specific selective sweep was located on chromosome 13. The same region also harbored the strongest genome-wide association (GWA) signal for susceptibility to the periodic fever syndrome (p_raw = 2.3×10⁻⁶, p_genome = 0.01). Dense targeted resequencing revealed two partially overlapping duplications, 14.3 Kb and 16.1 Kb in size, unique to Shar-Pei and upstream of the Hyaluronic Acid Synthase 2 (HAS2) gene. HAS2 encodes the rate-limiting enzyme synthesizing hyaluronan (HA), a major component of the skin. HA is up-regulated and accumulates in the thickened skin of Shar-Pei. A high copy number of the 16.1 Kb duplication was associated with an increased expression of HAS2 as well as the periodic fever syndrome (p<0.0001). When fragmented, HA can act as a trigger of the innate immune system and stimulate sterile fever and inflammation. The strong selection for the skin phenotype therefore appears to enrich for a pleiotropic mutation predisposing these dogs to a periodic fever syndrome. The identification of HA as a major risk factor for this canine disease raises the potential of this glycosaminoglycan as a risk factor for human periodic fevers and as an important driver of chronic inflammation.

Shar-Pei dogs have two unique features: a breed defining “wrinkled” skin phenotype and a genetic disorder called Familial Shar-Pei Fever (FSF). The wrinkled phenotype is strongly selected for and is the result of excessive hyaluronan (HA) deposited in the skin. HA is a molecule that may behave in a pro-inflammatory manner and create a “danger signal” by being analogous to molecules on the surface of pathogens. FSF is characterized by unprovoked episodes of fever and/or inflammation and resembles several human autoinflammatory syndromes. Here we show that the two features are connected and have the same genetic origin, a regulatory mutation located close to a HA synthesizing gene (HAS2). The mutation is a 16.1 Kb duplication, the copy number of which correlates with HAS2 expression and disease. We suggest that the large amount of HA responsible for the skin condition predisposes to sterile fever and inflammation. HAS2 was previously not known to associate with autoinflammatory disease, and this finding is of wide interest since approximately 60% of human patients with periodic fever syndrome remain genetically unexplained. This investigation also demonstrates how strong artificial selection may affect not only desired and selected phenotypes, but also the health of domestic animals.

Meta-analysis of genome-wide association studies in celiac disease and rheumatoid arthritis identifies fourteen non-HLA shared loci

Epidemiology and candidate gene studies indicate a shared genetic basis for celiac disease (CD) and rheumatoid arthritis (RA), but the extent of this sharing has not been systematically explored. Previous studies demonstrate that 6 of the established non-HLA CD and RA risk loci (out of 26 loci for each disease) are shared between both diseases. We hypothesized that there are additional shared risk alleles and that combining genome-wide association study (GWAS) data from each disease would increase power to identify these shared risk alleles. We performed a meta-analysis of two published GWAS on CD (4,533 cases and 10,750 controls) and RA (5,539 cases and 17,231 controls). After genotyping the top associated SNPs in 2,169 CD cases and 2,255 controls, and 2,845 RA cases and 4,944 controls, 8 additional SNPs demonstrated P<5×10⁻⁸ in a combined analysis of all 50,266 samples, including four SNPs that have not been previously confirmed in either disease: rs10892279 near the DDX6 gene (P_combined = 1.2×10⁻¹²), rs864537 near CD247 (P_combined = 2.2×10⁻¹¹), rs2298428 near UBE2L3 (P_combined = 2.5×10⁻¹⁰), and rs11203203 near UBASH3A (P_combined = 1.1×10⁻⁸). We also confirmed that 4 gene loci previously established in either CD or RA are associated with the other autoimmune disease at combined P<5×10⁻⁸ (SH2B3, 8q24, STAT4, and TRAF1-C5). From the 14 shared gene loci, 7 SNPs showed a genome-wide significant effect on expression of one or more transcripts in the linkage disequilibrium (LD) block around the SNP. These associations implicate antigen presentation and T-cell activation as a shared mechanism of disease pathogenesis and underscore the utility of cross-disease meta-analysis for identification of genetic risk factors with pleiotropic effects between two clinically distinct diseases.

Celiac disease (CD) and rheumatoid arthritis (RA) are two autoimmune diseases characterized by distinct clinical features but increased co-occurrence in families and individuals. Genome-wide association studies (GWAS) performed in CD and RA have identified the HLA region and 26 non-HLA genetic risk loci in each disease. Of the 26 CD and 26 RA risk loci, previous studies have shown that six are shared between the two diseases. In this study we aimed to identify additional shared risk alleles and, in doing so, gain more insight into shared disease pathogenesis. We first empirically investigated the distribution of putative risk alleles from GWAS across both diseases (after removing known risk loci for both diseases). We found that CD risk alleles are non-randomly distributed in the RA GWAS (and vice versa), indicating that CD risk alleles have an increased prior probability of being associated with RA (and vice versa). Next, we performed a GWAS meta-analysis to search for shared risk alleles by combing the RA and CD GWAS, performing both directional and opposite allelic effect analyses, followed by replication testing in independent case-control datasets in both diseases. In addition to the already established six non-HLA shared risk loci, we observed statistically robust associations at eight SNPs, thereby increasing the number of shared non-HLA risk loci to fourteen. Finally, we used gene expression studies and pathway analysis tools to identify the plausible candidate genes in the fourteen associated loci. We observed remarkable overrepresentation of T-cell signaling molecules among the shared genes.

Effective sample size: Quick estimation of the effect of related samples in genetic case-control association analyses

Affected relatives are essential for pedigree linkage analysis, however, they cause a violation of the independent sample assumption in case-control association studies. To avoid the correlation between samples, a common practice is to take only one affected sample per pedigree in association analysis. Although several methods exist in handling correlated samples, they are still not widely used in part because these are not easily implemented, or because they are not widely known. We advocate the effective sample size method as a simple and accessible approach for case-control association analysis with correlated samples. This method modifies the chi-square test statistic, p-value, and 95% confidence interval of the odds-ratio by replacing the apparent number of allele or genotype counts with the effective ones in the standard formula, without the need for specialized computer programs. We present a simple formula for calculating effective sample size for many types of relative pairs and relative sets. For allele frequency estimation, the effective sample size method captures the variance inflation exactly. For genotype frequency, simulations showed that effective sample size provides a satisfactory approximation. A gene which is previously identified as a type 1 diabetes susceptibility locus, the interferon-induced helicase gene (IFIH1), is shown to be significantly associated with rheumatoid arthritis when the effective sample size method is applied. This significant association is not established if only one affected sib per pedigree were used in the association analysis. Relationship between the effective sample size method and other methods - the generalized estimation equation, variance of eigenvalues for correlation matrices, and genomic controls - are discussed.

Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci

To identify new genetic risk factors for rheumatoid arthritis, we conducted a genome-wide association study meta-analysis of 5,539 autoantibody-positive individuals with rheumatoid arthritis (cases) and 20,169 controls of European descent, followed by replication in an independent set of 6,768 rheumatoid arthritis cases and 8,806 controls. Of 34 SNPs selected for replication, 7 new rheumatoid arthritis risk alleles were identified at genome-wide significance (P < 5 x 10(-8)) in an analysis of all 41,282 samples. The associated SNPs are near genes of known immune function, including IL6ST, SPRED2, RBPJ, CCR6, IRF5 and PXK. We also refined associations at two established rheumatoid arthritis risk loci (IL2RA and CCL21) and confirmed the association at AFF3. These new associations bring the total number of confirmed rheumatoid arthritis risk loci to 31 among individuals of European ancestry. An additional 11 SNPs replicated at P < 0.05, many of which are validated autoimmune risk alleles, suggesting that most represent genuine rheumatoid arthritis risk alleles.

A novel mutation in TNFRSF1A associated with overlapping features of tumor necrosis factor receptor-associated periodic syndrome and hyper-IgD syndrome

We describe a 10-year-old child with a novel mutation, c.352A>G/p.Thr118Ala (T89A) in the tumour necrosis factor receptor superfamily 1A (TNFRSF1A) gene. The patient presented with periodic fevers beginning at 2 years of age. He had overlapping clinical and laboratory features of tumour necrosis factor receptor-associated periodic syndrome (TRAPS) and hyper-IgD syndrome (HIDS). This patient expands the clinical and genetic spectrum of TRAPS.

Data for Genetic Analysis Workshop 16 Problem 1, association analysis of rheumatoid arthritis data

For Genetic Analysis Workshop 16 Problem 1, we provided data for genome-wide association analysis of rheumatoid arthritis. Single-nucleotide polymorphism (SNP) genotype data were provided for 868 cases and 1194 controls that had been assayed using an Illumina 550 k platform. In addition, phenotypic data were provided from genotyping DRB1 alleles, which were classified according to the rheumatoid arthritis shared epitope, levels of anti-cyclic citrullinated peptide, and levels of rheumatoid factor IgM. Several questions could be addressed using the data, including analysis of genetic associations using single SNPs or haplotypes, as well as gene-gene and genetic analysis of SNPs for qualitative and quantitative factors.

Clinical features and functional significance of the P369S/R408Q variant in pyrin, the familial Mediterranean fever protein

OBJECTIVES

Familial Mediterranean fever (FMF) is caused by mutations in MEFV, which encodes pyrin. The nature of substitutions P369S and R408Q in exon 3 remains unclear. Exon 3 encoding pyrin's B-box domain is necessary for interactions with proline serine threonine phosphatase interacting protein 1 (PSTPIP1). The aim was to characterise the phenotype of patients with these substitutions and to determine their functional significance.

METHODS

A database of genetic tests undertaken at the US National Institutes of Health was interrogated. Symptoms and signs were classified according to Tel-Hashomer criteria. Coimmunoprecipitation techniques were employed to determine the variants' effects on pyrin/PSTPIP1 interactions.

RESULTS

A total of 40 symptomatic and 4 asymptomatic family members with these substitutions were identified. P369S and R408Q were found in cis, and cosegregated in all patients sequenced. Clinical details were available on 22 patients. In all, 5 patients had symptoms and signs fulfilling a clinical diagnosis of FMF, and 15 received colchicine. In patients not achieving the criteria, trials of anti-tumour necrosis factor (TNF) agents resulted in partial or no benefit; resolution of symptoms was noted in those receiving anakinra. The carrier frequency was higher in the patient cohort than in controls but was not statistically significant. Coimmunoprecipitation studies demonstrated that these pyrin variants did not affect binding to PSTPIP1.

CONCLUSIONS

P369S/R408Q substitutions are associated with a highly variable phenotype, and are infrequently associated with typical FMF symptoms, however a trial of colchicine is warranted in all. Functional and modelling studies suggest that these substitutions do not significantly affect pyrin's interaction with PSTPIP1. This study highlights the need for caution in interpreting genetic tests in patients with atypical symptoms.

Familial Mediterranean fever with a single MEFV mutation: where is the second hit?

OBJECTIVE

Familial Mediterranean fever (FMF) has traditionally been considered an autosomal-recessive disease; however, it has been observed that a substantial number of patients with clinical FMF possess only 1 demonstrable MEFV mutation. The purpose of this study was to perform an extensive search for a second MEFV mutation in 46 patients diagnosed clinically as having FMF and carrying only 1 high-penetrance FMF mutation.

METHODS

MEFV and other candidate genes were sequenced by standard capillary electrophoresis. In 10 patients, the entire 15-kb MEFV genomic region was resequenced using hybridization-based chip technology. MEFV gene expression levels were determined by quantitative reverse transcription-polymerase chain reaction. Pyrin protein levels were examined by Western blotting.

RESULTS

A second MEFV mutation was not identified in any of the patients who were screened. Haplotype analysis did not identify a common haplotype that might be associated with the transmission of a second FMF allele. Western blots did not demonstrate a significant difference in pyrin levels between patients with a single mutation and those with a double mutation; however, FMF patients of both types showed higher protein expression as compared with controls and with non-FMF patients with active inflammation. Screening of genes encoding pyrin-interacting proteins identified rare mutations in a small number of patients, suggesting the possibility of digenic inheritance.

CONCLUSION

Our data underscore the existence of a significant subset of FMF patients who are carriers of only 1 MEFV mutation and demonstrate that complete MEFV sequencing is not likely to yield a second mutation. Screening for the set of the most common mutations and detection of a single mutation appears to be sufficient in the presence of clinical symptoms for the diagnosis of FMF and the initiation of a trial of colchicine.

REL, encoding a member of the NF-kappaB family of transcription factors, is a newly defined risk locus for rheumatoid arthritis

We conducted a genome-wide association study of rheumatoid arthritis in 2,418 cases and 4,504 controls from North America and identified an association at the REL locus, encoding c-Rel, on chromosome 2p13 (rs13031237, P = 6.01 x 10(-10)). Replication in independent case-control datasets comprising 2,604 cases and 2,882 controls confirmed this association, yielding an allelic OR = 1.25 (P = 3.08 x 10(-14)) for marker rs13031237 and an allelic OR = 1.21 (P = 2.60 x 10(-11)) for marker rs13017599 in the combined dataset. The combined dataset also provides definitive support for associations at both CTLA4 (rs231735; OR = 0.85; P = 6.25 x 10(-9)) and BLK (rs2736340; OR = 1.19; P = 5.69 x 10(-9)). c-Rel is an NF-kappaB family member with distinct functional properties in hematopoietic cells, and its association with rheumatoid arthritis suggests disease pathways that involve other recently identified rheumatoid arthritis susceptibility genes including CD40, TRAF1, TNFAIP3 and PRKCQ.

Functional phosphodiesterase 11A mutations may modify the risk of familial and bilateral testicular germ cell tumors

Inactivating germline mutations in phosphodiesterase 11A (PDE11A) have been implicated in adrenal tumor susceptibility. PDE11A is highly expressed in endocrine steroidogenic tissues, especially the testis, and mice with inactivated Pde11a exhibit male infertility, a known testicular germ cell tumor (TGCT) risk factor. We sequenced the PDE11A gene-coding region in 95 patients with TGCT from 64 unrelated kindreds. We identified 8 nonsynonymous substitutions in 20 patients from 15 families: four (R52T, F258Y, G291R, and V820M) were newly recognized, three (R804H, R867G, and M878V) were functional variants previously implicated in adrenal tumor predisposition, and one (Y727C) was a known polymorphism. We compared the frequency of these variants in our patients to unrelated controls that had been screened and found negative for any endocrine diseases: only the two previously reported variants, R804H and R867G, known to be frequent in general population, were detected in these controls. The frequency of all PDE11A-gene variants (combined) was significantly higher among patients with TGCT (P = 0.0002), present in 19% of the families of our cohort. Most variants were detected in the general population, but functional studies showed that all these mutations reduced PDE activity, and that PDE11A protein expression was decreased (or absent) in TGCT samples from carriers. This is the first demonstration of the involvement of a PDE gene in TGCT, although the cyclic AMP signaling pathway has been investigated extensively in reproductive organ function and their diseases. In conclusion, we report that PDE11A-inactivating sequence variants may modify the risk of familial and bilateral TGCT.

An autoinflammatory disease with deficiency of the interleukin-1-receptor antagonist

BACKGROUND

Autoinflammatory diseases manifest inflammation without evidence of infection, high-titer autoantibodies, or autoreactive T cells. We report a disorder caused by mutations of IL1RN, which encodes the interleukin-1-receptor antagonist, with prominent involvement of skin and bone.

METHODS

We studied nine children from six families who had neonatal onset of sterile multifocal osteomyelitis, periostitis, and pustulosis. Response to empirical treatment with the recombinant interleukin-1-receptor antagonist anakinra in the first patient prompted us to test for the presence of mutations and changes in proteins and their function in interleukin-1-pathway genes including IL1RN.

RESULTS

We identified homozygous mutations of IL1RN in nine affected children, from one family from Newfoundland, Canada, three families from The Netherlands, and one consanguineous family from Lebanon. A nonconsanguineous patient from Puerto Rico was homozygous for a genomic deletion that includes IL1RN and five other interleukin-1-family members. At least three of the mutations are founder mutations; heterozygous carriers were asymptomatic, with no cytokine abnormalities in vitro. The IL1RN mutations resulted in a truncated protein that is not secreted, thereby rendering cells hyperresponsive to interleukin-1beta stimulation. Patients treated with anakinra responded rapidly.

CONCLUSIONS

We propose the term deficiency of the interleukin-1-receptor antagonist, or DIRA, to denote this autosomal recessive autoinflammatory disease caused by mutations affecting IL1RN. The absence of interleukin-1-receptor antagonist allows unopposed action of interleukin-1, resulting in life-threatening systemic inflammation with skin and bone involvement. (ClinicalTrials.gov number, NCT00059748.)

Genetic risk factors for rheumatoid arthritis differ in Caucasian and Korean populations

OBJECTIVE

Recent studies have identified a number of novel rheumatoid arthritis (RA) susceptibility loci in Caucasian populations. The aim of this study was to determine whether the genetic variants at 4q27, 6q23, CCL21, TRAF1/C5, and CD40 identified in Caucasians are also associated with RA in a Korean case-control collection. We also comprehensively evaluated the genetic variation within PTPN22, a well-established autoimmune disease-associated gene.

METHODS

We designed an experiment to thoroughly evaluate the PTPN22 linkage disequilibrium region, using tag single-nucleotide polymorphisms (SNPs) and disease-associated SNPs at 5 RA-associated loci recently identified in Caucasians, in 1,128 Korean patients with RA and 1,022 ethnically matched control subjects. We also resequenced the PTPN22 gene to seek novel coding variants that might be contributing to disease in this population.

RESULTS

None of the susceptibility loci identified in Caucasian patients with RA contributed significantly to disease in Koreans. Although tag SNPs covering the PTPN22 linkage disequilibrium block were polymorphic, they did not reveal any disease association, and resequencing did not identify any new common coding region variants in this population. The 6q23 and 4q27 SNPs assayed were nonpolymorphic in this population, and the TRAF1/C5, CD40, and CCL21 SNPs did not show any evidence for association with RA in this population of Korean patients.

CONCLUSION

The genetic risk factors for RA are different in Caucasian and Korean patients. Although patients of different ethnic groups share the HLA region as a major genetic risk locus, most other genes shown to be significantly associated with disease in Caucasians appear not to play a role in Korean patients with RA.

The chromosome 7q region association with rheumatoid arthritis in females in a British population is not replicated in a North American case-control series

OBJECTIVE

The single-nucleotide polymorphism (SNP) rs11761231 on chromosome 7q has been reported to be sexually dimorphic marker for rheumatoid arthritis (RA) susceptibility in a British population. We sought to replicate this finding and to better characterize susceptibility alleles in the region in a North American population.

METHODS

DNA from 2 North American collections of RA patients and controls (1,605 cases and 2,640 controls) was genotyped for rs11761231 and 16 additional chromosome 7q tag SNPs using Sequenom iPlex assays. Association tests were performed for each collection and also separately, contrasting male cases with male controls and female cases with female controls. Principal components analysis (EigenStrat) was used to determine association with RA before and after adjusting for population stratification in the subset of the samples for which there were whole-genome SNP data (772 cases and 1,213 controls).

RESULTS

We failed to replicate an association of the 7q region with RA. Initially, rs11761231 showed evidence for association with RA in the North American Rheumatoid Arthritis Consortium (NARAC) collection (P=0.0073), and rs11765576 showed association with RA in both the NARAC (P=0.038) and RA replication (P = 0.0013) collections. These markers also exhibited sex differentiation. However, in the whole-genome subset, neither SNP showed significant association with RA after correction for population stratification.

CONCLUSION

While 2 SNPs on chromosome 7q appeared to be associated with RA in a North American cohort, the significance of this finding did not withstand correction for population substructure. Our results emphasize the need to carefully account for population structure to avoid false-positive disease associations.