Genetic and microbial factors modulating the ubiquitin proteasome system in inflammatory bowel disease

OBJECTIVE: Altered microbiota composition, changes in immune responses and impaired intestinal barrier functions are observed in IBD. Most of these features are controlled by proteases and their inhibitors to maintain gut homeostasis. Unrestrained or excessive proteolysis can lead to pathological gastrointestinal conditions. The aim was to validate the identified protease IBD candidates from a previously performed systematic review through a genetic association study and functional follow-up. DESIGN: We performed a genetic association study in a large multicentre cohort of patients with Crohn's disease (CD) and UC from five European IBD referral centres in a total of 2320 CD patients, 2112 UC patients and 1796 healthy controls. Subsequently, we did an extensive functional assessment of the candidate genes to explore their causality in IBD pathogenesis. RESULTS: Ten single nucleotide polymorphisms (SNPs) in four genes were significantly associated with CD: CYLD, USP40, APEH and USP3. CYLD was the most significant gene with the intronically located rs12324931 the strongest associated SNP (pFDR=1.74e-17, OR=2.24 (1.83 to 2.74)). Five SNPs in four genes were significantly associated with UC: USP40, APEH, DAG1 and USP3. CYLD, as well as some of the other associated genes, is part of the ubiquitin proteasome system (UPS). We therefore determined if the IBD-associated adherent-invasive Escherichia coli (AIEC) can modulate the UPS functioning. Infection of intestinal epithelial cells with the AIEC LF82 reference strain modulated the UPS turnover by reducing poly-ubiquitin conjugate accumulation, increasing 26S proteasome activities and decreasing protein levels of the NF-κB regulator CYLD. This resulted in IκB-α degradation and NF-κB activation. This activity was very important for the pathogenicity of AIEC since decreased CYLD resulted in increased ability of AIEC LF82 to replicate intracellularly. CONCLUSIONS: Our results reveal the UPS, and CYLD specifically, as an important contributor to IBD pathogenesis, which is favoured by both genetic and microbial factors.

Genetic Analyses Reveal a Role for Vitamin D Insufficiency in HCV-Associated Hepatocellular Carcinoma Development

Background Vitamin D insufficiency has been associated with the occurrence of various types of cancer, but causal relationships remain elusive. We therefore aimed to determine the relationship between genetic determinants of vitamin D serum levels and the risk of developing hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC). Methodology/Principal Findings Associations between CYP2R1, GC, and DHCR7 genotypes that are determinants of reduced 25-hydroxyvitamin D (25[OH]D3) serum levels and the risk of HCV-related HCC development were investigated for 1279 chronic hepatitis C patients with HCC and 4325 without HCC, respectively. The well-known associations between CYP2R1 (rs1993116, rs10741657), GC (rs2282679), and DHCR7 (rs7944926, rs12785878) genotypes and 25(OH)D3 serum levels were also apparent in patients with chronic hepatitis C. The same genotypes of these single nucleotide polymorphisms (SNPs) that are associated with reduced 25(OH)D3 serum levels were found to be associated with HCV-related HCC (P = 0.07 [OR = 1.13, 95% CI = 0.99–1.28] for CYP2R1, P = 0.007 [OR = 1.56, 95% CI = 1.12–2.15] for GC, P = 0.003 [OR = 1.42, 95% CI = 1.13–1.78] for DHCR7; ORs for risk genotypes). In contrast, no association between these genetic variations and liver fibrosis progression rate (P>0.2 for each SNP) or outcome of standard therapy with pegylated interferon-α and ribavirin (P>0.2 for each SNP) was observed, suggesting a specific influence of the genetic determinants of 25(OH)D3 serum levels on hepatocarcinogenesis. Conclusions/Significance Our data suggest a relatively weak but functionally relevant role for vitamin D in the prevention of HCV-related hepatocarcinogenesis.

Bovine prion protein gene (PRNP) promoter polymorphisms modulate PRNP expression and may be responsible for differences in bovine spongiform encephalopathy susceptibility

The susceptibility of humans to the variant Creutzfeldt-Jakob disease is greatly influenced by polymorphisms within the human prion protein gene (PRNP). Similar genetic differences exist in sheep, in which PRNP polymorphisms modify the susceptibility to scrapie. However, the known coding polymorphisms within the bovine PRNP gene have little or no effect on bovine spongiform encephalopathy (BSE) susceptibility in cattle. We have recently found a tentative association between PRNP promoter polymorphisms and BSE susceptibility in German cattle (Sander, P., Hamann, H., Pfeiffer, I., Wemheuer, W., Brenig, B., Groschup, M., Ziegler, U., Distl, O., and Leeb, T. (2004) Neurogenetics 5, 19-25). A plausible hypothesis explaining this observation could be that the bovine PRNP promoter polymorphisms cause changes in PRNP expression that might be responsible for differences in BSE incubation time and/or BSE susceptibility. To test this hypothesis, we performed a functional promoter analysis of the different bovine PRNP promoter alleles by reporter gene assays in vitro and by measuring PRNP mRNA levels in calves with different PRNP genotypes in vivo. Two variable sites, a 23-bp insertion/deletion (indel) polymorphism containing a RP58-binding site and a 12-bp indel polymorphism containing an SP1-binding site, were investigated. Band shift assays indicated differences in transcription factor binding to the different alleles at the two polymorphisms. Reporter gene assays demonstrated an interaction between the two postulated transcription factors and lower expression levels of the ins/ins allele compared with the del/del allele. The in vivo data revealed substantial individual variation of PRNP expression in different tissues. In intestinal lymph nodes, expression levels differed between the different PRNP genotypes.

Polymorphisms within the canine MLPH gene are associated with dilute coat color in dogs

BACKGROUND Pinschers and other dogs with coat color dilution show a characteristic pigmentation phenotype. The fur colors are a lighter shade, e.g. silvery grey (blue) instead of black and a sandy color (Isabella fawn) instead of red or brown. In some dogs the coat color dilution is sometimes accompanied by hair loss and recurrent skin inflammation, the so called color dilution alopecia (CDA) or black hair follicular dysplasia (BHFD). In humans and mice a comparable pigmentation phenotype without any documented hair loss is caused by mutations within the melanophilin gene (MLPH). RESULTS We sequenced the canine MLPH gene and performed a mutation analysis of the MLPH exons in 6 Doberman Pinschers and 5 German Pinschers. A total of 48 sequence variations was identified within and between the breeds. Three families of dogs showed co-segregation for at least one polymorphism in an MLPH exon and the dilute phenotype. No single polymorphism was identified in the coding sequences or at splice sites that is likely to be causative for the dilute phenotype of all dogs examined. In 18 German Pinschers a mutation in exon 7 (R199H) was consistently associated with the dilute phenotype. However, as this mutation was present in homozygous state in four dogs of other breeds with wildtype pigmentation, it seems unlikely that this mutation is truly causative for coat color dilution. In Doberman Pinschers as well as in Large Munsterlanders with BHFD, a set of single nucleotide polymorphisms (SNPs) around exon 2 was identified that show a highly significant association to the dilute phenotype. CONCLUSION This study provides evidence that coat color dilution is caused by one or more mutations within or near the MLPH gene in several dog breeds. The data on polymorphisms that are strongly associated with the dilute phenotype will allow the genetic testing of Pinschers to facilitate the breeding of dogs with defined coat colors and to select against Large Munsterlanders carrying BHFD.

Genetic mechanisms of Na+, K+-ATPase regulation in hypertension

The spontaneously hypertensive rat (SHR) is a model of essential hypertension. During the early development of hypertension, the SHR demonstrates increased proximal tubule (PT) Na+ reabsorption. I hypothesized that the increased PT Na+ reabsorption exhibited by the young SHR was due to altered sub-cellular distribution of Na+, K +-ATPase compared to the normotensive Wistar Kyoto (WKY). The hypothesis is supported, herein, by observations of greater Na+, K +-ATPase α 1 abundance in PT plasma membrane and lower abundance in late endosomes of 4wk SHR despite no difference in total PT α 1 abundance. There is a greater amount of Ser-18 unphosphorylated α 1 in the 4wk SHR PT. Total PT Na+, K+-ATPase γ abundance is greater in SHR at 4wk and 16wk but γ abundance in plasma membrane is greater only at 4wk. The phosphatase, calcineurin, was chosen for study because it is involved in the stimulation of Na+, K +-ATPase. No difference in calcineurin coding sequence, expression, or activity was observed in SHR. Gene expression arrays were next used to find candidate genes involved in the regulation of Na+, K +-ATPase. The first candidate analyzed was soluble epoxide hydrolase (sEH). The gene encoding sEH (EPHX2) showed lower expression in SHR. There was also a reduction in sEH protein abundance but there was no correlation between protein abundance and blood pressure in F2 progeny. Two EPHX2 alleles were identified, an ancestral allele and a variant allele containing four polymorphisms. sEH activity was greater in animals carrying the variant allele but the inheritance of the variant allele did not correlate with blood pressure. Gene expression arrays also led to the examination of genes involved in redox balance/Na+, K+-ATPase regulation. A pattern of lower expression of genes involved in reactive radical detoxification in SHR was discerned. Six transcription factor binding sites were identified that occurred more often in these genes. Three transcription factors that bind to the HNF1 site were expressed at lower levels in SHR. This points to the HNF1 transcriptional complex as an important trans-acting regulator of a wide range of genes involved in altered redox balance in SHR. ^

Are functional and genetic components of platelets linked to coronary artery disease: A case-control study

Coronary artery disease (CAD) is a multifactorial disease process involving behavioral, inflammatory, clinical, thrombotic, and genetic components. Previous epidemiologic studies focused on identifying behavioral and demographic risk factors of CAD, but none focused on platelets. Current platelet literature lacks the known effects of platelet function and platelet receptor polymorphisms on CAD. This case-control analysis addressed these issues by analyzing data collected for a previous study. Cases were individuals who had undergone CABG and thus had been diagnosed with CAD, while the controls were volunteers presumed to be CAD free. The platelet function variables analyzed included fibrinogen Von Willebrand Factor activity (VWF), shear-induced platelet aggregation (SIPA), sCD40L, and mean platelet volume; and the platelet polymorphisms studied included PIA, α2 807, Ko, Kozak, and VNTR. Univariate analysis found fibrinogen, VWF, SIPA, and PIA to be independent risk factors of CAD. Logistic regression was used to build a predictive model for CAD using the platelet function and platelet polymorphism data adjusted for age, sex, race, and current smoking status. A model containing only platelet polymorphisms and their respective receptor densities, found polymorphisms within GPIbα to be associated with CAD, yielding an 86% (95% C.I. 0.97–3.55) increased risk with the presence of at least 1 polymorphism in Ko, Kozak, or VNTR. Another model included both platelet function and platelet polymorphism data. Fibrinogen, the receptor density of GPIbα, and the polymorphism in GPIa-IIa (α2 807) were all associated with CAD with odds ratios of 1.10, 1.04, and 2.30 for fibrinogen (10mg/dl increase), GPIbα receptors (1 MFI increase), and GPIa-IIa, respectively. In addition, risk estimates and 99% confidence intervals adjusted for race were calculated to determine if the presence of a platelet receptor polymorphism was associated with CAD. The results were as follows: PIA (1.64, 0.74–3.65); α2 807 (1.35, 0.77–2.37); Ko (1.71, 0.70–4.16); Kozak (1.17, 0.54–2.52); and VNTR (1.24, 0.52–2.91). Although not statistically significant, all platelet polymorphisms were associated with an increased risk for CAD. These exploratory findings indicate that platelets do appear to have a role in atherosclerosis and that anti-platelet drugs targeting GPI-IIa and GPIbα may be better treatment candidates for individuals with CAD. ^

Nonlinear tests for genetic studies of complex disease

Linkage and association studies are major analytical tools to search for susceptibility genes for complex diseases. With the availability of large collection of single nucleotide polymorphisms (SNPs) and the rapid progresses for high throughput genotyping technologies, together with the ambitious goals of the International HapMap Project, genetic markers covering the whole genome will be available for genome-wide linkage and association studies. In order not to inflate the type I error rate in performing genome-wide linkage and association studies, multiple adjustment for the significant level for each independent linkage and/or association test is required, and this has led to the suggestion of genome-wide significant cut-off as low as 5 × 10 −7. Almost no linkage and/or association study can meet such a stringent threshold by the standard statistical methods. Developing new statistics with high power is urgently needed to tackle this problem. This dissertation proposes and explores a class of novel test statistics that can be used in both population-based and family-based genetic data by employing a completely new strategy, which uses nonlinear transformation of the sample means to construct test statistics for linkage and association studies. Extensive simulation studies are used to illustrate the properties of the nonlinear test statistics. Power calculations are performed using both analytical and empirical methods. Finally, real data sets are analyzed with the nonlinear test statistics. Results show that the nonlinear test statistics have correct type I error rates, and most of the studied nonlinear test statistics have higher power than the standard chi-square test. This dissertation introduces a new idea to design novel test statistics with high power and might open new ways to mapping susceptibility genes for complex diseases. ^

Identification of genetic variation influencing apolipoprotein E levels: Follow-up of genome-wide linkage scans in the GENOA study

Apolipoprotein E (ApoE) plays a major role in the metabolism of high density and low density lipoproteins (HDL and LDL). Its common protein isoforms (E2, E3, E4) are risk factors for coronary artery disease (CAD) and explain between 16 to 23% of the inter-individual variation in plasma apoE levels. Linkage analysis has been completed for plasma apoE levels in the GENOA study (Genetic Epidemiology Network of Atherosclerosis). After stratification of the population by lipoprotein levels and body mass index (BMI) to create more homogeneity with regard to biological context for apoE levels, Hispanic families showed significant linkage on chromosome 17q for two strata (LOD=2.93 at 104 cM for a low cholesterol group, LOD=3.04 at 111 cM for a low cholesterol, high HDLC group). Replication of 17q linkage was observed for apoB and apoE levels in the unstratified Hispanic and African-American populations, and for apoE levels in African-American families. Replication of this 17q linkage in different populations and strata provides strong support for the presence of gene(s) in this region with significant roles in the determination of inter-individual variation in plasma apoE levels. Through a positional and functional candidate gene approach, ten genes were identified in the 17q linked region, and 62 polymorphisms in these genes were genotyped in the GENOA families. Association analysis was performed with FBAT, GEE, and variance-component based tests followed by conditional linkage analysis. Association studies with partial coverage of TagSNPs in the gene coding for apolipoprotein H (APOH) were performed, and significant results were found for 2 SNPs (APOH_20951 and APOH_05407) in the Hispanic low cholesterol strata accounting for 3.49% of the inter-individual variation in plasma apoE levels. Among the other candidate genes, we identified a haplotype block in the ACE1 gene that contains two major haplotypes associated with apoE levels as well as total cholesterol, apoB and LDLC levels in the unstratified Hispanic population. Identifying genes responsible for the remaining 60% of inter-individual variation in plasma apoE level, will yield new insights into the understanding of genetic interactions involved in the lipid metabolism, and a more precise understanding of the risk factors leading to CAD. ^

DNA nucleotide excision repair gene single nucleotide polymorphisms and hereditary nonpolyposis colorectal cancer

Hereditary nonpolyposis colorectal cancer (HNPCC) is an autosomal dominant disease caused by germline mutations in DNA mismatch repair(MMR) genes. The nucleotide excision repair(NER) pathway plays a very important role in cancer development. We systematically studied interactions between NER and MMR genes to identify NER gene single nucleotide polymorphism (SNP) risk factors that modify the effect of MMR mutations on risk for cancer in HNPCC. We analyzed data from polymorphisms in 10 NER genes that had been genotyped in HNPCC patients that carry MSH2 and MLH1 gene mutations. The influence of the NER gene SNPs on time to onset of colorectal cancer (CRC) was assessed using survival analysis and a semiparametric proportional hazard model. We found the median age of onset for CRC among MMR mutation carriers with the ERCC1 mutation was 3.9 years earlier than patients with wildtype ERCC1(median 47.7 vs 51.6, log-rank test p=0.035). The influence of Rad23B A249V SNP on age of onset of HNPCC is age dependent (likelihood ratio test p=0.0056). Interestingly, using the likelihood ratio test, we also found evidence of genetic interactions between the MMR gene mutations and SNPs in ERCC1 gene(C8092A) and XPG/ERCC5 gene(D1104H) with p-values of 0.004 and 0.042, respectively. An assessment using tree structured survival analysis (TSSA) showed distinct gene interactions in MLH1 mutation carriers and MSH2 mutation carriers. ERCC1 SNP genotypes greatly modified the age onset of HNPCC in MSH2 mutation carriers, while no effect was detected in MLH1 mutation carriers. Given the NER genes in this study play different roles in NER pathway, they may have distinct influences on the development of HNPCC. The findings of this study are very important for elucidation of the molecular mechanism of colon cancer development and for understanding why some mutation carriers of the MSH2 and MLH1 gene develop CRC early and others never develop CRC. Overall, the findings also have important implications for the development of early detection strategies and prevention as well as understanding the mechanism of colorectal carcinogenesis in HNPCC. ^

Characterization of glucocorticoid receptor (NR3C1) gene polymorphisms and a family-based association study with hypertension

Hypertension is usually defined as having values of systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg. Hypertension is one of the main adverse effects of glucocorticoid on the cardiovascular system. Glucocorticoids are essential hormones, secreted from adrenal glands in circadian fashion. Glucocorticoid's effect on blood pressure is conveyed by the glucocorticoid receptor (NR3C1), an omnipresent nuclear transcription factor. Although polymorphisms in this gene have long been implicated to be a causal factor for cardiovascular diseases such as hypertension, no study has yet thoroughly interrogated the gene's polymorphisms for their effect on blood pressure levels. Therefore, I have first resequenced ∼30 kb of the gene, encompassing all exons, promoter regions, 5'/3' UTRs as well as at least 1.5 kb of the gene's flanking regions from 114 chromosome 5 monosomic cell lines, comprised of three major American ethnic groups—European American, African American and Mexican American. I observed 115 polymorphisms and 14 common molecularly phased haplotypes. A subset of markers was chosen for genotyping study populations of GENOA (Genetic Epidemiology Network of Atherosclerosis; 1022 non-Hispanic whites, 1228 African Americans and 954 Mexican Americans). Since these study populations include sibships, the family-based association test was performed on 4 blood pressure-related quantitative variables—pulse, systolic blood pressure, diastolic blood pressure and mean arterial pressure. Using these analyses, multiple correlated SNPs are significantly protective against high systolic blood pressure in non-Hispanic whites, which includes rsb198, a SNP formerly associated with beneficial body compositions. Haplotype association analysis also supports this finding and all p-values remained significant after permutation tests. I therefore conclude that multiple correlated SNPs on the gene may confer protection against high blood pressure in non-Hispanic whites. ^

Immune and genetic risk factors in glioma

Despite extensive research, the etiology of adult glioma remains largely unknown. We sought to further explore the role of immune and genetic factors in glioma etiology using data from the Harris County Brain Tumor Study and the first U.S. genome-wide association study of glioma. First, using a case-control study design, we examined the association between adult glioma risk and surrogates of the timing and frequency of common early childhood infections, birth order and sibship size, respectively. We found that each one-unit increase in birth order was associated with a 12% decreased risk of glioma development in adulthood (OR=0.88, 95% CI=0.81-0.96); however, sibship size was not associated with adult glioma risk (OR=0.96, 95% CI=0.91-1.02). Second, we used a multi-strategic approach to explore the relationships between glioma risk, history of asthma/allergies, and 23 functional SNPs in 11 inflammation genes. We found three inflammation gene SNPs to be significantly associated with glioma risk (COX2/PTGS2 rs20417 [OR=1.41]; IL10 rs1800896 [OR=1.57]; and IL13 rs20541 [OR=0.39]). Joint effects analysis of the risk-conferring alleles of these three SNPs revealed a trend of increasing risk with increasing number of adverse alleles among those without asthma/allergies (p<0.0001). Finally, we conducted a case-only study to explore pairwise SNP-SNP interactions in immune-related pathways among a population of 1304 non-Hispanic white glioma cases. After correction for multiple comparisons, we found 279 significant SNP-SNP interactions among polymorphisms of immune-related genes, many of which have not been previously examined. Our results, cumulatively, suggest an important role for immune and genetic factors in glioma etiology and provide several new hypotheses for future studies.^

Single nucleotide polymorphisms (SNPs) associated with TGF-beta pathway and their significance in systemic sclerosis - A multilevel analysis

Systemic sclerosis (SSc) or Scleroderma is a complex disease and its etiopathogenesis remains unelucidated. Fibrosis in multiple organs is a key feature of SSc and studies have shown that transforming growth factor-β (TGF-β) pathway has a crucial role in fibrotic responses. For a complex disease such as SSc, expression quantitative trait loci (eQTL) analysis is a powerful tool for identifying genetic variations that affect expression of genes involved in this disease. In this study, a multilevel model is described to perform a multivariate eQTL for identifying genetic variation (SNPs) specifically associated with the expression of three members of TGF-β pathway, CTGF, SPARC and COL3A1. The uniqueness of this model is that all three genes were included in one model, rather than one gene being examined at a time. A protein might contribute to multiple pathways and this approach allows the identification of important genetic variations linked to multiple genes belonging to the same pathway. In this study, 29 SNPs were identified and 16 of them located in known genes. Exploring the roles of these genes in TGF-β regulation will help elucidate the etiology of SSc, which will in turn help to better manage this complex disease. ^

Genetic variants in the mammalian target of rapamycin (mTOR) signaling pathway as predictors of survival and clinical response in women with ovarian cancer

Background. The mTOR pathway is commonly altered in human tumors and promotes cell survival and proliferation. Preliminary evidence suggests this pathway's involvement in chemoresistance to platinum and taxanes, first line therapy for epithelial ovarian cancer. A pathway-based approach was used to identify individual germline single nucleotide polymorphisms (SNPs) and cumulative effects of multiple genetic variants in mTOR pathway genes and their association with clinical outcome in women with ovarian cancer. ^ Methods. The case-series was restricted to 319 non-Hispanic white women with high grade ovarian cancer treated with surgery and platinum-based chemotherapy. 135 SNPs in 20 representative genes in the mTOR pathway were genotyped. Hazard ratios (HRs) for death and Odds ratios (ORs) for failure to respond to primary therapy were estimated for each SNP using the multivariate Cox proportional hazards model and multivariate logistic regression model, respectively, while adjusting for age, stage, histology and treatment sequence. A survival tree analysis of SNPs with a statistically significant association (p<0.05) was performed to identify higher order gene-gene interactions and their association with overall survival. ^ Results. There was no statistically significant difference in survival by tumor histology or treatment regimen. The median survival for the cohort was 48.3 months. Seven SNPs were significantly associated with decreased survival. Compared to those with no unfavorable genotypes, the HR for death increased significantly with the increasing number of unfavorable genotypes and women in the highest risk category had HR of 4.06 (95% CI 2.29–7.21). The survival tree analysis also identified patients with different survival patterns based on their genetic profiles. 13 SNPs on five different genes were found to be significantly associated with a treatment response, defined as no evidence of disease after completion of primary therapy. Rare homozygous genotype of SNP rs6973428 showed a 5.5-fold increased risk compared to the wild type carrying genotypes. In the cumulative effect analysis, the highest risk group (individuals with ≥8 unfavorable genotypes) was significantly less likely to respond to chemotherapy (OR=8.40, 95% CI 3.10–22.75) compared to the low risk group (≤4 unfavorable genotypes). ^ Conclusions. A pathway-based approach can demonstrate cumulative effects of multiple genetic variants on clinical response to chemotherapy and survival. Therapy targeting the mTOR pathway may modify outcome in select patients.^

Molecular epidemiological studies on genetic predispositions to squamous cell carcinoma of the head and neck

Two molecular epidemiological studies were conducted to examine associations between genetic variation and risk of squamous cell carcinoma of the head and neck (SCCHN). In the first study, we hypothesized that genetic variation in p53 response elements (REs) may play roles in the etiology of SCCHN. We selected and genotyped five polymorphic p53 REs as well as a most frequently studied p53 codon 72 (Arg72Pro, rs1042522) polymorphism in 1,100 non-Hispanic White SCCHN patients and 1,122 age-and sex-matched cancer-free controls recruited at The University of Texas M. D. Anderson Cancer Center. In multivariate logistic regression analysis with adjustment for age, sex, smoking and drinking status, marital status and education level, we observed that the EOMES rs3806624 CC genotype had a significant effect of protection against SCCHN risk (adjusted odds ratio= 0.79, 95% confidence interval =0.64–0.98), compared with the -838TT+CT genotypes. Moreover, a significantly increased risk associated with the combined genotypes of p53 codon 72CC and EOMES -838TT+CT was observed, especially in the subgroup of non-oropharyneal cancer patients. The values of false-positive report probability were also calculated for significant findings. In the second study, we assessed the association between SCCHN risk and four potential regulatory single nucleotide polymorphisms (SNPs) of DEC1 (deleted in esophageal cancer 1) gene, a candidate tumor suppressor gene for esophageal cancer. After adjustment for age, sex, and smoking and drinking status, the variant -606CC (i.e., -249CC) homozygotes had a significantly reduced SCCHN risk (adjusted odds ratio = 0.71, 95% confidence interval = 0.52–0.99), compared with the -606TT homozygotes. Stratification analyses showed that a reduced risk associated with the -606CC genotype was more pronounced in subgroups of non-smokers, non-drinkers, younger subjects (defined as ≤ 57 years), carriers of TP53 Arg/Arg (rs1042522) genotype, patients with oropharyngeal cancer or late-stage SCCHN. Further in silico analysis revealed that the -249 T-to-C change led to a gain of a transcription factor binding site. Additional functional analysis showed that the -249T-to-C change significantly enhanced transcriptional activity of the DEC1 promoter and the DNA-protein binding activity. We conclude that the DEC1 promoter -249 T>C (rs2012775) polymorphism is functional, modulating susceptibility to SCCHN among non-Hispanic Whites. Additional large-scale, preferably population-based studies are needed to validate our findings.^

Genetic Predictors of Clinical Outcomes in Non-Small Cell Lung Cancer Patients

Lung cancer is the leading cause of cancer-related mortality in the US. Emerging evidence has shown that host genetic factors can interact with environmental exposures to influence patient susceptibility to the diseases as well as clinical outcomes, such as survival and recurrence. We aimed to identify genetic prognostic markers for non-small cell lung cancer (NSCLC), a major (85%) subtype of lung cancer, and also in other subgroups. With the fast evolution of genotyping technology, genetic association studies have went through candidate gene approach, to pathway-based approach, to the genome wide association study (GWAS). Even in the era of GWAS, pathway-based approach has its own advantages on studying cancer clinical outcomes: it is cost-effective, requiring a smaller sample size than GWAS easier to identify a validation population and explore gene-gene interactions. In the current study, we adopted pathway-based approach focusing on two critical pathways - miRNA and inflammation pathways. MicroRNAs (miRNA) post-transcriptionally regulate around 30% of human genes. Polymorphisms within miRNA processing pathways and binding sites may influence patients’ prognosis through altered gene regulation. Inflammation plays an important role in cancer initiation and progression, and also has shown to impact patients’ clinical outcomes. We first evaluated 240 single nucleotide polymorphisms (SNPs) in miRNA biogenesis genes and predicted binding sites in NSCLC patients to determine associations with clinical outcomes in early-stage (stage I and II) and late-stage (stage III and IV) lung cancer patients, respectively. First, in 535 early-stage patients, after correcting multiple comparisons, FZD4:rs713065 (hazard ratio [HR]:0.46, 95% confidence interval [CI]:0.32-0.65) showed a significant inverse association with survival in early stage surgery-only patients. SP1:rs17695156 (HR:2.22, 95% CI:1.44-3.41) and DROSHA:rs6886834 (HR:6.38, 95% CI:2.49-16.31) conferred increased risk of progression in the all patients and surgery-only populations, respectively. FAS:rs2234978 was significantly associated with improved survival in all patients (HR:0.59, 95% CI:0.44-0.77) and in the surgery plus chemotherapy populations (HR:0.19, 95% CI:0.07-0.46).. Functional genomics analysis demonstrated that this variant creates a miR-651 binding site resulting in altered miRNA regulation of FAS, providing biological plausibility for the observed association. We then analyzed these associations in 598 late-stage patients. After multiple comparison corrections, no SNPs remained significant in the late stage group, while the top SNP NAT1:rs15561 (HR=1.98, 96%CI=1.32-2.94) conferred a significantly increased risk of death in the chemotherapy subgroup. To test the hypothesis that genetic variants in the inflammation-related pathways may be associated with survival in NSCLC patients, we first conducted a three-stage study. In the discovery phase, we investigated a comprehensive panel of 11,930 inflammation-related SNPs in three independent lung cancer populations. A missense SNP (rs2071554) in HLA-DOB was significantly associated with poor survival in the discovery population (HR: 1.46, 95% CI: 1.02-2.09), internal validation population (HR: 1.51, 95% CI: 1.02-2.25), and external validation (HR: 1.52, 95% CI: 1.01-2.29) population. Rs2900420 in KLRK1 was significantly associated with a reduced risk for death in the discovery (HR: 0.76, 95% CI: 0.60-0.96) and internal validation (HR: 0.77, 95% CI: 0.61-0.99) populations, and the association reached borderline significance in the external validation population (HR: 0.80, 95% CI: 0.63-1.02). We also evaluated these inflammation-related SNPs in NSCLC patients in never smokers. Lung cancer in never smokers has been increasingly recognized as distinct disease from that in ever-smokers. A two-stage study was performed using a discovery population from MD Anderson (411 patients) and a validation population from Mayo Clinic (311 patients). Three SNPs (IL17RA:rs879576, BMP8A:rs698141, and STK:rs290229) that were significantly associated with survival were validated (pCD74:rs1056400 and CD38:rs10805347) were borderline significant (p=0.08) in the Mayo Clinic population. In the combined analysis, IL17RA:rs879576 resulted in a 40% reduction in the risk for death (p=4.1 × 10-5 [p=0.61, heterogeneity test]). We also validated a survival tree created in MD Anderson population in the Mayo Clinic population. In conclusion, our results provided strong evidence that genetic variations in specific pathways that examined (miRNA and inflammation pathways) influenced clinical outcomes in NSCLC patients, and with further functional studies, the novel loci have potential to be translated into clinical use.