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.

ROLE OF STAT3 IN KERATINOCYTE STEM CELLS DURING SKIN TUMORIGENESIS

STATs play crucial roles in a wide variety of biological functions, including development, proliferation, differentiation, migration and in cancer development. In the present study, we examined the impact of Stat3 deletion or activation on behavior of keratinocytes, including keratinocyte stem cells (KSCs). Deletion of Stat3 specifically in the bulge region of the hair follicle using K15.CrePR1 X Stat3fl/fl mice led to decreased tumor development by altering survival of bulge region KSCs. To further understand the role of KSCs in skin tumorigenesis, K5.Stat3C transgenic (Tg) mice which express a constitutively active/dimerized form of Stat3 called Stat3C via the bovine keratin 5 (K5) promoter were studied. The number of CD34 and α6 integrin positive cells was significantly reduced in Tg mice as compared to non-transgenic (NTg) littermates. There was a concomitant increase in the progenitor populations (Lgr-6, Lrig-1 and Sca-1) in the Tg mice vs. the stem cell population (CD34 and Keratin15). To investigate the mechanism underlying the increase in the progenitor population at the expense of bulge region KSCs we examined if Stat3C expression was involved in inducing migration of the bulge region KSCs. There was altered β-catenin and α6-integrin expression in the hair follicles of Tg mice, which may have contributed to reduced adhesive interactions between the epithelial cells and the basement membrane facilitating migration out of the niche. To further study the effect of Stat3 on differentiation of keratinocytes we analyzed the epidermal keratinocytes in K5.Cre X Stat3fl/fl mice. There was an increase in the expression of epidermal differentiation markers in the Stat3 knockout mice. These data suggest that deletion of Stat3 in the epidermis and hair follicle induced differentiation in these cells. Preliminary studies done with the BK5.Stat3C mouse model suggests that multiple hair follicle stem/progenitor populations may be involved in skin tumor development and progression in this model of skin tumorigenesis. Overall, these data suggest that Stat3 plays an important role in differentiation as well as migration of keratinocytes and that these effects may play a role during epithelial carcinogenesis.

Targeting AKT/mTOR Signaling Pathways During Murine Skin Tumor Promotion and the Impact of Dietary Energy Balance Manipulation

The prevalence of obesity has continued to rise over the last several decades in the United States lending to overall increases in risk for chronic diseases including many types of cancer. In contrast, reduction in energy consumption via calorie restriction (CR) has been shown to be a potent inhibitor of carcinogenesis across a broad range of species and tumor types. Previous data has demonstrated differential signaling through Akt and mTOR via the IGF-1R and other growth factor receptors across the diet-induced obesity (DIO)/CR spectrum. Furthermore, mTORC1 is known to be regulated directly via nutrient availability, supporting its role in the link between epithelial carcinogenesis and diet-induced obesity. In an effort to better understand the importance of mTORC1 in the context of both positive and negative energy balance during epithelial carcinogenesis, we have employed the use of specific pharmacological inhibitors, rapamycin (mTORC1 inhibitor) and metformin (AMPK activator) to target mTORC1 or various components of this pathway during skin tumor promotion. Two-stage skin carcinogenesis studies demonstrated that mTORC1 inhibition via rapamycin, metformin or combination treatments greatly inhibited skin tumor development in normal, overweight and obese mice. Furthermore, mechanisms by which these chemopreventive agents may be exerting their anti-tumor effects were explored. In addition, the effect of these compounds on the epidermal proliferative response was analyzed and drastic decreases in epidermal hyperproliferation and hyperplasia were found. Rapamycin also inhibited dermal inflammatory cell infiltration in a dose-dependent manner. Both compounds also blocked or attenuated TPA-induced signaling through epidermal mTORC1 as well as several downstream targets. In addition, inhibition of this pathway by metformin appeared to be, at least in part, dependent on AMPK activation in the skin. Overall, the data indicate that pharmacological strategies targeting this pathway offset the tumor-enhancing effects of DIO and may serve as possible CR mimetics. They suggest that mTORC1 contributes significantly to the process of skin tumor promotion, specifically during dietary energy balance effects. Exploiting the mechanistic information underlying dietary energy balance responsive pathways will help translate decades of research into effective strategies for prevention of epithelial carcinogenesis.

Adaptive clinical trial design for targeted therapy development

The development of targeted therapy involve many challenges. Our study will address some of the key issues involved in biomarker identification and clinical trial design. In our study, we propose two biomarker selection methods, and then apply them in two different clinical trial designs for targeted therapy development. In particular, we propose a Bayesian two-step lasso procedure for biomarker selection in the proportional hazards model in Chapter 2. In the first step of this strategy, we use the Bayesian group lasso to identify the important marker groups, wherein each group contains the main effect of a single marker and its interactions with treatments. In the second step, we zoom in to select each individual marker and the interactions between markers and treatments in order to identify prognostic or predictive markers using the Bayesian adaptive lasso. In Chapter 3, we propose a Bayesian two-stage adaptive design for targeted therapy development while implementing the variable selection method given in Chapter 2. In Chapter 4, we proposed an alternate frequentist adaptive randomization strategy for situations where a large number of biomarkers need to be incorporated in the study design. We also propose a new adaptive randomization rule, which takes into account the variations associated with the point estimates of survival times. In all of our designs, we seek to identify the key markers that are either prognostic or predictive with respect to treatment. We are going to use extensive simulation to evaluate the operating characteristics of our methods.^

Instrumental variable method for the causal relationship between soluble receptor for advanced glycation end-products (sRAGE) and risk of pancreatic cancer

This thesis project is motivated by the potential problem of using observational data to draw inferences about a causal relationship in observational epidemiology research when controlled randomization is not applicable. Instrumental variable (IV) method is one of the statistical tools to overcome this problem. Mendelian randomization study uses genetic variants as IVs in genetic association study. In this thesis, the IV method, as well as standard logistic and linear regression models, is used to investigate the causal association between risk of pancreatic cancer and the circulating levels of soluble receptor for advanced glycation end-products (sRAGE). Higher levels of serum sRAGE were found to be associated with a lower risk of pancreatic cancer in a previous observational study (255 cases and 485 controls). However, such a novel association may be biased by unknown confounding factors. In a case-control study, we aimed to use the IV approach to confirm or refute this observation in a subset of study subjects for whom the genotyping data were available (178 cases and 177 controls). Two-stage IV method using generalized method of moments-structural mean models (GMM-SMM) was conducted and the relative risk (RR) was calculated. In the first stage analysis, we found that the single nucleotide polymorphism (SNP) rs2070600 of the receptor for advanced glycation end-products (AGER) gene meets all three general assumptions for a genetic IV in examining the causal association between sRAGE and risk of pancreatic cancer. The variant allele of SNP rs2070600 of the AGER gene was associated with lower levels of sRAGE, and it was neither associated with risk of pancreatic cancer, nor with the confounding factors. It was a potential strong IV (F statistic = 29.2). However, in the second stage analysis, the GMM-SMM model failed to converge due to non- concaveness probably because of the small sample size. Therefore, the IV analysis could not support the causality of the association between serum sRAGE levels and risk of pancreatic cancer. Nevertheless, these analyses suggest that rs2070600 was a potentially good genetic IV for testing the causality between the risk of pancreatic cancer and sRAGE levels. A larger sample size is required to conduct a credible IV analysis.^

Modulation of cytochrome P450 enzyme activity and PAH -induced skin carcinogenesis by naturally occurring coumarins

The present study was designed to determine the potential anticarcinogenic activity of naturally occurring coumarins and their mechanism of action. The results indicated that several naturally occurring coumarins including bergamottin, coriandrin, imperatorin, isopimpinellin, and ostruthin, to which humans are routinely exposed in the diet, were effective inhibitors and/or inactivators of CYP1A1-mediated ethoxyresorufin-O-dealkylase (EROD) or CYP2B1-mediated pentoxyresorufin-O-dealkylase (PROD) in mouse liver microsomes. In addition, bergamottin and corandrin were also found to be inhibitors of purified human P450 1A1 in vitro. Further studies with coriandrin revealed that this compound was a mechanism-based inactivator of P450 1A1 and covalently bound to the P450 1A1 apoprotein. In cultured mouse keratinocytes, bergamottin and coriandrin effectively inhibited the B(a) P metabolism and significantly decreased covalent binding of B(a) P and DMBA to keratinocyte DNA and anti-diol-epoxide-DNA adducts derived from both B(a) P and DMBA in keratinocytes. The data from in vivo experiments showed that bergamottin and coriandrin were potent inhibitors of covalent binding of B (a) P to epidermal DNA and the formation of (+) anti BPDE-DNA adduct, whereas imperatorin and isopimpinellin were more potent inhibitors of covalent binding of DMBA to epidermal DNA. The ability of coumarins to inhibit covalent binding of B (a) P to DNA in mouse epidermis was positively correlated with their inhibitory effect P450 1A1 in vitro, while the inhibitory effect of coumarins on covalent binding of DMBA to epidermal DNA was positively correlated with their inhibitory effects on P450 2B1 and negatively to their inhibitory activity toward P450 1A1. The data from tumor experiments indicated that bergamottin, ostruthin, and coriandrin inhibited tumor initiation by B (a) P in a two-stage carcinogenesis protocol. Bergamottin was most effective in this regard and produced a dose dependent inhibition of papilloma formation in these experiments. In addition, imperatorin was an effective inhibitor of skin tumorigenesis induced by DMBA in SENCAR mouse skin using both a two-stage and a complete carcinogenesis protocol. At dose levels higher than those effective against DMBA, imperatorin also inhibited tumor initiation by B (a) P. The results to date demonstrate that several naturally occurring coumarins possess the ability to block tumor initiation and tumorigenesis by PAHs such as B (a) P and DMBA through inhibition of the P450s involved in the metabolic activation of these hydrocarbons. A working model for the involvement of specific P450s in the metabolic activation of these two PAHs was proposed. ^