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.

A Model for Family Preservation Case Assessment

The outcomes of family preservation practice have been researched and debated. The effectiveness of family preservation is still inconclusive and many of the findings may only be inferred to specific situations. Few studies have addressed the assessment techniques or outcome factors from a qualitative perspective. This article synthesizes current literature, research and practice, and proposes a practice framework with questioning techniques to assist practitioners in assessing the strengths and characteristics of a family, and making decisions on whether or not familybased services are appropriate for the family. Two actual cases are presented to illustrate how the worker can benefit from having the assessment data derived from this model.

THE EXPRESSION PROFILE OF A GWAS DETECTED GENE (NINJ2) IN THE BRAIN OF A MOUSE MODEL OF ISCHEMIC STROKE

Stroke is the third leading cause of death and a major debilitating disease in the United States. Multiple factors, including genetic factors, contribute to the development of the disease. Genome-wide association studies (GWAS) have contributed to the identification of genetic loci influencing risk for complex diseases, such as stroke. In 2010, a GWAS of incident stroke was performed in four large prospective cohorts from the USA and Europe and identified an association of two Single Nucleotide Polymorphisms (SNPs) on chromosome 12p13 with a greater risk of ischemic stroke in individuals of European and African-American ancestry. These SNPs are located 11 Kb upstream of the nerve injury-induced gene 2, Ninjurin2 (NINJ2), suggesting that this gene may be involved in stroke pathogenesis. NINJ2 is a cell adhesion molecule induced in the distal glial cells from a sciatic-nerve injury at 7-days after injury. In an effort to ascribe a possible role of NINJ2 in stroke, we have assessed changes in the level of gene and protein expression of NINJ2 following a time-course from a transiently induced middle cerebral artery ischemic stroke in mice brains. We report an increase in the gene expression of NINJ2 in the ischemic and peri-infarct (ipsilateral) cortical tissues at 7 and 14-days after stroke. We also report an increase in the protein expression of NINJ2 in the cortex of both the ipsilateral and contralateral cortical tissues at the same time-points. We conclude that the expression of NINJ2 is regulated by an ischemic stroke in the cortex and is consistent with NINJ2 being involved in the recovery time-points of stroke.

Detection and sequence analysis of a nickel-induced mutation in a retroviral model system

We have developed a novel way to assess the mutagenicity of environmentally important metal carcinogens, such as nickel, by creating a positive selection system based upon the conditional expression of a retroviral transforming gene. The target gene is the v-mos gene in MuSVts110, a murine retrovirus possessing a growth temperature dependent defect in expression of the transforming gene due to viral RNA splicing. In normal rat kidney cells infected with MuSVts110 (6m2 cells), splicing of the MuSVts110 RNA to form the mRNA from which the transforming protein, p85$\sp{\rm gag-mos}$, is translated is growth-temperature dependent, occurring at 33 C and below but not at 39 C and above. This splicing "defect" is mediated by cis-acting viral sequences. Nickel chloride treatment of 6m2 cells followed by growth at 39 C, allowed the selection of "revertant" cells which constitutively express p85$\sp{\rm gag-mos}$ due to stable changes in the viral RNA splicing phenotype, suggesting that nickel, a carcinogen whose mutagenicity has not been well established, could induce mutations in mammalian genes. We also show by direct sequencing of PCR-amplified integrated MuSVts110 DNA from a 6m2 nickel-revertant cell line that the nickel-induced mutation affecting the splicing phenotype is a cis-acting 70-base duplication of a region of the viral DNA surrounding the 3$\sp\prime$ splice site. These findings provide the first example of the molecular basis for a nickel-induced DNA lesion and establish the mutagenicity of this potent carcinogen. ^

Differential gene expression in multistep tumorigenesis using an {\it in vitro\/} human cell model

Using a human terato-carcinoma cell line, PA-1, the functional role of the oncogenes and tumor suppressor gene involved in the multistep process of carcinogenesis have been analyzed. The expression of AP-2 was strongly correlated with the susceptibility to ras transformation. The differential responsiveness to growth factors between stage 1 ras resistant cells and stage 2 ras susceptible cells was observed, indicating that the ability of stage 2 cells to respond to the mutated ras oncogenes in transformation correlated with the ability to be stimulated by certain growth factors. Using differential screening of cDNA libraries, a number of differentially expressed cDNA clones was isolated. One of those, clone 12, is overexpressed in ras transformed stage 3 cells. The amino acid sequence of clone 12 is almost identical to a mouse LLrep3 gene that was growth-regulated, and 78% similar to a yeast ribosomal protein S4. These results suggest that the S4 gene may be involved in regulation of growth. Clone 9 is expressed in stage 1 ras resistant cells (3.5-kb and 3.0-kb transcripts) but the expression of this clone in stage 2 ras susceptible cells and stage 3 ras-transformed cells is greatly diminished. The expression of this cDNA clone was increased to at least five fold in ras resistant cells and nontumorigenic hybrids treated with retinoic acid but not increased in retinoic acid treated ras susceptible cells, ras transformed cells and the tumorigenic segregants. Partial sequence of this clone showed no homology to the sequences in Genbank. These findings suggest that clone 9 could be a suppressor gene or the genes that are involved in the biochemical pathway of tumor suppression or neurogenic differentiation. The apparent pleiotropic effect of the loss of this suppressor gene function support Harris' proposal that tumor suppressor genes regulate differentiation. The tumor suppressor gene may act as negative regulator of tumor growth by controlling gene expression in differentiation. ^

A convolution model for energy transport in a therapeutic fast neutron beam

A three-dimensional model has been proposed that uses Monte Carlo and fast Fourier transform convolution techniques to calculate the dose distribution from a fast neutron beam. This method transports scattered neutrons and photons in the forward, lateral, and backward directions and protons, electrons, and positrons in the forward and lateral directions by convolving energy spread kernels with initial interaction available energy distributions. The primary neutron and photon spectrums have been derived from narrow beam attenuation measurements. The positions and strengths of the effective primary neutron, scattered neutron, and photon sources have been derived from dual ion chamber measurements. The size of the effective primary neutron source has been measured using a copper activation technique. Heterogeneous tissue calculations require a weighted sum of two convolutions for each component since the kernels must be invariant for FFT convolution. Comparisons between calculations and measurements were performed for several water and heterogeneous phantom geometries. ^

Quantitative MRI of spinal cord injury in a rat model: Correlative studies

Serial quantitative and correlative studies of experimental spinal cord injury (SCI) in rats were conducted using three-dimensional magnetic resonance imaging (MRI). Correlative measures included morphological histopathology, neurobehavioral measures of functional deficit, and biochemical assays for N-acetyl-aspartate (NAA), lactate, pyruvate, and ATP. A spinal cord injury device was characterized and provided a reproducible injury severity. Injuries were moderate and consistent to within $\pm$20% (standard deviation). For MRI, a three-dimensional implementation of the single spin-echo FATE (Fast optimum angle, short TE) pulse sequence was used for rapid acquisition, with a 128 x 128 x 32 (x,y,z) matrix size and a 0.21 x 0.21 x 1.5 mm resolution. These serial studies revealed a bimodal characteristic in the evolution in MRI pathology with time. Early and late phases of SCI pathology were clearly visualized in $T\sb2$-weighted MRI, and these corresponded to specific histopathological changes in the spinal cord. Centralized hypointense MRI regions correlated with evidence of hemorrhagic and necrotic tissue, while surrounding hyperintense regions represented edema or myelomalacia. Unexpectedly, $T\sb2$-weighted MRI pathology contrast at 24 hours after injury appeared to subside before peaking at 72 hours after injury. This change is likely attributable to ongoing secondary injury processes, which may alter local $T\sb2$ values or reduce the natural anisotropy of the spinal cord. MRI, functional, and histological measures all indicated that 72 hours after injury was the temporal maximum for quantitative measures of spinal cord pathology. Thereafter, significant improvement was seen only in neurobehavioral scores. Significant correlations were found between quantitated MRI pathology and histopathology. Also, NAA and lactate levels correlated with behavioral measures of the level of function deficit. Asymmetric (rostral/caudal) changes in NAA and lactate due to injury indicate that rostral and caudal segments from the injury site are affected differently by the injury. These studies indicate that volumetric quantitation of MRI pathology from $T\sb2$-weighted images may play an important role in early prediction of neurologic deficit and spinal cord pathology. The loss of $T\sb2$ contrast at 24 hours suggests MR may be able to detect certain delayed mechanisms of secondary injury which are not resolved by histopathology or other radiological modalities. Furthermore, in vivo proton magnetic resonance spectroscopy (MRS) studies of SCI may provide a valuable addition source of information about changes in regional spinal cord lactate and NAA levels, which are indicative of local metabolic and pathological changes. ^

Involvement of heparin-like glycosaminoglycans and expression of a novel heparan sulfate binding protein (p24) during human placentation and in a model for human implantation

Previous studies in our laboratory have indicated that heparan sulfate proteoglycans (HSPGs) play an important role in murine embryo implantation. To investigate the potential function of HSPGs in human implantation, two human cell lines (RL95 and JAR) were selected to model uterine epithelium and embryonal trophectoderm, respectively. A heterologous cell-cell adhesion assay showed that initial binding between JAR and RL95 cells is mediated by cell surface glycosaminoglycans (GAG) with heparin-like properties, i.e., heparan sulfate and dermatan sulfate. Furthermore, a single class of highly specific, protease-sensitive heparin/heparan sulfate binding sites exist on the surface of RL95 cells. Three heparin binding, tryptic peptide fragments were isolated from RL95 cell surfaces and their amino termini partially sequenced. Reverse transcription-polymerase chain reaction (RT-PCR) generated 1 to 4 PCR products per tryptic peptide. Northern blot analysis of RNA from RL95 cells using one of these RT-PCR products identified a 1.2 Kb mRNA species (p24). The amino acid sequence predicted from the cDNA sequence contains a putative heparin-binding domain. A synthetic peptide representing this putative heparin binding domain was used to generate a rabbit polyclonal antibody (anti-p24). Indirect immunofluorescence studies on RL95 and JAR cells as well as binding studies of anti-p24 to intact RL95 cells demonstrate that p24 is distributed on the cell surface. Western blots of RL95 membrane preparations identify a 24 kDa protein (p24) highly enriched in the 100,000 g pellet plasma membrane-enriched fraction. p24 eluted from membranes with 0.8 M NaCl, but not 0.6 M NaCl, suggesting that it is a peripheral membrane component. Solubilized p24 binds heparin by heparin affinity chromatography and $\sp{125}$I-heparin binding assays. Furthermore, indirect immunofluorescence studies indicate that cytotrophoblast of floating and attached villi of the human fetal-maternal interface are recognized by anti-p24. The study also indicates that the HSPG, perlecan, accumulates where chorionic villi are attached to uterine stroma and where p24-expressing cytotrophoblast penetrate the stroma. Collectively, these data indicate that p24 is a cell surface membrane-associated heparin/heparan sulfate binding protein found in cytotrophoblast, but not many other cell types of the fetal-maternal interface. Furthermore, p24 colocalizes with HSPGs in regions of cytotrophoblast invasion. These observations are consistent with a role for HSPGs and HSPG binding proteins in human trophoblast-uterine cell interactions. ^

Genetic analysis of tumor progression susceptibility in the mouse skin model

In the field of chemical carcinogenesis the use of animal models has proved to be a useful tool in dissecting the multistage process of tumor formation. In this regard the outbred SENCAR mouse has been the strain of choice in the analysis of skin carcinogenesis given its high sensitivity to the chemically induced acquisition of premalignant lesions, papillomas, and the later progression of these lesions into squamous cell carcinomas (SCC).^ The derivation of an inbred strain from the SENCAR stock called SSIN, that in spite of a high sensitivity to the development of papillomas lack the ability to transform these premalignant lesions into SCC, suggested that tumor promotion and progression were under the genetic control of different sets of genes.^ In the present study the nature of susceptibility to tumor progression was investigated. Analysis of F1 hybrids between the outbred SENCAR and SSIN mice suggested that there is at least one dominant gene responsible for susceptibility to tumor progression.^ Later development of another inbred strain from the outbred SENCAR stock, that had sensitivity to both tumor promotion and progression, allowed the formulation of a more accurate genetic model. Using this newly derived line, SENCAR B/Pt. and SSIN it was determined that there is one dominant tumor progression susceptibility gene. Linkage analysis showed that this gene maps to mouse chromosome 14 and it was possible to narrow the region to a 16 cM interval.^ In order to better characterize the nature of the progression susceptibility differences between these two strains, their proliferative pattern was investigated. It was found that SENCAR B/Pt, have an enlarged proliferative compartment with overexpression of cyclin D1, p16 and p21. Further studies showed an aberrant overexpression of TGF-$\beta$ in the susceptible strain, an increase in apoptosis, p53 protein accumulation and early loss of connexin 26. These results taken together suggest that papillomas in the SENCAR B/Pt. mice have higher proliferation and may have an increase in genomic instability, these two factors would contribute to a higher sensitivity to tumor progression. ^

A LAPLACE TRANSFORM PAIR MODEL TO DETERMINE BREMSSTRAHLUNG SPECTRA FROM ATTENUATION DATA (RADIATION, X-RAY, SPECTRAL MODELING)

Every x-ray attenuation curve inherently contains all the information necessary to extract the complete energy spectrum of a beam. To date, attempts to obtain accurate spectral information from attenuation data have been inadequate.^ This investigation presents a mathematical pair model, grounded in physical reality by the Laplace Transformation, to describe the attenuation of a photon beam and the corresponding bremsstrahlung spectral distribution. In addition the Laplace model has been mathematically extended to include characteristic radiation in a physically meaningful way. A method to determine the fraction of characteristic radiation in any diagnostic x-ray beam was introduced for use with the extended model.^ This work has examined the reconstructive capability of the Laplace pair model for a photon beam range of from 50 kVp to 25 MV, using both theoretical and experimental methods.^ In the diagnostic region, excellent agreement between a wide variety of experimental spectra and those reconstructed with the Laplace model was obtained when the atomic composition of the attenuators was accurately known. The model successfully reproduced a 2 MV spectrum but demonstrated difficulty in accurately reconstructing orthovoltage and 6 MV spectra. The 25 MV spectrum was successfully reconstructed although poor agreement with the spectrum obtained by Levy was found.^ The analysis of errors, performed with diagnostic energy data, demonstrated the relative insensitivity of the model to typical experimental errors and confirmed that the model can be successfully used to theoretically derive accurate spectral information from experimental attenuation data. ^

AUGMENTATION OF RECONSTITUTION OF GUT ASSOCIATED LYMPHOID TISSUE IN THE SYNGENEIC MURINE BONE MARROW TRANSPLANTATION MODEL

Gut was studied as a prototypical mucosal membrane in the murine BDF-1 syngeneic bone marrow transplant model. Measures of jejunal intraepithelial lymphocytes (IELs) and crypt cells were obtained by standard techniques and a method of quantifying gut lamina propria plasma cells (PCs) was developed. The degree of ablation of gut PCs and IELs after 900 rads total body irradiation with ('60)Co, and their repopulation effected by transplantation with 2.0 x 10('5) or 1.0 x 10('6) bone marrow cells demonstrated a prolonged period of profound depression in population levels of these cells which was not reflected by the extent of damage sustained to the epithelium. Differences in the depopulation and recovery of gut PCs and IELs revealed a tendency towards initial differentiation of effector cells. A positive dose response to high bone marrow cell innocula was obtained. Subsequent studies determined that gut IEL and PC repopulation was potentiated by the addition of IELs or buffy coat cells (BCs) to the bone marrow transplant. A method of isolating 1.4 - 4.0 x 10('7) viable IELs per gram of murine small bowel was devised employing intralumenal hyaluronidase digestion of the epithelial layer and centrifugation of the resulting suspension through discontinuous Percoll gradients. Irradiated mice received 2.0 x 10('5) bone marrow cells along with an equal number of IELs or BCs. The extent and duration of depression of numbers of IELs and PCs was markedly reduced by the addition of the IEL isolate to the transplantation innocula, and to a lesser degree by the addition of BCs. The augmentation of repopuation far exceeded that expected by simple lodging of cells suggesting that the additionally transplanted cells contained a subpopulation of mucosal membrane lymphoid stem cells or helper cells. Correlation analysis of PC versus IEL levels suggests a possible feedback mechanism governing the relative size of their populations. Normal ratios of IgA, IgM, and IgG bearing PCs was maintained post transplantation with all of the regimens. ^

Enhanced in vivo angiogenesis within a model tissue engineered construct using biodegradable microspheres containing encapsulated VEGF

Introduction. Tissue engineering techniques offer a potential means to develop a tissue engineered construct (TEC) for the treatment of tissue and organ deficiencies. However, a lack of adequate vascularization is a limiting factor in the development of most viable engineered tissues. Vascular endothelial growth factor (VEGF) could aid in the development of a viable vascular network within TECs. The long-term goals of this research are to develop clinically relevant, appropriately vascularized TECs for use in humans. This project tested the hypothesis that the delivery of VEGF via controlled release from biodegradable microspheres would increase the vascular density and rate of angiogenesis within a model TEC. ^ Materials and methods. Biodegradable VEGF-encapsulated microspheres were manufactured using a novel method entitled the Solid Encapsulation/Single Emulsion/Solvent Extraction technique. Using a PLGA/PEG polymer blend, microspheres were manufactured and characterized in vitro. A model TEC using fibrin was designed for in vivo tissue engineering experimentation. At the appropriate timepoint, the TECs were explanted, and stained and quantified for CD31 using a novel semi-automated thresholding technique. ^ Results. In vitro results show the microspheres could be manufactured, stored, degrade, and release biologically active VEGF. The in vivo investigations revealed that skeletal muscle was the optimal implantation site as compared to dermis. In addition, the TECs containing fibrin with VEGF demonstrated significantly more angiogenesis than the controls. The TECs containing VEGF microspheres displayed a significant increase in vascular density by day 10. Furthermore, TECs containing VEGF microspheres had a significantly increased relative rate of angiogenesis from implantation day 5 to day 10. ^ Conclusions. A novel technique for producing microspheres loaded with biologically active proteins was developed. A defined concentration of microspheres can deliver a quantifiable level of VEGF with known release kinetics. A novel model TEC for in vivo tissue engineering investigations was developed. VEGF and VEGF microspheres stimulate angiogenesis within the model TEC. This investigation determined that biodegradable rhVEGF 165-encapsulated microspheres increased the vascular density and relative rate of angiogenesis within a model TEC. Future applications could include the incorporation of microvascular fragments into the model TEC and the incorporation of specific tissues, such as fat or bone. ^

CATCH across Texas: Building a model for the diffusion of effective school-based health promotion programs

The impact of health promotion programs is related to both program effectiveness and the extent to which the program is implemented among the target population. The purpose of this dissertation was to describe the development and evaluation of a school-based program diffusion intervention designed to increase the rate of dissemination and adoption of the Child and Adolescent Trial for Cardiovascular Health, or CATCH program (recently renamed the Coordinated Approach to Child Health). ^ The first study described the process by which schools across the state of Texas spontaneously began to adopt the CATCH program after it was tested and proven effective in a multi-site randomized efficacy trial. A survey of teachers and administrator representatives of all schools on record that purchased the CATCH program, but were not involved in the efficacy trial, was used to find out who brought CATCH into the schools, how they garnered support for its adoption, why they decided to adopt the program, and what was involved in deciding to adopt. ^ The second study described how the Intervention Mapping framework guided the planning, development and implementation of a program for the diffusion of CATCH. An iterative process was used to integrate theory, literature, the experience of project staff and data from the target population into a meaningful set of program determinants and performance objectives. Proximal program objectives were specified and translated into both media and interpersonal communication strategies for program diffusion. ^ The third study assessed the effectiveness of the diffusion program in a case-comparison design. Three of the twenty Education Service Center regions in Texas were chosen, selected based on similar demographic criteria, and were followed for adoption of the CATCH curriculum. One of these regions received the full media and interpersonal channel intervention; a second received a reduced media-only intervention, and a third received no intervention. Results suggested the use of the interpersonal channels with media follow-up is an effective means to facilitate program dissemination and adoption. The media-alone condition was not effective in facilitating program adoption. ^

Mixed effects model in negative exponential curve

Despite many researches on development in education and psychology, not often is the methodology tested with real data. A major barrier to test the growth model is that the design of study includes repeated observations and the nature of the growth is nonlinear. The repeat measurements on a nonlinear model require sophisticated statistical methods. In this study, we present mixed effects model in a negative exponential curve to describe the development of children's reading skills. This model can describe the nature of the growth on children's reading skills and account for intra-individual and inter-individual variation. We also apply simple techniques including cross-validation, regression, and graphical methods to determine the most appropriate curve for data, to find efficient initial values of parameters, and to select potential covariates. We illustrate with an example that motivated this research: a longitudinal study of academic skills from grade 1 to grade 12 in Connecticut public schools. ^