A regulatory variant at 19p13.3 is associated with primary biliary cholangitis risk and ARID3A expression

Genome-wide association studies have identified 19p13.3 locus associated with primary biliary cholangitis (PBC). Here we aim to identify causative variant(s) and initiate efforts to define the mechanism by which the 19p13.3 locus variant(s) contributes to the pathogenesis of PBC. A genome-wide meta-analysis of 1931 PBC subjects and 7852 controls in two Han Chinese cohorts confirms the strong association between 19p13.3 locus and PBC. By integrating functional annotations, luciferase reporter assay and allele-specific chromatin immunoprecipitation, we prioritize rs2238574, an AT-Rich Interaction Domain 3A (ARID3A) intronic variant, as a potential causal variant at 19p13.3 locus. The risk allele of rs2238574 shows higher binding affinity of transcription factors, leading to an increased enhancer activity in myeloid cells. Genome-editing demonstrates the regulatory effect of rs2238574 on ARID3A expression through allele-specific enhancer activity. Furthermore, knock-down of ARID3A inhibits myeloid differentiation and activation pathway, and overexpression of the gene has the opposite effect. Finally, we find ARID3A expression and rs2238574 genotypes linked to disease severity in PBC. Our work provides several lines of evidence that a non-coding variant regulates ARID3A expression, presenting a mechanistic basis for association of 19p13.3 locus with the susceptibility to PBC.

Genetic admixture and cardiovascular disease risk in postmenopausal Hispanic women

BACKGROUND

Hispanics are a heterogeneous population with differences in the prevalence of cardiovascular disease (CVD) and its related risk factors among ethnic sub-groups. This study evaluated the association of genetic admixture and CVD in self-identified Hispanic women from the Women's Health Initiative (WHI).

METHODS

Data came from the WHI Observational Study and the Clinical Trial Components conducted among postmenopausal women. The CVD outcomes included coronary heart disease (CHD) and stroke. The proportions of European (EUR), sub-Saharan African (AFR), and Amerindian (AMI) admixture were estimated using 92 ancestry-informative markers. Cox regression models were used to assess the relationship between genetic admixture and CVD adjusting for age, lifestyle risk factors, known risk factors, and neighborhood socioeconomic status.

RESULTS

Among 5195 participants EUR ancestry was associated with a lower CHD risk after adjusting for age (HR 0.41, p = 0.02), and in the fully adjusted model (HR 0.40, p = 0.03). AFR ancestry was associated with a higher CHD risk after adjusting for age (HR 2.91, p = 0.03), but it only showed a trend in in the fully adjusted model (HR 2.46, p = 0.10). AMI ancestry was not statistically significantly associated with CHD and none of the genetic admixture proportions were statistically significantly associated with stroke (p > 0.05).

CONCLUSION

EUR ancestry was associated with a lower risk of CHD in Hispanic women. This highlights the need to account for genetic admixture in future CVD studies to consider different heritage groups to understand the role that genetic, neighborhood socioeconomic status, and environmental factors contribute to CVD health disparities in Hispanic women.

Genome-wide meta-analysis identifies susceptibility loci for autoimmune hepatitis type 1

BACKGROUND AND AIMS

Autoimmune hepatitis (AIH) is a rare and chronic autoimmune liver disease. While genetic factors are believed to play a crucial role in the etiopathogenesis of AIH, our understanding of these genetic risk factors is still limited. In this study, we aimed to identify susceptibility loci to further understand the pathogenesis of this disease.

APPROACH AND RESULTS

We conducted a case-control association study of 1,622 Chinese patients with AIH type 1 and 10,466 population controls from two independent cohorts. A meta-analysis was performed to ascertain variants associated with AIH type 1. A single-nucleotide polymorphism within the human leukocyte antigen (HLA) region showed the strongest association with AIH (rs6932730: OR = 2.32; p = 9.21 × 10^-73 ). The meta-analysis also identified two non-HLA loci significantly associated with AIH: CD28/CTLA4/ICOS on 2q33.3 (rs72929257: OR = 1.31; p = 2.92 × 10^-9 ) and SYNPR on 3p14.2 (rs6809477: OR = 1.25; p = 5.48 × 10^-9 ). In silico annotation, reporter gene assays, and CRISPR activation experiments identified a distal enhancer at 2q33.3 that regulated expression of CTLA4. In addition, variants near STAT1/STAT4 (rs11889341: OR = 1.24; p = 1.34 × 10^-7 ), LINC00392 (rs9564997: OR = 0.81; p = 2.53 × 10^-7 ), IRF8 (rs11117432: OR = 0.72; p = 6.10 × 10^-6 ), and LILRA4/LILRA5 (rs11084330: OR = 0.65; p = 5.19 × 10^-6 ) had suggestive association signals with AIH.

CONCLUSIONS

Our study identifies two novel loci (CD28/CTLA4/ICOS and SYNPR) exceeding genome-wide significance and suggests four loci as potential risk factors. These findings highlight the importance of costimulatory signaling and neuro-immune interaction in the pathogenesis of AIH.

Increased sensitivity of gp210 autoantibody detection using a newly designed gp210 antigen

OBJECTIVES

The detection of autoantibody to glycoprotein 210 (gp210 Ab) against a 15 amino-acid peptide epitope by enzyme-linked immunosorbent assay (ELISA) has been widely used in the diagnosis of primary biliary cholangitis (PBC). However, this small peptide antigen presents spatial limitations for antibody access, which reduces the sensitivity of autoantibody detection. A recombinant gp210 antigen was constructed for increased sensitivity in antibody detection is described here.

METHODS

The gp210 C terminal 18 amino acid coding sequence was ligated to the modified C-terminal 108 amino acid coding sequence of human serum albumin (mHSA108) and produced as a recombinant gp210 antigen mHSA108-gp210-C18. Measurements of gp210 Ab using the gp210 C-terminal 25 amino acid peptide (gp210-C25) and mHSA108-gp210-C18 by in-house ELISA were compared. ELISAs with mHSA108-gp210-C18 and commercial INOVA kit for gp210 Ab detection were also compared in PBC patients and healthy controls. The correlation between the two assays was analyzed and their efficiency in diagnosing was compared.

RESULTS

Of 86 PBC samples, 35 (40.70%) and 44 (52.33%) positive samples were detected for anti-gp210 Ab using gp210-C25 and mHSA108-gp210-C18, respectively. Of 252 samples from PBC, 114 (45.24%) were positive for mHSA108-gp210-C18 ELISA whereas 94 (37.3%) for commercial ELISA (INOVA). All positive samples detected with commercial ELISA kit were also tested positive in mHSA108-gp210-C18 ELISA. Among 374 patients with other autoimmune diseases, anti-gp210 Ab were detected by mHSA108-gp210-C18 ELISA in 0.95% systemic lupus erythematosus (SLE) patients (2/210), 13.04% rheumatoid arthritis (RA) patients (13/97), and 1.47% of Sjögren's Syndrome (SS) patients (1/67).

CONCLUSIONS

Compared to the gp210 peptide antigen, the sensitivity of the ELISA system using mHSA108-gp210-C18 antigen was improved. The novel gp210 antigen could be useful for screening patients known to be at increased risk of developing PBC.

An international genome-wide meta-analysis of primary biliary cholangitis: Novel risk loci and candidate drugs

BACKGROUNDS & AIMS

Primary biliary cholangitis (PBC) is a chronic liver disease in which autoimmune destruction of the small intrahepatic bile ducts eventually leads to cirrhosis. Many patients have inadequate response to licensed medications, motivating the search for novel therapies. Previous genome-wide association studies (GWAS) and meta-analyses (GWMA) of PBC have identified numerous risk loci for this condition, providing insight into its aetiology. We undertook the largest GWMA of PBC to date, aiming to identify additional risk loci and prioritise candidate genes for in silico drug efficacy screening.

METHODS

We combined new and existing genotype data for 10,516 cases and 20,772 controls from 5 European and 2 East Asian cohorts.

RESULTS

We identified 56 genome-wide significant loci (20 novel) including 46 in European, 13 in Asian, and 41 in combined cohorts; and a 57th genome-wide significant locus (also novel) in conditional analysis of the European cohorts. Candidate genes at newly identified loci include FCRL3, INAVA, PRDM1, IRF7, CCR6, CD226, and IL12RB1, which each play key roles in immunity. Pathway analysis reiterated the likely importance of pattern recognition receptor and TNF signalling, JAK-STAT signalling, and differentiation of T helper (TH)1 and TH17 cells in the pathogenesis of this disease. Drug efficacy screening identified several medications predicted to be therapeutic in PBC, some of which are well-established in the treatment of other autoimmune disorders.

CONCLUSIONS

This study has identified additional risk loci for PBC, provided a hierarchy of agents that could be trialled in this condition, and emphasised the value of genetic and genomic approaches to drug discovery in complex disorders.

LAY SUMMARY

Primary biliary cholangitis (PBC) is a chronic liver disease that eventually leads to cirrhosis. In this study, we analysed genetic information from 10,516 people with PBC and 20,772 healthy individuals recruited in Canada, China, Italy, Japan, the UK, or the USA. We identified several genetic regions associated with PBC. Each of these regions contains several genes. For each region, we used diverse sources of evidence to help us choose the gene most likely to be involved in causing PBC. We used these 'candidate genes' to help us identify medications that are currently used for treatment of other conditions, which might also be useful for treatment of PBC.

X Chromosome Contribution to the Genetic Architecture of Primary Biliary Cholangitis

BACKGROUND & AIMS

Genome-wide association studies in primary biliary cholangitis (PBC) have failed to find X chromosome (chrX) variants associated with the disease. Here, we specifically explore the chrX contribution to PBC, a sexually dimorphic complex autoimmune disease.

METHODS

We performed a chrX-wide association study, including genotype data from 5 genome-wide association studies (from Italy, United Kingdom, Canada, China, and Japan; 5244 case patients and 11,875 control individuals).

RESULTS

Single-marker association analyses found approximately 100 loci displaying P < 5 × 10-4, with the most significant being a signal within the OTUD5 gene (rs3027490; P = 4.80 × 10-6; odds ratio [OR], 1.39; 95% confidence interval [CI], 1.028-1.88; Japanese cohort). Although the transethnic meta-analysis evidenced only a suggestive signal (rs2239452, mapping within the PIM2 gene; OR, 1.17; 95% CI, 1.09-1.26; P = 9.93 × 10-8), the population-specific meta-analysis showed a genome-wide significant locus in East Asian individuals pointing to the same region (rs7059064, mapping within the GRIPAP1 gene; P = 6.2 × 10-9; OR, 1.33; 95% CI, 1.21-1.46). Indeed, rs7059064 tags a unique linkage disequilibrium block including 7 genes: TIMM17B, PQBP1, PIM2, SLC35A2, OTUD5, KCND1, and GRIPAP1, as well as a superenhancer (GH0XJ048933 within OTUD5) targeting all these genes. GH0XJ048933 is also predicted to target FOXP3, the main T-regulatory cell lineage specification factor. Consistently, OTUD5 and FOXP3 RNA levels were up-regulated in PBC case patients (1.75- and 1.64-fold, respectively).

CONCLUSIONS

This work represents the first comprehensive study, to our knowledge, of the chrX contribution to the genetics of an autoimmune liver disease and shows a novel PBC-related genome-wide significant locus.

Fine mapping of the MHC region identifies major independent variants associated with Han Chinese primary biliary cholangitis

The genetic association of primary biliary cholangitis with major histocompatibility complex (MHC) has been widely confirmed among different ethnicities. To map specific MHC region variants associated with PBC in a Han Chinese cohort, we imputed HLA antigens and amino acids (AA) in 1126 PBC cases and 1770 healthy control subjects using a Han-MHC reference database. We demonstrate that HLA-DRB1 and/or HLA-DQB1 contributed the strongest signals, and that HLA-DPB1 was a separate independent locus. Regression analyses with classical HLA alleles indicate that HLA-DQB1*03:01 or HLA-DQβ1-Pro55, HLA-DPB1*17:01 or HLA-DPβ1-Asp84 and HLA-DRB1*08:03 could largely explain MHC association with PBC. Forward stepwise regression analyses with HLA amino acid variants localize the major signals to HLA-DRβ1-Ala74, HLA-DQβ1-Pro55 and HLA-DPβ1-Asp84. Electrostatic potential calculations implicated AA variations at HLA-DQβ1 position 55 and HLA-DPβ1 position 84 as critical to peptide binding properties. Furthermore, although several critical Han Chinese AA variants differed from those shown in European populations, the predicted effects on antigen binding are likely to be very similar or identical and underlie the major component of MHC association with PBC.

Genome-wide Association Studies of Specific Antinuclear Autoantibody Subphenotypes in Primary Biliary Cholangitis

Anti-nuclear antibodies to speckled 100 kDa (sp100) and glycoprotein 210 (gp210) are specific serologic markers of primary biliary cholangitis (PBC) of uncertain/controversial clinical or prognostic significance. To study the genetic determinants associated with sp100 and gp210 autoantibody subphenotypes, we performed a genome-wide association analysis of 930 PBC cases based on their autoantibody status, followed by a replication study in 1,252 PBC cases. We confirmed single-nucleotide polymorphisms rs492899 (P = 3.27 × 10-22 ; odds ratio [OR], 2.90; 95% confidence interval [CI], 2.34-3.66) and rs1794280 (P = 5.78 × 10-28 ; OR, 3.89; 95% CI, 3.05-4.96) in the human major histocompatibility complex (MHC) region associated with the sp100 autoantibody. However, no genetic variant was identified as being associated with the gp210 autoantibody. To further define specific classical human leukocyte antigen (HLA) alleles or amino acids associated with the sp100 autoantibody, we imputed 922 PBC cases (211 anti-sp100-positive versus 711 negative cases) using a Han Chinese MHC reference database. Conditional analysis identified that HLA-DRβ1-Asn77/Arg74, DRβ1-Ser37, and DPβ1-Lys65 were major determinants for sp100 production. For the classical HLA alleles, the strongest association was with DRB1*03:01 (P = 1.51 × 10-9 ; OR, 2.97; 95% CI, 2.06-4.29). Regression analysis with classical HLA alleles identified DRB1*03:01, DRB1*15:01, DRB1*01, and DPB1*03:01 alleles can explain most of the HLA association with sp100 autoantibody. Conclusion: This study indicated significant genetic predisposition to the sp100 autoantibody, but not the gp210 autoantibody, subphenotype in PBC patients. Additional studies will be necessary to determine if these findings have clinical significance to PBC pathogenesis and/or therapeutics.

Genetic variation at the CD28 locus and its impact on expansion of pro-inflammatory CD28 negative T cells in healthy individuals

The CD28 locus is associated with susceptibility to a variety of autoimmune and immune-mediated inflammatory diseases including primary sclerosing cholangitis (PSC). Previously, we linked the CD28 pathway in PSC disease pathology and found that vitamin D could maintain CD28 expression. Here, we assessed whether the PSC-associated CD28 risk variant A (rs7426056) affects CD28 expression and T cell function in healthy individuals (n = 14 AA, n = 14 AG, n = 14 GG). Homozygotes for the PSC disease risk allele (AA) showed significantly lower CD28 mRNA expression ex-vivo than either GG or AG (p < 0.001) in total peripheral blood mononuclear cells. However, the CD28 risk variant alone was not sufficient to explain CD28 protein loss on CD4⁺ T cells. All genotypes responded equally to vitamin D as indicated by induction of a regulatory phenotype and an increased anti-inflammatory/pro-inflammatory cytokine ratio. A genotypic effect on response to TNFα stimuli was detected, which was inhibited by vitamin D. Together our results show: (a) an altered gene expression in carriers of the susceptible CD28 variant, (b) no differences in protein levels on CD4⁺ T cells, and (c) a protective effect of the variant upon CD28 protein loss on CD4⁺ T cells under inflammatory conditions.

A genome-wide association study identifies six novel risk loci for primary biliary cholangitis

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.

Immune-Array Analysis in Sporadic Inclusion Body Myositis Reveals HLA-DRB1 Amino Acid Heterogeneity Across the Myositis Spectrum

OBJECTIVE

Inclusion body myositis (IBM) is characterized by a combination of inflammatory and degenerative changes affecting muscle. While the primary cause of IBM is unknown, genetic factors may influence disease susceptibility. To determine genetic factors contributing to the etiology of IBM, we conducted the largest genetic association study of the disease to date, investigating immune-related genes using the Immunochip.

METHODS

A total of 252 Caucasian patients with IBM were recruited from 11 countries through the Myositis Genetics Consortium and compared with 1,008 ethnically matched controls. Classic HLA alleles and amino acids were imputed using SNP2HLA.

RESULTS

The HLA region was confirmed as the most strongly associated region in IBM (P = 3.58 × 10^-33 ). HLA imputation identified 3 independent associations (with HLA-DRB1*03:01, DRB1*01:01, and DRB1*13:01), although the strongest association was with amino acid positions 26 and 11 of the HLA-DRB1 molecule. No association with anti-cytosolic 5'-nucleotidase 1A-positive status was found independent of HLA-DRB1*03:01. There was no association of HLA genotypes with age at onset of IBM. Three non-HLA regions reached suggestive significance, including the chromosome 3 p21.31 region, an established risk locus for autoimmune disease, where a frameshift mutation in CCR5 is thought to be the causal variant.

CONCLUSION

This is the largest, most comprehensive genetic association study to date in IBM. The data confirm that HLA is the most strongly associated region and identifies novel amino acid associations that may explain the risk in this locus. These amino acid associations differentiate IBM from polymyositis and dermatomyositis and may determine properties of the peptide-binding groove, allowing it to preferentially bind autoantigenic peptides. A novel suggestive association within the chromosome 3 p21.31 region suggests a role for CCR5.

The OncoArray Consortium: A Network for Understanding the Genetic Architecture of Common Cancers

BACKGROUND

Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits.

METHODS

The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background.

RESULTS

The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis.

CONCLUSIONS

Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures.

IMPACT

Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126-35. ©2016 AACR.

FastPop: a rapid principal component derived method to infer intercontinental ancestry using genetic data

BACKGROUND

Identifying subpopulations within a study and inferring intercontinental ancestry of the samples are important steps in genome wide association studies. Two software packages are widely used in analysis of substructure: Structure and Eigenstrat. Structure assigns each individual to a population by using a Bayesian method with multiple tuning parameters. It requires considerable computational time when dealing with thousands of samples and lacks the ability to create scores that could be used as covariates. Eigenstrat uses a principal component analysis method to model all sources of sampling variation. However, it does not readily provide information directly relevant to ancestral origin; the eigenvectors generated by Eigenstrat are sample specific and thus cannot be generalized to other individuals.

RESULTS

We developed FastPop, an efficient R package that fills the gap between Structure and Eigenstrat. It can: 1, generate PCA scores that identify ancestral origins and can be used for multiple studies; 2, infer ancestry information for data arising from two or more intercontinental origins. We demonstrate the use of FastPop using 2318 SNP markers selected from the genome based on high variability among European, Asian and West African (African) populations. We conducted an analysis of 505 Hapmap samples with European, African or Asian ancestry along with 19661 additional samples of unknown ancestry. The results from FastPop are highly consistent with those obtained by Structure across the 19661 samples we studied. The correlations of the results between FastPop and Structure are 0.99, 0.97 and 0.99 for European, African and Asian ancestry scores, respectively. Compared with Structure, FastPop is more efficient as it finished ancestry inference for 19661 samples in 16 min compared with 21-24 h required by Structure. FastPop also provided scores based on SNP weights so the scores of reference population can be applied to other studies provided the same set of markers are used. We also present application of the method for studying four continental populations (European, Asian, African, and Native American).

CONCLUSIONS

We developed an algorithm that can infer ancestries on data involving two or more intercontinental origins. It is efficient for analyzing large datasets. Additionally the PCA derived scores can be applied to multiple data sets to ensure the same ancestry analysis is applied to all studies.

A Common Variant in CLDN14 is Associated with Primary Biliary Cirrhosis and Bone Mineral Density

Primary biliary cirrhosis (PBC), a chronic autoimmune liver disease, has been associated with increased incidence of osteoporosis. Intriguingly, two PBC susceptibility loci identified through genome-wide association studies are also involved in bone mineral density (BMD). These observations led us to investigate the genetic variants shared between PBC and BMD. We evaluated 72 genome-wide significant BMD SNPs for association with PBC using two European GWAS data sets (n = 8392), with replication of significant findings in a Chinese cohort (685 cases, 1152 controls). Our analysis identified a novel variant in the intron of the CLDN14 gene (rs170183, P_fdr = 0.015) after multiple testing correction. The three associated variants were followed-up in the Chinese cohort; one SNP rs170183 demonstrated consistent evidence of association in diverse ethnic populations (P_combined = 2.43 × 10⁻⁵). Notably, expression quantitative trait loci (eQTL) data revealed that rs170183 was correlated with a decline in CLDN14 expression in both lymphoblastoid cell lines and T cells (P_adj = 0.003 and 0.016, respectively). In conclusion, our study identified a novel PBC susceptibility variant that has been shown to be strongly associated with BMD, highlighting the potential of pleiotropy to improve gene discovery.

Multiple genetic variants associated with primary biliary cirrhosis in a Han Chinese population

Multiple genome-wide association studies of primary biliary cirrhosis (PBC) in both European and Japanese ancestries have shown significant associations of many genetic loci contributing to the susceptibility to PBC. Major differences in susceptibility loci between these two population groups were observed. In this study, we examined whether the most significant loci observed in either European and/or Japanese cohorts are associated with PBC in a Han Chinese population. In 1070 PBC patients and 1198 controls, we observed highly significant associations at CD80 (rs2293370, P = 2.67 × 10(-8)) and TNFSF15 (rs4979462, P = 3.86 × 10(-8)) and significant associations at 17q12-21 (rs9303277), PDGFB (rs715505), NF-κB1 (rs7665090), IL12RB2 (rs11209050), and STAT4 (rs7574865; all corrected P values <0.01). However, no association was observed for POU2AF1 (rs4938534), IL12A (rs485499 and rs2366408), IL7R (rs6897932), CXCR5 (rs715412), SOCS1 (rs725613), and TNFRSF1A (rs1800693). STAT4 (rs7574865) was strongly associated after additional control samples were analyzed. Our study is the first large-scale genetic analysis in a Han Chinese PBC cohort. These results do not only reflect that Han Chinese PBC patients share common genetic susceptibility genes with both their Japanese and European counterparts but also suggest a distinctly different genetic susceptibility profile.

Genome-Wide Association Study of Late-Onset Myasthenia Gravis: Confirmation of TNFRSF11A and Identification of ZBTB10 and Three Distinct HLA Associations

To investigate the genetics of late-onset myasthenia gravis (LOMG), we conducted a genome-wide association study imputation of>6 million single nucleotide polymorphisms (SNPs) in 532 LOMG cases (anti-acetylcholine receptor [AChR] antibody positive; onset age≥50 years) and 2,128 controls matched for sex and population substructure. The data confirm reported TNFRSF11A associations (rs4574025, P = 3.9 × 10^-7, odds ratio [OR] 1.42) and identify a novel candidate gene, ZBTB10, achieving genome-wide significance (rs6998967, P = 8.9 × 10^-10, OR 0.53). Several other SNPs showed suggestive significance including rs2476601 (P = 6.5 × 10^-6, OR 1.62) encoding the PTPN22 R620W variant noted in early-onset myasthenia gravis (EOMG) and other autoimmune diseases. In contrast, EOMG-associated SNPs in TNIP1 showed no association in LOMG, nor did other loci suggested for EOMG. Many SNPs within the major histocompatibility complex (MHC) region showed strong associations in LOMG, but with smaller effect sizes than in EOMG (highest OR ~2 versus ~6 in EOMG). Moreover, the strongest associations were in opposite directions from EOMG, including an OR of 0.54 for DQA1*05:01 in LOMG (P = 5.9 × 10^-12) versus 2.82 in EOMG (P = 3.86 × 10^-45). Association and conditioning studies for the MHC region showed three distinct and largely independent association peaks for LOMG corresponding to (a) MHC class II (highest attenuation when conditioning on DQA1), (b) HLA-A and (c) MHC class III SNPs. Conditioning studies of human leukocyte antigen (HLA) amino acid residues also suggest potential functional correlates. Together, these findings emphasize the value of subgrouping myasthenia gravis patients for clinical and basic investigations and imply distinct predisposing mechanisms in LOMG.

The genetics of human autoimmune disease: A perspective on progress in the field and future directions

Progress in defining the genetics of autoimmune disease has been dramatically enhanced by large scale genetic studies. Genome-wide approaches, examining hundreds or for some diseases thousands of cases and controls, have been implemented using high throughput genotyping and appropriate algorithms to provide a wealth of data over the last decade. These studies have identified hundreds of non-HLA loci as well as further defining HLA variations that predispose to different autoimmune diseases. These studies to identify genetic risk loci are also complemented by progress in gene expression studies including definition of expression quantitative trait loci (eQTL), various alterations in chromatin structure including histone marks, DNase I sensitivity, repressed chromatin regions as well as transcript factor binding sites. Integration of this information can partially explain why particular variations can alter proclivity to autoimmune phenotypes. Despite our incomplete knowledge base with only partial definition of hereditary factors and possible functional connections, this progress has and will continue to facilitate a better understanding of critical pathways and critical changes in immunoregulation. Advances in defining and understanding functional variants potentially can lead to both novel therapeutics and personalized medicine in which therapeutic approaches are chosen based on particular molecular phenotypes and genomic alterations.

Genome-Wide Association and Trans-ethnic Meta-Analysis for Advanced Diabetic Kidney Disease: Family Investigation of Nephropathy and Diabetes (FIND)

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.

Type 2 diabetes is the most common cause of severe kidney disease worldwide and diabetic kidney disease (DKD) associates with premature death. Individuals of non-European ancestry have the highest burden of type 2 DKD; hence understanding the causes of DKD remains critical to reducing health disparities. Family studies demonstrate that genes regulate the onset and progression of DKD; however, identifying these genes has proven to be challenging. The Family Investigation of Diabetes and Nephropathy consortium (FIND) recruited a large multi-ethnic collection of individuals with type 2 diabetes with and without kidney disease in order to detect genes associated with DKD. FIND discovered and replicated a DKD-associated genetic locus on human chromosome 6q25.2 (rs955333) between the SCAF8 and CNKSR genes. Findings were supported by significantly different expression of genes in this region from kidney tissue of subjects with, versus without DKD. The present findings identify a novel kidney disease susceptibility locus in individuals with type 2 diabetes which is consistent across subjects of differing ancestries. In addition, FIND results provide a rich catalogue of genetic variation in DKD patients for future research on the genetic architecture regulating this common and devastating disease.

TYK2 protein-coding variants protect against rheumatoid arthritis and autoimmunity, with no evidence of major pleiotropic effects on non-autoimmune complex traits

Despite the success of genome-wide association studies (GWAS) in detecting a large number of loci for complex phenotypes such as rheumatoid arthritis (RA) susceptibility, the lack of information on the causal genes leaves important challenges to interpret GWAS results in the context of the disease biology. Here, we genetically fine-map the RA risk locus at 19p13 to define causal variants, and explore the pleiotropic effects of these same variants in other complex traits. First, we combined Immunochip dense genotyping (n = 23,092 case/control samples), Exomechip genotyping (n = 18,409 case/control samples) and targeted exon-sequencing (n = 2,236 case/controls samples) to demonstrate that three protein-coding variants in TYK2 (tyrosine kinase 2) independently protect against RA: P1104A (rs34536443, OR = 0.66, P = 2.3x10^-21), A928V (rs35018800, OR = 0.53, P = 1.2x10^-9), and I684S (rs12720356, OR = 0.86, P = 4.6x10^-7). Second, we show that the same three TYK2 variants protect against systemic lupus erythematosus (SLE, P_omnibus = 6x10^-18), and provide suggestive evidence that two of the TYK2 variants (P1104A and A928V) may also protect against inflammatory bowel disease (IBD; P_omnibus = 0.005). Finally, in a phenome-wide association study (PheWAS) assessing >500 phenotypes using electronic medical records (EMR) in >29,000 subjects, we found no convincing evidence for association of P1104A and A928V with complex phenotypes other than autoimmune diseases such as RA, SLE and IBD. Together, our results demonstrate the role of TYK2 in the pathogenesis of RA, SLE and IBD, and provide supporting evidence for TYK2 as a promising drug target for the treatment of autoimmune diseases.

Relationship of pain and ancestry in African American women

BACKGROUND

African Americans are reported to be more sensitive to pain than European Americans. Pain sensitivity has been shown to be genetically linked in animal models and is likely to be in humans.

METHODS

Exactly, 11,239 self-identified African American post-menopausal women enrolled in the Women's Health Initiative had percentage African ancestry determined by ancestry informative markers, "Pain Construct" measurements and covariate information. They answered five questions about specific types and location of pain, such as joint, neck, low back, headache and urinary. They also answered two questions which were used to derive a "Pain Construct", a measure of general pain scored on a scale of 1-100. Associations were tested in linear regression models adjusting for age, self-reported medical conditions, neighbourhood socio-economic status, education and depression.

RESULTS

In the unadjusted model of the five specific types of pain measures, greater pain perception was associated with a higher proportion of African ancestry. However, some of the specific types of pain measures were no longer associated with African ancestry after adjustment for other study covariates. The Pain Construct was statistically significantly associated with African ancestry in both the unadjusted [β = -0.132, 95% confidence interval (CI) = -099 to -0.164; r = -0.075, 95% CI -0.056 to -0.093] and the adjusted models (β = -0.069 95% CI = -0.04 to -0.10).

CONCLUSIONS

Greater African ancestry was associated with higher levels of self-reported pain, although this accounted for only a minor fraction of the overall variation in the Pain Construct.

Genome-wide analysis of DNA methylation, copy number variation, and gene expression in monozygotic twins discordant for primary biliary cirrhosis

Primary biliary cirrhosis (PBC) is an uncommon autoimmune disease with a homogeneous clinical phenotype that reflects incomplete disease concordance in monozygotic (MZ) twins. We have taken advantage of a unique collection consisting of genomic DNA and mRNA from peripheral blood cells of female MZ twins (n = 3 sets) and sisters of similar age (n = 8 pairs) discordant for disease. We performed a genome-wide study to investigate differences in (i) DNA methylation (using a custom tiled four-plex array containing tiled 50-mers 19,084 randomly chosen methylation sites), (ii) copy number variation (CNV) (with a chip including markers derived from the 1000 Genomes Project, all three HapMap phases, and recently published studies), and/or (iii) gene expression (by whole-genome expression arrays). Based on the results obtained from these three approaches we utilized quantitative PCR to compare the expression of candidate genes. Importantly, our data support consistent differences in discordant twins and siblings for the (i) methylation profiles of 60 gene regions, (ii) CNV of 10 genes, and (iii) the expression of 2 interferon-dependent genes. Quantitative PCR analysis showed that 17 of these genes are differentially expressed in discordant sibling pairs. In conclusion, we report that MZ twins and sisters discordant for PBC manifest particular epigenetic differences and highlight the value of the epigenetic study of twins.

Genome-wide analysis of DNA methylation, copy number variation, and gene expression in monozygotic twins discordant for primary biliary cirrhosis

Primary biliary cirrhosis (PBC) is an uncommon autoimmune disease with a homogeneous clinical phenotype that reflects incomplete disease concordance in monozygotic (MZ) twins. We have taken advantage of a unique collection consisting of genomic DNA and mRNA from peripheral blood cells of female MZ twins (n = 3 sets) and sisters of similar age (n = 8 pairs) discordant for disease. We performed a genome-wide study to investigate differences in (i) DNA methylation (using a custom tiled four-plex array containing tiled 50-mers 19,084 randomly chosen methylation sites), (ii) copy number variation (CNV) (with a chip including markers derived from the 1000 Genomes Project, all three HapMap phases, and recently published studies), and/or (iii) gene expression (by whole-genome expression arrays). Based on the results obtained from these three approaches we utilized quantitative PCR to compare the expression of candidate genes. Importantly, our data support consistent differences in discordant twins and siblings for the (i) methylation profiles of 60 gene regions, (ii) CNV of 10 genes, and (iii) the expression of 2 interferon-dependent genes. Quantitative PCR analysis showed that 17 of these genes are differentially expressed in discordant sibling pairs. In conclusion, we report that MZ twins and sisters discordant for PBC manifest particular epigenetic differences and highlight the value of the epigenetic study of twins.

Relationship between glaucoma and admixture in postmenopausal African American women

OBJECTIVE

To investigate the association between African admixture and glaucoma prevalence among African American women.

DESIGN, SETTING, PARTICIPANTS

Participants included 11616 African American women from the Women's Health Initiative Study (WHI) for whom admixture information was available and included 2548 who self-reported a diagnosis of glaucoma.

MAIN OUTCOME MEASURES

Glaucoma.

RESULTS

Significant association was observed between self-identified glaucoma status and admixture. However, this association was not significant in a model that included neighborhood socioeconomic status (NSES), hypertension, diabetes and body mass index (BMI). Self-identified glaucoma status was associated with diabetes that persisted after adjustment for admixture, NSES, hypertension, and BMI. Lower NSES was also associated with higher glaucoma risk but this association was marginal in the fully adjusted model and neither hypertension nor BMI showed association. When glaucoma status was limited to those reporting use or no use of appropriate ophthalmologic medication, no associations were observed in any of the models.

CONCLUSION

This study failed to find an independent association of glaucoma status and African admixture and these findings suggest that the higher frequency glaucoma in African Americans may be largely due to other factors.

Association study of nicotinic acetylcholine receptor genes identifies a novel lung cancer susceptibility locus near CHRNA1 in African-Americans

Studies in European and East Asian populations have identified lung cancer susceptibility loci in nicotinic acetylcholine receptor (nAChR) genes on chromosome 15q25.1 which also appear to influence smoking behaviors. We sought to determine if genetic variation in nAChR genes influences lung cancer susceptibly in African-Americans, and evaluated the association of these cancer susceptibility loci with smoking behavior. A total of 1308 African-Americans with lung cancer and 1241 African-American controls from three centers were genotyped for 378 single nucleotide polymorphisms (SNPs) spanning the sixteen human nAChR genes. Associations between SNPs and the risk of lung cancer were estimated using logistic regression, adjusted for relevant covariates. Seven SNPs in three nAChR genes were significantly associated with lung cancer at a strict Bonferroni-corrected level, including a novel association on chromosome 2 near the promoter of CHRNA1 (rs3755486: OR = 1.40, 95% CI = 1.18-1.67, P = 1.0 x 10-4). Association analysis of an additional 305 imputed SNPs on 2q31.1 supported this association. Publicly available expression data demonstrated that the rs3755486 risk allele correlates with increased CHRNA1 gene expression. Additional SNP associations were observed on 15q25.1 in genes previously associated with lung cancer, including a missense variant in CHRNA5 (rs16969968: OR = 1.60, 95% CI = 1.27-2.01, P = 5.9 x 10-5). Risk alleles on 15q25.1 also correlated with an increased number of cigarettes smoked per day among the controls. These findings identify a novel lung cancer risk locus on 2q31.1 which correlates with CHRNA1 expression and replicate previous associations on 15q25.1 in African-Americans.

Socioeconomic status and lung cancer: unraveling the contribution of genetic admixture

OBJECTIVES

We examined the relationship between genetic ancestry, socioeconomic status (SES), and lung cancer among African Americans and Latinos.

METHODS

We evaluated SES and genetic ancestry in a Northern California lung cancer case-control study (1998-2003) of African Americans and Latinos. Lung cancer case and control participants were frequency matched on age, gender, and race/ethnicity. We assessed case-control differences in individual admixture proportions using the 2-sample t test and analysis of covariance. Logistic regression models examined associations among genetic ancestry, socioeconomic characteristics, and lung cancer.

RESULTS

Decreased Amerindian ancestry was associated with higher education among Latino control participants and greater African ancestry was associated with decreased education among African lung cancer case participants. Education was associated with lung cancer among both Latinos and African Americans, independent of smoking, ancestry, age, and gender. Genetic ancestry was not associated with lung cancer among African Americans.

CONCLUSIONS

Findings suggest that socioeconomic factors may have a greater impact than genetic ancestry on lung cancer among African Americans. The genetic heterogeneity and recent dynamic migration and acculturation of Latinos complicate recruitment; thus, epidemiological analyses and findings should be interpreted cautiously.

Relationship between hysterectomy and admixture in African American women

OBJECTIVE

Most studies suggest that hysterectomies are more common in African American women than in other ethnic groups. To assess this ethnic surgical disparity in a novel way, our main goal was to determine whether admixture (the proportion of sub-Saharan African or European origin in individuals) is associated with hysterectomy frequency in African American women in the Women's Health Initiative.

STUDY DESIGN

In this retrospective study, we used ancestry informative single nucleotide polymorphisms to estimate admixture proportions in >10,000 African American women from the Women's Health Initiative. Logistic regression models were used to assess the association between admixture and self-reported history of hysterectomy with and without controls for relevant covariates. Multinomial logistic regression models were used to assess the association between admixture and self-reported age of hysterectomy. We also considered other potential risk factors (adiposity, hypertension, and education) for hysterectomy accounting for admixture.

RESULTS

African admixture was a strong risk factor after the adjustment for multiple covariates (odds ratio, 1.85; P < .0001). The admixture risk for hysterectomy was highest for those procedures that were performed in the 35-39 age range (odds ratio, 3.08; P < .0001) and least evident in oldest ages (≥45 years old). Our analyses also suggest that adiposity, hypertension, and education were associated independently with hysterectomy in this population group.

CONCLUSION

These results suggest that higher African admixture is associated with higher frequencies of hysterectomy and that genetic studies that specifically target African American women and diseases that are associated with hysterectomy may be especially useful in understanding the pathogenesis and underlying cause of this disparity in health outcome.

Association of DXA-derived bone mineral density and fat mass with African ancestry

CONTEXT

Both genes and environment have been implicated in determining the complex body composition phenotypes in individuals of European ancestry; however, few studies have been conducted in other race/ethnic groups.

OBJECTIVE

We conducted a genome-wide admixture mapping study in an attempt to localize novel genomic regions associated with genetic ancestry.

SETTING/PARTICIPANTS

We selected a sample of 842 African-American women from the Women's Health Initiative single nucleotide polymorphism (SNP) Health Association Resource for whom several dual-energy X-ray absorptiometry (DXA)-derived bone mineral density (BMD) and fat mass phenotypes were available.

METHODS

We derived both global and local ancestry estimates for each individual from Affymetrix 6.0 data and analyzed the correlation of DXA phenotypes with global African ancestry. For each phenotype, we examined the association of local genetic ancestry (number of African ancestral alleles at each marker) and each DXA phenotype at 570 282 markers across the genome in additive models with adjustment for important covariates.

RESULTS

We identified statistically significant correlations of whole-body fat mass, trunk fat mass, and all 6 measures of BMD with a proportion of African ancestry. Genome-wide (admixture) significance for femoral neck BMD was achieved across 2 regions ∼3.7 MB and 0.3 MB on chromosome 19q13; similarly, total hip and intertrochanter BMD were associated with local ancestry in these regions. Trunk fat was the most significant fat mass phenotype showing strong, but not genomewide significant associations on chromosome Xp22.

CONCLUSIONS

Our results suggest that genomic regions in postmenopausal African-American women contribute to variance in BMD and fat mass existence and warrant further study.

Pathway-based analysis of primary biliary cirrhosis genome-wide association studies

Genome-wide association studies (GWAS) have successfully identified several loci associated with primary biliary cirrhosis (PBC) risk. Pathway analysis complements conventional GWAS analysis. We applied the recently developed linear combination test for pathways to datasets drawn from independent PBC GWAS in Italian and Canadian subjects. Of the Kyoto Encyclopedia of Genes and Genomes and BioCarta pathways tested, 25 pathways in the Italian dataset (449 cases, 940 controls) and 26 pathways in the Canadian dataset (530 cases, 398 controls) were associated with PBC susceptibility (P<0.05). After correcting for multiple comparisons, only the eight most significant pathways in the Italian dataset had FDR <0.25 with tumor necrosis factor/stress-related signaling emerging as the top pathway (P=7.38 × 10⁻⁴, FDR=0.18). Two pathways, phosphatidylinositol signaling and hedgehog signaling, were replicated in both datasets (P<0.05), and subjected to two additional complementary pathway tests. Both pathway signals remained significant in the Italian dataset on modified gene set enrichment analysis (P<0.05). In both GWAS, variants nominally associated with PBC were significantly overrepresented in the phosphatidylinositol pathway (Fisher exact P<0.05). These results point to established and novel pathway-level associations with inherited predisposition to PBC that, on further independent replication and functional validation, may provide fresh insights into PBC etiology.

European genetic ancestry is associated with a decreased risk of lupus nephritis

OBJECTIVE

African Americans, East Asians, and Hispanics with systemic lupus erythematosus (SLE) are more likely to develop renal disease than are SLE patients of European descent. This study was undertaken to investigate whether European genetic ancestry protects against the development of lupus nephritis, with the aim of exploring the genetic and socioeconomic factors that might explain this effect.

METHODS

This was a cross-sectional study of SLE patients from a multiethnic case collection. Participants were genotyped for 126 single-nucleotide polymorphisms (SNPs) informative for ancestry. A subset of participants was also genotyped for 80 SNPs in 14 candidate genes for renal disease in SLE. Logistic regression was used to test the association between European ancestry and renal disease. Analyses were adjusted for continental ancestries, socioeconomic status (SES), and candidate genes.

RESULTS

Participants (n = 1,906) had, on average, 62.4% European, 15.8% African, 11.5% East Asian, 6.5% Amerindian, and 3.8% South Asian ancestry. Among the participants, 656 (34%) had renal disease. A 10% increase in the proportion of European ancestry estimated in each participant was associated with a 15% reduction in the odds of having renal disease, after adjustment for disease duration and sex (odds ratio 0.85, 95% confidence interval 0.82-0.87; P = 1.9 × 10(-30) ). Adjustment for other genetic ancestries, measures of SES, or SNPs in the genes most associated with renal disease (IRF5 [rs4728142], BLK [rs2736340], STAT4 [rs3024912], and HLA-DRB1*0301 and DRB1*1501) did not substantively alter this relationship.

CONCLUSION

European ancestry is protective against the development of renal disease in SLE, an effect that is independent of other genetic ancestries, candidate risk alleles, and socioeconomic factors.

Fine-mapping of the 5p15.33, 6p22.1-p21.31, and 15q25.1 regions identifies functional and histology-specific lung cancer susceptibility loci in African-Americans

BACKGROUND

Genome-wide association studies of European and East Asian populations have identified lung cancer susceptibility loci on chromosomes 5p15.33, 6p22.1-p21.31, and 15q25.1. We investigated whether these regions contain lung cancer susceptibly loci in African-Americans and refined previous association signals by using the reduced linkage disequilibrium observed in African-Americans.

METHODS

1,308 African-American cases and 1,241 African-American controls from 3 centers were genotyped for 760 single-nucleotide polymorphisms (SNP) spanning 3 regions, and additional SNP imputation was carried out. Associations between polymorphisms and lung cancer risk were estimated using logistic regression, stratified by tumor histology where appropriate.

RESULTS

The strongest associations were observed on 15q25.1 in/near CHRNA5, including a missense substitution [rs16969968: OR, 1.57; 95% confidence interval (CI), 1.25-1.97; P, 1.1 × 10(-4)) and variants in the 5'-UTR. Associations on 6p22.1-p21.31 were histology specific and included a missense variant in BAT2 associated with squamous cell carcinoma (rs2736158: OR, 0.64; 95% CI, 0.48-0.85; P, 1.82 × 10(-3)). Associations on 5p15.33 were detected near TERT, the strongest of which was rs2735940 (OR, 0.82; 95% CI, 0.73-0.93; P, 1.1 × 10(-3)). This association was stronger among cases with adenocarcinoma (OR, 0.75; 95% CI, 0.65-0.86; P, 8.1 × 10(-5)).

CONCLUSIONS

Polymorphisms in 5p15.33, 6p22.1-p21.31, and 15q25.1 are associated with lung cancer in African-Americans. Variants on 5p15.33 are stronger risk factors for adenocarcinoma and variants on 6p21.33 associated only with squamous cell carcinoma.

IMPACT

Results implicate the BAT2, TERT, and CHRNA5 genes in the pathogenesis of specific lung cancer histologies.

Unraveling multiple MHC gene associations with systemic lupus erythematosus: model choice indicates a role for HLA alleles and non-HLA genes in Europeans

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1(∗)03:01, HLA-DRB1(∗)08:01, and HLA-DQA1(∗)01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).

Risk for myasthenia gravis maps to a (151) Pro→Ala change in TNIP1 and to human leukocyte antigen-B*08

OBJECTIVE

The objective of this study is to comprehensively define the genetic basis of early onset myasthenia gravis (EOMG).

METHODS

We have carried out a 2-stage genome-wide association study on a total of 649 North European EOMG patients. Cases were matched 1:4 with controls of European ancestry. We performed imputation and conditional analyses across the major histocompatibility complex, as well as in the top regions of association outside the human leukocyte antigen (HLA) region.

RESULTS

We observed the strongest association in the HLA class I region at rs7750641 (p = 1.2 × 10(-92) ; odds ratio [OR], 6.25). By imputation and conditional analyses, HLA-B*08 proves to be the major associated allele (p = 2.87 × 10(-113) ; OR, 6.41). In addition to the expected association with PTPN22 (rs2476601; OR, 1.71; p = 8.2 × 10(-10) ), an imputed coding variant (rs2233290) at position 151 (Pro→Ala) in the TNFAIP3-interacting protein 1, TNIP1, confers even stronger risk than PTPN22 (OR, 1.91; p = 3.2 × 10(-10) ).

INTERPRETATION

The association at TNIP1 in EOMG implies disease mechanisms involving ubiquitin-dependent dysregulation of NF-κB signaling. The localization of the major HLA signal to the HLA-B*08 allele suggests that CD8(+) T cells may play a key role in disease initiation or pathogenesis.

Immunochip analyses identify a novel risk locus for primary biliary cirrhosis at 13q14, multiple independent associations at four established risk loci and epistasis between 1p31 and 7q32 risk variants

To further characterize the genetic basis of primary biliary cirrhosis (PBC), we genotyped 2426 PBC patients and 5731 unaffected controls from three independent cohorts using a single nucleotide polymorphism (SNP) array (Immunochip) enriched for autoimmune disease risk loci. Meta-analysis of the genotype data sets identified a novel disease-associated locus near the TNFSF11 gene at 13q14, provided evidence for association at six additional immune-related loci not previously implicated in PBC and confirmed associations at 19 of 22 established risk loci. Results of conditional analyses also provided evidence for multiple independent association signals at four risk loci, with haplotype analyses suggesting independent SNP effects at the 2q32 and 16p13 loci, but complex haplotype driven effects at the 3q25 and 6p21 loci. By imputing classical HLA alleles from this data set, four class II alleles independently contributing to the association signal from this region were identified. Imputation of genotypes at the non-HLA loci also provided additional associations, but none with stronger effects than the genotyped variants. An epistatic interaction between the IL12RB2 risk locus at 1p31and the IRF5 risk locus at 7q32 was also identified and suggests a complementary effect of these loci in predisposing to disease. These data expand the repertoire of genes with potential roles in PBC pathogenesis that need to be explored by follow-up biological studies.

Relationship between gallbladder surgery and ethnic admixture in African American and Hispanic American women

OBJECTIVES

The objective of this study was to investigate whether differences in admixture in African American and Hispanic American adult women are associated with differences in gallbladder surgery.

METHODS

Gallbladder surgery history on entry to the Women's Health Initiative's (WHI) study was used as a dichotomous outcome measure for this study. The proportion of European, sub-Saharan African, and Amerindian (AMI) admixture was estimated for 10,841 African American and 4,620 Hispanic American women in WHI using 92 ancestry informative markers. Logistic regression analyses assessed the relationship between admixture and gallbladder surgery in WHI women (enrollment at ages >50, mean age 61 years) with or without adjusting for multiple covariates, including measures of adiposity, parity, alcohol use, and education.

RESULTS

There was a significant positive association between AMI admixture and the frequency of gallbladder surgery in Hispanic Americans. The odds ratio (OR) and 95% confidence intervals (CIs) for AMI admixture group was OR=2.97, CI=2.01-4.38, P<10(-4). Although there were strong positive associations with parity and adiposity, and negative associations with alcohol consumption and education, accounting for these covariates did not remove the admixture association (OR=2.46, CI=1.62-3.73). In contrast, the effect of African admixture was nearly indistinguishable from that of the European admixture, both of which were protective in the Hispanic American group, and African admixture had a marginal association with decreased gallbladder surgery in the African American group. Measures of adiposity were associated with increased risk for gallbladder surgery and remained significant after accounting for admixture and each of the other covariates. Education level and alcohol use were associated with decreased risk for gallbladder disease. However, after accounting for the other covariates these variably remained significant.

CONCLUSIONS

AMI admixture is strongly associated with gallbladder surgery in women, even after adjustment for selected risk factors for cholelithiasis. Additional studies to ascertain the specific genetic risk factors underlying these associations are warranted.

Classical HLA-DRB1 and DPB1 alleles account for HLA associations with primary biliary cirrhosis

Susceptibility to primary biliary cirrhosis (PBC) is strongly associated with HLA region polymorphisms. To determine if associations can be explained by classical HLA determinants we studied Italian 676 cases and 1440 controls with genotyped with dense single nucleotide polymorphisms (SNPs) for which classical HLA alleles and amino acids were imputed. Although previous genome-wide association studies and our results show stronger SNP associations near DQB1, we demonstrate that the HLA signals can be attributed to classical DRB1 and DPB1 genes. Strong support for the predominant role of DRB1 is provided by our conditional analyses. We also demonstrate an independent association of DPB1. Specific HLA-DRB1 genes (*08, *11 and *14) account for most of the DRB1 association signal. Consistent with previous studies, DRB1*08 (p = 1.59 × 10⁻¹¹) was the strongest predisposing allele where as DRB1*11 (p = 1.42 × 10⁻¹⁰) was protective. Additionally DRB1*14 and the DPB1 association (DPB1*03:01) (p = 9.18 × 10⁻⁷) were predisposing risk alleles. No signal was observed in the HLA class 1 or class 3 regions. These findings better define the association of PBC with HLA and specifically support the role of classical HLA-DRB1 and DPB1 genes and alleles in susceptibility to PBC.

Relationship between diabetes risk and admixture in postmenopausal African-American and Hispanic-American women

AIMS/HYPOTHESIS

Type 2 diabetes is more prevalent in African-Americans (AFAs) and Hispanic-Americans (HAs) than in European-Americans. We assessed whether continental admixture was correlated with diabetes risk in these high-risk groups.

METHODS

We estimated the proportion of sub-Saharan African (AFR), Amerindian (AMI) and European admixture using 92 ancestry-informative marker genotypes in 16,476 AFA and HA women from the Women's Health Initiative. Cox regression models were used to examine the association between admixture and diabetes risk, with and without accounting for socioeconomic status (SES) and adiposity measurements.

RESULTS

AFR admixture was significantly associated with diabetes risk in AFA women when adjusting for entry age, neighbourhood SES and BMI or waist/hip ratio (WHR) (all p < 0.0001). In HA women, AMI admixture had significant associations with diabetes risk that remained significant after adjustment for SES and BMI (all p < 0.0005). In both AFAs and HAs, SES showed significant negative associations while BMI or WHR had significant positive associations with diabetes risk, with and without adjustment for genetic admixture.

CONCLUSIONS/INTERPRETATION

In AFAs, admixture, SES and BMI/WHR each independently contribute to diabetes risk after accounting for each of the other factors; in HAs, admixture, SES and BMI each independently contribute to diabetes risk after accounting for each of the other factors, whereas admixture is not significantly associated with diabetes risk after accounting for SES and WHR. The findings emphasise the importance of considering both genetic and environmental causes in the aetiology of type 2 diabetes.

Herpesvirus saimiri encodes a new cytokine, IL-17, which binds to a novel cytokine receptor

Herpesvirus Saimiri gene 13 (HVS13) exhibits 57% identity with the predicted sequence of a T cell-derived molecule termed CTLA8. Recombinant HVS13 and CTLA8 stimulate transcriptional factor NF-kappaB activity and Interleukin-6 (IL-6) secretion in fibroblasts, and costimulate T cell proliferation. An HVS13.Fc fusion protein was used to isolate a cDNA encoding a novel receptor that also binds CTLA8. This receptor is unrelated to previously identified cytokine receptor families. A recombinant soluble receptor inhibited T cell proliferation and IL-2 production induced by PHA, concanavalin A (conA), and anti-TCR MAb. These results define CTLA8 and HVS13 as novel cytokines that bind to a novel cytokine receptor. We propose to call these molecules IL-17, vIL-17, and IL-17R, respectively.

Genetic variants of TSLP and asthma in an admixed urban population

BACKGROUND

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.

OBJECTIVES

To determine whether genetic variants in TSLP are associated with asthma in an urban admixed population. METHODOLOGY AND MAIN RESULTS: Ten tag-SNPs in the TSLP gene were analyzed for association with asthma using 387 clinically diagnosed asthmatic cases and 212 healthy controls from an urban admixed population. One SNP (rs1898671) showed nominally significant association with asthma (odds ratio (OR) = 1.50; 95% confidence interval (95% CI): 1.09-2.05, p = 0.01) after adjusting for age, BMI, income, education and population stratification. Association results were consistent using two different approaches to adjust for population stratification. When stratified by smoking status, the same SNP showed a significantly increased risk associated with asthma in ex-smokers (OR = 2.00, 95% CI: 1.04-3.83, p = 0.04) but not significant in never-smokers (OR = 1.34; 95% CI: 0.93-1.94, p = 0.11). Haplotype-specific score test indicated that an elevated risk for asthma was associated with a specific haplotype of TSLP involving SNP rs1898671 (OR = 1.58, 95% CI: 1.10-2.27, p = 0.01). Association of this SNP with asthma was confirmed in an independent large population-based cohort consortium study (OR = 1.15, 95% CI: 1.07-1.23, p = 0.0003) and the results stratified by smoking status were also validated (ex-smokers: OR = 1.21, 95% CI: 1.08-1.34, p = 0.003; never-smokers: OR = 1.06, 95% CI: 0.94-1.17, p = 0.33).

CONCLUSIONS

Genetic variants in TSLP may contribute to asthma susceptibility in admixed urban populations with a gene and environment interaction.

Genome-wide association study identifies novel loci predisposing to cutaneous melanoma

We performed a multistage genome-wide association study of melanoma. In a discovery cohort of 1804 melanoma cases and 1026 controls, we identified loci at chromosomes 15q13.1 (HERC2/OCA2 region) and 16q24.3 (MC1R) regions that reached genome-wide significance within this study and also found strong evidence for genetic effects on susceptibility to melanoma from markers on chromosome 9p21.3 in the p16/ARF region and on chromosome 1q21.3 (ARNT/LASS2/ANXA9 region). The most significant single-nucleotide polymorphisms (SNPs) in the 15q13.1 locus (rs1129038 and rs12913832) lie within a genomic region that has profound effects on eye and skin color; notably, 50% of variability in eye color is associated with variation in the SNP rs12913832. Because eye and skin colors vary across European populations, we further evaluated the associations of the significant SNPs after carefully adjusting for European substructure. We also evaluated the top 10 most significant SNPs by using data from three other genome-wide scans. Additional in silico data provided replication of the findings from the most significant region on chromosome 1q21.3 rs7412746 (P = 6 × 10(-10)). Together, these data identified several candidate genes for additional studies to identify causal variants predisposing to increased risk for developing melanoma.

Stronger bone correlates with African admixture in African-American women

Osteoporosis affects all populations, but the risk for low bone density and fracture varies greatly by self-reported race and ethnicity. In this article, the relationship between measured percent African admixture and hip structural geometry, estimated from the hip structural analysis (HSA) program, was examined in a subcohort of the Women's Health Initiative (WHI). The study sample included 793 African-American women and 8559 non-Hispanic white women. All the participants were postmenopausal, between the ages of 50 and 79 years, at the time of recruitment and were followed for up to 9 years. Bone density and hip geometry were assessed using dual-energy X-ray absorptiometry. African admixture was measured for African Americans using genetic ancestry informative markers. Multiple regression and mixed-effects models were used for cross-sectional and longitudinal analyses, respectively. Covariates assessed from questionnaires and physical measurements were included in the analysis to control for possible confounding effects. The study results show significant correlations between percent of African admixture and HSA assessments. In comparison with non-Hispanic white women, significantly greater bone strength, as indicated by higher hip bone mineral density and stronger hip geometry, in women with higher African admixture was observed. However, women with higher percent African admixture had larger reductions in bone strength than non-Hispanic white women during the follow-up.

New approaches to disease mapping in admixed populations

Admixed populations such as African Americans and Hispanic Americans are often medically underserved and bear a disproportionately high burden of disease. Owing to the diversity of their genomes, these populations have both advantages and disadvantages for genetic studies of complex phenotypes. Advances in statistical methodologies that can infer genetic contributions from ancestral populations may yield new insights into the aetiology of disease and may contribute to the applicability of genomic medicine to these admixed population groups.

The rs4774 CIITA missense variant is associated with risk of systemic lupus erythematosus

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor required for human leukocyte antigens (HLA) class II MHC-restricted antigen presentation. MHC genes, including the HLA class II DRB1*03:01 allele, are strongly associated with systemic lupus erythematosus (SLE). Recently the rs4774 CIITA missense variant (+1632G/C) was reported to be associated with susceptibility to multiple sclerosis. In the current study, we investigated CIITA, DRB1*03:01 and risk of SLE using a multi-stage analysis. In stage 1, 9 CIITA variants were tested in 658 cases and 1363 controls (N=2021). In stage 2, rs4774 was tested in 684 cases and 2938 controls (N=3622). We also performed a meta-analysis of the pooled 1342 cases and 4301 controls (N=5643). In stage 1, rs4774(*)C was associated with SLE (odds ratio (OR)=1.24, 95% confidence interval (95% CI)=1.07-1.44, P=4.2 × 10(-3)). Similar results were observed in stage 2 (OR=1.16, 95% CI=1.02-1.33, P=8.5 × 10(-3)) and the meta-analysis of the combined data set (OR=1.20, 95% CI=1.09-1.33, P(meta)=2.5 × 10(-4)). In all three analyses, the strongest evidence for association between rs4774(*)C and SLE was present in individuals who carried at least one copy of DRB1*03:01 (P(meta)=1.9 × 10(-3)). Results support a role for CIITA in SLE, which appears to be stronger in the presence of DRB1*03:01.

Evidence for malaria selection of a CR1 haplotype in Sardinia

Complement receptor 1 (CR1) levels have been associated with malarial susceptibility and/or severity of the disease in different population groups, and CR1 is a receptor for Plasmodium falciparum. In this study, multiple CR1 single-nucleotide polymorphisms (SNPs) showed strong evidence of population differentiation between Sardinian and other European ethnic groups. Cross population algorithms comparing haplotype structure and differences in haplotype and allele frequency distribution provided additional support for natural selection of CR1 in Sardinia. The predominant Sardinian CR1 haplotype included SNPs that are associated with decreased CR1 levels in Europeans and other population groups. Previous studies have shown that the SNPs within the dominant Sardinian haplotype have a significantly higher frequency in a malaria endemic compared with non-endemic regions in India. Together with the historical evidence of the prevalence of malaria in Sardinia, these data support the role of malaria leading to positive selection of this CR1 haplotype in Sardinia.

Replicated association of 17q12-21 with susceptibility of primary biliary cirrhosis in a Japanese cohort

To examine the genetics of susceptibility to primary biliary cirrhosis (PBC), genome-wide association studies GWAS have been performed in patients of European ancestry and have shown the significant associations of IL12-related pathways, SPIB, IRF5-TNPO3, and 17q12-21. We tested whether these findings could be extended to a Japanese cohort, 303 Japanese PBC and 298 controls. We failed to detect significant associations at IL12A (rs574808, rs1075498) and IL12RB2 (rs3790567). There was no genetic variance at IRF5-TNPO3 (rs10488631) in Japanese. A single nucleotide polymorphism (SNP) at SPIB (rs3745516) reached nominal significance, but the corrected P value did not reach significance. For the 17q12-21 region, two SNPs had nominally significant associations [GSDMB (rs2305480, P = 0.022) and ZPBP2 (rs11557467, P = 0.021)] and we noted a significant P value at a SNP in IKZF3 (rs939327, P = 0.0024, P(c) = 0.017) after correction for multiple comparisons. Thus, these results indicate a haplotype on 17q12-21 with a similar association in Japanese and European PBC.

Relationship between adiposity and admixture in African-American and Hispanic-American women

Objective

To investigate whether differences in admixture in African American (AFA) and Hispanic American (HA) adult women are associated with adiposity and adipose distribution.

Design

The proportion of European, sub– Saharan African and Amerindian admixture was estimated for AFA and HA women in the Women's Heath Initiative using 92 ancestry informative markers. Analyses assessed the relationship between admixture and adiposity indices.

Subjects

11712 AFA and 5088 HA self– identified post– menopausal women.

Results

There was a significant positive association between body mass index (BMI) and African admixture when BMI was considered as a continuous variable, and age, education, physical activity, parity, family income and smoking were included covariates (p < 10^{− 4}). A dichotomous model (upper and lower BMI quartiles) showed that African admixture was associated with a high odds ratio [OR = 3.27 (for 100% admixture compared to 0% admixture), 95% confidence interval (CI) 2.08 – 5.15]. For HA there was no association between BMI and admixture. In contrast, when waist to hip ratio (WHR) was used as a measure of adipose distribution, there was no significant association between WHR and admixture in AFA but there was a strong association in HA (p<10^{− 4}; OR Amerindian admixture = 5.93, CI = 3.52 – 9.97).

Conclusion

These studies show that 1) African admixture is associated with BMI in AFA women; 2) Amerindian admixture is associated with WHR but not BMI in HA women; and 3) it may be important to consider different measurements of adiposity and adipose distribution in different ethnic population groups.

ALOX5 gene variants affect eicosanoid production and response to fish oil supplementation

The objective of this study was to determine whether 5-lipoxygenase (ALOX5) gene variants associated with cardiovascular disease affect eicosanoid production by monocytes. The study was a randomized, double-masked, parallel intervention trial with fish oil (5.0 g of fish oil daily, containing 2.0 g of eicosapentaenoic acid [EPA] and 1.0 g of docosahexaenoic acid [DHA]) or placebo oil (5.0 g of corn/soy mixture). A total of 116 subjects (68% female, 20-59 years old) of African American ancestry enrolled, and 98 subjects completed the study. Neither ALOX5 protein nor arachidonic acid-derived LTB4, LTD4, and LTE4 varied by genotype, but 5-hydroxyeicosatetraenoate (5-HETE), 6-trans-LTB4, 5-oxo-ETE, 15-HETE, and 5,15-diHETE levels were higher in subjects homozygous for the ALOX5 promoter allele containing five Sp1 element tandem repeats ("55" genotype) than in subjects with one deletion (d) (three or four repeats) and one common ("d5" genotype) allele or with two deletion ("dd") alleles. The EPA-derived metabolites 5-HEPE and 15-HEPE and the DHA-derived metabolite 17-HDoHE had similar associations with genotype and increased with supplementation; 5-HEPE and 15-HEPE increased, and 5-oxo-ETE decreased to a greater degree in the 55 than in the other genotypes. This differential eicosanoid response is consistent with the previously observed interaction of these variants with dietary intake of omega-3 fatty acids in predicting cardiovascular disease risk.