A BAYESIAN APPROACH TO DOSE-RESPONSE ASSESSMENT AND DRUG-DRUG INTERACTION ANALYSIS: APPLICATION TO IN VITRO STUDIES

The considerable search for synergistic agents in cancer research is motivated by the therapeutic benefits achieved by combining anti-cancer agents. Synergistic agents make it possible to reduce dosage while maintaining or enhancing a desired effect. Other favorable outcomes of synergistic agents include reduction in toxicity and minimizing or delaying drug resistance. Dose-response assessment and drug-drug interaction analysis play an important part in the drug discovery process, however analysis are often poorly done. This dissertation is an effort to notably improve dose-response assessment and drug-drug interaction analysis. The most commonly used method in published analysis is the Median-Effect Principle/Combination Index method (Chou and Talalay, 1984). The Median-Effect Principle/Combination Index method leads to inefficiency by ignoring important sources of variation inherent in dose-response data and discarding data points that do not fit the Median-Effect Principle. Previous work has shown that the conventional method yields a high rate of false positives (Boik, Boik, Newman, 2008; Hennessey, Rosner, Bast, Chen, 2010) and, in some cases, low power to detect synergy. There is a great need for improving the current methodology. We developed a Bayesian framework for dose-response modeling and drug-drug interaction analysis. First, we developed a hierarchical meta-regression dose-response model that accounts for various sources of variation and uncertainty and allows one to incorporate knowledge from prior studies into the current analysis, thus offering a more efficient and reliable inference. Second, in the case that parametric dose-response models do not fit the data, we developed a practical and flexible nonparametric regression method for meta-analysis of independently repeated dose-response experiments. Third, and lastly, we developed a method, based on Loewe additivity that allows one to quantitatively assess interaction between two agents combined at a fixed dose ratio. The proposed method makes a comprehensive and honest account of uncertainty within drug interaction assessment. Extensive simulation studies show that the novel methodology improves the screening process of effective/synergistic agents and reduces the incidence of type I error. We consider an ovarian cancer cell line study that investigates the combined effect of DNA methylation inhibitors and histone deacetylation inhibitors in human ovarian cancer cell lines. The hypothesis is that the combination of DNA methylation inhibitors and histone deacetylation inhibitors will enhance antiproliferative activity in human ovarian cancer cell lines compared to treatment with each inhibitor alone. By applying the proposed Bayesian methodology, in vitro synergy was declared for DNA methylation inhibitor, 5-AZA-2'-deoxycytidine combined with one histone deacetylation inhibitor, suberoylanilide hydroxamic acid or trichostatin A in the cell lines HEY and SKOV3. This suggests potential new epigenetic therapies in cell growth inhibition of ovarian cancer cells.

Development of combined micro -preparation, separation, and mass spectrometry methods and application in the microdialysis study of neuropeptide metabolism

Two sets of mass spectrometry-based methods were developed specifically for the in vivo study of extracellular neuropeptide biochemistry. First, an integrated micro-concentration/desalting/matrix-addition device was constructed for matrix-assisted laser desorption ionization mass spectrometry (MALDI MS) to achieve attomole sensitivity for microdialysis samples. Second, capillary electrophoresis (CE) was incorporated into the above micro-liquid chromatography (LC) and MALDI MS system to provide two-dimensional separation and identification (i.e. electrophoretic mobility and molecular mass) for the analysis of complex mixtures. The latter technique includes two parts of instrumentation: (1) the coupling of a preconcentration LC column to the inlet of a CE capillary, and (2) the utilization of a matrix-precoated membrane target for continuous CE effluent deposition and for automatic MALDI MS analysis (imaging) of the CE track.^ Initial in vivo data reveals a carboxypeptidase A (CPA) activity in rat brain involved in extracellular neurotensin metabolism. Benzylsuccinic acid, a CPA inhibitor, inhibited neurotensin metabolite NT1-12 formation by 70%, while inhibitors of other major extracellular peptide metabolizing enzymes increased NT1-12 formation. CPA activity has not been observed in previous in vitro experiments. Next, the validity of the methodology was demonstrated in the detection and structural elucidation of an endogenous neuropeptide, (L)VV-hemorphin-7, in rat brain upon ATP stimulation. Finally, the combined micro-LC/CE/MALDI MS was used in the in vivo metabolic study of peptide E, a mu-selective opioid peptide with 25 amino acid residues. Profiles of 88 metabolites were obtained, their identity being determined by their mass-to-charge ratio and electrophoretic mobility. The results indicate that there are several primary cleavage sites in vivo for peptide E in the release of its enkephalin-containing fragments. ^

Genome-wide gene-gene interaction analysis for cardiovascular disease

Numerous studies have been carried out to try to better understand the genetic predisposition for cardiovascular disease. Although it is widely believed that multifactorial diseases such as cardiovascular disease is the result from effects of many genes which working alone or interact with other genes, most genetic studies have been focused on identifying of cardiovascular disease susceptibility genes and usually ignore the effects of gene-gene interactions in the analysis. The current study applies a novel linkage disequilibrium based statistic for testing interactions between two linked loci using data from a genome-wide study of cardiovascular disease. A total of 53,394 single nucleotide polymorphisms (SNPs) are tested for pair-wise interactions, and 8,644 interactions are found to be significant with p-values less than 3.5×10-11. Results indicate that known cardiovascular disease susceptibility genes tend not to have many significantly interactions. One SNP in the CACNG1 (calcium channel, voltage-dependent, gamma subunit 1) gene and one SNP in the IL3RA (interleukin 3 receptor, alpha) gene are found to have the most significant pair-wise interactions. Findings from the current study should be replicated in other independent cohort to eliminate potential false positive results.^

Assessing the Influence of Brain-based Instructional Methods and Its Impact on Systemic Change for Learning Communities

Background: As scholars who prepare future school leaders to be innovative instructional leaders for their learning communities, we are on the verge of a curriculum design revolution. The application of brain research findings promotes educational reform efforts to systemically change the way in which children experience school. However, most educators, school leaders, board members, and policy makers are ill prepared to reconsider the implications for assessment, pedagogy, school climate, daily schedules, and use of technology. This qualitative study asked future school leaders to reconsider how school leadership preparedness programs prepared them to become instructional leaders for the 21st century. The findings from this study will enhance the field of school leadership, challenging the current emphasis placed on standardized testing, traditional school calendars, assessments, monocultural instructional methods, and meeting the needs of diverse learning communities. [See PDF for complete abstract]

THE INTERACTION BETWEEN INTERNAL AND EXTERNAL INFORMATION ON RELATIONAL DATA SEARCH

People often use tools to search for information. In order to improve the quality of an information search, it is important to understand how internal information, which is stored in user’s mind, and external information, represented by the interface of tools interact with each other. How information is distributed between internal and external representations significantly affects information search performance. However, few studies have examined the relationship between types of interface and types of search task in the context of information search. For a distributed information search task, how data are distributed, represented, and formatted significantly affects the user search performance in terms of response time and accuracy. Guided by UFuRT (User, Function, Representation, Task), a human-centered process, I propose a search model, task taxonomy. The model defines its relationship with other existing information models. The taxonomy clarifies the legitimate operations for each type of search task of relation data. Based on the model and taxonomy, I have also developed prototypes of interface for the search tasks of relational data. These prototypes were used for experiments. The experiments described in this study are of a within-subject design with a sample of 24 participants recruited from the graduate schools located in the Texas Medical Center. Participants performed one-dimensional nominal search tasks over nominal, ordinal, and ratio displays, and searched one-dimensional nominal, ordinal, interval, and ratio tasks over table and graph displays. Participants also performed the same task and display combination for twodimensional searches. Distributed cognition theory has been adopted as a theoretical framework for analyzing and predicting the search performance of relational data. It has been shown that the representation dimensions and data scales, as well as the search task types, are main factors in determining search efficiency and effectiveness. In particular, the more external representations used, the better search task performance, and the results suggest the ideal search performance occurs when the question type and corresponding data scale representation match. The implications of the study lie in contributing to the effective design of search interface for relational data, especially laboratory results, which are often used in healthcare activities.

Dissecting the Interaction between p53 and TRIM24

Dissecting the Interaction of p53 and TRIM24 Aundrietta DeVan Duncan Supervisory Professor, Michelle Barton, Ph.D. p53, the “guardian of the genome”, plays an important role in multiple biological processes including cell cycle, angiogenesis, DNA repair and apoptosis. Because it is mutated in over 50% of cancers, p53 has been widely studied in established cancer cell lines. However, little is known about the function of p53 in a normal cell. We focused on characterizing p53 in normal cells and during differentiation. Our lab recently identified a novel binding partner of p53, Tripartite Motif 24 protein (TRIM24). TRIM24 is a member of the TRIM family of proteins, defined by their conserved RING, B-box, and coiled coil domains. Specifically, TRIM24 is a member of the TIF1 subfamily, which is characterized by PHD and Bromo domains in the C-terminus. Between the Coiled-coil and PHD domain is a linker region, 437 amino acids in length. This linker region houses important functions of TRIM24 including it’s site of interaction with nuclear receptors. TRIM24 is an E3-ubiquitin ligase, recently discovered to negatively regulate p53 by targeting it for degradation. Though it is known that Trim24 and p53 interact, it is not known if the interaction is direct and what effect this interaction has on the function of TRIM24 and p53. My study aims to elucidate the specific interaction domains of p53 and TRIM24. To determine the specific domains of p53 required for interaction with TRIM24, we performed co-immuoprecipitation (Co-IP) with recombinant full-length Flag-tagged TRIM24 protein and various deletion constructs of in vitro translated GST-p53, as well as the reverse. I found that TRIM24 binds both the carboxy terminus and DNA binding domain of p53. Furthermore, my results show that binding is altered when post-translational modifications of p53 are present, suggesting that the interaction between p53 and TRIM24 may be affected by these post-translational modifications. To determine the specific domains of TRIM24 required for p53 interaction, we performed GST pull-downs with in vitro translated, Flag-TRIM24 protein constructs and recombinant GST-p53 protein purified from E. coli. We found that the Linker region is sufficient for interaction of p53 and TRIM24. Taken together, these data indicate that the interaction between p53 and TRIM24 does occur in vitro and that interaction may be influenced by post-translational modifications of the proteins.

REGULATION OF CYTOSKELETAL ASSEMBLY: A STUDY OF THE INTERACTION OF FILAMIN AND F-ACTIN

Filamin is a high molecular weight (2 x 250,000) actin crosslinking protein found in a wide variety of cells and tissues. The most striking feature of filamin is its ability to crosslink F-actin filaments and cause ATP-independent gelation and contraction of F-actin solutions. The gelation of actin filaments by filamin involves binding to actin and crosslinking of the filaments by filamin self-association. In order to understand the role of filamin-actin interactions in the regulation of cytoskeletal assembly, two approaches were used. First, the structural relationship between self-association and actin-binding was examined using proteolytic fragments of filamin. Treatment of filamin with papain generated two major fragments, 90Kd and 180Kd. Upon incubation of the papain digest with F-actin and centrifugation at 100,000 x g, only the 180Kd fragment co-sedimented with F-actin. The binding of the 180Kd fragment, P180, was similar to native filamin in its sensitivity to ionic strength. Analytical gel filtration studies indicated that, unlike native filamin, P180 was monomeric and did not self-associate. Thermolysin treatment of P180 produced a 170Kd fragment, PT170, which no longer bound and co-sedimented with F-actin. These results suggested that filamin contained a discrete actin-binding domain. In order to locate the actin-binding domain, affinity purified antibodies to the papain and thermolysin sensitive regions of filamin were used in conjunction with filamin fragments generated by digestion with S. aureus V8 protease and elastase. The results indicated that the papain and thermolysin cleavage sites were close together, and, most likely, within 10Kd of one another. Taken together, these data suggest that filamin contains a discrete, internal actin-binding domain. The second approach was to use the non-crosslinking fragment P180 to develop a quantitative assay of filamin-actin binding. The binding of ('14)C-carboxyalkylated P180 was examined using the co-sedimentation assay. ('14)C-P180 binding to actin was equivalent to that of unlabelled P180 and exhibited comparable sensitivity of binding to changes in ionic strength. Within 5 min. of incubation the process had reached equilibrium. The specificity of binding was shown by the lack of binding of ('14)C-PT170. The binding of ('14)C-P180 was found to be a reversible and saturable process, with a K(,d) of 2 x 10('-7) M. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

A meta-analysis: Obesity and colorectal cancer screening

Of cancer death, colorectal cancer death ranks second in the United States. Obesity is an important risk factor for colorectal cancer (1). Early detection of colorectal cancer when it is localized can effectively reduce mortality of colorectal cancer and increase survival time of patients if they are treated. Also, previous studies showed that obese women were more likely to delay breast cancer screening and cervical cancer screening than normal weight women (2-5). However, results from prior studies demonstrating the relationship between obesity and colorectal cancer screening are not consistent. This research was done to conduct a meta-analysis of previous cross-sectional studies selected from the Medline database and to evaluate the association between obesity and colorectal cancer screening. While the odds ratio was not statistically different from one, the results from this meta-analysis under the random effects model showed that obese people are slightly less likely to have colorectal cancer screening compared to normal weight individuals (OR,0.93;95% CI 0.75-1.15). This meta-analysis was particularly sensitive to one individual study (6) and the effect of obesity on colorectal cancer screening was statistically significant (OR, 0.87; 95% CI, 0.81-0.92) after removing Heo's study. Further systematic studies focused on whether the effect of obesity on colorectal cancer screening is limited to women only are suggested. ^

Doctor-patient communication in the medical interaction: Context, implications, and practice

Effective communication; whether from an interpersonal, mass media, or global perspective, is a critical component in public health. It is an essential conduit in increasing public awareness of available health resources, potential health hazards and related disease prevention strategies, and in delivering better health care. Within this context, available literature asserts doctor-patient communication as central to healthcare delivery. It has been shown to affect patient health outcomes, satisfaction with care, adherence to treatment recommendations, and even understanding of medical information. While research supports the essential imperative of interventions aimed at teaching doctors and patients the communication skills necessary for a successful and meaningful medical interaction, most interventions to date, focus on teaching these communication skills to doctors and seem to rely, largely, on mass media for providing patients with the information needed to increase communication efficacy. This study sought to fill a significant gap in the doctor-patient communication literature by reviewing the context of the doctor-patient exchange in the medical interaction, the implications of this exchange in resulting care of the patient, and the potential improvements to practice through interventions aimed at improving the communication exchange. Closing with an evaluation of a patient-centered communication intervention, the “How to Talk to Your Doctor” (HTTTYD) program that combines previously identified optimal strategies for improving communication between doctors and patients, this study examined the patients’ perspective of their potential as better communicators in the medical interaction. ^ Specific Aims, Hypotheses or Questions (Aim I) To examine the context of health communication within a public health framework and its relation to health care delivery. (Aim II) To review doctor-patient communication as a central focus within health care delivery and the resulting implications to patient care. (Aim III) To assess the utility of interventions to improve doctor-patient communication. Specifically, to evaluate the effectiveness of a patient-centered community education intervention, the “How to Talk to Your Doctor” (HTTTYD) program, aimed at improving patient communication efficacy.^

INTERACTION BETWEEN BRK AND HER2 IN BREAST CANCER

INTERACTION BETWEEN BRK AND HER2 IN BREAST CANCER Midan Ai, Ph.D. Supervisory Professor: Zhen Fan, M.D. Breast tumor kinase (Brk) is a nonreceptor protein-tyrosine kinase that is highly expressed in approximately two thirds of breast cancers but is not detectable or is expressed at very low levels in normal mammary epithelium. Brk plays important roles in promoting proliferation, survival, invasion, and metastasis of breast cancer cells, but the mechanism(s) of which remain largely unknown. Recent studies showed that Brk is frequently co-overexpressed with human epidermal growth factor receptor-2 (HER2) and is physically associated with HER2 in breast cancer. The mechanism needs to be determined. In my studies, I found that high expression of HER2 is correlated with high expression of Brk in breast cancer cell lines. Silencing HER2 expression via RNA interference in HER2 over-expressed breast cancer cells resulted in Brk protein decrease and overexpression of HER2 in HER2 low-expressed breast cancer cells up-regulated Brk expression. The mechanism study indicated that overexpression of HER2 increased Brk protein stability. Brk was degraded through a Ca2+-dependent protease pathway involving calpain and HER2 stimulated Brk expression via inhibiting calpain activity. Calpastatin is a calpain endogenous inhibitor and the calpain-calpastatin system has been implicated in a number of cell physiological functions. HER2 restrained calpain activation via up-regulating calpastatin expression and HER2 downstream signaling, MAPK pathway, was involved in the regulation. Furthermore, silencing of Brk expression by RNA interference in HER2-overexpressing breast cancer cells decreased HER2-mediated cell proliferation, survival, invasion/metastasis potential and increased cell sensitivity to HER2 kinase inhibitor, lapatinib, treatment, indicating that Brk plays important roles in regulating and mediating the oncogenic functions of HER2. The Stat3 pathway played important roles in Brk mediated cell survival and invasion/metastasis in the context of HER2-overexpressing breast cancer cells. However, transgenic mice with inducible expression of constitutively active Brk (CA) in the mammary epithelium failed to develop malignant change in the mammary glands after Brk induction for 15 months which indicated that expression of Brk protein alone was not sufficiently to induce spontaneous breast tumor. Bitransgenic mice with co-expression of HER2/neu and inducible expression of Brk in the mammary epithelium developed multifocal mammary tumors, but there were no significant difference in the tumor occurring time, tumor size, tumor weight and tumor multiplicity between the mouse group with co-expression of Brk and HER2/neu and the mouse group with HER2/neu expression only.

Gene-gene interaction and gene-pathway analysis of genome-wide association for schizophrenia

Schizophrenia (SZ) is a complex disorder with high heritability and variable phenotypes that has limited success in finding causal genes associated with the disease development. Pathway-based analysis is an effective approach in investigating the molecular mechanism of susceptible genes associated with complex diseases. The etiology of complex diseases could be a network of genetic factors and within the genes, interaction may occur. In this work we argue that some genes might be of small effect that by itself are neither sufficient nor necessary to cause the disease however, their effect may induce slight changes to the gene expression or affect the protein function, therefore, analyzing the gene-gene interaction mechanism within the disease pathway would play crucial role in dissecting the genetic architecture of complex diseases, making the pathway-based analysis a complementary approach to GWAS technique. ^ In this study, we implemented three novel linkage disequilibrium based statistics, the linear combination, the quadratic, and the decorrelation test statistics, to investigate the interaction between linked and unlinked genes in two independent case-control GWAS datasets for SZ including participants of European (EA) and African (AA) ancestries. The EA population included 1,173 cases and 1,378 controls with 729,454 genotyped SNPs, while the AA population included 219 cases and 288 controls with 845,814 genotyped SNPs. We identified 17,186 interacting gene-sets at significant level in EA dataset, and 12,691 gene-sets in AA dataset using the gene-gene interaction method. We also identified 18,846 genes in EA dataset and 19,431 genes in AA dataset that were in the disease pathways. However, few genes were reported of significant association to SZ. ^ Our research determined the pathways characteristics for schizophrenia through the gene-gene interaction and gene-pathway based approaches. Our findings suggest insightful inferences of our methods in studying the molecular mechanisms of common complex diseases.^

The roles of electrostatic charge pairing in the interaction between cytochrome P-450 and NADPH-cytochrome P-450 reductase

The cytochrome P450 enzyme catalysis requires two electrons transferred from NADPH-cytochrome P450 reductase (reductase) to P450. Electrostatic charge-pairing has been proposed to be one of the major forces in the interaction between P450 and reductase. In order to obtain further insight into the molecular basis for the protein interaction, I used two methods, chemical modification and specific anti-peptide antibodies, to study the involvement and importance of charged amino acid residues. Acetylation of lysine residues of P450c and P450b by acetic anhydride dramatically inhibited the reductase-supported P450c-dependent ethoxycoumarin hydroxylation activity, but P450 activity supported by cumene hydroperoxide is relatively unchanged. The modification of lysine residues of P450c and P450b did not grossly disturb the protein conformation as revealed by several spectral studies. This differential effect of lysine modification on the P450 activity in the system reconstituted with reductase versus the system supported by cumene hydroperoxide suggested an important role for P450 lysine residues in the interaction with reductase. Using $\rm\sp{14}C$-acetic anhydride, P450 lysine residues were labelled and further identified on P450c and P450b. Those lysine residues are at position 97, 271, 279, and 407 for P450c, and 251, 384, 422, 433, and 473 for P450b. Alignment of those identified lysine residues on P450c and P450b with amino acid residues identified in other studies indicated those residues reside in three major sequence areas. Modification of arginine residues of P450b by phenylglyoxal and 2, 3-butanedione have no significant effect on P450 activity either supported by NADPH and reductase or supported by cumene hydroperoxide. Further studies using $\rm\sp{14}C$-phenylglyoxal reveals that no incorporation of phenylglyoxal into P450b was found. These results demonstrated a predominant role of lysine residues of P450 in the electrostatic interaction with reductase. To understand the protein binding sites on each of P450 and reductase, I generated three anti-peptide antibodies against regions on reductase and five anti-peptide antibodies against five putative reductase binding sites on P450c. These anti-peptide antibodies were affinity purified and characterized on ELISA and by Western blot analysis. Inhibition experiments using these antibodies demonstrated that regions 109-120 and 204-220 of reductase are probably the two major binding sites for P450. The association of reductase with cytochromes P450 and cytochrome c may rely on different mechanisms. The data from experiments using anti-peptide (P450c) antibodies supports the important role of P450c lysine residues 271/279 and 458/460 in the interaction with reductase. ^

Secondary data analysis of a pilot study of use of an instructional video for dietary assessment

Introduction. Food frequency questionnaires (FFQ) are used study the association between dietary intake and disease. An instructional video may potentially offer a low cost, practical method of dietary assessment training for participants thereby reducing recall bias in FFQs. There is little evidence in the literature of the effect of using instructional videos on FFQ-based intake. Objective. This analysis compared the reported energy and macronutrient intake of two groups that were randomized either to watch an instructional video before completing an FFQ or to view the same instructional video after completing the same FFQ. Methods. In the parent study, a diverse group of students, faculty and staff from Houston Community College were randomized to two groups, stratified by ethnicity, and completed an FFQ. The "video before" group watched an instructional video about completing the FFQ prior to answering the FFQ. The "video after" group watched the instructional video after completing the FFQ. The two groups were compared on mean daily energy (Kcal/day), fat (g/day), protein (g/day), carbohydrate (g/day) and fiber (g/day) intakes using descriptive statistics and one-way ANOVA. Demographic, height, and weight information was collected. Dietary intakes were adjusted for total energy intake before the comparative analysis. BMI and age were ruled out as potential confounders. Results. There were no significant differences between the two groups in mean daily dietary intakes of energy, total fat, protein, carbohydrates and fiber. However, a pattern of higher energy intake and lower fiber intake was reported in the group that viewed the instructional video before completing the FFQ compared to those who viewed the video after. Discussion. Analysis of the difference between reported intake of energy and macronutrients showed an overall pattern, albeit not statistically significant, of higher intake in the video before versus the video after group. Application of instructional videos for dietary assessment may require further research to address the validity of reported dietary intakes in those who are randomized to watch an instructional video before reporting diet compared to a control groups that does not view a video.^

Motivational interviewing and feedback for college drinkers: Handling missing values with complete case analysis and multiple imputation

Objective: In this secondary data analysis, three statistical methodologies were implemented to handle cases with missing data in a motivational interviewing and feedback study. The aim was to evaluate the impact that these methodologies have on the data analysis. ^ Methods: We first evaluated whether the assumption of missing completely at random held for this study. We then proceeded to conduct a secondary data analysis using a mixed linear model to handle missing data with three methodologies (a) complete case analysis, (b) multiple imputation with explicit model containing outcome variables, time, and the interaction of time and treatment, and (c) multiple imputation with explicit model containing outcome variables, time, the interaction of time and treatment, and additional covariates (e.g., age, gender, smoke, years in school, marital status, housing, race/ethnicity, and if participants play on athletic team). Several comparisons were conducted including the following ones: 1) the motivation interviewing with feedback group (MIF) vs. the assessment only group (AO), the motivation interviewing group (MIO) vs. AO, and the intervention of the feedback only group (FBO) vs. AO, 2) MIF vs. FBO, and 3) MIF vs. MIO.^ Results: We first evaluated the patterns of missingness in this study, which indicated that about 13% of participants showed monotone missing patterns, and about 3.5% showed non-monotone missing patterns. Then we evaluated the assumption of missing completely at random by Little's missing completely at random (MCAR) test, in which the Chi-Square test statistic was 167.8 with 125 degrees of freedom, and its associated p-value was p=0.006, which indicated that the data could not be assumed to be missing completely at random. After that, we compared if the three different strategies reached the same results. For the comparison between MIF and AO as well as the comparison between MIF and FBO, only the multiple imputation with additional covariates by uncongenial and congenial models reached different results. For the comparison between MIF and MIO, all the methodologies for handling missing values obtained different results. ^ Discussions: The study indicated that, first, missingness was crucial in this study. Second, to understand the assumptions of the model was important since we could not identify if the data were missing at random or missing not at random. Therefore, future researches should focus on exploring more sensitivity analyses under missing not at random assumption.^