The SNP rs3128965 of HLA-DPB1 as a Genetic Marker of the AERD Phenotype

Background

Two common clinical syndromes of acetylsalicylic acid (aspirin) hypersensitivity, aspirin-exacerbated respiratory disease (AERD) and aspirin-exacerbated cutaneous disease (AECD), were subjected to a genome-wide association study to identify strong genetic markers for aspirin hypersensitivity in a Korean population.

Methods

A comparison of SNP genotype frequencies on an Affymetrix Genome-Wide Human SNP array of 179 AERD patients and 1989 healthy normal control subjects (NC) revealed SNPs on chromosome 6 that were associated with AERD, but not AECD. To validate the association, we enrolled a second cohort comprising AERD (n = 264), NC (n = 238) and disease-control (aspirin tolerant asthma; ATA, n = 387) groups.

Results

The minor genotype frequency (AG or AA) of a particular SNP, rs3128965, in the HLA-DPB1 region was higher in the AERD group compared to the ATA or NC group (P = 0.001, P = 0.002, in a co-dominant analysis model, respectively). Comparison of rs3128965 alleles with the clinical features of asthmatics revealed that patients harboring the A allele had increased bronchial hyperresponsiveness to inhaled aspirin and methacholine, and higher 15-HETE levels, than those without the A allele (P = 0.039, 0.037, and 0.004, respectively).

Conclusions

This implies the potential of rs3128965 as a genetic marker for diagnosis and prediction of the AERD phenotype.     

Elevation of Eosinophil-Derived Neurotoxin in Plasma of the Subjects with Aspirin-Exacerbated Respiratory Disease: A Possible Peripheral Blood Protein Biomarker

Aspirin-exacerbated respiratory disease (AERD) remains widely underdiagnosed in asthmatics, primarily due to insufficient awareness of the relationship between aspirin ingestion and asthma exacerbation. The identification of aspirin hypersensitivity is therefore essential to avoid serious aspirin complications. The goal of the study was to develop plasma biomarkers to predict AERD. We identified differentially expressed genes in peripheral blood mononuclear cells (PBMC) between subjects with AERD and those with aspirin-tolerant asthma (ATA). The genes were matched with the secreted protein database (http://spd.cbi.pku.edu.cn/) to select candidate proteins in the plasma. Plasma levels of the candidate proteins were then measured in AERD (n = 40) and ATA (n = 40) subjects using an enzyme-linked immunosorbent assay (ELISA). Target genes were validated as AERD biomarkers using an ROC curve analysis. From 175 differentially expressed genes (p-value <0.0001) that were queried to the secreted protein database, 11 secreted proteins were retrieved. The gene expression patterns were predicted as elevated for 7 genes and decreased for 4 genes in AERD as compared with ATA subjects. Among these genes, significantly higher levels of plasma eosinophil-derived neurotoxin (RNASE2) were observed in AERD as compared with ATA subjects (70(14.62∼311.92) µg/ml vs. 12(2.55∼272.84) µg/ml, p-value <0.0003). Based on the ROC curve analysis, the AUC was 0.74 (p-value = 0.0001, asymptotic 95% confidence interval [lower bound: 0.62, upper bound: 0.83]) with 95% sensitivity, 60% specificity, and a cut-off value of 27.15 µg/ml. Eosinophil-derived neurotoxin represents a novel biomarker to distinguish AERD from ATA.     

Modeling of Environmental Effects in Genome-Wide Association Studies Identifies SLC2A2 and HP as Novel Loci Influencing Serum Cholesterol Levels

Genome-wide association studies (GWAS) have identified 38 larger genetic regions affecting classical blood lipid levels without adjusting for important environmental influences. We modeled diet and physical activity in a GWAS in order to identify novel loci affecting total cholesterol, LDL cholesterol, HDL cholesterol, and triglyceride levels. The Swedish (SE) EUROSPAN cohort (N _SE = 656) was screened for candidate genes and the non-Swedish (NS) EUROSPAN cohorts (N _NS = 3,282) were used for replication. In total, 3 SNPs were associated in the Swedish sample and were replicated in the non-Swedish cohorts. While SNP rs1532624 was a replication of the previously published association between CETP and HDL cholesterol, the other two were novel findings. For the latter SNPs, the p-value for association was substantially improved by inclusion of environmental covariates: SNP rs5400 (p _{SE,unadjusted} = 3.6×10⁻⁵, p _SE,adjusted = 2.2×10⁻⁶, p _{NS,unadjusted} = 0.047) in the SLC2A2 (Glucose transporter type 2) and rs2000999 (p _{SE,unadjusted} = 1.1×10⁻³, p _SE,adjusted = 3.8×10⁻⁴, p _{NS,unadjusted} = 0.035) in the HP gene (Haptoglobin-related protein precursor). Both showed evidence of association with total cholesterol. These results demonstrate that inclusion of important environmental factors in the analysis model can reveal new genetic susceptibility loci.     

A Major Histocompatibility Class I Locus Contributes to Multiple Sclerosis Susceptibility Independently from HLA-DRB1*15:01

Background

In Northern European descended populations, genetic susceptibility for multiple sclerosis (MS) is associated with alleles of the human leukocyte antigen (HLA) Class II gene DRB1. Whether other major histocompatibility complex (MHC) genes contribute to MS susceptibility is controversial.

Methodology/Principal Findings

A case control analysis was performed using 958 single nucleotide polymorphisms (SNPs) spanning the MHC assayed in two independent datasets. The discovery dataset consisted of 1,018 cases and 1,795 controls and the replication dataset was composed of 1,343 cases and 1,379 controls. The most significantly MS-associated SNP in the discovery dataset was rs3135391, a Class II SNP known to tag the HLA-DRB1*15:01 allele, the primary MS susceptibility allele in the MHC (O.R. = 3.04, p<1×10⁻⁷⁸). To control for the effects of the HLA-DRB1*15:01 haplotype, case control analysis was performed adjusting for this HLA-DRB1*15:01 tagging SNP. After correction for multiple comparisons (false discovery rate = .05) 52 SNPs in the Class I, II and III regions were significantly associated with MS susceptibility in both datasets using the Cochran Armitage trend test. The discovery and replication datasets were merged and subjects carrying the HLA-DRB1*15:01 tagging SNP were excluded. Association tests showed that 48 of the 52 replicated SNPs retained significant associations with MS susceptibility independently of the HLA-DRB1*15:01 as defined by the tagging SNP. 20 Class I SNPs were associated with MS susceptibility with p-values ≤1×10⁻⁸. The most significantly associated SNP was rs4959039, a SNP in the downstream un-translated region of the non-classical HLA-G gene (Odds ratio 1.59, 95% CI 1.40, 1.81, p = 8.45×10⁻¹³) and is in linkage disequilibrium with several nearby SNPs. Logistic regression modeling showed that this SNP''s contribution to MS susceptibility was independent of the Class II and Class III SNPs identified in this screen.

Conclusions

A MHC Class I locus contributes to MS susceptibility independently of the HLA-DRB1*15:01 haplotype.     

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.     

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.     

Common Polymorphisms Influencing Serum Uric Acid Levels Contribute to Susceptibility to Gout,but Not to Coronary Artery Disease

Background

Recently, a large meta-analysis including over 28,000 participants identified nine different loci with association to serum uric acid (UA) levels. Since elevated serum UA levels potentially cause gout and are a possible risk factor for coronary artery disease (CAD) and myocardial infarction (MI), we performed two large case-control association analyses with participants from the German MI Family Study. In the first study, we assessed the association of the qualitative trait gout and ten single nucleotide polymorphisms (SNP) markers that showed association to UA serum levels. In the second study, the same genetic polymorphisms were analyzed for association with CAD.

Methods and Findings

A total of 683 patients suffering from gout and 1,563 healthy controls from the German MI Family Study were genotyped. Nine SNPs were identified from a recently performed genome-wide meta-analysis on serum UA levels (rs12129861, rs780094, rs734553, rs2231142, rs742132, rs1183201, rs12356193, rs17300741 and rs505802). Additionally, the marker rs6855911 was included which has been associated with gout in our cohort in a previous study. SNPs rs734553 and rs6855911, located in SLC2A9, and SNP rs2231142, known to be a missense polymorphism in ABCG2, were associated with gout (p = 5.6*10⁻⁷, p = 1.1*10⁻⁷, and p = 1.3*10⁻³, respectively). Other SNPs in the genes PDZK1, GCKR, LRRC16A, SLC17A1-SLC17A3, SLC16A9, SLC22A11 and SLC22A12 failed the significance level. None of the ten markers were associated with risk to CAD in our study sample of 1,473 CAD cases and 1,241 CAD-free controls.

Conclusion

SNP markers in SLC2A9 and ABCG2 genes were found to be strongly associated with the phenotype gout. However, not all SNP markers influencing serum UA levels were also directly associated with the clinical manifestation of gout in our study sample. In addition, none of these SNPs showed association with the risk to CAD in the German MI Family Study.     

Identification of a Functional Genetic Variant at 16q12.1 for Breast Cancer Risk: Results from the Asia Breast Cancer Consortium

Genetic factors play an important role in the etiology of breast cancer. We carried out a multi-stage genome-wide association (GWA) study in over 28,000 cases and controls recruited from 12 studies conducted in Asian and European American women to identify genetic susceptibility loci for breast cancer. After analyzing 684,457 SNPs in 2,073 cases and 2,084 controls in Chinese women, we evaluated 53 SNPs for fast-track replication in an independent set of 4,425 cases and 1,915 controls of Chinese origin. Four replicated SNPs were further investigated in an independent set of 6,173 cases and 6,340 controls from seven other studies conducted in Asian women. SNP rs4784227 was consistently associated with breast cancer risk across all studies with adjusted odds ratios (95% confidence intervals) of 1.25 (1.20−1.31) per allele (P = 3.2×10⁻²⁵) in the pooled analysis of samples from all Asian samples. This SNP was also associated with breast cancer risk among European Americans (per allele OR  = 1.19, 95% CI  = 1.09−1.31, P = 1.3×10⁻⁴, 2,797 cases and 2,662 controls). SNP rs4784227 is located at 16q12.1, a region identified previously for breast cancer risk among Europeans. The association of this SNP with breast cancer risk remained highly statistically significant in Asians after adjusting for previously-reported SNPs in this region. In vitro experiments using both luciferase reporter and electrophoretic mobility shift assays demonstrated functional significance of this SNP. These results provide strong evidence implicating rs4784227 as a functional causal variant for breast cancer in the locus 16q12.1 and demonstrate the utility of conducting genetic association studies in populations with different genetic architectures.     

Extracellular Water and Blood Pressure in Adults with Growth Hormone (GH) Deficiency: A Genotype-Phenotype Association Study

Objectives

Growth hormone deficiency (GHD) in adults is associated with decreased extracellular water volume (ECW). In response to GH replacement therapy (GHRT), ECW increases and blood pressure (BP) reduces or remains unchanged. Our primary aim was to study the association between polymorphisms in genes related to renal tubular function with ECW and BP before and 1 year after GHRT. The ECW measures using bioimpedance analysis (BIA) and bioimpedance spectroscopy (BIS) were validated against a reference method, the sodium bromide dilution method (Br⁻).

Design and Methods

Using a candidate gene approach, fifteen single-nucleotide polymorphisms (SNPs) in nine genes with known impact on renal tubular function (AGT, SCNN1A, SCNN1G, SLC12A1, SLC12A3, KCNJ1, STK39, WNK1 and CASR) were genotyped and analyzed for associations with ECW and BP at baseline and with their changes after 1 year of GHRT in 311 adult GHD patients. ECW was measured with the Br⁻, BIA, and BIS.

Results

Both BIA and BIS measurements demonstrated similar ECW results as the reference method. At baseline, after adjustment for sex and BMI, SNP rs2291340 in the SLC12A1 gene was associated with ECW volume in GHD patients (p = 0.039). None of the SNPs influenced the ECW response to GHRT. One SNP in the SLC12A3 gene (rs11643718; p = 0.024) and three SNPs in the SCNN1G gene [rs5723 (p = 0.02), rs5729 (p = 0.016) and rs13331086 (p = 0.035)] were associated with the inter-individual differences in BP levels at baseline. A polymorphism in the calcium-sensing receptor (CASR) gene (rs1965357) was associated with changes in systolic BP after GHRT (p = 0.036). None of these associations remained statistically significant when corrected for multiple testing.

Conclusion

The BIA and BIS are as accurate as Br⁻ to measure ECW in GHD adults before and during GHRT. Our study provides the first evidence that individual polymorphisms may have clinically relevant effects on ECW and BP in GHD adults.     

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.     

Fine Mapping and Functional Analysis of the Multiple Sclerosis Risk Gene CD6

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2^nd SRCR domain with susceptibility to MS (P _{max(T) permutation} = 1×10⁻⁴). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. – CD4⁺ naïve cells, P = 0.0001; CD8⁺ naïve cells, P<0.0001; CD4⁺ and CD8⁺ central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4⁺ and CD8⁺ T cells.     

Evaluation of 41 Candidate Gene Variants for Obesity in the EPIC-Potsdam Cohort by Multi-Locus Stepwise Regression

Objective

Obesity has become a leading preventable cause of morbidity and mortality in many parts of the world. It is thought to originate from multiple genetic and environmental determinants. The aim of the current study was to introduce haplotype-based multi-locus stepwise regression (MSR) as a method to investigate combinations of unlinked single nucleotide polymorphisms (SNPs) for obesity phenotypes.

Methods

In 2,122 healthy randomly selected men and women of the EPIC-Potsdam cohort, the association between 41 SNPs from 18 obesity-candidate genes and either body mass index (BMI, mean = 25.9 kg/m², SD = 4.1) or waist circumference (WC, mean = 85.2 cm, SD = 12.6) was assessed. Single SNP analyses were done by using linear regression adjusted for age, sex, and other covariates. Subsequently, MSR was applied to search for the ‘best’ SNP combinations. Combinations were selected according to specific AIC_c and p-value criteria. Model uncertainty was accounted for by a permutation test.

Results

The strongest single SNP effects on BMI were found for TBC1D1 rs637797 (β = −0.33, SE = 0.13), FTO rs9939609 (β = 0.28, SE = 0.13), MC4R rs17700144 (β = 0.41, SE = 0.15), and MC4R rs10871777 (β = 0.34, SE = 0.14). All these SNPs showed similar effects on waist circumference. The two ‘best’ six-SNP combinations for BMI (global p-value = 3.45⋅10^–6 and 6.82⋅10^–6) showed effects ranging from −1.70 (SE = 0.34) to 0.74 kg/m² (SE = 0.21) per allele combination. We selected two six-SNP combinations on waist circumference (global p-value = 7.80⋅10^–6 and 9.76⋅10^–6) with an allele combination effect of −2.96 cm (SE = 0.76) at maximum. Additional adjustment for BMI revealed 15 three-SNP combinations (global p-values ranged from 3.09⋅10^–4 to 1.02⋅10^–2). However, after carrying out the permutation test all SNP combinations lost significance indicating that the statistical associations might have occurred by chance.

Conclusion

MSR provides a tool to search for risk-related SNP combinations of common traits or diseases. However, the search process does not always find meaningful SNP combinations in a dataset.     

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.     

Novel Genetic Loci Identified for the Pathophysiology of Childhood Obesity in the Hispanic Population

Genetic variants responsible for susceptibility to obesity and its comorbidities among Hispanic children have not been identified. The VIVA LA FAMILIA Study was designed to genetically map childhood obesity and associated biological processes in the Hispanic population. A genome-wide association study (GWAS) entailed genotyping 1.1 million single nucleotide polymorphisms (SNPs) using the Illumina Infinium technology in 815 children. Measured genotype analysis was performed between genetic markers and obesity-related traits i.e., anthropometry, body composition, growth, metabolites, hormones, inflammation, diet, energy expenditure, substrate utilization and physical activity. Identified genome-wide significant loci: 1) corroborated genes implicated in other studies (MTNR1B, ZNF259/APOA5, XPA/FOXE1 (TTF-2), DARC, CCR3, ABO); 2) localized novel genes in plausible biological pathways (PCSK2, ARHGAP11A, CHRNA3); and 3) revealed novel genes with unknown function in obesity pathogenesis (MATK, COL4A1). Salient findings include a nonsynonymous SNP (rs1056513) in INADL (p = 1.2E-07) for weight; an intronic variant in MTNR1B associated with fasting glucose (p = 3.7E-08); variants in the APOA5-ZNF259 region associated with triglycerides (p = 2.5-4.8E-08); an intronic variant in PCSK2 associated with total antioxidants (p = 7.6E-08); a block of 23 SNPs in XPA/FOXE1 (TTF-2) associated with serum TSH (p = 5.5E-08 to 1.0E-09); a nonsynonymous SNP (p = 1.3E-21), an intronic SNP (p = 3.6E-13) in DARC identified for MCP-1; an intronic variant in ARHGAP11A associated with sleep duration (p = 5.0E-08); and, after adjusting for body weight, variants in MATK for total energy expenditure (p = 2.7E-08) and in CHRNA3 for sleeping energy expenditure (p = 6.0E-08). Unprecedented phenotyping and high-density SNP genotyping enabled localization of novel genetic loci associated with the pathophysiology of childhood obesity.     

Genome-Wide Scan Identifies Variant in TNFSF13 Associated with Serum IgM in a Healthy Chinese Male Population

IgM provides a first line of defense during microbial infections. Serum IgM levels are detected routinely in clinical practice. And IgM is a genetically complex trait. We conducted a two-stage genome-wide association study (GWAS) to identify genetic variants affecting serum IgM levels in a Chinese population of 3495, including 1999 unrelated subjects in the first stage and 1496 independent individuals in the second stage. Our data show that a common single nucleotide polymorphism (SNP), rs11552708 located in the TNFSF13 gene was significantly associated with IgM levels (p = 5.00×10⁻⁷ in first stage, p = 1.34×10⁻³ in second stage, and p = 4.22×10⁻⁹ when combined). Besides, smoking was identified to be associated with IgM levels in both stages (P<0.05), but there was no significant interaction between smoking and the identified SNP (P>0.05). It is suggested that TNFSF13 may be a susceptibility gene affecting serum IgM levels in Chinese male population.     

Results of Genome-Wide Analyses on Neurodevelopmental Phenotypes at Four-Year Follow-Up following Cardiac Surgery in Infancy

Background

Adverse neurodevelopmental sequelae are reported among children who undergo early cardiac surgery to repair congenital heart defects (CHD). APOE genotype has previously been determined to contribute to the prediction of these outcomes. Understanding further genetic causes for the development of poor neurobehavioral outcomes should enhance patient risk stratification and improve both prevention and treatment strategies.

Methods

We performed a prospective observational study of children who underwent cardiac surgery before six months of age; this included a neurodevelopmental evaluation between their fourth and fifth birthdays. Attention and behavioral skills were assessed through parental report utilizing the Attention Deficit-Hyperactivity Disorder-IV scale preschool edition (ADHD-IV), and Child Behavior Checklist (CBCL/1.5-5), respectively. Of the seven investigated, three neurodevelopmental phenotypes met genomic quality control criteria. Linear regression was performed to determine the effect of genome-wide genetic variation on these three neurodevelopmental measures in 316 subjects.

Results

This genome-wide association study identified single nucleotide polymorphisms (SNPs) associated with three neurobehavioral phenotypes in the postoperative children ADHD-IV Impulsivity/Hyperactivity, CBCL/1.5-5 PDPs, and CBCL/1.5-5 Total Problems. The most predictive SNPs for each phenotype were: a LGALS8 intronic SNP, rs4659682, associated with ADHD-IV Impulsivity (P = 1.03×10⁻⁶); a PCSK5 intronic SNP, rs2261722, associated with CBCL/1.5-5 PDPs (P = 1.11×10⁻⁶); and an intergenic SNP, rs11617488, 50 kb from FGF9, associated with CBCL/1.5-5 Total Problems (P = 3.47×10⁻⁷). 10 SNPs (3 for ADHD-IV Impulsivity, 5 for CBCL/1.5-5 PDPs, and 2 for CBCL/1.5-5 Total Problems) had p<10⁻⁵.

Conclusions

No SNPs met genome-wide significance for our three neurobehavioral phenotypes; however, 10 SNPs reached a threshold for suggestive significance (p<10⁻⁵). Given the unique nature of this cohort, larger studies and/or replication are not possible. Studies to further investigate the mechanisms through which these newly identified genes may influence neurodevelopment dysfunction are warranted.     

Association of Single Nucleotide Polymorphisms in the Lens Epithelium-Derived Growth Factor (LEDGF/p75) with HIV-1 Infection Outcomes in Brazilian HIV-1+ Individuals

The lens epithelium-derived growth factor p75 (LEDGF/p75), coded by the PSIP1 gene, is an important host co-factor that interacts with HIV-1 integrase to target integration of viral cDNA into active genes. The aim of this study was to investigate the association of SNPs in the PSIP1 gene with disease outcome in HIV-1 infected patients. We performed a genetic association study in a cohort of 171 HIV-1 seropositive Brazilian individuals classified as rapid progressors (RP, n = 69), typical progressors (TP, n = 79) and long-term nonprogressors (LTNP, n = 23). The exonic SNP rs61744944 and 9 tag SNPs were genotyped. A group of 192 healthy subjects was analyzed to determine the frequency of SNPs and haplotypes in the general population. Linkage disequilibrium (LD) analyses indicated that the SNPs analyzed were not in high LD (r²<0.8). Logistic regression models suggested that patients carrying the T allele rs61744944 (472L) were more likely to develop a LTNP phenotype (OR = 4.98; p = 0.05) as compared to TP group. The same trend was observed when LTNPs were compared to the RP group (OR = 3.26). Results of haplotype analyses reinforced this association, since the OR values obtained for the haplotype carrying allele T at rs61744944 also reflected an association with LTNP status (OR = 6.05; p = 0.08 and OR = 3.44; p = 0.12 for comparisons to TP and RP, respectively). The rare missense variations Ile436Ser and Thr473Ile were not identified in the patients enrolled in this study. Gene expression analyses showed lower LEDGF/p75 mRNA levels in peripheral blood mononuclear cells obtained from HIV-1 infected individuals. However, these levels were not influenced by any of the SNPs investigated. In spite of the limited number of LTNPs, these data suggest that the PSIP1 gene could be associated with the outcome of HIV-1 infection. Further analyses of this gene may guide the identification of causative variants to help predict disease course.     

Changes in Uric Acid Levels following Bariatric Surgery Are Not Associated with SLC2A9 Variants in the Swedish Obese Subjects Study

Context and Objective

Obesity and SLC2A9 genotype are strong determinants of uric acid levels. However, data on SLC2A9 variants and weight loss induced changes in uric acid levels are missing. We examined whether the changes in uric acid levels two- and ten-years after weight loss induced by bariatric surgery were associated with SLC2A9 single nucleotide polymorphisms (SNPs) in the Swedish Obese Subjects study.

Methods

SNPs (N = 14) identified by genome-wide association studies and exonic SNPs in the SLC2A9 gene locus were genotyped. Cross-sectional associations were tested before (N = 1806), two (N = 1664) and ten years (N = 1201) after bariatric surgery. Changes in uric acid were compared between baseline and Year 2 (N = 1660) and years 2 and 10 (N = 1172). A multiple testing corrected threshold of P = 0.007 was used for statistical significance.

Results

Overall, 11 of the 14 tested SLC2A9 SNPs were significantly associated with cross-sectional uric acid levels at all three time points, with rs13113918 showing the strongest association at each time point (R² = 3.7−5.2%, 3.9×10⁻²²≤p≤7.7×10⁻¹¹). One SNP (rs737267) showed a significant association (R² = 0.60%, P = 0.002) with change in uric acid levels from baseline to Year 2, as common allele homozygotes (C/C, N = 957) showed a larger decrease in uric acid (−61.4 µmol/L) compared to minor allele carriers (A/X: −51.7 µmol/L, N = 702). No SNPs were associated with changes in uric acid from years 2 to 10.

Conclusions

SNPs in the SLC2A9 locus contribute significantly to uric acid levels in obese individuals, and the associations persist even after considerable weight loss due to bariatric surgery. However, we found little evidence for an interaction between genotype and weight change on the response of uric acid to bariatric surgery over ten years. Thus, the fluctuations in uric acid levels among the surgery group appear to be driven by the weight losses and gains, independent of SLC2A9 genotypes.     

Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium

Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3×10⁻⁶⁴) and lower levels of eicosapentaenoic acid (EPA, p = 5×10⁻⁵⁸) and docosapentaenoic acid (DPA, p = 4×10⁻¹⁵⁴). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2×10⁻¹²) and DPA (p = 1×10⁻⁴³) and lower docosahexaenoic acid (DHA, p = 1×10⁻¹⁵). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1×10⁻⁸). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.