An empirical evaluation of the Random Forests classifier models for variable selection in a large-scale lung cancer case-control study

Random Forests™ is reported to be one of the most accurate classification algorithms in complex data analysis. It shows excellent performance even when most predictors are noisy and the number of variables is much larger than the number of observations. In this thesis Random Forests was applied to a large-scale lung cancer case-control study. A novel way of automatically selecting prognostic factors was proposed. Also, synthetic positive control was used to validate Random Forests method. Throughout this study we showed that Random Forests can deal with large number of weak input variables without overfitting. It can account for non-additive interactions between these input variables. Random Forests can also be used for variable selection without being adversely affected by collinearities. ^ Random Forests can deal with the large-scale data sets without rigorous data preprocessing. It has robust variable importance ranking measure. Proposed is a novel variable selection method in context of Random Forests that uses the data noise level as the cut-off value to determine the subset of the important predictors. This new approach enhanced the ability of the Random Forests algorithm to automatically identify important predictors for complex data. The cut-off value can also be adjusted based on the results of the synthetic positive control experiments. ^ When the data set had high variables to observations ratio, Random Forests complemented the established logistic regression. This study suggested that Random Forests is recommended for such high dimensionality data. One can use Random Forests to select the important variables and then use logistic regression or Random Forests itself to estimate the effect size of the predictors and to classify new observations. ^ We also found that the mean decrease of accuracy is a more reliable variable ranking measurement than mean decrease of Gini. ^

The effect of specific carotenoids on lung cancer risk

A population-based case-comparison study of histologically confirmed lung cancer among white male and female residents of six Texas coastal counties was conducted from 1979-1982. Dietary information as well as information concerning smoking, alcohol consumption, occupational and residential exposures, and family history of cancer was obtained from 149 living cases and 359 comparison subjects.^ These data support the findings of previous studies that reported a protective association of total carotene intake with lung cancer (OR = 4.07, CI = 3.36-4.78), and no association for total Vitamin A or retinol. Of six specific carotenoids examined, these data reveal a statistically significant protective effect for alpha carotene (OR = 3.58, CI = 2.85-4.31), and an elevated, although non-significant effect for beta carotene (OR = 1.46, CI = 0.61-2.31). No consistent significant effect was found for cryptoxanthin, other xanthins, lutein or lycopene.^ Similar results were found for both males and females, and for both squamous cell and adenocarcinoma subtypes, although loss of power resulting from stratification may have rendered the celltype specific results more imprecise. These results should be considered with caution until confirmed by other studies, however, they suggest the importance of evaluating specific carotenoids in future diet-lung cancer investigations. ^

PIM KINASE INHIBITION: SINGLE AGENT AND COMBINATION APPROACHES IN B-CELL LYMPHOMA

Proviral integration site for Moloney murine leukemia virus (Pim) kinases are Ser/Thr/Tyr kinases. They modulate B-cell development but become oncoproteins and promote cancer development once overexpressed. Containing three isoforms, Pim-1, -2 and -3 are known to phosphorylate various substrates that regulate transcription, translation, cell cycle, and survival pathways in both hematological and solid tumors. Mantle cell lymphoma (MCL) is an aggressive B-cell lymphoma. Elevated Pim kinase levels are common in MCL, and it negatively correlates with patient outcome. SGI-1776 is a small molecule inhibitor selective for Pim-1/-3. We hypothesize that SGI-1776 treatment in MCL will inhibit Pim kinase function, and inhibition of downstream substrates phosphorylation will disrupt transcriptional, translational, and cell cycle processes while promoting apoptosis. SGI-1776 treatment induced moderate to high levels of apoptosis in four MCL cell lines (JeKo-1, Mino, SP-53 and Granta-519) and peripheral blood mononuclear cells (PBMCs) from MCL patients. Phosphorylation of transcription and translation regulators, c-Myc and 4E-BP1 declined in both model systems. Additionally, levels of short-lived Mcl-1 mRNA and protein also decreased and correlated with decline of global RNA synthesis. Collectively, our investigations highlight Pim kinases as viable drug targets in MCL and emphasize their roles in transcriptional and translational regulation. We further investigated a combination strategy using SGI-1776 with bendamustine, an FDA-approved DNA-damaging alkylating agent for treating non-Hodgkin’s lymphoma. We hypothesized this combination will enhance SGI-1776-induced transcription and translation inhibition, while promoting bendamustine-triggered DNA damage and inducing additive to synergistic cytotoxicity in B-cell lymphoma. Bendamustine alone resulted in moderate levels of apoptosis induction in MCL cell lines (JeKo-1 and Mino), and in MCL and splenic marginal zone lymphoma (a type of B-cell lymphoma) primary cells. An additive effect in cell killing was observed when combined with SGI-1776. Expectedly, SGI-1776 effectively decreased global RNA and protein synthesis levels, while bendamustine significantly inhibited DNA synthesis and generated DNA damage response. In combination, intensified inhibitory effects in DNA, RNA and protein syntheses were observed. Together, these data suggested feasibility of using Pim kinase inhibitor in combination with chemotherapeutic agents such as bendamustine in B-cell lymphoma, and provided foundation of their mechanism of actions in lymphoma cells.

Functional studies of *Ras and Bcl-2 oncoproteins in keratinocyte homeostasis and multistep skin carcinogenesis

To assess the effect of deregulated Ha-ras and bcl-2, individually and in combination on epidermal keratinocyte homeostasis and during multistep skin carcinogenesis, we generated skin-specific transgenic mice and keratinocyte transfectants constitutively expressing oncogenic Ha-ras and bcl-2 proteins. The deregulated Ha-ras and bcl-2 expression contributing to homeostatic imbalances in the skin had an additive effect on the probability of tumor development. They were also cooperative in incidence, growth, and latency of tumor formation, and they exhibited synergistic cooperation in malignant transformation of benign papillomas. To explain the homeostatic imbalances by Ha-ras and bcl-2 overexpression in the skin, we investigated the three major cellular processes of proliferation, cell death, and differentiation. Epidermal expression of Bcl-2 retarded keratinocyte proliferation in the epidermis of neonatal mice compared with results for control littermates. Constitutive expression of Ha-ras increased keratinocyte proliferation, and co-expression of bcl-2 modestly suppressed the ras-mediated abnormal proliferation of neonatal keratinocytes. Bcl-2 proteins in keratinocytes protected UV-treated cells from apoptotic cell death regardless of oncogenic ras expression in both non-neoplastic neonatal epidermis and human keratinocyte cell lines. The spontaneous apoptotic index (AI) was also lower in papillomas constitutively expressing bcl-2 compared with the ones that developed in control mice. Ras-overexpressing epidermis, including that in ras/bcl-2 double transgenic mice, had abnormal differentiation patterns compared with controls. The oncogenic ras protein had alterations in both epidermal distribution and the extent of cytokeratin 14 and involucrin expression. Abnormal expression of the hyperproliferation marker cytokeratin 6 and modest down regulation of cytokeratin 1 were also detected. Late appearance of filaggrin was another abnormal phenotype of the ras-expressing epidermis. Overexpression of bcl-2 had no effect on epidermal differentiation. Together, these findings suggest that constitutive expression of oncogenic Ha-ras and bcl-2 are important determinants of epidermal proliferation, viability and differentiation. In summary, our results demonstrated that the disruption of epidermal homeostasis by overexpressed ras and bcl-2 predisposes to hyperplastic growth of the epidermis and to papilloma development and that these proteins with distinct mechanisms for oncogenesis are functionally synergistic for malignant transformation of chemically induced skin carcinogenesis. ^

Predictors of survival in Alzheimer's disease

Alzheimer's disease (AD) is associated with greater mortality and reduced survival among individuals with Alzheimer's disease as compared to those without dementia. It is uncertain how these survival estimates change when the clinical signs and/or symptoms of comorbid conditions are present in individuals' with Alzheimer's disease. Cardiovascular risk factors such as hypertension, hyperlipidemia, congestive heart failure, coronary artery disease, and diabetes mellitus are common conditions in the aged population. Independently, these factors influence mortality and may have an additive effect on reduced survival in an individual with concomitant Alzheimer's disease. The bulk of the evidence from previous research efforts suggests an association between vascular co-morbidities and Alzheimer's disease incidence, but their role in survival remains to be elucidated. The objective of this proposed study was to examine the effects of cardiovascular comorbidities on the survival experience of individuals with probable Alzheimer's disease in order to identify prognostic factors for life expectancy following onset of disease. This study utilized data from the Baylor College of Medicine Alzheimer's Disease Center (ADC) longitudinal study of Alzheimer's disease and other memory disorders. Individuals between the ages of 55-69, 70-79, and ≥80 had a median survival from date of onset of 9.2 years, 8.0 years, and 7.2 years, respectively (p<0.001) and 5.5 years, 4.3 years, and 3.4 years from diagnosis. Sex was the strongest predictor of death from onset of AD, with females having a 30 percent lower risk compared to males. These findings further support the notion that age (both from onset and from diagnosis) and sex are the strongest predictors of survival among those with AD. ^

A STUDY ON THE EPIDEMIOLOGY OF POSITIVE PHYSICAL HEALTH (NHANES I)

HANES 1 detailed sample data were used to operationalize a definition of health in the absence of disease and to describe and compare the characteristics of the normal (healthy) group versus an abnormal (unhealthy) group.^ Parallel screening gave a 3.8 percent prevalence proportion of physical health, with a female:male ratio of 2:1 and younger ages in the healthy group. Statistically significant Mantel-Haenszel gender-age-adjusted odds ratios (MHOR) were estimated among abnormal non-migrants (1.53), skilled workers/unemployed (1.76), annual family incomes of less than $10,000 (1.56), having ever smoked (1.58), and started smoking before 18 years of age (1.58). Significant MHOR were also found for abnormals for health promoting measures: non-iodized salt use (1.94), needed dental care (1.91); and for fair to poor perceived health (4.28), perceiving health problems (2.52), and low energy level (1.68). Significant protective effects for much to moderate recreational exercise (MHOR 0.42) and very active to moderate non-recreational activity (MHOR 0.49) were also obtained. Covariance analysis additive models detected statistically significant higher mean values for abnormals than normals for serum magnesium, hemoglobin, hematocrit, urinary creatinine, and systolic and diastolic blood pressures, and lower values for abnormals than normals for serum iron. No difference was detected for serum cholesterol. Significant non-additive joint effects were found for body mass index.^ The results suggest positive physical health can be measured with cross-sectional survey data. Gender differentials, and associations between ecologic, socioeconomic, hazardous risk factors, health promoting activities and physical health are in general agreement with published findings on studies of morbidity. Longitudinal prospective studies are suggested to establish the direction of the associations and to enhance present knowledge of health and its promoting factors. ^

TRIM24-REGULATED ESTROGEN RESPONSE IS DEPENDENT ON SPECIFIC HISTONE MODIFICATIONS IN BREAST CANCER CELLS

In this dissertation, I discovered that function of TRIM24 as a co-activator of ERα-mediated transcriptional activation is dependent on specific histone modifications in tumorigenic human breast cancer-derived MCF7 cells. In the first part, I proved that TRIM24-PHD finger domain, which recognizes unmethylated histone H3 lysine K4 (H3K4me0), is critical for ERα-regulated transcription. Therefore, when LSD1-mediated demethylation of H3K4 is inhibited, activation of TRIM24-regulated ERα target genes is greatly impaired. Importantly, I demonstrated that TRIM24 and LSD1 are cyclically recruited to estrogen responsive elements (EREs) in a time-dependent manner upon estrogen induction, and depletion of their expression exert corresponding time-dependent effect on target gene activation. I also identified that phosphorylation of histone H3 threonine T6 disrupts TRIM24 from binding to the chromatin and from activating ERα-regulated targets. In the second part, I revealed that TRIM24 depletion has additive effect to LSD1 inhibitor- and Tamoxifen-mediated reduction in survival and proliferation in breast cancer cells.

THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.

Molecular and evolutionary genetics of the X -linked visual pigment genes in humans and New World monkeys

Normal humans have one red and at least one green visual pigment genes. These genes are tightly linked as tandem repeats on the X chromosome and each of them has six exons. There is only one X-linked visual pigment gene in New World monkeys (NWMs) but the locus has three polymorphic alleles encoding red, yellow and green visual pigments, respectively. The spectral properties of the squirrel monkey and the marmoset (both NWMs) have been studied and partial sequences of the three alleles are available. To study the evolutionary history of these X-linked opsin genes in humans and NWMs, coding and intron sequences of the three squirrel monkey alleles and the three marmoset alleles were amplified by PCR followed by subcloning and sequencing. Introns 2 and 4 of the human red and green pigment genes were also sequenced. The results obtained are as follows: (1) The sequences of introns 2 and 4 of the human red and green opsin genes are significantly more similar between the two genes than are coding sequences, contrary to the usual situation where coding regions are better conserved in evolution than are introns. The high similarities in the two introns are probably due to recent gene conversion events during evolution of the human lineage. (2) Phylogenetic analysis of both intron and exon sequences indicates that the phylogenetic tree of the available primate opsin genes is the same as the species tree. The two human genes were derived from a gene duplication event after the divergence of the human and NWM lineages. The three alleles in each of the two NWM species diverged after the split of the two NWMs but have persisted in the population for at least 5 million years. (3) Allelic gene conversion might have occurred between the three squirrel monkey alleles. (4) A model of additive effect of hydroxyl-bearing amino acids on spectral tuning is proposed by treating some unknown variables as groups. Under the assumption that some residues have no effect, it is found that at least five amino acid residues, at positions 178 (3 nm), 180 (5 nm), 230 ($-$4 nm), 277 (9 nm) and 285 (13 nm), have linear spectral tuning effects. (5) Adaptive evolution of the opsin genes to different spectral peaks was observed at four residues that are important for spectral tuning. ^

THE EFFECT OF PERSON-ENVIRONMENT FIT ON THE PSYCHOLOGICAL WELL-BEING OF PSYCHIATRIC AIDES WORKING IN A STATE HOSPITAL (COPING RESOURCES, SOCIAL SUPPORT, OCCUPATIONAL STRESS, LOCUS-OF-CONTROL)

This study addressed two purposes: (1) to determine the effect of person-environment fit on the psychological well-being of psychiatric aides and (2) to determine what role the coping resources of social support and control have on the above relationship. Two hundred and ten psychiatric aides working in a state hospital in Texas responded to a questionnaire pertaining to these issues.^ Person-environment fit, as a measure of occupational stress, was assessed through a modified version of the Work Environment Scale (WES). The WES subscales used in this study were: involvement, autonomy, job pressure, job clarity, and physical comfort. Psychological well-being was measured with the General Well-Being Schedule which was developed by the National Center for Health Statistics. Co-worker and supervisor support were measured through the WES and finally, control was assessed through Rotter's Locus of Control Scale.^ The results of this study were as follows: (1) all person-environment (p-e) dimensions appeared to have linear relationships with psychological well-being; (2) the p-e fit - well-being relationship did not appear to be confounded by demographic factors; (3) all p-e fit dimensions were significantly related to well-being except for autonomy; (4) p-e fit was more strongly related to well-being than the environmental measure alone; (5) supervisor support and non-work related support were found to have additive effects on the relationship between p-e fit and well-being, however no interaction or buffering effects were observed; (6) locus of control was found to have additive effects in the prediction of well-being and showed interactive effects with work pressure, involvement and physical comfort; and (7) the testing of the overall study model which included many of the components mentioned above yielded an R('2) = .27.^ Implications of these findings are discussed, future research suggested and applications proposed. ^