Analysis of a functional BTNL2 polymorphism in type 1 diabetes, rheumatoid arthritis, and systemic lupus erythematosus

Missense Variant rs28362680 in BTNL2 Reduces Risk of Coronary Heart Disease

Background

The pathological basis of coronary heart disease (CHD) is atherosclerosis. BTNL2 can inhibit the activation of T cells. We aimed to explore the association between BTNL2 genetic variants and CHD risk in the southern Chinese Han population.

Methods

We recruited 1419 participants to perform an association analysis between missense variants in BTNL2 and CHD risk through SNPStats online software. Genotyping of all candidate SNPs were completed by the Agena MassARRAY. In addition, we used false-positive report probability analysis to detect whether the positive findings were noteworthy observations. We also used Haploview 4.2 software and SNPStats online software to conduct the haplotype analysis and analysis of linkage disequilibrium (LD). Finally, the interaction of SNP-SNP in CHD risk was evaluated by multi-factor dimensionality reduction (MDR).

Results

The results showed that BTNL2-rs35624343, -rs117896888, -rs41441651, -rs41417449, -rs28362680 and -rs2076523 were significantly associated with the CHD susceptibility. Especially for BTNL2-rs28362680, the allele A (OR = 0.68, p < 0.0001), genotype AA (OR = 0.40, p = 0.001) or GA (OR = 0.68, p < 0.0001) were associated with the reducing CHD risk. And -rs28362680 significantly reduced the CHD risk under all genetic models (dominant: OR = 0.64, p < 0.0001; recessive: OR = 0.47, p = 0.003; overdominant: OR = 0.73, p = 0.004; log-additive: OR = 0.66, p < 0.0001). And -rs28362680 was also closely associated with CHD risk reduction in all stratified analyses (age, gender, smoking, drinking, hypertension and diabetes). In addition, haplotype analysis showed that the “C_rs117896888C_rs41441651T_rs41417449A_rs28362680” (OR = 0.65, p < 0.0001) and “G_rs117896888T_rs41441651C_rs41417449A_rs28362680” (OR = 0.68, p = 0.013) may reduce CHD risk.

Conclusion

Missense variants (rs35624343, rs117896888, rs41441651, rs41417449, rs28362680, rs2076523) may be protective factors for the CHD risk. In particular, there were sufficient evidences that BTNL2-rs28362680 can protective CHD risk.

B7 immune checkpoint family members as putative therapeutics in autoimmune disease: An updated overview

Autoimmune diseases, especially among young people in the US, are one of the leading causes of morbidity and death. The immune responses are the fundamental pathogenicity of autoimmune disorders. The equilibrium between stimulatory and inhibitory signals is critical for the stimulation, migration, survival, and T cell-related immune responses. The B7 family can substantially regulate T cell-mediated immune responses. Nevertheless, recent breakthroughs in immune checkpoint blockade in cancer immunotherapy have facilitated autoimmune diseases, especially among the prone populations. In the current study, we tried to concisely review the role of the B7 family in regulating immune reactions and the influence of immune checkpoint inhibitors on autoimmunity development.

First-appearing islet autoantibodies for type 1 diabetes in young children: maternal life events during pregnancy and the child's genetic risk

AIMS/HYPOTHESIS

Psychological stress has long been considered a possible trigger of type 1 diabetes, although prospective studies examining the link between psychological stress or life events during pregnancy and the child's type 1 diabetes risk are rare. The objective of this study was to examine the association between life events during pregnancy and first-appearing islet autoantibodies (IA) in young children, conditioned by the child's type 1 diabetes-related genetic risk.

METHODS

The IA status of 7317 genetically at-risk The Environmental Determinants of Diabetes in the Young (TEDDY) participants was assessed every 3 months from 3 months to 4 years, and bi-annually thereafter. Reports of major life events during pregnancy were collected at study inception when the child was 3 months of age and placed into one of six categories. Life events during pregnancy were examined for association with first-appearing insulin (IAA) (N = 222) or GAD (GADA) (N = 209) autoantibodies in the child until 6 years of age using proportional hazard models. Relative excess risk due to interaction (RERI) by the child's HLA-DR and SNP profile was estimated.

RESULTS

Overall, 65% of mothers reported a life event during pregnancy; disease/injury (25%), serious interpersonal (28%) and job-related (25%) life events were most common. The association of life events during pregnancy differed between IAA and GADA as the first-appearing autoantibody. Serious interpersonal life events correlated with increased risk of GADA-first only in HLA-DR3 children with the BACH2-T allele (HR 2.28, p < 0.0001), an additive interaction (RERI 1.87, p = 0.0004). Job-related life events were also associated with increased risk of GADA-first among HLA-DR3/4 children (HR 1.53, p = 0.04) independent of serious interpersonal life events (HR 1.90, p = 0.002), an additive interaction (RERI 1.19, p = 0.004). Job-related life events correlated with reduced risk of IAA-first (HR 0.55, p = 0.004), particularly in children with the BTNL2-GG allele (HR 0.48; 95% CI 0.31, 0.76).

CONCLUSIONS/INTERPRETATION

Specific life events during pregnancy are differentially related to IAA vs GADA as first-appearing IA and interact with different HLA and non-HLA genetic factors, supporting the concept of different endotypes underlying type 1 diabetes. However, the mechanisms underlying these associations remain to be discovered. Life events may be markers for other yet-to-be-identified factors important to the development of first-appearing IA.

Lossless integration of multiple electronic health records for identifying pleiotropy using summary statistics

Increasingly, clinical phenotypes with matched genetic data from bio-bank linked electronic health records (EHRs) have been used for pleiotropy analyses. Thus far, pleiotropy analysis using individual-level EHR data has been limited to data from one site. However, it is desirable to integrate EHR data from multiple sites to improve the detection power and generalizability of the results. Due to privacy concerns, individual-level patients' data are not easily shared across institutions. As a result, we introduce Sum-Share, a method designed to efficiently integrate EHR and genetic data from multiple sites to perform pleiotropy analysis. Sum-Share requires only summary-level data and one round of communication from each site, yet it produces identical test statistics compared with that of pooled individual-level data. Consequently, Sum-Share can achieve lossless integration of multiple datasets. Using real EHR data from eMERGE, Sum-Share is able to identify 1734 potential pleiotropic SNPs for five cardiovascular diseases.

DNA methylation profiles are associated with complex regional pain syndrome after traumatic injury

Factors contributing to development of complex regional pain syndrome (CRPS) are not fully understood. This study examined possible epigenetic mechanisms that may contribute to CRPS after traumatic injury. DNA methylation profiles were compared between individuals developing CRPS (n = 9) and those developing non-CRPS neuropathic pain (n = 38) after undergoing amputation following military trauma. Linear Models for Microarray (LIMMA) analyses revealed 48 differentially methylated cytosine-phosphate-guanine dinucleotide (CpG) sites between groups (unadjusted P's < 0.005), with the top gene COL11A1 meeting Bonferroni-adjusted P < 0.05. The second largest differential methylation was observed for the HLA-DRB6 gene, an immune-related gene linked previously to CRPS in a small gene expression study. For all but 7 of the significant CpG sites, the CRPS group was hypomethylated. Numerous functional Gene Ontology-Biological Process categories were significantly enriched (false discovery rate-adjusted q value <0.15), including multiple immune-related categories (eg, activation of immune response, immune system development, regulation of immune system processes, and antigen processing and presentation). Differentially methylated genes were more highly connected in human protein-protein networks than expected by chance (P < 0.05), supporting the biological relevance of the findings. Results were validated in an independent sample linking a DNA biobank with electronic health records (n = 126 CRPS phenotype, n = 19,768 non-CRPS chronic pain phenotype). Analyses using PrediXcan methodology indicated differences in the genetically determined component of gene expression in 7 of 48 genes identified in methylation analyses (P's < 0.02). Results suggest that immune- and inflammatory-related factors might confer risk of developing CRPS after traumatic injury. Validation findings demonstrate the potential of using electronic health records linked to DNA for genomic studies of CRPS.

Genome-Wide Association Study of Lean Nonalcoholic Fatty Liver Disease Suggests Human Leukocyte Antigen as a Novel Candidate Locus

Nonalcoholic fatty liver disease (NAFLD) is supposed to manifest its metabolic phenotype in the liver, but it is common to have lean individuals diagnosed with NAFLD, known as lean NAFLD. We conducted a two-stage analysis to identify NAFLD-associated loci in Japanese patients. In stage I, 275 metabolically healthy normal-weight patients with NAFLD were compared with 1,411 non-NAFLD controls adjusted for age, sex, and alcohol consumption by a genome-wide association study (GWAS). In stage II, human leukocyte antigen (HLA) in chromosome 6 (chr6) (P = 6.73E-08), microRNA (MIR) MIR548F3 in chr7 (P = 4.25E-07), myosin light chain 2 (MYL2) in chr12 (P = 4.39E-07), and glycoprotein precursor (GPC)6 in chr13 (P = 5.43E-07), as suggested by the GWAS, were assessed by single nucleotide polymorphism (SNP) association analysis of whole NAFLD against non-NAFLD in 9,726 members of the general population. A minor allele of the secondary lead SNP in chr6, rs2076529, was significantly associated (odds ratio [OR], 1.19; 95% confidence interval [CI], 1.11-1.28; P = 2.10E-06) and the lead SNP in chr7 was weakly associated (OR 1.15; 95% CI, 1.04-1.27; P = 6.19E-03) with increased NAFLD risk. Imputation-based typing of HLA showed a significant difference in the distribution of HLA-B, HLA-DR-beta chain 1 (DRB1), and HLA-DQ-beta chain 1 (DQB1) alleles in lean NAFLD GWAS. Next-generation sequence-based typing of HLA in 5,649 members of the general population replicated the significant difference of HLA-B allele distribution and the significant increase of the HLA-B*54:01 allele in whole NAFLD. Fecal metagenomic analysis of 3,420 members of the general population showed significant dissimilarity in beta-diversity analysis of rs2076529 and HLA-B*54:01 allele carriers from noncarriers. Veillonellaceae was increased but Verrucomicrobia was decreased in rs2076529 minor allele and HLA-B*54:01 allele carriers as in NAFLD. Conclusion: HLA was identified as a novel locus associated with NAFLD susceptibility, which might be affected by the alteration of gut microbiota.

Integration of single nucleotide variants and whole-genome DNA methylation profiles for classification of rheumatoid arthritis cases from controls

This study evaluated the use of multiomics data for classification accuracy of rheumatoid arthritis (RA). Three approaches were used and compared in terms of prediction accuracy: (1) whole-genome prediction (WGP) using SNP marker information only, (2) whole-methylome prediction (WMP) using methylation profiles only, and (3) whole-genome/methylome prediction (WGMP) with combining both omics layers. The number of SNP and of methylation sites varied in each scenario, with either 1, 10, or 50% of these preselected based on four approaches: randomly, evenly spaced, lowest p value (genome-wide association or epigenome-wide association study), and estimated effect size using a Bayesian ridge regression (BRR) model. To remove effects of high levels of pairwise linkage disequilibrium (LD), SNPs were also preselected with an LD-pruning method. Five Bayesian regression models were studied for classification, including BRR, Bayes-A, Bayes-B, Bayes-C, and the Bayesian LASSO. Adjusting methylation profiles for cellular heterogeneity within whole blood samples had a detrimental effect on the classification ability of the models. Overall, WGMP using Bayes-B model has the best performance. In particular, selecting SNPs based on LD-pruning with 1% of the methylation sites selected based on BRR included in the model, and fitting the most significant SNP as a fixed effect was the best method for predicting disease risk with a classification accuracy of 0.975. Our results showed that multiomics data can be used to effectively predict the risk of RA and identify cases in early stages to prevent or alter disease progression via appropriate interventions.

The Role of Immune Checkpoint Receptors in Regulating Immune Reactivity in Lupus

Immune checkpoint receptors with co-stimulatory and co-inhibitory signals are important modulators for the immune system. However, unrestricted co-stimulation and/or inadequate co-inhibition may cause breakdown of self-tolerance, leading to autoimmunity. Systemic lupus erythematosus (SLE) is a complex multi-organ disease with skewed and dysregulated immune responses interacting with genetics and the environment. The close connections between co-signaling pathways and SLE have gradually been established in past research. Also, the recent success of immune checkpoint blockade in cancer therapy illustrates the importance of the co-inhibitory receptors in cancer immunotherapy. Moreover, immune checkpoint blockade could result in substantial immune-related adverse events that mimic autoimmune diseases, including lupus. Together, immune checkpoint regulators represent viable immunotherapeutic targets for the treatment of both autoimmunity and cancer. Therefore, it appears reasonable to treat SLE by restoring the out-of-order co-signaling axis or by manipulating collateral pathways to control the pathogenic immune responses. Here, we review the current state of knowledge regarding the relationships between SLE and the co-signaling pathways of T cells, B cells, dendritic cells, and neutrophils, and highlight their potential clinical implications. Current clinical trials targeting the specific co-signaling axes involved in SLE help to advance such knowledge, but further in-depth exploration is still warranted.

BTNL2-Ig Protein Attenuates Type 1 Diabetes in Non-Obese Diabetic (NOD) Mice

Type 1 diabetes (T1D) is a T cell-mediated autoimmune disease in which insulin-producing β-cells are destroyed. Although butyrophilin-like 2 (BTNL2) has been shown to be a negative T cell regulator in vitro, its ability to inhibit T cell responses in vivo has not been determined. In this study, the effect of a recombinant BTNL2-IgG2a Fc (rBTNL2-Ig) fusion protein on T1D development in vivo is determined. It is shown here that in vivo administration of rBTNL2-Ig ameliorates T1D in non-obese diabetic (NOD) mice. This is associated with the ability of rBTNL2-Ig to inhibit the proliferation, activation, and inflammatory cytokine production from autoreactive T cells in vivo. In addition, rBTNL2-Ig treatment increases the generation of regulatory T cells. The results suggest that targeting the BTNL2 pathway has the potential to be used in the prevention and treatment of patients with T1D.

Genetic Contribution to the Divergence in Type 1 Diabetes Risk Between Children From the General Population and Children From Affected Families

The risk for autoimmunity and subsequently type 1 diabetes is 10-fold higher in children with a first-degree family history of type 1 diabetes (FDR children) than in children in the general population (GP children). We analyzed children with high-risk HLA genotypes (n = 4,573) in the longitudinal TEDDY birth cohort to determine how much of the divergent risk is attributable to genetic enrichment in affected families. Enrichment for susceptible genotypes of multiple type 1 diabetes-associated genes and a novel risk gene, BTNL2, was identified in FDR children compared with GP children. After correction for genetic enrichment, the risks in the FDR and GP children converged but were not identical for multiple islet autoantibodies (hazard ratio [HR] 2.26 [95% CI 1.6-3.02]) and for diabetes (HR 2.92 [95% CI 2.05-4.16]). Convergence varied depending upon the degree of genetic susceptibility. Risks were similar in the highest genetic susceptibility group for multiple islet autoantibodies (14.3% vs .12.7%) and diabetes (4.8% vs. 4.1%) and were up to 5.8-fold divergent for children in the lowest genetic susceptibility group, decreasing incrementally in GP children but not in FDR children. These findings suggest that additional factors enriched within affected families preferentially increase the risk of autoimmunity and type 1 diabetes in lower genetic susceptibility strata.

Comparative study for haplotype block partitioning methods - Evidence from chromosome 6 of the North American Rheumatoid Arthritis Consortium (NARAC) dataset

Haplotype-based methods compete with “one-SNP-at-a-time” approaches on being preferred for association studies. Chromosome 6 contains most of the known genetic biomarkers for rheumatoid arthritis (RA) disease. Therefore, chromosome 6 serves as a benchmark for the haplotype methods testing. The aim of this study is to test the North American Rheumatoid Arthritis Consortium (NARAC) dataset to find out if haplotype block methods or single-locus approaches alone can sufficiently provide the significant single nucleotide polymorphisms (SNPs) associated with RA. In addition, could we be satisfied with only one method of the haplotype block methods for partitioning chromosome 6 of the NARAC dataset? In the NARAC dataset, chromosome 6 comprises 35,574 SNPs for 2,062 individuals (868 cases, 1,194 controls). Individual SNP approach and three haplotype block methods were applied to the NARAC dataset to identify the RA biomarkers. We employed three haplotype partitioning methods which are confidence interval test (CIT), four gamete test (FGT), and solid spine of linkage disequilibrium (SSLD). P-values after stringent Bonferroni correction for multiple testing were measured to assess the strength of association between the genetic variants and RA susceptibility. Moreover, the block size (in base pairs (bp) and number of SNPs included), number of blocks, percentage of uncovered SNPs by the block method, percentage of significant blocks from the total number of blocks, number of significant haplotypes and SNPs were used to compare among the three haplotype block methods. Individual SNP, CIT, FGT, and SSLD methods detected 432, 1,086, 1,099, and 1,322 associated SNPs, respectively. Each method identified significant SNPs that were not detected by any other method (Individual SNP: 12, FGT: 37, CIT: 55, and SSLD: 189 SNPs). 916 SNPs were discovered by all the three haplotype block methods. 367 SNPs were discovered by the haplotype block methods and the individual SNP approach. The P-values of these 367 SNPs were lower than those of the SNPs uniquely detected by only one method. The 367 SNPs detected by all the methods represent promising candidates for RA susceptibility. They should be further investigated for the European population. A hybrid technique including the four methods should be applied to detect the significant SNPs associated with RA for chromosome 6 of the NARAC dataset. Moreover, SSLD method may be preferred for its favored benefits in case of selecting only one method.

In vivo administration of recombinant BTNL2-Fc fusion protein ameliorates graft-versus-host disease in mice

Although hematopoietic stem cell transplantation (HSCT) has been widely used in the treatment of many diseases, graft-versus-host disease (GVHD) remains a major complication after allogeneic HSCT. Butyrophilin-like 2 (BTNL2) protein has been reported to have the ability to inhibit T cell proliferation in vitro; its ability to inhibit T cell responses in vivo has not been determined. We show here that in vivo administration of recombinant BTNL2-IgG2a Fc (rBTNL2-Ig) fusion protein ameliorates GVHD in mice. This is related to the ability of rBTNL2-Ig to inhibit T cell proliferation, activation and Th1/Th17 cytokine production in vivo. Furthermore, rBTNL2-Ig treatment increases the generation of regulatory T cells. Our results suggest that rBTNL2-Ig has the potential to be used in the prevention and treatment of patients with GVHD.

Deep Resequencing of Ulcerative Colitis-Associated Genes Identifies Novel Variants in Candidate Genes in the Korean Population

BACKGROUND

Genome-wide association studies and meta-analyses have revealed the genetic background of ulcerative colitis (UC) by identifying common variants. However, these variants do not fully explain the disease variance in UC. To identify novel variants, we performed deep resequencing of UC-associated genes in Korean UC patients and subsequently investigated the functional roles of identified susceptibility genes.

METHODS

We performed targeted deep resequencing of 108 genes in 24 Korean UC patients and then performed association analysis with data from 126 healthy controls. We validated these variants using 2-stage replication studies including 793 UC patients and 783 controls. We performed in silico and pathway analyses and functional analyses.

RESULTS

The combined analysis including 2 replication studies identified 6 novel susceptibility loci and reconfirmed 10 previously reported loci. Among the novel single nucleotide variants (SNVs), rs10035653 in C5orf55 (P = 2.08 × 10-3; OR = 1.50), rs41417449 in BTNL2 (P = 1.27 × 10-2; OR = 1.32), rs3117099 in HCG23 (P = 9.98 × 10-6; OR = 1.40), rs7192 in HLA-DRA (P = 6.95 × 10-9; OR = 1.57), and rs3744246 in ORMDL3 (P = 2.21 × 10-2; OR = 1.21) were identified as causal variants, whereas rs713669 in IL17REL (P = 2.69 × 10-2; OR = 0.84) as a protective variant for UC. When correcting multiple testing, 3 novel SNVs (rs41417449 in BTNL2, rs3744246 in ORMDL3, and rs713669 in IL17REL) and 4 previously reported SNVs did not reach a statistical significance. Functional study suggested that SNVs of BTNL2 and C5orf55 exacerbated the inflammatory response both in vitro and in vivo.

CONCLUSIONS

This study identified 3 novel susceptibility loci and validated 6 previously reported SNVs for UC through deep resequencing in Koreans and revealed the functional roles of BTNL2 and C5orf55.

Prediction of Disordered Regions and Their Roles in the Anti-Pathogenic and Immunomodulatory Functions of Butyrophilins

Butyrophilins (BTNs) are a group of the moonlighting proteins, some members of which are secreted in milk. They constitute a large family of structurally similar type 1 transmembrane proteins from the immunoglobulin superfamily. Although the founding member of this family is related to lactation, participating in the secretion, formation and stabilization of milk fat globules, it may also have a cell surface receptor function. Generally, the BTN family members are known to modulate co-stimulatory responses, T cell selection, differentiation, and cell fate determination. Polymorphism of these genes was shown to be associated with the pathology of several human diseases. Despite their biological significance, structural information on human butyrophilins is rather limited. Based on their remarkable multifunctionality, butyrophilins seem to belong to the category of moonlighting proteins, which are known to contain intrinsically disordered protein regions (IDPRs). However, the disorder status of human BTNs was not systematically investigated as of yet. The goal of this study is to fill this gap and to evaluate peculiarities of intrinsic disorder predisposition of the members of human BTN family, and to find if they have IDPRs that can be attributed to the multifunctionality of these important proteins.

Butyrophilins: an important new element of resistance

Butyrophilins belonging to the immunoglobulin superfamily are new immune system regulators because they are present on lymphocytes, dendritic cells, monocytes, macrophages, neutrophils and eosinophils, and they exert a stimulatory and (or) inhibitory effect on them. The role of butyrophilins is associated and results from their similarity to the regulatory B7 protein family involved in the modulation of immune phenomena. Butyrophilins are glycoproteins built of two extracellular immunoglobulin domains, stabilized with disulfide bonds: constant IgC, and variable IgV and a transmembrane region. Most of these proteins contain a conserved domain encoded by a single exon – B30.2, also referred to as PRYSPRY. In humans, the family of butyrophilins includes 7 butyrophilin proteins, 5 butyrophilin-like proteins and the SKINT-like factor. Butyrophilins have been also demonstrated to play a role in various infections, e.g. tuberculosis or diseases that include sarcoidosis, systemic lupus erythematosus, rheumatoid arthritis, genetic metabolic diseases, ulcerative colitis, cancer and kidney disease.

The BTNL2 G16071A gene polymorphism increases granulomatous disease susceptibility: A meta-analysis including FPRP test of 8710 participants

OBJECTIVE

The butyrophilin-like 2 (BTNL2) G16071A gene polymorphism has been implicated in the susceptibility to granulomatous diseases, but the results were inconclusive. The objective of the current study was to precisely explore the relationship between BTNL2 G16071A gene polymorphism and granulomatous disease susceptibility by the meta-analysis including false-positive report probability (FPRP) test.

METHODS

A systematic literature search in the PubMed, Embase, and Wanfang databases, China National Knowledge Internet, and commercial Internet search engines was conducted to identify studies published up to April 1, 2016. The odds ratio (OR) with 95% confidence interval (CI) was used to assess the effect size. Statistical analysis was conducted using the STATA 12.0 software and FPRP test sheet.

RESULTS

In total, all 4324 cases and 4386 controls from 14 eligible studies were included in the current meta-analysis. By the overall meta-analysis, we found a significant association between BTNL2 G16071A gene polymorphism and granulomatous disease susceptibility (A vs G: OR = 1.25, 95% CI = 1.07-1.45, P = 0.005). The meta-regression analyses showed that a large proportion of the between-study heterogeneity was significantly attributed to the ethnicity (A vs G, P = 0.013) and the types of granulomatous diseases (A vs G, P = 0.002). By the subgroup meta-analysis, the BTNL2 G16071A gene polymorphism was associated with granulomatous disease susceptibility in Caucasians (A vs G: OR = 1.37, 95% CI = 1.18-1.58, P < 0.001). Moreover, a significant relationship between the BTNL2 G16071A gene polymorphism and sarcoidosis susceptibility (A vs G: OR = 1.52, 95% CI = 1.39-1.66, P < 0.001) was found. However, to avoid the "false-positive report," we further investigated the significant associations observed in the present meta-analysis by the FPRP test. Interestingly, the results of FPRP test indicated that the BTNL2 G16071A gene polymorphism was truly associated with sarcoidosis susceptibility (A vs G, FPRP < 0.001). Additionally, the FPRP test confirmed that the BTNL2 G16071A gene polymorphism was associated only with granulomatous disease susceptibility among Caucasians (A vs G, FPRP < 0.001) at the level of a prior probability, which was 0.001.

CONCLUSION

The meta-analysis indicated that BTNL2 G16071A gene polymorphism may as a likelihood factor contributed to granulomatous disease susceptibility, especially increasing the sarcoidosis susceptibility. In addition, the polymorphism may be greatly associated with likelihood of granulomatous diseases among Caucasians.

Constraints on Biological Mechanism from Disease Comorbidity Using Electronic Medical Records and Database of Genetic Variants

Patterns of disease co-occurrence that deviate from statistical independence may represent important constraints on biological mechanism, which sometimes can be explained by shared genetics. In this work we study the relationship between disease co-occurrence and commonly shared genetic architecture of disease. Records of pairs of diseases were combined from two different electronic medical systems (Columbia, Stanford), and compared to a large database of published disease-associated genetic variants (VARIMED); data on 35 disorders were available across all three sources, which include medical records for over 1.2 million patients and variants from over 17,000 publications. Based on the sources in which they appeared, disease pairs were categorized as having predominant clinical, genetic, or both kinds of manifestations. Confounding effects of age on disease incidence were controlled for by only comparing diseases when they fall in the same cluster of similarly shaped incidence patterns. We find that disease pairs that are overrepresented in both electronic medical record systems and in VARIMED come from two main disease classes, autoimmune and neuropsychiatric. We furthermore identify specific genes that are shared within these disease groups.

Genetic Susceptibility to Vitiligo: GWAS Approaches for Identifying Vitiligo Susceptibility Genes and Loci

Vitiligo is an autoimmune disease with a strong genetic component, characterized by areas of depigmented skin resulting from loss of epidermal melanocytes. Genetic factors are known to play key roles in vitiligo through discoveries in association studies and family studies. Previously, vitiligo susceptibility genes were mainly revealed through linkage analysis and candidate gene studies. Recently, our understanding of the genetic basis of vitiligo has been rapidly advancing through genome-wide association study (GWAS). More than 40 robust susceptible loci have been identified and confirmed to be associated with vitiligo by using GWAS. Most of these associated genes participate in important pathways involved in the pathogenesis of vitiligo. Many susceptible loci with unknown functions in the pathogenesis of vitiligo have also been identified, indicating that additional molecular mechanisms may contribute to the risk of developing vitiligo. In this review, we summarize the key loci that are of genome-wide significance, which have been shown to influence vitiligo risk. These genetic loci may help build the foundation for genetic diagnosis and personalize treatment for patients with vitiligo in the future. However, substantial additional studies, including gene-targeted and functional studies, are required to confirm the causality of the genetic variants and their biological relevance in the development of vitiligo.

Novel Immune Check-Point Regulators in Tolerance Maintenance

The great success of anti-cytotoxic lymphocyte antigen 4 (CTLA4) and anti-programed cell death protein 1 (PD1) in cancer treatment has encouraged more effort in harnessing the immune response through immunomodulatory molecules in various diseases. The immunoglobulin (Ig) super family comprises the majority of immunomodulatory molecules. Discovery of novel Ig super family members has brought novel insights into the function of different immune cells in tolerance maintenance. In this review, we discuss the function of newly identified B7 family molecules, B7-H4 and V-domain Ig Suppressor of T cell Activation (VISTA), and the butyrophilin/butyrophilin-like family members. We discuss the current stages of immunomodulatory molecules in clinical trials of organ transplantation. The potential of engaging the novel Ig superfamily members in tolerance maintenance is also discussed. We conclude with the challenges remaining to manipulate these molecules in the immune response.

Knowledge-based analysis of genetic associations of rheumatoid arthritis to inform studies searching for pleiotropic genes: a literature review and network analysis

Introduction

Pleiotropy describes the genetic effect of a single gene on multiple phenotypic traits. Gene variants directly affect the normal processes of a series of physiological and biochemical reactions, and therefore cause a variety of diseases traits to be changed accordingly. Moreover, a shared genetic susceptibility mechanism may exist between different diseases. Therefore, shared genes, with pleiotropic effects, are important to understand the sharing pathogenesis and hence the mechanisms underlying comorbidity.

Methods

In this study, we proposed combining genome-wide association studies (GWAS) and public knowledge databases to search for potential pleiotropic genes associated with rheumatoid arthritis (RA) and eight other related diseases. Here, a GWAS-based network analysis is used to recognize risk genes significantly associated with RA. These RA risk genes are re-extracted as potential pleiotropic genes if they have been proved to be susceptible genes for at least one of eight other diseases in the OMIM or PubMed databases.

Results

In total, we extracted 116 potential functional pleiotropic genes for RA and eight other diseases, including five hub pleiotropic genes, BTNL2, HLA-DRA, NOTCH4, TNXB, and C6orf10, where BTNL2, NOTCH4, and C6orf10 are novel pleiotropic genes identified by our analysis.

Conclusions

This study demonstrates that pleiotropy is a common property of genes associated with disease traits. Our results ascertained the shared genetic risk profiles that predisposed individuals to RA and other diseases, which could have implications for identification of molecular targets for drug development, and classification of diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0715-1) contains supplementary material, which is available to authorized users.

The T-Cell Inhibitory Molecule Butyrophilin-Like 2 Is Up-regulated in Mild Plasmodium falciparum Infection and Is Protective During Experimental Cerebral Malaria

Plasmodium falciparum infection can result in severe disease that is associated with elevated inflammation and vital organ dysfunction; however, malaria-endemic residents gain protection from lethal outcomes and manifest only mild symptoms during infection. To characterize host responses associated with this more effective antimalarial response, we characterized whole-blood transcriptional profiles in Rwandan adults during a mild malaria episode and compared them with findings from a convalescence sample. We observed transcriptional up-regulation in many pathways, including type I interferon, interferon γ, complement activation, and nitric oxide during malaria infection, which provide benchmarks of mild disease physiology. Transcripts encoding negative regulators of T-cell activation, such as programmed death ligand 1 (PD-L1), programmed death 1 ligand 2 (PD-L2), and the butyrophilin family member butyrophilin-like 2 (BTNL2) were also increased. To support an important functional role for BTNL2 during malaria infection, we studied chimeric mice reconstituted with BTNL2(-/-) or wild-type hematopoietic cells that were inoculated with Plasmodium berghei ANKA, a murine model of cerebral malaria. We found that BTNL2(-/-) chimeric mice had a significant decrease in survival compared with wild-type counterparts. Collectively these data characterize the immune responses associated with mild malaria and uncover a novel role for BTNL2 in the host response to malaria.

The NeST long ncRNA controls microbial susceptibility and epigenetic activation of the interferon-γ locus

Long noncoding RNAs (lncRNAs) are increasingly appreciated as regulators of cell-specific gene expression. Here, an enhancer-like lncRNA termed NeST (nettoie Salmonella pas Theiler's [cleanup Salmonella not Theiler's]) is shown to be causal for all phenotypes conferred by murine viral susceptibility locus Tmevp3. This locus was defined by crosses between SJL/J and B10.S mice and contains several candidate genes, including NeST. The SJL/J-derived locus confers higher lncRNA expression, increased interferon-γ (IFN-γ) abundance in activated CD8(+) T cells, increased Theiler's virus persistence, and decreased Salmonella enterica pathogenesis. Transgenic expression of NeST lncRNA alone was sufficient to confer all phenotypes of the SJL/J locus. NeST RNA was found to bind WDR5, a component of the histone H3 lysine 4 methyltransferase complex, and to alter histone 3 methylation at the IFN-γ locus. Thus, this lncRNA regulates epigenetic marking of IFN-γ-encoding chromatin, expression of IFN-γ, and susceptibility to a viral and a bacterial pathogen.

Genetic testing in diffuse parenchymal lung disease

Diffuse parenchymal lung diseases (DPLD) represent a diverse group of disorders affecting the distal lung parenchyma, specifically the tissue and spaces surrounding the alveoli, which may be filled with inflammatory cells, proliferating fibroblasts or established fibrosis, often leading to architectural distortion and impaired gas exchange. While the underlying pathogenetic mechanisms are known or inferred for some DPLD (such as sarcoidosis, silicosis, drug reactions and collagen vascular diseases), the pathogenesis of the majority of these entities - particularly those characterized by progressive fibrosis - is poorly understood.

Several lines of evidence indicate that the development of pulmonary fibrosis is genetically determined. They include: 1. familial clustering; 2. the occurrence of pulmonary fibrosis in the context of rare inherited disorders; 3. substantial variability in the development of pulmonary fibrosis amongst individuals exposed to organic or inorganic dusts; 4. difference in susceptibility to fibrogenic stimuli amongst inbred strains of mice.

This review focuses on idiopathic pulmonary fibrosis (IPF) and sarcoidosis, the two most common DPLD and the two entities for which there is stronger evidence of a genetic predisposition, although how aberrant genes interact with each other and with environmental factors, such as smoking in IPF and infectious agents in sarcoidosis, in determining disease susceptibility and clinical phenotypes is largely unknown. Finally, we discuss practical issues and implications for both patients and physicians of recent advances in the genetics of sarcoidosis and IPF.

Clinical profile of coexisting conditions in type 1 diabetes mellitus patients

AIMS

Type 1 diabetes mellitus (T1DM) is associated with various genetic and autoimmune diseases implicated in its etiopathogenesis. We hereby profile the clinical association of such diseases among patients from our center.

METHODS

Consecutive patients of T1DM presenting to department of Endocrinology from May 1997 to December 2011 were retrospectively analyzed in context of associated clinical profile.

RESULTS

Among 260 patients diagnosed as T1DM, 21 (8%) had hypothyroidism, 4 (1.5%) had hyperthyroidism and 2 (0.7%) had primary adrenal insufficiency. Eighteen patients (7%) had celiac disease, 9 (3.5%) had Turner's syndrome, 5 patients (1.9%) had Klinefelter's syndrome, whereas Down's syndrome and Noonan's syndrome was present in 2 and 1 patients (0.7%) respectively. One patient had Wolframs' syndrome and 1 patients had myasthenia gravis. Systemic lupus erythematosus and rheumatoid arthritis were present in 3 and 1 patients respectively. Total of 5 patients with cerebral palsy, 4 cases with deaf mutism, 4 cases with acute psychosis and 16 patients with depression were noted. Mean age of study patients was 20.8±9.8 years (range, 3-23 years).

CONCLUSION

Various conditions including genetic (Down, Turner, Noonan, and Klinefelter's), autoimmune (thyroid and adrenal disorders, myasthenia gravis, SLE, rheumatoid arthritis) and central nervous system diseases were the associated diseases encountered in our patients. Routine screening is required for early diagnosis and treatment of associated co morbidities.

Butyrophilins: an emerging family of immune regulators

Butyrophilins (Btns) and butyrophilin-like (Btnl) molecules are emerging as novel regulators of immune responses in mice and humans. Several clues point to their probable importance: many of the genes are located within the MHC; they are structurally related to B7-co-stimulatory molecules; they are functionally implicated in T cell inhibition and in the modulation of epithelial cell-T cell interactions; and they are genetically associated with inflammatory diseases. Nonetheless, initial immersion into the current literature can uncover confusion over even basic information such as gene names and expression patterns, and seemingly conflicting data regarding the biological activities of different family members. This review addresses each of these issues, concluding with the attractive potential of Btn and Btnl molecules to act as specific attenuators of tissue-associated inflammatory responses.

The BTNL2 A allele variant is frequent in Danish patients with sarcoidosis

BACKGROUND

The butyrophilin-like 2 (BTNL2) gene is located on chromosome 6p21.3 close to the HLA-class II genes. An association has been reported between sarcoidosis and a single nucleotide polymorphism in BTNL2, rs2076530, also termed the A allele.

OBJECTIVES

To evaluate whether patients with sarcoidosis carry the A allele more frequently than healthy subjects.

METHODS

The series comprised 87 ethnic Danes with sarcoidosis and 113 healthy control subjects. Analysis of rs2076530 was performed by Taqman assay, polymerase chain reaction and sequencing of genomic DNA.

RESULTS

Sarcoidosis patients had a higher frequency of the A allele than controls (73.9% vs 55.8%) (P < 0.025). The frequency of GG, GA and AA genotype was 5.7%, 40.2% and 54.0% in patients vs 16.0%, 56.6% and 27.4% in controls (P < 0.001). The AA genotype was associated with increased risk of sarcoidosis in both a dominant [odds ratio (OR) 3.1; 95% confidence interval (CI) 1.1-8.7; P < 0.03] and a recessive model (OR 3.1; 95% CI 1.72-5.61; P < 0.001). Population attributable fraction for disease was 50% in a dominant model and 25% in a recessive model.

CONCLUSIONS

The BTNL2 A allele variant occurs with a high frequency in Danish patients with sarcoidosis and the AA genotype is associated with a ∼threefold higher risk of sarcoidosis than the GG genotype. Our results should encourage future studies on the interrelationship between the BTNL2 protein and granuloma formation in sarcoidosis.

The role of bioinformatics in studying rheumatic and autoimmune disorders

In the past decade, the availability and abundance of individual-level molecular data, such as gene expression, proteomics and sequence data, has enabled the use of integrative computational approaches to pose and answer novel questions about disease. In this article, we discuss several examples of applications of bioinformatics techniques to study autoimmune and rheumatic disorders. We focus our discussion on how integrative techniques can be applied to analyze gene expression and genetic variation data across different diseases, and discuss the implications of such analyses. We also outline current challenges and future directions of these approaches. We show that integrative computational methods are essential for translational research and provide a powerful opportunity to improve human health by refining the current knowledge about diagnostics, therapeutics and mechanisms of disease pathogenesis.

Genome-wide analysis identifies a quantitative trait locus in the MHC class II region associated with generalized vitiligo age of onset

Generalized vitiligo is a common autoimmune disease in which acquired patchy depigmentation of skin, hair, and mucous membranes results from loss of melanocytes from involved areas. Previous genetic analyses have focused on vitiligo susceptibility, and have identified a number of genes involved in disease risk. Age of onset of generalized vitiligo also involves a substantial genetic component, but has not previously been studied systematically. In this study, we report a genome-wide association study of vitiligo age of onset in 1,339 generalized vitiligo patients, with replication in an independent cohort of 677 cases. We identified a quantitative trait locus for vitiligo age of onset in the major histocompatibility complex (MHC) class II region, located near c6orf10-BTNL2 (rs7758128; P=8.14 × 10(-11)), a region that is also associated with generalized vitiligo susceptibility. In contrast, there was no association of vitiligo age of onset with any other MHC or non-MHC loci that are associated with vitiligo susceptibility. These findings highlight the differing roles played by genes involved in vitiligo susceptibility versus vitiligo age of onset, and illustrate that genome-wide analyses can be used to identify genes involved in quantitative aspects of disease natural history, as well as disease susceptibility per se.

Butyrophilin-like 2 in pulmonary sarcoidosis: a factor for susceptibility and progression?

The aims of this study were to assess the association of BTNL2G16071A with the course of pulmonary sarcoidosis and to verify the association with disease predisposition. In addition, the linkage between BTNL2G16071A and certain human leukocyte antigen (HLA)-DRB1/DQB1 types was investigated. In a retrospective case-control study BTNL2G16071A, HLA-DQB1, and HLA-DRB1 were typed in 632 sarcoidosis patients. These patients were classified into 304 patients with persistent sarcoidosis and 328 patients with nonpersistent sarcoidosis. The BTNL2 16071A variant allele was present significantly more often in patients with persistent disease (92.4%; 281/304) compared with patients demonstrating a nonpersistent course (86.6%; 284/328; odds ratio (OR) = 1.89 with 95% confidence interval (95% CI) 1.11-3.22). Furthermore, BTNL2 16071A variant allele carriers have an increased risk (OR = 1.85, 95% CI 1.19-2.88) of developing sarcoidosis. Moreover, the strong linkage between variant allele and HLA-DRB1*15 presence (OR = 8.43, 95% CI 3.02-23.5) was confirmed. The presence of a BTNL2G16071A variant allele almost doubles the risk of progressing to persistent pulmonary sarcoidosis in addition to increasing the risk of developing sarcoidosis. Presumably, these increased risks are caused by the strong linkage between BTNL2G16071A and DRB1*15. The choice between determining BTNL2G16071A SNP or the HLA-DRB1 type depends on the ability and/or availability to perform either test.

Association of a polymorphism of BTN2A1 with chronic kidney disease in individuals with or without hypertension or diabetes mellitus

Hypertension and diabetes mellitus are important risk factors for chronic kidney disease (CKD). We previously showed that the C→T polymorphism (rs6929846) of BTN2A1 was significantly associated with myocardial infarction. The purpose of the present study was to examine an association of rs6929846 of BTN2A1 with CKD in individuals with or without hypertension or diabetes mellitus, thereby contributing to the personalized prevention of CKD in such individuals separately. The study population comprised 7,542 unrelated individuals, including 2,289 subjects with CKD [estimated glomerular filtration rate (eGFR) <60 ml/min/1.73 m(2)] and 5,253 controls (eGFR ≥60 ml/min/1.73 m(2)) with or without hypertension or diabetes mellitus. The Chi-square test, a multivariable logistic regression analysis with adjustment for covariates, as well as a stepwise forward selection procedure revealed that the C→T polymorphism (rs6929846) of BTN2A1 was significantly associated with CKD in normotensive individuals, in diabetic individuals and in individuals with hypertension and diabetes mellitus, or without either condition, with the T allele representing a risk factor for CKD. Stratification of subjects based on hypertension or diabetes mellitus may thus be important in order to achieve personalized prevention of CKD with the use of genetic information.

Association of a polymorphism of BTN2A1 with myocardial infarction in East Asian populations

OBJECTIVE

We have performed a genome-wide association study (GWAS) to identify genetic variants that confer susceptibility to myocardial infarction (MI) in Japanese and Korean populations.

METHODS

A total of 17,447 Japanese or Korean individuals from four independent subject panels was examined. Japanese subject panels A, B, and C comprised 134 individuals with MI and 137 controls, 1431 individuals with MI and 3161 controls, and 643 individuals with MI and 1347 controls, respectively, whereas the Korean population comprised 1880 individuals with MI and 8714 controls. A GWAS for MI was performed in Japanese subject panel A with the use of the Affymetrix GeneChip Human Mapping 500K Array Set.

RESULTS

Seventy single nucleotide polymorphisms (SNPs) significantly (P<1.0×10(-7)) associated with MI by the GWAS were examined further in Japanese subject panel B, revealing two SNPs (rs6929846 of BTN2A1, rs2569512 of ILF3) to be significantly (P<0.0007) associated with MI. The rs6929846 SNP of BTN2A1, but not rs2569512 of ILF3, was also significantly associated with MI in Japanese subject panel C. However, the association of neither rs6929846 nor rs2569512 with MI was replicated in the Korean population.

CONCLUSION

BTN2A1 may be a susceptibility gene for MI in Japanese individuals.

Gatekeepers of intestinal inflammation

The intestine is subjected to a barrage of insults from food, bacterial flora, and pathogens. Despite this constant antigenic challenge, the mucosal tissues lining the intestinal tract remain largely under control. The mechanisms regulating the homeostatic balance in the gut have been investigated for many years by many groups, but the precise nature of the regulatory control remains elusive. In this review, we provide an overview of pathways proposed to be involved in dampening the inflammatory response and maintaining the homeostatic balance in the intestine, and how these pathways may be disrupted in ulcerative colitis and Crohn's disease.

Gene hunting of the Genetic Analysis Workshop 16 rheumatoid arthritis data using rough set theory

We propose to use the rough set theory to identify genes affecting rheumatoid arthritis risk from the data collected by the North American Rheumatoid Arthritis Consortium. For each gene, we employ generalized dynamic reducts in the rough set theory to select a subset of single-nucleotide polymorphisms (SNPs) to represent the genetic information from this gene. We then group the study subjects into different clusters based on their genotype similarity at the selected markers. Statistical association between disease status and cluster membership is then studied to identify genes associated with rheumatoid arthritis. Based on our proposed approach, we are able to identify a number of statistically significant genes associated with rheumatoid arthritis. Aside from genes on chromosome 6, our identified genes include known disease-associated genes such as PTPN22 and TRAF1. In addition, our list contains other biologically plausible genes, such as ADAM15 and AGPAT2. Our findings suggest that ADAM15 and AGPAT2 may contribute to a genetic predisposition through abnormal angiogenesis and adipose tissue.

Detecting significant single-nucleotide polymorphisms in a rheumatoid arthritis study using random forests

Random forest is an efficient approach for investigating not only the effects of individual markers on a trait but also the effect of the interactions among the markers in genetic association studies. This approach is especially appealing for the analysis of genome-wide data, such as those obtained from gene expression/single-nucleotide polymorphism (SNP) array experiments in which the number of candidate genes/SNPs is vast. We applied this approach to the Genetic Analysis Workshop 16 Problem 1 data to identify SNPs that contribute to rheumatoid arthritis. The random forest computed a raw importance score for each SNP marker, where higher importance score suggests higher level of association between the marker and the trait. The significance level of the association was determined empirically by repeatedly reapplying the random forest on randomly generated data under the null hypothesis that no association exists between the markers and the trait. Using random forest, we were able to identify 228 significant SNPs (at the genome-wide significant level of 0.05) across the whole genome, over two-thirds of which are located on chromosome 6, especially clustered in the region of 6p21 containing the human leukocyte antigen (HLA) genes, such as gene HLA-DRB1 and HLA-DRA. Further analysis of this region indicates a strong association to the rheumatoid arthritis status.

Conditional analysis of the major histocompatibility complex in rheumatoid arthritis

We performed a whole-genome association study of rheumatoid arthritis susceptibility using Illumina 550k single-nucleotide polymorphism (SNP) genotypes of 868 cases and 1194 controls from the North American Rheumatoid Arthritis Consortium (NARAC). Structured association analysis with adjustment for potential population stratification yielded 200 SNPs with p < 1 x 10-8 for association with RA, all of which were on chromosome 6 in a 2.7-Mb region of the major histocompatibility complex (MHC). Given the extensive linkage equilibrium in the region and known risk of HLA-DRB1 alleles, we then applied conditional analyses to ascertain independent signals for RA susceptibility among these 200 candidate SNPs. Conditional analyses incorporating risk categories of the HLA-DRB1 "shared epitope" revealed three SNPs having independent associations with RA (conditional p < 0.001). This supports the presence of significant effects on RA susceptibility in the MHC in addition to the shared epitope.

Genome-wide association study of rheumatoid arthritis by a score test based on wavelet transformation

Background

We have conducted a genome-wide association study on the Genetic Analysis Workshop (GAW) 16 rheumatoid arthritis data using a multilocus score test based on wavelet transform proposed recently by the authors. The wavelet-based test automatically adjusts for the amount of noise suppressed from the data. The power of the test is also increased by using the genetic information contained in the spatial ordering of single-nucleotide polymorphisms on a chromosome.

Results

After adjusting for the effect of population stratification, the test identified some previously discovered rheumatoid arthritis susceptibility loci (HLA-DRB1 and rs3761847) as well as some loci (rs2076530 and rs3130340) known to have association with sarcoidosis and bone mineral density. It was previously reported that patients with rheumatoid arthritis have elevated prevalence of sarcoidosis and have reduced bone mass.

Conclusion

This new test provides a useful tool in genome-wide association studies.

Case-control genome-wide association study of rheumatoid arthritis from Genetic Analysis Workshop 16 using penalized orthogonal-components regression-linear discriminant analysis

Currently, genome-wide association studies (GWAS) are conducted by collecting a massive number of SNPs (i.e., large p) for a relatively small number of individuals (i.e., small n) and associations are made between clinical phenotypes and genetic variation one single-nucleotide polymorphism (SNP) at a time. Univariate association approaches like this ignore the linkage disequilibrium between SNPs in regions of low recombination. This results in a low reliability of candidate gene identification. Here we propose to improve the case-control GWAS approach by implementing linear discriminant analysis (LDA) through a penalized orthogonal-components regression (POCRE), a newly developed variable selection method for large p small n data. The proposed POCRE-LDA method was applied to the Genetic Analysis Workshop 16 case-control data for rheumatoid arthritis (RA). In addition to the two regions on chromosomes 6 and 9 previously associated with RA by GWAS, we identified SNPs on chromosomes 10 and 18 as potential candidates for further investigation.

Autoimmune disease classification by inverse association with SNP alleles

With multiple genome-wide association studies (GWAS) performed across autoimmune diseases, there is a great opportunity to study the homogeneity of genetic architectures across autoimmune disease. Previous approaches have been limited in the scope of their analysis and have failed to properly incorporate the direction of allele-specific disease associations for SNPs. In this work, we refine the notion of a genetic variation profile for a given disease to capture strength of association with multiple SNPs in an allele-specific fashion. We apply this method to compare genetic variation profiles of six autoimmune diseases: multiple sclerosis (MS), ankylosing spondylitis (AS), autoimmune thyroid disease (ATD), rheumatoid arthritis (RA), Crohn's disease (CD), and type 1 diabetes (T1D), as well as five non-autoimmune diseases. We quantify pair-wise relationships between these diseases and find two broad clusters of autoimmune disease where SNPs that make an individual susceptible to one class of autoimmune disease also protect from diseases in the other autoimmune class. We find that RA and AS form one such class, and MS and ATD another. We identify specific SNPs and genes with opposite risk profiles for these two classes. We furthermore explore individual SNPs that play an important role in defining similarities and differences between disease pairs. We present a novel, systematic, cross-platform approach to identify allele-specific relationships between disease pairs based on genetic variation as well as the individual SNPs which drive the relationships. While recognizing similarities between diseases might lead to identifying novel treatment options, detecting differences between diseases previously thought to be similar may point to key novel disease-specific genes and pathways.

Yin-Yang of costimulation: crucial controls of immune tolerance and function

SUMMARY

In addition to signals from the T-cell receptor complex, it has been recognized for many years that a 'second' signal, most notably from CD28, is also important in T-cell activation. In the recent years, many new members of CD28 family as well as the molecules that share structural homology to CD28 ligands CD80 and CD86 have been discovered. Interestingly, some of these proteins function to dampen T-cell activation and regulate the induction of T-cell tolerance. Therefore, positive and negative costimulation are the two sides of the coin to fine tune T-cell receptor signaling to determine the outcome of T-cell receptor engagement-tolerance versus function.

Genetic impact of a butyrophilin-like 2 (BTNL2) gene variation on specific IgE responsiveness to Dermatophagoides farinae (Der f) in Japanese

BACKGROUND

Dermatophagoides farinae (Der f) is one of the most frequently implicated allergens in several allergic diseases. Several genome-wide screens have identified a linkage between chromosome 6p21 and mite-specific IgE responsiveness. Butyrophilin-like 2 (BTNL2) is a member of the immunoglobulin superfamily and, on the basis of its homology to B7-1, has been implicated as a costimulatory molecule involved in T-cell activation. BTNL2 resides in the HLA region on chromosome 6p21, and significant associations between BTNL2 gene polymorphisms and several inflammatory diseases have been reported.

OBJECTIVE

The aim of this study was to examine whether BTNL2 gene polymorphisms are associated with specific IgE responses to Der f.

METHODS

Three single nucleotide polymorphisms (SNPs), including 2 coding SNPs and 1 intron SNP, were studied. One of the coding SNPs was the rs2076530 A > G, which has a functional consequence. A total of 863 unrelated Japanese subjects (447 positive and 416 negative for IgE to Der f) were recruited for a case-control study.

RESULTS

Controlling for gender, age, smoking, and the presence of asthma, multiple logistic regression analyses showed that homozygosity of the rs2076530 A allele, which has been reported to be a risk allele for sarcoidosis, was associated with a risk of sensitization towards Der f (Odds ratio; 1.55, p = 0.0060).

CONCLUSIONS

Although an association which may be due to the linkage disequilibrium with other genes in 6p21 needs to be ruled out, the present findings suggest that the BTNL2 gene might be one of the candidate genes that is responsible for the pathogenesis of Der f-specific IgE responsiveness.

Three microsatellites from the T1DGC MHC data set show highly significant association with type 1 diabetes, independent of the HLA-DRB1, -DQA1 and -DQB1 genes

AIM

The aim of this study was to test the microsatellites in the Type 1 Diabetes Genetics Consortium major histocompatibility complex (MHC) data set for association with type 1 diabetes (T1D) independent of the HLA-DRB1, -DQA1 and -DQB1 genes.

METHODS

The data set was edited to contain only one affected child per family, and broad ethnic subgroups were defined. Genotypes for HLA-DRB1, -DQA1 and -DQB1 were replaced by a haplotype code spanning all three loci, with phase inferred based on common haplotypes. The final data set contained 8190 samples in 2301 families, 59 microsatellites and the DRB1-DQA1-DQB1 haplotype code. Statistical analyses consisted of conditional logistic regression and haplotype estimations and linkage disequilibrium calculations.

RESULTS

The data set was screened using a main effects test approach adjusted for DRB1-DQA1-DQB1, and significant results tested for validity. After these procedures, four markers remained significant at the Bonferroni-corrected threshold: D6S2773 (p = 0.00014), DG6S185 (p = 0.00015), DG6S398 (p = 0.00043) and D6S2998 (p = 0.00015). These results were supported by allelic tests conditioned on DRB1-DQA1-DQB1 haplotypes, except for DG6S185, which may contain artefacts.

CONCLUSIONS

We have identified three microsatellites that mark additional risk factors for T1D at highly significant levels in the MHC. Further analyses are needed to establish the relationship with other possible genetic determinants in this region.

Protein tyrosine phosphatases in autoimmunity

Protein tyrosine phosphatases (PTPs) are important regulators of many cellular functions and a growing number of PTPs have been implicated in human disease conditions, such as developmental defects, neoplastic disorders, and immunodeficiency. Here, we review the involvement of PTPs in human autoimmunity. The leading examples include the allelic variant of the lymphoid tyrosine phosphatase (PTPN22), which is associated with multiple autoimmune diseases, and mutations that affect the exon-intron splicing of CD45 (PTPRC). We also find it likely that additional PTPs are involved in susceptibility to autoimmune and inflammatory diseases. Finally, we discuss the possibility that PTPs regulating the immune system may serve as therapeutic targets.

Analysis of BTNL2 genetic polymorphisms in British and Dutch patients with sarcoidosis

Sarcoidosis is a heterogeneous disorder, both phenotypically and genetically. Two independent studies have recently shown that a functional polymorphism within butyrophilin-like 2 (BTNL2) gene predisposes to sarcoidosis independently of the human leukocyte antigen (HLA)-DRB1 alleles. However, in both studies, data analysis was not stratified by Löfgren's syndrome, a clinically and genetically distinct sarcoidosis subset. BTNL2, potentially encoding an immune coreceptor, is adjacent and in linkage disequilibrium (LD) with HLA-DRB1. We investigated six BTNL2 variants, including the functional rs2076530 (G > A), as well as HLA-DRB1 alleles, by sequence-specific primers-polymerase chain reaction, in 288 patients and 446 controls from two European countries. In the patient group as a whole, the HLA-DRB1*14 [odds ratio (OR) = 3.1, P(c) = 0.0003], DRB1*12 (OR = 2.5, P(c) = 0.003), and BTNL2 rs2076530 A allele (OR = 1.49, P(c) = 0.002) were all associated with disease susceptibility. However, after exclusion of patients presenting with Löfgren's syndrome and after adjusting for HLA-DRB1 alleles, the association between BTNL2 rs2076530 A and disease disappeared (P = 0.23). By contrast, both HLA-DRB1*14 and DRB1*12 remained strongly significant (OR = 3.60, P < 0.0001 and OR = 3.03, P = 0.003, respectively). BTNL2 haplotype 4, tagged by the rs2076530 G allele, also remained associated with non-Löfgren sarcoidosis after adjusting for HLA-DRB1 alleles (OR 0.37, P = 0.016). In summary, HLA-DRB1*14, DRB1*12, and BTNL2 haplotype 4--but not rs2076530 A--are associated with non-Löfgren sarcoidosis. However, the tight LD across the HLA complex makes it difficult to identify the precise location of the susceptibility locus/i. Larger sample sets from different ethnic groups, finer mapping, and more robust LD analyses across the HLA region are needed.

BTNL2 allele associations with chronic beryllium disease in HLA-DPB1*Glu69-negative individuals

Butyrophilin-like 2 (BTNL2) polymorphisms have been associated with sarcoidosis. We hypothesized that BTNL2 variants might confer a human leukocyte antigen (HLA)-independent risk for chronic beryllium disease (CBD), a granulomatous lung disease with similar clinical, radiological, and pathological features to sarcoidosis. Genomic DNA was obtained from CBD (n= 88), beryllium sensitized (BeS, n= 86), and beryllium exposed nondiseased control subjects (Be-exp, n= 173). Six BTNL2 polymorphisms, HLA-DPB1, DRB1, and DQB1 alleles were determined by sequence-specific primer-PCR. All BTNL2 polymorphisms were in Hardy-Weinberg equilibrium. No significant differences were found between BTNL2 polymorphisms or haplotypes and CBD, BeS, or Be-exp. In HLA-DPB1*Glu69-negative subjects (n= 10 CBD, n= 13 BeS, n= 102 Be-exp), DRB1*13 and BTNL2 rs3117099TT homozygosity were increased in CBD (70% and 40%, respectively) vs Be-exp (16%, P= 0.001 and 2.9%, P= 0.001, respectively). The BTNL2 rs3117099T-HLA-DRB1*13 combination was significantly increased in CBD (50%) compared with Be-exp (6.9%, P= 0.001). In conclusion, both DRB1*13 and rs3117099TT homozygosity are associated with CBD in *Glu69-negative subjects, while DPB1*Glu69 is associated with CBD and BeS compared with Be-exp. As a result of the small sample size and strong linkage disequilibrium between DRB1*13-DQB1*0603/4/9 and the BTNL2 rs3117099T allele, it is difficult to assess the primary association in DPB1*Glu69-negative CBD cases.