Dysregulated long non-coding RNA in Sjögren's disease impacts both interferon and adaptive immune responses.

OBJECTIVE: Sjögren's disease (SjD) is an autoimmune disease characterised by inflammatory destruction of exocrine glands. Patients with autoantibodies to Ro/SSA (SjDRo+) exhibit more severe disease. Long non-coding RNAs (lncRNAs) are a functionally diverse class of non-protein-coding RNAs whose role in autoimmune disease pathology has not been well characterised. METHODS: Whole blood RNA-sequencing (RNA-seq) was performed on SjD cases (n=23 Ro/SSA negative (SjDRo-); n=27 Ro/SSA positive (SjDRo+) and healthy controls (HCs; n=27). Bioinformatics and pathway analyses of differentially expressed (DE) transcripts (log2 fold change ≥2 or ≤0.5; padj<0.05) were used to predict lncRNA function. LINC01871 was characterised by RNA-seq analyses of HSB-2 cells with CRISPR-targeted LINC01871 deletion (LINC01871-/ -) and in vitro stimulation assays. RESULTS: Whole blood RNA-seq revealed autoantibody-specific transcription profiles and disproportionate downregulation of DE transcripts in SjD cases relative to HCs. Sixteen DE lncRNAs exhibited correlated expression with the interferon (IFN)-regulated gene, RSAD2, in SjDRo+ (r≥0.65 or ≤-0.6); four antisense lncRNAs exhibited IFN-regulated expression in immune cell lines. LINC01871 was upregulated in all SjD cases. RNA-seq and pathway analyses of LINC01871-/ - cells implicated roles in cytotoxic function, differentiation and IFNγ induction. LINC01871 was induced by IFNγ in a myeloid cell line and regulated by calcineurin/NFAT pathway and T cell receptor (TCR) signalling in primary human T cells. CONCLUSION: LINC01871 influences expression of many immune cell genes and growth factors, is IFNγ inducible, and regulated by calcineurin signalling and TCR ligand engagement. Altered LINC01871 expression may influence the dysregulated T cell inflammatory pathways implicated in SjD.

Variants on the UBE2L3/YDJC Autoimmune Disease Risk Haplotype Increase UBE2L3 Expression by Modulating CCCTC-Binding Factor and YY1 Binding.

OBJECTIVE: Genetic variants spanning UBE2L3 are associated with increased expression of the UBE2L3-encoded E2 ubiquitin-conjugating enzyme H7 (UbcH7), which facilitates activation of proinflammatory NF-κB signaling and susceptibility to autoimmune diseases. We undertook this study to delineate how genetic variants carried on the UBE2L3/YDJC autoimmune risk haplotype function to drive hypermorphic UBE2L3 expression. METHODS: We used bioinformatic analyses, electrophoretic mobility shift assays, and luciferase reporter assays to identify and functionally characterize allele-specific effects of risk variants positioned in chromatin accessible regions of immune cells. Chromatin conformation capture with quantitative polymerase chain reaction (3C-qPCR), chromatin immunoprecipitation (ChIP)-qPCR, and small interfering RNA (siRNA) knockdown assays were performed on patient-derived Epstein-Barr virus-transformed B cells homozygous for the UBE2L3/YDJC nonrisk or risk haplotype to determine if the risk haplotype increases UBE2L3 expression by altering the regulatory chromatin architecture in the region. RESULTS: Of the 7 prioritized variants, 5 demonstrated allele-specific increases in nuclear protein binding affinity and regulatory activity. High-throughput sequencing of chromosome conformation capture coupled with ChIP (HiChIP) and 3C-qPCR uncovered a long-range interaction between the UBE2L3 promoter (rs140490, rs140491, rs11089620) and the downstream YDJC promoter (rs3747093) that was strengthened in the presence of the UBE2L3/YDJC risk haplotype, and correlated with the loss of CCCTC-binding factor (CTCF) and gain of YY1 binding at the risk alleles. Depleting YY1 by siRNA disrupted the long-range interaction between the 2 promoters and reduced UBE2L3 expression. CONCLUSION: The UBE2L3/YDJC autoimmune risk haplotype increases UBE2L3 expression through strengthening a YY1-mediated interaction between the UBE2L3 and YDJC promoters.

Role of Systemic Lupus Erythematosus Risk Variants With Opposing Functional Effects as a Driver of Hypomorphic Expression of TNIP1 and Other Genes Within a Three-Dimensional Chromatin Network.

OBJECTIVE: Genetic variants in the region of tumor necrosis factor-induced protein 3-interacting protein 1 (TNIP1) are associated with autoimmune disease and reduced TNIP1 gene expression. The aim of this study was to define the functional genetic mechanisms driving TNIP1 hypomorphic expression imparted by the systemic lupus erythematosus-associated TNIP1 H1 risk haplotype. METHODS: Dual luciferase expression and electrophoretic mobility shift assays were used to evaluate the allelic effects of 11 risk variants on enhancer function and nuclear protein binding in immune cell line models (Epstein-Barr virus [EBV]-transformed human B cells, Jurkat cells, and THP-1 cells), left in a resting state or stimulated with phorbol 12-myristate 13-acetate/ionomycin. HiChIP was used to define the regulatory 3-dimensional (3-D) chromatin network of the TNIP1 haplotype by detecting in situ long-range DNA contacts associated with H3K27ac-marked chromatin in EBV B cells. Then, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used to determine the expression of genes within the 3-D chromatin network. RESULTS: Bioinformatics analyses of 50 single-nucleotide polymorphisms on the TNIP1 H1 risk haplotype identified 11 non-protein-coding variants with a high likelihood of influencing TNIP1 gene expression. Eight variants in EBV B cells, 5 in THP-1 cells, and 2 in Jurkat cells exhibited various allelic effects on enhancer activation, resulting in a cumulative suppressive effect on TNIP1 expression (net effect of risk variants -7.14 fold, -6.80 fold, and -2.44 fold, respectively; n > 3). Specifically, in EBV B cells, only 2 variants (rs10057690 and rs13180950) exhibited allele-specific loss of both enhancer activity and nuclear protein binding (each P < 0.01 relative to nonrisk alleles). In contrast, the rs10036748 risk allele reduced binding affinities of the transcriptional repressors basic helix-loop-helix family member 40/differentially expressed in chondrocytes 1 (bHLHe40/DEC1) (P < 0.05 relative to nonrisk alleles) and CREB-1 (P not significant) in EBV B cells, resulting in a gain of enhancer activity (P < 0.05). HiChIP and qRT-PCR analyses revealed that overall transcriptional repression of the TNIP1 haplotype extended to the neighboring genes DCTN4 and GMA2, both of which also showed decreased expression in the presence of the TNIP1 risk haplotype (P < 0.001 and P < 0.01, respectively, relative to the nonrisk haplotype); notably, it was found that these genes share a 3-D chromatin network. CONCLUSION: Hypomorphic TNIP1 expression results from the combined concordant and opposing effects of multiple risk variants carried on the TNIP1 risk haplotype, with the strongest regulatory effect in B lymphoid lineage cells. Furthermore, the TNIP1 risk haplotype effect extends to neighboring genes within a shared chromatin network.

Enhancer histone-QTLs are enriched on autoimmune risk haplotypes and influence gene expression within chromatin networks.

Genetic variants can confer risk to complex genetic diseases by modulating gene expression through changes to the epigenome. To assess the degree to which genetic variants influence epigenome activity, we integrate epigenetic and genotypic data from lupus patient lymphoblastoid cell lines to identify variants that induce allelic imbalance in the magnitude of histone post-translational modifications, referred to herein as histone quantitative trait loci (hQTLs). We demonstrate that enhancer hQTLs are enriched on autoimmune disease risk haplotypes and disproportionately influence gene expression variability compared with non-hQTL variants in strong linkage disequilibrium. We show that the epigenome regulates HLA class II genes differently in individuals who carry HLA-DR3 or HLA-DR15 haplotypes, resulting in differential 3D chromatin conformation and gene expression. Finally, we identify significant expression QTL (eQTL) x hQTL interactions that reveal substructure within eQTL gene expression, suggesting potential implications for functional genomic studies that leverage eQTL data for subject selection and stratification.

Single-cell analysis of glandular T cell receptors in Sjögren's syndrome.

CD4+ T cells predominate in salivary gland (SG) inflammatory lesions in Sjögren's syndrome (SS). However, their antigen specificity, degree of clonal expansion, and relationship to clinical disease features remain unknown. We used multiplex reverse-transcriptase PCR to amplify paired T cell receptor α (TCRα) and ß transcripts of single CD4+CD45RA- T cells from SG and peripheral blood (PB) of 10 individuals with primary SS, 9 of whom shared the HLA DR3/DQ2 risk haplotype. TCRα and ß sequences were obtained from a median of 91 SG and 107 PB cells per subject. The degree of clonal expansion and frequency of cells expressing two productively rearranged α genes were increased in SG versus PB. Expanded clones from SG exhibited complementary-determining region 3 (CDR3) sequence similarity both within and among subjects, suggesting antigenic selection and shared antigen recognition. CDR3 similarities were shared among expanded clones from individuals discordant for canonical Ro and La autoantibodies, suggesting recognition of alternative SG antigen(s). The extent of SG clonal expansion correlated with reduced saliva production and increased SG fibrosis, linking expanded SG T cells with glandular dysfunction. Knowledge of paired TCRα and ß sequences enables further work toward identification of target antigens and development of novel therapies.

Regulatory polymorphisms modulate the expression of HLA class II molecules and promote autoimmunity.

Targeted sequencing of sixteen SLE risk loci among 1349 Caucasian cases and controls produced a comprehensive dataset of the variations causing susceptibility to systemic lupus erythematosus (SLE). Two independent disease association signals in the HLA-D region identified two regulatory regions containing 3562 polymorphisms that modified thirty-seven transcription factor binding sites. These extensive functional variations are a new and potent facet of HLA polymorphism. Variations modifying the consensus binding motifs of IRF4 and CTCF in the XL9 regulatory complex modified the transcription of HLA-DRB1, HLA-DQA1 and HLA-DQB1 in a chromosome-specific manner, resulting in a 2.5-fold increase in the surface expression of HLA-DR and DQ molecules on dendritic cells with SLE risk genotypes, which increases to over 4-fold after stimulation. Similar analyses of fifteen other SLE risk loci identified 1206 functional variants tightly linked with disease-associated SNPs and demonstrated that common disease alleles contain multiple causal variants modulating multiple immune system genes.

Forward genetic screening identifies a small molecule that blocks Toxoplasma gondii growth by inhibiting both host- and parasite-encoded kinases.

The simultaneous targeting of host and pathogen processes represents an untapped approach for the treatment of intracellular infections. Hypoxia-inducible factor-1 (HIF-1) is a host cell transcription factor that is activated by and required for the growth of the intracellular protozoan parasite Toxoplasma gondii at physiological oxygen levels. Parasite activation of HIF-1 is blocked by inhibiting the family of closely related Activin-Like Kinase (ALK) host cell receptors ALK4, ALK5, and ALK7, which was determined in part by use of an ALK4,5,7 inhibitor named SB505124. Besides inhibiting HIF-1 activation, SB505124 also potently blocks parasite replication under normoxic conditions. To determine whether SB505124 inhibition of parasite growth was exclusively due to inhibition of ALK4,5,7 or because the drug inhibited a second kinase, SB505124-resistant parasites were isolated by chemical mutagenesis. Whole-genome sequencing of these mutants revealed mutations in the Toxoplasma MAP kinase, TgMAPK1. Allelic replacement of mutant TgMAPK1 alleles into wild-type parasites was sufficient to confer SB505124 resistance. SB505124 independently impacts TgMAPK1 and ALK4,5,7 signaling since drug resistant parasites could not activate HIF-1 in the presence of SB505124 or grow in HIF-1 deficient cells. In addition, TgMAPK1 kinase activity is inhibited by SB505124. Finally, mice treated with SB505124 had significantly lower tissue burdens following Toxoplasma infection. These data therefore identify SB505124 as a novel small molecule inhibitor that acts by inhibiting two distinct targets, host HIF-1 and TgMAPK1.

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.

Allelic heterogeneity in NCF2 associated with systemic lupus erythematosus (SLE) susceptibility across four ethnic populations.

Recent reports have associated NCF2, encoding a core component of the multi-protein NADPH oxidase (NADPHO), with systemic lupus erythematosus (SLE) susceptibility in individuals of European ancestry. To identify ethnicity-specific and -robust variants within NCF2, we assessed 145 SNPs in and around the NCF2 gene in 5325 cases and 21 866 controls of European-American (EA), African-American (AA), Hispanic (HS) and Korean (KR) ancestry. Subsequent imputation, conditional, haplotype and bioinformatic analyses identified seven potentially functional SLE-predisposing variants. Association with non-synonymous rs17849502, previously reported in EA, was detected in EA, HS and AA (P(EA) = 1.01 × 10(-54), PHS = 3.68 × 10(-10), P(AA) = 0.03); synonymous rs17849501 was similarly significant. These SNPs were monomorphic in KR. Novel associations were detected with coding variants at rs35937854 in AA (PAA = 1.49 × 10(-9)), and rs13306575 in HS and KR (P(HS) = 7.04 × 10(-7), P(KR) = 3.30 × 10(-3)). In KR, a 3-SNP haplotype was significantly associated (P = 4.20 × 10(-7)), implying that SLE predisposing variants were tagged. Significant SNP-SNP interaction (P = 0.02) was detected between rs13306575 and rs17849502 in HS, and a dramatically increased risk (OR = 6.55) with a risk allele at each locus. Molecular modeling predicts that these non-synonymous mutations could disrupt NADPHO complex assembly. The risk allele of rs17849501, located in a conserved transcriptional regulatory region, increased reporter gene activity, suggesting in vivo enhancer function. Our results not only establish allelic heterogeneity within NCF2 associated with SLE, but also emphasize the utility of multi-ethnic cohorts to identify predisposing variants explaining additional phenotypic variance ('missing heritability') of complex diseases like SLE.

Use of next-generation DNA sequencing to analyze genetic variants in rheumatic disease.

Next-generation DNA sequencing has revolutionized the field of genetics and genomics, providing researchers with the tools to efficiently identify novel rare and low frequency risk variants, which was not practical with previously available methodologies. These methods allow for the sequence capture of a specific locus or small genetic region all the way up to the entire six billion base pairs of the diploid human genome. Rheumatic diseases are a huge burden on the US population, affecting more than 46 million Americans. Those afflicted suffer from one or more of the more than 100 diseases characterized by inflammation and loss of function, mainly of the joints, tendons, ligaments, bones, and muscles. While genetics studies of many of these diseases (for example, systemic lupus erythematosus, rheumatoid arthritis, and inflammatory bowel disease) have had major successes in defining their genetic architecture, causal alleles and rare variants have still been elusive. This review describes the current high-throughput DNA sequencing methodologies commercially available and their application to rheumatic diseases in both case­control as well as family-based studies.

Infinium assay for large-scale SNP genotyping applications.

Genotyping variants in the human genome has proven to be an efficient method to identify genetic associations with phenotypes. The distribution of variants within families or populations can facilitate identification of the genetic factors of disease. Illumina's panel of genotyping BeadChips allows investigators to genotype thousands or millions of single nucleotide polymorphisms (SNPs) or to analyze other genomic variants, such as copy number, across a large number of DNA samples. These SNPs can be spread throughout the genome or targeted in specific regions in order to maximize potential discovery. The Infinium assay has been optimized to yield high-quality, accurate results quickly. With proper setup, a single technician can process from a few hundred to over a thousand DNA samples per week, depending on the type of array. This assay guides users through every step, starting with genomic DNA and ending with the scanning of the array. Using propriety reagents, samples are amplified, fragmented, precipitated, resuspended, hybridized to the chip, extended by a single base, stained, and scanned on either an iScan or Hi Scan high-resolution optical imaging system. One overnight step is required to amplify the DNA. The DNA is denatured and isothermally amplified by whole-genome amplification; therefore, no PCR is required. Samples are hybridized to the arrays during a second overnight step. By the third day, the samples are ready to be scanned and analyzed. Amplified DNA may be stockpiled in large quantities, allowing bead arrays to be processed every day of the week, thereby maximizing throughput.

BAC-pool sequencing and analysis of large segments of A12 and D12 homoeologous chromosomes in upland cotton.

Although new and emerging next-generation sequencing (NGS) technologies have reduced sequencing costs significantly, much work remains to implement them for de novo sequencing of complex and highly repetitive genomes such as the tetraploid genome of Upland cotton (Gossypium hirsutum L.). Herein we report the results from implementing a novel, hybrid Sanger/454-based BAC-pool sequencing strategy using minimum tiling path (MTP) BACs from Ctg-3301 and Ctg-465, two large genomic segments in A12 and D12 homoeologous chromosomes (Ctg). To enable generation of longer contig sequences in assembly, we implemented a hybrid assembly method to process ~35x data from 454 technology and 2.8-3x data from Sanger method. Hybrid assemblies offered higher sequence coverage and better sequence assemblies. Homology studies revealed the presence of retrotransposon regions like Copia and Gypsy elements in these contigs and also helped in identifying new genomic SSRs. Unigenes were anchored to the sequences in Ctg-3301 and Ctg-465 to support the physical map. Gene density, gene structure and protein sequence information derived from protein prediction programs were used to obtain the functional annotation of these genes. Comparative analysis of both contigs with Arabidopsis genome exhibited synteny and microcollinearity with a conserved gene order in both genomes. This study provides insight about use of MTP-based BAC-pool sequencing approach for sequencing complex polyploid genomes with limited constraints in generating better sequence assemblies to build reference scaffold sequences. Combining the utilities of MTP-based BAC-pool sequencing with current longer and short read NGS technologies in multiplexed format would provide a new direction to cost-effectively and precisely sequence complex plant genomes.

An enhancer element harboring variants associated with systemic lupus erythematosus engages the TNFAIP3 promoter to influence A20 expression.

Functional characterization of causal variants present on risk haplotypes identified through genome-wide association studies (GWAS) is a primary objective of human genetics. In this report, we evaluate the function of a pair of tandem polymorphic dinucleotides, 42 kb downstream of the promoter of TNFAIP3, (rs148314165, rs200820567, collectively referred to as TT>A) recently nominated as causal variants responsible for genetic association of systemic lupus erythematosus (SLE) with tumor necrosis factor alpha inducible protein 3 (TNFAIP3). TNFAIP3 encodes the ubiquitin-editing enzyme, A20, a key negative regulator of NF-κB signaling. A20 expression is reduced in subjects carrying the TT>A risk alleles; however, the underlying functional mechanism by which this occurs is unclear. We used a combination of electrophoretic mobility shift assays (EMSA), mass spectrometry (MS), reporter assays, chromatin immunoprecipitation-PCR (ChIP-PCR) and chromosome conformation capture (3C) EBV transformed lymphoblastoid cell lines (LCL) from individuals carrying risk and non-risk TNFAIP3 haplotypes to characterize the effect of TT>A on A20 expression. Our results demonstrate that the TT>A variants reside in an enhancer element that binds NF-κB and SATB1 enabling physical interaction of the enhancer with the TNFAIP3 promoter through long-range DNA looping. Impaired binding of NF-κB to the TT>A risk alleles or knockdown of SATB1 expression by shRNA, inhibits the looping interaction resulting in reduced A20 expression. Together, these data reveal a novel mechanism of TNFAIP3 transcriptional regulation and establish the functional basis by which the TT>A risk variants attenuate A20 expression through inefficient delivery of NF-κB to the TNFAIP3 promoter. These results provide critical functional evidence supporting a direct causal role for TT>A in the genetic predisposition to SLE.

RBPJ mutations identified in two families affected by Adams-Oliver syndrome.

Through exome resequencing, we identified two unique mutations in recombination signal binding protein for immunoglobulin kappa J (RBPJ) in two independent families affected by Adams-Oliver syndrome (AOS), a rare multiple-malformation disorder consisting primarily of aplasia cutis congenita of the vertex scalp and transverse terminal limb defects. These identified mutations link RBPJ, the primary transcriptional regulator for the Notch pathway, with AOS, a human genetic disorder. Functional assays confirmed impaired DNA binding of mutated RBPJ, placing it among other notch-pathway proteins altered in human genetic syndromes.

WITHDRAWN: A female autoimmunity gene exists: DDX3X.

This article has been withdrawn at the request of the author(s) and/or editor. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at http://www.elsevier.com/locate/withdrawalpolicy.

Association of two independent functional risk haplotypes in TNIP1 with systemic lupus erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified >30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting NF-κB signaling. The present study was undertaken to investigate the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway. METHODS: We analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog TNIP2, in case-control populations of diverse ethnic origins. TNIP1, TNIP2, and TAX1BP1 were fine-mapped in a total of 8,372 SLE cases and 7,492 healthy controls from European-ancestry, African American, Hispanic, East Asian, and African American Gullah populations. Levels of TNIP1 messenger RNA (mRNA) and ABIN1 protein in Epstein-Barr virus-transformed human B cell lines were analyzed by quantitative reverse transcription-polymerase chain reaction and Western blotting, respectively. RESULTS: We found significant associations between SLE and genetic variants within TNIP1, but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified 2 independent risk haplotypes within TNIP1 in individuals of European ancestry that were also present in African American and Hispanic populations. Levels of TNIP1 mRNA and ABIN1 protein were reduced among subjects with these haplotypes, suggesting that they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression. CONCLUSION: Our results confirm the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.

Detection of members of the Secoviridae in the Tallgrass Prairie Preserve, Osage County, Oklahoma, USA.

Viruses are most frequently discovered because they cause disease. To expand knowledge of plant-associated viruses beyond these narrow constraints, non-cultivated plants of the Tallgrass Prairie of the United States were systematically surveyed for evidence of viruses. This report discusses putative viruses of the family Secoviridae identified by the survey. Sequence analysis suggests the presence of at least six viruses in the study site, including Bean pod mottle virus, Maize chlorotic dwarf virus, three previously undescribed viruses within the subfamily Comovirinae and one unclassifiable virus.

Identification of IRF8, TMEM39A, and IKZF3-ZPBP2 as susceptibility loci for systemic lupus erythematosus in a large-scale multiracial replication study.

Systemic lupus erythematosus (SLE) is a chronic heterogeneous autoimmune disorder characterized by the loss of tolerance to self-antigens and dysregulated interferon responses. The etiology of SLE is complex, involving both heritable and environmental factors. Candidate-gene studies and genome-wide association (GWA) scans have been successful in identifying new loci that contribute to disease susceptibility; however, much of the heritable risk has yet to be identified. In this study, we sought to replicate 1,580 variants showing suggestive association with SLE in a previously published GWA scan of European Americans; we tested a multiethnic population consisting of 7,998 SLE cases and 7,492 controls of European, African American, Asian, Hispanic, Gullah, and Amerindian ancestry to find association with the disease. Several genes relevant to immunological pathways showed association with SLE. Three loci exceeded the genome-wide significance threshold: interferon regulatory factor 8 (IRF8; rs11644034; p(meta-Euro) = 2.08 × 10(-10)), transmembrane protein 39A (TMEM39A; rs1132200; p(meta-all) = 8.62 × 10(-9)), and 17q21 (rs1453560; p(meta-all) = 3.48 × 10(-10)) between IKAROS family of zinc finger 3 (AIOLOS; IKZF3) and zona pellucida binding protein 2 (ZPBP2). Fine mapping, resequencing, imputation, and haplotype analysis of IRF8 indicated that three independent effects tagged by rs8046526, rs450443, and rs4843869, respectively, were required for risk in individuals of European ancestry. Eleven additional replicated effects (5 × 10(-8) < p(meta-Euro) < 9.99 × 10(-5)) were observed with CFHR1, CADM2, LOC730109/IL12A, LPP, LOC63920, SLU7, ADAMTSL1, C10orf64, OR8D4, FAM19A2, and STXBP6. The results of this study increase the number of confirmed SLE risk loci and identify others warranting further investigation.

The genomics of autoimmune disease in the era of genome-wide association studies and beyond.

Recent advances in the field of genetics have dramatically changed our understanding of autoimmune disease. Candidate gene and, more recently, genome-wide association (GWA) studies have led to an explosion in the number of loci and pathways known to contribute to autoimmune phenotypes. Since the 1970s, researchers have known that several alleles in the MHC region play a role in the pathogenesis of many autoimmune diseases. More recent work has identified numerous risk loci involving both the innate and adaptive immune responses. However, much remains to be learned about the heritability of autoimmune conditions. Most regions found through GWA scans have yet to isolate the association to the causal allele(s) responsible for conferring disease risk. A role for rare variants (allele frequencies of <1%) has begun to emerge. Future research will use next-generation sequencing (NGS) technology to comprehensively evaluate the human genome for risk variants. Whole-transcriptome sequencing is now possible, which will provide much more detailed gene expression data. The dramatic drop in the cost and time required to sequence the entire human genome will ultimately make it possible for this technology to be used as a clinical diagnostic tool.