Genetic Variations in Key MicroRNA Processing Genes and Risk of Head and Neck Cancer: A Case-Control Study in Chinese Population

MicroRNAs (miRNAs) have been reported to play a key role in oncogenesis. Genetic variations in miRNA processing genes and miRNA binding sites may affect the biogenesis of miRNA and the miRNA-mRNA interactions, hence promoting tumorigenesis. In the present study, we hypothesized that potentially functional polymorphisms in miRNA processing genes may contribute to head and neck cancer (HNC) susceptibility. To test this hypothesis, we genotyped three SNPs at miRNA binding sites of miRNA processing genes (rs1057035 in 3′UTR of DICER, rs3803012 in 3′UTR of RAN and rs10773771 in 3′UTR of HIWI) with a case-control study including 397 HNC cases and 900 controls matched by age and sex in Chinese. Although none of three SNPs was significantly associated with overall risk of HNC, rs1057035 in 3′UTR of DICER was associated with a significantly decreased risk of oral cancer (TC/CC vs. TT: adjusted OR  = 0.65, 95% CI  = 0.46–0.92). Furthermore, luciferase activity assay showed that rs1057035 variant C allele led to significantly lower expression levels as compared to the T allele, which may be due to the relatively high inhibition of hsa-miR-574-3p on DICER mRNA. These findings indicated that rs1057035 located at 3′UTR of DICER may contribute to the risk of oral cancer by affecting the binding of miRNAs to DICER. Large-scale and well-designed studies are warranted to validate our findings.     

MicroRNA Related Polymorphisms and Breast Cancer Risk

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88–0.96), rs1052532 (OR 0.97; 95% CI: 0.95–0.99), rs10719 (OR 0.97; 95% CI: 0.94–0.99), rs4687554 (OR 0.97; 95% CI: 0.95–0.99, and rs3134615 (OR 1.03; 95% CI: 1.01–1.05) located in the 3′ UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.     

Rare BRCA1 haplotypes including 3′UTR SNPs associated with breast cancer risk

Genetic markers identifying women at an increased risk of developing breast cancer exist, yet the majority of inherited risk remains elusive. While numerous BRCA1 coding sequence mutations are associated with breast cancer risk, BRCA1 mutations account for less then 5% of breast cancer risk. Since 3′ untranslated region (3′UTR) polymorphisms disrupting microRNA (miRNA) binding can be functional and can act as genetic markers of cancer risk, we tested the hypothesis that such polymorphisms in the 3′UTR of BRCA1 and haplotypes containing these functional polymorphisms may be associated with breast cancer risk. We sequenced the BRCA1 3′UTR from breast cancer patients to identify miRNA disrupting polymorphisms. We further evaluated haplotypes of this region including the identified 3′UTR variants in a large population of controls and breast cancer patients (n = 221) with known breast cancer subtypes and ethnicities. We identified three 3′UTR variants in BRCA1 that are polymorphic in breast cancer populations, and haplotype analysis including these variants revealed that breast cancer patients harbor five rare haplotypes not generally found among controls (9.50% for breast cancer chromosomes, 0.11% for control chromosomes, p = 0.0001). Three of these rare haplotypes contain the rs8176318 BRCA1 3′UTR functional variant. These haplotypes are not biomarkers for BRCA1 coding mutations, as they are found rarely in BRCA1 mutant breast cancer patients (1/129 patients = 0.78%). These rare BRCA1 haplotypes and 3′UTR SNPs may represent new genetic markers of breast cancer risk.Key words: BRCA1, haplotype, microRNA, SNP, 3′UTR, breast cancer, triple negative breast cancer     

Association of Impulsivity and Polymorphic MicroRNA-641 Target Sites in the SNAP-25 Gene

Impulsivity is a personality trait of high impact and is connected with several types of maladaptive behavior and psychiatric diseases, such as attention deficit hyperactivity disorder, alcohol and drug abuse, as well as pathological gambling and mood disorders. Polymorphic variants of the SNAP-25 gene emerged as putative genetic components of impulsivity, as SNAP-25 protein plays an important role in the central nervous system, and its SNPs are associated with several psychiatric disorders. In this study we aimed to investigate if polymorphisms in the regulatory regions of the SNAP-25 gene are in association with normal variability of impulsivity. Genotypes and haplotypes of two polymorphisms in the promoter (rs6077690 and rs6039769) and two SNPs in the 3′ UTR (rs3746544 and rs1051312) of the SNAP-25 gene were determined in a healthy Hungarian population (N = 901) using PCR–RFLP or real-time PCR in combination with sequence specific probes. Significant association was found between the T–T 3′ UTR haplotype and impulsivity, whereas no association could be detected with genotypes or haplotypes of the promoter loci. According to sequence alignment, the polymorphisms in the 3′ UTR of the gene alter the binding site of microRNA-641, which was analyzed by luciferase reporter system. It was observed that haplotypes altering one or two nucleotides in the binding site of the seed region of microRNA-641 significantly increased the amount of generated protein in vitro. These findings support the role of polymorphic SNAP-25 variants both at psychogenetic and molecular biological levels.     

A Study on Genetic Variants of Fibroblast Growth Factor Receptor 2 (FGFR2) and the Risk of Breast Cancer from North India

Genome-Wide Association Studies (GWAS) have identified Fibroblast growth factor receptor 2 (FGFR2) as a candidate gene for breast cancer with single nucleotide polymorphisms (SNPs) located in intron 2 region as the susceptibility loci strongly associated with the risk. However, replicate studies have often failed to extrapolate the association to diverse ethnic regions. This hints towards the existing heterogeneity among different populations, arising due to differential linkage disequilibrium (LD) structures and frequencies of SNPs within the associated regions of the genome. It is therefore important to revisit the previously linked candidates in varied population groups to unravel the extent of heterogeneity. In an attempt to investigate the role of FGFR2 polymorphisms in susceptibility to the risk of breast cancer among North Indian women, we genotyped rs2981582, rs1219648, rs2981578 and rs7895676 polymorphisms in 368 breast cancer patients and 484 healthy controls by Polymerase chain reaction-Restriction fragment length polymorphism (PCR-RFLP) assay. We observed a statistically significant association with breast cancer risk for all the four genetic variants (P<0.05). In per-allele model for rs2981582, rs1219648, rs7895676 and in dominant model for rs2981578, association remained significant after bonferroni correction (P<0.0125). On performing stratified analysis, significant correlations with various clinicopathological as well as environmental and lifestyle characteristics were observed. It was evident that rs1219648 and rs2981578 interacted with exogenous hormone use and advanced clinical stage III (after Bonferroni correction, P<0.000694), respectively. Furthermore, combined analysis on these four loci revealed that compared to women with 0–1 risk loci, those with 2–4 risk loci had increased risk (OR = 1.645, 95%CI = 1.152–2.347, P = 0.006). In haplotype analysis, for rs2981578, rs2981582 and rs1219648, risk haplotype (GTG) was associated with a significantly increased risk compared to the common (ACA) haplotype (OR = 1.365, 95% CI = 1.086–1.717, P = 0.008). Our results suggest that intron 2 SNPs of FGFR2 may contribute to genetic susceptibility of breast cancer in North India population.     

Genetic Variation in DROSHA 3’UTR Regulated by hsa-miR-27b Is Associated with Bladder Cancer Risk

Purpose

miRNAs can regulate the biological processes, including differentiation, proliferation and apoptosis. DICER and DROSHA are two members of RNase III family, playing pivotal roles in the pathway of miRNAs biogenesis. In this study, we hypothesized that genetic variations of the DICER and DROSHA genes were associated with the bladder cancer risk.

Experimental Design

We performed a case-control study of 685 bladder cancer cases and 730 controls to investigate the association between the seven functional SNPs of DICER and DROSHA genes and bladder cancer risk. We then evaluated the functionality of the important SNPs.

Results

We found that rs10719T>C polymorphism located in 3’ untranslated region (UTR) of DROSHA gene was associated with the increased risk of bladder cancer. Stratified analysis suggested that rs10719TC/CC genotype can increase risk of bladder cancer among male patients (Adjusted OR = 1.34, 95% CI = 1.05-1.70, P = 0.018), and ever smokers (1.56, 1.14-2.14, 0.006), compared with TT genotype. Furthermore, DROSHA rs10719T>C polymorphism was predicted to regulate the binding activity of hsa-miR-27a/b. Luciferase reported gene assay confirmed that rs10719 T to G substitution disrupted the binding site for hsa-miR-27b, resulting the increased levels of DROSHA protein.

Conclusions

Taken together, these findings suggested that DROSHA rs10719T>C polymorphism may be associated with bladder cancer risk in a Chinese population, and hsa-miR-27b can influence the expression of DROSHA protein by binding with 3’UTR.     

Type II Transmembrane Serine Protease Gene Variants Associate with Breast Cancer

Type II transmembrane serine proteases (TTSPs) are related to tumor growth, invasion, and metastasis in cancer. Genetic variants in these genes may alter their function, leading to cancer onset and progression, and affect patient outcome. Here, 464 breast cancer cases and 370 controls were genotyped for 82 single-nucleotide polymorphisms covering eight genes. Association of the genotypes was estimated against breast cancer risk, breast cancer–specific survival, and survival in different treatment groups, and clinicopathological variables. SNPs in TMPRSS3 (rs3814903 and rs11203200), TMPRSS7 (rs1844925), and HGF (rs5745752) associated significantly with breast cancer risk (P _trend = 0.008–0.042). SNPs in TMPRSS1 (rs12151195 and rs12461158), TMPRSS2 (rs2276205), TMPRSS3 (rs3814903), and TMPRSS7 (rs2399403) associated with prognosis (P = 0.004–0.046). When estimating the combined effect of the variants, the risk of breast cancer was higher with 4–5 alleles present compared to 0–2 alleles (P = 0.0001; OR, 2.34; 95% CI, 1.39–3.94). Women with 6–8 survival-associating alleles had a 3.3 times higher risk of dying of breast cancer compared to women with 1–3 alleles (P = 0.001; HR, 3.30; 95% CI, 1.58–6.88). The results demonstrate the combined effect of variants in TTSPs and their related genes in breast cancer risk and patient outcome. Functional analysis of these variants will lead to further understanding of this gene family, which may improve individualized risk estimation and development of new strategies for treatment of breast cancer.     

Genetic Variations in the Flanking Regions of miR-101-2 Are Associated with Increased Risk of Breast Cancer

Genetic variants in human microRNA (miRNA) genes may alter mature miRNA processing and/or target selection, and likely contribute to cancer susceptibility and disease progression. Previous studies have suggested that miR-101 may play important roles in the development of cancer by regulating key tumor-associated genes. However, the role of single nucleotide polymorphisms (SNPs) of miR-101 in breast cancer susceptibility remains unclear. In this study, we genotyped 11 SNPs of the miR-101 genes (including miR-101-1 and miR-101-2) in a case-control study of 1064 breast cancer cases and 1073 cancer-free controls. The results revealed that rs462480 and rs1053872 in the flank regions of pre-miR-101-2 were significantly associated with increased risk of breast cancer (rs462480 AC/CC vs AA: adjusted OR = 1.182, 95% CI: 1.030–1.357, P = 0.017; rs1053872 CG/GG vs CC: adjusted OR = 1.179, 95% CI: 1.040–1.337, P = 0.010). However, the remaining 9 SNPs were not significantly associated with risk of breast cancer. Additionally, combined analysis of the two high-risk SNPs revealed that subjects carrying the variant genotypes of rs462480 and rs1053872 had increased risk of breast cancer in a dose-response manner (P _trend = 0.002). Compared with individuals with “0–1” risk allele, those carrying “2–4” risk alleles had 1.29-fold risk of breast cancer. In conclusion, these findings suggested that the SNPs rs462480 and rs1053872 residing in miR-101-2 gene may have a solid impact on genetic susceptibility to breast cancer, which may improve our understanding of the potential contribution of miRNA SNPs to cancer pathogenesis.     

Polymorphisms in a Putative Enhancer at the 10q21.2 Breast Cancer Risk Locus Regulate NRBF2 Expression

Genome-wide association studies have identified SNPs near ZNF365 at 10q21.2 that are associated with both breast cancer risk and mammographic density. To identify the most likely causal SNPs, we fine mapped the association signal by genotyping 428 SNPs across the region in 89,050 European and 12,893 Asian case and control subjects from the Breast Cancer Association Consortium. We identified four independent sets of correlated, highly trait-associated variants (iCHAVs), three of which were located within ZNF365. The most strongly risk-associated SNP, rs10995201 in iCHAV1, showed clear evidence of association with both estrogen receptor (ER)-positive (OR = 0.85 [0.82–0.88]) and ER-negative (OR = 0.87 [0.82–0.91]) disease, and was also the SNP most strongly associated with percent mammographic density. iCHAV2 (lead SNP, chr10: 64,258,684:D) and iCHAV3 (lead SNP, rs7922449) were also associated with ER-positive (OR = 0.93 [0.91–0.95] and OR = 1.06 [1.03–1.09]) and ER-negative (OR = 0.95 [0.91–0.98] and OR = 1.08 [1.04–1.13]) disease. There was weaker evidence for iCHAV4, located 5′ of ADO, associated only with ER-positive breast cancer (OR = 0.93 [0.90–0.96]). We found 12, 17, 18, and 2 candidate causal SNPs for breast cancer in iCHAVs 1–4, respectively. Chromosome conformation capture analysis showed that iCHAV2 interacts with the ZNF365 and NRBF2 (more than 600 kb away) promoters in normal and cancerous breast epithelial cells. Luciferase assays did not identify SNPs that affect transactivation of ZNF365, but identified a protective haplotype in iCHAV2, associated with silencing of the NRBF2 promoter, implicating this gene in the etiology of breast cancer.     

Association Analysis of GABRB3 Promoter Variants with Heroin Dependence

GABRB3 encoding the β3 subunit of GABAA receptor has been implicated in multiple neuropsychiatric disorders, including substance abuse. Previous studies reported that SNPs at the 5′ regulatory region of GABRB3 could regulate GABRB3 gene expression and associated with childhood absence epilepsy (CAE). The study aimed to investigate whether SNPs at the 5′ regulatory region of GABRB3 were associated with heroin dependence in our population. We first re-sequenced 1.5 kb of the 5′regulatory region of GABRB3 gene to examine the SNP profile in the genomic DNA of 365 control subjects. Then, we conducted a case-control association analysis between 576 subjects with heroin dependence (549 males, 27 females) and 886 controls (472 males, 414 females) by genotyping the rs4906902 as a tag SNP. We also conducted a reporter gene assay to assess the promoter activity of two major haplotypes derived from SNPs at this region. We detected 3 common SNPs (rs4906902, rs8179184 and rs20317) at this region that had strong pair-wise linkage disequilibrium. The C allele of rs4906902 was found to be associated with increased risk of heroin dependence (odds ratio:1.27, p = 0.002). Two major haplotypes (C-A-G and T-G-C) derived from these 3 SNPs accounted for 99% of this sample, and reporter gene activity assay showed that haplotype C-A-G that contained the C allele of the tag SNP rs4906902 had higher activity than haplotype T-G-C. Our data suggest that GABRB3 might be associated with heroin dependence, and increased expression of GABRB3 might contribute to the pathogenesis of heroin dependence.     

Fine-Mapping Angiotensin-Converting Enzyme Gene: Separate QTLs Identified for Hypertension and for ACE Activity

Angiotensin-converting enzyme (ACE) has been implicated in multiple biological system, particularly cardiovascular diseases. However, findings associating ACE insertion/deletion polymorphism with hypertension or other related traits are inconsistent. Therefore, in a two-stage approach, we aimed to fine-map ACE in order to narrow-down the function-specific locations. We genotyped 31 single nucleotide polymorphisms (SNPs) of ACE from 1168 individuals from 305 young-onset (age ≤40) hypertension pedigrees, and found four linkage disequilibrium (LD) blocks. A tag-SNP, rs1800764 on LD block 2, upstream of and near the ACE promoter, was significantly associated with young-onset hypertension (p = 0.04). Tag-SNPs on all LD blocks were significantly associated with ACE activity (p-value: 10^–16 to <10^–33). The two regions most associated with ACE activity were found between exon13 and intron18 and between intron 20 and 3′UTR, as revealed by measured haplotype analysis. These two major QTLs of ACE activity and the moderate effect variant upstream of ACE promoter for young-onset hypertension were replicated by another independent association study with 842 subjects.     

Association analysis of miRNA-related genetic polymorphisms in miR-143/145 and KRAS with colorectal cancer susceptibility and survival

Background: There is accumulating evidence of aberrant expression of miR-143 and miR-145 and their target gene KRAS in colorectal cancer (CRC). We hypothesize that single nucleotide polymorphisms (SNPs) within or near mRNA–microRNA (miRNA) binding sites may affect miRNA/target gene interaction, resulting in differential mRNA/protein expression and promoting the development and progression of CRC. Methods: We conducted a case–control study of 507 patients with CRC recruited from a tertiary hospital and 497 population-based controls to assess the association of genetic polymorphisms in miR-143/145 and the KRAS 3′ untranslated region (3′UTR) with susceptibility to CRC and patients’ survival. In addition, genetic variations of genomic regions located from 500 bp upstream to 500 bp downstream of the miR-143/miR-145 gene and the 3′UTR of KRAS were selected for analysis using the Haploview and HaploReg software. Results: Using publicly available expression profiling data, we found that miR-143/145 and KRAS expression were all reduced in rectal cancer tissue compared with adjacent non-neoplastic large intestinal mucosa. The rs74693964 C/T variant located 65 bp downstream of miR-145 genomic regions was observed to be associated with susceptibility to CRC (adjusted odds ratio (OR): 2.414, 95% CI: 1.385–4.206). Cumulative effects of miR-143 and miR-145 on CRC risk were observed (P_trend=0.03). Patients having CRC carrying variant genotype TT of KRAS rs712 had poorer survival (log-rank P=0.044, adjusted hazard ratio (HR): 4.328, 95% CI: 1.236–15.147). Conclusions: Our results indicate that miRNA-related polymorphisms in miR-143/145 and KRAS are likely to be deleterious and represent potential biomarkers for susceptibility to CRC and patients’ survival.     

A Newfound Association between MDC1 Functional Polymorphism and Lung Cancer Risk in Chinese

Mediator of DNA damage checkpoint protein 1 (MDC1) plays an early and core role in Double-Strand Break Repair (DDR) and ataxia telangiectasia-mutated (ATM) mediated response to DNA double-strand breaks (DSBs), and thus involves the pathogenesis of several DNA damage-related diseases such as cancer. We hypothesized that the single nucleotide polymorphisms (SNPs) of MDC1 which have potencies on affecting MDC1 expression or function were associated with risk of lung cancer. In a two-stage case-control study, we tested the association between 5 putatively functional SNPs of MDC1 and lung cancer risk in a southern Chinese population, and validated the promising association in an eastern Chinese population. We found the SNP rs4713354A>C that is located in the 5′-untranslated region of MDC1 was significantly associated with lung cancer risk in both populations (P = 0.024), with an odds ratio as 1.23(95% confidence interval  = 1.35–1.26) for the rs4713354C (CA+CC) genotypes compared to the rs4713354AA genotype. However, no significant association was observed between other SNPs and lung cancer risk. The gene-based analysis rested with these SNPs suggested the MDC1 as a susceptible gene for lung cancer (P = 0.009). Moreover, by querying the gene expression database, we further found that the rs4713354C genotypes confer a significantly lower mRNA expression of MDC1 than the rs4713354AA genotype in 260 cases of lymphoblastoid cells (P = 0.002). Our data suggested that the SNP rs4713354A>C of MDC1 may be a functional genetic biomarker for susceptibility to lung cancer in Chinese.     

Integrative Approach Detected Association between Genetic Variants of microRNA Binding Sites of TLRs Pathway Genes and OSCC Susceptibility in Chinese Han Population

Oral squamous cell carcinoma (OSCC) is a leading malignancy worldwide; the overall 5-year survival rate is approximately 50%. A variety of proteins in Toll-like receptors (TLRs) pathway have been related with the risk of OSCC. However, the influence of genetic variations in TLRs pathway genes on OSCC susceptibility is unclear. Previous studies mainly focused on the coding region of genes, while the UTR region remains unstudied. In the current study, a bioinformatics approach was performed to select candidate single nucleotide polymorphisms (SNPs) on microRNA binding sites of TLRs pathway genes related with OSCC. After screening 90 OSCC related TLRs pathway genes, 16 SNPs were selected for genotyping. We found that rs5030486, the polymorphisms on 3′ UTR of TRAF6, was significantly associated with OSCC risk. AG genotype of TRAF6 was strongly associated with a decreased risk of OSCC (OR = 0.252; 95% CI = 0.106, 0.598; p = 0.001). In addition, AG genotype was also related with a reduced risk of OSCC progression both in univariable analysis (HR = 0.303, 95% CI = 0.092, 0.995) and multivariable analysis (HR = 0.272, 95% CI = 0.082, 0.903). Furthermore, after detecting the mRNA expression level of TRAF6 in 24 OSCC patients, we found that TRAF6 expression level was significantly different between patients carrying different genotypes at locus rs5030486 (p = 0.013), indicating that rs5030486 of TRAF6 might contribute to OSCC risk by altering TRAF6 expression level. In general, these data indicated that SNP rs5030486 could be a potential bio-marker for OSCC risk and our results might provide new insights into the association of polymorphisms within the non-coding area of genes with cancers.     

Polymorphisms in ERCC1 and XPF Genes and Risk of Gastric Cancer in an Eastern Chinese Population

Background

Inherited functional single nucleotide polymorphisms (SNPs) in DNA repair genes may alter DNA repair capacity and thus contribute to cancer risk.

Methods

Three ERCC1 functional SNPs (rs2298881C>A, rs3212986C>A and rs11615G>A) and two XPF/ERCC4 functional SNPs (rs2276466C>G and rs6498486A>C) were genotyped for 1125 gastric adenocarcinoma cases and 1196 cancer-free controls by Taqman assays. Odds ratios (OR) and 95% confidence intervals (CI) were used to estimate risk associations, and false-positive report probabilities (FPRP) were calculated for assessing significant findings.

Results

ERCC1 rs2298881C and rs11615A variant genotypes were associated with increased gastric cancer risk (adjusted OR = 1.33, 95% CI = 1.05–1.67 for rs2298881 AC/CC and adjusted OR = 1.23, 95% CI = 1.05–1.46 for rs11615 AG/AA, compared with their common genotype AA and GG, respectively). Patients with 2–3 ERCC1 risk genotypes had significant increased risk (adjusted OR = 1.56, 95% CI = 1.27–1.93), compared with those with 0–1 ERCC1 risk genotypes, and this risk was more significantly in subgroups of never drinkers, non-gastric cardia adenocarcinoma (NGCA) and clinical stage I+II. All these risks were not observed for XPF SNPs.

Conclusions

These findings suggest that functional ERCC1 SNPs may contribute to risk of gastric cancer. Larger and well-designed studies with different ethnic populations are needed to validate our findings.     

Genetic Variants in Inflammation-Related Genes Are Associated with Radiation-Induced Toxicity Following Treatment for Non-Small Cell Lung Cancer

Treatment of non-small cell lung cancer (NSCLC) with radiotherapy or chemoradiotherapy is often accompanied by the development of esophagitis and pneumonitis. Identifying patients who might be at increased risk for normal tissue toxicity would help in determination of the optimal radiation dose to avoid these events. We profiled 59 single nucleotide polymorphisms (SNPs) from 37 inflammation-related genes in 173 NSCLC patients with stage IIIA/IIIB (dry) disease who were treated with definitive radiation or chemoradiation. For esophagitis risk, nine SNPs were associated with a 1.5- to 4-fold increase in risk, including three PTGS2 (COX2) variants: rs20417 (HR:1.93, 95% CI:1.10–3.39), rs5275 (HR:1.58, 95% CI:1.09–2.27), and rs689470 (HR:3.38, 95% CI:1.09–10.49). Significantly increased risk of pneumonitis was observed for patients with genetic variation in the proinflammatory genes IL1A, IL8, TNF, TNFRSF1B, and MIF. In contrast, NOS3:rs1799983 displayed a protective effect with a 45% reduction in pneumonitis risk (HR:0.55, 95% CI:0.31–0.96). Pneumonitis risk was also modulated by polymorphisms in anti-inflammatory genes, including genetic variation in IL13. rs20541 and rs180925 each resulted in increased risk (HR:2.95, 95% CI:1.14–7.63 and HR:3.23, 95% CI:1.03–10.18, respectively). The cumulative effect of these SNPs on risk was dose-dependent, as evidenced by a significantly increased risk of either toxicity with an increasing number of risk genotypes (P<0.001). These results suggest that genetic variations among inflammation pathway genes may modulate the development of radiation-induced toxicity and, ultimately, help in identifying patients who are at an increased likelihood for such events.