Significant eQTL Research Articles

Genetic Variants in the Bone Morphogenic Protein Gene Family Modify the Association between Residential Exposure to Traffic and Peripheral Arterial Disease.

There is a growing literature indicating that genetic variants modify many of the associations between environmental exposures and clinical outcomes, potentially by increasing susceptibility to these exposures. However, genome-scale investigations of these interactions have been rarely performed particularly in the case of air pollution exposures. We performed race-stratified genome-wide gene-environment interaction association studies on European-American (EA, N = 1623) and African-American (AA, N = 554) cohorts to investigate the joint influence of common single nucleotide polymorphisms (SNPs) and residential exposure to traffic (“traffic exposure”)—a recognized vascular disease risk factor—on peripheral arterial disease (PAD). Traffic exposure was estimated via the distance from the primary residence to the nearest major roadway, defined as the nearest limited access highways or major arterial. The rs755249-traffic exposure interaction was associated with PAD at a genome-wide significant level (P = 2.29x10-8) in European-Americans. Rs755249 is located in the 3’ untranslated region of BMP8A, a member of the bone morphogenic protein (BMP) gene family. Further investigation revealed several variants in BMP genes associated with PAD via an interaction with traffic exposure in both the EA and AA cohorts; this included interactions with non-synonymous variants in BMP2, which is regulated by air pollution exposure. The BMP family of genes is linked to vascular growth and calcification and is a novel gene family for the study of PAD pathophysiology. Further investigation of BMP8A using the Genotype Tissue Expression Database revealed multiple variants with nominally significant (P < 0.05) interaction P-values in our EA cohort were significant BMP8A eQTLs in tissue types highlight relevant for PAD such as rs755249 (tibial nerve, eQTL P = 3.6x10-6) and rs1180341 (tibial artery, eQTL P = 5.3x10-6). Together these results reveal a novel gene, and possibly gene family, associated with PAD via an interaction with traffic air pollution exposure. These results also highlight the potential for interactions studies, particularly at the genome scale, to reveal novel biology linking environmental exposures to clinical outcomes.

IFORM/eQTL: an ultrahigh-dimensional platform for inferring the global genetic architecture of gene transcripts.

Knowledge about how changes in gene expression are encoded by expression quantitative trait loci (eQTLs) is a key to construct the genotype-phenotype map for complex traits or diseases. Traditional eQTL mapping is to associate one transcript with a single marker at a time, thereby limiting our inference about a complete picture of the genetic architecture of gene expression. Here, we implemented an ultrahigh-dimensional variable selection model to build a computing platform that can systematically scan main effects and interaction effects among all possible loci and identify a set of significant eQTLs modulating differentiation and function of gene expression. This platform, named iFORM/eQTL, was assembled by forward-selection-based procedures to tackle complex covariance structures of gene-gene interactions. iFORM/eQTL can particularly discern the role of cis-QTLs, trans-QTLs and their epistatic interactions in gene expression. Results from the reanalysis of a published genetic and genomic data set through iFORM/eQTL gain new discoveries on the genetic origin of gene expression differentiation in Caenorhabditis elegans, which could not be detected by a traditional one-locus/one-transcript analysis approach.

Chromatin interactions and candidate genes at ten prostate cancer risk loci.

Genome-wide association studies have identified more than 100 common single nucleotide polymorphisms (SNPs) that are associated with prostate cancer risk. However, the vast majority of these SNPs lie in noncoding regions of the genome. To test whether these risk SNPs regulate their target genes through long-range chromatin interactions, we applied capture-based 3C sequencing technology to investigate possible cis-interactions at ten prostate cancer risk loci in six cell lines. We identified significant physical interactions between risk regions and their potential target genes including CAPG at 2p11.2, C2orf43 at 2p24.1, RFX6 at 6q22.1, NFASC at 1q32.1, MYC at 8q24.1 and AGAP7P at 10q11.23. Most of the interaction peaks were co-localized to regions of active histone modification and transcription factor binding sites. Expression quantitative trait locus (eQTL) analysis showed suggestive eQTL signals at rs1446669, rs699664 and rs1078004 for CAPG (p < 0.004), rs13394027 for C2orf43 (p = 2.25E-27), rs10993994 and rs4631830 for AGAP7P (p < 8.02E-5). Further analysis revealed an enhancer activity at genomic region surrounding rs4631830 which was expected to disrupt HOXB-like DNA binding affinity. This study identifies a set of candidate genes and their potential regulatory variants, and provides additional evidence showing the role of long-range chromatin interactions in prostate cancer etiology.

Alzheimer's Disease Risk Polymorphisms Regulate Gene Expression in the ZCWPW1 and the CELF1 Loci.

Late onset Alzheimer’s disease (LOAD) is a genetically complex and clinically heterogeneous disease. Recent large-scale genome wide association studies (GWAS) have identified more than twenty loci that modify risk for AD. Despite the identification of these loci, little progress has been made in identifying the functional variants that explain the association with AD risk. Thus, we sought to determine whether the novel LOAD GWAS single nucleotide polymorphisms (SNPs) alter expression of LOAD GWAS genes and whether expression of these genes is altered in AD brains. The majority of LOAD GWAS SNPs occur in gene dense regions under large linkage disequilibrium (LD) blocks, making it unclear which gene(s) are modified by the SNP. Thus, we tested for brain expression quantitative trait loci (eQTLs) between LOAD GWAS SNPs and SNPs in high LD with the LOAD GWAS SNPs in all of the genes within the GWAS loci. We found a significant eQTL between rs1476679 and PILRB and GATS, which occurs within the ZCWPW1 locus. PILRB and GATS expression levels, within the ZCWPW1 locus, were also associated with AD status. Rs7120548 was associated with MTCH2 expression, which occurs within the CELF1 locus. Additionally, expression of several genes within the CELF1 locus, including MTCH2, were highly correlated with one another and were associated with AD status. We further demonstrate that PILRB, as well as other genes within the GWAS loci, are most highly expressed in microglia. These findings together with the function of PILRB as a DAP12 receptor supports the critical role of microglia and neuroinflammation in AD risk.

Association of Forced Vital Capacity with the Developmental Gene NCOR2.

BackgroundForced Vital Capacity (FVC) is an important predictor of all-cause mortality in the absence of chronic respiratory conditions. Epidemiological evidence highlights the role of early life factors on adult FVC, pointing to environmental exposures and genes affecting lung development as risk factors for low FVC later in life. Although highly heritable, a small number of genes have been found associated with FVC, and we aimed at identifying further genetic variants by focusing on lung development genes.MethodsPer-allele effects of 24,728 SNPs in 403 genes involved in lung development were tested in 7,749 adults from three studies (NFBC1966, ECRHS, EGEA). The most significant SNP for the top 25 genes was followed-up in 46,103 adults (CHARGE and SpiroMeta consortia) and 5,062 children (ALSPAC). Associations were considered replicated if the replication p-value survived Bonferroni correction (p<0.002; 0.05/25), with a nominal p-value considered as suggestive evidence. For SNPs with evidence of replication, effects on the expression levels of nearby genes in lung tissue were tested in 1,111 lung samples (Lung eQTL consortium), with further functional investigation performed using public epigenomic profiling data (ENCODE).ResultsNCOR2-rs12708369 showed strong replication in children (p = 0.0002), with replication unavailable in adults due to low imputation quality. This intronic variant is in a strong transcriptional enhancer element in lung fibroblasts, but its eQTL effects could not be tested due to low imputation quality in the eQTL dataset. SERPINE2-rs6754561 replicated at nominal level in both adults (p = 0.036) and children (p = 0.045), while WNT16-rs2707469 replicated at nominal level only in adults (p = 0.026). The eQTL analyses showed association of WNT16-rs2707469 with expression levels of the nearby gene CPED1. We found no statistically significant eQTL effects for SERPINE2-rs6754561.ConclusionsWe have identified a new gene, NCOR2, in the retinoic acid signalling pathway pointing to a role of vitamin A metabolism in the regulation of FVC. Our findings also support SERPINE2, a COPD gene with weak previous evidence of association with FVC, and suggest WNT16 as a further promising candidate.

Identification of candidate genes for prostate cancer-risk SNPs utilizing a normal prostate tissue eQTL data set.

Multiple studies have identified loci associated with the risk of developing prostate cancer but the associated genes are not well studied. Here we create a normal prostate tissue-specific eQTL data set and apply this data set to previously identified prostate cancer (PrCa)-risk SNPs in an effort to identify candidate target genes. The eQTL data set is constructed by the genotyping and RNA sequencing of 471 samples. We focus on 146 PrCa-risk SNPs, including all SNPs in linkage disequilibrium with each risk SNP, resulting in 100 unique risk intervals. We analyse cis-acting associations where the transcript is located within 2 Mb (±1 Mb) of the risk SNP interval. Of all SNP–gene combinations tested, 41.7% of SNPs demonstrate a significant eQTL signal after adjustment for sample histology and 14 expression principal component covariates. Of the 100 PrCa-risk intervals, 51 have a significant eQTL signal and these are associated with 88 genes. This study provides a rich resource to study biological mechanisms underlying genetic risk to PrCa.

Comprehensively Evaluating cis-Regulatory Variation in the Human Prostate Transcriptome by Using Gene-Level Allele-Specific Expression

The identification of cis-acting regulatory variation in primary tissues has the potential to elucidate the genetic basis of complex traits and further our understanding of transcriptomic diversity across cell types. Expression quantitative trait locus (eQTL) association analysis using RNA sequencing (RNA-seq) data can improve upon the detection of cis-acting regulatory variation by leveraging allele-specific expression (ASE) patterns in association analysis. Here, we present a comprehensive evaluation of cis-acting eQTLs by analyzing RNA-seq gene-expression data and genome-wide high-density genotypes from 471 samples of normal primary prostate tissue. Using statistical models that integrate ASE information, we identified extensive cis-eQTLs across the prostate transcriptome and found that approximately 70% of expressed genes corresponded to a significant eQTL at a gene-level false-discovery rate of 0.05. Overall, cis-eQTLs were heavily concentrated near the transcription start and stop sites of affected genes, and effects were negatively correlated with distance. We identified multiple instances of cis-acting co-regulation by using phased genotype data and discovered 233 SNPs as the most strongly associated eQTLs for more than one gene. We also noted significant enrichment (25/50, p = 2E-5) of previously reported prostate cancer risk SNPs in prostate eQTLs. Our results illustrate the benefit of assessing ASE data in cis-eQTL analyses by showing better reproducibility of prior eQTL findings than of eQTL mapping based on total expression alone. Altogether, our analysis provides extensive functional context of thousands of SNPs in prostate tissue, and these results will be of critical value in guiding studies examining disease of the human prostate.

Expression quantitative trait loci (eQTL) mapping in Puerto Rican children.

BackgroundExpression quantitative trait loci (eQTL) have been identified using tissue or cell samples from diverse human populations, thus enhancing our understanding of regulation of gene expression. However, few studies have attempted to identify eQTL in racially admixed populations such as Hispanics.MethodsWe performed a systematic eQTL study to identify regulatory variants of gene expression in whole blood from 121 Puerto Rican children with (n = 63) and without (n = 58) asthma. Genome-wide genotyping was conducted using the Illumina Omni2.5M Bead Chip, and gene expression was assessed using the Illumina HT-12 microarray. After completing quality control, we performed a pair-wise genome analysis of ~15 K transcripts and ~1.3 M SNPs for both local and distal effects. This analysis was conducted under a regression framework adjusting for age, gender and principal components derived from both genotypic and mRNA data. We used a false discovery rate (FDR) approach to identify significant eQTL signals, which were next compared to top eQTL signals from existing eQTL databases. We then performed a pathway analysis for our top genes.ResultsWe identified 36,720 local pairs in 3,391 unique genes and 1,851 distal pairs in 446 unique genes at FDR <0.05, corresponding to unadjusted P values lower than 1.5x10-4 and 4.5x10-9, respectively. A significant proportion of genes identified in our study overlapped with those identified in previous studies. We also found an enrichment of disease-related genes in our eQTL list.ConclusionsWe present results from the first eQTL study in Puerto Rican children, who are members of a unique Hispanic cohort disproportionately affected with asthma, prematurity, obesity and other common diseases. Our study confirmed eQTL signals identified in other ethnic groups, while also detecting additional eQTLs unique to our study population. The identified eQTLs will help prioritize findings from future genome-wide association studies in Puerto Ricans.

Genetic control of gene expression at novel and established chronic obstructive pulmonary disease loci.

Genetic risk loci have been identified for a wide range of diseases through genome-wide association studies (GWAS), but the relevant functional mechanisms have been identified for only a small proportion of these GWAS-identified loci. By integrating results from the largest current GWAS of chronic obstructive disease (COPD) with expression quantitative trait locus (eQTL) analysis in whole blood and sputum from 121 subjects with COPD from the ECLIPSE Study, this analysis identifies loci that are simultaneously associated with COPD and the expression of nearby genes (COPD eQTLs). After integrative analysis, 19 COPD eQTLs were identified, including all four previously identified genome-wide significant loci near HHIP, FAM13A, and the 15q25 and 19q13 loci. For each COPD eQTL, fine mapping and colocalization analysis to identify causal shared eQTL and GWAS variants identified a subset of sites with moderate-to-strong evidence of harboring at least one shared variant responsible for both the eQTL and GWAS signals. Transcription factor binding site (TFBS) analysis confirms that multiple COPD eQTL lead SNPs disrupt TFBS, and enhancer enrichment analysis for loci with the strongest colocalization signals showed enrichment for blood-related cell types (CD3 and CD4+ T cells, lymphoblastoid cell lines). In summary, integrative eQTL and GWAS analysis confirms that genetic control of gene expression plays a key role in the genetic architecture of COPD and identifies specific blood-related cell types as likely participants in the functional pathway from GWAS-associated variant to disease phenotype.

A meta-analysis of gene expression quantitative trait loci in brain.

Current catalogs of brain expression quantitative trait loci (eQTL) are incomplete and the findings do not replicate well across studies. All existing cortical eQTL studies are small and emphasize the need for a meta-analysis. We performed a meta-analysis of 424 brain samples across five studies to identify regulatory variants influencing gene expression in human cortex. We identified 3584 genes in autosomes and chromosome X with false discovery rate q<0.05 whose expression was significantly associated with DNA sequence variation. Consistent with previous eQTL studies, local regulatory variants tended to occur symmetrically around transcription start sites and the effect was more evident in studies with large sample sizes. In contrast to random SNPs, we observed that significant eQTLs were more likely to be near 5'-untranslated regions and intersect with regulatory features. Permutation-based enrichment analysis revealed that SNPs associated with schizophrenia and bipolar disorder were enriched among brain eQTLs. Genes with significant eQTL evidence were also strongly associated with diseases from OMIM (Online Mendelian Inheritance in Man) and the NHGRI (National Human Genome Research Institute) genome-wide association study catalog. Surprisingly, we found that a large proportion (28%) of ~1000 autosomal genes encoding proteins needed for mitochondrial structure or function were eQTLs (enrichment P-value=1.3 × 10−9), suggesting a potential role for common genetic variation influencing the robustness of energy supply in brain and a possible role in the etiology of some psychiatric disorders. These systematically generated eQTL information should be a valuable resource in determining the functional mechanisms of brain gene expression and the underlying biology of associations with psychiatric disorders.

Genetic variability in the regulation of gene expression in ten regions of the human brain.

Germ-line genetic control of gene expression occurs via expression quantitative trait loci (eQTLs). We present a large, exon-specific eQTL data set covering ten human brain regions. We found that cis-eQTL signals (within 1 Mb of their target gene) were numerous, and many acted heterogeneously among regions and exons. Co-regulation analysis of shared eQTL signals produced well-defined modules of region-specific co-regulated genes, in contrast to standard coexpression analysis of the same samples. We report cis-eQTL signals for 23.1% of catalogued genome-wide association study hits for adult-onset neurological disorders. The data set is publicly available via public data repositories and via http://www.braineac.org/. Our study increases our understanding of the regulation of gene expression in the human brain and will be of value to others pursuing functional follow-up of disease-associated variants.

Coeliac disease-associated polymorphisms influence thymic gene expression.

Significant associations between coeliac disease (CD) and single nucleotide polymorphisms (SNPs) distributed over 40 genetic regions have been established. The majority of these SNPs are non-coding and 20 SNPs were, by expression quantitative trait loci (eQTL) analysis, found to harbour cis regulatory potential in peripheral blood mononuclear cells (PBMC). Almost all regions contain genes with an immunological relevant function, of which many act in the same biological pathways. One such pathway is T-cell development in the thymus, a pathway previously not explored in CD pathogenesis. The aim of our study was to explore the regulatory potential of the CD-associated SNPs (n=50) by eQTL analysis in thymic tissue from 42 subjects. In total, 43 nominal significant (P<0.05) eQTLs were found within 24 CD-associated chromosomal regions, corresponding to 27 expression-altering SNPs (eSNPs) and 40 probes (eProbes) that represents 39 unique genes (eGenes). Nine significant probe-SNP pairs (corresponding to 8 eSNPs and 7 eGenes) overlapped with previous findings in PBMC (rs12727642-PARK7, rs296547-DDX59, rs917997-IL18RAP, rs842647-AHSA2, rs13003464-AHSA2, rs6974491-ELMO1, rs2074404-NSF (two independent probes) and rs2298428-UBE2L3). When compared across more tissues, we found that 14 eQTLs could represent potentially novel thymus-specific eQTLs. This implies that CD risk polymorphisms could affect gene regulation in thymus.

720 Rare Variants in CRC Susceptibility Genes Do Not Predispose to Serrated Polyposis Syndrome

G A A b st ra ct s both Crohn's disease (CD) and ulcerative colitis (UC), were significant eQTLs for INPP5E and UBA7, respectively. The CD-associated SNP rs2549794 was an eQTL for ERAP2, while 2 UC-associated SNPs rs661946 and rs9847710 were eQTLs for NXPE1, REXO2 (both associated with rs661946), and SFMBT1, respectively. Finally, the CRC-associated SNP rs3802842 was noted to be a significant eQTL for C11orf93. Conclusion: eQTL mapping in the human colon reveals a large number of significant cis eQTLs. A number of SNPs associated with CRC and IBD were noted to be eQTLs in the colon. Further studies are needed to understand the functional mechanisms by which these SNPs impact disease risk.

Mapping of hepatic expression quantitative trait loci (eQTLs) in a Han Chinese population

BackgroundElucidating the genetic basis underlying hepatic gene expression variability is of importance to understand the aetiology of the disease and variation in drug metabolism. To date, no genome-wide expression quantitative...

Expression Quantitative Trait Loci Analysis of Carotenoid Metabolism-related Genes in Citrus

The expression profiles of carotenoid metabolism pathway genes are one of the important factors that regulate carotenoid content and composition in the fruit of Citrus cultivars. To identify regulatory factors of gene expressions in the carotenoid pathway, expression quantitative trait loci (eQTL) analysis was applied using a hybrid population. In the analyzed population, significant eQTLs for phytoene synthase (PSY), phytoene desaturase (PDS), ζ-carotene desaturase (ZDS), β-ring hydroxylase (HYb), and zeaxanthin epoxidase (ZEP) were detected. Among these eQTLs, eQTLs for HYb and ZEP were located on their responsible gene loci, respectively. This result indicates that the expression level of HYb and ZEP was influenced by cis-elements in their up-stream regions. Conversely, eQTLs for PDS and ZDS were located in different loci from those of the responsible genes, indicating that the loci were trans-regulating factors. It also suggested that expression levels of PDS and ZDS were regulated by a common transcription factor because their eQTLs were co-localized. Significant eQTL of lycopene β-cyclase (LCYb) and 9-cis-epoxycarotenoid dioxygenase (NCED) were not detected in this population. The major factors to control gene expression depend on each carotenoid metabolism gene, and their combination might cause complex carotenoid metabolism and extend the variation of content and composition in citrus varieties. eQTL analysis was demonstrated as a powerful tool to identify cis- or trans-regulation for carotenoid metabolism genes, as well as in generating functionally defined markers for marker-assisted selection to increase the carotenoid content of citrus in breeding programs.

Vitamin D-responsive SGPP2 variants associated with lung cell expression and lung function

BackgroundVitamin D is associated with lung health in epidemiologic studies, but mechanisms mediating observed associations are poorly understood. This study explores mechanisms for an effect of vitamin D in lung through an in vivo gene expression study, an expression quantitative trait loci (eQTL) analysis in lung tissue, and a population-based cohort study of sequence variants.MethodsMicroarray analysis investigated the association of gene expression in small airway epithelial cells with serum 25(OH)D in adult non-smokers. Sequence variants in candidate genes identified by the microarray were investigated in a lung tissue eQTL database, and also in relation to cross-sectional pulmonary function in the Health, Aging, and Body Composition (Health ABC) study, stratified by race, with replication in the Framingham Heart Study (FHS).Results13 candidate genes had significant differences in expression by serum 25(OH)D (nominal p < 0.05), and a genome-wide significant eQTL association was detected for SGPP2. In Health ABC, SGPP2 SNPs were associated with FEV1 in both European- and African-Americans, and the gene-level association was replicated in European-American FHS participants. SNPs in 5 additional candidate genes (DAPK1, FSTL1, KAL1, KCNS3, and RSAD2) were associated with FEV1 in Health ABC participants.ConclusionsSGPP2, a sphingosine-1-phosphate phosphatase, is a novel vitamin D-responsive gene associated with lung function. The identified associations will need to be followed up in further studies.

Expression quantitative trait locus mapping across water availability environments reveals contrasting associations with genomic features in Arabidopsis.

The regulation of gene expression is crucial for an organism's development and response to stress, and an understanding of the evolution of gene expression is of fundamental importance to basic and applied biology. To improve this understanding, we conducted expression quantitative trait locus (eQTL) mapping in the Tsu-1 (Tsushima, Japan) × Kas-1 (Kashmir, India) recombinant inbred line population of Arabidopsis thaliana across soil drying treatments. We then used genome resequencing data to evaluate whether genomic features (promoter polymorphism, recombination rate, gene length, and gene density) are associated with genes responding to the environment (E) or with genes with genetic variation (G) in gene expression in the form of eQTLs. We identified thousands of genes that responded to soil drying and hundreds of main-effect eQTLs. However, we identified very few statistically significant eQTLs that interacted with the soil drying treatment (GxE eQTL). Analysis of genome resequencing data revealed associations of several genomic features with G and E genes. In general, E genes had lower promoter diversity and local recombination rates. By contrast, genes with eQTLs (G) had significantly greater promoter diversity and were located in genomic regions with higher recombination. These results suggest that genomic architecture may play an important a role in the evolution of gene expression.

Refining Susceptibility Loci of Chronic Obstructive Pulmonary Disease with Lung eqtls

Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of mortality worldwide. Recent genome-wide association studies (GWAS) have identified robust susceptibility loci associated with COPD. However, the mechanisms mediating the risk conferred by these loci remain to be found. The goal of this study was to identify causal genes/variants within susceptibility loci associated with COPD. In the discovery cohort, genome-wide gene expression profiles of 500 non-tumor lung specimens were obtained from patients undergoing lung surgery. Blood-DNA from the same patients were genotyped for 1,2 million SNPs. Following genotyping and gene expression quality control filters, 409 samples were analyzed. Lung expression quantitative trait loci (eQTLs) were identified and overlaid onto three COPD susceptibility loci derived from GWAS; 4q31 (HHIP), 4q22 (FAM13A), and 19q13 (RAB4B, EGLN2, MIA, CYP2A6). Significant eQTLs were replicated in two independent datasets (n = 363 and 339). SNPs previously associated with COPD and lung function on 4q31 (rs1828591, rs13118928) were associated with the mRNA expression of HHIP. An association between mRNA expression level of FAM13A and SNP rs2045517 was detected at 4q22, but did not reach statistical significance. At 19q13, significant eQTLs were detected with EGLN2. In summary, this study supports HHIP, FAM13A, and EGLN2 as the most likely causal COPD genes on 4q31, 4q22, and 19q13, respectively. Strong lung eQTL SNPs identified in this study will need to be tested for association with COPD in case-control studies. Further functional studies will also be needed to understand the role of genes regulated by disease-related variants in COPD.

THU0002 Fine mapping and expression studies of the 12Q13-14 locus associated with rheumatoid arthritis

BackgroundWith evidence initially from The Wellcome Trust Case Control Consortium genome wide association study (GWAS), a region on 12q13 has been convincingly and reproducibly associated with rheumatoid arthritis (RA) and...

Abstract 3626: Inflammation-related genetic variants predict radiation-induced toxicity following definitive radiotherapy for non-small cell lung cancer.

Abstract Definitive radiotherapy combined with chemotherapy is associated with enhanced locoregional control and clinical outcomes in non-small cell lung cancer (NSCLC) patients who are not candidates for surgery. However, radiation-induced pneumonitis and esophagitis are common dose-limiting inflammatory conditions that occur in these patients. These events can compromise the effect of the overall treatment and reduce quality of life in these patients. It is not clear what factors contribute to differences in the development of toxicities among patients with same clinical or dosimetric parameters. In this study, we conducted a large scale pathway-based analysis to identify inflammation-related genetic variants associated with radiation-induced pneumonitis or esophagitis. In the initial discovery phase, 11,930 inflammation-related genetic variants were genotyped in 243 stage I-III NSCLC patients who were treated with definitive radiotherapy at MD Anderson Cancer Center. A fast-track validation was performed for the top candidate SNPs in an additional 212 NSCLC cases following the same criteria. A total of 14 SNPs and 12 SNPs were validated for pneumonitis and esophagitis, respectively. Meta-analysis was used to summarize the results from two phases. The most significant SNPs for esophagitis and pneumonitis mapped to chromosome 18p (OR=9.05, 95%CI=2.68-30.59, P-meta=3.90×10−4) and 14q (OR=4.75, 95%CI=2.17-10.41, P-meta =9.61×10−5), respectively. To explore the potential underlying mechanisms, we analyzed radiation sensitivity and eQTL relationships of the identified SNPs and their host genes in 277 ethnically defined human lymphoblastoid cell lines. Among the 24 host genes for the 26 identified SNPs, three genes (DDX58, GZMB, CDC2) showed significant association with radiation responses in these cell lines (P&amp;lt;0.05). Moreover, a significant eQTL relationship was identified between CDC2 expression and the genotype of the identified CDC2 SNP (P=0.022). Interestingly, this CDC2 SNP is predicted to create a new miRNA binding site. Luciferase reporter assays to test the influence of this SNP on miRNA regulation is currently ongoing. Our results provide evidence that inflammation-related genetic variants contribute to the development of radiation-induced pneumonitis and esophagitis. The results also suggest that gene expression changes through altered miRNA binding may be an underlying mechanism for these events. Citation Format: Xia Pu, Liewei Wang, Michelle AT Hildebrandt, Yuanqing Ye, Joe Y. Chang, Charles Lu, Heath D. Skinner, Nifang Niu, Gregory D. Jenkins, Ritsuko Komaki, John D. Minna, Jack A. Roth, Richard M. Weinshilboum, Xifeng Wu. Inflammation-related genetic variants predict radiation-induced toxicity following definitive radiotherapy for non-small cell lung cancer. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3626. doi:10.1158/1538-7445.AM2013-3626