Inherited and noninherited sources of variation in metastasis and the growth rate of tumors: Implications for tumor evolution and therapy

I have undertaken measurements of the genetic (or inherited) and nongenetic (or noninherited) components of the variability of metastasis formation and tumor diameter doubling time in more than 100 metastatic lines from each of three murine tumors (sarcoma SANH, sarcoma SA4020, and hepatocarcinoma HCA-I) syngeneic to C3Hf/Kam mice. These lines were isolated twice from lung metastases and analysed immediately thereafter to obtain the variance to spontaneous lung metastasis and tumor diameter doubling time. Additional studies utilized cells obtained from within 4 passages of isolation. Under the assumption that no genetic differences in metastasis formation or diameter doubling time existed among the cells of a given line, the variance within a line would estimate nongenetic variation. The variability derived from differences between lines would represent genetic origin. The estimates of the genetic contribution to the variation of metastasis and tumor diameter doubling time were significantly greater than zero, but only in the metastatic lines of tumor SANH was genetic variation the major source of metastatic variability (contributing 53% of the variability). In the tumor cell lines of SA4020 and HCA-I, however, the contribution of nongenetic factors predominated over genetic factors in the variability of the number of metastasis and tumor diameter doubling time. A number of other parameters examined, such as DNA content, karyotype, and selection and variance analysis with passage in vivo, indicated that genetic differences existed within the cell lines and that these differences were probably created by genetic instability. The mean metastatic propensity of the lines may have increased somewhat during their isolation and isotransplantation, but the variance was only slightly affected, if at all. Analysis of the DNA profiles of the metastatic lines of SA4020 and HCA-I revealed differences between these lines and their primary parent tumors, but not among the SANH lines and their parent tumor. Furthermore, there was a direct correlation between the extent of genetic influence on metastasis formation and the ability of the tumor cells to develop resistance to cisplatinum. Thus although nongenetic factors might predominate in contributing to metastasis formation, it is probably genetic variation and genetic instability that cause the progression of tumor cells to a more metastatic phenotype and leads to the emergence of drug resistance. ^

The evolution and assessment of pharmaceutical care: Structure, processes, and outcomes

The evolution of pharmaceutical care is identified through a complete review of the literature published in the American Journal of Health-System Pharmacy, the sole comprehensive publication of institutional pharmacy practice. The evolution is categorized according to characteristics of structure (organizational structure, the role of the pharmacist), process (drug delivery systems, formulary management, acquiring drug products, methods to impact drug therapy decisions), and outcomes (cost of drug delivery, cost of drug acquisition and use, improved safety, improved health outcomes) recorded from the 1950s through the 1990s. While significant progress has been made in implementing basic drug distribution systems, levels of pharmacy involvement with direct patient care is still limited.^ A new practice framework suggests enhanced direct patient care involvement through increase in the efficiency and effectiveness of traditional pharmacy services. Recommendations advance internal and external organizational structure relationships that position pharmacists to fully use their unique skills and knowledge to impact drug therapy decisions and outcomes. Specific strategies facilitate expansion of the breadth and scope of each process component in order to expand the depth of integration of pharmacy and pharmaceutical care within the broad healthcare environment. Economic evaluation methods formally evaluate the impact of both operational and clinical interventions.^ Outcome measurements include specific recommendations and methods to increase efficiency of drug acquisition, emphasizing pharmacists' roles that impact physician prescribing decisions. Effectiveness measures include those that improve safety of drug distribution systems, decrease the potential of adverse drug therapy events, and demonstrate that pharmaceutical care can significantly contribute to improvement in overall health status.^ The implementation of the new framework is modeled on a case study at the M.D. Anderson Cancer Center. The implementation of several new drug distribution methods facilitated the redeployment of personnel from distributive functions to direct patient care activities with significant personnel and drug cost reduction. A cost-benefit analysis illustrates that framework process enhancements produced a benefit-to-cost ratio of 7.9. In addition, measures of effectiveness demonstrated significant levels of safety and enhanced drug therapy outcomes. ^

BMP-SIGNALING REGULATES A COMMON TRANSCRIPTIONAL PROGRAM TO CONTROL FACIAL FORM AND SKELETAL MORPHOGENESIS

Much of the craniofacial skeleton, such as the skull vault, mandible and midface, develops through direct, intramembranous ossification of the cranial neural crest (CNC) derived progenitor cells. Bmp-signaling plays critical roles in normal craniofacial development, and Bmp4 deficiency results in craniofacial abnormalities, such as cleft lip and palate. We performed an in depth analysis of Bmp4, a critical regulator of development, disease, and evolution, in the CNC. Conditional Bmp4 overexpression, using a tetracycline regulated Bmp4 gain of function allele, resulted in facial form changes that were most dramatic after an E10.5 Bmp4 induction. Expression profiling uncovered a signature of Bmp4 induced genes (BIG) composed predominantly of transcriptional regulators controlling self-renewal, osteoblast differentiation, and negative Bmp autoregulation. The complimentary experiment, CNC inactivation of Bmp2, Bmp4, and Bmp7, resulted in complete or partial loss of multiple CNC derived skeletal elements revealing a critical requirement for Bmp-signaling in membranous bone and cartilage development. Importantly, the BIG signature was reduced in Bmp loss of function mutants indicating similar Bmp-regulated target genes underlying facial form modulation and normal skeletal morphogenesis. Chromatin immunoprecipitation (ChIP) revealed a subset of the BIG signature, including Satb2, Smad6, Hand1, Gadd45g and Gata3 that was bound by Smad1/5 in the developing mandible revealing direct, Smad-mediated regulation. These data indicate that Bmp-signaling regulates craniofacial skeletal development and facial form by balancing self-renewal and differentiation pathways in CNC progenitors.

Population and molecular evolution studies on the Melanocortin 1 receptor locus implicate its role in human pigmentation variation

Human pigmentation is a complex trait with the observed variation caused by the varied production of eumelanin (brown/black melanins) and phaeomelanin (red/yellow melanins) by the melanocytes. The melanocortin 1 receptor (MC1R), a G protein-coupled receptor expressed in the melanocytes, is a regulator eu- and phaeomelanin synthesis, and MC1R mutations causing skin and coat color changes are known in many mammals. To understand the role of MC1R in human pigmentation variation, I have sequenced the MC1R gene in 121 individuals sampled from world populations. In addition, I have sequenced the MC1R gene in common and pygmy chimpanzees, gorilla, orangutan, and baboon to study the evolution of MC1R and to infer the ancestral human MC1R sequence. The ancestral MC1R sequence is observed in all 25 African individuals studied, but at lower frequencies in the other populations examined, especially in East and Southeast Asians. The Arg163Gln variant is absent in the Africans studied, almost absent in Europeans, and at a low frequency in Indians, but is at an exceptionally high frequency (70%) in East and Southeast Asians. To further evaluate the role of MC1R variants in human pigmentation variation, I have combined these molecular evolution and population studies with functional assays on MC1R variants and primate MC1Rs. ^

Population genetics of VNTR loci and their applications in evolutionary studies

Variable number of tandem repeats (VNTR) are genetic loci at which short sequence motifs are found repeated different numbers of times among chromosomes. To explore the potential utility of VNTR loci in evolutionary studies, I have conducted a series of studies to address the following questions: (1) What are the population genetic properties of these loci? (2) What are the mutational mechanisms of repeat number change at these loci? (3) Can DNA profiles be used to measure the relatedness between a pair of individuals? (4) Can DNA fingerprint be used to measure the relatedness between populations in evolutionary studies? (5) Can microsatellite and short tandem repeat (STR) loci which mutate stepwisely be used in evolutionary analyses?^ A large number of VNTR loci typed in many populations were studied by means of statistical methods developed recently. The results of this work indicate that there is no significant departure from Hardy-Weinberg expectation (HWE) at VNTR loci in most of the human populations examined, and the departure from HWE in some VNTR loci are not solely caused by the presence of population sub-structure.^ A statistical procedure is developed to investigate the mutational mechanisms of VNTR loci by studying the allele frequency distributions of these loci. Comparisons of frequency distribution data on several hundreds VNTR loci with the predictions of two mutation models demonstrated that there are differences among VNTR loci grouped by repeat unit sizes.^ By extending the ITO method, I derived the distribution of the number of shared bands between individuals with any kinship relationship. A maximum likelihood estimation procedure is proposed to estimate the relatedness between individuals from the observed number of shared bands between them.^ It was believed that classical measures of genetic distance are not applicable to analysis of DNA fingerprints which reveal many minisatellite loci simultaneously in the genome, because the information regarding underlying alleles and loci is not available. I proposed a new measure of genetic distance based on band sharing between individuals that is applicable to DNA fingerprint data.^ To address the concern that microsatellite and STR loci may not be useful for evolutionary studies because of the convergent nature of their mutation mechanisms, by a theoretical study as well as by computer simulation, I conclude that the possible bias caused by the convergent mutations can be corrected, and a novel measure of genetic distance that makes the correction is suggested. In summary, I conclude that hypervariable VNTR loci are useful in evolutionary studies of closely related populations or species, especially in the study of human evolution and the history of geographic dispersal of Homo sapiens. (Abstract shortened by UMI.) ^

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

INDUCTION OF DNA STRAND BREAKS AND CROSS-LINKS BY THE NEW ANTICANCER AGENT 3,6-DIAZIRIDINYL-2,5-BIS(CARBOETHOXYAMINO) -1,4-BENZOQUINONE

The DNA breakage effect of the anticancer agent 3,6-diaziridinyl-2,5-bis(carboethoxyamino)-1,4-benzoquinone (AZQ, NSC-182986) on bacteriophage PM2 DNA was investigated using agarose gel electrophoresis. AZQ caused both single-stranded and double-stranded breaks after reduction with NaBH(,4), but it was not active in the native state. At 120 (mu)M, it degraded 50% of the closed circular form I DNA into 40% form II DNA (single-stranded break) and 10% form III DNA (double-stranded break). It produced a dose-response breakage between 1 (mu)M and 320 (mu)M. The DNA breakage exhibited a marked pH dependency. At 320 (mu)M, AZQ degraded 80% and 60% of form I DNA at pH 4 and 10 respectively, but none between pH 6 to 8. The DNA breakage at physiologic pH was greatly enhanced when 10 (mu)M cupric sulfate was included in the incubation mixture. The DNA strand scission was inhibited by catalase, glutathione, KI, histidine, Tiron, and DABCO. These results suggest that the DNA breakage may be caused by active oxygen metabolites including hydroxyl free radical. The bifunctional cross-linking activity of reduced AZQ on isolated calf thymus DNA was investigated by ethidium fluorescence assay. The cross-linking activity exhibited a similar pH dependency; highest in acidic and alkaline pH, inactive under neutral conditions. Using the alkaline elution method, we found that AZQ induced DNA single-stranded breaks in Chinese hamster ovary cells treated with 50 (mu)M of AZQ for 2 hr. The single-stranded break frequencies in rad equivalents were 17 with 50 (mu)M and 140 with 100 (mu)M of AZQ. In comparison, DNA cross-links appeared in cells treated with only 1 to 25 (mu)M of AZQ for 2 hr. The cross-linking frequencies in rad equivalents were 39 and 90 for 1 and 5 (mu)M of AZQ, respectively. Both DNA-DNA and DNa-protein cross-links were induced by AZQ in CHO cells as revealed by the proteinas K digestion assay. DNA cross-links increased within the first 4 hr of incubation in drug-free medium and slightly decreased by 12 hr, and most of the cross-links disappeared after cells were allowed to recovered for 24 hr.^ By electrochemical analysis, we found that AZQ was more readily reduced at acidic pH. However, incubation of AZQ with NaBH(,4) at pH 7.8 or 10, but not at 4, produced superoxide anion. The opening of the aziridinyl rings of AZQ at pH 4 was faster in the presence of NaBH(,4) than in its absence; no ring-opening was detected at pH 7.8 regardless of the inclusion of NaBH(,4). . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^