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

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

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

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

The transcriptional regulation of the {\it neu\/} gene: Examples of activation, repression, and cell-specific expression

The neu oncogene encodes a growth factor receptor-like protein, p185, with an intrinsic tyrosine kinase activity. A single point mutation, an A to T transversion resulting in an amino acid substitution from valine to glutamic acid, in the transmembrane domain of the rat neu gene was found to be responsible for the transforming and tumorigenic phenotype of the cells that carry it. In contrast, the human proto-neu oncogene is frequently amplified in tumors and cell lines derived from tumors and the human neu gene overexpression/amplification in breast and ovarian cancers is known to correlate with poor patient prognosis. Examples of the human neu gene overexpression in the absence of gene amplification have been observed, which may suggest the significant role of the transcriptional and/or post-transcriptional control of the neu gene in the oncogenic process. However, little is known about the transcriptional mechanisms which regulate the neu gene expression. In this study, three examples are presented to demonstrate the positive and negative control of the neu gene expression.^ First, by using band shift assays and methylation interference analyses, I have identified a specific protein-binding sequence, AAGATAAAACC ($-$466 to $-$456), that binds a specific trans-acting factor termed RVF (for EcoRV factor on the neu promoter). The RVF-binding site is required for maximum transcriptional activity of the rat neu promoter. This same sequence is also found in the corresponding regions of both human and mouse neu promoters. Furthermore, this sequence can enhance the CAT activity driven by a minimum promoter of the thymidine kinase gene in an orientation-independent manner, and thus it behaves as an enhancer. In addition, Southwestern (DNA-protein) blot analysis using the RVF-binding site as a probe points to a 60-kDa polypeptide as a potential candidate for RVF.^ Second, it has been reported that the E3 region of adenovirus 5 induces down-regulation of epidermal growth factor (EGF) receptor through endocytosis. I found that the human neu gene product, p185, (an EGF receptor-related protein) is also down-regulated by adenovirus 5, but via a different mechanism. I demonstrate that the adenovirus E1a gene is responsible for the repression of the human neu gene at the transcriptional level.^ Third, a differential expression of the neu gene has been found in two cell model systems: between the mouse fibroblast Swiss-Webster 3T3 (SW3T3) and its variant NR-6 cells; and between the mouse liver tumor cell line, Hep1-a, and the mouse pancreas tumor cell line, 266-6. Both NR-6 and 266-6 cell lines are not able to express the neu gene product, p185. I demonstrate that, in both cases, the transcriptional repression of the neu gene may account for the lack of the p185 expression in these two cell lines. ^

Characterization of the MLC box: A conserved {\it cis\/}-acting DNA element present in the chicken MLC1s/1c gene

The slow/cardiac alkali myosin light chain (MLC1s/1c) is a member of a multigene family whose protein products are essential for activation of the myosin ATPase. In the adult, the MLC1s/1c isoform is expressed in both cardiac and slow-twitch skeletal muscles, while it is expressed by all skeletal muscles during development.^ To elucidate the molecular mechanisms that underlie the transcriptional regulation of MLC1s/1c gene expression, the immediate 5$\sp\prime$ flanking region of the gene was isolated and shown to be capable of directing reporter gene expression. Analysis of this region revealed a 110 bp muscle-specific enhancer that includes a myocyte-specific enhancer-binding factor 2 (MEF-2) site, E-boxes, which are potential binding sites for the basic-helix-loop-helix proteins such as MyoD, and a MLC box. The focus of the thesis was to identify the role of the MLC box in expression of the MLC1s/1c gene.^ The MLC box is a member of the family of CArG box containing cis-acting DNA elements. Mutagenesis showed that the MLC box is necessary, but not sufficient, for the expression of a reporter gene linked to the 5$\sp\prime$ flanking region of the MLC1s/1c gene. Linker scanner and site-directed mutagenesis identified a number of potential sites within the 110 bp muscle-specific enhancer that may cooperate with the MLC box. These are the MEF-2 site, the E-box site, and a 10 bp element located upstream of the MEF-2 site that does not have sequence similarity with any known cis-acting element. The MLC box is capable of binding to factors present in muscle nuclear extracts, as well as to human recombinant serum response factor (SRF). Binding of SRF to the MLC box was correlated with the ability of the 5$\sp\prime$ flanking region of the MLC1s/1c gene to drive reporter gene expression. Results suggest a model in which binding of SRF to the MLC box activates expression of the MLC1s/1c gene while binding of the factors present in the nuclear extracts suppresses the expression of the gene. (Abstract shortened with permission of author.) ^

Alterations in the ras oncogene and the p53 tumor suppressor gene in ultraviolet radiation-induced skin carcinogenesis

A rapid increase of the ultraviolet radiation (UVR)-related skin cancer incidence has attracted more and more public attention during the last few decades. Prevention and treatment of UVR-related skin cancer has become an important public health issue in the United States. Recent studies indicate that mutations in ras and/or p53 genes may be involved in UVR-induced skin tumor development but the precise molecular mechanism remains unclear. In this study, alterations of H-ras and p53 genes were investigated in different stages of carcinogenesis in a chronic UVR (solar simulator) exposure-induced Sencar mouse skin carcinogenesis model in order to clarify the role of the alterations of these genes during the skin carcinogenesis process and to further understand the mechanisms by which UVR causes skin cancer.^ Positive ras-p21 staining in cell membranes and cytosol were detected in 18/33 (55%) of squamous cell carcinomas (SCCs), but were not detected in UV-exposed skin, papillomas, or spindle cell tumors (SCTs). Positive staining of the malignant progression marker K13 was found in 17/33 (52%) of SCCs only. A significant positive correlation was observed between the K13 and the ras-p21 expression. Polymerase chain reaction (PCR)-based single strand conformation polymorphism (SSCP) analysis and gene sequencing analysis revealed three point mutations, one (codon 56) in UV-exposed non-tumor bearing skin and the other two (codons 21 and 13) in SCCs. No UV-specific mutation patterns were found.^ Positive p53 nuclear staining was found in 10/37 (27%) of SCCs and 12/24 (50%) of SCTs, but was not detected in normal skin or papillomas. PCR-based SSCP and sequencing analysis revealed eight point mutations in exons 5 and 6 (four in SCTs, two in SCCs, and two in UV-exposed skin) including six C-T or C-A transitions. Four of the mutations occurred at a dipyrimidine (CC) sequence. The pattern of the mutations indicated that the mutagenic lesions were induced by UVR.^ These results indicate that overexpression of ras-p21 in conjunction with aberrant expression of K13 occurred frequently in UVR-induced SCCs in Sencar mouse skin. The point mutation in the H-ras gene appeared to be a rare event in UVR skin carcinogenesis and may not be responsible for overexpression of ras-p21. UVR-induced P53 gene alteration is a frequent event in UVR-induced SCCs and later stage SCT tumors in Sencar mice skin, suggesting the p53 gene mutation plays an important role in skin tumor malignant progression. ^

Transcriptional regulation of the rat {\it mdr1b\/} gene

The expression of P-glycoproteins encoded by the mdr gene family is associated with the emergence of multidrug-resistance phenotype in animal cells. This gene family includes two members, MDR1 and MDR2, in humans, and three members, mdr1a, mdr1b, and mdr2, in rodents. Among them, the rat mdr1b is known to be highly activated during hepatocarcinogenesis, and its expression is sensitive to the treatment with growth factors, cytotoxic drugs, as well as other physical or chemical stresses. It is believed that the transcriptional regulation plays an important role in above events, however little has been known about mechanisms involved.^ To elucidate how mdr1b expression is regulated, we isolated the genomic sequence of the rat mdr1b and functionally dissected its 5$\prime$ promoter region. Our results demonstrated that: (1) the transcription start site of the rat mdr1b is identical to that of the murine mdr1b homologue; (2) a palindromic sequence from bp $-$189 to $-$180 bp is essential for the basal promoter function of the rat mdr1b, and binds to a specific protein that appears to be a novel transcription factor implicated in the regulation of the rat mdr1b expression; (3) a NF-$\kappa$B-binding site from bp $-$167 to $-$159 is also involved in the basal promoter function. The p65/p50 subunits of the NF-$\kappa$B and raf-1 kinase are implicated in the insulin-inducible promoter activity of the mdr1b, suggesting the important role of NF-$\kappa$B in the regulation of the mdr1b by growth factors; (4) a p53-binding site from bp $-$199 to $-$180 is not only essential for the basal promoter activity but also responsible for the induction of mdr1b by cytotoxic agents. In addition, we provided evidence showing that endogenous mdr1b expression can be modulated by wild-type p53. On the basis of these findings, a model of transcriptional regulation of the rat mdr1b was proposed. ^

Mechanism of molecular regulation of nuclear -encoded mitochondrial gene expression by electrical stimulation in the cultured neonatal rat cardiac myocytes

To answer the question whether increased energy demand resulting from myocyte hypertrophy and enhanced $\beta$-myosin heavy chain mRNA, contractile protein synthesis and assembly leads to mitochondrial proliferation and differentiation, we set up an electrical stimulation model of cultured neonatal rat cardiac myocytes. We describe, as a result of increased contractile activity, increased mitochondrial profiles, cytochrome oxidase mRNA, and activity, as well as a switch in mitochondrial carnitine palmitoyltransferase-I (CPT-I) from the liver to muscle isoform. We investigate physiological pathways that lead to accumulation of gene transcripts for nuclear encoded mitochondrial proteins in the heart. Cardiomyocytes were stimulated for varying times up to 72 hr in serum-free culture. The mRNA contents for genes associated with transcriptional activation (c-fos, c-jun, junB, nuclear respiratory factor 1 (Nrf-1)), mitochondrial proliferation (cytochrome c (Cyt c), cytochrome oxidase), and mitochondrial differentiation (carnitine palmitonyltransferase I (CPT-I) isoforms) were measured. The results establish a temporal pattern of mRNA induction beginning with c-fos (0.25-3 hr) and followed by c-jun (0.5-3 hr), junB (0.5-6 hr), NRF-1 (1-12 hr), Cyt c (12-72 hr), cytochrome c oxidase (12-72 hr). Induction of the latter was accompanied by a marked decrease in the liver-specific CPT-I mRNA. Electrical stimulation increased c-fos, $\beta$-myosin heavy chain, and Cyt c promoter activities. These increases coincided with a rise in their respective endogenous gene transcripts. NRF-1, cAMP response element (CRE), and Sp-1 site mutations within the Cyt c promoter reduced luciferase expression in both stimulated and nonstimulated myocytes. Mutations in the Nrf-1 and CRE sites inhibited the induction by electrical stimulation or by transfection of c-jun into non-paced cardiac myocytes whereas mutation of the Sp-1 site maintained or increased the fold induction. This is consistent with the appearance of NRF-1 and fos/jun mRNAs prior to that of Cyt c. Overexpression of c-jun by transfection also activates the Nrf-1 and Cyt c mRNA sequentially. Electrical stimulation of cardiac myocytes activates the c-Jun-N-terminal kinase so that the fold-activation of the cyt c promoter is increased by pacing when either c-jun or c-fos/c-jun are cotransfected. We have identified physical association of Nrf-1 protein with the Nrf-1 enhancer element and of c-Jun with the CRE binding sites on the Cyt c promoter. This is the first demonstration that induction of Nrf-1 and c-Jun by pacing of cardiac myocytes directly mediates Cyt c gene expression and mitochondrial proliferation in response to hypertrophic stimuli in the heart.^ Subsequent to gene activation pathways that lead to mitochondrial proliferation, we observed an isoform switch in CPT-I from the liver to muscle mRNA. We have found that the half-life for the muscle CPT-I is not affected by electrical stimulation, but electrical decrease the T1/2 in the liver CPT-I by greater than 50%. This suggests that the liver CPT-I switch to muscle isoform is due to (1) a decrease in T1/2 of liver CPT-I and (2) activation of muscle CPT-Itranscripts by electrical stimulation. (Abstract shortened by UMI.) ^

Genetic integration of a nuclear -encoded mitochondrial gene with cardiac function

Cardiovascular disease (CVD) is the leading cause of death in the United States. One manifestation of CVD known to increase mortality is an enlarged, or hypertrophic heart. Hypertrophic cardiomyocytes adapt to increased contractile demand at the genetic level with a re-emergence of the fetal gene program and a downregulation of fatty acid oxidation genes with concomitant increased reliance on glucose-based metabolism. To understand the transcriptional regulatory pathways that implement hypertrophic directives we analyzed the upstream promoter region of the muscle specific isoform of the nuclear-encoded mitochondrial gene, carnitine palmitoyltransferase-1β (CPT-1β) in cultured rat neonatal cardiac myocytes. This enzyme catalyzes the rate-limiting step of fatty acid entry into β-oxidation and is downregulated in cardiac hypertrophy and failure, making it an attractive model for the study of hypertrophic gene regulation and metabolic adaptations. We demonstrate that the muscle-enriched transcription factors GATA-4 and SRF synergistically activate CPT-1β; moreover, DNA binding to cognate sites and intact protein structure are required. This mechanism coordinates upregulation of energy generating processes with activation of the energy consuming contractile promoter for cardiac α-actin. We hypothesized that fatty acid or glucose responsive transcription factors may also regulate CPT-1β. Oleate weakly stimulates CPT-1β activity; in contrast, the glucose responsive Upstream Stimulatory Factors (USF) dramatically depresses the CPT-1β reporter. USF regulates CPT-1β through a novel physical interaction with the cofactor PGC-1 and abrogation of MEF2A/PGC-1 synergistic stimulation. In this way, USF can inversely regulate metabolic gene programs and may play a role in the shift of metabolic substrate preference seen in hypertrophy. Failing hearts have elevated expression of the nuclear hormone receptor COUP-TF. We report that COUP-TF significantly suppresses reporter transcription independent of DNA binding and specific interactions with GATA-4, Nkx2.5 or USF. In summary, CPT-1β transcriptional regulation integrates mitochondrial gene expression with two essential cardiac functions: contraction and metabolic substrate oxidation. ^

Regulation of tyrosine hydroxylase gene expression by mRNA stabilization

Tyrosine hydroxylase (TH) expression increases in adrenal chromaffin cells treated with the nicotinic agonist, dimethylphenylpiperazinium (DMPP; 1 μM). We are using this response as a model of the changes in TH level that occur during increased cholinergic neural activity. Here we report a 4-fold increase in TH mRNA half-life in DMPP-treated chromaffin cells that is apparent when using a pulse-chase analysis to measure TH mRNA half-life. No increase is apparent using actinomycin D to measure half-life, indicating a requirement for ongoing transcription. Characterization of protein binding to the TH 3′UTR using RNA electro-mobility shift assays show the presence of two complexes both of which are increased by DMPP-treatment. The faster migrating complex (FMC) increases 2.5-fold and the slower migrating complex (SMC) increases 1.5-fold. Separation of UV crosslinked RNA-protein complexes on SDS polyacrylamide gels shows FMC to contain a single protein whereas SMC contains two proteins. Northwesterns yielded similar results. Transfection studies reveal an increase in expression of the full-length TH transcript due to DMPP-treatment similar to that of endogenous TH mRNA. This finding suggests the increased expression is due primarily to mRNA stabilization. Transfection of luciferase reporter constructs containing regions of the TH 3′UTR reveal only the full-length 3′UTR influenced the expression level of reporter transcripts. ^

Gain-of-function of the p53R172H mutation in a mouse model of Li -Fraumeni syndrome

Missense mutations in the p53 tumor-suppressor gene are the most common alterations of p53 in somatic tumors and in patients with Li-Fraumeni syndrome. p53 missense mutations occur in the DNA binding region and disrupt the ability of p53 to activate transcription. In vitro studies have shown that some p53 missense mutants have a gain-of-function or dominant-negative activity. ^ The p53 175 Arg-to-His (p53 R175H) mutation in humans has been shown to have dominant-negative and gain-of-function properties in vitro. This mutation is observed in the germline of individuals with Li-Fraumeni syndrome. To accurately model Li-Fraumeni syndrome and to examine the mechanistic nature of a gain-of-function missense mutation on in vivo tumorigenesis, we generated and characterized a mouse with the corresponding mutation, p53 R172H. p53R172H homozygous and heterozygous mice developed similar tumor spectra and survival curves as p53 −/− and p53+/− mice, respectively. However, tumors in p53+/R172H mice metastasized to various organs with high frequency, suggesting a gain-of-function phenotype by p53R172H in vivo. Mouse embryonic fibroblasts (MEFs) from p53R172H mice also showed gain-of-function phenotypes in cell proliferation, DNA synthesis, and transformation potential, while cells from p53+/− and p53−/− mice did not. ^ To mechanistically characterize the gain-of-function phenotype of the p53R172H mutant, the role of p53 family members, p63 and p73, was analyzed. Disruption of p63 and p73 by siRNAs in p53 −/− MEFs increased transformation potential and reinitiated DNA synthesis to levels observed in p53R172H/R172H cells. Additionally, p63 and p73 were bound and functionally inactivated by p53R172H in metastatic p53 R172H tumor-derived cell lines, indicating a role for the p53 family members in the gain-of-function phenotype. This study provides in vivo evidence for the gain-of-function effect of p53 missense mutations and more accurately models the Li-Fraumeni syndrome. ^

Breast cancer -specific activation of topoisomerase IIalpha and its application in the tumor -targeting gene therapy

To ensure the success of systemic gene therapy, it is critical to enhance the tumor specificity and activity of the promoter. In the current study, we identified the breast cancer-specific activity of the topoisomerase IIα promoter. We further showed that cdk2 and cyclin A activate topoisomerase IIα promoter in a breast cancer-specific manner. An element containing an inverted CCAAT box (ICB) was shown to respond this signaling. When the ICB-harboring topoisomerase IIα minimal promoter was linked with an enhancer sequence from the cytomegalovirus immediate early gene promoter (CMV promoter), this composite promoter, CT90, exhibited activity comparable to or higher than the CMV promoter in breast cancer cells in vitro and in vivo, yet expresses much lower activity in normal cell lines and normal organs than the CMV promoter. A CT90-driven construct expressing BikDD, a potent pro-apoptotic gene, was shown to selectively kill breast cancer cells in vitro and to suppress mammary tumor development in an animal model of intravenously administrated, liposome-delivered gene therapy. Expression of BikDD was readily detectable in the tumors but not in the normal organs of CT90-BikDD-treated animals. Finally, we demonstrated that CT90-BikDD treatment potentially enhanced the sensitivity of breast cancer cells to chemotherapeutic agents, especially doxorubicin and taxol. The results indicate that liposomal CT90-BikDD is a novel and effective systemic breast cancer-targeting gene therapy, and its combination with chemotherapy may further improve the current adjuvant therapy for breast cancer. ^