To discover new rheumatoid arthritis (RA) risk loci, we systematically examined 370 SNPs from 179 independent loci with P < 0.001 in a published meta-analysis of RA genome-wide association studies (GWAS) of 3,393 cases and 12,462 controls. We used Gene Relationships Across Implicated Loci (GRAIL), a computational method that applies statistical text mining to PubMed abstracts, to score these 179 loci for functional relationships to genes in 16 established RA disease loci. We identified 22 loci with a significant degree of functional connectivity. We genotyped 22 representative SNPs in an independent set of 7,957 cases and 11,958 matched controls. Three were convincingly validated: CD2-CD58 (rs11586238, P = 1 x 10(-6) replication, P = 1 x 10(-9) overall), CD28 (rs1980422, P = 5 x 10(-6) replication, P = 1 x 10(-9) overall) and PRDM1 (rs548234, P = 1 x 10(-5) replication, P = 2 x 10(-8) overall). An additional four were replicated (P < 0.0023): TAGAP (rs394581, P = 0.0002 replication, P = 4 x 10(-7) overall), PTPRC (rs10919563, P = 0.0003 replication, P = 7 x 10(-7) overall), TRAF6-RAG1 (rs540386, P = 0.0008 replication, P = 4 x 10(-6) overall) and FCGR2A (rs12746613, P = 0.0022 replication, P = 2 x 10(-5) overall). Many of these loci are also associated to other immunologic diseases.

Case-control genetic studies of complex human diseases can be confounded by population stratification. This issue can be addressed using panels of ancestry informative markers (AIMs) that can provide substantial population substructure information. Previously, we described a panel of 128 SNP AIMs that were designed as a tool for ascertaining the origins of subjects from Europe, Sub-Saharan Africa, Americas, and East Asia.

Accounting for population genetic substructure is important in reducing type 1 errors in genetic studies of complex disease. As efforts to understand complex genetic disease are expanded to different continental populations the understanding of genetic substructure within these continents will be useful in design and execution of association tests. In this study, population differentiation (Fst) and Principal Components Analyses (PCA) are examined using >200 K genotypes from multiple populations of East Asian ancestry. The population groups included those from the Human Genome Diversity Panel [Cambodian, Yi, Daur, Mongolian, Lahu, Dai, Hezhen, Miaozu, Naxi, Oroqen, She, Tu, Tujia, Naxi, Xibo, and Yakut], HapMap [ Han Chinese (CHB) and Japanese (JPT)], and East Asian or East Asian American subjects of Vietnamese, Korean, Filipino and Chinese ancestry. Paired Fst (Wei and Cockerham) showed close relationships between CHB and several large East Asian population groups (CHB/Korean, 0.0019; CHB/JPT, 00651; CHB/Vietnamese, 0.0065) with larger separation with Filipino (CHB/Filipino, 0.014). Low levels of differentiation were also observed between Dai and Vietnamese (0.0045) and between Vietnamese and Cambodian (0.0062). Similarly, small Fst's were observed among different presumed Han Chinese populations originating in different regions of mainland of China and Taiwan (Fst's <0.0025 with CHB). For PCA, the first two PC's showed a pattern of relationships that closely followed the geographic distribution of the different East Asian populations. PCA showed substructure both between different East Asian groups and within the Han Chinese population. These studies have also identified a subset of East Asian substructure ancestry informative markers (EASTASAIMS) that may be useful for future complex genetic disease association studies in reducing type 1 errors and in identifying homogeneous groups that may increase the power of such studies.

Primary biliary cholangitis (PBC) is a chronic autoimmune liver disease with progressive cholestasis and liver fibrosis. Similar to human patients with PBC, p40-/-IL-2Rα-/- mice spontaneously develop severe autoimmune cholangitis. Although there has been considerable work on immune regulation and autoimmunity, there is a relative paucity of work directed at the functional implications of the key peritoneal cavity (PC) B cell subset, coined B1a cells in PBC. We used flow cytometry and high-resolution microarrays to study the qualitative and quantitative characteristics of B cells, particularly B1a cells, in the PC of p40-/-IL-2Rα-/- mice compared to controls. Importantly, B1a cell proliferation was markedly lower as the expression of Ki67 decreased. Meanwhile, the apoptosis level was much higher. These lead to a reduction of B1a cells in the PC of p40-/-IL-2Rα-/- mice compared to controls. In contrast, there was a dramatic increase of CD4+ and CD8+ T cells accompanied by elevated production of IFN-γ. In addition, we found a negative correlation between the frequency of B1a cells and the presence of autoreactive CD8+ T cells in both liver and PC of p40-/-IL-2Rα-/- mice. From a functional perspective, B cells from p40-/-IL-2Rα-/- mice downregulated IL-10 production and CTLA-4 expression, leading to loss of B cell regulatory function. We suggest that the dysfunction of B1a cells in the PC in this murine model of autoimmune cholangitis results in defective regulatory function. This highlights a new potential therapeutic target in PBC.

Using the Immunochip custom SNP array, which was designed for dense genotyping of 186 loci identified through genome-wide association studies (GWAS), we analyzed 11,475 individuals with rheumatoid arthritis (cases) of European ancestry and 15,870 controls for 129,464 markers. We combined these data in a meta-analysis with GWAS data from additional independent cases (n = 2,363) and controls (n = 17,872). We identified 14 new susceptibility loci, 9 of which were associated with rheumatoid arthritis overall and five of which were specifically associated with disease that was positive for anticitrullinated peptide antibodies, bringing the number of confirmed rheumatoid arthritis risk loci in individuals of European ancestry to 46. We refined the peak of association to a single gene for 19 loci, identified secondary independent effects at 6 loci and identified association to low-frequency variants at 4 loci. Bioinformatic analyses generated strong hypotheses for the causal SNP at seven loci. This study illustrates the advantages of dense SNP mapping analysis to inform subsequent functional investigations.

Primary biliary cirrhosis (PBC) is considered a model autoimmune disease, with the most highly directed and specific autoantibody in both murine and human autoimmunity, the anti-mitochondrial autoantibody (AMA). However, therapeutic advances in this disease have lagged behind. Herein we have taken advantage of our unique model of murine PBC in which mice immunized with 2-octynoic acid coupled to BSA (2OA-BSA), a compound identified by quantitative structure activity relationships (QSAR) of human AMA binding, develop an intense inflammatory cholangitis with striking similarities to humans with PBC. In particular, we have constructed several unique gene-deleted mice, including mice deleted of IL-12p40, IL-12p35, IFN-γ, IL-23p19, IL-17A, IL-17F and IL-22, immunized these animals with 2OA-BSA and followed the natural history of immunopathology to identify key pathways that might provide clues for successful therapy. Our data indicate that whereas both IL-12/Th1 and IL-23/Th17 are involved in cholangitis, it is the IL-12/Th1 signaling pathway that elicits pathology. In fact, deletion of IFN-γ prevents disease and suppresses autoantibodies. Importantly, deletion of the Th17 cytokines IL-17A and IL-22, but not IL-17F, reduces biliary damage; IL-17A-knockout mice have reduced levels of anti-mitochondrial antibody. We further demonstrate that the production of IFN-γ is significantly decreased in the liver of IL-23p19(-/-), IL-17A(-/-) and IL-22(-/-) mice compared with controls. However, the ability of T cells to produce IFN-γ was not affected in Th17 cytokine-deficient mice. Our data indicate that a deficient Th17 pathway suppresses the accumulation of IFN-γ producing cells in liver during the early phase of cholangitis. In conclusion, whereas IFN-γ has a pivotal role in the early events involved in the pathogenesis of autoimmune cholangitis induced by 2OA-BSA, the IL-23/Th17 pathway potentiates the effects of IL-12/IFN-γ-mediated immunopathology.

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD.

Primary biliary cirrhosis (PBC) is a classical autoimmune liver disease for which effective immunomodulatory therapy is lacking. Here we perform meta-analyses of discovery data sets from genome-wide association studies of European subjects (n=2,764 cases and 10,475 controls) followed by validation genotyping in an independent cohort (n=3,716 cases and 4,261 controls). We discover and validate six previously unknown risk loci for PBC (Pcombined<5 × 10(-8)) and used pathway analysis to identify JAK-STAT/IL12/IL27 signalling and cytokine-cytokine pathways, for which relevant therapies exist.

The aryl hydrocarbon receptor repressor (AhRR) is known to repress aryl hydrocarbon receptor (AhR) signaling, but very little is known regarding the role of the AhRR in vivo.

Thymic stromal lymphopoietin (TSLP), an IL7-like cytokine produced by bronchial epithelial cells is upregulated in asthma and induces dendritic cell maturation supporting a Th2 response. Environmental pollutants, including tobacco smoke and diesel exhaust particles upregulate TSLP suggesting that TSLP may be an interface between environmental pollution and immune responses in asthma. Since asthma is prevalent in urban communities, variants in the TSLP gene may be important in asthma susceptibility in these populations.

The minor allele of the R620W missense single-nucleotide polymorphism (SNP) (rs2476601) in the hematopoietic-specific protein tyrosine phosphatase gene, PTPN22, has been associated with multiple autoimmune diseases, including rheumatoid arthritis (RA). These genetic data, combined with biochemical evidence that this SNP affects PTPN22 function, suggest that this phosphatase is a key regulator of autoimmunity. To determine whether other genetic variants in PTPN22 contribute to the development of RA, we sequenced the coding regions of this gene in 48 white North American patients with RA and identified 15 previously unreported SNPs, including 2 coding SNPs in the catalytic domain. We then genotyped 37 SNPs in or near PTPN22 in 475 patients with RA and 475 individually matched controls (sample set 1) and selected a subset of markers for replication in an additional 661 patients with RA and 1,322 individually matched controls (sample set 2). Analyses of these results predict 10 common (frequency >1%) PTPN22 haplotypes in white North Americans. The sole haplotype found to carry the previously identified W620 risk allele was strongly associated with disease in both sample sets, whereas another haplotype, identical at all other SNPs but carrying the R620 allele, showed no association. R620W, however, does not fully explain the association between PTPN22 and RA, since significant differences between cases and controls persisted in both sample sets after the haplotype data were stratified by R620W. Additional analyses identified two SNPs on a single common haplotype that are associated with RA independent of R620W, suggesting that R620W and at least one additional variant in the PTPN22 gene region influence RA susceptibility.

For admixture mapping studies in Mexican Americans (MAM), we define a genomewide single-nucleotide-polymorphism (SNP) panel that can distinguish between chromosomal segments of Amerindian (AMI) or European (EUR) ancestry. These studies used genotypes for >400,000 SNPs, defined in EUR and both Pima and Mayan AMI, to define a set of ancestry-informative markers (AIMs). The use of two AMI populations was necessary to remove a subset of SNPs that distinguished genotypes of only one AMI subgroup from EUR genotypes. The AIMs set contained 8,144 SNPs separated by a minimum of 50 kb with only three intermarker intervals >1 Mb and had EUR/AMI FST values >0.30 (mean FST = 0.48) and Mayan/Pima FST values <0.05 (mean FST < 0.01). Analysis of a subset of these SNP AIMs suggested that this panel may also distinguish ancestry between EUR and other disparate AMI groups, including Quechuan from South America. We show, using realistic simulation parameters that are based on our analyses of MAM genotyping results, that this panel of SNP AIMs provides good power for detecting disease-associated chromosomal segments for genes with modest ethnicity risk ratios. A reduced set of 5,287 SNP AIMs captured almost the same admixture mapping information, but smaller SNP sets showed substantial drop-off in admixture mapping information and power. The results will enable studies of type 2 diabetes, rheumatoid arthritis, and other diseases among which epidemiological studies suggest differences in the distribution of ancestry-associated susceptibility.

Mucosal-associated invariant T (MAIT) cells are novel innate-like T cells constituting a significant proportion of circulating and hepatic T cells. Herein, we extensively examine the phenotypical and functional alterations of MAIT cells and their regulation in a cohort of 56 patients with Primary Biliary Cholangitis (PBC) and 53 healthy controls (HC). Additionally alterations of MAIT cells were assessed before and after UDCA treatment. Finally the localization of MAIT cell in liver was examined using specific tetramer staining and the underlying mechanisms of these alterations in PBC were explored. Our data demonstrated that the frequency and number of circulating MAIT cells were decreased, whereas hepatic MAIT cells were increased in PBC compared to HC. Moreover, circulating MAIT cells were more activated in PBC than HC, reflected by elevated expression levels of granzyme B. Six months of UDCA treatment significantly attenuated the circulating MAIT cells differences in PBC. Of note, the expression levels of IL-7 were significantly increased in both plasma and liver from PBC as compared to HC, which promoted the production of inflammatory cytokines and granzyme B by inducing signal transduction and activation of transcription 5 (STAT5) phosphorylation in MAIT cells. Finally, cholic acid, one of the major bile acids in liver, upregulated IL-7 expression in hepatocyte cell line L02 by inducing Farnesoid X Receptor (FXR) binding to the IL-7 promoter. Hence MAIT cells are activated and enriched in the liver of PBC. Cholic acid-induced IL-7 production in hepatocytes plays a critical role in regulating MAIT cell function, highlighting that hepatocytes may bridge cholangiocyte injury and innate immunity through a bile acid signaling pathway.

The immune and nervous systems are two major organ systems responsible for host defense and memory. Both systems achieve memory and learning that can be retained, retrieved, and utilized for decades. Here, we report the surprising discovery that peripheral sensory neurons of the dorsal root ganglia (DRGs) of immunized mice contain antigen-specific antibodies. Using a combination of rigorous molecular genetic analyses, transgenic mice, and adoptive transfer experiments, we demonstrate that DRGs do not synthesize these antigen-specific antibodies, but rather sequester primarily IgG1 subtype antibodies. As revealed by RNA-seq and targeted quantitative PCR (qPCR), dorsal root ganglion (DRG) sensory neurons harvested from either naïve or immunized mice lack enzymes (i.e., RAG1, RAG2, AID, or UNG) required for generating antibody diversity and, therefore, cannot make antibodies. Additionally, transgenic mice that express a reporter fluorescent protein under the control of Igγ1 constant region fail to express Ighg1 transcripts in DRG sensory neurons. Furthermore, neural sequestration of antibodies occurs in mice rendered deficient in neuronal Rag2, but antibody sequestration is not observed in DRG sensory neurons isolated from mice that lack mature B cells [e.g., Rag1 knock out (KO) or μMT mice]. Finally, adoptive transfer of Rag1-deficient bone marrow (BM) into wild-type (WT) mice or WT BM into Rag1 KO mice revealed that antibody sequestration was observed in DRG sensory neurons of chimeric mice with WT BM but not with Rag1-deficient BM. Together, these results indicate that DRG sensory neurons sequester and retain antigen-specific antibodies released by antibody-secreting plasma cells. Coupling this work with previous studies implicating DRG sensory neurons in regulating antigen trafficking during immunization raises the interesting possibility that the nervous system collaborates with the immune system to regulate antigen-mediated responses.

Primary biliary cholangitis (PBC) is an autoimmune liver disease with a strong hereditary component. Here, we report a genome-wide association study that included 1,122 PBC cases and 4,036 controls of Han Chinese descent, with subsequent replication in a separate cohort of 907 PBC cases and 2,127 controls. Our results show genome-wide association of 14 PBC risk loci including previously identified 6p21 (HLA-DRA and DPB1), 17q12 (ORMDL3), 3q13.33 (CD80), 2q32.3 (STAT1/STAT4), 3q25.33 (IL12A), 4q24 (NF-κB) and 22q13.1 (RPL3/SYNGR1). We also identified variants in IL21, IL21R, CD28/CTLA4/ICOS, CD58, ARID3A and IL16 as novel PBC risk loci. These new findings and histochemical studies showing enhanced expression of IL21 and IL21R in PBC livers (particularly in the hepatic portal tracks) support a disease mechanism in which the deregulation of the IL21 signalling pathway, in addition to CD4 T-cell activation and T-cell co-stimulation are critical components in the development of PBC.

Endometriosis is a chronic and underdiagnosed disease which affects 5-10% of women of childbearing age and is characterized by growth of endometrial tissue outside of the uterus, most often in the peritoneal cavity. Delay in diagnosis is a major problem for management of this disorder, and treatment is often not initiated until the disease has progressed for many years. Although the exact etiology of endometriosis remains unknown, retrograde menstruation is recognized as a common underlying factor leading to the deposit of menstrual effluent (ME) into the peritoneal cavity. Differences in the cellular biology and genetics of the cells within ME are therefore likely to explain why endometriosis develops in only a subset of women.

Can a previously defined relationship between sperm capacitation and the probability of a man generating pregnancy within three cycles, prospectively predict male fertility in diverse clinical settings? A second study asked, what is the prevalence of impaired sperm fertilizing ability in men questioning their fertility (MQF), and does this relate to traditional semen analysis metrics?

Regulatory B cells (Breg) are considered as immunosuppressive cells. Different subsets of Breg cells have been identified both in human beings and in mice. However, there is a lack of unique markers to identify Breg cells, and the heterogeneity of Breg cells in different organs needs to be further illuminated. In this study, we performed high-throughput single-cell RNA sequencing (scRNA-seq) and single-cell B-cell receptor sequencing (scBCR-seq) of B cells from the murine spleen, liver, mesenteric lymph nodes, bone marrow, and peritoneal cavity to better define the phenotype of these cells. Breg cells were identified based on the expression of immunosuppressive genes and IL-10-producing B (B10) cell-related genes, to define B10 and non-B10 subsets in Breg cells based on the score of the B10 gene signatures. Moreover, we characterized 19 common genes significantly expressed in Breg cells, including Fcrl5, Zbtb20, Ccdc28b, Cd9, and Ptpn22, and further analyzed the transcription factor activity in defined Breg cells. Last, a BCR analysis was used to determine the clonally expanded clusters and the relationship of Breg cells across different organs. We demonstrated that Atf3 may potentially modulate the function of Breg cells as a transcription factor and that seven organ-specific subsets of Breg cells are found. Depending on gene expression and functional modules, non-B10 Breg cells exhibited activated the TGF-β pathway, thus suggesting that non-B10 Breg cells have specific immunosuppressive properties different from conventional B10 cells. In conclusion, our work provides new insights into Breg cells and illustrates their transcriptional profiles and BCR repertoire in different organs under physiological conditions.

Current treatment guidelines recommend ursodeoxycholic acid (UDCA) as the first-line treatment for new-diagnosed primary biliary cholangitis (PBC) patients. However, up to 40% patients are insensitive to UDCA monotherapy, and evaluation of UDCA response at 12 months may result in long period of ineffective treatment. We aimed to develop a new criterion to reliably identify non-response patients much earlier.

Invariant natural killer T (iNKT) cells play complex roles in bridging innate and adaptive immunity by engaging with glycolipid antigens presented by CD1d. Our earlier work suggested that iNKT cells were involved in the initiation of the original loss of tolerance in primary biliary cirrhosis (PBC). To address this issue in more detail and, in particular, to focus on whether iNKT cells activated by a Th2-biasing agonist (2s,3s,4r)-1-O-(α-D-galactopyranosyl)-N-tetracosanoyl-2-amino-1,3,4-nonanetriol (OCH), can influence the development of PBC in a xenobiotic-induced PBC murine model. Groups of mice were treated with either OCH or, as a control, α-galactosylceramide (α-GalCer) and thence serially followed for cytokine production, markers of T cell activation, liver histopathology and anti-mitochondrial antibody responses. Further, additional groups of CD1d deleted mice were similarly studied. Our data indicate that administration of OCH has a dramatic influence with exacerbation of portal inflammation and hepatic fibrosis similar to mice treated with α-GalCer. Further, iNKT cell deficient CD1d knockout mice have decreased inflammatory portal cell infiltrates and reduced anti-mitochondrial antibody responses. We submit that activation of iNKT cells can occur via overlapping and/or promiscuous pathways and highlight the critical role of innate immunity in the natural history of autoimmune cholangitis. These data have implications for humans with PBC and emphasize that therapeutic strategies must focus not only on suppressing adaptive responses, but also innate immunity.