Genome Wide Analysis of Narcolepsy in China Implicates Novel Immune Loci and Reveals Changes in Association Prior to Versus After the 2009 H1N1 Influenza Pandemic

Previous studies in narcolepsy, an autoimmune disorder affecting hypocretin (orexin) neurons and recently associated with H1N1 influenza, have demonstrated significant associations with five loci. Using a well-characterized Chinese cohort, we refined known associations in TRA@ and P2RY11-DNMT1 and identified new associations in the TCR beta (TRB@; rs9648789 max P = 3.7×10⁻⁹ OR 0.77), ZNF365 (rs10995245 max P = 1.2×10⁻¹¹ OR 1.23), and IL10RB-IFNAR1 loci (rs2252931 max P = 2.2×10⁻⁹ OR 0.75). Variants in the Human Leukocyte Antigen (HLA)- DQ region were associated with age of onset (rs7744020 P = 7.9×10⁻⁹ beta −1.9 years) and varied significantly among cases with onset after the 2009 H1N1 influenza pandemic compared to previous years (rs9271117 P = 7.8×10⁻¹⁰ OR 0.57). These reflected an association of DQB1*03:01 with earlier onset and decreased DQB1*06:02 homozygosity following 2009. Our results illustrate how genetic association can change in the presence of new environmental challenges and suggest that the monitoring of genetic architecture over time may help reveal the appearance of novel triggers for autoimmune diseases.     

SORL1 Is Genetically Associated with Late-Onset Alzheimer’s Disease in Japanese,Koreans and Caucasians

To discover susceptibility genes of late-onset Alzheimer’s disease (LOAD), we conducted a 3-stage genome-wide association study (GWAS) using three populations: Japanese from the Japanese Genetic Consortium for Alzheimer Disease (JGSCAD), Koreans, and Caucasians from the Alzheimer Disease Genetic Consortium (ADGC). In Stage 1, we evaluated data for 5,877,918 genotyped and imputed SNPs in Japanese cases (n = 1,008) and controls (n = 1,016). Genome-wide significance was observed with 12 SNPs in the APOE region. Seven SNPs from other distinct regions with p-values <2×10⁻⁵ were genotyped in a second Japanese sample (885 cases, 985 controls), and evidence of association was confirmed for one SORL1 SNP (rs3781834, P = 7.33×10⁻⁷ in the combined sample). Subsequent analysis combining results for several SORL1 SNPs in the Japanese, Korean (339 cases, 1,129 controls) and Caucasians (11,840 AD cases, 10,931 controls) revealed genome wide significance with rs11218343 (P = 1.77×10⁻⁹) and rs3781834 (P = 1.04×10⁻⁸). SNPs in previously established AD loci in Caucasians showed strong evidence of association in Japanese including rs3851179 near PICALM (P = 1.71×10⁻⁵) and rs744373 near BIN1 (P = 1.39×10⁻⁴). The associated allele for each of these SNPs was the same as in Caucasians. These data demonstrate for the first time genome-wide significance of LOAD with SORL1 and confirm the role of other known loci for LOAD in Japanese. Our study highlights the importance of examining associations in multiple ethnic populations.     

Genome-Wide Association Study Link Novel Loci to Endometriosis

Endometriosis is a common gynecological condition with complex etiology defined by the presence of endometrial glands and stroma outside the womb. Endometriosis is a common cause of both cyclic and chronic pelvic pain, reduced fertility, and reduced quality-of-life. Diagnosis and treatment of endometriosis is, on average, delayed by 7–10 years from the onset of symptoms. Absence of a timely and non-invasive diagnostic tool is presently the greatest barrier to the identification and treatment of endometriosis. Twin and family studies have documented an increased relative risk in families. To identify genetic factors that contribute to endometriosis we conducted a two-stage genome-wide association study (GWAS) of a European cohort including 2,019 surgically confirmed endometriosis cases and 14,471 controls. Three of the SNPs we identify associated at P<5×10⁻⁸ in our combined analysis belong to two loci: LINC00339-WNT4 on 1p36.12 (rs2235529; P = 8.65×10⁻⁹, OR = 1.29, CI = 1.18–1.40) and RND3-RBM43 on 2q23.3 (rs1519761; P = 4.70×10⁻⁸, OR = 1.20, Cl = 1.13–1.29, and rs6757804; P = 4.05×10⁻⁸, OR = 1.20, Cl = 1.13–1.29). Using an adjusted Bonferoni significance threshold of 4.51×10⁻⁷ we identify two additional loci in our meta-analysis that associate with endometriosis:, RNF144B-ID4 on 6p22.3 (rs6907340; P = 2.19×10⁻⁷, OR = 1.20, Cl = 1.12–1.28), and HNRNPA3P1-LOC100130539 on 10q11.21 (rs10508881; P = 4.08×10⁻⁷, OR = 1.19, Cl = 1.11–1.27). Consistent with previously suggested associations to WNT4 our study implicate a 150 kb region around WNT4 that also include LINC00339 and CDC42. A univariate analysis of documented infertility, age at menarche, and family history did not show allelic association with these SNP markers. Clinical data from patients in our study reveal an average delay in diagnosis of 8.4 years and confirm a strong correlation between endometriosis severity and infertility (n = 1182, P<0.001, OR = 2.18). This GWAS of endometriosis was conducted with high diagnostic certainty in cases, and with stringent handling of population substructure. Our findings broaden the understanding of the genetic factors that play a role in endometriosis.     

A Genome-Wide Association Study of Circulating Galectin-3

Galectin-3 is a lectin involved in fibrosis, inflammation and proliferation. Increased circulating levels of galectin-3 have been associated with various diseases, including cancer, immunological disorders, and cardiovascular disease. To enhance our knowledge on galectin-3 biology we performed the first genome-wide association study (GWAS) using the Illumina HumanCytoSNP-12 array imputed with the HapMap 2 CEU panel on plasma galectin-3 levels in 3,776 subjects and follow-up genotyping in an additional 3,516 subjects. We identified 2 genome wide significant loci associated with plasma galectin-3 levels. One locus harbours the LGALS3 gene (rs2274273; P = 2.35×10⁻¹⁸⁸) and the other locus the ABO gene (rs644234; P = 3.65×10⁻⁴⁷). The variance explained by the LGALS3 locus was 25.6% and by the ABO locus 3.8% and jointly they explained 29.2%. Rs2274273 lies in high linkage disequilibrium with two non-synonymous SNPs (rs4644; r² = 1.0, and rs4652; r² = 0.91) and wet lab follow-up genotyping revealed that both are strongly associated with galectin-3 levels (rs4644; P = 4.97×10⁻⁴⁶⁵ and rs4652 P = 1.50×10⁻⁴²¹) and were also associated with LGALS3 gene-expression. The origins of our associations should be further validated by means of functional experiments.     

A Genome-Wide Association Study Reveals Variants in ARL15 that Influence Adiponectin Levels

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P≤5×10⁻⁸). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P≤0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2×10⁻¹⁹ for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9×10⁻⁸, n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5×10⁻⁶, n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2×10⁻³, n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.     

A Large Scale Gene-Centric Association Study of Lung Function in Newly-Hired Female Cotton Textile Workers with Endotoxin Exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P _trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.     

Multiple Independent Genetic Factors at NOS1AP Modulate the QT Interval in a Multi-Ethnic Population

Extremes of electrocardiographic QT interval are associated with increased risk for sudden cardiac death (SCD); thus, identification and characterization of genetic variants that modulate QT interval may elucidate the underlying etiology of SCD. Previous studies have revealed an association between a common genetic variant in NOS1AP and QT interval in populations of European ancestry, but this finding has not been extended to other ethnic populations. We sought to characterize the effects of NOS1AP genetic variants on QT interval in the multi-ethnic population-based Dallas Heart Study (DHS, n = 3,072). The SNP most strongly associated with QT interval in previous samples of European ancestry, rs16847548, was the most strongly associated in White (P = 0.005) and Black (P = 3.6×10⁻⁵) participants, with the same direction of effect in Hispanics (P = 0.17), and further showed a significant SNP × sex-interaction (P = 0.03). A second SNP, rs16856785, uncorrelated with rs16847548, was also associated with QT interval in Blacks (P = 0.01), with qualitatively similar results in Whites and Hispanics. In a previously genotyped cohort of 14,107 White individuals drawn from the combined Atherosclerotic Risk in Communities (ARIC) and Cardiovascular Health Study (CHS) cohorts, we validated both the second locus at rs16856785 (P = 7.63×10⁻⁸), as well as the sex-interaction with rs16847548 (P = 8.68×10⁻⁶). These data extend the association of genetic variants in NOS1AP with QT interval to a Black population, with similar trends, though not statistically significant at P<0.05, in Hispanics. In addition, we identify a strong sex-interaction and the presence of a second independent site within NOS1AP associated with the QT interval. These results highlight the consistent and complex role of NOS1AP genetic variants in modulating QT interval.     

Genome-Wide Association Study Identifies Major Loci for Carcass Weight on BTA14 in Hanwoo (Korean Cattle)

This genome-wide association study (GWAS) was conducted to identify major loci that are significantly associated with carcass weight, and their effects, in order to provide increased understanding of the genetic architecture of carcass weight in Hanwoo. This genome-wide association study identified one major chromosome region ranging from 23 Mb to 25 Mb on chromosome 14 as being associated with carcass weight in Hanwoo. Significant Bonferroni-corrected genome-wide associations (P<1.52×10⁻⁶) were detected for 6 Single Nucleotide Polymorphic (SNP) loci for carcass weight on chromosome 14. The most significant SNP was BTB-01280026 (P = 4.02×10⁻¹¹), located in the 25 Mb region on Bos taurus autosome 14 (BTA14). The other 5 significant SNPs were Hapmap27934-BTC-065223 (P = 4.04×10⁻¹¹) in 25.2 Mb, BTB-01143580 (P = 6.35×10⁻¹¹) in 24.3 Mb, Hapmap30932-BTC-011225 (P = 5.92×10⁻¹⁰) in 24.8 Mb, Hapmap27112-BTC-063342 (P = 5.18×10⁻⁹) in 25.4 Mb, and Hapmap24414-BTC-073009 (P = 7.38×10⁻⁸) in 25.4 Mb, all on BTA 14. One SNP (BTB-01143580; P = 6.35×10⁻¹¹) lies independently from the other 5 SNPs. The 5 SNPs that lie together showed a large Linkage disequilibrium (LD) block (block size of 553 kb) with LD coefficients ranging from 0.53 to 0.89 within the block. The most significant SNPs accounted for 6.73% to 10.55% of additive genetic variance, which is quite a large proportion of the total additive genetic variance. The most significant SNP (BTB-01280026; P = 4.02×10⁻¹¹) had 16.96 kg of allele substitution effect, and the second most significant SNP (Hapmap27934-BTC-065223; P = 4.04×10⁻¹¹) had 18.06 kg of effect on carcass weight, which correspond to 44% and 47%, respectively, of the phenotypic standard deviation for carcass weight in Hanwoo cattle. Our results demonstrated that carcass weight was affected by a major Quantitative Trait Locus (QTL) with a large effect and by many SNPs with small effects that are normally distributed.     

Novel Association of HK1 with Glycated Hemoglobin in a Non-Diabetic Population: A Genome-Wide Evaluation of 14,618 Participants in the Women's Genome Health Study

Type 2 diabetes is a leading cause of morbidity and mortality. While genetic variants have been found to influence the risk of type 2 diabetes mellitus, relatively few studies have focused on genes associated with glycated hemoglobin, an index of the mean blood glucose concentration of the preceding 8–12 weeks. Epidemiologic studies and randomized clinical trials have documented the relationship between glycated hemoglobin levels and the development of long-term complications in diabetes; moreover, higher glycated hemoglobin levels in the subdiabetic range have been shown to predict type 2 diabetes risk and cardiovascular disease. To examine the common genetic determinants of glycated hemoglobin levels, we performed a genome-wide association study that evaluated 337,343 SNPs in 14,618 apparently healthy Caucasian women. The results show that glycated hemoglobin levels are associated with genetic variation at the GCK (rs730497; P = 2.8×10⁻¹²), SLC30A8 (rs13266634; P = 9.8×10⁻⁸), G6PC2 (rs1402837; P = 6.8×10⁻¹⁰), and HK1 (rs7072268; P = 6.4×10⁻⁹) loci. While associations at the GCK, SLC30A8, and G6PC2 loci are confirmatory, the findings at HK1 are novel. We were able to replicate this novel association in an independent validation sample of 455 additional non-diabetic men and women. HK1 encodes the enzyme hexokinase, the first step in glycolysis and a likely candidate for the control of glucose metabolism. This observed genetic association between glycated hemoglobin levels and HK1 polymorphisms paves the way for further studies of the role of HK1 in hemoglobin glycation, glucose metabolism, and diabetes.     

Common Genetic Variation Near the Phospholamban Gene Is Associated with Cardiac Repolarisation: Meta-Analysis of Three Genome-Wide Association Studies

To identify loci affecting the electrocardiographic QT interval, a measure of cardiac repolarisation associated with risk of ventricular arrhythmias and sudden cardiac death, we conducted a meta-analysis of three genome-wide association studies (GWAS) including 3,558 subjects from the TwinsUK and BRIGHT cohorts in the UK and the DCCT/EDIC cohort from North America. Five loci were significantly associated with QT interval at P<1×10⁻⁶. To validate these findings we performed an in silico comparison with data from two QT consortia: QTSCD (n = 15,842) and QTGEN (n = 13,685). Analysis confirmed the association between common variants near NOS1AP (P = 1.4×10⁻⁸³) and the phospholamban (PLN) gene (P = 1.9×10⁻²⁹). The most associated SNP near NOS1AP (rs12143842) explains 0.82% variance; the SNP near PLN (rs11153730) explains 0.74% variance of QT interval duration. We found no evidence for interaction between these two SNPs (P = 0.99). PLN is a key regulator of cardiac diastolic function and is involved in regulating intracellular calcium cycling, it has only recently been identified as a susceptibility locus for QT interval. These data offer further mechanistic insights into genetic influence on the QT interval which may predispose to life threatening arrhythmias and sudden cardiac death.     

Endometriosis Is Associated with Rare Copy Number Variants

Endometriosis is a complex gynecological condition that affects 6–10% of women in their reproductive years and is defined by the presence of endometrial glands and stroma outside the uterus. Twin, family, and genome-wide association (GWA) studies have confirmed a genetic role, yet only a small part of the genetic risk can be explained by SNP variation. Copy number variants (CNVs) account for a greater portion of human genetic variation than SNPs and include more recent mutations of large effect. CNVs, likely to be prominent in conditions with decreased reproductive fitness, have not previously been examined as a genetic contributor to endometriosis. Here we employ a high-density genotyping microarray in a genome-wide survey of CNVs in a case-control population that includes 2,126 surgically confirmed endometriosis cases and 17,974 population controls of European ancestry. We apply stringent quality filters to reduce the false positive rate common to many CNV-detection algorithms from 77.7% to 7.3% without noticeable reduction in the true positive rate. We detected no differences in the CNV landscape between cases and controls on the global level which showed an average of 1.92 CNVs per individual with an average size of 142.3 kb. On the local level we identify 22 CNV-regions at the nominal significance threshold (P<0.05), which is greater than the 8.15 CNV-regions expected based on permutation analysis (P<0.001). Three CNV''s passed a genome-wide P-value threshold of 9.3×10⁻⁴; a deletion at SGCZ on 8p22 (P = 7.3×10⁻⁴, OR = 8.5, Cl = 2.3–31.7), a deletion in MALRD1 on 10p12.31 (P = 5.6×10⁻⁴, OR = 14.1, Cl = 2.7–90.9), and a deletion at 11q14.1 (P = 5.7×10⁻⁴, OR = 33.8, Cl = 3.3–1651). Two SNPs within the 22 CNVRs show significant genotypic association with endometriosis after adjusting for multiple testing; rs758316 in DPP6 on 7q36.2 (P = 0.0045) and rs4837864 in ASTN2 on 9q33.1 (P = 0.0002). Together, the CNV-loci are detected in 6.9% of affected women compared to 2.1% in the general population.     

Susceptibility Loci Associated with Specific and Shared Subtypes of Lymphoid Malignancies

The genetics of lymphoma susceptibility reflect the marked heterogeneity of diseases that comprise this broad phenotype. However, multiple subtypes of lymphoma are observed in some families, suggesting shared pathways of genetic predisposition to these pathologically distinct entities. Using a two-stage GWAS, we tested 530,583 SNPs in 944 cases of lymphoma, including 282 familial cases, and 4,044 public shared controls, followed by genotyping of 50 SNPs in 1,245 cases and 2,596 controls. A novel region on 11q12.1 showed association with combined lymphoma (LYM) subtypes. SNPs in this region included rs12289961 near LPXN, (P_LYM = 3.89×10⁻⁸, OR = 1.29) and rs948562 (P_LYM = 5.85×10⁻⁷, OR = 1.29). A SNP in a novel non-HLA region on 6p23 (rs707824, P_NHL = 5.72×10⁻⁷) was suggestive of an association conferring susceptibility to lymphoma. Four SNPs, all in a previously reported HLA region, 6p21.32, showed genome-wide significant associations with follicular lymphoma. The most significant association with follicular lymphoma was for rs4530903 (P_FL = 2.69×10⁻¹², OR = 1.93). Three novel SNPs near the HLA locus, rs9268853, rs2647046, and rs2621416, demonstrated additional variation contributing toward genetic susceptibility to FL associated with this region. Genes implicated by GWAS were also found to be cis-eQTLs in lymphoblastoid cell lines; candidate genes in these regions have been implicated in hematopoiesis and immune function. These results, showing novel susceptibility regions and allelic heterogeneity, point to the existence of pathways of susceptibility to both shared as well as specific subtypes of lymphoid malignancy.     

Association of SNPs of CD40 Gene with Multiple Sclerosis in Russians

Multiple sclerosis (MS) is a serious, incurable neurological disease. In 2009, the ANZgene studies detected the suggestive association of located upstream of CD40 gene in chromosome 20q13 (p = 1.3×10⁻⁷). Identification of the causal variant(s) in the CD40 locus leads to a better understanding of the mechanism underlying the development of autoimmune pathologies. We determined the genotypes of rs6074022, rs1883832, rs1535045, and rs11086996 in patients with MS (n = 1684) and in the control group (n = 879). Two SNPs were significantly associated with MS: rs6074022 (additive model C allele OR = 1.27, 95% CI = [1.12–1.45], p = 3×10⁻⁴) and rs1883832 (additive model T allele OR = 1.20, 95% CI = [1.05–1.38], p = 7×10⁻³). In the meta-analysis of our results and the results of four previous studies, we obtain the association p-value of 2.34×10⁻¹², which confirmed the association between MS and rs6074022 at a genome-wide significant level. Next, we demonstrated that the model including rs6074022 only sufficiently described the association. From our analysis, we can speculate that the association between rs1883832 and MS was induced by LD, whereas rs6074022 was a marker in stronger LD with the functional variant or was the functional variant itself. Our results indicated that the functional variants were located in the upstream region of the gene CD40 and were in higher LD with rs6074022 than LD with rs1883832.     

Genome-Wide Association Study for Serum Complement C3 and C4 Levels in Healthy Chinese Subjects

Complement C3 and C4 play key roles in the main physiological activities of complement system, and their deficiencies or over-expression are associated with many clinical infectious or immunity diseases. A two-stage genome-wide association study (GWAS) was performed for serum levels of C3 and C4. The first stage was conducted in 1,999 healthy Chinese men, and the second stage was performed in an additional 1,496 subjects. We identified two SNPs, rs3753394 in CFH gene and rs3745567 in C3 gene, that are significantly associated with serum C3 levels at a genome-wide significance level (P = 7.33×10⁻¹¹ and P = 1.83×10⁻⁹, respectively). For C4, one large genomic region on chromosome 6p21.3 is significantly associated with serum C4 levels. Two SNPs (rs1052693 and rs11575839) were located in the MHC class I area that include HLA-A, HLA-C, and HLA-B genes. Two SNPs (rs2075799 and rs2857009) were located 5′ and 3′ of C4 gene. The other four SNPs, rs2071278, rs3763317, rs9276606, and rs241428, were located in the MHC class II region that includes HLA-DRA, HLA-DRB, and HLA-DQB genes. The combined P-values for those eight SNPs ranged from 3.19×10⁻²² to 5.62×10⁻⁹⁷. HBsAg-positive subjects have significantly lower C3 and C4 protein concentrations compared with HBsAg-negative subjects (P<0.05). Our study is the first GWAS report which shows genetic components influence the levels of complement C3 and C4. Our significant findings provide novel insights of their related autoimmune, infectious diseases, and molecular mechanisms.     

PTPN22 Association in Systemic Lupus Erythematosus (SLE) with Respect to Individual Ancestry and Clinical Sub-Phenotypes

Protein tyrosine phosphatase non-receptor type 22 (PTPN22) is a negative regulator of T-cell activation associated with several autoimmune diseases, including systemic lupus erythematosus (SLE). Missense rs2476601 is associated with SLE in individuals with European ancestry. Since the rs2476601 risk allele frequency differs dramatically across ethnicities, we assessed robustness of PTPN22 association with SLE and its clinical sub-phenotypes across four ethnically diverse populations. Ten SNPs were genotyped in 8220 SLE cases and 7369 controls from in European-Americans (EA), African-Americans (AA), Asians (AS), and Hispanics (HS). We performed imputation-based association followed by conditional analysis to identify independent associations. Significantly associated SNPs were tested for association with SLE clinical sub-phenotypes, including autoantibody profiles. Multiple testing was accounted for by using false discovery rate. We successfully imputed and tested allelic association for 107 SNPs within the PTPN22 region and detected evidence of ethnic-specific associations from EA and HS. In EA, the strongest association was at rs2476601 (P = 4.7×10⁻⁹, OR = 1.40 (95% CI = 1.25–1.56)). Independent association with rs1217414 was also observed in EA, and both SNPs are correlated with increased European ancestry. For HS imputed intronic SNP, rs3765598, predicted to be a cis-eQTL, was associated (P = 0.007, OR = 0.79 and 95% CI = 0.67–0.94). No significant associations were observed in AA or AS. Case-only analysis using lupus-related clinical criteria revealed differences between EA SLE patients positive for moderate to high titers of IgG anti-cardiolipin (aCL IgG >20) versus negative aCL IgG at rs2476601 (P = 0.012, OR = 1.65). Association was reinforced when these cases were compared to controls (P = 2.7×10⁻⁵, OR = 2.11). Our results validate that rs2476601 is the most significantly associated SNP in individuals with European ancestry. Additionally, rs1217414 and rs3765598 may be associated with SLE. Further studies are required to confirm the involvement of rs2476601 with aCL IgG.     

A Genome-Wide Association Study of Total Serum and Mite-Specific IgEs in Asthma Patients

Immunoglobulin E (IgE) is one of the central players in asthma and allergic diseases. Although the serum IgE level, a useful endophenotype, is generally increased in patients with asthma, genetic factors influencing IgE regulation in asthma are still not fully understood. To identify the genetic variations associated with total serum and mite-specific IgEs in asthmatics, a genome-wide association study (GWAS) of 657,366 single nucleotide polymorphisms (SNPs) was performed in 877 Korean asthmatics. This study found that several new genes might be associated with total IgE in asthmatics, such as CRIM1 (rs848512, P = 1.18×10⁻⁶; rs711254, P = 6.73×10⁻⁶), ZNF71 (rs10404342, P = 7.60×10⁻⁶), TLN1 (rs4879926, P = 7.74×10⁻⁶), and SYNPO2 (rs1472066, P = 8.36×10⁻⁶; rs1038770, P = 8.66×10⁻⁶). Regarding the association of specific IgE to house dust mites, it was observed that intergenic SNPs nearby to OPRK1 and LOC730217 might be associated with Dermatophagoides pteronyssinus (D.p.) and Dermatophagoides farinae (D.f.) in asthmatics, respectively. In further pathway analysis, the phosphatidylinositol signaling system and adherens junction pathways were estimated to play a role in the regulation of total IgE levels in asthma. Although functional evaluations and replications of these results in other populations are needed, this GWAS of serum IgE in asthmatics could facilitate improved understanding of the role of the newly identified genetic variants in asthma and its related phenotypes.     

Genome-Wide Association Study in BRCA1 Mutation Carriers Identifies Novel Loci Associated with Breast and Ovarian Cancer Risk

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7×10⁻⁸, HR = 1.14, 95% CI: 1.09–1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4×10⁻⁸, HR = 1.27, 95% CI: 1.17–1.38) and 4q32.3 (rs4691139, P = 3.4×10⁻⁸, HR = 1.20, 95% CI: 1.17–1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific association. The 17q21.31 locus was also associated with ovarian cancer risk in 8,211 BRCA2 carriers (P = 2×10⁻⁴). These loci may lead to an improved understanding of the etiology of breast and ovarian tumors in BRCA1 carriers. Based on the joint distribution of the known BRCA1 breast cancer risk-modifying loci, we estimated that the breast cancer lifetime risks for the 5% of BRCA1 carriers at lowest risk are 28%–50% compared to 81%–100% for the 5% at highest risk. Similarly, based on the known ovarian cancer risk-modifying loci, the 5% of BRCA1 carriers at lowest risk have an estimated lifetime risk of developing ovarian cancer of 28% or lower, whereas the 5% at highest risk will have a risk of 63% or higher. Such differences in risk may have important implications for risk prediction and clinical management for BRCA1 carriers.     

Polymorphisms in the TNFA and IL6 Genes Represent Risk Factors for Autoimmune Thyroid Disease

Background

Autoimmune thyroid disease (AITD) comprises diseases including Hashimoto''s thyroiditis and Graves'' disease, both characterized by reactivity to autoantigens causing, respectively, inflammatory destruction and autoimmune stimulation of the thyroid-stimulating hormone receptor. AITD is the most common thyroid disease and the leading form of autoimmune disease in women. Cytokines are key regulators of the immune and inflammatory responses; therefore, genetic variants at cytokine-encoding genes are potential risk factors for AITD.

Methods

Polymorphisms in the IL6-174 G/C (rs1800795), TNFA-308 G/A (rs1800629), IL1B-511 C/T (rs16944), and IFNGR1-56 T/C (rs2234711) genes were assessed in a case-control study comprising 420 Hashimoto''s thyroiditis patients, 111 Graves'' disease patients and 735 unrelated controls from Portugal. Genetic variants were discriminated by real-time PCR using TaqMan SNP genotyping assays.

Results

A significant association was found between the allele A in TNFA-308 G/A and Hashimoto''s thyroiditis, both in the dominant (OR = 1.82, CI = 1.37–2.43, p-value = 4.4×10⁻⁵) and log-additive (OR = 1.64, CI = 1.28–2.10, p-value = 8.2×10⁻⁵) models. The allele C in IL6-174 G/C is also associated with Hashimoto''s thyroiditis, however, only retained significance after multiple testing correction in the log-additive model (OR = 1.28, CI = 1.06–1.54, p-value = 8.9×10⁻³). The group with Graves'' disease also registered a higher frequency of the allele A in TNFA-308 G/A compared with controls both in the dominant (OR = 1.85, CI = 1.19–2.87, p-value = 7.0×10⁻³) and log-additive (OR = 1.69, CI = 1.17–2.44, p-value = 6.6×10⁻³) models. The risk for Hashimoto''s thyroiditis and Graves'' disease increases with the number of risk alleles (OR for two risk alleles is, respectively, 2.27 and 2.59).

Conclusions

This study reports significant associations of genetic variants in TNFA and IL6 with the risk for AITD, highlighting the relevance of polymorphisms in inflammation-related genes in the etiopathogenesis of AITD.