Automatic quantitation of multiple sclerosis lesions on MR images

A two-pronged approach for the automatic quantitation of multiple sclerosis (MS) lesions on magnetic resonance (MR) images has been developed. This method includes the design and use of a pulse sequence for improved lesion-to-tissue contrast (LTC) and seeks to identify and minimize the sources of false lesion classifications in segmented images. The new pulse sequence, referred to as AFFIRMATIVE (Attenuation of Fluid by Fast Inversion Recovery with MAgnetization Transfer Imaging with Variable Echoes), improves the LTC, relative to spin-echo images, by combining Fluid-Attenuated Inversion Recovery (FLAIR) and Magnetization Transfer Contrast (MTC). In addition to acquiring fast FLAIR/MTC images, the AFFIRMATIVE sequence simultaneously acquires fast spin-echo (FSE) images for spatial registration of images, which is necessary for accurate lesion quantitation. Flow has been found to be a primary source of false lesion classifications. Therefore, an imaging protocol and reconstruction methods are developed to generate "flow images" which depict both coherent (vascular) and incoherent (CSF) flow. An automatic technique is designed for the removal of extra-meningeal tissues, since these are known to be sources of false lesion classifications. A retrospective, three-dimensional (3D) registration algorithm is implemented to correct for patient movement which may have occurred between AFFIRMATIVE and flow imaging scans. Following application of these pre-processing steps, images are segmented into white matter, gray matter, cerebrospinal fluid, and MS lesions based on AFFIRMATIVE and flow images using an automatic algorithm. All algorithms are seamlessly integrated into a single MR image analysis software package. Lesion quantitation has been performed on images from 15 patient volunteers. The total processing time is less than two hours per patient on a SPARCstation 20. The automated nature of this approach should provide an objective means of monitoring the progression, stabilization, and/or regression of MS lesions in large-scale, multi-center clinical trials. ^

An Innovative Phase I Trial Design Allowing for the Identification of Multiple Potential Maximum Tolerated Doses with Combination Therapy of Targeted Agents

Treatment for cancer often involves combination therapies used both in medical practice and clinical trials. Korn and Simon listed three reasons for the utility of combinations: 1) biochemical synergism, 2) differential susceptibility of tumor cells to different agents, and 3) higher achievable dose intensity by exploiting non-overlapping toxicities to the host. Even if the toxicity profile of each agent of a given combination is known, the toxicity profile of the agents used in combination must be established. Thus, caution is required when designing and evaluating trials with combination therapies. Traditional clinical design is based on the consideration of a single drug. However, a trial of drugs in combination requires a dose-selection procedure that is vastly different than that needed for a single-drug trial. When two drugs are combined in a phase I trial, an important trial objective is to determine the maximum tolerated dose (MTD). The MTD is defined as the dose level below the dose at which two of six patients experience drug-related dose-limiting toxicity (DLT). In phase I trials that combine two agents, more than one MTD generally exists, although all are rarely determined. For example, there may be an MTD that includes high doses of drug A with lower doses of drug B, another one for high doses of drug B with lower doses of drug A, and yet another for intermediate doses of both drugs administered together. With classic phase I trial designs, only one MTD is identified. Our new trial design allows identification of more than one MTD efficiently, within the context of a single protocol. The two drugs combined in our phase I trial are temsirolimus and bevacizumab. Bevacizumab is a monoclonal antibody targeting the vascular endothelial growth factor (VEGF) pathway which is fundamental for tumor growth and metastasis. One mechanism of tumor resistance to antiangiogenic therapy is upregulation of hypoxia inducible factor 1α (HIF-1α) which mediates responses to hypoxic conditions. Temsirolimus has resulted in reduced levels of HIF-1α making this an ideal combination therapy. Dr. Donald Berry developed a trial design schema for evaluating low, intermediate and high dose levels of two drugs given in combination as illustrated in a recently published paper in Biometrics entitled “A Parallel Phase I/II Clinical Trial Design for Combination Therapies.” His trial design utilized cytotoxic chemotherapy. We adapted this design schema by incorporating greater numbers of dose levels for each drug. Additional dose levels are being examined because it has been the experience of phase I trials that targeted agents, when given in combination, are often effective at dosing levels lower than the FDA-approved dose of said drugs. A total of thirteen dose levels including representative high, intermediate and low dose levels of temsirolimus with representative high, intermediate, and low dose levels of bevacizumab will be evaluated. We hypothesize that our new trial design will facilitate identification of more than one MTD, if they exist, efficiently and within the context of a single protocol. Doses gleaned from this approach could potentially allow for a more personalized approach in dose selection from among the MTDs obtained that can be based upon a patient’s specific co-morbid conditions or anticipated toxicities.

DRUG METABOLISM IN LIVER TUMORS: PURIFICATION AND CHARACTERIZATION OF NADPH-CYTOCHROME-P450 REDUCTASE AND THE RESOLUTION OF MULTIPLE FORMS OF CYTOCHROME P-450

The hydroxylation of N- and O-methyl drugs and a polycyclic hydrocarbon has been demonstrated in microsomes prepared from two transplantable Morris hepatomas (i.e., 7288C. t.c. and 5123 t.c.(H). The hydroxylation rates of the drug benzphetamine and the polycyclic hydrocarbon benzo {(alpha)} pyrene by tumor microsomes were inducible 2 to 3-fold and 2-fold, respectively by pretreatment of rats with phenobarbital/hydrocortisone. Hepatoma 5123t.c.(h) microsomal hydroxylation activities were more inducible after these pretreatments than hepatoma 7288C.t.c. Two chemotherapeutic drugs (cyclophosphamide and isophosphamide) were shown to be mutagenic after activation by the tumor hemogenate with the TA100 strain of Salmonella typhimurium bacteria. NADPH-cytochrome P-450 was purified from phenobarbital/hydrocortisone treated rat hepatoma 5123t.c.(H) microsomes 353-fold with a specific activity 63.6 nmol of cytochrome c reduced per min per mg of protein. The purified enzyme, has an apparent molecular weight of 79,500 daltons, and contained an equal molar ratio of FMN and FAD, with a total flavin content of 16.4 nmol per mg of protein. The purified enzyme also catalyzed electron transfer to artificial electron acceptors with the K(,m) values of the hepatoma reductase similar to those of purified liver reductase. The K(,m) value of the hepatoma reductase (13 uM) for NADPH was similar to that of purified liver reductase (5.0 uM). In addition the purified hepatoma reductase was immunochemically similar to the liver reductase.^ Hepatoma cytochrome P-450, the hemeprotein component of the hepatoma microsomes of rats pretreated with phenobarbital/hydrocortisone. The resolution of the six forms was achieved by the DE-53 ion-exchange chromatography, and further purified by hydroxyapatite. The six different fractions that contained P-450 activity, had specific contents from 0.47 to 1.75 nmol of cytochrome P-450 per mg of protein, and indicated a 2 to 9-fold purification as compared to the original microsomes. In addition, difference spectra, molecular weights and immunological results suggest there are at least six different forms of cytochrome P-450 in hepatoma 5123 t.c.(H). ^

Variables affecting the stability of multiple cloning sites and hairpin structures in retroviral vectors

Retroviruses are RNA viruses that replicate through a double-stranded DNA intermediate. The viral enzyme reverse transcriptase copies the retroviral genomic RNA into this DNA intermediate through the process of reverse transcription. Many variables can affect the fidelity of reverse transcriptase during reverse transcription, including specific sequences within the retroviral genome. ^ Previous studies have observed that multiple cloning sites (MCS) and sequences predicted to form stable hairpin structures are hotspots for deletion during retroviral replication. The studies described in this dissertation were performed to elucidate the variables that affect the stability of MCS and hairpin structures in retroviral vectors. Two series of retroviral vectors were constructed and characterized in these studies. ^ Spleen necrosis virus-based vectors were constructed containing separate MCS insertions of varying length, orientation, and symmetry. The only MCS that was a hotspot for deletion formed a stable hairpin structure. Upon more detailed study, the MCS previously reported as a hotspot for deletion was found to contain a tandem linker insertion that formed a hairpin structure. Murine leukemia virus-based vectors were constructed containing separate sequence insertions of either inverted repeat symmetry (122IR) that could form a hairpin structure, or little symmetry (122c) that would form a less stable structure. These insertions were made into either the neomycin resistance marker ( neo) or the hygromycin resistance marker (hyg) of the vector. 122c was stable in both neo and hyg, while 122IR was preferentially deleted in neo and was remarkably unstable in hyg. ^ These results suggest that MCS are hotspots for deletion in retroviral vectors if they can form hairpin structures, and that hairpin structures can be highly unstable at certain locations in retroviral vectors. This information may contribute to improved design of retroviral vectors for such uses as human gene therapy, and will contribute to a greater understanding of the basic science of retroviral reverse transcription. ^

Application of stellar photometry to the analysis of microarray images

Improvements in the analysis of microarray images are critical for accurately quantifying gene expression levels. The acquisition of accurate spot intensities directly influences the results and interpretation of statistical analyses. This dissertation discusses the implementation of a novel approach to the analysis of cDNA microarray images. We use a stellar photometric model, the Moffat function, to quantify microarray spots from nylon microarray images. The inherent flexibility of the Moffat shape model makes it ideal for quantifying microarray spots. We apply our novel approach to a Wilms' tumor microarray study and compare our results with a fixed-circle segmentation approach for spot quantification. Our results suggest that different spot feature extraction methods can have an impact on the ability of statistical methods to identify differentially expressed genes. We also used the Moffat function to simulate a series of microarray images under various experimental conditions. These simulations were used to validate the performance of various statistical methods for identifying differentially expressed genes. Our simulation results indicate that tests taking into account the dependency between mean spot intensity and variance estimation, such as the smoothened t-test, can better identify differentially expressed genes, especially when the number of replicates and mean fold change are low. The analysis of the simulations also showed that overall, a rank sum test (Mann-Whitney) performed well at identifying differentially expressed genes. Previous work has suggested the strengths of nonparametric approaches for identifying differentially expressed genes. We also show that multivariate approaches, such as hierarchical and k-means cluster analysis along with principal components analysis, are only effective at classifying samples when replicate numbers and mean fold change are high. Finally, we show how our stellar shape model approach can be extended to the analysis of 2D-gel images by adapting the Moffat function to take into account the elliptical nature of spots in such images. Our results indicate that stellar shape models offer a previously unexplored approach for the quantification of 2D-gel spots. ^

Obesity and risk of multiple myeloma: A case-control study

Introduction. Several studies have reported a positive association of body mass index (BMI) with multiple myeloma; however, the period of adulthood where BMI is most important remains unclear. In addition, it is well known that body fat is associated with both sex-steroid hormone storage and with increasing insulin levels; therefore, it was hypothesized that the association between obesity and multiple myeloma may be attributed to increased aromatization of androgen in adipose tissue. Objective. The overall objective of this case-control study was to determine whether multiple myeloma cases had higher BMI and greater adult weight gain relative to healthy controls. In addition, we tested the hypothesis that hormone replacement therapy use among women will further increase the association between BMI and risk of multiple myeloma. This study used data from a pilot case-control study at M.D. Anderson Cancer Center (MDACC), entitled Etiology of multiple myeloma, directed by Dr. Sara Strom and Dr. Sergio Giralt. Methods. The pilot study recruited a total of 122 cases of histopathologically confirmed multiple myeloma from MDACC. Controls (n=183) were selected from a database of random digit dialing controls accrued in the Department of Epidemiology at MDACC and were frequency matched to the cases on age (±5 years), gender, and race/ethnicity. Demographic and risk factor information were obtained from all participants who completed a self-administered questionnaire. Items included in the questionnaire include demographic information, height and weight at age 25, 40 and current/diagnosis, medical history, family history of cancer, smoking and alcohol use. Statistical analysis. Initial descriptive analysis included Student's t-test and Pearson's chi-squared tests. Odds ratios and 95% confidence intervals were calculated to quantify the association between the variables of interest and multiple myeloma. A multivariable model will be developed using unconditional logistic regression. Results. MM cases were 1.79 times (95% CI=0.99-3.32) more likely to have been overweight or obese (BMI > 25 kg/m2) at age 25 relative to healthy controls after controlling for age, gender, race/ethnicty, education and family history of cancer. Being overweight or obese at age 40 was not significantly associated with mutliple myeloma risk (OR=1.42, 95% CI=0.86-2.34) nor was being overweight or obses at diagnosis (OR=1.43, 95% CI=0.78, 2.63). We observed a statistically significant 2-fold increased odds of multiple myeloma in individuals who gained more than 4.7 kg during between 25 and 40 years (OR=1.97, 95% CI=1.15-3.39). When assessing HRT as a modifier of the BMI and multiple myeloma association among women (N=123), no association between obesity and MM status was observed among women who have never used HRT (OR=0.60, 95% CI=0.23-1.61; n=73). Yet among women who have ever used HRT (n=50), being overweight or obese was associated with an increase in MM risk (OR=2. 93, 95% CI=0.81-10.6) after adjusting for age; however, the association was not statistically significant. Significance. This study provides further evidence that increased BMI increases the risk of multiple myeloma. Furthermore, among women, HRT use may modify risk of disease. ^

Safety, pharmacokinetics, and antiretroviral activity of multiple doses of ibalizumab (formerly TNX-355), an anti-CD4 monoclonal antibody in HIV-1 infected adults

Context: Despite tremendous strides in HIV treatment over the past decade, resistance remains a major problem. A growing number of patients develop resistance and require new therapies to suppress viral replication. ^ Objective: To assess the safety of multiple administrations of the anti-CD4 receptor (anti-CD4) monoclonal antibody ibalizumab given as intravenous (IV) infusions, in three dosage regimens, in subjects infected with human immunodeficiency virus (HIV-1). ^ Design: Phase 1, multi-center, open-label, randomized clinical trial comparing the safety, pharmacokinetics and antiviral activity of three dosages of ibalizumab. ^ Setting: Six clinical trial sites in the United States. ^ Participants: A total of twenty-two HIV-positive patients on no anti-retroviral therapy or a stable failing regimen. ^ Intervention: Randomized to one of two treatment groups in Arms A and B followed by non-randomized enrollment in Arm C. Patients randomized to Arm A received 10 mg/kg of ibalizumab every 7 days, for a total of 10 doses; patients randomized to Arm B received a total of six doses of ibalizumab; a single loading dose of 10 mg/kg on Day 1 followed by five maintenance doses of 6 mg/kg every 14 days, starting at Week 1. Patients assigned to Arm C received 25 mg/kg of ibalizumab every 14 days for a total of 5 doses. All patients were followed for safety for an additional 7 to 8 weeks. ^ Main Outcome Measures: Clinical and laboratory assessments of safety and tolerability of multiple administrations of ibalizumab in HIV-infected patients. Secondary measures of efficacy include HIV-1 RNA (viral load) measurements. ^ Results: 21 patients were treatment-experienced and 1 was naïve to HIV therapy. Six patients were failing despite therapy and 15 were on no current HIV treatment. Mean baseline viral load (4.78 log 10; range 3.7-5.9) and CD4+ cell counts (332/μL; range 89-494) were similar across cohorts. Mean peak decreases in viral load from baseline of 0.99 log10(1.11 log10, and 0.96 log 10 occurred by Wk 2 in Cohorts A, B and C, respectively. Viral loads decreased by >1.0 log10 in 64%; 4 patients viral loads were suppressed to < 400 copies/mL. Viral loads returned towards baseline by Week 9 with reduced susceptibility to ibalizumab. CD4+ cell counts rose transiently and returned toward baseline. Maximum median elevations above BL in CD4+ cell counts for Cohorts A, B and C were +257, +198 and +103 cells/μL, respectively and occurred within 3 Wks in 16 of 22 subjects. The half-life of ibalizumab was 3-3.5 days and elimination was characteristic of capacity-limited kinetics. Administration of ibalizumab was well tolerated. Four serious adverse events were reported during the study. None of these events were related to study drug. Headache, nausea and cough were the most frequently reported treatment emergent adverse events and there were no laboratory abnormalities related to study drug. ^ Conclusions: Ibalizumab administered either weekly or bi-weekly was safe, well tolerated, and demonstrated antiviral activity. Further studies with ibalizumab in combination with standard antiretroviral treatments are warranted.^

Comparison of multiple comparison methods for identifying differential gene expression in simulated and real papillary thyroid cancer microarray data

The difficulty of detecting differential gene expression in microarray data has existed for many years. Several correction procedures try to avoid the family-wise error rate in multiple comparison process, including the Bonferroni and Sidak single-step p-value adjustments, Holm's step-down correction method, and Benjamini and Hochberg's false discovery rate (FDR) correction procedure. Each multiple comparison technique has its advantages and weaknesses. We studied each multiple comparison method through numerical studies (simulations) and applied the methods to the real exploratory DNA microarray data, which detect of molecular signatures in papillary thyroid cancer (PTC) patients. According to our results of simulation studies, Benjamini and Hochberg step-up FDR controlling procedure is the best process among these multiple comparison methods and we discovered 1277 potential biomarkers among 54675 probe sets after applying the Benjamini and Hochberg's method to PTC microarray data.^

Generalized fuzzy clustering for segmentation of multi-spectral magnetic resonance images.

An integrated approach for multi-spectral segmentation of MR images is presented. This method is based on the fuzzy c-means (FCM) and includes bias field correction and contextual constraints over spatial intensity distribution and accounts for the non-spherical cluster's shape in the feature space. The bias field is modeled as a linear combination of smooth polynomial basis functions for fast computation in the clustering iterations. Regularization terms for the neighborhood continuity of intensity are added into the FCM cost functions. To reduce the computational complexity, the contextual regularizations are separated from the clustering iterations. Since the feature space is not isotropic, distance measure adopted in Gustafson-Kessel (G-K) algorithm is used instead of the Euclidean distance, to account for the non-spherical shape of the clusters in the feature space. These algorithms are quantitatively evaluated on MR brain images using the similarity measures.

Multiple site phosphorylation of tyrosine hydroxylase: A potential role for protein kinase C in the regulation of catecholamine synthesis

Tyrosine hydroxylase (TH), the initial and rate limiting enzyme in the catecholaminergic biosynthetic pathway, is phosphorylated on multiple serine residues by multiple protein kinases. Although it has been demonstrated that many protein kinases are capable of phosphorylating and activating TH in vitro, it is less clear which protein kinases participate in the physiological regulation of catecholamine synthesis in situ. These studies were designed to determine if protein kinase C (PK-C) plays such a regulatory role.^ Stimulation of intact bovine adrenal chromaffin cells with phorbol esters results in stimulation of catecholamine synthesis, tyrosine hydroxylase phosphorylation and activation. These responses are both time and concentration dependent, and are specific for those phorbol ester analogues which activate PK-C. RP-HPLC analysis of TH tryptic phosphopeptides indicate that PK-C phosphorylates TH on three putative sites. One of these (pepetide 6) is the same as that phosphorylated by both cAMP-dependent protein kinase (PK-A) and calcium/calmodulin-dependent protein kinase (CaM-K). However, two of these sites (peptides 4 and 7) are unique, and, to date, have not been shown to be phosphorylated by any other protein kinase. These peptides correspond to those which are phosphorylated with a slow time course in response to stimulation of chromaffin cells with the natural agonist acetylcholine. The activation of TH produced by PK-C is most closely correlated with the phosphorylation of peptide 6. But, as evident from pH profiles of tyrosine hydroxylase activity, phosphorylation of peptides 4 and 7 affect the expression of the activation produced by phosphorylation of peptide 6.^ These data support a role for PK-C in the control of TH activity, and suggest a two stage model for the physiological regulation of catecholamine synthesis by phosphorylation in response to cholinergic stimulation. An initial fast response, which appears to be mediated by CaM-K, and a slower, sustained response which appears to be mediated by PK-C. In addition, the multiple site phosphorylation of TH provides a mechanism whereby the regulation of catecholamine synthesis appears to be under the control of multiple protein kinases, and allows for the convergence of multiple, diverse physiological and biochemical signals. ^

Involvement of human chromosome 17 in the control of the metastatic phenotype in melanoma

Fusion of nonmetastatic murine melanoma K1735 C19H cells with metastatic human melanoma A375 C15N cells resulted in a hybrid (termed H7) which was highly metastatic in a nude mouse model. The H7 hybrid retained chromosome 17 as the sole intact human chromosome in the cell. A lung bioassay showed that the K1735 C19H cells were present in the lungs of nude mice with s.c. tumors, yet at 6-weeks after tumor resection, no cells remained in the lung and therefore did not form lung metastases. Examination of various phenotypic properties such as in vivo and in vitro growth demonstrated that phenotypically the H7 hybrid was most like the K1735 C19H cell line except for its metastatic ability. In contrast the H7 hybrid cells containing single or multiple copies of human chromosome 17 with a point mutation at codon 249 (arg-gly) of the p53 gene, readily formed lung metastases. A plasmid containing the human p53 from the H7 hybrid and four other contructs with mutations at codon 143 (val-arg), 175 (arg-his), 249 (arg-ser) and 273 (arg-his) were transfected into K1735 C19H cells. K1735 C19H cells expressing human p53 genes with mutations at codons 249, both the arg-ser mutation and the mutation from the H7 hybrid and 273 produced significantly more lung metastases.^ In vitro assays demonstrated that responses to various cytotoxic and DNA damaging agents varied with the presence of mutant p53 and with the type of agent used. When cultured in mouse lung-conditioned medium, the K1735 C19H cell line was growth-inhibited, while cells containing a mutant human p53 (either on the whole chromosome 17, as in the H7 hybrid cells or from a stably transfected construct) were growth stimulated. Western blot analysis of lung-conditioned media derived from either 6-month or 15-month old mice has detected high levels of soluble Fas ligand in the medium from older animals. Comparison of the levels of Fas receptor on the K1735 C19H cell line and the H7 hybrid were almost identical, but 50% of the K1735 C19H cells were killed in the presence of anti-Fas antibody as opposed to 7% of the H7 hybrid cells. The growth-inhibitory effects of the lung-conditioned medium on the K1735 C19H cells were abrogated by coculture with Fas-Fc, which competes with the Fas ligand for receptor binding. Growth-inhibition of the K1735 C19H was 54% when cultured in 60 $\mu$g/0.2 ml lung-conditioned medium and a control Fc, with only 9% inhibition in 60 $\mu$g/0.2 ml lung-conditioned medium and Fas-Fc. Growth of the H7 cells and K1735 C19H cells transfected with various mutant human p53 genes were unchanged by the presence of either the control Fc or the Fas-Fc. This indicates that the presence of human chromosome 17, and mutant p53 in part protects the cells from Fas:Fas ligand induced apoptosis, and allows the growth of lung metastases. ^

CYTOLOGICAL, ULTRASTRUCTURAL, AND BIOCHEMICAL STUDIES OF THE STRUCTURE OF PREMATURELY CONDENSED CHROMOSOMES

The phenomenon of premature chromosome condensation, resulting from fusion between mitotic and interphase cells, includes dissolution of the interphase nuclear framework, thus allowing a direct visualization of interphase chromosomes. Light microscope morphology of prematurely condensed chromosomes (PCC) from synchronized HeLa cells supports the model of an interphase "chromosome condensation cycle". PCC are increasingly attenuated as cells progress through G(,1). A maximum degree of decondensation is observed at active sites of DNA replication during S phase, and a condensed morphology is rapidly resumed following completion of replication of a chromosome segment.^ To permit ultrastructural and biochemical studies of PCC, a procedure was developed to induce premature chromosome condensation at high frequency. This was achieved by polyethylene glycol (PEG)-mediated fusion of a dense monolayer of mitotic and interphase cells induced by centrifugation onto lectin-coated culture dishes. Using this method, PCC induction frequencies of 60-90% are routinely obtained.^ Scanning electron microscope analysis of PCC spreads revealed that the extension of PCC during progression through G(,1) is accompanied by a transition of the basic 30 nm chromatin fiber from tightly packed looping fibers to extended longitudinal fibers. Sites of active DNA replication is S-PCC were indicated to be organized a single longitudinal fibers. Following replication of a chromosome segment, a rapid reorganization from the extended longitudinal fiber to packed looping fibers occurs. The postreplication maturation process appears to include the assembly of a chromosome core consisting of multiple longitudinal fibers.^ The role of histone H1 phosphorylation in PCC formation was investigated by acidurea polyacrylamide gel electrophoresis of total histone extracted from metaphase chromosomes and PCC following high frequency fusion. This investigation failed to demonstrate an extensive phosphorylation of H1 associated with PCC formation. However, significant dephosphorylation of superphosphorylated metaphase chromosome H1 was observed, indicating that interphase H1-phosphatase activity is dominant over metaphase H1 kinase activity. These observations provide evidence against models suggesting a role for H1 superphosphorylation in triggering mitotic condensation of chromosomes. ^

Multiple imputation for missing continuous and dichotomous data: The effects of patterns of missingness and variable correlations on type I error with small samples

The purpose of this study is to investigate the effects of predictor variable correlations and patterns of missingness with dichotomous and/or continuous data in small samples when missing data is multiply imputed. Missing data of predictor variables is multiply imputed under three different multivariate models: the multivariate normal model for continuous data, the multinomial model for dichotomous data and the general location model for mixed dichotomous and continuous data. Subsequent to the multiple imputation process, Type I error rates of the regression coefficients obtained with logistic regression analysis are estimated under various conditions of correlation structure, sample size, type of data and patterns of missing data. The distributional properties of average mean, variance and correlations among the predictor variables are assessed after the multiple imputation process. ^ For continuous predictor data under the multivariate normal model, Type I error rates are generally within the nominal values with samples of size n = 100. Smaller samples of size n = 50 resulted in more conservative estimates (i.e., lower than the nominal value). Correlation and variance estimates of the original data are retained after multiple imputation with less than 50% missing continuous predictor data. For dichotomous predictor data under the multinomial model, Type I error rates are generally conservative, which in part is due to the sparseness of the data. The correlation structure for the predictor variables is not well retained on multiply-imputed data from small samples with more than 50% missing data with this model. For mixed continuous and dichotomous predictor data, the results are similar to those found under the multivariate normal model for continuous data and under the multinomial model for dichotomous data. With all data types, a fully-observed variable included with variables subject to missingness in the multiple imputation process and subsequent statistical analysis provided liberal (larger than nominal values) Type I error rates under a specific pattern of missing data. It is suggested that future studies focus on the effects of multiple imputation in multivariate settings with more realistic data characteristics and a variety of multivariate analyses, assessing both Type I error and power. ^