Cell surface $\beta$-1,4-galactosyltransferase function during invasion and metastasis

Metastasis is the complex process of tumor cell spread which is responsible for the majority of cancer-related deaths. Metastasis necessitates complex phenotypic changes, many of which are mediated by changes in the activities of cell surface molecules. One of these is cell surface $\beta$1,4-galactosyltransferase (GalTase), which is elevated on more highly metastatic cells. In this study, both molecular and biochemical methods were used to perturb and manipulate cell surface GalTase levels on K1735 murine melanoma cell lines in order to examine its function in metastasis.^ As expected, highly metastatic K1735 variants have higher cell surface GalTase than poorly metastatic variants. Stably transfected K1735 cell lines that overexpress surface GalTase were created. These cell lines were assayed for metastatic ability using an invasion chamber with Matrigel-coated filter inserts. Cells with increased surface GalTase were uniformly more invasive than neo transfected controls. With multiple cell lines, there was a direct correlation (r = 0.918) between surface GalTase activity and invasiveness. Homologous recombination was used to create K1735 cells with decreased levels of surface GalTase. These cells were uniformly less invasive than controls. Cell surface GalTase was inhibited using two different biochemical strategies. In both cases, inhibition of surface GalTase led to a decrease in in vivo metastatic ability of K1735 cells. This is the first direct experimental evidence addressing GalTase function in metastasis. These data provide several lines of independent evidence which show that cell surface GalTase facilitates invasion and metastasis. ^

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Middle managers' attitudes toward employee assistance programs in the corporate setting

Several studies have shown that successful Employee Assistance Programs (EAPs) have strong management endorsement. Strong management endorsement is defined as positive support in utilizing EAP services for themselves and their employees. This study focuses solely on middle management as opposed to upper or general management support. The study further examines success or lack of success of an EAP by the utilization rate defined as the number of employees over a year period who access EAP services.^ A analytical cross-sectional design was used to compare and observe differences between two groups of middle managers (utilizers and nonutilizers). Middle manager data was collected through a mail questionnaire. The study focused on identifying predictors that influence middle managers' utilization rate specifically: attitude toward EAPs, EAP knowledge level, attitude toward mental health professionals, age, gender, years worked as a middle manager, education level, training, and other possible predictors of utilization. The overall hypothesis states middle manager utilizers of EAP services have more positive attitudes and a better understanding of their EAP than middle management nonutilizers.^ As predicted, nonparametric bivariate results showed significant differences between the two groups. Middle managers in the utilization group (n = 473) tended to show more positive attitudes toward their EAP and mental health professionals and demonstrated greater EAP knowledge compared to the nonutilization group (n = 154). These findings support past studies on variables that influence EAP utilization rates.^ Further variables found to influence middle management utilization were identified by multivariate logistic regression results. These variable were gender (female supervisors), educational levels of employees supervised (employees with lower levels of education), number of employees supervised (greater the number supervised, more likely to utilize), managerial EAP training (trained supervisors) and awareness that problems do influence an employee's productivity.^ These findings strengthen the assertion that middle management's attitudes, as well as other variables may influence utilization. Study findings add new information about important variables specifically influencing middle management who utilize EAPs. An understanding of these variables is essential in developing competent EAP program training and orientation programs for middle managers. ^

Phototaxis signaling by the membrane complex of archaeal sensory rhodopsins and their transducers

Sensory rhodopsins I and II (SRI and SRII) are visual pigment-like phototaxis receptors in the archaeon Halobacterium salinarum. The receptor proteins each consist of a single polypeptide that folds into 7 $\alpha$-helical membrane-spanning segments forming an internal pocket where the chromophore retinal is bound. They transmit signals to their tightly bound transducer proteins, HtrI and HtrII, respectively, which in turn control a phosphotransfer pathway modulating the flagellar motors. SRI-HtrI mediates attractant responses to orange-light and repellent responses to UV light, while SRII-HtrII mediates repellent response to blue light. Experiments were designed to analyze the molecular processes in the SR-Htr complexes responsible for receptor activation, which previously had been shown by our laboratory to involve proton transfer reactions of the retinylidene Schiff base in the photoactive site, transfer of signals from receptor to transducer, and signaling specificity by the receptor-transducer complex.^ Site-directed mutagenesis and laser-flash kinetic spectroscopy revealed that His-166 in SRI (i) plays a role in the proton transfers both to and from the Schiffbase, either as a structurally critical residue or possibly as a direct participant, (ii) is involved in the modulation of SIU photoreaction kinetics by HtrI, and (iii) modulates the pKa of Asp-76, an important residue in the photoactive site, through a long-distance electrostatic interaction. Computerized cell tracking and motion analysis demonstrated that (iv) His-166 is crucial in phototaxis signaling: a spectrum of substitutions either eliminate signaling or greatly perturb the activation process that produces attractant and repellent signaling states of the receptor.^ The signaling states of SRI are communicated to HtrI, whose oligomeric structure and conformational changes were investigated by engineered sulfhydryl probes. It was found that signaling by the SRI-HtrI complex involves reversible conformational changes within a preexisting HtrI dimer, which is likely accomplished through a slight winding or unwinding of the two HtrT monomers via their loose coiled coil association. To elucidate which domains of the Htr dimers confer specificity for interaction with SRI or SRII, chimeras of HtrI and HtrII were constructed. The only determinant needed for functional and specific interaction with SRI or SRII was found to be the four transmembrane segments of the HtrI or HtrII dimers, respectively. The entire cytoplasmic parts of HtrI and HtrII, which include the functionally important signaling and adaptation domains, were interchangeable.^ These observations support a model in which SRI and SRII undergo conformational changes coupled to light-induced proton transfers in their photoactive sites, and that lateral helix-helix interactions with their cognate transducers' 4-helix bundle in the membrane relay these conformational changes into different states of the Htr proteins which regulate the down-stream phosphotransfer pathway. ^

Conditional expression of Sry, and Wnt4 by gene-targeting: Studies in sexual and skeletal development

Sry and Wnt4 cDNAs were individually introduced into the ubiquitously-expressed Rosa26 ( R26) locus by gene targeting in embryonic stem (ES) cells to create a conditional gene expression system in mice. In the targeted alleles, expression of these cDNAs should be blocked by a neomycin resistance selection cassette that is flanked by loxP sites. Transgene expression should be activated after the blocking cassette is deleted by Cre recombinase. ^ To test this conditional expression system, I have bred R26-stop- Sry and R26-stop-Wnt4 heterozygotes with a MisRII-Cre mouse line that expresses Cre in the gonads of both sexes. Analysis of these two types of bigenic heterozygotes indicated that their gonads developed normally like those of wild types. However, one XX R26-Sry/R26-Sry; MisR2-Cre/+ showed epididymis-like structures resembling those of males. In contrast, only normal phenotypes were observed in XY R26-Wnt4/R26-Wnt4; MisR2-Cre /+ mice. To interpret these results, I have tested for Cre recombinase activity by Southern blot and transcription of the Sry and Wnt4 transgenes by RT-PCR. Results showed that bigenic mutants had insufficient activation of the transgenes in their gonads at E12.5 and E13.5. Therefore, the failure to observe mutant phenotypes may have resulted from low activity of MisR2-Cre recombination at the appropriate time. ^ Col2a1-Cre transgenic mice express Cre in differentiating chondrocytes. R26-Wnt4; Col2a1-Cre bigenic heterozygous mice were found to exhibit a dramatic alteration in growth presumably caused by Wnt4 overexpression during chondrogenesis. R26-Wnt4; Col2a1-Cre mice exhibited dwarfism beginning approximately 10 days after birth. In addition, they also had craniofacial abnormalities, and had delayed ossification of the lumbar vertebrate and pelvic bones. Histological analysis of the growth plates of R26-Wnt4; Col2a1-Cre mice revealed less structural organization and a delay in onset of the primary and secondary ossification centers. Molecular studies confirmed that overexpression of Wnt4 causes decreased proliferation and early maturation of chondrocytes. In addition, R26-Wnt4; Col2a1-Cre mice had decreased expression of vascular endothelial growth factor (VEGF), suggesting that defects in vascularization may contribute to the dwarf phenotype. Finally, 9-month-old R26-Wnt4; Col2a1-Cre mice had significantly more fat cells in the marrow cavities of their metaphysis long bones, implying that long-term overexpression of Wnt4may cause bone marrow pathologies. In conclusion, Wnt4 was activated by Col2a1-Cre recombinase and was overexpressed in the growth plate, resulting in aberrant proliferation and differentiation of chondrocytes, and ultimately leads to dwarfism in mice. ^

Elucidation of the role of 53BP1 in DNA double strand break (DSB) repair and tumor suppression

The protein p53 binding protein one (53BP1) was discovered in a yeast two-hybrid screen that used the DNA binding domain of p53 as bait. Cloning of full-length 53BP1 showed that this protein contains several protein domains which help make up the protein, which include two tandem BRCT domains and a amino-terminal serine/glutamine cluster domain (SCD). These are two protein domains are often seen in factors that are involved in the cellular response to DNA damage and control of cell cycle checkpoints and we hypothesize that 53BP1 is involved in the cellular response to DNA damage. In support of this hypothesis we observe that 53BP1 is phosphorylated and undergoes a dramatic nuclear re-localization in response to DNA damaging agents. 53BP1 also interacts with several factors that are important in the cellular response to DNA damage, such as the BRCA1 tumor suppressor, ATM and Rad3 related (ATR), and the phosphorylated version of the histone variant H2AX. Mice deficient in 53BP1 display increased sensitivity ionizing radiation (IR), a DNA damaging agent that introduces DNA double strand breaks (DSBs). In addition, 53BP1-deficient mice do not properly undergo the process of class switch recombination (CSR). We also observe that when a defect in 53BP1 is combined with a defect in p53; the resulting mice have an increased rate of formation of spontaneous tumors, notably the formation of B and T lineage lymphomas. The T lineage tumors arise by two distinct mechanisms: one driven by defects in cell cycle regulation and a second driven by defects in the ability to repair DNA DSBs. The B lineage tumors arise by the inability to repair DNA damage and over-expression of the oncogene c-myc. ^ With these observations, we conclude that not only does 53BP1 function in the cellular response to DNA damage, but it also works in concert with p53 to suppress tumor formation. ^

An ecologic comparison study of the impact of economic disadvantage on Texas Assessment of Knowledge and Skills performance, graduation rates, and readiness for higher education for students attending public and charter schools in Texas from 2004 to 2006

Purpose. No Child Left Behind aimed to "improve the academic achievement of the disadvantaged." The primary research question considered how academic achievement of those from economic disadvantage compared to those not from disadvantage? ^ Economically disadvantaged students can potentially have added academic disadvantage. Research shows low academic achievement can potentially result in drug abuse, youth violence, and teen pregnancy. ^ Methods. To compare the student populations, measures included TAKS results and academic indicator data collected by the Texas Education Agency. ^ Results. T-test analyses showed a significant difference between the economically and non-economically disadvantaged student populations in meeting the TAKS passing standard, graduation, and preparation for higher education.^ Conclusions. The achievement gap between students remained as indicated by the Texas testing program. More research and time are needed to observe if the desired impact on those from economic disadvantage will be reflected by academic achievement data.^

The association between maternal use of spermicides, condoms, intra-uterine devices or progesterone and major structural birth defects

Birth defects occur in 1 of every 33 babies born in the United States, and are the leading cause of infant death. Mothers using contraceptives that become pregnant may continue to use their contraceptives after their first missed menstrual period, thus exposing their baby in utero to the contraceptive product. Progesterone is also sometimes prescribed during the first trimester of pregnancy to mothers with a history of miscarriages or infertility problems. To ensure the safety of these products, it is important to investigate whether there is an increased occurrence of babies born with birth defects to mothers using various contraceptive methods or progesterone in early pregnancy. Using data from the National Birth Defects Prevention Study (NBDPS), an ongoing multi-state, population based case-control study, this study assessed maternal exposures to IUDs, spermicides, condoms and progesterone in early pregnancy. ^ Progesterone used for threatened miscarriage during the first three months of pregnancy was associated with an increased occurrence of hypoplastic left heart (adjusted odds ratios (OR) 2.24, 95% CI 1.13-4.21), perimembranous ventricular septal defects (OR 1.64, 95% CI 1.10-2.41), septal associations (OR 2.52, 95% CI 1.45-4.24), esophageal atresia (OR 1.82, 95% CI 1.04-3.08), and hypospadias (OR 2.12, 95% CI 1.41-3.18). Mothers using progesterone for injectable contraception had increased (OR > 2.5), but insignificant odds ratios for anencephaly, septal associations, small intestinal atresias and omphalocel. Progesterone used for fertility was not associated with an increased occurrence of any birth defects examined. ^ Mothers using progesterone for fertility assistance and threatened miscarriage were very similar with respect to their demographics and pregnancy history. They also both reported similar types of progesterone. Thus, if progesterone was a causal risk factor for birth defects we would have expected to observe similar increases in risk among mothers using progesterone for both indications. Because we predominantly observed increased associations among mothers using progesterone for threatened miscarriage but not fertility assistance, it is possible the increased associations we observed were confounded by indication (i.e. progesterone was administered for vaginal bleeding which occurred as a sequelae to the formation of a congenital anomaly. ^ No significant increased associations were observed between maternal spermicide use during pregnancy and 26 of 27 types of structural malformations. While multiple statistical tests were performed we observed first trimester maternal spermicide use to be associated with a significant increased occurrence of perimembranous ventricular septal defects (OR 2.21, 95% CI 1.16-4.21). A decreased occurrence (OR < 1.0) was observed for several categories of birth defects among mothers who conceived in the first cycle after discontinuing the use of spermicides (22 of 28) or male condoms (23 of 33). ^ Overall the percent of IUD use was similar between mothers of controls and mothers of all cases in aggregate (crude OR 1.05, 95% CI 0.61-1.84). Power was limited to detect significant associations between IUD use and birth defects, however mothers using an IUD in the month immediately prior to conception or during pregnancy were not associated with an increase of birth defects. Limb defects and amniotic band sequence previously reported to be associated with IUD use during pregnancy were not found to occur among any mothers reporting the use of an IUD during pregnancy.^

Yoga off the mat: How far and to whom do the benefits extend?

At least 15 million American adults have participated in yoga at least once in their lifetime (Saper, Eisenberg, Davis, Culpepper, & Phillips, 2004). The field of yoga research is relatively new in the United States, and the majority of studies have concentrated on yoga's effect on measures of physiology (cardiovascular disease, diabetes, obesity) or psychological measures of stress and anxiety. This review attempted to identify studies that had been conducted measuring a different set of outcome measures, specifically violence, trauma, eating, and other behavioral disorders. In 9 of 10 studies reviewed, researchers observed statistically significant effects of yoga interventions. Effects were most evident within multifaceted studies that combined intensive yoga practice with group discussion and training. Only one group (Mitchell, Mazzeo, Rausch, & Cooke, 2007) failed to observe any significant differences between yoga practice groups and control groups. Effects were seen in both sexes, although a majority of the studies were aimed specifically at women. All studies were limited by small sample size and lack of follow-up data. Future research should seek to increase sample size, to diversify recruitment to allow for the randomization of treatment and control groups, and to allow for long-term follow-up.^

A descriptive study on Hurricanes Katrina and Rita evacuees participating in the Harris County Housing Authority Disaster Housing Assistance Program

Natural disasters occur in various forms such as hurricanes, tsunamis, earthquakes, outbreaks, etc. The most unsettling aspect of a natural disaster is that it can strike at any moment. Over the past decade, our society has experienced an alarming increase of natural disasters. How to expeditiously respond and recover from natural disasters has become a precedent question for public health officials. To date, the most recent natural disaster was the January 12, 2010 earthquake in Haiti; however the most memorable was that of Hurricane Katrina (“Haiti Earthquake”, 2010). ^ This study provides insight on the need to develop a National Disaster Response and Recovery Program which effectively responds to natural disasters. The specific aims of this paper were to (1) observe the government’s role on federal, state and local levels in assisting Hurricanes Katrina and Rita evacuees, (2) assess the prevalence of needs among Hurricanes Katrina and Rita families participating in the Disaster Housing Assistance Program (DHAP) and (3) describe the level of progress towards “self sufficiency” for the DHAP families receiving case management social services. ^ Secondary data from a cross-sectional “Needs Assessment” questionnaire were analyzed. The questionnaire was administered initially and again six months later (follow-up) by H.A.U.L. case managers. The “Needs Assessment” questionnaire collected data regarding participants’ education, employment, transportation, child care, health resources, income, permanent housing and disability needs. Case managers determined the appropriate level of social services required for each family based on the data collected from the “Needs Assessment” questionnaire. ^ Secondary data provided by the H.A.U.L. were analyzed to determine the prevalence of needs among the DHAP families. In addition, differences measured between the initial and follow-up (at six months) questionnaires were analyzed to determine statistical significance between case management services provided and prevalence of needs among the DHAP families from initial to 6 months later at follow-up. The data analyzed describe the level of progress made by these families to achieve program “self sufficiency” (see Appendix A). Disaster assistance programs which first address basic human needs; then socioeconomic needs may offer an essential tool in aiding disaster affected communities quickly recover from natural disasters. ^

Volunteering: Demographic characteristics, motives and potential benefits from service

Volunteering is intricately woven into the fabric of our society. In 2009 alone, approximately 63.4 million Americans participated in volunteer activities, collectively donating over 8.1 billion service-hours (Corporation for National and Community Service [CNCS], 2010). Each service-hour is determined by the U.S. Bureau of Labor Statistics (2010) to be valued at $20.85/hr which translates to a national savings of $169 billion. Thus, we can clearly observe the significance of volunteer contribution to the overall benefit of society. In addition, there is now evidence that voluntary service may also benefit the actual volunteer, especially individuals who are 65+ years. As we reach 2020 this elderly class, composed of nearly 13 million (CNCS, 2010) Americans, will be of much consequence. Their potential to contribute in community-related efforts may save the U.S. billions in labor costs, and may also help reduce healthcare-related expenditures if volunteering proves to be a protective factor. In this literature review, we set out to explore the potential relationship between volunteer participation and increased mental and physical wellness. We also examined volunteer demographic characteristics and common motives for engaging in service-related activities. Analysis showed that volunteer work often combined low-impact physical activity and mental satisfaction from serving others, resulting in overall health benefit. Demographic characteristics displayed were consistent with previous studies and found that a majority of volunteers were female, White, married status, having received college degree or higher, employed, middle-high SES. In addition, age was seen to be a key characteristic in forecasting volunteer motivation and self-reported perceived health benefits.^

FZD6, MATN2 AND SLC25A32, POSSIBLE CANDIDATE GENES IN NONSYNDROMIC CLEFT LIP AND PALATE

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect with a multifactorial etiology. Despite decades of research, the genetic underpinnings of NSCLP still remain largely unexplained. A genome wide association study (GWAS) of a large NSCLP African American family with seven affected individuals across three generations found evidence for linkage at 8q21.3-24.12 (LOD = 2.98). This region contained three biologically relevant candidate genes: Frizzled-6 (FZD6) (LOD = 2.8), Matrilin-2 (MATN2) (LOD = 2.3), and Solute Carrier Family 25, Member 32 (SLC26A32) (LOD = 1.6). Sequencing of the coding regions and the 5’ and 3’ UTRs of these genes in two affected family members identified a rare intronic variant, rs138557689 (c.-153+432A>C), in FZD6. The rs138557689/C allele segregated with the NSCLP phenotype; in silico analysis predicted and EMSA analysis showed that the 138557689/C allele creates new DNA binding sites. FZD6 is part of the WNT pathway, which is involved in craniofacial development, including midface development and upper lip fusion. Our novel findings suggest that an alteration in FZD6 gene regulation may perturb this tightly controlled biological pathway and in turn contribute to the development of NSCLP in this family. Studies are underway to further define how the rs138557689/C variant affects expression of FZD6.

CONFORMATIONAL DYNAMICS OF K-RAS AND H-RAS PROTEINS: IS THERE FUNCTIONAL SPECIFICITY AT THE CATALYTIC DOMAIN?

Ras proteins serve as crucial signaling modulators in cell proliferation through their ability to hydrolyze GTP and exist in a GTP “on” state and GTP “off” state. There are three different human Ras isoforms: H-ras, N-ras and K-ras (4A and 4B). Although their sequence identity is very high at the catalytic domain, these isoforms differ in their ability to activate different effectors and hence different signaling pathways. Much of the previous work on this topic has attributed this difference to the hyper variable region of Ras proteins, which contains most of the sequence variance among the isoforms and encodes specificity for differential distribution in the membrane. However, we hypothesize that sequence variation on lobe II of Ras catalytic domain alters dynamics and leads to differential preference for different effectors or modulators. In this work, we used all atom molecular dynamics to analyze the dynamics in the catalytic domain of H-ras and K-ras. We have also analyzed the dynamics of a transforming mutant of H-ras and K-ras and further studied the dynamics of an effectorselective mutant of H-ras. Collectively we have determined that wild type K-ras is more dynamic than H-ras and that the structure of the effector binding loop more closely resembles that of the T35S Raf-selective mutant, possibly giving us a new view and insight into the v mode of effector specificity. Furthermore we have determined that specific mutations at the same location perturb the conformational equilibrium differently in H-ras and K-ras and that an enhanced oncogenic potential may arise from different structural perturbations for each point mutation of a specific isoform.

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.