Exome sequencing identifies novel rheumatoid arthritis-susceptible variants in the BTNL2

Low-frequency and rare genetic variants associated with rheumatoid arthritis risk

Rheumatoid arthritis (RA) has an estimated heritability of nearly 50%, which is particularly high in seropositive RA. HLA alleles account for a large proportion of this heritability, in addition to many common single-nucleotide polymorphisms with smaller individual effects. Low-frequency and rare variants, such as those captured by next-generation sequencing, can also have a large role in heritability in some individuals. Rare variant discovery has informed the development of drugs such as inhibitors of PCSK9 and Janus kinases. Some 34 low-frequency and rare variants are currently associated with RA risk. One variant (19:10352442G>C in TYK2) was identified in five separate studies, and might therefore represent a promising therapeutic target. Following a set of best practices in future studies, including studying diverse populations, using large sample sizes, validating RA and serostatus, replicating findings, adjusting for other variants and performing functional assessment, could help to ensure the relevance of identified variants. Exciting opportunities are now on the horizon for genetics in RA, including larger datasets and consortia, whole-genome sequencing and direct applications of findings in the management, and especially treatment, of RA.

The B7:CD28 family and friends: Unraveling coinhibitory interactions

Immune responses must be tightly regulated to ensure both optimal protective immunity and tolerance. Costimulatory pathways within the B7:CD28 family provide essential signals for optimal T cell activation and clonal expansion. They provide crucial inhibitory signals that maintain immune homeostasis, control resolution of inflammation, regulate host defense, and promote tolerance to prevent autoimmunity. Tumors and chronic pathogens can exploit these pathways to evade eradication by the immune system. Advances in understanding B7:CD28 pathways have ushered in a new era of immunotherapy with effective drugs to treat cancer, autoimmune diseases, infectious diseases, and transplant rejection. Here, we discuss current understanding of the mechanisms underlying the coinhibitory functions of CTLA-4, PD-1, PD-L1:B7-1 and PD-L2:RGMb interactions and less studied B7 family members, including HHLA2, VISTA, BTNL2, and BTN3A1, as well as their overlapping and unique roles in regulating immune responses, and the therapeutic potential of these insights.

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.

Application of NGS Technology in Understanding the Pathology of Autoimmune Diseases

NGS technologies have transformed clinical diagnostics and broadly used from neonatal emergencies to adult conditions where the diagnosis cannot be made based on clinical symptoms. Autoimmune diseases reveal complicate molecular background and traditional methods could not fully capture them. Certainly, NGS technologies meet the needs of modern exploratory research, diagnostic and pharmacotherapy. Therefore, the main purpose of this review was to briefly present the application of NGS technology used in recent years in the understanding of autoimmune diseases paying particular attention to autoimmune connective tissue diseases. The main issues are presented in four parts: (a) panels, whole-genome and -exome sequencing (WGS and WES) in diagnostic, (b) Human leukocyte antigens (HLA) as a diagnostic tool, (c) RNAseq, (d) microRNA and (f) microbiome. Although all these areas of research are extensive, it seems that epigenetic impact on the development of systemic autoimmune diseases will set trends for future studies on this area.

kTWAS: integrating kernel machine with transcriptome-wide association studies improves statistical power and reveals novel genes

The power of genotype-phenotype association mapping studies increases greatly when contributions from multiple variants in a focal region are meaningfully aggregated. Currently, there are two popular categories of variant aggregation methods. Transcriptome-wide association studies (TWAS) represent a set of emerging methods that select variants based on their effect on gene expressions, providing pretrained linear combinations of variants for downstream association mapping. In contrast to this, kernel methods such as sequence kernel association test (SKAT) model genotypic and phenotypic variance use various kernel functions that capture genetic similarity between subjects, allowing nonlinear effects to be included. From the perspective of machine learning, these two methods cover two complementary aspects of feature engineering: feature selection/pruning and feature aggregation. Thus far, no thorough comparison has been made between these categories, and no methods exist which incorporate the advantages of TWAS- and kernel-based methods. In this work, we developed a novel method called kernel-based TWAS (kTWAS) that applies TWAS-like feature selection to a SKAT-like kernel association test, combining the strengths of both approaches. Through extensive simulations, we demonstrate that kTWAS has higher power than TWAS and multiple SKAT-based protocols, and we identify novel disease-associated genes in Wellcome Trust Case Control Consortium genotyping array data and MSSNG (Autism) sequence data. The source code for kTWAS and our simulations are available in our GitHub repository (https://github.com/theLongLab/kTWAS).

Genetic Variation and Cardiovascular Risk Factors: A Cohort Study on Migrants from the Former Soviet Union and a Native German Population

Resettlers are a large migrant group of more than 2 million people in Germany who migrated mainly from the former Soviet Union to Germany after 1989. We sought to compare the distribution of the major risk factors for cardiovascular disease (CVD) and to investigate the overall genetic differences in a study population which consisted of resettlers and native (autochthone) Germans. This was a joint analysis of two cohort studies which were performed in the region of Augsburg, Bavaria, Germany, with 3363 native Germans and 363 resettlers. Data from questionnaires and physical examinations were used to compare the risk factors for cardiovascular diseases between the resettlers and native Germans. A population-based genome-wide association analysis was performed in order to identify the genetic differences between the two groups. The distribution of the major risk factors for CVD differed between the two groups. The resettlers lead a less active lifestyle. While female resettlers smoked less than their German counterparts, the men showed similar smoking behavior. SNPs from three genes (BTNL2, DGKB, TGFBR3) indicated a difference in the two populations. In other studies, these genes have been shown to be associated with CVD, rheumatoid arthritis and osteoporosis, respectively.

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.

Germline genetic patterns underlying familial rheumatoid arthritis, systemic lupus erythematosus and primary Sjögren's syndrome highlight T cell-initiated autoimmunity

OBJECTIVES

Familial aggregation of primary Sjögren's syndrome (pSS), systemic lupus erythematosus (SLE), rheumatoid arthritis (RA) and co-aggregation of these autoimmune diseases (ADs) (also called familial autoimmunity) is well recognised. However, the genetic predisposition variants that explain this clustering remains poorly defined.

METHODS

We used whole-exome sequencing on 31 families (9 pSS, 11 SLE, 6 RA and 5 mixed autoimmunity), followed by heterozygous filtering and cosegregation analysis of a family-focused approach to document rare variants predicted to be pathogenic by in silico analysis. Potential importance in immune-related processes, gene ontology, pathway enrichment and overlap analyses were performed to prioritise gene sets.

RESULTS

A range from 1 to 50 rare possible pathogenic variants, including 39 variants in immune-related genes across SLE, RA and pSS families, were identified. Among this gene set, regulation of T cell activation (p=4.06×10-7) and T cell receptor (TCR) signalling pathway (p=1.73×10-6) were particularly concentrated, including PTPRC (CD45), LCK, LAT-SLP76 complex genes (THEMIS, LAT, ITK, TEC, TESPA1, PLCL1), DGKD, PRKD1, PAK2 and NFAT5, shared across 14 SLE, RA and pSS families. TCR-interactive genes P2RX7, LAG3, PTPN3 and LAX1 were also detected. Overlap analysis demonstrated that the antiviral immunity gene DUS2 variant cosegregated with SLE, RA and pSS phenotypes in an extended family, that variants in the TCR-pathway genes CD45, LCK and PRKD1 occurred independently in three mixed autoimmunity families, and that variants in CD36 and VWA8 occurred in both RA-pSS and SLE-pSS families.

CONCLUSIONS

Our preliminary results define common genetic characteristics linked to familial pSS, SLE and RA and highlight rare genetic variations in TCR signalling pathway genes which might provide innovative molecular targets for therapeutic interventions for those three ADs.

Primary Immunodeficiency with Severe Multi-Organ Immune Dysregulation

Polyglandular autoimmune syndrome type 1, also known as autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED), is a rare primary immunodeficiency disorder with multi-organ involvement. Besides for being predisposed to severe life-threatening infections, patients with APECED are also prone to organ impairment secondary to severe autoimmunity. As this is an autosomal recessive disorder, a biallelic mutation in the AIRE gene is responsible for APECED. The author presents a case of APECED with a single AIRE mutation. Whole exome sequencing identified a mutation in the BTNL2 gene that the author suggests may have contributed to the patient's presentation.

Phosphorylation-related SNPs influence lipid levels and rheumatoid arthritis risk by altering gene expression and plasma protein levels

Objectives

Phosphorylation-related single-nucleotide polymorphisms (phosSNPs) are missense SNPs that may influence protein phosphorylation. The aim of this study was to evaluate the effect of phosSNPs on lipid levels and RA.

Methods

We examined the association of phosSNPs with lipid levels and RA in large-scale genome-wide association studies (GWAS) and performed random sampling and fgwas analyses to determine whether the phosSNPs associated with lipid levels and RA were significantly enriched. Furthermore, we performed QTL analysis and Mendelian randomization analysis to obtain additional evidence to be associated with the identified phosSNPs and genes.

Results

We found 483 phosSNPs for lipid levels and 243 phosSNPs for RA in the GWAS loci (P &lt; 1.0 × 10−5). SNPs associated with high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, Total cholesterol (TC) and RA were significantly enriched with phosSNPs. Almost all of the identified phosSNPs showed expression quantitative trait loci (eQTL) effects. A total of 48 protein QTLs and 9 metabolite QTLs were found. The phosSNP rs3184504 (p.Trp262Arg) at SH2B3 was significantly associated with RA, SH2B3 expression level, and plasma levels of high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, TC, hypoxanthine and 80 proteins, including beta-2-microglobulin. SH2B3 was differentially expressed between RA cases and controls in peripheral blood mononuclear cells and synovial tissues. Mendelian randomization analysis showed that SH2B3 expression level was significantly associated with TC level and RA. Plasma beta-2-microglobulin level was causally associated with high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, TC levels and RA.

Conclusion

The findings suggested that phosSNPs may play important roles in lipid metabolism and the pathological mechanisms of RA. PhosSNPs may influence lipid levels and RA risk by altering gene expression and plasma protein levels.

Clustering of immune-mediated diseases in sarcoidosis

PURPOSE OF REVIEW

Sarcoidosis is an immune-mediated disease of unknown cause. Immune-mediated diseases appear to cluster in patients and in families. We review what is known on this topic for sarcoidosis, and what factors may underlie disease clustering.

RECENT FINDINGS

In populations of patients with sarcoidosis, relative risk estimates of Sjögren's syndrome, systemic lupus erythematosus, autoimmune hepatitis, ankylosing spondylitis, multiple sclerosis (MS), celiac disease, autoimmune thyroid disease, and ulcerative colitis, varied between 2.1 and 11.6. In relatives of patients with sarcoidosis, relative risk estimates varied between 1.3 and 5.8 for sarcoidosis, MS, celiac disease, type 1 diabetes, Graves' disease, rheumatoid arthritis, Crohn's disease, and ulcerative colitis. Shared risk loci in key immunological pathways provide evidence for a contribution to development of multiple diseases. Identical changes in the immune status, epigenetic alterations, and environmental triggers have been detected in several diseases, and drug-induced disease is likely responsible for a small portion of co-occurring disease.

SUMMARY

Clustering of sarcoidosis and other immune-mediated diseases in patients and in their relatives occurs for sarcoidosis, MS, celiac disease, Graves' disease, and ulcerative colitis. Further research is needed to substantiate causal links and risk estimates in patients and their relatives.

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.

T-cell immunology in sarcoidosis: Disruption of a delicate balance between helper and regulatory T-cells

PURPOSE OF REVIEW

Although the aetiology of sarcoidosis is not yet completely understood, immunological changes within the T-cell compartment are characteristic for an exaggerated antigen-driven immune response. In this review, we describe the most recent findings on T-cell subset responses and regulation in sarcoidosis. We discuss how future immunological research can advance the field to unravel pathobiological mechanisms of this intriguingly complex disease.

RECENT FINDINGS

Research into the field of T-cell plasticity has recently challenged the long-held T helper type 1 (Th1) paradigm in sarcoidosis and striking parallels with autoimmune disorders and common variable immunodeficiency were recognized. For instance, it was demonstrated that Th17.1-cells rather than Th1-cells are responsible for the exaggerated IFN-γ production in pulmonary sarcoidosis. Furthermore, impaired regulatory T-cell function and alterations within the expression of co-inhibitory receptors that control T-cell responses, such as PD-1, CTLA-4 and BTNL2, raise new questions regarding T-cell regulation in pulmonary sarcoidosis.

SUMMARY

It becomes increasingly clear that Th17(.1)-cells and regulatory T-cells are key players in sarcoidosis T-cell immunology. New findings on plasticity and co-inhibitory receptor expression by these subsets help build a more comprehensive model for T-cell regulation in sarcoidosis and will finally shed light on the potential of new treatment modalities.

Differentially Methylated DNA Regions in Monozygotic Twin Pairs Discordant for Rheumatoid Arthritis: An Epigenome-Wide Study

Objectives

In an explorative epigenome-wide association study (EWAS) to search for gene independent, differentially methylated DNA positions and regions (DMRs) associated with rheumatoid arthritis (RA) by studying monozygotic (MZ) twin pairs discordant for RA.

Methods

Genomic DNA was isolated from whole blood samples from 28 MZ twin pairs discordant for RA. DNA methylation was measured using the HumanMethylation450 BeadChips. Smoking, anti-cyclic citrullinated peptide antibodies, and immunosuppressive treatment were included as covariates. Pathway analysis was performed using GREAT.

Results

Smoking was significantly associated with hypomethylation of a DMR overlapping the promoter region of the RNF5 and the AGPAT1, which are implicated in inflammation and autoimmunity, whereas DMARD treatment induced hypermethylation of the same region. Additionally, the promotor region of both S100A6 and EFCAB4B were hypomethylated, and both genes have previously been associated with RA. We replicated several candidate genes identified in a previous EWAS in treatment-naïve RA singletons. Gene-set analysis indicated the involvement of immunologic signatures and cancer-related pathways in RA.

Conclusion

We identified several differentially methylated regions associated with RA, which may represent environmental effects or consequences of the disease and plausible biological pathways pertinent to the pathogenesis of RA.

The ontogeny of Butyrophilin-like (Btnl) 1 and Btnl6 in murine small intestine

Murine Butyrophilin-like (Btnl) 1 and Btnl6 are primarily restricted to intestinal epithelium where they regulate the function of intraepithelial T lymphocytes. We recently demonstrated that Btnl1 and Btnl6 can form an intra-family heterocomplex and that the Btnl1-Btnl6 complex selectively expands Vγ7Vδ4 TCR IELs. To define the regulation of Btnl expression in the small intestine during ontogeny we examined the presence of Btnl1 and Btnl6 in the small bowel of newborn to 4-week-old mice. Although RNA expression of Btnl1 and Btnl6 was detected in the small intestine at day 0, Btnl1 and Btnl6 protein expression was substantially delayed and was not detectable in the intestinal epithelium until the mice reached 2–3 weeks of age. The markedly elevated Btnl protein level at week 3 coincided with a significant increase of γδ TCR IELs, particularly those bearing the Vγ7Vδ4 receptor. This was not dependent on gut microbial colonization as mice housed in germ-free conditions had normal Btnl protein levels. Taken together, our data show that the expression of Btnl1 and Btnl6 is delayed in the murine neonatal gut and that the appearance of the Btnl1 and Btnl6 proteins in the intestinal mucosa associates with the expansion of Vγ7Vδ4 TCR IELs.

Contribution of a Non-classical HLA Gene, HLA-DOA, to the Risk of Rheumatoid Arthritis

Despite the progress in human leukocyte antigen (HLA) causal variant mapping, independent localization of major histocompatibility complex (MHC) risk from classical HLA genes is challenging. Here, we conducted a large-scale MHC fine-mapping analysis of rheumatoid arthritis (RA) in a Japanese population (6,244 RA cases and 23,731 controls) population by using HLA imputation, followed by a multi-ethnic validation study including east Asian and European populations (n = 7,097 and 23,149, respectively). Our study identified an independent risk of a synonymous mutation at HLA-DOA, a non-classical HLA gene, on anti-citrullinated protein autoantibody (ACPA)-positive RA risk (p = 1.4 × 10(-9)), which demonstrated a cis-expression quantitative trait loci (cis-eQTL) effect on HLA-DOA expression. Trans-ethnic comparison revealed different linkage disequilibrium (LD) patterns in HLA-DOA and HLA-DRB1, explaining the observed HLA-DOA variant risk heterogeneity among ethnicities, which was most evident in the Japanese population. Although previous HLA fine-mapping studies have identified amino acid polymorphisms of the classical HLA genes as driving genetic susceptibility to disease, our study additionally identifies the dosage contribution of a non-classical HLA gene to disease etiology. Our study contributes to the understanding of HLA immunology in human diseases and suggests the value of incorporating additional ancestry in MHC fine-mapping.

Deep resequencing of 131 Crohn's disease associated genes in pooled DNA confirmed three reported variants and identified eight novel variants

OBJECTIVE

Genome wide association studies (GWAS) and meta-analyses for Crohn's disease (CD) have not fully explained the heritability of CD, suggesting that additional loci are yet to be found and that the known loci may contain high effect rare risk variants that have thus far gone undetected by GWAS. While the cost of deep sequencing remains too high to analyse many samples, targeted sequencing of pooled DNA samples allows the efficient and cost effective capture of all variations in a target region.

DESIGN

We performed pooled sequencing in 500 Korean CD cases and 1000 controls to evaluate the coding exon and 5' and 3' untranslated regions of 131 CD associated genes. The identified genetic variants were validated using genotyping in an independent set of 500 CD cases and 1000 controls.

RESULTS

Pooled sequencing identified 30 common/low single nucleotide variants (SNVs) in 12 genes and 3 rare SNVs in 3 genes. Our results confirmed a significant association of CD with the following previously reported risk loci: rs3810936 in TNFSF15 (OR=1.83, p<2.2×10(-16)), rs76418789 in IL23R (OR=0.47, p=1.14×10(-8)) and rs2241880 in ATG16L1 (OR=1.30, p=5.28×10(-6)). In addition, novel loci were identified in TNFSF8 (rs3181374, OR=1.53, p=1.03×10(-14)), BTNL2 (rs28362680, OR=1.47, p=9.67×10(-11)), HLA-DQA2 (rs3208181, OR=1.36, p=4.66×10(-6)), STAT3 (rs1053004, OR=1.29, p=2.07×10(-5)), NFKBIA (rs2273650, OR=0.80, p=3.93×10(-4)), NKX2-3 (rs888208, OR=0.82, p=6.37×10(-4)) and DNAH12 (rs4462937, OR=1.13, p=3.17×10(-2)). A novel rare SNV, rs200735402 in CARD9, was shown to have a protective effect (OR=0.09, p=5.28×10(-5)).

CONCLUSIONS

Our deep resequencing of 131 CD associated genes confirmed 3 reported risk loci and identified 8 novel risk loci for CD in Koreans, providing new insights into the genetic architecture of CD.

Altered expression of Butyrophilin (BTN) and BTN-like (BTNL) genes in intestinal inflammation and colon cancer

Several Butyrophilin (BTN) and Btn‐like (BTNL) molecules control T lymphocyte responses, and are genetically associated with inflammatory disorders and cancer. In this study, we present a comprehensive expression analysis of human and murine BTN and BTNL genes in conditions associated with intestinal inflammation and cancer. Using real‐time PCR, expression of human BTN and BTNL genes was analyzed in samples from patients with ulcerative colitis, irritable bowel syndrome, and colon tumors. Expression of murine Btn and Btnl genes was examined in mouse models of spontaneous colitis (Muc2^−/−) and intestinal tumorigenesis (Apc^Min/+). Our analysis indicates a strong association of several of the human genes with ulcerative colitis and colon cancer; while especially BTN1A1, BTN2A2, BTN3A3, and BTNL8 were significantly altered in inflammation, colonic tumors exhibited significantly decreased levels of BTNL2, BTNL3, BTNL8, and BTNL9 as compared to unaffected tissue. Colonic inflammation in Muc2^−/− mice significantly down‐regulated the expression of particularly Btnl1, Btnl4, and Btnl6 mRNA, and intestinal polyps derived from Apc^Min/+ mice displayed altered levels of Btn1a1, Btn2a2, and Btnl1 transcripts. Thus, our data present an association of BTN and BTNL genes with intestinal inflammation and cancer and represent a valuable resource for further studies of this gene family.

Analysis of genomic copy number variation in equine recurrent airway obstruction (heaves)

We explored the involvement of genomic copy number variants (CNVs) in susceptibility to recurrent airway obstruction (RAO), or heaves-an asthmalike inflammatory disease in horses. Analysis of 16 RAO-susceptible (cases) and six RAO-resistant (control) horses on a custom-made whole-genome 400K equine tiling array identified 245 CNV regions (CNVRs), 197 previously known and 48 new, distributed on all horse autosomes and the X chromosome. Among the new CNVRs, 30 were exclusively found in RAO cases and were further analyzed by quantitative PCR, including additional cases and controls. Suggestive association (P = 0.03; corrected P = 0.06) was found between RAO and a loss on chromosome 5 involving NME7, a gene necessary for ciliary functions in lungs and involved in primary ciliary dyskinesia in humans. The CNVR could be a potential marker for RAO susceptibility but needs further study in additional RAO cohorts. Other CNVRs were not associated with RAO, although several involved genes of interest, such as SPI2/SERPINA1 from the serpin gene family, which are associated with chronic obstructive pulmonary disease and asthma in humans. The SPI2/SERPINA1 CNVR showed striking variation among horses, but it was not significantly different between RAO cases and controls. The findings provide baseline information on the relationship between CNVs and RAO susceptibility. Discovery of new CNVs and the use of a larger population of RAO-affected and control horses are needed to shed more light on their significance in modulating this complex and heterogeneous disease.

Murine Butyrophilin-Like 1 and Btnl6 Form Heteromeric Complexes in Small Intestinal Epithelial Cells and Promote Proliferation of Local T Lymphocytes

To date, few molecular conduits mediating the cross-talk between intestinal epithelial cells and intraepithelial lymphocytes (IELs) have been described. We recently showed that butyrophilin-like (Btnl) 1 can attenuate the epithelial response to activated IELs, resulting in reduced production of proinflammatory mediators, such as IL-6 and CXCL1. We here report that like Btnl1, murine Btnl6 expression is primarily confined to the intestinal epithelium. Although Btnl1 can exist in a cell surface-expressed homomeric form, we found that it additionally forms heteromeric complexes with Btnl6, and that the engagement of Btnl1 is a prerequisite for surface expression of Btnl6 on intestinal epithelial cells. In an IEL-epithelial cell coculture system, enforced epithelial cell expression of Btnl1 significantly enhanced the proliferation of IELs in the absence of exogenous activation. The effect on proliferation was dependent on the presence of IL-2 or IL-15 and restricted to IELs upregulating CD25. In the γδ T-cell subset, the Btnl1-Btnl6 complex, but not Btnl1, specifically elevated the proliferation of IELs bearing the Vγ7Vδ4 receptor. Thus, our results show that murine epithelial cell-specific Btnl proteins can form intrafamily heterocomplexes and suggest that the interaction between Btnl proteins and IELs regulates the expansion of IELs in the intestinal mucosa.

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.

Aggregation of rare/low-frequency variants of the mitochondria respiratory chain-related proteins in rheumatoid arthritis patients

Exome sequencings were conducted using 59 patients having rheumatoid arthritis (RA) and 93 controls. After stepwise filtering, 107 genes showed less than 0.05 of P-values by gene-burden tests. Among 107 genes, NDUFA7 which is a subunit of the complex I in the mitochondrial respiratory chain was selected for further analysis based on previous reports. A case-control study was performed on the three single-nucleotide variants (SNVs) of NDUFA7 with 432 cases and 432 controls. An association was observed between NDUFA7 and RA with severe erosive arthritis. These results together with previous reports suggested the involvement of reactive oxygen species (ROS) in the pathogenesis of RA. In the next step, four SNVs from three genes related to the mitochondrial respiratory chain were selected, which is a major source of ROS, and conducted a case-control study. An association was observed based on a pathway-burden test comprising NDUFA7, SDHAF2, SCO1 and ATP5O: P=1.56E-04, odds ratio=2.16, 95% confidence interval=1.43-3.28. Previous reports suggested the involvement of ROS in the pathogenesis of RA. The aggregation of SNVs in the mitochondria respiratory chain suggests the pivotal role of those SNVs in the pathogenesis of RA with severe erosive arthritis.

Pooled sequencing of 531 genes in inflammatory bowel disease identifies an associated rare variant in BTNL2 and implicates other immune related genes

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn’s disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10⁻¹⁰, OR = 2.3[95% CI = 1.75–3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (<1%) and low frequency (1–5%) variants in 3 additional genes showed suggestive association (p<0.005) with either an increased risk (ARIH2 c.338-6C>T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis.

Crohn’s disease and ulcerative colitis are two forms of inflammatory bowel disease which cause chronic inflammation of the gastrointestinal tract. Common genetic variants in more than 160 regions of the human genome have been associated with an altered risk of these disorders, but leave much of the estimated genetic contribution to disease risk unexplained. We sought to establish whether rare genetic variants which alter the structure or function of the proteins encoded by genes also contribute to disease susceptibility. We used high throughput DNA sequencing to screen over 500 genes for such variants in nearly 500 patients and controls, and validated interesting variants in about 10,000 patients and 7,000 controls. We detected association of a limited number of rare variants from coding regions with disease, suggesting that they do not account for a large proportion of genetic susceptibility. However, they highlight the involvement of genes of potential importance in the development of inflammatory bowel disease, including those involved in the activation of immune cells, the regulation of immune response genes, and the degradation of proteins in cells.

Clozapine-induced agranulocytosis is associated with rare HLA-DQB1 and HLA-B alleles

Clozapine is a particularly effective antipsychotic medication but its use is curtailed by the risk of clozapine-induced agranulocytosis/granulocytopenia (CIAG), a severe adverse drug reaction occurring in up to 1% of treated individuals. Identifying genetic risk factors for CIAG could enable safer and more widespread use of clozapine. Here we perform the largest and most comprehensive genetic study of CIAG to date by interrogating 163 cases using genome-wide genotyping and whole-exome sequencing. We find that two loci in the major histocompatibility complex are independently associated with CIAG: a single amino acid in HLA-DQB1 (126Q) (P=4.7 × 10(-14), odds ratio (OR)=0.19, 95% confidence interval (CI)=0.12-0.29) and an amino acid change in the extracellular binding pocket of HLA-B (158T) (P=6.4 × 10(-10), OR=3.3, 95% CI=2.3-4.9). These associations dovetail with the roles of these genes in immunogenetic phenotypes and adverse drug responses for other medications, and provide insight into the pathophysiology of CIAG.

BTNL2 gene SNPs as a contributing factor to sarcoidosis pathogenesis in a cohort of Greek patients

Sarcoidosis is a multisystemic granulomatous disease of unknown etiology that primarily affects adults between the ages of 20 and 40 years old. It is characterized by the activation of Th1 lymphocytes resulting in the production of inflammatory cytokines and the formation of noncaseating epithelioid cell granulomas in affected tissues. The lungs and lymphatic system are the ones most frequently affected. The disease usually presents spontaneous remission in the first two years and, in a few patients, the disease progresses to pulmonary fibrosis or other fatal complications depending on the affected organ. The pathogenesis of sarcoidosis is still not clearly defined, and is considered an interaction between the environment and risk alleles in many genes.

The present case control study consisted of 146 Greek patients with sarcoidosis and 90 healthy volunteers from the same ethnic group. The coding and neighboring intronic regions of the BTNL2 gene were sequenced and risk alleles were compared amongst the two groups. Thirty-seven different variants were detected from which 12 were synonymous substitutions and 25 non-synonymous. With the help of in silico tools (SIFT, PolyPhen, PROVEAN, PMut and EX_SKIP), 13 variants were classified as possible pathological risk variants including 4 novel ones. The most common risk variants contributing to phenotypic modulation of sarcoidosis were p.S360G and p.S334L, with the latter contributing to a more severe disease stage with extra-pulmonary manifestations such as skin granulomas and relapses being more common.

•

Thirteen pathological risk variants were discovered in BTNL2 gene in sarcoidosis patients.

•

The common risk variant p.S360G was found in cases (37.67%) and controls (30%).

•

The 2nd most common variant p.S334L was found in cases (16.43%) and controls (12.22%).

•

Six novel variants were detected with 4 deemed pathological.

•

Coinheritance of common and novel variants affected the final clinical phenotype.

Genetic analysis in young-age-of-onset Graves' disease reveals new susceptibility loci

CONTEXT

Genetic and environmental factors play an essential role in the pathogenesis of Graves' Disease (GD). Children with GD have less exposure time to environmental factors and therefore are believed to harbor stronger genetic susceptibility than adults.

OBJECTIVE

The aim of the study was to identify susceptibility loci that predispose to GD in patients with young-age-of-onset (YAO) GD.

SETTING AND DESIGN

One hundred six patients with YAO GD (onset <30 y) and 855 healthy subjects were studied. Cases and controls were genotyped using the Illumina Infinium Immunochip, designed to genotype 196,524 polymorphisms. Case control association analyses were performed using the PLINK computer package. Ingenuity Pathway Analysis program (QIAGEN) was used to carry out pathway analyses.

RESULTS

Immunochip genetic association analysis identified 30 single-nucleotide polymorphisms in several genes that were significantly associated with YAO GD, including major histocompatibility complex class I and class II genes, BTNL2, NOTCH4, TNFAIP3, and CXCR4. Candidate gene analysis revealed that most of the genes previously shown to be associated with adult-onset GD were also associated with YAO GD. Pathway analysis demonstrated that antigen presentation, T-helper cell differentiation, and B cell development were the major pathways contributing to the pathogenesis of YAO GD.

CONCLUSIONS

Genetic analysis identified novel susceptibility loci in YAO GD adding a new dimension to the understanding of GD etiology.

Diversity of extended HLA-DRB1 haplotypes in the Finnish population

The Major Histocompatibility Complex (MHC, 6p21) codes for traditional HLA and other host response related genes. The polymorphic HLA-DRB1 gene in MHC Class II has been associated with several complex diseases. In this study we focus on MHC haplotype structures in the Finnish population. We explore the variability of extended HLA-DRB1 haplotypes in relation to the other traditional HLA genes and a selected group of MHC class III genes. A total of 150 healthy Finnish individuals were included in the study. Subjects were genotyped for HLA alleles (HLA-A, -B, -DRB1, -DQB1, and -DPB1). The polymorphism of TNF, LTA, C4, BTNL2 and HLA-DRA genes was studied with 74 SNPs (single nucleotide polymorphism). The C4A and C4B gene copy numbers and a 2-bp silencing insertion at exon 29 in C4A gene were analysed with quantitative genomic realtime-PCR. The allele frequencies for each locus were calculated and haplotypes were constructed using both the traditional HLA alleles and SNP blocks. The most frequent Finnish A∼B∼DR -haplotype, uncommon in elsewhere in Europe, was A*03∼B*35∼DRB1*01∶01. The second most common haplotype was a common European ancestral haplotype AH 8.1 (A*01∼B*08∼DRB1*03∶01). Extended haplotypes containing HLA-B, TNF block, C4 and HLA-DPB1 strongly increased the number of HLA-DRB1 haplotypes showing variability in the extended HLA-DRB1 haplotype structures. On the contrary, BTNL2 block and HLA-DQB1 were more conserved showing linkage with the HLA-DRB1 alleles. We show that the use of HLA-DRB1 haplotypes rather than single HLA-DRB1 alleles is advantageous when studying the polymorphisms and LD patters of the MHC region. For disease association studies the HLA-DRB1 haplotypes with various MHC markers allows us to cluster haplotypes with functionally important gene variants such as C4 deficiency and cytokines TNF and LTA, and provides hypotheses for further assessment. Our study corroborates the importance of studying population-specific MHC haplotypes.

Germline missense variants in the BTNL2 gene are associated with prostate cancer susceptibility

BACKGROUND

Rare, inherited mutations account for 5% to 10% of all prostate cancer cases. However, to date, few causative mutations have been identified.

METHODS

To identify rare mutations for prostate cancer, we conducted whole-exome sequencing (WES) in multiple kindreds (n = 91) from 19 hereditary prostate cancer (HPC) families characterized by aggressive or early-onset phenotypes. Candidate variants (n = 130) identified through family- and bioinformatics-based filtering of WES data were then genotyped in an independent set of 270 HPC families (n = 819 prostate cancer cases; n = 496 unaffected relatives) for replication. Two variants with supportive evidence were subsequently genotyped in a population-based case-control study (n = 1,155 incident prostate cancer cases; n = 1,060 age-matched controls) for further confirmation. All participants were men of European ancestry.

RESULTS

The strongest evidence was for two germline missense variants in the butyrophilin-like 2 (BTNL2) gene (rs41441651, p.Asp336Asn and rs28362675, p.Gly454Cys) that segregated with affection status in two of the WES families. In the independent set of 270 HPC families, 1.5% (rs41441651; P = 0.0032) and 1.2% (rs28362675; P = 0.0070) of affected men, but no unaffected men, carried a variant. Both variants were associated with elevated prostate cancer risk in the population-based study (rs41441651: OR, 2.7; 95% CI, 1.27-5.87; P = 0.010; rs28362675: OR, 2.5; 95% CI, 1.16-5.46; P = 0.019).

CONCLUSIONS

Results indicate that rare BTNL2 variants play a role in susceptibility to both familial and sporadic prostate cancer.

IMPACT

Results implicate BTNL2 as a novel prostate cancer susceptibility gene.