Tuberculosis infection and lung adenocarcinoma: Mendelian randomization and pathway analysis of genome-wide association study data from never-smoking Asian women

We investigated whether genetic susceptibility to tuberculosis (TB) influences lung adenocarcinoma development among never-smokers using TB genome-wide association study (GWAS) results within the Female Lung Cancer Consortium in Asia. Pathway analysis with the adaptive rank truncated product method was used to assess the association between a TB-related gene-set and lung adenocarcinoma using GWAS data from 5512 lung adenocarcinoma cases and 6277 controls. The gene-set consisted of 31 genes containing known/suggestive associations with genetic variants from previous TB-GWAS. Subsequently, we followed-up with Mendelian Randomization to evaluate the association between TB and lung adenocarcinoma using three genome-wide significant variants from previous TB-GWAS in East Asians. The TB-related gene-set was associated with lung adenocarcinoma (p = 0.016). Additionally, the Mendelian Randomization showed an association between TB and lung adenocarcinoma (OR = 1.31, 95% CI: 1.03, 1.66, p = 0.027). Our findings support TB as a causal risk factor for lung cancer development among never-smoking Asian women.     

An analysis of single nucleotide polymorphisms of 125 DNA repair genes in the Texas genome-wide association study of lung cancer with a replication for the XRCC4 SNPs

DNA repair genes are important for maintaining genomic stability and limiting carcinogenesis. We analyzed all single nucleotide polymorphisms (SNPs) of 125 DNA repair genes covered by the Illumina HumanHap300 (v1.1) BeadChips in a previously conducted genome-wide association study (GWAS) of 1154 lung cancer cases and 1137 controls and replicated the top-hits of XRCC4 SNPs in an independent set of 597 cases and 611 controls in Texas populations. We found that six of 20 XRCC4 SNPs were associated with a decreased risk of lung cancer with a P-value of 0.01 or lower in the discovery dataset, of which the most significant SNP was rs10040363 (P for allelic test=4.89 x 10??). Moreover, the data in this region allowed us to impute a potentially functional SNP rs2075685 (imputed P for allelic test=1.3 x 10?3). A luciferase reporter assay demonstrated that the rs2075685G>T change in the XRCC4 promoter increased expression of the gene. In the replication study of rs10040363, rs1478486, rs9293329, and rs2075685, however, only rs10040363 achieved a borderline association with a decreased risk of lung cancer in a dominant model (adjusted OR=0.80, 95% CI=0.62-1.03 and P=0.079). In the final combined analysis of both the Texas GWAS discovery and replication datasets, the strength of the association was increased for rs10040363 (adjusted OR=0.77, 95% CI=0.66-0.89, P(dominant)=5 x 10?? and P for trend=5 x 10??) and rs1478486 (adjusted OR=0.82, 95% CI=0.71-0.94, P(dominant)=6 x 10?3 and P for trend=3.5 x 10?3). Finally, we conducted a meta-analysis of these XRCC4 SNPs with available data from published GWA studies of lung cancer with a total of 12,312 cases and 47,921 controls, in which none of these XRCC4 SNPs was associated with lung cancer risk. It appeared that rs2075685, although associated with increased expression of a reporter gene and lung cancer risk in the Texas populations, did not have an effect on lung cancer risk in other populations. This study underscores the importance of replication using published data in larger populations.     

Genome-wide association study on bipolar disorder in the Bulgarian population

Bipolar disorder is a severe psychiatric disorder influenced by environmental and genetic factors. Genetic studies have implicated many variants in the disease's etiology but only few have been successfully replicated. We conducted a genome-wide association study (GWAS) on bipolar disorder in the Bulgarian population followed by a replication study of the top 100 single nucleotide polymorphisms (SNPs) showing the smallest P values. The GWAS was performed on 188 bipolar disorder patients and 376 control subjects genotyped on the Illumina 550 platform. The replication study was conducted on 122 patients and 328 controls. Although our study did not show any association P value that achieved genome-wide significance, and none of the top 100 SNPs reached the Bonferroni-corrected P value in the replication study, the plausible involvement of some variants cannot be entirely discarded. Three polymorphisms, rs8099939 [P = 2.12 × 10(-6), odds ratio (OR) = 1.95, 95% confidence interval (CI) = 1.43-2.67] in GRIK5, rs6122972 (P = 3.11 × 10(-6), OR = 2.02, 95% CI = 1.46-2.80) in PARD6B and rs2289700 (P = 9.14 × 10(-6), OR = 2.13, 95% CI = 1.53-2.95) in CTSH remained associated at a similar level after Mantel-Haenszel test for combining the results from the genome-wide and replication studies. A modest association was also detected for SNP rs1012053 (GWAS P = 4.50 × 10(-2)) in DGKH, which has already been reported as the most significant variant in a previous genome-wide scan on bipolar disorder. However, further studies using larger datasets are needed to identify variants with smaller effects that contribute to the risk of bipolar disorder.     

Effects of cis and trans regulatory variations on the expression divergence of heat shock response genes between yeast strains

Recently, two genome-wide association studies in Asia identified gene polymorphisms known as rs4488809, rs9816619 in TP63 and rs2131877, rs952481 in C3orf21. It has been proposed that these polymorphisms are susceptibility loci for non-small cell lung cancer (NSCLC) development among Japanese and Korean populations. We ask whether susceptibility to NSCLC is limited to the Chinese population or whether the environment also affects genetic polymorphisms. We conducted a matched case-control study to explore this question. Results show that polymorphism of TP63 was not associated with NSCLC development, whereas variant genotypes of C3orf21 were nominally associated with a reduced risk of lung adenocarcinoma (OR=0.619, 95% CI=0.390-0.976). These results strongly suggest that environmental agents interact with human genetic polymorphism independent of ethnic background. In addition, the C3orf21 gene may be a potential susceptibility marker for lung adenocarcinoma independent of ethnic background and environmental agents.     

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.     

Identification of a systemic lupus erythematosus susceptibility locus at 11p13 between PDHX and CD44 in a multiethnic study

Systemic lupus erythematosus (SLE) is considered to be the prototypic autoimmune disease, with a complex genetic architecture influenced by environmental factors. We sought to replicate a putative association at 11p13 not yet exceeding genome-wide significance (p < 5 × 10(-8)) identified in a genome-wide association study (GWAS). Our GWA scan identified two intergenic SNPs located between PDHX and CD44 showing suggestive evidence of association with SLE in cases of European descent (rs2732552, p = 0.004, odds ratio [OR] = 0.78; rs387619, p = 0.003, OR = 0.78). The replication cohort consisted of >15,000 subjects, including 3562 SLE cases and 3491 controls of European ancestry, 1527 cases and 1811 controls of African American (AA) descent, and 1265 cases and 1260 controls of Asian origin. We observed robust association at both rs2732552 (p = 9.03 × 10(-8), OR = 0.83) and rs387619 (p = 7.7 × 10(-7), OR = 0.83) in the European samples with p(meta) = 1.82 × 10(-9) for rs2732552. The AA and Asian SLE cases also demonstrated association at rs2732552 (p = 5 × 10(-3), OR = 0.81 and p = 4.3 × 10(-4), OR = 0.80, respectively). A meta-analysis of rs2732552 for all racial and ethnic groups studied produced p(meta) = 2.36 × 10(-13). This locus contains multiple regulatory sites that could potentially affect expression and functions of CD44, a cell-surface glycoprotein influencing immunologic, inflammatory, and oncologic phenotypes, or PDHX, a subunit of the pyruvate dehydrogenase complex.     

Use of a multiethnic approach to identify rheumatoid- arthritis-susceptibility loci, 1p36 and 17q12

We have previously shown that rheumatoid arthritis (RA) risk alleles overlap between different ethnic groups. Here, we utilize a multiethnic approach to show that we can effectively discover RA risk alleles. Thirteen putatively associated SNPs that had not yet exceeded genome-wide significance (p < 5 × 10(-8)) in our previous RA genome-wide association study (GWAS) were analyzed in independent sample sets consisting of 4,366 cases and 17,765 controls of European, African American, and East Asian ancestry. Additionally, we conducted an overall association test across all 65,833 samples (a GWAS meta-analysis plus the replication samples). Of the 13 SNPs investigated, four were significantly below the study-wide Bonferroni corrected p value threshold (p < 0.0038) in the replication samples. Two SNPs (rs3890745 at the 1p36 locus [p = 2.3 × 10(-12)] and rs2872507 at the 17q12 locus [p = 1.7 × 10(-9)]) surpassed genome-wide significance in all 16,659 RA cases and 49,174 controls combined. We used available GWAS data to fine map these two loci in Europeans and East Asians, and we found that the same allele conferred risk in both ethnic groups. A series of bioinformatic analyses identified TNFRSF14-MMEL1 at the 1p36 locus and IKZF3-ORMDL3-GSDMB at the 17q12 locus as the genes most likely associated with RA. These findings demonstrate empirically that a multiethnic approach is an effective strategy for discovering RA risk loci, and they suggest that combining GWASs across ethnic groups represents an efficient strategy for gaining statistical power.     

Quantitative Assessment of the Influence of TP63 Gene Polymorphisms and Lung Cancer Risk: Evidence Based on 93,751 Subjects

Background

Several genome-wide association studies on lung cancer (LC) have reported similar findings of a new susceptibility locus, 3q28. After that, a number of studies reported that the rs10937405, and rs4488809 polymorphism in chromosome 3q28 has been implicated in LC risk. However, the studies have yielded contradictory results.

Methods

PubMed, ISI web of science, EMBASE and the Chinese National Knowledge Infrastructure databases were systematically searched to identify relevant studies. Data were abstracted independently by two reviewers. A meta-analysis was performed to examine the association between rs10937405, rs4488809 polymorphism at 3q28 and susceptibility to LC. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated. Heterogeneity and publication bias were also tested.

Results

A total of 9 studies including 35,961 LC cases and 57,790 controls were involved in this meta-analysis. An overall random-effects per-allele OR of1.19 (95% CI: 1.14–1.25; P<10⁻⁵) and 1.19 (95% CI: 1.13–1.25; P<10⁻⁵) was found for the rs10937405 and rs4488809 polymorphism respectively. Similar results were also observed using dominant or recessive genetic model. After stratified by ethnicity, significant associations were found among East Asians (per-allele OR = 1.22, 95% CI: 1.17–1.27; P<10⁻⁵); whereas no significant associations were found among Caucasians for rs10937405. In the sub-group analysis by sample size, significantly increased risks were found for these polymorphisms in all genetic models. When analyzed according to histological type, the effects of rs10937405, and rs4488809 at 3q28 on the risk of lung cancer were significant mostly for lung adenocarcinoma.

Conclusions

Our findings demonstrated that rs10937405-G allele and rs4488809-G allele might be risk-conferring factors for the development of lung cancer, especially for East Asian populations.     

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.     

GWAS of follicular lymphoma reveals allelic heterogeneity at 6p21.32 and suggests shared genetic susceptibility with diffuse large B-cell lymphoma

Non-Hodgkin lymphoma (NHL) represents a diverse group of hematological malignancies, of which follicular lymphoma (FL) is a prevalent subtype. A previous genome-wide association study has established a marker, rs10484561 in the human leukocyte antigen (HLA) class II region on 6p21.32 associated with increased FL risk. Here, in a three-stage genome-wide association study, starting with a genome-wide scan of 379 FL cases and 791 controls followed by validation in 1,049 cases and 5,790 controls, we identified a second independent FL-associated locus on 6p21.32, rs2647012 (OR(combined) = 0.64, P(combined) = 2 × 10(-21)) located 962 bp away from rs10484561 (r(2)<0.1 in controls). After mutual adjustment, the associations at the two SNPs remained genome-wide significant (rs2647012:OR(adjusted) = 0.70, P(adjusted) = 4 × 10(-12); rs10484561:OR(adjusted) = 1.64, P(adjusted) = 5 × 10(-15)). Haplotype and coalescence analyses indicated that rs2647012 arose on an evolutionarily distinct haplotype from that of rs10484561 and tags a novel allele with an opposite (protective) effect on FL risk. Moreover, in a follow-up analysis of the top 6 FL-associated SNPs in 4,449 cases of other NHL subtypes, rs10484561 was associated with risk of diffuse large B-cell lymphoma (OR(combined) = 1.36, P(combined) = 1.4 × 10(-7)). Our results reveal the presence of allelic heterogeneity within the HLA class II region influencing FL susceptibility and indicate a possible shared genetic etiology with diffuse large B-cell lymphoma. These findings suggest that the HLA class II region plays a complex yet important role in NHL.     

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.     

Common genetic variants on 3q28 contribute to non-small cell lung cancer susceptibility: evidence from 10 case-control studies

The association between common variations (rs10937405, rs4488809) on 3q28 and lung cancer has been widely evaluated in various ethnic groups, since it was first identified through genome-wide association approach. However, the results have been inconclusive. To derive a more precise estimation of the relationship and the effect of factors that might modify the risk, we performed this meta-analysis. The random-effects model was applied, addressing heterogeneity and publication bias. A total of 10 articles involving 36,221 cases and 58,108 controls were included. Overall, the summary per-allele OR of 1.19 (95 % CI 1.14–1.25, P < 10^?5) and 1.17 (95 % CI 1.10–1.23, P < 10^?5) was found for the rs10937405 and rs4488809 polymorphisms, respectively. Significant results were also observed in heterozygous and homozygous when compared with wild genotype for these polymorphisms. Significant results were found in East Asians when stratified by ethnicity, whereas no significant associations were found among Caucasians. After stratifying by sample size, study design, control source and sex, significant associations were also obtained. In addition, our data indicate that these polymorphisms are involved in lung cancer susceptibility and confer its effect primarily in lung adenocarcinoma when stratified by histological subtype. Furthermore, significant associations were also detected both never-smokers and smokers for these polymorphisms. In conclusion, this meta-analysis demonstrated that rs10937405 and rs4488809 are a risk factor associated with increased non-small cell lung cancer susceptibility, particularly for East Asian populations.     

Genome-wide association analysis of incident coronary heart disease (CHD) in African Americans: a short report

African Americans have the highest rate of mortality due to coronary heart disease (CHD). Although multiple loci have been identified influencing CHD risk in European-Americans using a genome-wide association (GWAS) approach, no GWAS of incident CHD has been reported for African Americans. We performed a GWAS for incident CHD events collected during 19 years of follow-up in 2,905 African Americans from the Atherosclerosis Risk in Communities (ARIC) study. We identified a genome-wide significant SNP (rs1859023, MAF = 31%) located at 7q21 near the PFTK1 gene (HR = 0.57, 95% CI 0.46 to 0.69, p = 1.86×10(-08)), which replicated in an independent sample of over 8,000 African American women from the Women's Health Initiative (WHI) (HR = 0.81, 95% CI 0.70 to 0.93, p = 0.005). PFTK1 encodes a serine/threonine-protein kinase, PFTAIRE-1, that acts as a cyclin-dependent kinase regulating cell cycle progression and cell proliferation. This is the first finding of incident CHD locus identified by GWAS in African Americans.     

Genome-wide association study identifies a single major locus contributing to survival into old age; the APOE locus revisited

By studying the loci that contribute to human longevity, we aim to identify mechanisms that contribute to healthy aging. To identify such loci, we performed a genome-wide association study (GWAS) comparing 403 unrelated nonagenarians from long-living families included in the Leiden Longevity Study (LLS) and 1670 younger population controls. The strongest candidate SNPs from this GWAS have been analyzed in a meta-analysis of nonagenarian cases from the Rotterdam Study, Leiden 85-plus study, and Danish 1905 cohort. Only one of the 62 prioritized SNPs from the GWAS analysis (P<1×10(-4) ) showed genome-wide significance with survival into old age in the meta-analysis of 4149 nonagenarian cases and 7582 younger controls [OR=0.71 (95% CI 0.65-0.77), P=3.39 × 10(-17) ]. This SNP, rs2075650, is located in TOMM40 at chromosome 19q13.32 close to the apolipoprotein E (APOE) gene. Although there was only moderate linkage disequilibrium between rs2075650 and the ApoE ε4 defining SNP rs429358, we could not find an APOE-independent effect of rs2075650 on longevity, either in cross-sectional or in longitudinal analyses. As expected, rs429358 associated with metabolic phenotypes in the offspring of the nonagenarian cases from the LLS and their partners. In addition, we observed a novel association between this locus and serum levels of IGF-1 in women (P=0.005). In conclusion, the major locus determining familial longevity up to high age as detected by GWAS was marked by rs2075650, which tags the deleterious effects of the ApoE ε4 allele. No other major longevity locus was found.     

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.     

A Genome-wide Association Study of Lung Cancer Identifies a Region of Chromosome 5p15 Associated with Risk for Adenocarcinoma

Three genetic loci for lung cancer risk have been identified by genome-wide association studies (GWAS), but inherited susceptibility to specific histologic types of lung cancer is not well established. We conducted a GWAS of lung cancer and its major histologic types, genotyping 515,922 single-nucleotide polymorphisms (SNPs) in 5739 lung cancer cases and 5848 controls from one population-based case-control study and three cohort studies. Results were combined with summary data from ten additional studies, for a total of 13,300 cases and 19,666 controls of European descent. Four studies also provided histology data for replication, resulting in 3333 adenocarcinomas (AD), 2589 squamous cell carcinomas (SQ), and 1418 small cell carcinomas (SC). In analyses by histology, rs2736100 (TERT), on chromosome 5p15.33, was associated with risk of adenocarcinoma (odds ratio [OR] = 1.23, 95% confidence interval [CI] = 1.13–1.33, p = 3.02 × 10⁻⁷), but not with other histologic types (OR = 1.01, p = 0.84 and OR = 1.00, p = 0.93 for SQ and SC, respectively). This finding was confirmed in each replication study and overall meta-analysis (OR = 1.24, 95% CI = 1.17–1.31, p = 3.74 × 10⁻¹⁴ for AD; OR = 0.99, p = 0.69 and OR = 0.97, p = 0.48 for SQ and SC, respectively). Other previously reported association signals on 15q25 and 6p21 were also refined, but no additional loci reached genome-wide significance. In conclusion, a lung cancer GWAS identified a distinct hereditary contribution to adenocarcinoma.     

Genetic Variants at 6p21.1 and 7p15.3 Are Associated with Risk of Multiple Cancers in Han Chinese

Cancer susceptibility loci identified in reported genome-wide association studies (GWAS) are often tumor-specific; however, evidence of pleiotropy of some genes/loci has also been observed and biologically plausible. We hypothesized that there are important regions in the genome harboring genetic variants associated with risk of multiple types of cancer. In the current study, we attempted to map genetic variants that have consistent effects on risk of multiple cancers using our existing genome-wide scan data of lung cancer, noncardia gastric cancer, and esophageal squamous-cell carcinoma with overall 5,368 cases and 4,006 controls (GWAS stage), followed by a further evaluation in additional 9,001 cases with one of these cancer types and 11,436 controls (replication stage). Five variants satisfying the criteria of pleiotropy with p values from 1.10 × 10⁻⁸ to 8.96 × 10⁻⁶ for genome-wide scans of three cancer types were further evaluated in the replication stage. We found consistent associations of rs2494938 at 6p21.1 and rs2285947 at 7p15.3 with these three cancers in both GWAS and replication stages. In combined samples of GWAS and replication stages, the minor alleles of rs2494938 and rs2285947 were significantly associated with an increased risk of the cancers (odds ratio [OR] = 1.15, 95% confidence interval [CI], 1.10–1.19 and OR = 1.17, 95% CI, 1.12–1.21), with the p values being 1.20 × 10⁻¹² and 1.26 × 10⁻¹⁶, respectively, which are at a genome-wide significance level. Our findings highlight the potential importance of variants at 6p21.1 and 7p15.3 in the susceptibility to multiple cancers.     

A genome-wide association meta-analysis of circulating sex hormone-binding globulin reveals multiple Loci implicated in sex steroid hormone regulation

Sex hormone-binding globulin (SHBG) is a glycoprotein responsible for the transport and biologic availability of sex steroid hormones, primarily testosterone and estradiol. SHBG has been associated with chronic diseases including type 2 diabetes (T2D) and with hormone-sensitive cancers such as breast and prostate cancer. We performed a genome-wide association study (GWAS) meta-analysis of 21,791 individuals from 10 epidemiologic studies and validated these findings in 7,046 individuals in an additional six studies. We identified twelve genomic regions (SNPs) associated with circulating SHBG concentrations. Loci near the identified SNPs included SHBG (rs12150660, 17p13.1, p = 1.8 × 10(-106)), PRMT6 (rs17496332, 1p13.3, p = 1.4 × 10(-11)), GCKR (rs780093, 2p23.3, p = 2.2 × 10(-16)), ZBTB10 (rs440837, 8q21.13, p = 3.4 × 10(-09)), JMJD1C (rs7910927, 10q21.3, p = 6.1 × 10(-35)), SLCO1B1 (rs4149056, 12p12.1, p = 1.9 × 10(-08)), NR2F2 (rs8023580, 15q26.2, p = 8.3 × 10(-12)), ZNF652 (rs2411984, 17q21.32, p = 3.5 × 10(-14)), TDGF3 (rs1573036, Xq22.3, p = 4.1 × 10(-14)), LHCGR (rs10454142, 2p16.3, p = 1.3 × 10(-07)), BAIAP2L1 (rs3779195, 7q21.3, p = 2.7 × 10(-08)), and UGT2B15 (rs293428, 4q13.2, p = 5.5 × 10(-06)). These genes encompass multiple biologic pathways, including hepatic function, lipid metabolism, carbohydrate metabolism and T2D, androgen and estrogen receptor function, epigenetic effects, and the biology of sex steroid hormone-responsive cancers including breast and prostate cancer. We found evidence of sex-differentiated genetic influences on SHBG. In a sex-specific GWAS, the loci 4q13.2-UGT2B15 was significant in men only (men p = 2.5 × 10(-08), women p = 0.66, heterogeneity p = 0.003). Additionally, three loci showed strong sex-differentiated effects: 17p13.1-SHBG and Xq22.3-TDGF3 were stronger in men, whereas 8q21.12-ZBTB10 was stronger in women. Conditional analyses identified additional signals at the SHBG gene that together almost double the proportion of variance explained at the locus. Using an independent study of 1,129 individuals, all SNPs identified in the overall or sex-differentiated or conditional analyses explained ~15.6% and ~8.4% of the genetic variation of SHBG concentrations in men and women, respectively. The evidence for sex-differentiated effects and allelic heterogeneity highlight the importance of considering these features when estimating complex trait variance.     

BMD-associated variation at the Osterix locus is correlated with childhood obesity in females

Recent genome wide association studies (GWAS) have revealed a number of genetic variants robustly associated with bone mineral density (BMD) and/or osteoporosis. Evidence from epidemiological and clinical studies has shown an association between BMD and BMI, presumably as a consequence of bone loading. We investigated the 23 previously published BMD GWAS-derived loci in the context of childhood obesity by leveraging our existing genome-wide genotyped European American cohort of 1106 obese children (BMI ≥ 95th percentile) and 5997 controls (BMI < 95th percentile). Evidence of association was only observed at one locus, namely Osterix (SP7), with the G allele of rs2016266 being significantly over-represented among childhood obesity cases (P = 2.85 × 10(-3)). When restricting these analyses to each gender, we observed strong association between rs2016266 and childhood obesity in females (477 cases and 2867 controls; P = 3.56 × 10(-4)), which survived adjustment for all tests applied. However, no evidence of association was observed among males. Interestingly, Osterix is the only GWAS locus uncovered to date that has also been previously implicated in the determination of BMD in childhood. In conclusion, these findings indicate that a well established variant at the Osterix locus associated with increased BMD is also associated with childhood obesity primarily in females.     

Genome-wide association study identifies genetic loci associated with iron deficiency

The existence of multiple inherited disorders of iron metabolism in man, rodents and other vertebrates suggests genetic contributions to iron deficiency. To identify new genomic locations associated with iron deficiency, a genome-wide association study (GWAS) was performed using DNA collected from white men aged≥25 y and women≥50 y in the Hemochromatosis and Iron Overload Screening (HEIRS) Study with serum ferritin (SF)≤12 μg/L (cases) and iron replete controls (SF>100 μg/L in men, SF>50 μg/L in women). Regression analysis was used to examine the association between case-control status (336 cases, 343 controls) and quantitative serum iron measures and 331,060 single nucleotide polymorphism (SNP) genotypes, with replication analyses performed in a sample of 71 cases and 161 controls from a population of white male and female veterans screened at a US Veterans Affairs (VA) medical center. Five SNPs identified in the GWAS met genome-wide statistical significance for association with at least one iron measure, rs2698530 on chr. 2p14; rs3811647 on chr. 3q22, a known SNP in the transferrin (TF) gene region; rs1800562 on chr. 6p22, the C282Y mutation in the HFE gene; rs7787204 on chr. 7p21; and rs987710 on chr. 22q11 (GWAS observed P<1.51×10(-7) for all). An association between total iron binding capacity and SNP rs3811647 in the TF gene (GWAS observed P=7.0×10(-9), corrected P=0.012) was replicated within the VA samples (observed P=0.012). Associations with the C282Y mutation in the HFE gene also were replicated. The joint analysis of the HEIRS and VA samples revealed strong associations between rs2698530 on chr. 2p14 and iron status outcomes. These results confirm a previously-described TF polymorphism and implicate one potential new locus as a target for gene identification.