Dietary cholesterol increases paraoxonase 1 enzyme activity

HDL-associated paraoxonase 1 (PON1) activity has been consistently associated with cardiovascular and other diseases. Vitamins C and E intake have previously been positively associated with PON1 in a subset of the Carotid Lesion Epidemiology and Risk (CLEAR) cohort. The goal of this study was to replicate these findings and determine whether other nutrient intake affected PON1 activity. To predict nutrient and mineral intake values, 1,402 subjects completed a standardized food frequency survey of their dietary habits over the past year. Stepwise regression was used to evaluate dietary and covariate effects on PON1 arylesterase activity. Five dietary components, cholesterol (P < 2.0 × 10⁻¹⁶), alcohol (P = 8.51 × 10⁻⁸), vitamin C (P = 7.97 × 10⁻⁵), iron (P = 0.0026), and folic acid (0.037) were independently predictive of PON1 activity. Dietary cholesterol was positively associated and predicted 5.5% of PON1 activity, second in variance explained. This study presents a novel finding of dietary cholesterol, iron, and folic acid predicting PON1 activity in humans and confirms prior reported associations, including that with vitamin C. Identifying and understanding environmental factors that affect PON1 activity is necessary to understand its role and that of HDL in human disease.     

Novel common and rare genetic determinants of paraoxonase activity: FTO,SERPINA12, and ITGAL

HDL-associated paraoxonase 1 (PON1) activity is associated with cardiovascular and other human diseases. As the role of genetic variants outside of the PON gene cluster on PON1 activity is unknown, we sought to identify common and rare variants in such loci. We typed 33,057 variants on the CVD chip in 1,362 subjects to test for their effects on adjusted-PON1 activity. Three novel genes (FTO, ITGAL, and SERPINA12) and the PON gene cluster had SNPs associated with PON1 arylesterase (AREase) activity. These loci were carried forward for rare-variant analysis using Exome chip genotypes in an overlapping subset of 1,051 subjects using sequence kernel association testing. PON1 (P = 2.24 × 10⁻⁴), PON3 (P = 0.022), FTO (P = 0.019), and SERPINA12 (P = 0.039) had both common and rare variants associated with PON1 AREase. ITGAL variants were associated with PON1 activity when using weighted sequence kernel association testing (SKAT) analysis (P = 2.63 × 10⁻³). When adjusting for the initial common variants, SERPINA12 became marginally significant (P = 0.09), whereas all other findings remained significant (P < 0.05), suggesting independent rare-variant effects. We present novel findings that common and rare variants in FTO, SERPINA12, and ITGAL predict PON1 activity. These results further link PON1 to diabetes and inflammation and may inform the role of HDL in human disease.     

A common variant highly associated with plasma VEGFA levels also contributes to the variation of both LDL-C and HDL-C

Vascular endothelial growth factor A (VEGFA) is among the most-significant stimulators of angiogenesis. Its effect on cardiovascular diseases and on the variation of related risk factors such as lipid parameters is considered important, although as yet unclear. Recently, we identified four common variants (rs6921438, rs4416670, rs6993770, and rs10738760) that explain up to 50% of the heritability of plasma VEGFA levels. In the present study, we aimed at assessing the contribution of these variants to the variation of blood lipid levels (including apoE, triglycerides, total cholesterol, low- and high-density lipoprotein cholesterol levels (LDL-C and HDL-C)] in healthy subjects. The effect of these single-nucleotide polymorphisms (SNPs) on lipid levels was assessed using linear regression in discovery and replication samples (n = 1,006 and n = 1,145; respectively), followed by a meta-analysis. Their gene×gene and gene×environment interactions were also assessed. SNP rs6921438 was associated with HDL-C (β = −0.08 mmol/l, P_overall = 1.2 × 10⁻⁷) and LDL-C (β = 0.13 mmol/l, P_overall = 1.5 × 10⁻⁴). We also identified a significant association between the interaction rs4416670×hypertension and apoE variation (P_overall = 1.7 × 10⁻⁵). Therefore, our present study shows a common genetic regulation between VEGFA and cholesterol homeostasis molecules. The SNP rs6921438 is in linkage disequilibrium with variants located in an enhancer- and promoter-associated histone mark region and could have a regulatory effect in the expression of surrounding genes, including VEGFA.     

Genetic variation at the FADS1-FADS2 gene locus influences delta-5 desaturase activity and LC-PUFA proportions after fish oil supplement

Delta-5 and delta-6 desaturases (D5D and D6D) are key enzymes in endogenous synthesis of long-chain PUFAs. In this sample of healthy subjects (n = 310), genotypes of single nucleotide polymorphisms (SNPs) rs174537, rs174561, and rs3834458 in the FADS1-FADS2 gene cluster were strongly associated with proportions of LC-PUFAs and desaturase activities estimated in plasma and erythrocytes. Minor allele carriage associated with decreased activities of D5D (FADS1) (5.84 × 10⁻¹⁹ ≤ P ≤ 4.5 × 10⁻¹⁸) and D6D (FADS2) (6.05 × 10⁻⁸ ≤ P ≤ 4.20 × 10⁻⁷) was accompanied by increased substrate and decreased product proportions (0.05 ≤ P ≤ 2.49 × 10⁻¹⁶). The significance of haplotype association with D5D activity (P = 2.19 × 10⁻¹⁷) was comparable to that of single SNPs, but haplotype association with D6D activity (P = 3.39 × 10⁻²⁸) was much stronger. In a randomized controlled dietary intervention, increasing eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3) intake significantly increased D5D (P = 4.0 × 10⁻⁹) and decreased D6D activity (P = 9.16 × 10⁻⁶) after doses of 0.45, 0.9, and 1.8 g/day for six months. Interaction of rs174537 genotype with treatment was a determinant of D5D activity estimated in plasma (P = 0.05). In conclusion, different sites at the FADS1-FADS2 locus appear to influence D5D and D6D activity, and rs174537 genotype interacts with dietary EPA+DHA to modulate D5D.     

Inherited genetic variants associated with glucocorticoid sensitivity in leukaemia cells

Identification of genetic variants associated with glucocorticoids (GC) sensitivity of leukaemia cells may provide insight into potential drug targets and tailored therapy. In the present study, within 72 leukaemic cell lines derived from Japanese patients with B‐cell precursor acute lymphoblastic leukaemia (ALL), we conducted genome‐wide genotyping of single nucleotide polymorphisms (SNP) and attempted to identify genetic variants associated with GC sensitivity and NR3C1 (GC receptor) gene expression. IC50 measures for prednisolone (Pred) and dexamethasone (Dex) were available using an alamarBlue cell viability assay. IC50 values of Pred showed the strongest association with rs904419 (P = 4.34 × 10⁻⁸), located between the FRMD4B and MITF genes. The median IC50 values of prednisolone for cell lines with rs904419 AA (n = 13), AG (n = 31) and GG (n = 28) genotypes were 0.089, 0.139 and 297 µmol/L, respectively. For dexamethasone sensitivity, suggestive association was observed for SNP rs2306888 (P = 1.43 × 10⁻⁶), a synonymous SNP of the TGFBR3 gene. For NR3C1 gene expression, suggestive association was observed for SNP rs11982167 (P = 6.44 × 10⁻⁸), located in the PLEKHA8 gene. These genetic variants may affect GC sensitivity of ALL cells and may give rise to opportunities in personalized medicine for effective and safe chemotherapy in ALL patients.     

Rare coding variation in paraoxonase-1 is associated with ischemic stroke in the NHLBI Exome Sequencing Project

HDL-associated paraoxonase-1 (PON1) is an enzyme whose activity is associated with cerebrovascular disease. Common PON1 genetic variants have not been consistently associated with cerebrovascular disease. Rare coding variation that likely alters PON1 enzyme function may be more strongly associated with stroke. The National Heart, Lung, and Blood Institute Exome Sequencing Project sequenced the coding regions (exomes) of the genome for heart, lung, and blood-related phenotypes (including ischemic stroke). In this sample of 4,204 unrelated participants, 496 had verified, noncardioembolic ischemic stroke. After filtering, 28 nonsynonymous PON1 variants were identified. Analysis with the sequence kernel association test, adjusted for covariates, identified significant associations between PON1 variants and ischemic stroke (P = 3.01 × 10⁻³). Stratified analyses demonstrated a stronger association of PON1 variants with ischemic stroke in African ancestry (AA) participants (P = 5.03 × 10⁻³). Ethnic differences in the association between PON1 variants with stroke could be due to the effects of PON1_Val109Ile (overall P = 7.88 × 10⁻³; AA P = 6.52 × 10⁻⁴), found at higher frequency in AA participants (1.16% vs. 0.02%) and whose protein is less stable than the common allele. In summary, rare genetic variation in PON1 was associated with ischemic stroke, with stronger associations identified in those of AA. Increased focus on PON1 enzyme function and its role in cerebrovascular disease is warranted.     

Genetic Variants Associated with Lipid Profiles in Chinese Patients with Type 2 Diabetes

Dyslipidemia is a strong risk factor for cardiovascular disease among patients with type 2 diabetes (T2D). The aim of this study was to identify lipid-related genetic variants in T2D patients of Han Chinese ancestry. Among 4,908 Chinese T2D patients who were not taking lipid-lowering medications, single nucleotide polymorphisms (SNPs) in seven genes previously found to be associated with lipid traits in genome-wide association studies conducted in populations of European ancestry (ABCA1, GCKR, BAZ1B, TOMM40, DOCK7, HNF1A, and HNF4A) were genotyped. After adjusting for multiple covariates, SNPs in ABCA1, GCKR, BAZ1B, TOMM40, and HNF1A were identified as significantly associated with triglyceride levels in T2D patients (P < 0.05). The associations between the SNPs in ABCA1 (rs3890182), GCKR (rs780094), and BAZ1B (rs2240466) remained significant even after correction for multiple testing (P = 8.85×10⁻³, 7.88×10⁻⁷, and 2.03×10⁻⁶, respectively). BAZ1B (rs2240466) also was associated with the total cholesterol level (P = 4.75×10⁻²). In addition, SNP rs157580 in TOMM40 was associated with the low-density lipoprotein cholesterol level (P = 6.94×10⁻³). Our findings confirm that lipid-related genetic loci are associated with lipid profiles in Chinese patients with type 2 diabetes.     

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.     

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.     

Systematic Evaluation of Genetic Variants for Polycystic Ovary Syndrome in a Chinese Population

To date, eleven genome-wide significant (GWS) loci (P < 5×10⁻⁸) for polycystic ovary syndrome (PCOS) have been identified through genome-wide association studies (GWAS). Some of the risk loci have been selected for replications and validated in multiple ethnicities, however, few previous studies investigated all loci. Scanning all the GWAS variants would demonstrate a more informative profile of variance they explained. Thus, we analyzed all the 17 single nucleotide polymorphisms (SNPs) mapping to the 11 GWAS loci in an independent sample set of 800 Chinese subjects with PCOS and 1110 healthy controls systematically. Variants of rs3802457 in C9orf3 locus (P = 5.99×10⁻⁴) and rs13405728 in LHCGR locus (P = 3.73×10⁻⁴) were significantly associated with PCOS after the strict Bonferroni correction in our data set. The further haplotype analysis indicated that in the block of C9orf3 gene (rs4385527 and rs3802457), GA haplotype played a protective role in PCOS (8.7 vs 5.0, P = 9.85×10⁻⁶, OR = 0.548, 95%CI = 0.418–0.717), while GG haplotype was found suffering from an extraordinarily increased risk of PCOS (73.6% vs79.2%, P = 3.41×10⁻⁵, OR = 1.394, 95%CI = 1.191–1.632). Moreover, the directions of effects for all SNPs were consistent with previous GWAS reports (P = 1.53×10⁻⁵). Polygenic score analysis demonstrated that these 17 SNPs have a significant capacity on predicting case-control status in our samples (P = 7.17×10⁻⁹), meanwhile all these gathered 17 SNPs explained about 2.40% of variance. Our findings supported that C9orf3 and LHCGR loci variants were vital susceptibility of PCOS.     

A Genome-Wide Association Study Uncovers a Genetic Locus Associated with Thoracic-to-Hip Ratio in Koreans

The thoracic-to-hip circumference ratio (THR) is an anthropometric marker recently described as a predictor of type 2 diabetes. In this study, we performed a genome-wide association study (GWAS) followed by confirmatory analyses to identify genetic markers associated with THR. A total of 7,240 Korean subjects (4,988 for the discovery stage and 2,252 for the confirmatory analyses) were recruited for this study, and genome-wide single nucleotide polymorphism (SNP) genotyping of the initial 4,988 individuals was performed using Affymetrix Human SNP array 5.0. Linear regression analysis was then performed to adjust for the effects of age, sex, and current diabetes medication status on the THR of the study subjects. In the initial discovery stage, there was a statistically nominal association between minor alleles of SNP markers on chromosomes 4, 8, 10, and 12, and THR changes (p < 5.0 × 10⁻⁶). The subsequent confirmatory analyses of these markers, however, only detected a significant association between two SNPs in the HECTD4 gene and decreased THRs. Notably, this association was detected in male (rs11066280: p = 1.14 × 10⁻²; rs2074356: p = 1.10 × 10⁻²), but not in female subjects. Meanwhile, the combined results from the two analyses (initial and confirmatory) indicated that minor alleles of these two intronic variants exhibited a significant genome-wide association with decreased THR in the male subjects (n = 3,155; rs11066280: effect size = −0.008624, p = 6.19 × 10⁻⁹; rs2074356: effect size = −0.008762, p = 1.89 × 10⁻⁸). Furthermore, minor alleles of these two SNPs exhibited protective effects on patients’ risks for developing type 2 diabetes. In conclusion, we have identified two genetic variations in HECTD4 that are associated with THR, particularly in men.     

Genetic Determinants of Long-Term Changes in Blood Lipid Concentrations: 10-Year Follow-Up of the GLACIER Study

Recent genome-wide meta-analyses identified 157 loci associated with cross-sectional lipid traits. Here we tested whether these loci associate (singly and in trait-specific genetic risk scores [GRS]) with longitudinal changes in total cholesterol (TC) and triglyceride (TG) levels in a population-based prospective cohort from Northern Sweden (the GLACIER Study). We sought replication in a southern Swedish cohort (the MDC Study; N = 2,943). GLACIER Study participants (N = 6,064) were genotyped with the MetaboChip array. Up to 3,495 participants had 10-yr follow-up data available in the GLACIER Study. The TC- and TG-specific GRSs were strongly associated with change in lipid levels (β = 0.02 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻¹¹ for TC; β = 0.02 mmol/l per effect allele per decade follow-up, P = 5.0×10⁻⁵ for TG). In individual SNP analysis, one TC locus, apolipoprotein E (APOE) rs4420638 (β = 0.12 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻⁵), and two TG loci, tribbles pseudokinase 1 (TRIB1) rs2954029 (β = 0.09 mmol/l per effect allele per decade follow-up, P = 5.1×10⁻⁴) and apolipoprotein A-I (APOA1) rs6589564 (β = 0.31 mmol/l per effect allele per decade follow-up, P = 1.4×10⁻⁸), remained significantly associated with longitudinal changes for the respective traits after correction for multiple testing. An additional 12 loci were nominally associated with TC or TG changes. In replication analyses, the APOE rs4420638, TRIB1 rs2954029, and APOA1 rs6589564 associations were confirmed (P≤0.001). In summary, trait-specific GRSs are robustly associated with 10-yr changes in lipid levels and three individual SNPs were strongly associated with 10-yr changes in lipid levels.     

FADS genetic variants and ω-6 polyunsaturated fatty acid metabolism in a homogeneous island population

Long-chain polyunsaturated fatty acids (PUFA) orchestrate immunity and inflammation through their capacity to be converted to potent inflammatory mediators. We assessed associations of FADS gene cluster polymorphisms and fasting serum PUFA concentrations in a fully ascertained, geographically isolated founder population of European descent. Concentrations of 22 PUFAs were determined by gas chromatography, of which ten fatty acids and five ratios defining FADS1 and FADS2 activity were tested for genetic association against 16 single nucleotide polymorphisms (SNP) in 224 individuals. A cluster of SNPs in tight linkage disequilibrium in the FADS1 gene (rs174537, rs174545, rs174546, rs174553, rs174556, rs174561, rs174568, and rs99780) were strongly associated with arachidonic acid (AA) (P = 5.8 × 10⁻⁷ – 1.7 × 10⁻⁸) among other PUFAs, but the strongest associations were with the ratio measuring FADS1 activity in the ω-6 series (P = 2.11 × 10⁻¹³ – 1.8 × 10⁻²⁰). The minor allele across all SNPs was consistently associated with decreased ω-6 PUFAs, with the exception of dihomo-γ-linoleic acid (DHGLA), where the minor allele was consistently associated with increased levels. Our findings in a geographically isolated population with a homogenous dietary environment suggest that variants in the Δ-5 desaturase enzymatic step likely regulate the efficiency of conversion of medium-chain PUFAs to potentially inflammatory PUFAs, such as AA.     

Association analyses confirm five susceptibility loci for systemic lupus erythematosus in the Han Chinese population

IntroductionSystemic lupus erythematosus (SLE) is a multisystem autoimmune disease. Currently, numerous genetic loci of SLE have been confirmed. Here we try to further explore additional genes contributing to SLE susceptibility in this study.MethodsForty nine single nucleotide polymorphisms (SNPs) with moderate-risk for SLE in previous study were genotyped in a large-scale replication study with a total of 3,522 cases and 8,252 controls using the Sequenom Massarray system. Association analyses were performed using logistic regression with gender or sample cohorts as a covariate through PLINK 1.07 software.ResultsThis replication effort confirmed five reported SLE susceptibility loci reaching genome-wide levels of significance (P_meta <5.00 × 10⁻⁰⁸): TNFSF4 (rs1418190, odds ratio (OR) = 0.81, P_meta = 1.08 × 10⁻⁰⁸; rs4916219, OR = 0.80, P_meta = 7.77 × 10⁻⁰⁹), IRF8 (rs2934498, OR = 1.25, P_meta = 4.97 × 10⁻⁰⁹), miR-146a (rs2431697, OR = 0.69, P_meta = 1.15 × 10⁻²²), CD44 (rs2732547, OR = 0.82, P_meta = 1.55 × 10⁻¹¹), and TMEM39A (rs12494314, OR = 0.84, P_meta = 1.01 × 10⁻⁰⁹). Further logistic regression analysis indicated that the genetic effects within TNFSF4 detected in this study are independent from our previously reported signals.ConclusionsThis study increases the number of established susceptibility loci for SLE in Han Chinese population and highlights the contribution of multiple variants of modest effect. Although further studies will be required to identify the causal alleles within these loci, the findings make a significant step forward in our understanding of the genetic contribution to SLE in Chinese population.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0602-9) contains supplementary material, which is available to authorized users.     

Identification of an Interaction between VWF rs7965413 and Platelet Count as a Novel Risk Marker for Metabolic Syndrome: An Extensive Search of Candidate Polymorphisms in a Case-Control Study

Although many single nucleotide polymorphisms (SNPs) have been identified to be associated with metabolic syndrome (MetS), there was only a slight improvement in the ability to predict future MetS by the simply addition of SNPs to clinical risk markers. To improve the ability to predict future MetS, combinational effects, such as SNP—SNP interaction, SNP—environment interaction, and SNP—clinical parameter (SNP × CP) interaction should be also considered. We performed a case-control study to explore novel SNP × CP interactions as risk markers for MetS based on health check-up data of Japanese male employees. We selected 99 SNPs that were previously reported to be associated with MetS and components of MetS; subsequently, we genotyped these SNPs from 360 cases and 1983 control subjects. First, we performed logistic regression analyses to assess the association of each SNP with MetS. Of these SNPs, five SNPs were significantly associated with MetS (P < 0.05): LRP2 rs2544390, rs1800592 between UCP1 and TBC1D9, APOA5 rs662799, VWF rs7965413, and rs1411766 between MYO16 and IRS2. Furthermore, we performed multiple logistic regression analyses, including an SNP term, a CP term, and an SNP × CP interaction term for each CP and SNP that was significantly associated with MetS. We identified a novel SNP × CP interaction between rs7965413 and platelet count that was significantly associated with MetS [SNP term: odds ratio (OR) = 0.78, P = 0.004; SNP × CP interaction term: OR = 1.33, P = 0.001]. This association of the SNP × CP interaction with MetS remained nominally significant in multiple logistic regression analysis after adjustment for either the number of MetS components or MetS components excluding obesity. Our results reveal new insight into platelet count as a risk marker for MetS.     

A Large-Scale Genetic Analysis Reveals a Strong Contribution of the HLA Class II Region to Giant Cell Arteritis Susceptibility

We conducted a large-scale genetic analysis on giant cell arteritis (GCA), a polygenic immune-mediated vasculitis. A case-control cohort, comprising 1,651 case subjects with GCA and 15,306 unrelated control subjects from six different countries of European ancestry, was genotyped by the Immunochip array. We also imputed HLA data with a previously validated imputation method to perform a more comprehensive analysis of this genomic region. The strongest association signals were observed in the HLA region, with rs477515 representing the highest peak (p = 4.05 × 10⁻⁴⁰, OR = 1.73). A multivariate model including class II amino acids of HLA-DRβ1 and HLA-DQα1 and one class I amino acid of HLA-B explained most of the HLA association with GCA, consistent with previously reported associations of classical HLA alleles like HLA-DRB1^∗04. An omnibus test on polymorphic amino acid positions highlighted DRβ1 13 (p = 4.08 × 10⁻⁴³) and HLA-DQα1 47 (p = 4.02 × 10⁻⁴⁶), 56, and 76 (both p = 1.84 × 10⁻⁴⁵) as relevant positions for disease susceptibility. Outside the HLA region, the most significant loci included PTPN22 (rs2476601, p = 1.73 × 10⁻⁶, OR = 1.38), LRRC32 (rs10160518, p = 4.39 × 10⁻⁶, OR = 1.20), and REL (rs115674477, p = 1.10 × 10⁻⁵, OR = 1.63). Our study provides evidence of a strong contribution of HLA class I and II molecules to susceptibility to GCA. In the non-HLA region, we confirmed a key role for the functional PTPN22 rs2476601 variant and proposed other putative risk loci for GCA involved in Th1, Th17, and Treg cell function.     

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⁻⁷, odds ratio [OR] 1.42) and identify a novel candidate gene, ZBTB10, achieving genome-wide significance (rs6998967, P = 8.9 × 10⁻¹⁰, OR 0.53). Several other SNPs showed suggestive significance including rs2476601 (P = 6.5 × 10⁻⁶, 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⁻¹²) versus 2.82 in EOMG (P = 3.86 × 10⁻⁴⁵). 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.     

Genetic Analysis Identifies DDR2 as a Novel Gene Affecting Bone Mineral Density and Osteoporotic Fractures in Chinese Population

DDR2 gene, playing an essential role in regulating osteoblast differentiation and chondrocyte maturation, may influence bone mineral density (BMD) and osteoporosis, but the genetic variations actually leading to the association remain to be elucidated. Therefore, the aim of this study was to investigate whether the genetic variants in DDR2 are associated with BMD and fracture risk. This study was performed in three samples from two ethnicities, including 1,300 Chinese Han subjects, 700 Chinese Han subjects (350 with osteoporotic hip fractures and 350 healthy controls) and 2,286 US white subjects. Twenty-eight SNPs in DDR2 were genotyped and tested for associations with hip BMD and fractures. We identified 3 SNPs in DDR2 significantly associated with hip BMD in the Chinese population after multiple testing adjustments, which were rs7521233 (P = 1.06×10⁻⁴, β: −0.018 for allele C), rs7553831 (P = 1.30×10⁻⁴, β: −0.018 for allele T), and rs6697469 (P = 1.59×10⁻³, β: −0.015 for allele C), separately. These three SNPs were in high linkage disequilibrium. Haplotype analyses detected two significantly associated haplotypes, including one haplotype in block 2 (P = 9.54×10⁻⁴, β: −0.016) where these three SNPs located. SNP rs6697469 was also associated with hip fractures (P = 0.043, OR: 1.42) in the Chinese population. The effect on fracture risk was consistent with its association with lower BMD. However, in the white population, we didn’t observe significant associations with hip BMD. eQTL analyses revealed that SNPs associated with BMD also affected DDR2 mRNA expression levels in Chinese. Our findings, together with the prior biological evidence, suggest that DDR2 could be a new candidate for osteoporosis in Chinese population. Our results also reveal an ethnic difference, which highlights the need for further genetic studies in each ethnic group.