Pathway-Based Analysis Using Genome-wide Association Data from a Korean Non-Small Cell Lung Cancer Study

Pathway-based analysis, used in conjunction with genome-wide association study (GWAS) techniques, is a powerful tool to detect subtle but systematic patterns in genome that can help elucidate complex diseases, like cancers. Here, we stepped back from genetic polymorphisms at a single locus and examined how multiple association signals can be orchestrated to find pathways related to lung cancer susceptibility. We used single-nucleotide polymorphism (SNP) array data from 869 non-small cell lung cancer (NSCLC) cases from a previous GWAS at the National Cancer Center and 1,533 controls from the Korean Association Resource project for the pathway-based analysis. After mapping single-nucleotide polymorphisms to genes, considering their coding region and regulatory elements (±20 kbp), multivariate logistic regression of additive and dominant genetic models were fitted against disease status, with adjustments for age, gender, and smoking status. Pathway statistics were evaluated using Gene Set Enrichment Analysis (GSEA) and Adaptive Rank Truncated Product (ARTP) methods. Among 880 pathways, 11 showed relatively significant statistics compared to our positive controls (P_GSEA≤0.025, false discovery rate≤0.25). Candidate pathways were validated using the ARTP method and similarities between pathways were computed against each other. The top-ranked pathways were ABC Transporters (P_GSEA<0.001, P_ARTP = 0.001), VEGF Signaling Pathway (P_GSEA<0.001, P_ARTP = 0.008), G1/S Check Point (P_GSEA = 0.004, P_ARTP = 0.013), and NRAGE Signals Death through JNK (P_GSEA = 0.006, P_ARTP = 0.001). Our results demonstrate that pathway analysis can shed light on post-GWAS research and help identify potential targets for cancer susceptibility.     

The Influence of the CHIEF Pathway on Colorectal Cancer-Specific Mortality

Many components of the CHIEF (Convergence of Hormones, Inflammation, and Energy Related Factors) pathway could influence survival given their involvement in cell growth, apoptosis, angiogenesis, and tumor invasion stimulation. We used ARTP (Adaptive Rank Truncation Product) to test if genes in the pathway were associated with colorectal cancer-specific mortality. Colon cancer (n = 1555) and rectal cancer (n = 754) cases were followed over five years. Age, center, stage at diagnosis, and tumor molecular phenotype were considered when calculating ARTP p values. A polygenic risk score was used to summarize the magnitude of risk associated with this pathway. The JAK/STAT/SOC was significant for colon cancer survival (P_ARTP = 0.035). Fifteen genes (DUSP2, INFGR1, IL6, IRF2, JAK2, MAP3K10, MMP1, NFkB1A, NOS2A, PIK3CA, SEPX1, SMAD3, TLR2, TYK2, and VDR) were associated with colon cancer mortality (P_ARTP <0.05); JAK2 (P_ARTP  = 0.0086), PIK3CA (P_ARTP = 0.0098), and SMAD3 (P_ARTP = 0.0059) had the strongest associations. Over 40 SNPs were significantly associated with survival within the 15 significant genes (P_ARTP<0.05). SMAD3 had the strongest association with survival (HR_GG 2.46 95% CI 1.44,4.21 P_Ttrnd = 0.0002). Seven genes (IL2RA, IL8RA, IL8RB, IRF2, RAF1, RUNX3, and SEPX1) were significantly associated with rectal cancer (P_ARTP<0.05). The HR for colorectal cancer-specific mortality among colon cancer cases in the upper at-risk alleles group was 11.81 (95% CI 7.07, 19. 74) and was 10.99 (95% CI 5.30, 22.78) for rectal cancer. These results suggest that several genes in the CHIEF pathway are important for colorectal cancer survival; the risk associated with the pathway merits validation in other studies.     

Evidence for more than one Parkinson's disease-associated variant within the HLA region

Parkinson''s disease (PD) was recently found to be associated with HLA in a genome-wide association study (GWAS). Follow-up GWAS''s replicated the PD-HLA association but their top hits differ. Do the different hits tag the same locus or is there more than one PD-associated variant within HLA? We show that the top GWAS hits are not correlated with each other (0.00≤r²≤0.15). Using our GWAS (2000 cases, 1986 controls) we conducted step-wise conditional analysis on 107 SNPs with P<10⁻³ for PD-association; 103 dropped-out, four remained significant. Each SNP, when conditioned on the other three, yielded P_SNP1 = 5×10⁻⁴, P_SNP2 = 5×10⁻⁴, P_SNP3 = 4×10⁻³ and P_SNP4 = 0.025. The four SNPs were not correlated (0.01≤r²≤0.20). Haplotype analysis (excluding rare SNP2) revealed increasing PD risk with increasing risk alleles from OR = 1.27, P = 5×10⁻³ for one risk allele to OR = 1.65, P = 4×10⁻⁸ for three. Using additional 843 cases and 856 controls we replicated the independent effects of SNP1 (P_{conditioned-on-SNP4} = 0.04) and SNP4 (P_{conditioned-on-SNP1} = 0.04); SNP2 and SNP3 could not be replicated. In pooled GWAS and replication, SNP1 had OR_{conditioned-on-SNP4} = 1.23, P_{conditioned-on-SNP4} = 6×10⁻⁷; SNP4 had OR_{conditioned-on-SNP1} = 1.18, P_{conditioned-on-SNP1} = 3×10⁻³; and the haplotype with both risk alleles had OR = 1.48, P = 2×10⁻¹². Genotypic OR increased with the number of risk alleles an individual possessed up to OR = 1.94, P = 2×10⁻¹¹ for individuals who were homozygous for the risk allele at both SNP1 and SNP4. SNP1 is a variant in HLA-DRA and is associated with HLA-DRA, DRB5 and DQA2 gene expression. SNP4 is correlated (r² = 0.95) with variants that are associated with HLA-DQA2 expression, and with the top HLA SNP from the IPDGC GWAS (r² = 0.60). Our findings suggest more than one PD-HLA association; either different alleles of the same gene, or separate loci.     

Genome-wide association study for atopy and allergic rhinitis in a Singapore Chinese population

Allergic rhinitis (AR) is an atopic disease which affects about 600 million people worldwide and results from a complex interplay between genetic and environmental factors. However genetic association studies on known candidate genes yielded variable results. The aim of this study is to identify the genetic variants that influence predisposition towards allergic rhinitis in an ethnic Chinese population in Singapore using a genome-wide association study (GWAS) approach. A total of 4461 ethnic Chinese volunteers were recruited in Singapore and classified according to their allergic disease status. The GWAS included a discovery stage comparing 515 atopic cases (including 456 AR cases) and 486 non-allergic non-rhinitis (NANR) controls. The top SNPs were then validated in a replication cohort consisting of a separate 2323 atopic cases (including 676 AR cases) and 511 NANR controls. Two SNPs showed consistent association in both discovery and replication phases; MRPL4 SNP rs8111930 on 19q13.2 (OR = 0.69, P_combined = 4.46×10⁻⁰⁵) and BCAP SNP rs505010 on chromosome 10q24.1 (OR = 0.64, P_combined = 1.10×10⁻⁰⁴). In addition, we also replicated multiple associations within known candidates regions such as HLA-DQ and NPSR1 locus in the discovery phase. Our study suggests that MRPL4 and BCAP, key components of the HIF-1α and PI3K/Akt signaling pathways respectively, are two novel candidate genes for atopy and allergic rhinitis. Further study on these molecules and their signaling pathways would help in understanding of the pathogenesis of allergic rhinitis and identification of targets for new therapeutic intervention.     

A genome-wide association study of female sexual dysfunction

Background

Female sexual dysfunction (FSD) is an important but controversial problem with serious negative impact on women’s quality of life. Data from twin studies have shown a genetic contribution to the development and maintenance of FSD.

Methodology/Principal Findings

We performed a genome-wide association study (GWAS) on 2.5 million single-nucleotide polymorphisms (SNPs) in 1,104 female twins (25–81 years of age) in a population-based register and phenotypic data on lifelong sexual functioning. Although none reached conventional genome-wide level of significance (10×-8), we found strongly suggestive associations with the phenotypic dimension of arousal (rs13202860, P = 1.2×10⁻⁷; rs1876525, P = 1.2×10⁻⁷; and rs13209281 P = 8.3×10⁻⁷) on chromosome 6, around 500kb upstream of the locus HTR1E (5-hydroxytryptamine receptor 1E) locus, related to the serotonin brain pathways. We could not replicate previously reported candidate SNPs associated with FSD in the DRD4, 5HT2A and IL-1B loci.

Conclusions/Significance

We report the first GWAS of FSD symptoms in humans. This has pointed to several “risk alleles” and the implication of the serotonin and GABA pathways. Ultimately, understanding key mechanisms via this research may lead to new FSD treatments and inform clinical practice and developments in psychiatric nosology.     

Association analysis of single nucleotide polymorphisms at five loci: comparison between atopic dermatitis and asthma in the Chinese Han population

Atopic diseases, such as atopic dermatitis (AD) and asthma, are closely related to clinical phenotypes with hypersensitivity, and often share some similar genetic and pathogenic bases. Our recent GWAS identified three susceptibility gene/loci FLG (rs11204971 and rs3126085), 5q22.1 (rs10067777, rs7701890, rs13360927 and rs13361382) and 20q13.33 (rs6010620) to AD. The effect of these AD associated polymorphisms in asthma is so far unknown. To investigate whether AD relevant genetic variants is identical to asthma and reveal the differences in genetic factors between AD and asthma in Chinese Han population, seven AD associated single nucleotide polymorphisms (SNPs) as well as 3 other SNPs (rs7936562 and rs7124842 at 11q13.5 and rs4982958 at 14q11.2) from our previous AD GWAS were genotyped in 463 asthma patients and 985 controls using Sequenom MassArray system. We found rs4982958 at 14q11.2 was significantly associated with asthma (P = 3.04×10⁻⁴, OR = 0.73). We also detected one significant risk haplotype GGGA from the 4 SNPs (rs10067777, rs7701890, rs13360927 and rs13361382) at 5q22.1 in AD cases (P _correction = 3.60×10⁻¹⁰, OR = 1.26), and the haplotype was suggestive of risk in asthma cases in this study (P = 0.014, P _correction = 0.084, OR = 1.38). These SNPs (rs11204971, rs3126085, rs7936562, rs712484 and rs6010620) at AD susceptibility genes/loci FLG, 11q13.5 and 20q13.33 were not associated with asthma in this study. Our results further comfirmed that 14q11.2 was an important candidate locus for asthma and demonstrated that 5q22.1 might be shared by AD and asthma in Chinese Han population.     

Genome-wide association study identifies HLA-DP as a susceptibility gene for pediatric asthma in Asian populations

Asthma is a complex phenotype influenced by genetic and environmental factors. We conducted a genome-wide association study (GWAS) with 938 Japanese pediatric asthma patients and 2,376 controls. Single-nucleotide polymorphisms (SNPs) showing strong associations (P<1×10⁻⁸) in GWAS were further genotyped in an independent Japanese samples (818 cases and 1,032 controls) and in Korean samples (835 cases and 421 controls). SNP rs987870, located between HLA-DPA1 and HLA-DPB1, was consistently associated with pediatric asthma in 3 independent populations (P _combined = 2.3×10⁻¹⁰, odds ratio [OR] = 1.40). HLA-DP allele analysis showed that DPA1*0201 and DPB1*0901, which were in strong linkage disequilibrium, were strongly associated with pediatric asthma (DPA1*0201: P = 5.5×10⁻¹⁰, OR = 1.52, and DPB1*0901: P = 2.0×10⁻⁷, OR = 1.49). Our findings show that genetic variants in the HLA-DP locus are associated with the risk of pediatric asthma in Asian populations.     

CD3Z genetic polymorphism in immune response to hepatitis B vaccination in two independent Chinese populations

Vaccination against hepatitis B virus is an effective and routine practice that can prevent infection. However, vaccine-induced immunity to hepatitis B varies among individuals. CD4⁺ T helper cells, which play an important role in both cellular and humoral immunity, are involved in the immune response elicited by vaccination. Polymorphisms in the genes involved in stimulating the activation and proliferation of CD4⁺ T helper cells may influence the immune response to hepatitis B vaccination. In the first stage of the present study, a total of 111 single nucleotide polymorphisms (SNPs) in 17 genes were analyzed, using the iPLEX MassARRAY system, among 214 high responders and 107 low responders to hepatitis B vaccination. Three SNPs (rs12133337 and rs10918706 in CD3Z, rs10912564 in OX40L) were associated significantly with the immune response to hepatitis B vaccination (P = 0.008, 0.041, and 0.019, respectively). The three SNPs were analyzed further with the TaqMan-MGB or TaqMan-BHQ probe-based real-time polymerase chain reaction in another independent population, which included 1090 high responders and 636 low responders. The minor allele ‘C’ of rs12133337 continued to show an association with a lower response to hepatitis B vaccination (P = 0.033, odds radio = 1.28, 95% confidence interval = 1.01–1.61). Furthermore, in the stratified analysis for both the first and second populations, the association of the minor allele ‘C’ of rs12133337 with a lower response to hepatitis B vaccination was more prominent after individuals who were overweight or obese (body mass index ≥25 kg/m²) were excluded (1^st stage: P = 0.003, 2^nd stage: P = 0.002, P-combined = 9.47e-5). These findings suggest that the rs12133337 polymorphism in the CD3Z gene might affect the immune response to hepatitis B vaccination, and that a lower BMI might increase the contribution of the polymorphism to immunity to hepatitis B vaccination.     

Genome-wide gene-environment study identifies glutamate receptor gene GRIN2A as a Parkinson's disease modifier gene via interaction with coffee

Our aim was to identify genes that influence the inverse association of coffee with the risk of developing Parkinson''s disease (PD). We used genome-wide genotype data and lifetime caffeinated-coffee-consumption data on 1,458 persons with PD and 931 without PD from the NeuroGenetics Research Consortium (NGRC), and we performed a genome-wide association and interaction study (GWAIS), testing each SNP''s main-effect plus its interaction with coffee, adjusting for sex, age, and two principal components. We then stratified subjects as heavy or light coffee-drinkers and performed genome-wide association study (GWAS) in each group. We replicated the most significant SNP. Finally, we imputed the NGRC dataset, increasing genomic coverage to examine the region of interest in detail. The primary analyses (GWAIS, GWAS, Replication) were performed using genotyped data. In GWAIS, the most significant signal came from rs4998386 and the neighboring SNPs in GRIN2A. GRIN2A encodes an NMDA-glutamate-receptor subunit and regulates excitatory neurotransmission in the brain. Achieving P_2df = 10⁻⁶, GRIN2A surpassed all known PD susceptibility genes in significance in the GWAIS. In stratified GWAS, the GRIN2A signal was present in heavy coffee-drinkers (OR = 0.43; P = 6×10⁻⁷) but not in light coffee-drinkers. The a priori Replication hypothesis that “Among heavy coffee-drinkers, rs4998386_T carriers have lower PD risk than rs4998386_CC carriers” was confirmed: OR_Replication = 0.59, P_Replication = 10⁻³; OR_Pooled = 0.51, P_Pooled = 7×10⁻⁸. Compared to light coffee-drinkers with rs4998386_CC genotype, heavy coffee-drinkers with rs4998386_CC genotype had 18% lower risk (P = 3×10⁻³), whereas heavy coffee-drinkers with rs4998386_TC genotype had 59% lower risk (P = 6×10⁻¹³). Imputation revealed a block of SNPs that achieved P_2df<5×10⁻⁸ in GWAIS, and OR = 0.41, P = 3×10⁻⁸ in heavy coffee-drinkers. This study is proof of concept that inclusion of environmental factors can help identify genes that are missed in GWAS. Both adenosine antagonists (caffeine-like) and glutamate antagonists (GRIN2A-related) are being tested in clinical trials for treatment of PD. GRIN2A may be a useful pharmacogenetic marker for subdividing individuals in clinical trials to determine which medications might work best for which patients.     

SEPP1 Influences Breast Cancer Risk among Women with Greater Native American Ancestry: The Breast Cancer Health Disparities Study

Selenoproteins are a class of proteins containing a selenocysteine residue, many of which have been shown to have redox functions, acting as antioxidants to decrease oxidative stress. Selenoproteins have previously been associated with risk of various cancers and redox-related diseases. In this study we evaluated possible associations between breast cancer risk and survival and single nucleotide polymorphisms (SNPs) in the selenoprotein genes GPX1, GPX2, GPX3, GPX4, SELS, SEP15, SEPN1, SEPP1, SEPW1, TXNRD1, and TXNRD2 among Hispanic/Native American (2111 cases, 2597 controls) and non-Hispanic white (NHW) (1481 cases, 1586 controls) women in the Breast Cancer Health Disparities Study. Adaptive Rank Truncated Product (ARTP) analysis was used to determine both gene and pathway significance with these genes. The overall selenoprotein pathway P_ARTP was not significantly associated with breast cancer risk (P_ARTP = 0.69), and only one gene, GPX3, was of borderline significance for the overall population (P_ARTP =0.09) and marginally significant among women with 0-28% Native American (NA) ancestry (P_ARTP=0.06). The SEPP1 gene was statistically significantly associated with breast cancer risk among women with higher NA ancestry (P_ARTP=0.002) and contributed to a significant pathway among those women (P_ARTP=0.04). GPX1, GPX3, and SELS were associated with Estrogen Receptor-/Progesterone Receptor+ status (P_ARTP = 0.002, 0.05, and 0.01, respectively). Four SNPs (GPX3 rs2070593, rsGPX4 rs2074451, SELS rs9874, and TXNRD1 rs17202060) significantly interacted with dietary oxidative balance score after adjustment for multiple comparisons to alter breast cancer risk. GPX4 was significantly associated with breast cancer survival among those with the highest NA ancestry (P_ARTP = 0.05) only. Our data suggest that SEPP1 alters breast cancer risk among women with higher levels of NA ancestry.     

Both Baseline Clinical Factors and Genetic Polymorphisms Influence the Development of Severe Functional Status in Ankylosing Spondylitis

Functional severity in ankylosing spondylitis (AS) patients is variable and difficult to predict early. The aim of our study was to assess whether a combination of baseline clinical factors and genetic markers may predict the development of severe functional status in AS. We performed a cross-sectional association study on AS patients included in the Spanish National Registry of Spondyloarthropathies—REGISPONSER. Bath Ankylosing Spondylitis Functional Index (BASFI) was standardized by adjusting for disease duration since the first symptoms (BASFI/t). We considered as severe functional status the values of BASFI/t in the top of the 60th (p60), 65th (p65), 70th (p70), and 75th (p75) percentile. We selected 384 single nucleotide polymorphisms (SNPs) distributed in 190 genes to be analyzed. The study cohort included 456 patients with mean age 50.8(±10.5) years and with mean disease duration since first symptoms 24.7 (±10.1) years. Older age at disease onset and neck pain at baseline showed statistical significant association with severe BASFI/t. Polymorphisms associated in the allele frequencies test with severe BASFI/t in all classifications were: rs2542151 (p60 [P = .04], p65 [P = .04], p70 [P = .001] and p75 [P = .001]) and rs2254441 (p60 [P = .004], p65 [P = .02], p70 [P = .01] and p75 [P<.001]).. Genotype association, after adjustment for covariates, found an association in three of the four patients'' classifications for rs2542151 and in two of the classifications for rs2254441.Forward logistic regression did not identify any model with a good predictive power for severe functional development.In our study we identified clinical factors and 24 polymorphisms associated with development of severe functional status in AS patients. Validation of these results in other cohorts is required.     

Genetic polymorphisms in the hypothalamic pathway in relation to subsequent weight change--the DiOGenes study

Background

Single nucleotide polymorphisms (SNPs) in genes encoding the components involved in the hypothalamic pathway may influence weight gain and dietary factors may modify their effects.

Aim

We conducted a case-cohort study to investigate the associations of SNPs in candidate genes with weight change during an average of 6.8 years of follow-up and to examine the potential effect modification by glycemic index (GI) and protein intake.

Methods and Findings

Participants, aged 20–60 years at baseline, came from five European countries. Cases (‘weight gainers’) were selected from the total eligible cohort (n = 50,293) as those with the greatest unexplained annual weight gain (n = 5,584). A random subcohort (n = 6,566) was drawn with the intention to obtain an equal number of cases and noncases (n = 5,507). We genotyped 134 SNPs that captured all common genetic variation across the 15 candidate genes; 123 met the quality control criteria. Each SNP was tested for association with the risk of being a ‘weight gainer’ (logistic regression models) in the case-noncase data and with weight gain (linear regression models) in the random subcohort data. After accounting for multiple testing, none of the SNPs was significantly associated with weight change. Furthermore, we observed no significant effect modification by dietary factors, except for SNP rs7180849 in the neuromedin β gene (NMB). Carriers of the minor allele had a more pronounced weight gain at a higher GI (P = 2×10⁻⁷).

Conclusions

We found no evidence of association between SNPs in the studied hypothalamic genes with weight change. The interaction between GI and NMB SNP rs7180849 needs further confirmation.     

Genome-Wide Association Study of Gene by Smoking Interactions in Coronary Artery Calcification

Many GWAS have identified novel loci associated with common diseases, but have focused only on main effects of individual genetic variants rather than interactions with environmental factors (GxE). Identification of GxE interactions is particularly important for coronary heart disease (CHD), a major preventable source of morbidity and mortality with strong non-genetic risk factors. Atherosclerosis is the major cause of CHD, and coronary artery calcification (CAC) is directly correlated with quantity of coronary atherosclerotic plaque. In the current study, we tested for genetic variants influencing extent of CAC via interaction with smoking (GxS), by conducting a GxS discovery GWAS in Genetic Epidemiology Network of Arteriopathy (GENOA) sibships (N = 915 European Americans) followed by replication in Framingham Heart Study (FHS) sibships (N = 1025 European Americans). Generalized estimating equations accounted for the correlation within sibships in strata-specific groups of smokers and nonsmokers, as well as GxS interaction. Primary analysis found SNPs that showed suggestive associations (p≤10⁻⁵) in GENOA GWAS, but these index SNPs did not replicate in FHS. However, secondary analysis was able to replicate candidate gene regions in FHS using other SNPs (+/−250 kb of GENOA index SNP). In smoker and nonsmoker groups, replicated genes included TCF7L2 (p = 6.0×10⁻⁵) and WWOX (p = 4.5×10⁻⁶); and TNFRSF8 (p = 7.8×10⁻⁵), respectively. For GxS interactions, replicated genes included TBC1D4 (p = 6.9×10⁻⁵) and ADAMTS9 (P = 7.1×10⁻⁵). Interestingly, these genes are involved in inflammatory pathways mediated by the NF-κB axis. Since smoking is known to induce chronic and systemic inflammation, association of these genes likely reflects roles in CAC development via inflammatory pathways. Furthermore, the NF-κB axis regulates bone remodeling, a key physiological process in CAC development. In conclusion, GxS GWAS has yielded evidence for novel loci that are associated with CAC via interaction with smoking, providing promising new targets for future population-based and functional studies of CAC development.     

Different genes interact with particulate matter and tobacco smoke exposure in affecting lung function decline in the general population

Background

Oxidative stress related genes modify the effects of ambient air pollution or tobacco smoking on lung function decline. The impact of interactions might be substantial, but previous studies mostly focused on main effects of single genes.

Objectives

We studied the interaction of both exposures with a broad set of oxidative-stress related candidate genes and pathways on lung function decline and contrasted interactions between exposures.

Methods

For 12679 single nucleotide polymorphisms (SNPs), change in forced expiratory volume in one second (FEV₁), FEV₁ over forced vital capacity (FEV₁/FVC), and mean forced expiratory flow between 25 and 75% of the FVC (FEF_25-75) was regressed on interval exposure to particulate matter <10 µm in diameter (PM10) or packyears smoked (a), additive SNP effects (b), and interaction terms between (a) and (b) in 669 adults with GWAS data. Interaction p-values for 152 genes and 14 pathways were calculated by the adaptive rank truncation product (ARTP) method, and compared between exposures. Interaction effect sizes were contrasted for the strongest SNPs of nominally significant genes (p_interaction<0.05). Replication was attempted for SNPs with MAF>10% in 3320 SAPALDIA participants without GWAS.

Results

On the SNP-level, rs2035268 in gene SNCA accelerated FEV₁/FVC decline by 3.8% (p_interaction = 2.5×10⁻⁶), and rs12190800 in PARK2 attenuated FEV1 decline by 95.1 ml p_interaction = 9.7×10⁻⁸) over 11 years, while interacting with PM10. Genes and pathways nominally interacting with PM10 and packyears exposure differed substantially. Gene CRISP2 presented a significant interaction with PM10 (p_interaction = 3.0×10⁻⁴) on FEV₁/FVC decline. Pathway interactions were weak. Replications for the strongest SNPs in PARK2 and CRISP2 were not successful.

Conclusions

Consistent with a stratified response to increasing oxidative stress, different genes and pathways potentially mediate PM10 and tobacco smoke effects on lung function decline. Ignoring environmental exposures would miss these patterns, but achieving sufficient sample size and comparability across study samples is challenging.     

A comprehensive evaluation of potential lung function associated genes in the SpiroMeta general population sample

Rationale

Lung function measures are heritable traits that predict population morbidity and mortality and are essential for the diagnosis of chronic obstructive pulmonary disease (COPD). Variations in many genes have been reported to affect these traits, but attempts at replication have provided conflicting results. Recently, we undertook a meta-analysis of Genome Wide Association Study (GWAS) results for lung function measures in 20,288 individuals from the general population (the SpiroMeta consortium).

Objectives

To comprehensively analyse previously reported genetic associations with lung function measures, and to investigate whether single nucleotide polymorphisms (SNPs) in these genomic regions are associated with lung function in a large population sample.

Methods

We analysed association for SNPs tagging 130 genes and 48 intergenic regions (+/−10 kb), after conducting a systematic review of the literature in the PubMed database for genetic association studies reporting lung function associations.

Results

The analysis included 16,936 genotyped and imputed SNPs. No loci showed overall significant association for FEV₁ or FEV₁/FVC traits using a carefully defined significance threshold of 1.3×10⁻⁵. The most significant loci associated with FEV₁ include SNPs tagging MACROD2 (P = 6.81×10⁻⁵), CNTN5 (P = 4.37×10⁻⁴), and TRPV4 (P = 1.58×10⁻³). Among ever-smokers, SERPINA1 showed the most significant association with FEV₁ (P = 8.41×10⁻⁵), followed by PDE4D (P = 1.22×10⁻⁴). The strongest association with FEV₁/FVC ratio was observed with ABCC1 (P = 4.38×10⁻⁴), and ESR1 (P = 5.42×10⁻⁴) among ever-smokers.

Conclusions

Polymorphisms spanning previously associated lung function genes did not show strong evidence for association with lung function measures in the SpiroMeta consortium population. Common SERPINA1 polymorphisms may affect FEV₁ among smokers in the general population.     

ABCC5, a Gene That Influences the Anterior Chamber Depth,Is Associated with Primary Angle Closure Glaucoma

Anterior chamber depth (ACD) is a key anatomical risk factor for primary angle closure glaucoma (PACG). We conducted a genome-wide association study (GWAS) on ACD to discover novel genes for PACG on a total of 5,308 population-based individuals of Asian descent. Genome-wide significant association was observed at a sequence variant within ABCC5 (rs1401999; per-allele effect size = −0.045 mm, P = 8.17×10⁻⁹). This locus was associated with an increase in risk of PACG in a separate case-control study of 4,276 PACG cases and 18,801 controls (per-allele OR = 1.13 [95% CI: 1.06–1.22], P = 0.00046). The association was strengthened when a sub-group of controls with open angles were included in the analysis (per-allele OR = 1.30, P = 7.45×10⁻⁹; 3,458 cases vs. 3,831 controls). Our findings suggest that the increase in PACG risk could in part be mediated by genetic sequence variants influencing anterior chamber dimensions.     

GWAS identifies novel susceptibility loci on 6p21.32 and 21q21.3 for hepatocellular carcinoma in chronic hepatitis B virus carriers

Genome-wide association studies (GWAS) have recently identified KIF1B as susceptibility locus for hepatitis B virus (HBV)–related hepatocellular carcinoma (HCC). To further identify novel susceptibility loci associated with HBV–related HCC and replicate the previously reported association, we performed a large three-stage GWAS in the Han Chinese population. 523,663 autosomal SNPs in 1,538 HBV–positive HCC patients and 1,465 chronic HBV carriers were genotyped for the discovery stage. Top candidate SNPs were genotyped in the initial validation samples of 2,112 HBV–positive HCC cases and 2,208 HBV carriers and then in the second validation samples of 1,021 cases and 1,491 HBV carriers. We discovered two novel associations at rs9272105 (HLA-DQA1/DRB1) on 6p21.32 (OR = 1.30, P = 1.13×10⁻¹⁹) and rs455804 (GRIK1) on 21q21.3 (OR = 0.84, P = 1.86×10⁻⁸), which were further replicated in the fourth independent sample of 1,298 cases and 1,026 controls (rs9272105: OR = 1.25, P = 1.71×10⁻⁴; rs455804: OR = 0.84, P = 6.92×10⁻³). We also revealed the associations of HLA-DRB1*0405 and 0901*0602, which could partially account for the association at rs9272105. The association at rs455804 implicates GRIK1 as a novel susceptibility gene for HBV–related HCC, suggesting the involvement of glutamate signaling in the development of HBV–related HCC.