Genome-Wide Association Study of Breast Cancer in the Japanese Population

Breast cancer is the most common malignancy among women in worldwide including Japan. Several studies have identified common genetic variants to be associated with the risk of breast cancer. Due to the complex linkage disequilibrium structure and various environmental exposures in different populations, it is essential to identify variants associated with breast cancer in each population, which subsequently facilitate the better understanding of mammary carcinogenesis. In this study, we conducted a genome-wide association study (GWAS) as well as whole-genome imputation with 2,642 cases and 2,099 unaffected female controls. We further examined 13 suggestive loci (P<1.0×10⁻⁵) using an independent sample set of 2,885 cases and 3,395 controls and successfully validated two previously-reported loci, rs2981578 (combined P-value of 1.31×10⁻¹², OR = 1.23; 95% CI = 1.16–.30) on chromosome 10q26 (FGFR2), rs3803662 (combined P-value of 2.79×10⁻¹¹, OR = 1.21; 95% CI = 1.15–.28) and rs12922061 (combined P-value of 3.97×10⁻¹⁰, OR = 1.23; 95% CI = 1.15–.31) on chromosome 16q12 (TOX3-LOC643714). Weighted genetic risk score on the basis of three significantly associated variants and two previously reported breast cancer associated loci in East Asian population revealed that individuals who carry the most risk alleles in category 5 have 2.2 times higher risk of developing breast cancer in the Japanese population than those who carry the least risk alleles in reference category 1. Although we could not identify additional loci associated with breast cancer, our study utilized one of the largest sample sizes reported to date, and provided genetic status that represent the Japanese population. Further local and international collaborative study is essential to identify additional genetic variants that could lead to a better, accurate prediction for breast cancer.     

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.     

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.     

Association of Variants at UMOD with Chronic Kidney Disease and Kidney Stones—Role of Age and Comorbid Diseases

Chronic kidney disease (CKD) is a worldwide public health problem that is associated with substantial morbidity and mortality. To search for sequence variants that associate with CKD, we conducted a genome-wide association study (GWAS) that included a total of 3,203 Icelandic cases and 38,782 controls. We observed an association between CKD and a variant with 80% population frequency, rs4293393-T, positioned next to the UMOD gene (GeneID: 7369) on chromosome 16p12 (OR = 1.25, P = 4.1×10⁻¹⁰). This gene encodes uromodulin (Tamm-Horsfall protein), the most abundant protein in mammalian urine. The variant also associates significantly with serum creatinine concentration (SCr) in Icelandic subjects (N = 24,635, P = 1.3×10⁻²³) but not in a smaller set of healthy Dutch controls (N = 1,819, P = 0.39). Our findings validate the association between the UMOD variant and both CKD and SCr recently discovered in a large GWAS. In the Icelandic dataset, we demonstrate that the effect on SCr increases substantially with both age (P = 3.0×10⁻¹⁷) and number of comorbid diseases (P = 0.008). The association with CKD is also stronger in the older age groups. These results suggest that the UMOD variant may influence the adaptation of the kidney to age-related risk factors of kidney disease such as hypertension and diabetes. The variant also associates with serum urea (P = 1.0×10⁻⁶), uric acid (P = 0.0064), and suggestively with gout. In contrast to CKD, the UMOD variant confers protection against kidney stones when studied in 3,617 Icelandic and Dutch kidney stone cases and 43,201 controls (OR = 0.88, P = 5.7×10⁻⁵).     

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.     

Common Variants of FTO Are Associated with Childhood Obesity in a Cross-Sectional Study of 3,126 Urban Indian Children

Background

FTO variants are robustly associated with obesity and related traits in many population and shown to have variable impact during life course. Although studies have shown association of FTO variants with adiposity in adult Indian, its association in Indian children is yet to be confirmed.

Methods

Here we examined association of FTO variants (rs9939609 and rs8050136) with obesity and related anthropometric and biochemical traits in 3,126 Indian children (aged 11–17 years) including 2,230 normal-weight and 896 over-weight/obese children. We also compared effects observed in the present study with that observed in previous studies on South Asian adults and children of other ethnic groups.

Results

The variant rs9939609 showed significant association with risk of obesity [OR = 1.21, P = 2.5×10⁻³] and its measures BMI, weight, waist circumference and hip circumference [β range = 0.11 to 0.14 Z-score units; P range = 1.3×10⁻⁴ to 1.6×10⁻⁷] in children. The observed effect sizes in Indian children were similar to those reported for European children. Variant rs9939609 explained 0.88% of BMI variance in Indian children. The effect sizes of rs9939609 on BMI and WC were ∼2 fold higher in children than adults. Interestingly rs9939609 was also associated with serum levels of thyroid stimulating hormone (TSH) [β = 0.10 Z-score, P = 5.8×10⁻³]. The other variant rs8050136 was in strong linkage disequilibrium with rs9939609 (r² = 0.97) and provided similar association results.

Conclusion

The study provides first report of association of FTO variants with obesity and related anthropometric traits in Indian children with higher impact in children compared to adults. We also demonstrated association of FTO variant with serum levels of TSH, indicating putative influence of FTO in hypothalamic-pituitary-thyroid axis.     

Ancestry-Shift Refinement Mapping of the C6orf97-ESR1 Breast Cancer Susceptibility Locus

We used an approach that we term ancestry-shift refinement mapping to investigate an association, originally discovered in a GWAS of a Chinese population, between rs2046210[T] and breast cancer susceptibility. The locus is on 6q25.1 in proximity to the C6orf97 and estrogen receptor α (ESR1) genes. We identified a panel of SNPs that are correlated with rs2046210 in Chinese, but not necessarily so in other ancestral populations, and genotyped them in breast cancer case∶control samples of Asian, European, and African origin, a total of 10,176 cases and 13,286 controls. We found that rs2046210[T] does not confer substantial risk of breast cancer in Europeans and Africans (OR = 1.04, P = 0.099, and OR = 0.98, P = 0.77, respectively). Rather, in those ancestries, an association signal arises from a group of less common SNPs typified by rs9397435. The rs9397435[G] allele was found to confer risk of breast cancer in European (OR = 1.15, P = 1.2×10⁻³), African (OR = 1.35, P = 0.014), and Asian (OR = 1.23, P = 2.9×10⁻⁴) population samples. Combined over all ancestries, the OR was 1.19 (P = 3.9×10⁻⁷), was without significant heterogeneity between ancestries (P_het = 0.36) and the SNP fully accounted for the association signal in each ancestry. Haplotypes bearing rs9397435[G] are well tagged by rs2046210[T] only in Asians. The rs9397435[G] allele showed associations with both estrogen receptor positive and estrogen receptor negative breast cancer. Using early-draft data from the 1,000 Genomes project, we found that the risk allele of a novel SNP (rs77275268), which is closely correlated with rs9397435, disrupts a partially methylated CpG sequence within a known CTCF binding site. These studies demonstrate that shifting the analysis among ancestral populations can provide valuable resolution in association mapping.     

Association of the CTLA4 Gene with Graves' Disease in the Chinese Han Population

To determine whether genetic heterogeneity exists in patients with Graves'' disease (GD), the cytotoxic T-lymphocyte associated 4 (CTLA-4) gene, which is implicated a susceptibility gene for GD by considerable genetic and immunological evidence, was used for association analysis in a Chinese Han cohort recruited from various geographic regions. Our association study for the SNPs in the CTLA4 gene in 2640 GD patients and 2204 control subjects confirmed that CTLA4 is the susceptibility gene for GD in the Chinese Han population. Moreover, the logistic regression analysis in the combined Chinese Han cohort revealed that SNP rs231779 (allele frequencies p = 2.81×10⁻⁹, OR = 1.35, and genotype distributions p = 2.75×10⁻⁹, OR = 1.42) is likely the susceptibility variant for GD. Interestingly, the logistic regression analysis revealed that SNP rs35219727 may be the susceptibility variant to GD in the Shandong population; however, SNP, rs231779 in the CTLA4 gene probably independently confers GD susceptibility in the Xuzhou and southern China populations. These data suggest that the susceptibility variants of the CTLA4 gene varied between the different geographic populations with GD.     

IFNL4 ss469415590 Variant Shows Similar Performance to rs12979860 as Predictor of Response to Treatment against Hepatitis C Virus Genotype 1 or 4 in Caucasians

Objectives

The rs12979860 variant, linked to IL28B gene, predicts sustained viral response (SVR) to pegylated-interferon/ribavirin (pegIFN/RBV) therapy in Hepatitis C Virus genotype 1 or 4 (HCV-1/4)-infected patients. Recently, a functional variant, ss469415590, in linkage disequilibrium (LD) with rs12979860, has been discovered. Our objective was to assess the value of ss469415590 to predict SVR to pegIFN/RBV in Caucasian HCV-1/4-infected individuals and to compare its performance with that of rs12979860.

Methods

272 Caucasian HCV-1/4-infected patients who completed a course of pegIFN/RBV were genotyped for both rs12979860 and ss469415590 markers. Logistic regression models including factors with univariate association with SVR and each genetic marker were elaborated. The area under the receiver operating-characteristic curve (AUROC) was calculated for each model and both were compared.

Results

Both markers were in LD (r² = 0.82). For rs12979860, 66 (64.0%) CC versus 56 (33.1%) T allele carriers achieved SVR (Adjusted OR = 4.156, 95%CI = 2.388–7.232, p = 4.647×10⁻⁷). For ss469415590, 66 (66.0%) TT/TT versus 56 (32.5%) –G allele carriers (Adjusted OR = 4.783, 95%CI = 2.714–8.428, p = 6.153×10⁻⁸) achieved SVR. The AUROC of the model including rs12979860 was 0.742 (95%CI = 0.672–0.813) and of that based on ss469415590 was 0.756 (95%CI = 0.687–0.826) (p = 0.780).

Conclusions

The ss469415590 variant shows an equivalent performance to predict SVR to pegIFN/RBV than the rs2979860 in Caucasian HCV-1/4-infected patients.     

A Genome-Wide Association Study of Psoriasis and Psoriatic Arthritis Identifies New Disease Loci

A genome-wide association study was performed to identify genetic factors involved in susceptibility to psoriasis (PS) and psoriatic arthritis (PSA), inflammatory diseases of the skin and joints in humans. 223 PS cases (including 91 with PSA) were genotyped with 311,398 single nucleotide polymorphisms (SNPs), and results were compared with those from 519 Northern European controls. Replications were performed with an independent cohort of 577 PS cases and 737 controls from the U.S., and 576 PSA patients and 480 controls from the U.K.. Strongest associations were with the class I region of the major histocompatibility complex (MHC). The most highly associated SNP was rs10484554, which lies 34.7 kb upstream from HLA-C (P = 7.8×10⁻¹¹, GWA scan; P = 1.8×10⁻³⁰, replication; P = 1.8×10⁻³⁹, combined; U.K. PSA: P = 6.9×10⁻¹¹). However, rs2395029 encoding the G2V polymorphism within the class I gene HCP5 (combined P = 2.13×10⁻²⁶ in U.S. cases) yielded the highest ORs with both PS and PSA (4.1 and 3.2 respectively). This variant is associated with low viral set point following HIV infection and its effect is independent of rs10484554. We replicated the previously reported association with interleukin 23 receptor and interleukin 12B (IL12B) polymorphisms in PS and PSA cohorts (IL23R: rs11209026, U.S. PS, P = 1.4×10⁻⁴; U.K. PSA: P = 8.0×10⁻⁴; IL12B:rs6887695, U.S. PS, P = 5×10⁻⁵ and U.K. PSA, P = 1.3×10⁻³) and detected an independent association in the IL23R region with a SNP 4 kb upstream from IL12RB2 (P = 0.001). Novel associations replicated in the U.S. PS cohort included the region harboring lipoma HMGIC fusion partner (LHFP) and conserved oligomeric golgi complex component 6 (COG6) genes on chromosome 13q13 (combined P = 2×10⁻⁶ for rs7993214; OR = 0.71), the late cornified envelope gene cluster (LCE) from the Epidermal Differentiation Complex (PSORS4) (combined P = 6.2×10⁻⁵ for rs6701216; OR 1.45) and a region of LD at 15q21 (combined P = 2.9×10⁻⁵ for rs3803369; OR = 1.43). This region is of interest because it harbors ubiquitin-specific protease-8 whose processed pseudogene lies upstream from HLA-C. This region of 15q21 also harbors the gene for SPPL2A (signal peptide peptidase like 2a) which activates tumor necrosis factor alpha by cleavage, triggering the expression of IL12 in human dendritic cells. We also identified a novel PSA (and potentially PS) locus on chromosome 4q27. This region harbors the interleukin 2 (IL2) and interleukin 21 (IL21) genes and was recently shown to be associated with four autoimmune diseases (Celiac disease, Type 1 diabetes, Grave''s disease and Rheumatoid Arthritis).     

STK39 Polymorphism Is Associated with Essential Hypertension: A Systematic Review and Meta-Analysis

Background

A recent genome-wide association study identified STK39as a candidate gene for blood pressure (BP) in Europeans. Subsequently, several studies have attempted to replicate the association across different ethnic populations. However, the results have been inconsistent.

Objective and Methods

We performed a meta-analysis to elucidate the association between the STK39 rs3754777 polymorphism (or proxy) and hypertension. Published literature from PubMed and Embase databases were retrieved and pooled odds ratio (OR) with 95% confidence interval (CI) was calculated using fixed- or random-effects model.

Results

Using appropriate inclusion/exclusion criteria, we identified 10 studies that included 21, 863 hypertensive cases and 24, 480 controls from different ethnicities. The meta-analysis showed a significant association of STK39 rs3754777 variant with hypertension (OR = 1.10, 95%CI = 1.06–1.15, p = 7.95×10⁻⁶). Further subgroup analysis by ethnicity suggested that the association was significant in Europeans (OR = 1.08, 95% CI = 1.03–1.14, p = 0.002) and in East Asians (OR = 1.16, 95%CI = 1.07–1.25, p = 4.34×10⁻⁴), but not in Africans (OR = 1.01, 95%CI 0.80–1.27, p = 0.932). We further confirmed the positive association by sensitivity analysis. No publication bias was detected (Begg’s test, p = 0.721; Egger’s test, p = 0.744).

Conclusions

The present meta-analysis confirms the significant association of STK39 polymorphism with susceptibility to hypertension in Europeans and East Asians. Future studies should include gene–gene and gene–environment interactions to investigate the identified association.     

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.     

TNFAIP3 Gene Polymorphisms in a Chinese Han Population with Vogt–Koyanagi–Harada Syndrome

Background

This study was performed to evaluate the potential association of TNFAIP3 polymorphisms with Vogt–Koyanagi–Harada (VKH) disease in a Chinese Han population.

Methodology/Principal Findings

Five single-nucleotide polymorphisms (SNPs), rs10499194, rs610604, rs7753873, rs5029928 and rs9494885 of TNFAIP3 were genotyped in 834 VKH disease patients and 1415 healthy controls using a PCR-restriction fragment length polymorphism assay. An increased frequency of the C allele and CT genotype for rs9494885 were found in VKH patients in the Guangzhou and Chongqing cohorts (pc = 0.015, OR = 1.6, pc = 0.036, OR = 1.7; pc = 2.36×10−4, OR = 1.5, pc = 0.012, OR = 1.5, respectively). Meanwhile, a decreased frequency of the TT genotype for rs9494885 was observed in VKH patients in the Guangzhou and Chongqing cohorts (pc = 0.026, OR = 0.6, pc = 0.0074, OR = 0.7, respectively). The combined analysis showed that a significantly increased prevalence of the rs9494885 TC genotype and C allele were found in VKH disease patients compared with controls (pc = 2.26×10−5, OR = 1.7; pc = 1.09× 10−5, OR = 1.6, respectively). The frequency of the TT genotype of rs9494885 was markedly lower in VKH disease patients as compared with that in controls (pc = 1.12×10−5, OR = 0.6; p_c = 1.09×10⁻⁵, OR = 0.6, respectively). No association was found between rs10499194, rs610604, rs7753873 and rs5029928 polymorphisms and VKH disease. To our knowledge this is the first report describing the association of a TNFAIP3 gene polymorphism with VKH disease in a Chinese Han population.

Conclusions/Significance

The results suggest that the rs9494885 TC genotype and C allele may be predisposing factors to VKH disease, whereas the rs9494885 TT genotype and T allele may provide protection against this disease.     

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.     

Interactions between Urinary 4-tert-Octylphenol Levels and Metabolism Enzyme Gene Variants on Idiopathic Male Infertility

Octylphenol (OP) and Trichlorophenol (TCP) act as endocrine disruptors and have effects on male reproductive function. We studied the interactions between 4-tert-Octylphenol (4-t-OP), 4-n- Octylphenol (4-n-OP), 2,3,4-Trichlorophenol (2,3,4-TCP), 2,4,5-Trichlorophenol (2,4,5-TCP) urinary exposure levels and polymorphisms in selected xenobiotic metabolism enzyme genes among 589 idiopathic male infertile patients and 396 controls in a Han-Chinese population. Ultra high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) was used to measure alkylphenols and chlorophenols in urine. Polymorphisms were genotyped using the SNPstream platform and the Taqman method. Among four phenols that were detected, we found that only exposure to 4-t-OP increased the risk of male infertility (P _trend = 1.70×10⁻⁷). The strongest interaction was between 4-t-OP and rs4918758 in CYP2C9 (P _inter = 6.05×10⁻⁷). It presented a significant monotonic increase in risk estimates for male infertility with increasing 4-t-OP exposure levels among men with TC/CC genotype (low level compared with non-exposed, odds ratio (OR) = 2.26, 95% confidence intervals (CI) = 1.06, 4.83; high level compared with non-exposed, OR = 9.22, 95% CI = 2.78, 30.59), but no associations observed among men with TT genotype. We also found interactions between 4-t-OP and rs4986894 in CYP2C19, and between rs1048943 in CYP1A1, on male infertile risk (P _inter = 8.09×10⁻⁷, P _inter = 3.73×10⁻⁴, respectively).We observed notable interactions between 4-t-OP exposure and metabolism enzyme gene polymorphisms on idiopathic infertility in Han-Chinese men.     

Association of Aminoacyl-tRNA Synthetases Gene Polymorphisms with the Risk of Congenital Heart Disease in the Chinese Han Population

Aminoacyl-tRNA synthetases (ARSs) are in charge of cellular protein synthesis and have additional domains that function in a versatile manner beyond translation. Eight core ARSs (EPRS, MRS, QRS, RRS, IRS, LRS, KRS, DRS) combined with three nonenzymatic components form a complex known as multisynthetase complex (MSC).We hypothesize that the single-nucleotide polymorphisms (SNPs) of the eight core ARS coding genes might influence the susceptibility of sporadic congenital heart disease (CHD). Thus, we conducted a case-control study of 984 CHD cases and 2953 non-CHD controls in the Chinese Han population to evaluate the associations of 16 potentially functional SNPs within the eight ARS coding genes with the risk of CHD. We observed significant associations with the risk of CHD for rs1061248 [G/A; odds ratio (OR) = 0.90, 95% confidence interval (CI) = 0.81–0.99; P = 3.81×10⁻²], rs2230301 [A/C; OR = 0.73, 95%CI = 0.60–0.90, P = 3.81×10⁻²], rs1061160 [G/A; OR = 1.18, 95%CI = 1.06–1.31; P = 3.53×10⁻³] and rs5030754 [G/A; OR = 1.39, 95%CI = 1.11–1.75; P = 4.47×10⁻³] of EPRS gene. After multiple comparisons, rs1061248 conferred no predisposition to CHD. Additionally, a combined analysis showed a significant dosage-response effect of CHD risk among individuals carrying the different number of risk alleles (P _trend = 5.00×10⁻⁴). Compared with individuals with “0–2” risk allele, those carrying “3”, “4” or “5 or more” risk alleles had a 0.97-, 1.25- or 1.38-fold increased risk of CHD, respectively. These findings indicate that genetic variants of the EPRS gene may influence the individual susceptibility to CHD in the Chinese Han population.     

Leprosy and the Adaptation of Human Toll-Like Receptor 1

Leprosy is an infectious disease caused by the obligate intracellular pathogen Mycobacterium leprae and remains endemic in many parts of the world. Despite several major studies on susceptibility to leprosy, few genomic loci have been replicated independently. We have conducted an association analysis of more than 1,500 individuals from different case-control and family studies, and observed consistent associations between genetic variants in both TLR1 and the HLA-DRB1/DQA1 regions with susceptibility to leprosy (TLR1 I602S, case-control P = 5.7×10⁻⁸, OR = 0.31, 95% CI = 0.20–0.48, and HLA-DQA1 rs1071630, case-control P = 4.9×10⁻¹⁴, OR = 0.43, 95% CI = 0.35–0.54). The effect sizes of these associations suggest that TLR1 and HLA-DRB1/DQA1 are major susceptibility genes in susceptibility to leprosy. Further population differentiation analysis shows that the TLR1 locus is extremely differentiated. The protective dysfunctional 602S allele is rare in Africa but expands to become the dominant allele among individuals of European descent. This supports the hypothesis that this locus may be under selection from mycobacteria or other pathogens that are recognized by TLR1 and its co-receptors. These observations provide insight into the long standing host-pathogen relationship between human and mycobacteria and highlight the key role of the TLR pathway in infectious diseases.     

Association of CD247 Polymorphisms with Rheumatoid Arthritis: A Replication Study and a Meta-Analysis

Given the role of CD247 in the response of the T cells, its entailment in autoimmune diseases and in order to better clarify the role of this gene in RA susceptibility, we aimed to analyze CD247 gene variants previously associated with other autoimmune diseases (rs1052237, rs2056626 and rs864537) in a large independent European Caucasian population. However, no evidence of association was found for the analyzed CD247 single-nucleotide polymorphisms (SNPs) with RA and with the presence/absence of anti-cyclic citrullinated polypeptide. We performed a meta-analysis including previously published GWAS data from the rs864537 variant, revealing an overall genome-wide significant association between this CD247 SNP and RA with anti-CCP (OR = 0.90, CI 95% = 0.87–0.93, P_overall = 2.1×10⁻¹⁰). Our results show for first time a GWAS-level association between this CD247 polymorphism and RA risk.     

Contribution of Common PCSK1 Genetic Variants to Obesity in 8,359 Subjects from Multi-Ethnic American Population

Common PCSK1 variants (notably rs6232 and rs6235) have been shown to be associated with obesity in European, Asian and Mexican populations. To determine whether common PCSK1 variants contribute to obesity in American population, we conducted association analyses in 8,359 subjects using two multi-ethnic American studies: the Coronary Artery Risk Development in Young Adults (CARDIA) study and the Multi-Ethnic Study of Atherosclerosis (MESA). By evaluating the contribution of rs6232 and rs6235 in each ethnic group, we found that in European-American subjects from CARDIA, only rs6232 was associated with BMI (P = 0.006) and obesity (P = 0.018) but also increased the obesity incidence during the 20 years of follow-up (HR = 1.53 [1.07–2.19], P = 0.019). Alternatively, in African-American subjects from CARDIA, rs6235 was associated with BMI (P = 0.028) and obesity (P = 0.018). Further, by combining the two case-control ethnic groups from the CARDIA study in a meta-analysis, association between rs6235 and obesity risk remained significant (OR = 1.23 [1.05–1.45], P = 9.5×10⁻³). However, neither rs6232 nor rs6235 was associated with BMI or obesity in the MESA study. Interestingly, rs6232 was associated with BMI (P = 4.2×10⁻³) and obesity (P = 3.4×10⁻³) in the younger European-American group combining samples from the both studies [less than median age (53 years)], but not among the older age group (P = 0.756 and P = 0.935 for BMI and obesity, respectively). By combining all the case-control ethnic groups from CARDIA and MESA in a meta-analysis, we found no significant association for the both variants and obesity risk. Finally, by exploring the full PCSK1 locus, we observed that no variant remained significant after correction for multiple testing. These results indicate that common PCSK1 variants (notably rs6232 and rs6235) contribute modestly to obesity in multi-ethnic American population. Further, these results suggest that the association of rs6232 with obesity may be age-dependent in European-Americans. However, multiple replication studies in multi-ethnic American population are needed to confirm our findings.     

The PSEN1, p.E318G Variant Increases the Risk of Alzheimer's Disease in APOE-ε4 Carriers

The primary constituents of plaques (Aβ42/Aβ40) and neurofibrillary tangles (tau and phosphorylated forms of tau [ptau]) are the current leading diagnostic and prognostic cerebrospinal fluid (CSF) biomarkers for AD. In this study, we performed deep sequencing of APP, PSEN1, PSEN2, GRN, APOE and MAPT genes in individuals with extreme CSF Aβ42, tau, or ptau levels. One known pathogenic mutation (PSEN1 p.A426P), four high-risk variants for AD (APOE p.L46P, MAPT p.A152T, PSEN2 p.R62H and p.R71W) and nine novel variants were identified. Surprisingly, a coding variant in PSEN1, p.E318G (rs17125721-G) exhibited a significant association with high CSF tau (p = 9.2×10⁻⁴) and ptau (p = 1.8×10⁻³) levels. The association of the p.E318G variant with Aβ deposition was observed in APOE-ε4 allele carriers. Furthermore, we found that in a large case-control series (n = 5,161) individuals who are APOE-ε4 carriers and carry the p.E318G variant are at a risk of developing AD (OR = 10.7, 95% CI = 4.7–24.6) that is similar to APOE-ε4 homozygous (OR = 9.9, 95% CI = 7.2.9–13.6), and double the risk for APOE-ε4 carriers that do not carry p.E318G (OR = 3.9, 95% CI = 3.4–4.4). The p.E318G variant is present in 5.3% (n = 30) of the families from a large clinical series of LOAD families (n = 565) and exhibited a higher frequency in familial LOAD (MAF = 2.5%) than in sporadic LOAD (MAF = 1.6%) (p = 0.02). Additionally, we found that in the presence of at least one APOE-ε4 allele, p.E318G is associated with more Aβ plaques and faster cognitive decline. We demonstrate that the effect of PSEN1, p.E318G on AD susceptibility is largely dependent on an interaction with APOE-ε4 and mediated by an increased burden of Aβ deposition.