Hormone therapy and physical activity as predictors of lung cancer risk and survival: The California Teachers Study

Cigarette smoking is responsible for the majority of lung cancer cases worldwide; however, a proportion of never smokers still develop lung cancer over their lifetime, prompting investigation into additional factors that may modify lung cancer incidence, as well as mortality. Although hormone therapy (HT), physical activity (PA), and lung cancer have been previously examined, the associations remain unclear. This study investigated exposure to HT and PA that may modulate underlying mechanisms of lung cancer etiology and progression among women by using existing, de-identified data from the California Teachers Study (CTS).^ The CTS cohort, established in 1995–1996, has 133,479 active and retired female teachers and administrators, recruited through the California State Teachers Retirement System, and followed annually for cancer diagnosis, death, and change of address. Each woman enrolled in the CTS returned a questionnaire covering a wide variety of issues related to cancer risk and women's health, including recent and past HT use and physical activity, as well as active and environmental cigarette smoke exposure. Complete data to assess the associations between HT and lung cancer risk and survival were available for 60,592 postmenopausal women. Between 1995 and 2007, 727 of these women were diagnosed with invasive lung cancer; 441 of these died. Complete data to assess the associations between PA and lung cancer risk and survival were available for 118,513 women. Between 1995 and 2007, 853 of these women were diagnosed with invasive lung cancer; 516 of these died.^ After careful adjustment for smoking habits and other potential confounders, no measure of HT use was associated with lung cancer risk; however, any HT use (vs. no use) was associated with a decrease in lung-cancer-specific mortality. Specifically, among women who only used estrogen (E-only), decreases in lung cancer mortality were seen for recent use, but not for former use; no association was observed for estrogen plus progestin (E+P). Furthermore, among former users of HT, a statistically significant decrease in lung cancer mortality was observed for E-only use within 5 years prior to baseline, but not for E-only use >5 years prior to baseline. Neither long-term recreational PA nor recent recreational PA alone were associated with lung cancer risk; however, among women with a BMI<25 and ever smokers, high long-term moderate+strenuous PA was associated with a decrease in lung cancer risk. Women with non-local disease showed a decrease in lung cancer mortality associated with increasing duration of strenuous long-term activity, and 1.50-3.00 h/wk/y of recent moderate or recent strenuous PA. Long-term moderate PA was associated with decreased lung cancer mortality in never smokers, whereas recent moderate PA was associated with increased lung cancer mortality in current smokers. ^ Placing our findings in the context of the current literature, HT does not appear to be associated with lung cancer risk and previous studies reporting a protective effect of HT use on lung cancer risk may be subject to residual confounding by smoking. Looking at our findings regarding PA overall, the evidence still remains inconclusive regarding whether or not physical activity influence lung cancer risk or mortality. Our results suggest that recreational PA may associated with decreased lung cancer risk among women with BMI<25 and ever smoking-women; however, residual confounding by smoking should be strongly considered. To our knowledge, this is the first study to investigate lifetime recreational PA and lung cancer mortality among women. Our results contribute to the growing body of knowledge regarding non-smoking-related risk factors for lung cancer incidence and mortality among women. Given the potential clinical and interventional significance, further study and validation of these findings is warranted.^

Logistic regression with Markov chains as covariates

The tobacco-specific nitrosamine 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) is an obvious carcinogen for lung cancer. Since CBMN (Cytokinesis-blocked micronucleus) has been found to be extremely sensitive to NNK-induced genetic damage, it is a potential important factor to predict the lung cancer risk. However, the association between lung cancer and NNK-induced genetic damage measured by CBMN assay has not been rigorously examined. ^ This research develops a methodology to model the chromosomal changes under NNK-induced genetic damage in a logistic regression framework in order to predict the occurrence of lung cancer. Since these chromosomal changes were usually not observed very long due to laboratory cost and time, a resampling technique was applied to generate the Markov chain of the normal and the damaged cell for each individual. A joint likelihood between the resampled Markov chains and the logistic regression model including transition probabilities of this chain as covariates was established. The Maximum likelihood estimation was applied to carry on the statistical test for comparison. The ability of this approach to increase discriminating power to predict lung cancer was compared to a baseline "non-genetic" model. ^ Our method offered an option to understand the association between the dynamic cell information and lung cancer. Our study indicated the extent of DNA damage/non-damage using the CBMN assay provides critical information that impacts public health studies of lung cancer risk. This novel statistical method could simultaneously estimate the process of DNA damage/non-damage and its relationship with lung cancer for each individual.^

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

Intensive care unit utilization among patients with cancer at the end of life

Over the last 2 decades, survival rates in critically ill cancer patients have improved. Despite the increase in survival, the intensive care unit (ICU) continues to be a location where end-of-life care takes place. More than 20% of deaths in the United States occur after admission to an ICU, and as baby boomers reach the seventh and eighth decades of their lives, the volume of patients in the ICU is predicted to rise. The aim of this study was to evaluate intensive care unit utilization among patients with cancer who were at the end of life. End of life was defined using decedent and high-risk cohort study designs. The decedent study evaluated characteristics and ICU utilization during the terminal hospital stay among patients who died at The University of Texas MD Anderson Cancer Center during 2003-2007. The high-risk cohort study evaluated characteristics and ICU utilization during the index hospital stay among patients admitted to MD Anderson during 2003-2007 with a high risk of in-hospital mortality. Factors associated with higher ICU utilization in the decedent study included non-local residence, hematologic and non-metastatic solid tumor malignancies, malignancy diagnosed within 2 months, and elective admission to surgical or pediatric services. Having a palliative care consultation on admission was associated with dying in the hospital without ICU services. In the cohort of patients with high risk of in-hospital mortality, patients who went to the ICU were more likely to be younger, male, with newly diagnosed non-metastatic solid tumor or hematologic malignancy, and admitted from the emergency center to one of the surgical services. A palliative care consultation on admission was associated with a decreased likelihood of having an ICU stay. There were no differences in ethnicity, marital status, comorbidities, or insurance status between patients who did and did not utilize ICU services. Inpatient mortality probability models developed for the general population are inadequate in predicting in-hospital mortality for patients with cancer. The following characteristics that differed between the decedent study and high-risk cohort study can be considered in future research to predict risk of in-hospital mortality for patients with cancer: ethnicity, type and stage of malignancy, time since diagnosis, and having advance directives. Identifying those at risk can precipitate discussions in advance to ensure care remains appropriate and in accordance with the wishes of the patient and family.^

The local safety net, health insurance, preventable hospitalizations and costs: Harris County (Houston), Texas

The aim of this study was to examine the association between determinants of access to healthcare and preventable hospitalizations, based on Davidson et al.'s framework for evaluating the effects of individual and community determinants on access to healthcare. The study population consisted of the low income, non-elderly, hospitalized adults residing in Harris County, Texas in 2004. The objectives of this study were to examine the proportion of the variance in preventable hospitalizations at the ZIP-code level, to analyze the association between the proximity to the nearest safety net clinic and preventable hospitalizations, to examine how the safety net capacity relates to preventable hospitalizations, to compare the relative strength of the associations of health insurance and the proximity to the nearest safety net clinic with preventable hospitalizations, and to estimate and compare the costs of preventable hospitalizations in Harris County with the average cost in the literature. The data were collected from Texas Health Care Information Collection (2004), Census 2000, and Project Safety Net (2004). A total of 61,841 eligible individuals were included in the final data analysis. A random-intercept multi-level model was constructed with two different levels of data: the individual level and the ZIP-code level. The results of this study suggest that ZIP-code characteristics explain about two percent of the variance in preventable hospitalizations and safety net capacity was marginally significantly associated with preventable hospitalizations (p= 0.062). Proximity to the nearest safety net clinic was not related to preventable hospitalizations; however, health insurance was significantly associated with a decreased risk of preventable hospitalization. The average direct cost was $6,466 per preventable hospitalization, which is significantly different from reports in the literature. ^

A model for use by local public health departments to evaluate pandemic influenza plans

The events of the 1990's and early 2000's demonstrated the need for effective planning and response to natural and man-made disasters. One of those potential natural disasters is pandemic flu. Once defined, the CDC stated that program, or plan, effectiveness is improved through the process of program evaluation. (Centers for Disease Control and Prevention, 1999) Program evaluation should be accomplished not only periodically, but in the course of routine administration of the program. (Centers for Disease Control and Prevention, 1999) Accomplishing this task for a "rare, but significant event" is challenging. (Herbold, John R., PhD., 2008) To address this challenge, the RAND Corporation (under contract to the CDC) developed the "Facilitated Look-Backs" approach that was tested and validated at the state level. (Aledort et al., 2006).^ Nevertheless, no comprehensive and generally applicable pandemic influenza program evaluation tool or model is readily found for use at the local public health department level. This project developed such a model based on the "Facilitated Look-Backs" approach developed by RAND Corporation. (Aledort et al., 2006) Modifications to the RAND model included stakeholder additions, inclusion of all six CDC program evaluation steps, and suggestions for incorporating pandemic flu response plans in seasonal flu management implementation. Feedback on the model was then obtained from three LPHD's—one rural, one suburban, and one urban. These recommendations were incorporated into the final model. Feedback from the sites also supported the assumption that this model promotes the effective and efficient evaluation of both pandemic flu and seasonal flu response by reducing redundant evaluations of pandemic flu plans, seasonal flu plans, and funding requirement accountability. Site feedback also demonstrated that the model is comprehensive and flexible, so it can be adapted and applied to different LPHD needs and settings. It also stimulates evaluation of the major issues associated with pandemic flu planning. ^ The next phase in evaluating this model should be to apply it in a program evaluation of one or more LPHD's seasonal flu response that incorporates pandemic flu response plans.^

Bayesian and survival model for tumor relapse status and disease-specific survival, with applications to breast cancer

Breast cancer is the most common non-skin cancer and the second leading cause of cancer-related death in women in the United States. Studies on ipsilateral breast tumor relapse (IBTR) status and disease-specific survival will help guide clinic treatment and predict patient prognosis.^ After breast conservation therapy, patients with breast cancer may experience breast tumor relapse. This relapse is classified into two distinct types: true local recurrence (TR) and new ipsilateral primary tumor (NP). However, the methods used to classify the relapse types are imperfect and are prone to misclassification. In addition, some observed survival data (e.g., time to relapse and time from relapse to death)are strongly correlated with relapse types. The first part of this dissertation presents a Bayesian approach to (1) modeling the potentially misclassified relapse status and the correlated survival information, (2) estimating the sensitivity and specificity of the diagnostic methods, and (3) quantify the covariate effects on event probabilities. A shared frailty was used to account for the within-subject correlation between survival times. The inference was conducted using a Bayesian framework via Markov Chain Monte Carlo simulation implemented in softwareWinBUGS. Simulation was used to validate the Bayesian method and assess its frequentist properties. The new model has two important innovations: (1) it utilizes the additional survival times correlated with the relapse status to improve the parameter estimation, and (2) it provides tools to address the correlation between the two diagnostic methods conditional to the true relapse types.^ Prediction of patients at highest risk for IBTR after local excision of ductal carcinoma in situ (DCIS) remains a clinical concern. The goals of the second part of this dissertation were to evaluate a published nomogram from Memorial Sloan-Kettering Cancer Center, to determine the risk of IBTR in patients with DCIS treated with local excision, and to determine whether there is a subset of patients at low risk of IBTR. Patients who had undergone local excision from 1990 through 2007 at MD Anderson Cancer Center with a final diagnosis of DCIS (n=794) were included in this part. Clinicopathologic factors and the performance of the Memorial Sloan-Kettering Cancer Center nomogram for prediction of IBTR were assessed for 734 patients with complete data. Nomogram for prediction of 5- and 10-year IBTR probabilities were found to demonstrate imperfect calibration and discrimination, with an area under the receiver operating characteristic curve of .63 and a concordance index of .63. In conclusion, predictive models for IBTR in DCIS patients treated with local excision are imperfect. Our current ability to accurately predict recurrence based on clinical parameters is limited.^ The American Joint Committee on Cancer (AJCC) staging of breast cancer is widely used to determine prognosis, yet survival within each AJCC stage shows wide variation and remains unpredictable. For the third part of this dissertation, biologic markers were hypothesized to be responsible for some of this variation, and the addition of biologic markers to current AJCC staging were examined for possibly provide improved prognostication. The initial cohort included patients treated with surgery as first intervention at MDACC from 1997 to 2006. Cox proportional hazards models were used to create prognostic scoring systems. AJCC pathologic staging parameters and biologic tumor markers were investigated to devise the scoring systems. Surveillance Epidemiology and End Results (SEER) data was used as the external cohort to validate the scoring systems. Binary indicators for pathologic stage (PS), estrogen receptor status (E), and tumor grade (G) were summed to create PS+EG scoring systems devised to predict 5-year patient outcomes. These scoring systems facilitated separation of the study population into more refined subgroups than the current AJCC staging system. The ability of the PS+EG score to stratify outcomes was confirmed in both internal and external validation cohorts. The current study proposes and validates a new staging system by incorporating tumor grade and ER status into current AJCC staging. We recommend that biologic markers be incorporating into revised versions of the AJCC staging system for patients receiving surgery as the first intervention.^ Chapter 1 focuses on developing a Bayesian method to solve misclassified relapse status and application to breast cancer data. Chapter 2 focuses on evaluation of a breast cancer nomogram for predicting risk of IBTR in patients with DCIS after local excision gives the statement of the problem in the clinical research. Chapter 3 focuses on validation of a novel staging system for disease-specific survival in patients with breast cancer treated with surgery as the first intervention. ^