Genome-wide associated loci influencing interleukin (IL)-10, IL-1Ra, and IL-6 levels in African Americans

Interleukins (ILs) are key mediators of the immune response and inflammatory process. Plasma levels of IL-10, IL-1Ra, and IL-6 are associated with metabolic conditions, show large inter-individual variations, and are under strong genetic control. Therefore, elucidation of the genetic variants that influence levels of these ILs provides useful insights into mechanisms of immune response and pathogenesis of diseases. We conducted a genome-wide association study (GWAS) of IL-10, IL-1Ra, and IL-6 levels in 707 non-diabetic African Americans using 5,396,780 imputed and directly genotyped single nucleotide polymorphisms (SNPs) with adjustment for gender, age, and body mass index. IL-10 levels showed genome-wide significant associations (p < 5 × 10(-8)) with eight SNPs, the most significant of which was rs5743185 in the PMS1 gene (p = 2.30 × 10(-10)). We tested replication of SNPs that showed genome-wide significance in 425 non-diabetic individuals from West Africa, and successfully replicated rs17365948 in the YWHAZ gene (p = 0.02). IL-1Ra levels showed suggestive associations with two SNPs in the ASB3 gene (p = 2.55 × 10(-7)), ten SNPs in the IL-1 gene family (IL1F5, IL1F8, IL1F10, and IL1Ra, p = 1.04 × 10(-6) to 1.75 × 10(-6)), and 23 SNPs near the IL1A gene (p = 1.22 × 10(-6) to 1.63 × 10(-6)). We also successfully replicated rs4251961 (p = 0.009); this SNP was reported to be associated with IL-1Ra levels in a candidate gene study of Europeans. IL-6 levels showed genome-wide significant association with one SNP (RP11-314E23.1; chr6:133397598; p = 8.63 × 10(-9)). To our knowledge, this is the first GWAS on IL-10, IL-1Ra, and IL-6 levels. Follow-up of these findings may provide valuable insight into the pathobiology of IL actions and dysregulations in inflammation and human diseases.     

Adiponectin concentrations: a genome-wide association study

Adiponectin is associated with obesity and insulin resistance. To date, there has been no genome-wide association study (GWAS) of adiponectin levels in Asians. Here we present a GWAS of a cohort of Korean volunteers. A total of 4,001 subjects were genotyped by using a genome-wide marker panel in a two-stage design (979 subjects initially and 3,022 in a second stage). Another 2,304 subjects were used for follow-up replication studies with selected markers. In the discovery phase, the top SNP associated with mean log adiponectin was rs3865188 in CDH13 on chromosome 16 (p = 1.69 × 10(-15) in the initial sample, p = 6.58 × 10(-39) in the second genome-wide sample, and p = 2.12 × 10(-32) in the replication sample). The meta-analysis p value for rs3865188 in all 6,305 individuals was 2.82 × 10(-83). The association of rs3865188 with high-molecular-weight adiponectin (p = 7.36 × 10(-58)) was even stronger in the third sample. A reporter assay that evaluated the effects of a CDH13 promoter SNP in complete linkage disequilibrium with rs3865188 revealed that the major allele increased expression 2.2-fold. This study clearly shows that genetic variants in CDH13 influence adiponectin levels in Korean adults.     

Complement receptor 1 gene variants are associated with erythrocyte sedimentation rate

The erythrocyte sedimentation rate (ESR), a commonly performed test of the acute phase response, is the rate at which erythrocytes sediment in vitro in 1 hr. The molecular basis of erythrocyte sedimentation is unknown. To identify genetic variants associated with ESR, we carried out a genome-wide association study of 7607 patients in the Electronic Medical Records and Genomics (eMERGE) network. The discovery cohort consisted of 1979 individuals from the Mayo Clinic, and the replication cohort consisted of 5628 individuals from the remaining four eMERGE sites. A nonsynonymous SNP, rs6691117 (Val→IIe), in the complement receptor 1 gene (CR1) was associated with ESR (discovery cohort p = 7 × 10(-12), replication cohort p = 3 × 10(-14), combined cohort p = 9 × 10(-24)). We imputed 61 SNPs in CR1, and a "possibly damaging" SNP (rs2274567, His→Arg) in linkage disequilibrium (r(2) = 0.74) with rs6691117 was also associated with ESR (discovery p = 5 × 10(-11), replication p = 7 × 10(-17), and combined cohort p = 2 × 10(-25)). The two nonsynonymous SNPs in CR1 are near the C3b/C4b binding site, suggesting a possible mechanism by which the variants may influence ESR. In conclusion, genetic variation in CR1, which encodes a protein that clears complement-tagged inflammatory particles from the circulation, influences interindividual variation in ESR, highlighting an association between the innate immunity pathway and erythrocyte interactions.     

Large-scale genome-wide association study of Asian population reveals genetic factors in FRMD4A and other loci influencing smoking initiation and nicotine dependence

Diseases related to smoking are the second leading cause of death in the world. Cigarette smoking is a risk factor for several diseases such as cancer and cardiovascular and respiratory disorders. Despite increasing evidence of genetic determination, the susceptibility genes and loci underlying various aspects of smoking behavior are largely unknown. Moreover, almost all reported genome-wide association studies (GWASs) have been performed on samples of European origin, limiting the applicability of the results to other ethnic populations. In this first GWAS on smoking behavior in an Asian population, after analyzing 8,842 DNA samples from the Korea Association Resource project with 352,228 single nucleotide polymorphisms (SNPs) genotyped for each sample, we identified 8 SNPs significantly associated with smoking initiation (SI) and 4 with nicotine dependence (ND). Because of the current unavailability of an independent Asian smoking sample, we replicated the discoveries in independent samples of European-American and African-American origin. Of the 12 SNPs examined in the replicated samples, we identified two SNPs, in the regulator of G-protein signaling 17 gene (rs7747583, p value(meta)?=?6.40?×?10(-6); rs2349433, p value(meta)?=?5.57?×?10(-6)), associated with SI. Also, we found two SNPs significantly associated with ND; one in the FERM domain containing 4A (rs4424567, p value(meta)?=?2.30?×?10(-6)) and the other at 7q31.1 (rs848353, p value(meta)?=?9.16?×?10(-8)). These SNPs represent novel targets for examination of smoking behavior and warrant further investigation using independent samples.     

Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci

To identify genetic factors contributing to type 2 diabetes (T2D), we performed large-scale meta-analyses by using a custom ～50,000 SNP genotyping array (the ITMAT-Broad-CARe array) with ～2000 candidate genes in 39 multiethnic population-based studies, case-control studies, and clinical trials totaling 17,418 cases and 70,298 controls. First, meta-analysis of 25 studies comprising 14,073 cases and 57,489 controls of European descent confirmed eight established T2D loci at genome-wide significance. In silico follow-up analysis of putative association signals found in independent genome-wide association studies (including 8,130 cases and 38,987 controls) performed by the DIAGRAM consortium identified a T2D locus at genome-wide significance (GATAD2A/CILP2/PBX4; p = 5.7 × 10(-9)) and two loci exceeding study-wide significance (SREBF1, and TH/INS; p < 2.4 × 10(-6)). Second, meta-analyses of 1,986 cases and 7,695 controls from eight African-American studies identified study-wide-significant (p = 2.4 × 10(-7)) variants in HMGA2 and replicated variants in TCF7L2 (p = 5.1 × 10(-15)). Third, conditional analysis revealed multiple known and novel independent signals within five T2D-associated genes in samples of European ancestry and within HMGA2 in African-American samples. Fourth, a multiethnic meta-analysis of all 39 studies identified T2D-associated variants in BCL2 (p = 2.1 × 10(-8)). Finally, a composite genetic score of SNPs from new and established T2D signals was significantly associated with increased risk of diabetes in African-American, Hispanic, and Asian populations. In summary, large-scale meta-analysis involving a dense gene-centric approach has uncovered additional loci and variants that contribute to T2D risk and suggests substantial overlap of T2D association signals across multiple ethnic groups.     

Genome-wide association study of treatment refractory schizophrenia in Han Chinese

We report the first genome-wide association study of a joint analysis using 795 Han Chinese individuals with treatment-refractory schizophrenia (TRS) and 806 controls. Three loci showed suggestive significant association with TRS were identified. These loci include: rs10218843 (P = 3.04 × 10(-7)) and rs11265461 (P = 1.94 × 10(-7)) are adjacent to signaling lymphocytic activation molecule family member 1 (SLAMF1); rs4699030 (P = 1.94 × 10(-6)) and rs230529 (P = 1.74 × 10(-7)) are located in the gene nuclear factor of kappa light polypeptide gene enhancer in B-cells 1 (NFKB1); and rs13049286 (P = 3.05 × 10(-5)) and rs3827219 (P = 1.66 × 10(-5)) fall in receptor-interacting serine/threonine-protein kinase 4 (RIPK4). One isolated single nucleotide polymorphism (SNP), rs739617 (P = 3.87 × 10(-5)) was also identified to be associated with TRS. The -94delATTG allele (rs28362691) located in the promoter region of NFKB1 was identified by resequencing and was found to associate with TRS (P = 4.85 × 10(-6)). The promoter assay demonstrated that the -94delATTG allele had a significant lower promoter activity than the -94insATTG allele in the SH-SY5Y cells. This study suggests that rs28362691 in NFKB1 might be involved in the development of TRS.     

Genetic associations for activated partial thromboplastin time and prothrombin time, their gene expression profiles, and risk of coronary artery disease

Activated partial thromboplastin time (aPTT) and prothrombin time (PT) are clinical tests commonly used to screen for coagulation-factor deficiencies. One genome-wide association study (GWAS) has been reported previously for aPTT, but no GWAS has been reported for PT. We conducted a GWAS and meta-analysis to identify genetic loci for aPTT and PT. The GWAS for aPTT was conducted in 9,240 individuals of European ancestry from the Atherosclerosis Risk in Communities (ARIC) study, and the GWAS for PT was conducted in 2,583 participants from the Genetic Study of Three Population Microisolates in South Tyrol (MICROS) and the Lothian Birth Cohorts (LBC) of 1921 and 1936. Replication was assessed in 1,041 to 3,467 individuals. For aPTT, previously reported associations with KNG1, HRG, F11, F12, and ABO were confirmed. A second independent association in ABO was identified and replicated (rs8176704, p = 4.26 × 10(-24)). Pooling the ARIC and replication data yielded two additional loci in F5 (rs6028, p = 3.22 × 10(-9)) and AGBL1 (rs2469184, p = 3.61 × 10(-8)). For PT, significant associations were identified and confirmed in F7 (rs561241, p = 3.71 × 10(-56)) and PROCR/EDEM2 (rs2295888, p = 5.25 × 10(-13)). Assessment of existing gene expression and coronary artery disease (CAD) databases identified associations of five of the GWAS loci with altered gene expression and two with CAD. In summary, eight genetic loci that account for ～29% of the variance in aPTT and two loci that account for ～14% of the variance in PT were detected and supported by functional data.     

A genome-wide association analysis implicates SOX6 as a candidate gene for wrist bone mass

Osteoporosis is a highly heritable common bone disease leading to fractures that severely impair the life quality of patients. Wrist fractures caused by osteoporosis are largely due to the scarcity of wrist bone mass. Here we report the results of a genome-wide association study (GWAS) of wrist bone mineral density (BMD). We examined ～500000 SNP markers in 1000 unrelated homogeneous Caucasian subjects and found a novel allelic association with wrist BMD at rs11023787 in the SOX6 (SRY (sex determining region Y)-box 6) gene (P=9.00×10(-5)). Subjects carrying the C allele of rs11023787 in SOX6 had significantly higher mean wrist BMD values than those with the T allele (0.485:0.462 g cm(-2) for C allele vs. T allele carriers). For validation, we performed SOX6 association for BMD in an independent Chinese sample and found that SNP rs11023787 was significantly associated with wrist BMD in the Chinese sample (P=6.41×10(-3)). Meta-analyses of the GWAS scan and the replication studies yielded P-values of 5.20×10(-6) for rs11023787. Results of this study, together with the functional relevance of SOX6 in cartilage formation, support the SOX6 gene as an important gene for BMD variation.     

Meta-analysis of new genome-wide association studies of colorectal cancer risk

Colorectal cancer is the second leading cause of cancer death in developed countries. Genome-wide association studies (GWAS) have successfully identified novel susceptibility loci for colorectal cancer. To follow up on these findings, and try to identify novel colorectal cancer susceptibility loci, we present results for GWAS of colorectal cancer (2,906 cases, 3,416 controls) that have not previously published main associations. Specifically, we calculated odds ratios and 95% confidence intervals using log-additive models for each study. In order to improve our power to detect novel colorectal cancer susceptibility loci, we performed a meta-analysis combining the results across studies. We selected the most statistically significant single nucleotide polymorphisms (SNPs) for replication using ten independent studies (8,161 cases and 9,101 controls). We again used a meta-analysis to summarize results for the replication studies alone, and for a combined analysis of GWAS and replication studies. We measured ten SNPs previously identified in colorectal cancer susceptibility loci and found eight to be associated with colorectal cancer (p value range 0.02 to 1.8?×?10(-8)). When we excluded studies that have previously published on these SNPs, five SNPs remained significant at p?相似文献   

A genome-wide association study of total bilirubin and cholelithiasis risk in sickle cell anemia

Serum bilirubin levels have been associated with polymorphisms in the UGT1A1 promoter in normal populations and in patients with hemolytic anemias, including sickle cell anemia. When hemolysis occurs circulating heme increases, leading to elevated bilirubin levels and an increased incidence of cholelithiasis. We performed the first genome-wide association study (GWAS) of bilirubin levels and cholelithiasis risk in a discovery cohort of 1,117 sickle cell anemia patients. We found 15 single nucleotide polymorphisms (SNPs) associated with total bilirubin levels at the genome-wide significance level (p value <5 × 10(-8)). SNPs in UGT1A1, UGT1A3, UGT1A6, UGT1A8 and UGT1A10, different isoforms within the UGT1A locus, were identified (most significant rs887829, p = 9.08 × 10(-25)). All of these associations were validated in 4 independent sets of sickle cell anemia patients. We tested the association of the 15 SNPs with cholelithiasis in the discovery cohort and found a significant association (most significant p value 1.15 × 10(-4)). These results confirm that the UGT1A region is the major regulator of bilirubin metabolism in African Americans with sickle cell anemia, similar to what is observed in other ethnicities.     

Genetic variation in antioxidant enzymes and lung function

Not all cigarette smokers develop chronic obstructive pulmonary disease, and discovering susceptibility factors is an important research priority. The oxidative burden of smoking may overwhelm antioxidant defenses, and vulnerabilities may exist as a result of sequence variants in genes encoding antioxidant enzymes. This study explored the association between genetic variation in a network of antioxidant enzymes and lung phenotypes. Linear models evaluated single-locus marker associations in 2387 European American and African American participants in the Health, Aging, and Body Composition Study. After corrections were made for multiple comparisons, 15 statistically significant associations were identified, all of which were for SNP by smoking interactions. The most statistically significant findings were for genes encoding members of the isocitrate dehydrogenase gene family (IDH3A, IDH3B, IDH2). For rs6107100 (IDH3B) the variant genotype was associated with a difference of 6% in the FEV(1)/FVC ratio in African American current smokers, but the SNP had little or no association with FEV(1)/FVC in former and never smokers (nominal p(interaction)=5×10(-6)). A variant of the peroxiredoxin gene (rs9787810, PRDX5) was associated with lower percentage predicted FEV(1) and a lower ratio in European American current smokers, with little or no association in other smoking groups (nominal p(interaction)=0.0001 and 0.0003, respectively). The studied genes have not been reported in previous candidate gene association studies, and thus the findings suggest novel mechanisms and targets for future research and provide evidence for a contribution of sequence variation in genes encoding antioxidant enzymes to susceptibility in smokers.     

Four novel Loci (19q13, 6q24, 12q24, and 5q14) influence the microcirculation in vivo

There is increasing evidence that the microcirculation plays an important role in the pathogenesis of cardiovascular diseases. Changes in retinal vascular caliber reflect early microvascular disease and predict incident cardiovascular events. We performed a genome-wide association study to identify genetic variants associated with retinal vascular caliber. We analyzed data from four population-based discovery cohorts with 15,358 unrelated Caucasian individuals, who are members of the Cohort for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and replicated findings in four independent Caucasian cohorts (n = 6,652). All participants had retinal photography and retinal arteriolar and venular caliber measured from computer software. In the discovery cohorts, 179 single nucleotide polymorphisms (SNP) spread across five loci were significantly associated (p<5.0×10(-8)) with retinal venular caliber, but none showed association with arteriolar caliber. Collectively, these five loci explain 1.0%-3.2% of the variation in retinal venular caliber. Four out of these five loci were confirmed in independent replication samples. In the combined analyses, the top SNPs at each locus were: rs2287921 (19q13; p = 1.61×10(-25), within the RASIP1 locus), rs225717 (6q24; p?=?1.25×10(-16), adjacent to the VTA1 and NMBR loci), rs10774625 (12q24; p = 2.15×10(-13), in the region of ATXN2,SH2B3 and PTPN11 loci), and rs17421627 (5q14; p?=?7.32×10(-16), adjacent to the MEF2C locus). In two independent samples, locus 12q24 was also associated with coronary heart disease and hypertension. Our population-based genome-wide association study demonstrates four novel loci associated with retinal venular caliber, an endophenotype of the microcirculation associated with clinical cardiovascular disease. These data provide further insights into the contribution and biological mechanisms of microcirculatory changes that underlie cardiovascular disease.     

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ～2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)相似文献   

Association of common variants in TNFRSF13B, TNFSF13, and ANXA3 with serum levels of non-albumin protein and immunoglobulin isotypes in Japanese

We performed a genome-wide association study (GWAS) on levels of serum total protein (TP), albumin (ALB), and non-albumin protein (NAP). We analyzed SNPs on autosomal chromosomes using data from 9,103 Japanese individuals, followed by a replication study of 1,600 additional individuals. We confirmed the previously- reported association of GCKR on chromosome 2p23.3 with serum ALB (rs1260326, P(meta) = 3.1 × 10(-9)), and additionally identified the significant genome-wide association of rs4985726 in TNFRSF13B on 17p11.2 with both TP and NAP (P(meta) = 1.2 × 10(-14) and 7.1 × 10(-24), respectively). For NAP, rs3803800 and rs11552708 in TNFSF13 on 17p13.1 (P(meta) = 7.2 × 10(-15) and 7.5 × 10(-10), respectively) as well as rs10007186 on 4q21.2 near ANXA3 (P(meta) = 1.3 × 10(-9)) also indicated significant associations. Interestingly, TNFRSF13B and TNFSF13 encode a tumor necrosis factor (TNF) receptor and its ligand, which together constitute an important receptor-ligand axis for B-cell homeostasis and immunoglobulin production. Furthermore, three SNPs, rs4985726, rs3803800, and rs11552708 in TNFRSF13B and TNFSF13, were indicated to be associated with serum levels of IgG (P<2.3 × 10(-3)) and IgM (P<0.018), while rs3803800 and rs11552708 were associated with IgA (P<0.013). Rs10007186 in 4q21.2 was associated with serum levels of IgA (P = 0.036), IgM (P = 0.019), and IgE (P = 4.9 × 10(-4)). Our results should add interesting knowledge about the regulation of major serum components.     

Gene promoter methylation assayed in exhaled breath,with differences in smokers and lung cancer patients

Background

There is a need for new, noninvasive risk assessment tools for use in lung cancer population screening and prevention programs.

Methods

To investigate the technical feasibility of determining DNA methylation in exhaled breath condensate, we applied our previously-developed method for tag-adapted bisulfite genomic DNA sequencing (tBGS) for mapping of DNA methylation, and adapted it to exhaled breath condensate (EBC) from lung cancer cases and non-cancer controls. Promoter methylation patterns were analyzed in DAPK, RASSF1A and PAX5β promoters in EBC samples from 54 individuals, comprised of 37 controls [current- (n = 19), former- (n = 10), and never-smokers (n = 8)] and 17 lung cancer cases [current- (n = 5), former- (n = 11), and never-smokers (n = 1)].

Results

We found: (1) Wide inter-individual variability in methylation density and spatial distribution for DAPK, PAX5β and RASSF1A. (2) Methylation patterns from paired exhaled breath condensate and mouth rinse specimens were completely divergent. (3) For smoking status, the methylation density of RASSF1A was statistically different (p = 0.0285); pair-wise comparisons showed that the former smokers had higher methylation density versus never smokers and current smokers (p = 0.019 and p = 0.031). For DAPK and PAX5β, there was no such significant smoking-related difference. Underlying lung disease did not impact on methylation density for this geneset. (4) In case-control comparisons, CpG at -63 of DAPK promoter and +52 of PAX5β promoter were significantly associated with lung cancer status (p = 0.0042 and 0.0093, respectively). After adjusting for multiple testing, both loci were of borderline significance (p_adj = 0.054 and 0.031). (5) The DAPK gene had a regional methylation pattern with two blocks (1)~-215~-113 and (2) -84 ~+26); while similar in block 1, there was a significant case-control difference in methylation density in block 2 (p = 0.045); (6)Tumor stage and histology did not impact on the methylation density among the cases. (7) The results of qMSP applied to EBC correlated with the corresponding tBGS sequencing map loci.

Conclusion

Our results show that DNA methylation in exhaled breath condensate is detectable and is likely of lung origin. Suggestive correlations with smoking and lung cancer case-control status depend on individual gene and CpG site examined.     

Genome-wide copy number analysis uncovers a new HSCR gene: NRG3

Hirschsprung disease (HSCR) is a congenital disorder characterized by aganglionosis of the distal intestine. To assess the contribution of copy number variants (CNVs) to HSCR, we analysed the data generated from our previous genome-wide association study on HSCR patients, whereby we identified NRG1 as a new HSCR susceptibility locus. Analysis of 129 Chinese patients and 331 ethnically matched controls showed that HSCR patients have a greater burden of rare CNVs (p = 1.50 × 10(-5)), particularly for those encompassing genes (p = 5.00 × 10(-6)). Our study identified 246 rare-genic CNVs exclusive to patients. Among those, we detected a NRG3 deletion (p = 1.64 × 10(-3)). Subsequent follow-up (96 additional patients and 220 controls) on NRG3 revealed 9 deletions (combined p = 3.36 × 10(-5)) and 2 de novo duplications among patients and two deletions among controls. Importantly, NRG3 is a paralog of NRG1. Stratification of patients by presence/absence of HSCR-associated syndromes showed that while syndromic-HSCR patients carried significantly longer CNVs than the non-syndromic or controls (p = 1.50 × 10(-5)), non-syndromic patients were enriched in CNV number when compared to controls (p = 4.00 × 10(-6)) or the syndromic counterpart. Our results suggest a role for NRG3 in HSCR etiology and provide insights into the relative contribution of structural variants in both syndromic and non-syndromic HSCR. This would be the first genome-wide catalog of copy number variants identified in HSCR.