INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

Development and implementation of a dynamic heterogeneous proton equivalent anthropomorphic thorax phantom for the assessment of scanned proton beam therapy

DEVELOPMENT AND IMPLEMENTATION OF A DYNAMIC HETEROGENEOUS PROTON EQUIVALENT ANTHROPOMORPHIC THORAX PHANTOM FOR THE ASSESSMENT OF SCANNED PROTON BEAM THERAPY by James Leroy Neihart, B.S. APPROVED: ______________________________David Followill, Ph.D. ______________________________Peter Balter, Ph.D. ______________________________Narayan Sahoo, Ph.D. ______________________________Kenneth Hess, Ph.D. ______________________________Paige Summers, M.S. APPROVED: ____________________________ Dean, The University of Texas Graduate School of Biomedical Sciences at Houston DEVELOPMENT AND IMPLEMENTATION OF A DYNAMIC HETEROGENEOUS PROTON EQUIVALENT ANTHROPOMORPHIC THORAX PHANTOM FOR THE ASSESSMENT OF SCANNED PROTON BEAM THERAPY A THESIS Presented to the Faculty of The University of Texas Health Science Center at Houston andThe University of TexasMD Anderson Cancer CenterGraduate School of Biomedical Sciences in Partial Fulfillment of the Requirements for the Degree of MASTER OF SCIENCE by James Leroy Neihart, B.S. Houston, Texas Date of Graduation August, 2013 Acknowledgments I would like to acknowledge my advisory committee members, chair David Followill, Ph.D., Peter Balter, Ph.D, Narayan Sahoo, Ph.D., Kenneth Hess, Ph.D., Paige Summers M.S. and, for their time and effort contributed to this project. I would additionally like to thank the faculty and staff at the PTC-H and the RPC who assisted in many aspects of this project. Falk Pӧnisch, Ph.D. for his breath hold proton therapy treatment expertise, Matt Palmer and Jaques Bluett for proton dosimetry assistance, Matt Kerr for verification plan assistance, Carrie Amador, Nadia Hernandez, Trang Nguyen, Andrea Molineu, Lynda McDonald for TLD and film dosimetry assistance. Finally, I would like to thank my wife and family for their support and encouragement during my research and studies. Development and implementation of a dynamic heterogeneous proton equivalent anthropomorphic thorax phantom for the assessment of scanned proton beam therapy By: James Leroy Neihart, B.S. Chair of Advisory Committee: David Followill, Ph.D Proton therapy has been gaining ground recently in radiation oncology. To date, the most successful utilization of proton therapy is in head and neck cases as well as prostate cases. These tumor locations do not suffer from the resulting difficulties of treatment delivery as a result of respiratory motion. Lung tumors require either breath hold or motion tracking, neither of which have been assessed with an end-to-end phantom for proton treatments. Currently, the RPC does not have a dynamic thoracic phantom for proton therapy procedure assessment. Additionally, such a phantom could be an excellent means of assessing quality assurance of the procedures of proton therapy centers wishing to participate in clinical trials. An eventual goal of this phantom is to have a means of evaluating and auditing institutions for the ability to start clinical trials utilizing proton therapy procedures for lung cancers. Therefore, the hypothesis of this study is that a dynamic anthropomorphic thoracic phantom can be created to evaluate end-to-end proton therapy treatment procedures for lung cancer to assure agreement between the measured and calculated dose within 5% / 5 mm with a reproducibility of 2%. Multiple materials were assessed for thoracic heterogeneity equivalency. The phantom was designed from the materials found to be in greatest agreement. The phantom was treated in an end-to-end treatment four times, which included simulation, treatment planning and treatment delivery. Each treatment plan was delivered three times to assess reproducibility. The dose measured within the phantom was compared to that of the treatment plan. The hypothesis was fully supported for three of the treatment plans, but failed the reproducibility requirement for the most aggressive treatment plan.

Baseline quality of life as a predictor of exercise intervention adherence among endometrial cancer survivors

Background: The objective of this analysis is to test whether baseline quality of life (QOL) measurements, body mass index (BMI) and prior exercise behavior are significantly associated with (1) telephone counseling adherence, and (2) activity at the final assessment, in a physical activity promoting intervention among endometrial cancer survivors.^ Methods: One hundred endometrial cancer survivors not currently meeting physical activity guidelines completed baseline QOL and anthropometric assessments to measure general physical and mental health [Medical Outcomes Survey (SF-36)], sleep patterns and sleep quality [Pittsburgh Sleep Quality Index (PSQI)], perceived stress [Perceived Stress Scale (PSS)], cancer-specific concerns of long-term survivors [Quality of Life in Adult Cancer Survivors (QLACS)], and psychological distress [Brief Symptom Inventory-18 (BSI-18)]. Survivors were counseled by telephone during the 6-month intervention and their completion rate determined their adherence. The primary variables of interest included age, baseline BMI, baseline activity level, time since diagnosis, education, treatment received, and the SF-36 physical and mental component scores.^ Results: Final activity was most closely linked with baseline activity (p<.001) and less invasive surgery, being leaner and older, and experiencing less pain and more vitality. Telephone counseling was also predicted well by baseline activity, working less and having better overall cancer-related functioning.^ Conclusion: Above and beyond the QOL measures, baseline activity was the strongest predictor of both final activity and telephone counseling adherence. While education, surgery treatment type and bodily pain were important predictors for final exercise and employment status and cancer-related quality of life were important predictors for telephone counseling adherence, considering adaptive exercise interventions that focus heavily on engaging inactive participants may be a way to produce better exercise-related outcomes in the endometrial cancer survivor population.^

The risk prediction for SPM and its impact on the survival in patients with head and neck squamous cell carcinoma

Head and Neck Squamous Cell Carcinoma (HNSCC) is the sixth common malignancy in the world, with high rates of developing second primary malignancy (SPM) and moderately low survival rates. This disease has become an enormous challenge in the cancer research and treatments. For HNSCC patients, a highly significant cause of post-treatment mortality and morbidity is the development of SPM. Hence, assessment of predicting the risk for the development of SPM would be very helpful for patients, clinicians and policy makers to estimate the survival of patients with HNSCC. In this study, we built a prognostic model to predict the risk of developing SPM in patients with newly diagnosed HNSCC. The dataset used in this research was obtained from The University of Texas MD Anderson Cancer Center. For the first aim, we used stepwise logistic regression to identify the prognostic factors for the development of SPM. Our final model contained cancer site and overall cancer stage as our risk factors for SPM. The Hosmer-Lemeshow test (p-value= 0.15>0.05) showed the final prognostic model fit the data well. The area under the ROC curve was 0.72 that suggested the discrimination ability of our model was acceptable. The internal validation confirmed the prognostic model was a good fit and the final prognostic model would not over optimistically predict the risk of SPM. This model needs external validation by using large data sample size before it can be generalized to predict SPM risk for other HNSCC patients. For the second aim, we utilized a multistate survival analysis approach to estimate the probability of death for HNSCC patients taking into consideration of the possibility of SPM. Patients without SPM were associated with longer survival. These findings suggest that the development of SPM could be a predictor of survival rates among the patients with HNSCC.^

Microbiota and Microsatellite genotypes of Pacific oysters stemming from three oyster bed in the Northern Wadden Sea

Background: Studies of oyster microbiomes have revealed that a limited number of microbes, including pathogens, can dominate microbial communities in host tissues such as gills and gut. Much of the bacterial diversity however remains underexplored and unexplained, although environmental conditions and host genetics have been implicated. We used 454 next generation 16S rRNA amplicon sequencing of individually tagged PCR reactions to explore the diversity of bacterial communities in gill tissue of the invasive Pacific oyster Crassostrea gigas stemming from genetically differentiated beds under ambient outdoor conditions and after a multifaceted disturbance treatment imposing stress on the host. Results: While the gill associated microbial communities in oysters were dominated by few abundant taxa (i.e. Sphingomonas, Mycoplasma) the distribution of rare bacterial groups correlated to relatedness between the hosts under ambient conditions. Exposing the host to disturbance broke apart this relationship by removing rare phylotypes thereby reducing overall microbial diversity. Shifts in the microbiome composition in response to stress did not result in a net increase in genera known to contain potentially pathogenic strains. Conclusion: The decrease in microbial diversity and the disassociation between population genetic structure of the hosts and their associated microbiome suggest that disturbance (i.e. stress) may play a significant role for the assembly of the natural microbiome. Such community shifts may in turn also feed back on the course of disease and the occurrence of mass mortality events in oyster populations.

Carbon and Nitrogen concentrations and isotopic composition of ODP Site 170-1040 and Leg 205 sites

We determined the C and N concentrations and isotopic compositions of sediments in the prism sampled during Ocean Drilling Program Legs 170 and 205 offshore Costa Rica, with the goals of evaluating sediment sources and extents of diagenesis and identifying any effects of infiltrating fluids on the sedimentary C and N. The sediments from Leg 170 Site 1040 contain 0.85-1.96 wt% total organic carbon (TOC) with Vienna Peedee belemnite (VPDB) d13CVPDB from -26.3 per mil to -22.5 per mil, and 832-2221 ppm total nitrogen (TN) with d15Nair from +3.5 per mil to +6.6 per mil. Sediment TN concentrations and d15N values show dramatic downhole increases within the uppermost 130 m of the section and more gradual downhole decreases from 130 meters below seafloor (mbsf) to the base of the décollement at ~370 mbsf. Concentrations and isotopic compositions of TOC are relatively uniform within the entire section, showing some minor perturbation within the décollement zone. In the uppermost 100 m, upsection increases in TN concentrations at constant TOC concentrations produce significant increases in atomic TOC/TN ratios from ~8 to ~18. Carbonate (calcite) contents in the wedge sediments are generally low (<4 wt%), but the d13C and Vienna standard mean ocean water (VSMOW) d18OVSMOW values vary significantly from -26.1 per mil to +4.1 per mil and from +30.0 per mil to +35.3 per mil, respectively. Concentrations and isotopic compositions of TOC and TN for sediments from Leg 205 Sites 1254 and 1255 overlap well with C-N data for sediments from the same depth intervals obtained during Leg 170 at Site 1040.

Sedimentology and stable oxygen isotope record of the northern South China Sea

High-resolution sediment records from the South China Sea reveal a winter monsoon dominated glacial regime and a summer monsoon dominated Holocene regime during the last glacial cycle. A fundamental change between regimes occurred during deglaciation through a series of millennial reoccurrences of century-scale changes in the East Asian monsoon (EAM) climate. These abrupt events centered at 17.0, 15.9, 15.5, 14.7, 13.5, 13.9, 13.3, 12.1, 11.5, and 10.7 14C ka correlate well with the millennial-scale events in the Santa Barbara Basin and the Arabian Sea, i.e. a relationship between EAM and El Niño/Southern Oscillation systems. The abrupt increases in summer monsoon imply enhanced heat transport from low-latitude sea area to the midlatitude/high-latitude land area. The phase relationship between events of EAM and ice sheet may reflect a faster EAM response and a slower ice sheet response to the insolation change. A far-reaching conclusion is that the EAM might have triggered the Northern Hemisphere deglaciation.

CO2 and inorganic nutrient enrichment affect the performance of a calcifying green alga and its noncalcifying epiphyte

Ocean acidification studies in the past decade have greatly improved our knowledge of how calcifying organisms respond to increased surface ocean CO2 levels. It has become evident that, for many organisms, nutrient availability is an important factor that influences their physiological responses and competitive interactions with other species. Therefore, we tested how simulated ocean acidification and eutrophication (nitrate and phosphate enrichment) interact to affect the physiology and ecology of a calcifying chlorophyte macroalga (Halimeda opuntia (L.) J.V. Lamouroux) and its common noncalcifying epiphyte (Dictyota sp.) in a 4-week fully crossed multifactorial experiment. Inorganic nutrient enrichment (+NP) had a strong influence on all responses measured with the exception of net calcification. Elevated CO2 alone significantly decreased electron transport rates of the photosynthetic apparatus and resulted in phosphorus limitation in both species, but had no effect on oxygen production or respiration. The combination of CO2 and +NP significantly increased electron transport rates in both species. While +NP alone stimulated H. opuntia growth rates, Dictyota growth was significantly stimulated by nutrient enrichment only at elevated CO2, which led to the highest biomass ratios of Dictyota to Halimeda. Our results suggest that inorganic nutrient enrichment alone stimulates several aspects of H. opuntia physiology, but nutrient enrichment at a CO2 concentration predicted for the end of the century benefits Dictyota sp. and hinders its calcifying basibiont H. opuntia.

Biogenic and environmental fluctuations of temperature and pH in a novel mesocosm concept ("Kiel Outdoor Benthocosms") April 2013 to September 2014

The Kiel Outdoor Benthocosm infrastructure (Kiel, Germany,N 54°19.8'; E 010°09.0') allows combining natural in-situ fluctuations on all environmental variables with the controlled manipulation of treatment factors. The environmental fluctuations are admitted by a continuous flow-through of fjord water. The treatment is applied by delta-treatments which shift the mean of target variables (temperature and pH in this case) while maintaining the frequency and amplitude of natural fluctuations. The data presented here show the treatment levels and the continuously logged temperature and pH conditions in the experimental tanks. The dynamics of temperature and pH are driven by (i) in situ variability, (ii) the treatments imposed and (iii) the biology of the biota in the tanks. These contained macroalgal communities with associated mesograzers, mussels, and sea stars. The data set comprised 5 experimental runs: spring experiment (4.4.-19.6.2013), summer experiment 1 (4.7.-17.9.2013), autumn experiment (10.10-17.12.2013), winter experiment (16.1. - 1.4.2014), summer experiment 2 (10.7. - 26.9.2014).