Personalizing medicine for autoimmune and inflammatory diseases.

Therapies that target molecular pathways do not provide uniform benefits for all patients at present. New transformative therapies for autoimmune and inflammatory diseases require greater molecular understanding of patient subsets and the ability to personalize targeted therapies for each subset.

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.

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.

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.

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.

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.

A mutation in APP protects against Alzheimer's disease and age-related cognitive decline.

The prevalence of dementia in the Western world in people over the age of 60 has been estimated to be greater than 5%, about two-thirds of which are due to Alzheimer's disease. The age-specific prevalence of Alzheimer's disease nearly doubles every 5 years after age 65, leading to a prevalence of greater than 25% in those over the age of 90 (ref. 3). Here, to search for low-frequency variants in the amyloid-ß precursor protein (APP) gene with a significant effect on the risk of Alzheimer's disease, we studied coding variants in APP in a set of whole-genome sequence data from 1,795 Icelanders. We found a coding mutation (A673T) in the APP gene that protects against Alzheimer's disease and cognitive decline in the elderly without Alzheimer's disease. This substitution is adjacent to the aspartyl protease ß-site in APP, and results in an approximately 40% reduction in the formation of amyloidogenic peptides in vitro. The strong protective effect of the A673T substitution against Alzheimer's disease provides proof of principle for the hypothesis that reducing the ß-cleavage of APP may protect against the disease. Furthermore, as the A673T allele also protects against cognitive decline in the elderly without Alzheimer's disease, the two may be mediated through the same or similar mechanisms.

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.

Two functional lupus-associated BLK promoter variants control cell-type- and developmental-stage-specific transcription.

Efforts to identify lupus-associated causal variants in the FAM167A/BLK locus on 8p21 are hampered by highly associated noncausal variants. In this report, we used a trans-population mapping and sequencing strategy to identify a common variant (rs922483) in the proximal BLK promoter and a tri-allelic variant (rs1382568) in the upstream alternative BLK promoter as putative causal variants for association with systemic lupus erythematosus. The risk allele (T) at rs922483 reduced proximal promoter activity and modulated alternative promoter usage. Allelic differences at rs1382568 resulted in altered promoter activity in B progenitor cell lines. Thus, our results demonstrated that both lupus-associated functional variants contribute to the autoimmune disease association by modulating transcription of BLK in B cells and thus potentially altering immune responses.