Development of a Bayesian Joint Logistic Model to Better Study the Association between Haplotypes and Disease

In 2011, there will be an estimated 1,596,670 new cancer cases and 571,950 cancer-related deaths in the US. With the ever-increasing applications of cancer genetics in epidemiology, there is great potential to identify genetic risk factors that would help identify individuals with increased genetic susceptibility to cancer, which could be used to develop interventions or targeted therapies that could hopefully reduce cancer risk and mortality. In this dissertation, I propose to develop a new statistical method to evaluate the role of haplotypes in cancer susceptibility and development. This model will be flexible enough to handle not only haplotypes of any size, but also a variety of covariates. I will then apply this method to three cancer-related data sets (Hodgkin Disease, Glioma, and Lung Cancer). I hypothesize that there is substantial improvement in the estimation of association between haplotypes and disease, with the use of a Bayesian mathematical method to infer haplotypes that uses prior information from known genetics sources. Analysis based on haplotypes using information from publically available genetic sources generally show increased odds ratios and smaller p-values in both the Hodgkin, Glioma, and Lung data sets. For instance, the Bayesian Joint Logistic Model (BJLM) inferred haplotype TC had a substantially higher estimated effect size (OR=12.16, 95% CI = 2.47-90.1 vs. 9.24, 95% CI = 1.81-47.2) and more significant p-value (0.00044 vs. 0.008) for Hodgkin Disease compared to a traditional logistic regression approach. Also, the effect sizes of haplotypes modeled with recessive genetic effects were higher (and had more significant p-values) when analyzed with the BJLM. Full genetic models with haplotype information developed with the BJLM resulted in significantly higher discriminatory power and a significantly higher Net Reclassification Index compared to those developed with haplo.stats for lung cancer. Future analysis for this work could be to incorporate the 1000 Genomes project, which offers a larger selection of SNPs can be incorporated into the information from known genetic sources as well. Other future analysis include testing non-binary outcomes, like the levels of biomarkers that are present in lung cancer (NNK), and extending this analysis to full GWAS studies.

Household Hardships, Public Programs, and Their Associations with the Health and Development of Very Young Children: Insights from Children’s HealthWatch

America’s low-income families struggle to protect their children from multiple threats to their health and growth. Many research and advocacy groups explore the health and educational effects of food insecurity, but less is known about these effects on very young children. Children’s HealthWatch, a group of pediatric clinicians and public health researchers, has continuously collected data on the effects of food insecurity alone and in conjunction with other household hardships since 1998. The group’s peer reviewed research has shown that a number of economic risks at the household level, including food, housing and energy insecurity, tend to be correlated. These insecurities alone or in conjunction increase the risk that a young child will suffer various negative health consequences, including increases in lifetime hospitalizations, parental report of fair or poor health,1 or risk for developmental delays.2 Child food insecurity is an incremental risk indicator above and beyond the risk imposed by household-level food insecurity. The Children’sHealthwatch research also suggests public benefits programs modify some of these effects for families experiencing hardships. This empirical evidence is presented in a variety of public venues outside the usual scientific settings, such as congressional hearings, to support the needs of America’s most vulnerable population through policy change. Children’s HealthWatch research supports legislative solutions to food insecurity, including sustained funding for public programs and re-evaluation of the use of the Thrifty Food Plan as the basis of SNAP benefits calculations. Children’s HealthWatch is one of many models to support the American Academy of Pediatrics’ call to “stand up, speak up, and step up for children.”3 No isolated group or single intervention will solve child poverty or multiple hardships. However, working collaboratively each group has a role to play in supporting the health and well-being of young children and their families. 1. Cook JT, Frank DA, Berkowitz C, et al. Food insecurity is associated with adverse health outcomes among human infants and toddlers. J Nutr. 2004;134:1432-1438. 2. Rose-Jacobs R, Black MM, Casey PH, et al. Household food insecurity: associations with at-risk infant and toddler development. Pediatrics. 2008;121:65-72. 3. AAP leader says to stand up, speak up, and step up for child health [news release]. Boston, MA: American Academy of Pediatrics; October 11, 2008. http://www2.aap.org/pressroom/nce/nce08childhealth.htm. Accessed January 1, 2012.

Mouse models and quantitative trait loci of susceptibility to bleomycin- and radiation-induced pulmonary fibrosis

Radiotherapy involving the thoracic cavity and chemotherapy with the drug bleomycin are both dose limited by the development of pulmonary fibrosis. From evidence that there is variation in the population in susceptibility to pulmonary fibrosis, and animal data, it was hypothesized that individual variation in susceptibility to bleomycin-induced, or radiation-induced, pulmonary fibrosis is, in part, genetically controlled. In this thesis a three generation mouse genetic model of C57BL/6J (fibrosis prone) and C3Hf/Kam (fibrosis resistant) mouse strains and F1 and F2 (F1 intercross) progeny derived from the parental strains was developed to investigate the genetic basis of susceptibility to fibrosis. In the bleomycin studies the mice received 100 mg/kg (125 for females) of bleomycin, via mini osmotic pump. The animals were sacrificed at eight weeks following treatment or when their breathing rate indicated respiratory distress. In the radiation studies the mice were given a single dose of 14 or 16 Gy (Co$\sp{60})$ to the whole thorax and were sacrificed when moribund. The phenotype was defined as the percent of fibrosis area in the left lung as quantified with image analysis of histological sections. Quantitative trait loci (QTL) mapping was used to identify the chromosomal location of genes which contribute to susceptibility to bleomycin-induced pulmonary fibrosis in C57BL/6J mice compared to C3Hf/Kam mice and to determine if the QTL's which influence susceptibility to bleomycin-induced lung fibrosis in these progenitor strains could be implicated in susceptibility to radiation-induced lung fibrosis. For bleomycin, a genome wide scan revealed QTL's on chromosome 17, at the MHC, (LOD = 11.7 for males and 7.2 for females) accounting for approximately 21% of the phenotypic variance, and on chromosome 11 (LOD = 4.9), in male mice only, adding 8% of phenotypic variance. The bleomycin QTL on chromosome 17 was also implicated for susceptibility to radiation-induced fibrosis (LOD = 5.0) and contributes 7% of the phenotypic variance in the radiation study. In conclusion, susceptibility to both bleomycin-induced and radiation-induced pulmonary fibrosis are heritable traits, and are influenced by a genetic factor which maps to a genomic region containing the MHC. ^

Statistical characterization and development of summary measures of non-normal distributions of nuclear morphometry data from prostate cancer cells for use as prognostic indicators

Nuclear morphometry (NM) uses image analysis to measure features of the cell nucleus which are classified as: bulk properties, shape or form, and DNA distribution. Studies have used these measurements as diagnostic and prognostic indicators of disease with inconclusive results. The distributional properties of these variables have not been systematically investigated although much of the medical data exhibit nonnormal distributions. Measurements are done on several hundred cells per patient so summary measurements reflecting the underlying distribution are needed.^ Distributional characteristics of 34 NM variables from prostate cancer cells were investigated using graphical and analytical techniques. Cells per sample ranged from 52 to 458. A small sample of patients with benign prostatic hyperplasia (BPH), representing non-cancer cells, was used for general comparison with the cancer cells.^ Data transformations such as log, square root and 1/x did not yield normality as measured by the Shapiro-Wilks test for normality. A modulus transformation, used for distributions having abnormal kurtosis values, also did not produce normality.^ Kernel density histograms of the 34 variables exhibited non-normality and 18 variables also exhibited bimodality. A bimodality coefficient was calculated and 3 variables: DNA concentration, shape and elongation, showed the strongest evidence of bimodality and were studied further.^ Two analytical approaches were used to obtain a summary measure for each variable for each patient: cluster analysis to determine significant clusters and a mixture model analysis using a two component model having a Gaussian distribution with equal variances. The mixture component parameters were used to bootstrap the log likelihood ratio to determine the significant number of components, 1 or 2. These summary measures were used as predictors of disease severity in several proportional odds logistic regression models. The disease severity scale had 5 levels and was constructed of 3 components: extracapsulary penetration (ECP), lymph node involvement (LN+) and seminal vesicle involvement (SV+) which represent surrogate measures of prognosis. The summary measures were not strong predictors of disease severity. There was some indication from the mixture model results that there were changes in mean levels and proportions of the components in the lower severity levels. ^

Development of prostate cancer gene therapy

Osseous metastases account for most of the morbidity and mortality associated with prostate cancer, for which there are currently no effective therapies. In the skeletal metastatic environment, neoplastic prostatic epithelial cells interact in a bidirectional stimulatory manner with osteoblastic stromal cells. Similarly, the presence of osteoblastic cells is essential for the survival and maintenance of intraosseous prostate cancer cells. In this thesis, I have developed novel gene therapy strategies for the treatment of androgen-independent human prostate cancers in experimental animal models. First, Ad-CMV-p53, a recombinant adenovirus (Ad) containing p53 tumor suppressor gene driven by the universal cytomegalovirus promoter, was effective in inhibiting prostate cancer cell growth, and direct intratumoral injections of Ad-CMV-p53 resulted in tumor regression. Second, because prostate cancer cells as well as osteoblastic cells produce osteocalcin (OC), OC promoter mediated tissue/tumor specific toxic gene therapy is developed to interrupt stromal-epithelial communications by targeting both cell types. Ad-OC-TK, a recombinant Ad containing the herpes simplex virus thymidine kinase (TK) gene driven by the OC promoter, was generated to inhibit the growth of osteoblastic osteosarcoma with prodrug acyclovir (ACV). Ad-OC-TK/ACV also inhibited the growth of prostate cancer cells and suppressed the growth of subcutaneous and intraosseous prostate tumor. In order to combine treatment modalities to maximize tumor cell-kill with minimized host toxicities, Ad-OC-TK/ACV was applied in combination with low dose methotrexate to eradicate osteoblastic osteosarcoma. In targeting of micrometastatic disease, intravenous Ad-OC-TK/ACV treatment resulted in significant tumor nodule reduction and prolonged the survival of animals harboring osteosarcoma lung metastases without significant host toxicity. Ad-OC-TK is a rational choice for the treatment of prostate cancer skeletal metastasis because OC is uniformly detected in both primary and metastatic human prostate cancer specimens by immunohistochemistry. Ad-OC-TK/ACV inhibits the growth not only of prostate cancer cells but also of their supporting bone stromal cells. Targeting both prostate cancer epithelium and its supporting stroma may be most efficacious for the treatment of prostate cancer osseous metastases. ^

THE DESIGN AND DEVELOPMENT OF A SURVEY INSTRUMENT FOR THE ASSESSMENT OF OCCUPATIONAL STRESS

Using stress and coping as a unifying theoretical concept, a series of five models was developed in order to synthesize the survey questions and to classify information. These models identified the question, listed the research study, described measurements, listed workplace data, and listed industry and national reference data.^ A set of 38 instrument questions was developed within the five coping correlate categories. In addition, a set of 22 stress symptoms was also developed. The study was conducted within two groups, police and professors, on a large university campus. The groups were selected because their occupations were diverse, but they were a part of the same macroenvironment. The premise was that police officers would be more highly stressed than professors.^ Of a total study group of 80, there were 37 respondents. The difference in the mean stress responses was observable between the two groups. Not only were the responses similar within each group, but the stress level of response was also similar within each group. While the response to the survey instrument was good, only 3 respondents answered the stress symptom survey properly. It was determined that none of the 37 respondents believed that they were ill. This perception of being well was also evidenced by the grand mean of the stress scores of 2.76 (3.0 = moderate stress). This also caused fewer independent variables to be entered in the multiple regression model.^ The survey instrument was carefully designed to be universal. Universality is the ability to transcend occupational or regional definitions as applied to stress. It is the ability to measure responses within broad categories such as physiological, emotional, behavioral, social, and cognitive functions without losing the ability to measure the detail within the individual questions, or the relationships between questions and categories.^ Replication is much easier to achieve with standardized categories, questions, and measurement procedures such as those developed for the universal survey instrument. Because the survey instrument is universal it can be used as an analytical device, an assessment device, a basic tool for planning and a follow-up instrument to measure individual response to planned reductions in occupational stress. (Abstract shortened with permission of author.) ^

An Early-Warning System for Hypo-/Hyperglycemic Events Based on Fusion of Adaptive Prediction Models

Introduction: Early warning of future hypoglycemic and hyperglycemic events can improve the safety of type 1 diabetes mellitus (T1DM) patients. The aim of this study is to design and evaluate a hypoglycemia / hyperglycemia early warning system (EWS) for T1DM patients under sensor-augmented pump (SAP) therapy. Methods: The EWS is based on the combination of data-driven online adaptive prediction models and a warning algorithm. Three modeling approaches have been investigated: (i) autoregressive (ARX) models, (ii) auto-regressive with an output correction module (cARX) models, and (iii) recurrent neural network (RNN) models. The warning algorithm performs postprocessing of the models′ outputs and issues alerts if upcoming hypoglycemic/hyperglycemic events are detected. Fusion of the cARX and RNN models, due to their complementary prediction performances, resulted in the hybrid autoregressive with an output correction module/recurrent neural network (cARN)-based EWS. Results: The EWS was evaluated on 23 T1DM patients under SAP therapy. The ARX-based system achieved hypoglycemic (hyperglycemic) event prediction with median values of accuracy of 100.0% (100.0%), detection time of 10.0 (8.0) min, and daily false alarms of 0.7 (0.5). The respective values for the cARX-based system were 100.0% (100.0%), 17.5 (14.8) min, and 1.5 (1.3) and, for the RNN-based system, were 100.0% (92.0%), 8.4 (7.0) min, and 0.1 (0.2). The hybrid cARN-based EWS presented outperforming results with 100.0% (100.0%) prediction accuracy, detection 16.7 (14.7) min in advance, and 0.8 (0.8) daily false alarms. Conclusion: Combined use of cARX and RNN models for the development of an EWS outperformed the single use of each model, achieving accurate and prompt event prediction with few false alarms, thus providing increased safety and comfort.

Simulating Cortical Development as a Self Constructing Process: A Novel Multi-Scale Approach Combining Molecular and Physical Aspects

Current models of embryological development focus on intracellular processes such as gene expression and protein networks, rather than on the complex relationship between subcellular processes and the collective cellular organization these processes support. We have explored this collective behavior in the context of neocortical development, by modeling the expansion of a small number of progenitor cells into a laminated cortex with layer and cell type specific projections. The developmental process is steered by a formal language analogous to genomic instructions, and takes place in a physically realistic three-dimensional environment. A common genome inserted into individual cells control their individual behaviors, and thereby gives rise to collective developmental sequences in a biologically plausible manner. The simulation begins with a single progenitor cell containing the artificial genome. This progenitor then gives rise through a lineage of offspring to distinct populations of neuronal precursors that migrate to form the cortical laminae. The precursors differentiate by extending dendrites and axons, which reproduce the experimentally determined branching patterns of a number of different neuronal cell types observed in the cat visual cortex. This result is the first comprehensive demonstration of the principles of self-construction whereby the cortical architecture develops. In addition, our model makes several testable predictions concerning cell migration and branching mechanisms.

A Drosophila XPD model links cell cycle coordination with neuro-development and suggests links to cancer

XPD functions in transcription, DNA repair and in cell cycle control. Mutations in human XPD (also known as ERCC2) mainly cause three clinical phenotypes: xeroderma pigmentosum (XP), Cockayne syndrome (XP/CS) and trichothiodystrophy (TTD), and only XP patients have a high predisposition to developing cancer. Hence, we developed a fly model to obtain novel insights into the defects caused by individual hypomorphic alleles identified in human XP-D patients. This model revealed that the mutations that displayed the greatest in vivo UV sensitivity in Drosophila did not correlate with those that led to tumor formation in humans. Immunoprecipitations followed by targeted quantitative MS/MS analysis showed how different xpd mutations affected the formation or stability of different transcription factor IIH (TFIIH) subcomplexes. The XP mutants most clearly linked to high cancer risk, Xpd R683W and R601L, showed a reduced interaction with the core TFIIH and also an abnormal interaction with the Cdk-activating kinase (CAK) complex. Interestingly, these two XP alleles additionally displayed high levels of chromatin loss and free centrosomes during the rapid nuclear division phase of the Drosophila embryo. Finally, the xpd mutations showing defects in the coordination of cell cycle timing during the Drosophila embryonic divisions correlated with those human mutations that cause the neurodevelopmental abnormalities and developmental growth defects observed in XP/CS and TTD patients.

Avian area vasculosa and CAM as rapid in vivo pro-angiogenic and antiangiogenic models.

Angiogenesis, the development of new blood vessels from preexisting ones, is driven by coordinated signaling pathways governed by specific molecules, hemodynamic forces, and endothelial and periendothelial cells. The processes involve adhesion, migration, and survival machinery within the target endothelial and periendothelial cells. Factors that interfere with any of these processes may therefore influence angiogenesis either positively (pro-angiogenesis) or negatively (antiangiogenesis). The avian area vasculosa (AV) and the avian chorioallantoic membrane (CAM) are two useful tools for studying both angiogenesis and antiangiogenesis since they are amenable to both intravascular and topical administration of target, agents, are relatively rapid assays, and can be adapted very easily to study angiogenesis-dependent processes, such as tumor growth. Both models provide a physiological setting that permits investigation of pro-angiogenic and antiangiogenic agent interactions in vivo.

Early development of human lymphomas in a prostate cancer xenograft program using triple knock-out immunocompromised mice

BACKGROUND There is an urgent need for preclinical models of prostate cancer; however, clinically relevant patient-derived prostate cancer xenografts (PDXs) are demanding to establish. METHODS Sixty-seven patients who were undergoing palliative transurethral surgery or radical prostatectomy for histologically confirmed, clinically relevant prostate cancer were included in the study. Fresh prostate cancer tissue was identified by frozen analysis in 48 patients. The cancer tissue was transplanted subcutaneously and under the renal capsule of NSG and NOG mice supplemented with human testosterone. All growing PDXs were evaluated by histology and immunohistochemistry. RESULTS Early assessment of the animals at least three months after transplantation included 27/48 (56.3%) eligible PDX cohorts. PDX growth was detected in 10/27 (37%) mouse cohorts. Eight of the ten PDXs were identified as human donor derived lymphomas, including seven Epstein Barr virus (EBV)-positive diffuse large B-cell lymphomas and one EBV-negative peripheral T-cell lymphoma. One sample consisted of benign prostatic tissue, and one sample comprised a benign epithelial cyst. Prostate cancer was not detected in any of the samples. CONCLUSIONS Tumors that arise within the first three months after prostate cancer xenografting may represent patient-derived EBV-positive lymphomas in up to 80% of the early growing PDXs when using triple knockout NSG immunocompromised mice. Therefore, lymphoma should be excluded in prostate cancer xenografts that do not resemble typical prostatic adenocarcinoma. Prostate 9999: XX-XX, 2014. © 2015 Wiley Periodicals, Inc.

Change in physical activity after smoking cessation: the Coronary Artery Risk Development in Young Adults (CARDIA) study

AIMS To estimate physical activity trajectories for people who quit smoking, and compare them to what would have been expected had smoking continued. DESIGN, SETTING AND PARTICIPANTS A total of 5115 participants in the Coronary Artery Risk Development in Young Adults Study (CARDIA) study, a population-based study of African American and European American people recruited at age 18-30 years in 1985/6 and followed over 25 years. MEASUREMENTS Physical activity was self-reported during clinical examinations at baseline (1985/6) and at years 2, 5, 7, 10, 15, 20 and 25 (2010/11); smoking status was reported each year (at examinations or by telephone, and imputed where missing). We used mixed linear models to estimate trajectories of physical activity under varying smoking conditions, with adjustment for participant characteristics and secular trends. FINDINGS We found significant interactions by race/sex (P = 0.02 for the interaction with cumulative years of smoking), hence we investigated the subgroups separately. Increasing years of smoking were associated with a decline in physical activity in black and white women and black men [e.g. coefficient for 10 years of smoking: -0.14; 95% confidence interval (CI) = -0.20 to -0.07, P < 0.001 for white women]. An increase in physical activity was associated with years since smoking cessation in white men (coefficient 0.06; 95% CI = 0 to 0.13, P = 0.05). The physical activity trajectory for people who quit diverged progressively towards higher physical activity from the expected trajectory had smoking continued. For example, physical activity was 34% higher (95% CI = 18 to 52%; P < 0.001) for white women 10 years after stopping compared with continuing smoking for those 10 years (P = 0.21 for race/sex differences). CONCLUSIONS Smokers who quit have progressively higher levels of physical activity in the years after quitting compared with continuing smokers.

Postnatal Notch1 activation induces T‑cell malignancy in conditional and inducible mouse models

The Notch1 signaling pathway is essential for hematopoietic development. However, the effects of postnatal activation of Notch1 signaling on hematopoietic system is not yet fully understood. We previously generated ZEG‑IC‑Notch1 transgenic mice that have a floxed β‑geo/stop signal between a CMV promoter and intracellular domain of Notch1 (IC‑Notch1). Constitutively active IC‑Notch1 is silent until the introduction of Cre recombinase. In this study, endothelial/hematopoietic specific expression of IC‑Notch1 in double transgenic ZEG‑IC‑Notch1/Tie2‑Cre embryos induced embryonic lethality at E9.5 with defects in vascular system but not in hematopoietic system. Inducible IC‑Notch1 expression in adult mice was achieved by using tetracycline regulated Cre system. The ZEG‑IC‑Notch1/Tie2‑tTA/tet‑O‑Cre triple transgenic mice survived embryonic development when maintained on tetracycline. Post‑natal withdrawal of tetracycline induced expression of IC‑Notch1 transgene in hematopoietic cells of adult mice. The triple transgenic mice displayed extensive T‑cell infiltration in multiple organs and T‑cell malignancy of lymph nodes. In addition, the protein levels of p53 and alternative reading frame (ARF) were decreased in lymphoma‑like neoplasms from the triple transgenic mice while their mRNA expression remained unchanged, suggesting that IC‑Notch1 might repress ARF‑p53 pathway by a post‑transcriptional mechanism. This study demonstrated that activation of constitutive Notch1 signaling after embryonic development alters adult hematopoiesis and induces T‑cell malignancy.

Sample size considerations using mathematical models: an example with Chlamydia trachomatis infection and its sequelae pelvic inflammatory disease.

BACKGROUND The success of an intervention to prevent the complications of an infection is influenced by the natural history of the infection. Assumptions about the temporal relationship between infection and the development of sequelae can affect the predicted effect size of an intervention and the sample size calculation. This study investigates how a mathematical model can be used to inform sample size calculations for a randomised controlled trial (RCT) using the example of Chlamydia trachomatis infection and pelvic inflammatory disease (PID). METHODS We used a compartmental model to imitate the structure of a published RCT. We considered three different processes for the timing of PID development, in relation to the initial C. trachomatis infection: immediate, constant throughout, or at the end of the infectious period. For each process we assumed that, of all women infected, the same fraction would develop PID in the absence of an intervention. We examined two sets of assumptions used to calculate the sample size in a published RCT that investigated the effect of chlamydia screening on PID incidence. We also investigated the influence of the natural history parameters of chlamydia on the required sample size. RESULTS The assumed event rates and effect sizes used for the sample size calculation implicitly determined the temporal relationship between chlamydia infection and PID in the model. Even small changes in the assumed PID incidence and relative risk (RR) led to considerable differences in the hypothesised mechanism of PID development. The RR and the sample size needed per group also depend on the natural history parameters of chlamydia. CONCLUSIONS Mathematical modelling helps to understand the temporal relationship between an infection and its sequelae and can show how uncertainties about natural history parameters affect sample size calculations when planning a RCT.