Genome-wide association study identifies kallikrein 5 in type 2 inflammation-low asthma.

BACKGROUND: Clinical studies of type 2 (T2) cytokine-related neutralizing antibodies in asthma have identified a substantial subset of patients with low levels of T2 inflammation who do not benefit from T2 cytokine neutralizing antibody treatment. Non-T2 mechanisms are poorly understood in asthma but represent a redefined unmet medical need. OBJECTIVE: We sought to gain a better understanding of genetic contributions to T2-low asthma. METHODS: We utilized an unbiased genome-wide association study of patients with moderate to severe asthma stratified by T2 serum biomarker periostin. We also performed additional expression and biological analysis for the top genetic hits. RESULTS: We identified a novel protective single nucleotide polymorphism at chr19q13.41, which is selectively associated with T2-low asthma and establishes Kallikrein-related peptidase 5 (KLK5) as the causal gene mediating this association. Heterozygous carriers of the single nucleotide polymorphisms have reduced KLK5 expression. KLK5 is secreted by human bronchial epithelial cells and elevated in asthma bronchial alveolar lavage. T2 cytokines IL-4 and IL-13 downregulate KLK5 in human bronchial epithelial cells. KLK5, dependent on its catalytic function, induces epithelial chemokine/cytokine expression. Finally, overexpression of KLK5 in airway or lack of an endogenous KLK5 inhibitor, SPINK5, leads to spontaneous airway neutrophilic inflammation. CONCLUSION: Our data identify KLK5 to be the causal gene at a novel locus at chr19q13.41 associated with T2-low asthma.

ADAR and hnRNPC deficiency synergize in activating endogenous dsRNA-induced type I IFN responses.

Cytosolic double-stranded RNA (dsRNA) initiates type I IFN responses. Endogenous retroelements, notably Alu elements, constitute a source of dsRNA. Adenosine-to-inosine (A-to-I) editing by ADAR induces mismatches in dsRNA and prevents recognition by MDA5 and autoinflammation. To identify additional endogenous dsRNA checkpoints, we conducted a candidate screen in THP-1 monocytes and found that hnRNPC and ADAR deficiency resulted in synergistic induction of MDA5-dependent IFN responses. RNA-seq analysis demonstrated dysregulation of Alu-containing introns in hnRNPC-deficient cells via utilization of unmasked cryptic splice sites, including introns containing ADAR-dependent A-to-I editing clusters. These putative MDA5 ligands showed reduced editing in the absence of ADAR, providing a plausible mechanism for the combined effects of hnRNPC and ADAR. This study contributes to our understanding of the control of repetitive element-induced autoinflammation and suggests that patients with hnRNPC-mutated tumors might maximally benefit from ADAR inhibition-based immunotherapy.

STXBP6 and B3GNT6 Genes are Associated With Selective IgA Deficiency.

Immunoglobulin A Deficiency (IgAD) is a polygenic primary immune deficiency, with a strong genetic association to the human leukocyte antigen (HLA) region. Previous genome-wide association studies (GWAS) have identified five non-HLA risk loci (IFIH1, PVT1, ATG13-AMBRA1, AHI1 and CLEC16A). In this study, we investigated the genetic interactions between different HLA susceptibility haplotypes and non-MHC genes in IgAD. To do this, we stratified IgAD subjects and healthy controls based on HLA haplotypes (N = 10,993), and then performed GWAS to identify novel genetic regions contributing to IgAD susceptibility. After replicating previously published HLA risk haplotypes, we compared individuals carrying at least one HLA risk allele (HLA-B*08:01-DRB1*03:01-DQB1*02:01 or HLA-DRB1*07:01-DQB1*02:02 or HLA-DRB1*01-DQB1*05:01) with individuals lacking an HLA risk allele. Subsequently, we stratified subjects based on the susceptibility alleles/haplotypes and performed gene-based association analysis using 572,856 SNPs and 24,125 genes. A significant genome-wide association in STXBP6 (rs4097492; p = 7.63 × 10-9) was observed in the cohort carrying at least one MHC risk allele. We also identified a significant gene-based association for B3GNT6 (P Gene = 2.1 × 10-6) in patients not carrying known HLA susceptibility alleles. Our findings indicate that the etiology of IgAD differs depending on the genetic background of HLA susceptibility haplotypes.

Influence of genetic copy number variants of the human GLUT3 glucose transporter gene SLC2A3 on protein expression, glycolysis and rheumatoid arthritis risk: A genetic replication study.

OBJECTIVES: The gene encoding glucose transporter 3 (GLUT3, SLC2A3) is present in the human population at variable copy number. An overt disease phenotype of SLC2A3 copy number variants has not been reported; however, deletion of SLC2A3 has been previously reported to protect carriers from rheumatoid arthritis, implicating GLUT3 as a therapeutic target in rheumatoid arthritis. Here we aim to perform functional analysis of GLUT3 copy number variants in immune cells, and test the reported protective association of the GLUT3 copy number variants for rheumatoid arthritis in a genetic replication study. METHODS: Cells from genotyped healthy controls were analyzed for SLC2A3/GLUT3 expression and glycolysis capacity. We genotyped the SLC2A3 copy number variant in four independent cohorts of rheumatoid arthritis and controls and one cohort of multiple sclerosis and controls. RESULTS: Heterozygous deletion of SLC2A3 correlates directly with expression levels of GLUT3 and influences glycolysis rates in the human immune system. The frequency of the SLC2A3 copy number variant is not different between rheumatoid arthritis, multiple sclerosis and control groups. CONCLUSIONS: Despite a robust SLC2A3 gene copy number dependent phenotype, our study of large groups of rheumatoid arthritis cases and controls provides no evidence for rheumatoid arthritis disease protection in deletion carriers. These data emphasize the importance of well powered replication studies to confirm or refute genetic associations, particularly for relatively rare variants.

Previously reported placebo-response-associated variants do not predict patient outcomes in inflammatory disease Phase III trial placebo arms.

In clinical trials, a placebo response refers to improvement in disease symptoms arising from the psychological effect of receiving a treatment rather than the actual treatment under investigation. Previous research has reported genomic variation associated with the likelihood of observing a placebo response, but these studies have been limited in scope and have not been validated. Here, we analyzed whole-genome sequencing data from 784 patients undergoing placebo treatment in Phase III Asthma or Rheumatoid Arthritis trials to assess the impact of previously reported variation on patient outcomes in the placebo arms and to identify novel variants associated with the placebo response. Contrary to expectations based on previous reports, we did not observe any statistically significant associations between genomic variants and placebo treatment outcome. Our findings suggest that the biological origin of the placebo response is complex and likely to be variable between disease areas.

LACC1 Regulates TNF and IL-17 in Mouse Models of Arthritis and Inflammation.

Both common and rare genetic variants of laccase domain-containing 1 (LACC1, previously C13orf31) are associated with inflammatory bowel disease, leprosy, Behcet disease, and systemic juvenile idiopathic arthritis. However, the functional relevance of these variants is unclear. In this study, we use LACC1-deficient mice to gain insight into the role of LACC1 in regulating inflammation. Following oral administration of Citrobacter rodentium, LACC1 knockout (KO) mice had more severe colon lesions compared with wildtype (WT) controls. Immunization with collagen II, a collagen-induced arthritis (CIA) model, resulted in an accelerated onset of arthritis and significantly worse arthritis and inflammation in LACC1 KO mice. Similar results were obtained in a mannan-induced arthritis model. Serum and local TNF in CIA paws and C. rodentium colons were significantly increased in LACC1 KO mice compared with WT controls. The percentage of IL-17A-producing CD4+ T cells was elevated in LACC1 KO mice undergoing CIA as well as aged mice compared with WT controls. Neutralization of IL-17, but not TNF, prevented enhanced mannan-induced arthritis in LACC1 KO mice. These data provide new mechanistic insight into the function of LACC1 in regulating TNF and IL-17 during inflammatory responses. We hypothesize that these effects contribute to immune-driven pathologies observed in individuals carrying LACC1 variants.

Paired Immunoglobulin-like Type 2 Receptor Alpha G78R variant alters ligand binding and confers protection to Alzheimer's disease.

Paired Immunoglobulin-like Type 2 Receptor Alpha (PILRA) is a cell surface inhibitory receptor that recognizes specific O-glycosylated proteins and is expressed on various innate immune cell types including microglia. We show here that a common missense variant (G78R, rs1859788) of PILRA is the likely causal allele for the confirmed Alzheimer's disease risk locus at 7q21 (rs1476679). The G78R variant alters the interaction of residues essential for sialic acid engagement, resulting in >50% reduced binding for several PILRA ligands including a novel ligand, complement component 4A, and herpes simplex virus 1 (HSV-1) glycoprotein B. PILRA is an entry receptor for HSV-1 via glycoprotein B, and macrophages derived from R78 homozygous donors showed significantly decreased levels of HSV-1 infection at several multiplicities of infection compared to homozygous G78 macrophages. We propose that PILRA G78R protects individuals from Alzheimer's disease risk via reduced inhibitory signaling in microglia and reduced microglial infection during HSV-1 recurrence.

A meta-analysis of genome-wide association studies identifies 17 new Parkinson's disease risk loci.

Common variant genome-wide association studies (GWASs) have, to date, identified >24 risk loci for Parkinson's disease (PD). To discover additional loci, we carried out a GWAS comparing 6,476 PD cases with 302,042 controls, followed by a meta-analysis with a recent study of over 13,000 PD cases and 95,000 controls at 9,830 overlapping variants. We then tested 35 loci (P < 1 × 10-6) in a replication cohort of 5,851 cases and 5,866 controls. We identified 17 novel risk loci (P < 5 × 10-8) in a joint analysis of 26,035 cases and 403,190 controls. We used a neurocentric strategy to assign candidate risk genes to the loci. We identified protein-altering or cis-expression quantitative trait locus (cis-eQTL) variants in linkage disequilibrium with the index variant in 29 of the 41 PD loci. These results indicate a key role for autophagy and lysosomal biology in PD risk, and suggest potential new drug targets for PD.

Identifying and mitigating batch effects in whole genome sequencing data.

BACKGROUND: Large sample sets of whole genome sequencing with deep coverage are being generated, however assembling datasets from different sources inevitably introduces batch effects. These batch effects are not well understood and can be due to changes in the sequencing protocol or bioinformatics tools used to process the data. No systematic algorithms or heuristics exist to detect and filter batch effects or remove associations impacted by batch effects in whole genome sequencing data. RESULTS: We describe key quality metrics, provide a freely available software package to compute them, and demonstrate that identification of batch effects is aided by principal components analysis of these metrics. To mitigate batch effects, we developed new site-specific filters that identified and removed variants that falsely associated with the phenotype due to batch effect. These include filtering based on: a haplotype based genotype correction, a differential genotype quality test, and removing sites with missing genotype rate greater than 30% after setting genotypes with quality scores less than 20 to missing. This method removed 96.1% of unconfirmed genome-wide significant SNP associations and 97.6% of unconfirmed genome-wide significant indel associations. We performed analyses to demonstrate that: 1) These filters impacted variants known to be disease associated as 2 out of 16 confirmed associations in an AMD candidate SNP analysis were filtered, representing a reduction in power of 12.5%, 2) In the absence of batch effects, these filters removed only a small proportion of variants across the genome (type I error rate of 3%), and 3) in an independent dataset, the method removed 90.2% of unconfirmed genome-wide SNP associations and 89.8% of unconfirmed genome-wide indel associations. CONCLUSIONS: Researchers currently do not have effective tools to identify and mitigate batch effects in whole genome sequencing data. We developed and validated methods and filters to address this deficiency.

The Ox40/Ox40 Ligand Pathway Promotes Pathogenic Th Cell Responses, Plasmablast Accumulation, and Lupus Nephritis in NZB/W F1 Mice.

Ox40 ligand (Ox40L) locus genetic variants are associated with the risk for systemic lupus erythematosus (SLE); however, it is unclear how Ox40L contributes to SLE pathogenesis. In this study, we evaluated the contribution of Ox40L and its cognate receptor, Ox40, using in vivo agonist and antagonist approaches in the NZB × NZW (NZB/W) F1 mouse model of SLE. Ox40 was highly expressed on several CD4 Th cell subsets in the spleen and kidney of diseased mice, and expression correlated with disease severity. Treatment of aged NZB/W F1 mice with agonist anti-Ox40 mAbs potently exacerbated renal disease, which was accompanied by activation of kidney-infiltrating T cells and cytokine production. The agonist mAbs also induced activation and inflammatory gene expression in splenic CD4 T cells, including IFN-regulated genes, increased the number of follicular helper T cells and plasmablasts in the spleen, and led to elevated levels of serum IgM and enhanced renal glomerular IgM deposition. In a type I IFN-accelerated lupus model, treatment with an antagonist Ox40:Fc fusion protein significantly delayed the onset of severe proteinuria and improved survival. These data support the hypothesis that the Ox40/Ox40L pathway drives cellular and humoral autoimmune responses during lupus nephritis in NZB/W F1 mice and emphasize the potential clinical value of targeting this pathway in human lupus.

A Common Variant of IL-6R is Associated with Elevated IL-6 Pathway Activity in Alzheimer's Disease Brains.

The common p.D358A variant (rs2228145) in IL-6R is associated with risk for multiple diseases and with increased levels of soluble IL-6R in the periphery and central nervous system (CNS). Here, we show that the p.D358A allele leads to increased proteolysis of membrane bound IL-6R and demonstrate that IL-6R peptides with A358 are more susceptible to cleavage by ADAM10 and ADAM17. IL-6 responsive genes were identified in primary astrocytes and microglia and an IL-6 gene signature was increased in the CNS of late onset Alzheimer's disease subjects in an IL6R allele dependent manner. We conducted a screen to identify variants associated with the age of onset of Alzheimer's disease in APOE ɛ4 carriers. Across five datasets, p.D358A had a meta P = 3 ×10-4 and an odds ratio = 1.3, 95% confidence interval 1.12 -1.48. Our study suggests that a common coding region variant of the IL-6 receptor results in neuroinflammatory changes that may influence the age of onset of Alzheimer's disease in APOE ɛ4 carriers.

geneAttribution: trait agnostic identification of candidate genes associated with noncoding variation.

Motivation: We have developed geneAttribution, an R package that assigns candidate causal gene(s) to a risk variant identified by a genetic association study such as a GWAS. The method combines user-supplied functional annotation such as expression quantitative trait loci (eQTL) or Hi-C genome conformation data and reports the most likely candidate genes. In the absence of annotation data, geneAttribution relies on the distances between the genes and the input variant. Availability and Implementation: The package is freely available from http://www.bioconductor.org/ . A quick-start vignette is included with the package. Contact: wustera@gene.com.

Common variants at PVT1, ATG13-AMBRA1, AHI1 and CLEC16A are associated with selective IgA deficiency.

Selective immunoglobulin A deficiency (IgAD) is the most common primary immunodeficiency in Europeans. Our genome-wide association study (GWAS) meta-analysis of 1,635 patients with IgAD and 4,852 controls identified four new significant (P < 5 × 10-8) loci and association with a rare IFIH1 variant (p.Ile923Val). Peak new variants (PVT1, P = 4.3 × 10-11; ATG13-AMBRA1, P = 6.7 × 10-10; AHI1, P = 8.4 × 10-10; CLEC16A, P = 1.4 × 10-9) overlapped with autoimmune markers (3/4) and correlated with 21 putative regulatory variants, including expression quantitative trait loci (eQTLs) for AHI1 and DEXI and DNase hypersensitivity sites in FOXP3+ regulatory T cells. Pathway analysis of the meta-analysis results showed striking association with the KEGG pathway for IgA production (pathway P < 0.0001), with 22 of the 30 annotated pathway genes containing at least one variant with P ≤ 0.05 in the IgAD meta-analysis. These data suggest that a complex network of genetic effects, including genes known to influence the biology of IgA production, contributes to IgAD.

Thinking differently about lupus.

A search for the genetic causes of an autoimmune disease called systemic lupus erythematosus reveals a new twist on an old story.

Genetic association analyses implicate aberrant regulation of innate and adaptive immunity genes in the pathogenesis of systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a genetically complex autoimmune disease characterized by loss of immune tolerance to nuclear and cell surface antigens. Previous genome-wide association studies (GWAS) had modest sample sizes, reducing their scope and reliability. Our study comprised 7,219 cases and 15,991 controls of European ancestry, constituting a new GWAS, a meta-analysis with a published GWAS and a replication study. We have mapped 43 susceptibility loci, including ten new associations. Assisted by dense genome coverage, imputation provided evidence for missense variants underpinning associations in eight genes. Other likely causal genes were established by examining associated alleles for cis-acting eQTL effects in a range of ex vivo immune cells. We found an over-representation (n = 16) of transcription factors among SLE susceptibility genes. This finding supports the view that aberrantly regulated gene expression networks in multiple cell types in both the innate and adaptive immune response contribute to the risk of developing SLE.

Genetic Modifiers of Age at Onset in Carriers of the G206A Mutation in PSEN1 With Familial Alzheimer Disease Among Caribbean Hispanics.

IMPORTANCE: The present study identified potential genetic modifiers that may delay or accelerate age at onset of familial Alzheimer disease (AD) by examining age at onset in PSEN1 mutation carrier families, and further investigation of these modifiers may provide insight into the pathobiology of AD and potential therapeutic measures. OBJECTIVE: To identify genetic variants that modify age at onset of AD. DESIGN, SETTING, AND PARTICIPANTS: Using a subset of Caribbean Hispanic families that carry the PSEN1 p.G206A mutation, we performed a 2-stage genome study. The mutation carrier families from an ongoing genetic study served as a discovery set, and the cohort of those with LOAD served as a confirmation set. To identify candidate loci, we performed linkage analysis using 5 p.G206A carrier families (n = 56), and we also performed whole-exome association analysis using 31 p.G206A carriers from 26 families. To confirm the genetic modifiers identified from the p.G206A carrier families, we analyzed the GWAS data for 2888 elderly individuals with LOAD. All study participants were Caribbean Hispanics. MAIN OUTCOMES AND MEASURES: Age at onset of AD. RESULTS: Linkage analysis of AD identified the strongest linkage support at 4q35 (LOD [logarithm of odds] score, 3.69), and the GWAS of age at onset identified variants on 1p13.1, 2q13, 4q25, and 17p11. In the confirmation stage, genewise analysis identified SNX25, PDLIM3, and 3 SH3 domain genes (SORBS2, SH3RF3, and NPHP1) to be significantly associated with LOAD. Subsequent allelic association analysis confirmed SNX25, PDLIM3, and SORBS2 as genetic modifiers of age at onset of EOAD and LOAD and provided modest support for SH3RF3 and NPHP1. CONCLUSIONS AND RELEVANCE: Our 2-stage analysis revealed that SNX25, PDLIM3, and SORBS2 may serve as genetic modifiers of age at onset in both EOAD and LOAD.

Association of the interferon signature metric with serological disease manifestations but not global activity scores in multiple cohorts of patients with SLE.

OBJECTIVES: The interferon (IFN) signature (IS) in patients with systemic lupus erythematosus (SLE) includes over 100 genes induced by type I IFN pathway activation. We developed a method to quantify the IS using three genes-the IS metric (ISM)-and characterised the clinical characteristics of patients with SLE with different ISM status from multiple clinical trials. METHODS: Blood microarray expression data from a training cohort of patients with SLE confirmed the presence of the IS and identified surrogate genes. We assayed these genes in a quantitative PCR (qPCR) assay, yielding an ISM from the IS. The association of ISM status with clinical disease characteristics was assessed in patients with extrarenal lupus and lupus nephritis from four clinical trials. RESULTS: Three genes, HERC5, EPSTI and CMPK2, correlated well with the IS (p>0.96), and composed the ISM qPCR assay. Using the 95th centile for healthy control data, patients with SLE from different studies were classified into two ISM subsets-ISM-Low and ISM-High-that are longitudinally stable over 36 weeks. Significant associations were identified between ISM-High status and higher titres of anti-dsDNA antibodies, presence of anti extractable nuclear antigen autoantibodies, elevated serum B cell activating factor of the tumour necrosis factor family (BAFF) levels, and hypocomplementaemia. However, measures of overall clinical disease activity were similar for ISM-High and ISM-Low groups. CONCLUSIONS: The ISM is an IS biomarker that divides patients with SLE into two subpopulations-ISM-High and ISM-Low-with differing serological manifestations. The ISM does not distinguish between high and low disease activity, but may have utility in identifying patients more likely to respond to treatment(s) targeting IFN-α. CLINICALTRIALSGOV REGISTRATION NUMBER: NCT00962832.

Rarity of the Alzheimer disease-protective APP A673T variant in the United States.

IMPORTANCE: Recently, a rare variant in the amyloid precursor protein gene (APP) was described in a population from Iceland. This variant, in which alanine is replaced by threonine at position 673 (A673T), appears to protect against late-onset Alzheimer disease (AD). We evaluated the frequency of this variant in AD cases and cognitively normal controls to determine whether this variant will significantly contribute to risk assessment in individuals in the United States. OBJECTIVE: To determine the frequency of the APP A673T variant in a large group of elderly cognitively normal controls and AD cases from the United States and in 2 case-control cohorts from Sweden. DESIGN, SETTING, AND PARTICIPANTS: Case-control association analysis of variant APP A673T in US and Swedish white individuals comparing AD cases with cognitively intact elderly controls. Participants were ascertained at multiple university-associated medical centers and clinics across the United States and Sweden by study-specific sampling methods. They were from case-control studies, community-based prospective cohort studies, and studies that ascertained multiplex families from multiple sources. MAIN OUTCOMES AND MEASURES: Genotypes for the APP A673T variant were determined using the Infinium HumanExome V1 Beadchip (Illumina, Inc) and by TaqMan genotyping (Life Technologies). RESULTS: The A673T variant genotypes were evaluated in 8943 US AD cases, 10 480 US cognitively normal controls, 862 Swedish AD cases, and 707 Swedish cognitively normal controls. We identified 3 US individuals heterozygous for A673T, including 1 AD case (age at onset, 89 years) and 2 controls (age at last examination, 82 and 77 years). The remaining US samples were homozygous for the alanine (A673) allele. In the Swedish samples, 3 controls were heterozygous for A673T and all AD cases were homozygous for the A673 allele. We also genotyped a US family previously reported to harbor the A673T variant and found a mother-daughter pair, both cognitively normal at ages 72 and 84 years, respectively, who were both heterozygous for A673T; however, all individuals with AD in the family were homozygous for A673. CONCLUSIONS AND RELEVANCE: The A673T variant is extremely rare in US cohorts and does not play a substantial role in risk for AD in this population. This variant may be primarily restricted to Icelandic and Scandinavian populations.

A rare mutation in UNC5C predisposes to late-onset Alzheimer's disease and increases neuronal cell death.

We have identified a rare coding mutation, T835M (rs137875858), in the UNC5C netrin receptor gene that segregated with disease in an autosomal dominant pattern in two families enriched for late-onset Alzheimer's disease and that was associated with disease across four large case-control cohorts (odds ratio = 2.15, Pmeta = 0.0095). T835M alters a conserved residue in the hinge region of UNC5C, and in vitro studies demonstrate that this mutation leads to increased cell death in human HEK293T cells and in rodent neurons. Furthermore, neurons expressing T835M UNC5C are more susceptible to cell death from multiple neurotoxic stimuli, including ß-amyloid (Aß), glutamate and staurosporine. On the basis of these data and the enriched hippocampal expression of UNC5C in the adult nervous system, we propose that one possible mechanism in which T835M UNC5C contributes to the risk of Alzheimer's disease is by increasing susceptibility to neuronal cell death, particularly in vulnerable regions of the Alzheimer's disease brain.