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.

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 adaptive immune system uses several strategies to generate a repertoire of T cell receptors (TCR) with sufficient diversity to recognize the universe of potential pathogens. However, it remains unclear how differences in the T cell receptor (TCR) contribute to heterogeneity in T cell state. In this study, we used polychromatic flow cytometry to isolate highly pure CD4+/CD8+ naive and memory T cells, and applied deep sequencing to characterize corresponding TCR β-chain (TCRβ) complementary-determining region 3 (CDR3) repertoires. We find that shorter TCRβ CDR3s with fewer insertions were highly enriched during thymic selection. Antigen-experienced T cells (memory T cells) harbor shorter CDR3s vs. naive T cells. Moreover, the public TCRβ CDR3 clonotypes within cell subsets or interindividual tend to have shorter CDR3 length and a significantly larger size compared with "private" clonotypes. Taken together, shorter CDR3s highly enriched during thymic selection and antigen-driven selection, and further enriched in public T-cell responses. These results indicated that it may be evolutionary pressures drive short CDR3s to recognize most of antigen in nature.

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.

We previously reported that the G allele of rs3853839 at 3'untranslated region (UTR) of Toll-like receptor 7 (TLR7) was associated with elevated transcript expression and increased risk for systemic lupus erythematosus (SLE) in 9,274 Eastern Asians [P = 6.5×10(-10), odds ratio (OR) (95%CI) = 1.27 (1.17-1.36)]. Here, we conducted trans-ancestral fine-mapping in 13,339 subjects including European Americans, African Americans, and Amerindian/Hispanics and confirmed rs3853839 as the only variant within the TLR7-TLR8 region exhibiting consistent and independent association with SLE (Pmeta = 7.5×10(-11), OR = 1.24 [1.18-1.34]). The risk G allele was associated with significantly increased levels of TLR7 mRNA and protein in peripheral blood mononuclear cells (PBMCs) and elevated luciferase activity of reporter gene in transfected cells. TLR7 3'UTR sequence bearing the non-risk C allele of rs3853839 matches a predicted binding site of microRNA-3148 (miR-3148), suggesting that this microRNA may regulate TLR7 expression. Indeed, miR-3148 levels were inversely correlated with TLR7 transcript levels in PBMCs from SLE patients and controls (R(2) = 0.255, P = 0.001). Overexpression of miR-3148 in HEK-293 cells led to significant dose-dependent decrease in luciferase activity for construct driven by TLR7 3'UTR segment bearing the C allele (P = 0.0003). Compared with the G-allele construct, the C-allele construct showed greater than two-fold reduction of luciferase activity in the presence of miR-3148. Reduced modulation by miR-3148 conferred slower degradation of the risk G-allele containing TLR7 transcripts, resulting in elevated levels of gene products. These data establish rs3853839 of TLR7 as a shared risk variant of SLE in 22,613 subjects of Asian, EA, AA, and Amerindian/Hispanic ancestries (Pmeta  = 2.0×10(-19), OR = 1.25 [1.20-1.32]), which confers allelic effect on transcript turnover via differential binding to the epigenetic factor miR-3148.

Myeloid dendritic cells (mDC) are lost from blood in individuals with human immunodeficiency virus (HIV) infection but the mechanism for this loss and its relationship to disease progression are not known. We studied the mDC response in blood and lymph nodes of simian immunodeficiency virus (SIV)-infected rhesus macaques with different disease outcomes. Early changes in blood mDC number were inversely correlated with virus load and reflective of eventual disease outcome, as animals with stable infection that remained disease-free for more than one year had average increases in blood mDC of 200% over preinfection levels at virus set-point, whereas animals that progressed rapidly to AIDS had significant loss of mDC at this time. Short term antiretroviral therapy (ART) transiently reversed mDC loss in progressor animals, whereas discontinuation of ART resulted in a 3.5-fold increase in mDC over preinfection levels only in stable animals, approaching 10-fold in some cases. Progressive SIV infection was associated with increased CCR7 expression on blood mDC and an 8-fold increase in expression of CCL19 mRNA in lymph nodes, consistent with increased mDC recruitment. Paradoxically, lymph node mDC did not accumulate in progressive infection but rather died from caspase-8-dependent apoptosis that was reduced by ART, indicating that increased recruitment is offset by increased death. Lymph node mDC from both stable and progressor animals remained responsive to exogenous stimulation with a TLR7/8 agonist. These data suggest that mDC are mobilized in SIV infection but that an increase in the CCR7-CCL19 chemokine axis associated with high virus burden in progressive infection promotes exodus of activated mDC from blood into lymph nodes where they die from apoptosis. We suggest that inflamed lymph nodes serve as a sink for mDC through recruitment, activation and death that contributes to AIDS pathogenesis.

Recent genome-wide association studies (GWASs) conducted in Asian populations have identified novel risk loci for systemic lupus erythematosus (SLE). Here, we genotyped 10 single-nucleotide polymorphisms (SNPs) in eight such loci and investigated their disease associations in three independent Caucasian SLE case-control cohorts recruited from Sweden, Finland and the United States. The disease associations of the SNPs in ETS1, IKZF1, LRRC18-WDFY4, RASGRP3, SLC15A4, TNIP1 and 16p11.2 were replicated, whereas no solid evidence of association was observed for the 7q11.23 locus in the Caucasian cohorts. SLC15A4 was significantly associated with renal involvement in SLE. The association of TNIP1 was more pronounced in SLE patients with renal and immunological disorder, which is corroborated by two previous studies in Asian cohorts. The effects of all the associated SNPs, either conferring risk for or being protective against SLE, were in the same direction in Caucasians and Asians. The magnitudes of the allelic effects for most of the SNPs were also comparable across different ethnic groups. On the contrary, remarkable differences in allele frequencies between Caucasian and Asian populations were observed for all associated SNPs. In conclusion, most of the novel SLE risk loci identified by GWASs in Asian populations were also associated with SLE in Caucasian populations. We observed both similarities and differences with respect to the effect sizes and risk allele frequencies across ethnicities.

The transporter associated with antigen processing (TAP) translocates antigenic peptides from the cytosol into the lumen of the endoplasmic reticular and plays a critical role in the major histocompatibility complex (MHC) class I molecule-mediated antigenic presentation pathway. In this study, the porcine TAP1 gene was mapped to the pig chromosome 7 (SSC7) and was closely linked to the marker SSC2B02 (retention fraction=43%, LOD=15.18). Subcellular localization of TAP1 by transient transfection of PK15 cells indicated that the TAP1 protein might be located in the endoplasmic reticulum (ER) in pig kidney epithelial cells (PK-15). Gene expression analysis by semi-quantitative RT-PCR revealed that TAP1 was selectively expressed in some immune and immune-related tissues. Quantitative real-time PCR (qRT-PCR) analysis revealed that this gene was up-regulated after treatments that mimic viral and bacterial infection (polyriboinosinic-polyribocytidylic acid (poly(I:C)) and lipopolysaccharide (LPS), respectively). In addition, elevated TAP1 expression was detected after porcine reproductive and respiratory syndrome virus (PRRSV) infection in porcine white blood cells (WBCs). One single nucleotide polymorphism (SNP) in exon 3 of TAP1 was detected in a Landrace pig population by Bsp143I restriction enzyme digestion. Different genotypes of this SNP had significant associations (P<0.05) with the red blood cell distribution width (RDW) of 1-day-old (1 d) pigs (P=0.0168), the PRRSV antibody level (PRRSV Ab) (P=0.0445) and the absolute lymphocyte count (LYM#) (P=0.024) of 17 d pigs. Our results showed that the TAP1 gene might have important roles in swine immune responses, and these results provide useful information for further functional studies.

Systemic lupus erythematosus (SLE) is considered to be the prototypic autoimmune disease, with a complex genetic architecture influenced by environmental factors. We sought to replicate a putative association at 11p13 not yet exceeding genome-wide significance (p < 5 × 10(-8)) identified in a genome-wide association study (GWAS). Our GWA scan identified two intergenic SNPs located between PDHX and CD44 showing suggestive evidence of association with SLE in cases of European descent (rs2732552, p = 0.004, odds ratio [OR] = 0.78; rs387619, p = 0.003, OR = 0.78). The replication cohort consisted of >15,000 subjects, including 3562 SLE cases and 3491 controls of European ancestry, 1527 cases and 1811 controls of African American (AA) descent, and 1265 cases and 1260 controls of Asian origin. We observed robust association at both rs2732552 (p = 9.03 × 10(-8), OR = 0.83) and rs387619 (p = 7.7 × 10(-7), OR = 0.83) in the European samples with p(meta) = 1.82 × 10(-9) for rs2732552. The AA and Asian SLE cases also demonstrated association at rs2732552 (p = 5 × 10(-3), OR = 0.81 and p = 4.3 × 10(-4), OR = 0.80, respectively). A meta-analysis of rs2732552 for all racial and ethnic groups studied produced p(meta) = 2.36 × 10(-13). This locus contains multiple regulatory sites that could potentially affect expression and functions of CD44, a cell-surface glycoprotein influencing immunologic, inflammatory, and oncologic phenotypes, or PDHX, a subunit of the pyruvate dehydrogenase complex.

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.

Many lines of evidence implicate mitochondria in phenotypic variation: (a) rare mutations in mitochondrial proteins cause metabolic, neurological, and muscular disorders; (b) alterations in oxidative phosphorylation are characteristic of type 2 diabetes, Parkinson disease, Huntington disease, and other diseases; and (c) common missense variants in the mitochondrial genome (mtDNA) have been implicated as having been subject to natural selection for adaptation to cold climates and contributing to "energy deficiency" diseases today. To test the hypothesis that common mtDNA variation influences human physiology and disease, we identified all 144 variants with frequency >1% in Europeans from >900 publicly available European mtDNA sequences and selected 64 tagging single-nucleotide polymorphisms that efficiently capture all common variation (except the hypervariable D-loop). Next, we evaluated the complete set of common mtDNA variants for association with type 2 diabetes in a sample of 3,304 diabetics and 3,304 matched nondiabetic individuals. Association of mtDNA variants with other metabolic traits (body mass index, measures of insulin secretion and action, blood pressure, and cholesterol) was also tested in subsets of this sample. We did not find a significant association of common mtDNA variants with these metabolic phenotypes. Moreover, we failed to identify any physiological effect of alleles that were previously proposed to have been adaptive for energy metabolism in human evolution. More generally, this comprehensive association-testing framework can readily be applied to other diseases for which mitochondrial dysfunction has been implicated.

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.

Psoriasis is a common immune-mediated disorder that affects the skin, nails and joints. To identify psoriasis susceptibility loci, we genotyped 438,670 SNPs in 1,409 psoriasis cases and 1,436 controls of European ancestry. We followed up 21 promising SNPs in 5,048 psoriasis cases and 5,041 controls. Our results provide strong support for the association of at least seven genetic loci and psoriasis (each with combined P < 5 x 10(-8)). Loci with confirmed association include HLA-C, three genes involved in IL-23 signaling (IL23A, IL23R, IL12B), two genes that act downstream of TNF-alpha and regulate NF-kappaB signaling (TNIP1, TNFAIP3) and two genes involved in the modulation of Th2 immune responses (IL4, IL13). Although the proteins encoded in these loci are known to interact biologically, we found no evidence for epistasis between associated SNPs. Our results expand the catalog of genetic loci implicated in psoriasis susceptibility and suggest priority targets for study in other auto-immune disorders.

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.

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.

Experimental evidence indicated that interferon-inducible guanylate-binding proteins (GBPs) are similar with genes in the myxovirus (Mx) resistance protein subfamily of the large GTPases protein family and play important roles in the resistance to intracellular pathogens. As more diseases exert significant influence on pig industry, it is anticipated that more candidate disease-resistance genes could be found in future strategies aimed at improving genetic resistance to infectious diseases. In this study, we cloned cDNA sequences and analyzed the genomic structure of porcine GBP1 (poGBP1) and GBP2 (poGBP2). The two genes were mapped to SSC4q21-q23 and SSC4q24 by the SCHP panel respectively, further IMpRH panel analysis showed both genes were most closely linked to the marker SWR153. The deduced amino acid sequences of these two genes share the same three classical GTP-binding motifs at the amino terminus and are less conserved at the carboxyl termini except for a CaaX motif, compared with human and mouse counterparts. The reverse transcriptase-polymerase chain reaction (RT-PCR) revealed that poGBP1 and poGBP2 were both widely expressed in many tissues, and transient transfection indicated that poGBP1 and poGBP2 proteins were both located in cytoplasms within Pig Kidney Epithelial cells (PK15). Quantitative real-time PCR (Q-RT-PCR) analyses showed poGBP1 and poGBP2 had very similar expression patterns in PK15 cells at different time points after poly I:C stimulation, suggesting that ISRE (interferon-stimulated response element) plays a crucial role in the transcriptional regulation of these two genes. Four single nucleotide polymorphisms (SNPs) and three SNPs were detected in the cDNA sequences of poGBP1 and poGBP2, respectively. Association analyses revealed that the poGBP1 Eco81I and poGBP2 SspI polymorphisms both had significant associations (p<0.05) with red blood cell count (RBC), haemoglobin concentration (HGB) and hematocrit (HCT) of 17-day-old pigs.

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.

We applied nonparametric quantitative trait linkage analysis to two rheumatoid arthritis quantitative phenotypes, IgM rheumatoid factor (RF) and anti-cyclic citrullinated peptide autoantibody titer measurements, using 5700 genome-wide Illumina single-nucleotide polymorphism genotypes on 658 Caucasian North American Rheumatoid Arthritis Consortium families. Peak LOD scores for both quantitative traits were located in the human leukocyte antigen region 6p21 (15.8 and 13.8 for RF and anti-cyclic citrullinated peptide, respectively) followed by 11p12 (3.2 and 3.6). In addition, there were LOD scores of 3.2 on 2q32 for RF and 3.6 on 4q24 for anti-cyclic citrullinated peptide. The resulting linkage signals for both phenotypes are very similar to previous results for rheumatoid arthritis as a qualitative variable, with rheumatoid factor measurements being most closely aligned. Interestingly, anti-cyclic citrullinated peptide exhibits a stronger linkage peak on 2p14 than rheumatoid factor and rheumatoid arthritis, and stronger linkage on 4q24. Finally, we used ordered subset analyses to determine if sub-ranges of these two traits increased rheumatoid arthritis linkage signals; however, our analyses did not reveal significant effects of the quantitative traits on rheumatoid arthritis linkage signals in this population.

Droplet single-cell RNA-sequencing (dscRNA-seq) has enabled rapid, massively parallel profiling of transcriptomes. However, assessing differential expression across multiple individuals has been hampered by inefficient sample processing and technical batch effects. Here we describe a computational tool, demuxlet, that harnesses natural genetic variation to determine the sample identity of each droplet containing a single cell (singlet) and detect droplets containing two cells (doublets). These capabilities enable multiplexed dscRNA-seq experiments in which cells from unrelated individuals are pooled and captured at higher throughput than in standard workflows. Using simulated data, we show that 50 single-nucleotide polymorphisms (SNPs) per cell are sufficient to assign 97% of singlets and identify 92% of doublets in pools of up to 64 individuals. Given genotyping data for each of eight pooled samples, demuxlet correctly recovers the sample identity of >99% of singlets and identifies doublets at rates consistent with previous estimates. We apply demuxlet to assess cell-type-specific changes in gene expression in 8 pooled lupus patient samples treated with interferon (IFN)-β and perform eQTL analysis on 23 pooled samples.

Epigenetic modifications have previously been associated with rheumatoid arthritis (RA). In this study, we aimed to determine whether differential DNA methylation in peripheral blood cell subpopulations is associated with any of 4 clinical outcomes among RA patients.