A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus.

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

Screening for Functional Non-coding Genetic Variants Using Electrophoretic Mobility Shift Assay (EMSA) and DNA-affinity Precipitation Assay (DAPA).

Population and family-based genetic studies typically result in the identification of genetic variants that are statistically associated with a clinical disease or phenotype. For many diseases and traits, most variants are non-coding, and are thus likely to act by impacting subtle, comparatively hard to predict mechanisms controlling gene expression. Here, we describe a general strategic approach to prioritize non-coding variants, and screen them for their function. This approach involves computational prioritization using functional genomic databases followed by experimental analysis of differential binding of transcription factors (TFs) to risk and non-risk alleles. For both electrophoretic mobility shift assay (EMSA) and DNA affinity precipitation assay (DAPA) analysis of genetic variants, a synthetic DNA oligonucleotide (oligo) is used to identify factors in the nuclear lysate of disease or phenotype-relevant cells. For EMSA, the oligonucleotides with or without bound nuclear factors (often TFs) are analyzed by non-denaturing electrophoresis on a tris-borate-EDTA (TBE) polyacrylamide gel. For DAPA, the oligonucleotides are bound to a magnetic column and the nuclear factors that specifically bind the DNA sequence are eluted and analyzed through mass spectrometry or with a reducing sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) followed by Western blot analysis. This general approach can be widely used to study the function of non-coding genetic variants associated with any disease, trait, or phenotype.

Evaluation of the contribution of multiple DAMPs and DAMP receptors in cell death-induced sterile inflammatory responses.

When cells die by necrosis in vivo they stimulate an inflammatory response. It is thought that this response is triggered when the injured cells expose proinflammatory molecules, collectively referred to as damage associated molecular patterns (DAMPs), which are recognized by cells or soluble molecules of the innate or adaptive immune system. Several putative DAMPs and/or their receptors have been identified, but whether and how much they participate in responses in vivo is incompletely understood, and they have not previously been compared side-by-side in the same models. This study focuses on evaluating the contribution of multiple mechanisms that have been proposed to or potentially could participate in cell death-induced inflammation: The third component of complement (C3), ATP (and its receptor P2X7), antibodies, the C-type lectin receptor Mincle (Clec4e), and protease-activated receptor 2 (PAR2). We investigate the role of these factors in cell death-induced inflammation to dead cells in the peritoneum and acetaminophen-induced liver damage. We find that mice deficient in antibody, C3 or PAR2 have impaired inflammatory responses to dying cells. In contrast there was no reduction in inflammation to cell death in the peritoneum or liver of mice that genetically lack Mincle, the P2X7 receptor or that were treated with apyrase to deplete ATP. These results indicate that antibody, complement and PAR2 contribute to cell death-induced inflammation but that Mincle and ATP- P2X7 receptor are not required for this response in at least 2 different in vivo models.

The struggle to find reliable results in exome sequencing data: filtering out Mendelian errors.

Next Generation Sequencing studies generate a large quantity of genetic data in a relatively cost and time efficient manner and provide an unprecedented opportunity to identify candidate causative variants that lead to disease phenotypes. A challenge to these studies is the generation of sequencing artifacts by current technologies. To identify and characterize the properties that distinguish false positive variants from true variants, we sequenced a child and both parents (one trio) using DNA isolated from three sources (blood, buccal cells, and saliva). The trio strategy allowed us to identify variants in the proband that could not have been inherited from the parents (Mendelian errors) and would most likely indicate sequencing artifacts. Quality control measurements were examined and three measurements were found to identify the greatest number of Mendelian errors. These included read depth, genotype quality score, and alternate allele ratio. Filtering the variants on these measurements removed ~95% of the Mendelian errors while retaining 80% of the called variants. These filters were applied independently. After filtering, the concordance between identical samples isolated from different sources was 99.99% as compared to 87% before filtering. This high concordance suggests that different sources of DNA can be used in trio studies without affecting the ability to identify causative polymorphisms. To facilitate analysis of next generation sequencing data, we developed the Cincinnati Analytical Suite for Sequencing Informatics (CASSI) to store sequencing files, metadata (eg. relatedness information), file versioning, data filtering, variant annotation, and identify candidate causative polymorphisms that follow either de novo, rare recessive homozygous or compound heterozygous inheritance models. We conclude the data cleaning process improves the signal to noise ratio in terms of variants and facilitates the identification of candidate disease causative polymorphisms.

Lupus risk variants in the PXK locus alter B-cell receptor internalization.

Genome wide association studies have identified variants in PXK that confer risk for humoral autoimmune diseases, including systemic lupus erythematosus (SLE or lupus), rheumatoid arthritis and more recently systemic sclerosis. While PXK is involved in trafficking of epidermal growth factor Receptor (EGFR) in COS-7 cells, mechanisms linking PXK to lupus pathophysiology have remained undefined. In an effort to uncover the mechanism at this locus that increases lupus-risk, we undertook a fine-mapping analysis in a large multi-ancestral study of lupus patients and controls. We define a large (257kb) common haplotype marking a single causal variant that confers lupus risk detected only in European ancestral populations and spans the promoter through the 3' UTR of PXK. The strongest association was found at rs6445972 with P < 4.62 × 10(-10), OR 0.81 (0.75-0.86). Using stepwise logistic regression analysis, we demonstrate that one signal drives the genetic association in the region. Bayesian analysis confirms our results, identifying a 95% credible set consisting of 172 variants spanning 202 kb. Functionally, we found that PXK operates on the B-cell antigen receptor (BCR); we confirmed that PXK influenced the rate of BCR internalization. Furthermore, we demonstrate that individuals carrying the risk haplotype exhibited a decreased rate of BCR internalization, a process known to impact B cell survival and cell fate. Taken together, these data define a new candidate mechanism for the genetic association of variants around PXK with lupus risk and highlight the regulation of intracellular trafficking as a genetically regulated pathway mediating human autoimmunity.

The IL-1-dependent sterile inflammatory response has a substantial caspase-1-independent component that requires cathepsin C.

The sterile inflammatory response to cell death and irritant crystals is medically important because it causes disease. Although these stimuli are structurally distinct, they cause inflammation through a common pathway that requires the cytokine IL-1. In vitro, the inflammasome, and in particular its generation of active caspase-1, is absolutely required to produce bioactive IL-1ß. However, in this study, we report that caspase-1 is not required in vivo for much of the IL-1ß-dependent sterile inflammatory response. Furthermore, we find that cathepsin C, which controls the activity of a number of leukocyte serine proteases capable of processing IL-1ß, plays a major role in this caspase-1-independent pathway. Mice that are deficient in cathepsin C have reduced inflammatory responses to dying cells and silica crystals. In the absence of cathepsin C, caspase-1 becomes rate limiting such that mice doubly deficient in both of these proteases make little IL-1ß in vivo and have markedly attenuated inflammatory responses to the sterile stimuli. In contrast, these mutant mice generate normal inflammation in response to exogenous IL-1ß, indicating that cathepsin C and caspase-1 function upstream of IL-1ß, and, in their absence, all components of the pathway downstream of mature IL-1ß are intact.