A SIMULATION STUDY OF THE SIZE AND POWER OF SOME TWO-SAMPLE TESTS FOR UNEQUALLY CENSORED SURVIVAL DATA (RANK TESTS, LOGRANK, PROPORTIONAL HAZARDS, WILCOXON TESTS, EXPONENTIAL DISTRIBUTION)

Sizes and power of selected two-sample tests of the equality of survival distributions are compared by simulation for small samples from unequally, randomly-censored exponential distributions. The tests investigated include parametric tests (F, Score, Likelihood, Asymptotic), logrank tests (Mantel, Peto-Peto), and Wilcoxon-Type tests (Gehan, Prentice). Equal sized samples, n = 18, 16, 32 with 1000 (size) and 500 (power) simulation trials, are compared for 16 combinations of the censoring proportions 0%, 20%, 40%, and 60%. For n = 8 and 16, the Asymptotic, Peto-Peto, and Wilcoxon tests perform at nominal 5% size expectations, but the F, Score and Mantel tests exceeded 5% size confidence limits for 1/3 of the censoring combinations. For n = 32, all tests showed proper size, with the Peto-Peto test most conservative in the presence of unequal censoring. Powers of all tests are compared for exponential hazard ratios of 1.4 and 2.0. There is little difference in power characteristics of the tests within the classes of tests considered. The Mantel test showed 90% to 95% power efficiency relative to parametric tests. Wilcoxon-type tests have the lowest relative power but are robust to differential censoring patterns. A modified Peto-Peto test shows power comparable to the Mantel test. For n = 32, a specific Weibull-exponential comparison of crossing survival curves suggests that the relative powers of logrank and Wilcoxon-type tests are dependent on the scale parameter of the Weibull distribution. Wilcoxon-type tests appear more powerful than logrank tests in the case of late-crossing and less powerful for early-crossing survival curves. Guidelines for the appropriate selection of two-sample tests are given. ^

Effect of patient accrual patterns on power and size of log-rank test for time-to-event outcome in clinical trials

The determination of size as well as power of a test is a vital part of a Clinical Trial Design. This research focuses on the simulation of clinical trial data with time-to-event as the primary outcome. It investigates the impact of different recruitment patterns, and time dependent hazard structures on size and power of the log-rank test. A non-homogeneous Poisson process is used to simulate entry times according to the different accrual patterns. A Weibull distribution is employed to simulate survival times according to the different hazard structures. The current study utilizes simulation methods to evaluate the effect of different recruitment patterns on size and power estimates of the log-rank test. The size of the log-rank test is estimated by simulating survival times with identical hazard rates between the treatment and the control arm of the study resulting in a hazard ratio of one. Powers of the log-rank test at specific values of hazard ratio (≠1) are estimated by simulating survival times with different, but proportional hazard rates for the two arms of the study. Different shapes (constant, decreasing, or increasing) of the hazard function of the Weibull distribution are also considered to assess the effect of hazard structure on the size and power of the log-rank test. ^

NOVEL PHANTOMS AND POST-PROCESSING FOR DIFFUSION SPECTRUM IMAGING

High Angular Resolution Diffusion Imaging (HARDI) techniques, including Diffusion Spectrum Imaging (DSI), have been proposed to resolve crossing and other complex fiber architecture in the human brain white matter. In these methods, directional information of diffusion is inferred from the peaks in the orientation distribution function (ODF). Extensive studies using histology on macaque brain, cat cerebellum, rat hippocampus and optic tracts, and bovine tongue are qualitatively in agreement with the DSI-derived ODFs and tractography. However, there are only two studies in the literature which validated the DSI results using physical phantoms and both these studies were not performed on a clinical MRI scanner. Also, the limited studies which optimized DSI in a clinical setting, did not involve a comparison against physical phantoms. Finally, there is lack of consensus on the necessary pre- and post-processing steps in DSI; and ground truth diffusion fiber phantoms are not yet standardized. Therefore, the aims of this dissertation were to design and construct novel diffusion phantoms, employ post-processing techniques in order to systematically validate and optimize (DSI)-derived fiber ODFs in the crossing regions on a clinical 3T MR scanner, and develop user-friendly software for DSI data reconstruction and analysis. Phantoms with a fixed crossing fiber configuration of two crossing fibers at 90° and 45° respectively along with a phantom with three crossing fibers at 60°, using novel hollow plastic capillaries and novel placeholders, were constructed. T2-weighted MRI results on these phantoms demonstrated high SNR, homogeneous signal, and absence of air bubbles. Also, a technique to deconvolve the response function of an individual peak from the overall ODF was implemented, in addition to other DSI post-processing steps. This technique greatly improved the angular resolution of the otherwise unresolvable peaks in a crossing fiber ODF. The effects of DSI acquisition parameters and SNR on the resultant angular accuracy of DSI on the clinical scanner were studied and quantified using the developed phantoms. With a high angular direction sampling and reasonable levels of SNR, quantification of a crossing region in the 90°, 45° and 60° phantoms resulted in a successful detection of angular information with mean ± SD of 86.93°±2.65°, 44.61°±1.6° and 60.03°±2.21° respectively, while simultaneously enhancing the ODFs in regions containing single fibers. For the applicability of these validated methodologies in DSI, improvement in ODFs and fiber tracking from known crossing fiber regions in normal human subjects were demonstrated; and an in-house software package in MATLAB which streamlines the data reconstruction and post-processing for DSI, with easy to use graphical user interface was developed. In conclusion, the phantoms developed in this dissertation offer a means of providing ground truth for validation of reconstruction and tractography algorithms of various diffusion models (including DSI). Also, the deconvolution methodology (when applied as an additional DSI post-processing step) significantly improved the angular accuracy of the ODFs obtained from DSI, and should be applicable to ODFs obtained from the other high angular resolution diffusion imaging techniques.

Expression and function of $\beta$-tubulin isotypes in Chinese hamster ovary (CHO) cells

The availability of isotype specific antisera for $\beta$-tubulin, coupled with genetic and biochemical analysis, has allowed the determination of $\beta$-tubulin isotype expression and distribution in Chinese hamster ovary (CHO) cells. Using genetic manipulations involving selection for colcemid resistance followed by reversion and reselection for drug resistance, we have succeeded in isolating cell lines that exhibit three major and one minor $\beta$-tubulin spots by two-dimensional gel electrophoresis. In concert with isotype specific antibodies, analysis of these mutants demonstrates that CHO cells express two copies of isotype I, at least one copy of isotype IV, and very small amounts of isotype V. Their stoichiometry is approximately 1:1:0.7:0.2. All three isotypes assemble into both cytoplasmic and spindle microtubules, and are similar in their responses to cold, colcemid, and calcium induced depolymerization. They have comparable turnover rates and are equally sensitive to depression of synthesis upon colchicine treatment. These results suggest that $\beta$-tubulin isotypes are used interchangeably to assemble microtubule structures in CHO cells. However, of 18 colcemid resistant mutants with a demonstrable alteration in $\beta$-tubulin, all were found to have the alteration in isotype I, thus leaving open the possibility that subtle differences in isotype properties may exist. Under various conditions of the cell growth, the relative proportion of each expressed isotype does not significantly seem to change except in the early G1 phase of the cell cycle. At this time the synthesis of isotype V increases more than two fold relative to isotype I and IV, while at the same time, total $\beta$-tubulin synthesis is decreased about 60-70%. ^

Cell surface galactosyltransferase structure and function during mesenchymal cell migration

In order to more fully understand the function of surface GalTase on mesenchymal cells, anti-GalTase IgG was used to (a) examine the role of surface GalTase during mouse mesenchymal cell migration on laminin and fibronectin; (b) define the plasma membrane distribution of GalTase by indirect immunofluorescence on migrating cells; (c) quantitate the level of surface GalTase on migrating cells; and (d) determine whether GalTase is associated with the cytoskeleton.^ Results show that anti-GalTase IgG was able to inhibit migration (48-80% as compared to basal rate) when cells were migrating on laminin-containing matrices. Monovalent Fab fragments inhibited migration on laminin by 90% after 4 hours. On the other hand, anti-GalTase IgG had no effect on cells migrating on fibronectin. This illustrates the substrate specificity of GalTase mediated-migration. When anti-GalTase IgG was used to localize surface GalTase on cells migratory on laminin, the enzyme was restricted to the leading and trailing edges of the cell. Assays indicate that GalTase is elevated approximately 3-fold when cells are migrating on laminin-containing matrices as compared to migratory cells on plastic or fibronectin, or as compared to stationary cells on any substrate. Laminin appears to recruit GalTase from preexisting intracellular pools to the growing lamellipodia.^ Double-label indirect immunofluorescence studies indicate that there is an apparent co-localization between some of the surface GalTase and some actin filaments. This relationship was explored by extracting cells prelabeled with anti-GalTase IgG and quantitated by radiolabeled second antibodies. Results show that 79% of the surface GalTase is associated with the cytoskeleton (as judged by detergent insolubility) when monovalent antibodies (Fab) are used. However virtually all (80-100%) of the surface GalTase can be induced to associate with the cytoskeleton when cross-linked with bivalent antibodies. Furthermore, when cells in suspension are incubated with divalent antibodies, an additional 66% of the surface GalTase can be induced to associate with the cytoskeleton. The elevated levels of surface GalTase detectable on cells migrating on laminin also appear to be associated with the cytoskeleton.^ Several lines of evidence suggest that GalTase is associated with F-actin. Data suggest that laminin induces the expression of surface GalTase to the growing lamellipodia where it becomes associated with the cytoskeleton leading to cell spreading and migration. (Abstract shortened with permission of author.) ^

Cytochrome P450 in mammalian brain: Purification, anatomical distribution and regulation of the cytochrome P450 monooxygenase system in rat brain

The cytochrome P450 (P450) monooxygenase system plays a major role in metabolizing a wide variety of xenobiotic as well as endogenous compounds. In performing this function, it serves to protect the body from foreign substances. However, in a number of cases, P450 activates procarcinogens to cause harm. In most animals, the highest level of activity is found in the liver. Virtually all tissues demonstrate P450 activity, though, and the role of the P450 monooxygenase system in these other organs is not well understood. In this project I have studied the P450 system in rat brain; purifying NADPH-cytochrome P450 reductase (reductase) from that tissue. In addition, I have examined the distribution and regulation of expression of reductase and P450 in various anatomical regions of the rat brain.^ NADPH-cytochrome P450 reductase was purified to apparent homogeneity and cytochrome P450 partially purified from whole rat brain. Purified reductase from brain was identical to liver P450 reductase by SDS-PAGE and Western blot techniques. Kinetic studies utilizing cerebral P450 reductase reveal Km values in close agreement with those determined with enzyme purified from rat liver. Moreover, the brain P450 reductase was able to function successfully in a reconstituted microsomal system with partially purified brain cytochrome P450 and with purified hepatic P4501A1 as measured by 7-ethoxycoumarin and 7-ethoxyresorufin O-deethylation. These results indicate that the reductase and P450 components may interact to form a competent drug metabolism system in brain tissue.^ Since the brain is not a homogeneous organ, dependent upon the well orchestrated interaction of numerous parts, pathology in one nucleus may have a large impact upon its overall function. Hence, the anatomical distribution of the P450 monooxygenase system in brain is important in elucidating its function in that organ. Related to this is the regulation of P450 expression in brain. In order to study these issues female rats--both ovariectomized and not--were treated with a number of xenobiotic compounds and sex steroids. The brains from these animals were dissected into 8 discrete regions and the presence and relative level of message for P4502D and reductase determined using polymerase chain reaction. Results of this study indicate the presence of mRNA for reductase and P4502D isoforms throughout the rat brain. In addition, quantitative PCR has allowed the determination of factors affecting the expression of message for these enzymes. ^

Structure and function study of prostaglandin H synthase

Prostaglandin H synthase (PGHS) is a key enzyme in biosynthesis of prostaglandins, thromboxane, and prostacyclin. It has two activities, cyclooxygenase and peroxidase. "PGHS" means PGHS-1. A current hypothesis considers the cyclooxygenase reaction to be a free radical chain reaction, initiated by interaction of the synthase peroxidase with hydroperoxides leading to the production of a tyrosyl free radical. According to this hypothesis, tyrosyl residue(s) may play a key role in the cyclooxygenase reaction. Tetranitromethane (TNM) can relatively selectively nitrate tyrosines at pH 8.0. The effect of TNM on both cyclooxygenase activity and peroxidase activity has been examined: reaction of the synthase holoenzyme with TNM at pH 8.0 led to inactivation of both activities, with the cyclooxygenase activity being lost rapidly and completely, while the peroxidase activity was lost more slowly. Indomethacin, a non-steroidal anti-inflammatory agent, can protect the synthase from the inactivation of TNM. Amino acid analyses indicated that a loss of tyrosine and formation of nitrotyrosine residues occurred during reaction with TNM, and that TNM-reacted holoenzyme with $<$10% residual cyclooxygenase activity had about 2.0 nitrotyrosine/subunit.^ PGH synthase is known to be an endoplasmic reticulum membrane-associated protein. Antibodies directed at particular PGHS peptide segments and indirect immunofluorescence have been used to characterize the membrane topology of crucial portions of PGHS. PGHS was expressed in COS-1 cells transfected with the appropriate cDNA. Stably-transfected human endothelial cells were also used for the topology study. The cells were treated with streptolysin-O, which selectively permeabilizes the plasma membrane, or with saponin to achieve general membrane disruption, before incubation with the antipeptide antibodies. Bound antipeptide antibody was stained by FITC-labelled secondary antibody and visualized by fluorescence microscopy. With the antipeptide antibodies against residues 51-66, 156-170 or 377-390, there was a significant reticular and perinuclear pattern of staining in cells permeabilized with saponin but not in cells permeabilized with SLO alone. Antibodies directed against the endogenous C-terminal peptide or against residues 271-284 produced staining in cells permeabilized with saponin, and also in a lower, but significant fraction of cells permeabilized with SLO. Similar results were obtained when COS-1 cells expressing recombinant PGHS with a viral reporter peptide inserted at the C-terminus were stained with antibody against the reporter epitope.^ The PGHS C-terminal sequence is similar to that of the consensus KDEL ER retention signal. The potential function of the PGHS C-terminus segment in ER retention was examined by mutating this segment and analyzing the subcellular distribution of the mutants expressed in COS-1 cells. None of the mutants had an altered subcellular distribution, although some had greatly diminished the enzyme activities. (Abstract shortened by UMI.) ^

Effects of informative censoring on the hazard function: Application to the proportional hazards regression model

The standard analyses of survival data involve the assumption that survival and censoring are independent. When censoring and survival are related, the phenomenon is known as informative censoring. This paper examines the effects of an informative censoring assumption on the hazard function and the estimated hazard ratio provided by the Cox model.^ The limiting factor in all analyses of informative censoring is the problem of non-identifiability. Non-identifiability implies that it is impossible to distinguish a situation in which censoring and death are independent from one in which there is dependence. However, it is possible that informative censoring occurs. Examination of the literature indicates how others have approached the problem and covers the relevant theoretical background.^ Three models are examined in detail. The first model uses conditionally independent marginal hazards to obtain the unconditional survival function and hazards. The second model is based on the Gumbel Type A method for combining independent marginal distributions into bivariate distributions using a dependency parameter. Finally, a formulation based on a compartmental model is presented and its results described. For the latter two approaches, the resulting hazard is used in the Cox model in a simulation study.^ The unconditional survival distribution formed from the first model involves dependency, but the crude hazard resulting from this unconditional distribution is identical to the marginal hazard, and inferences based on the hazard are valid. The hazard ratios formed from two distributions following the Gumbel Type A model are biased by a factor dependent on the amount of censoring in the two populations and the strength of the dependency of death and censoring in the two populations. The Cox model estimates this biased hazard ratio. In general, the hazard resulting from the compartmental model is not constant, even if the individual marginal hazards are constant, unless censoring is non-informative. The hazard ratio tends to a specific limit.^ Methods of evaluating situations in which informative censoring is present are described, and the relative utility of the three models examined is discussed. ^

The assumptions and effects of conservative procedures in low -dose risk assessment for setting safety standards

Conservative procedures in low-dose risk assessment are used to set safety standards for known or suspected carcinogens. However, the assumptions upon which the methods are based and the effects of these methods are not well understood.^ To minimize the number of false-negatives and to reduce the cost of bioassays, animals are given very high doses of potential carcinogens. Results must then be extrapolated to much smaller doses to set safety standards for risks such as one per million. There are a number of competing methods that add a conservative safety factor into these calculations.^ A method of quantifying the conservatism of these methods was described and tested on eight procedures used in setting low-dose safety standards. The results using these procedures were compared by computer simulation and by the use of data from a large scale animal study.^ The method consisted of determining a "true safe dose" (tsd) according to an assumed underlying model. If one assumed that Y = the probability of cancer = P(d), a known mathematical function of the dose, then by setting Y to some predetermined acceptable risk, one can solve for d, the model's "true safe dose".^ Simulations were generated, assuming a binomial distribution, for an artificial bioassay. The eight procedures were then used to determine a "virtual safe dose" (vsd) that estimates the tsd, assuming a risk of one per million. A ratio R = ((tsd-vsd)/vsd) was calculated for each "experiment" (simulation). The mean R of 500 simulations and the probability R $<$ 0 was used to measure the over and under conservatism of each procedure.^ The eight procedures included Weil's method, Hoel's method, the Mantel-Byran method, the improved Mantel-Byran, Gross's method, fitting a one-hit model, Crump's procedure, and applying Rai and Van Ryzin's method to a Weibull model.^ None of the procedures performed uniformly well for all types of dose-response curves. When the data were linear, the one-hit model, Hoel's method, or the Gross-Mantel method worked reasonably well. However, when the data were non-linear, these same methods were overly conservative. Crump's procedure and the Weibull model performed better in these situations. ^

Expression and function of $\beta$-1,4-galactosyltransferase during wild-type and T/{\it t\/} -mutant spermatogenesis

In the mouse, gamete recognition is mediated in part by the binding of sperm surface $\beta$1,4 galactosyltransferase (GalTase) to specific oligosaccharide residues on the zona pellucida ZP3. The expression of GalTase on the sperm surface is regulated by alleles within the distal segment of the T/t complex and results in a haploid-specific increase in GalTase expression on spermatids and sperm from t-bearing males, suggesting that differences in sperm GalTase activity may contribute to t-sperm transmission ratio distortion. In this study, the expression of GalTase RNA during wild-type and T/t-mutant spermatogenesis was characterized and the role of GalTase was analyzed in transmission ratio distortion. It was found that spermatogenic cells predominantly express the long form of the GalTase RNA, which encodes the GalTase protein that is preferentially targeted to the cell surface in somatic cells. In wild-type testes, GalTase RNA accumulates during the maturation of primary spermatocytes, reaches peak levels prior to meiosis, and decreases and meiosis. GalTase RNA accumulates to similar levels during the maturation of +/t and t/t primary spermatocytes, but unlike wild-type, the level of GalTase RNA in t-spermatocytes remains elevated during meiotic division. Consequently, spermatids in t-mutant testes inherit higher levels of GalTase RNA than do wild-type spermatids, which likely accounts for the haploid-specific increase in surface GalTase activity characteristic of spermatids from t-bearing mice.^ The functional significance of the increased GalTase activity during t-sperm transmission ratio distortion was determined by examining the distribution of GalTase RNA and surface GalTase protein in haploid spermatids from +/t males. Results show that +- and t-spermatids have similar levels of both GalTase RNA and protein, indicating that transmission ratio distortion in +/t mice is not likely due to haploid-specific differences in sperm surface GalTase activity.^ The presence of GalTase on the surface of an early spermatogenic cells before it is required on the mature sperm to perform its function during gamete binding suggests a separate function for GalTase in Sertoli-germ cell adhesion. Studies indicate that cell surface GalTase partly mediates the initial adhesion of pachytene spermatocytes, but not haploid spermatids, to Sertoli cells. ^