High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

Application of cell line based genomic predictors to predict response to targeted therapies in breast cancer

Cancer cell lines can be treated with a drug and the molecular comparison of responders and non-responders may yield potential predictors that could be tested in the clinic. It is a bioinformatics challenge to apply the cell line-derived multivariable response predictors to patients who respond to therapy. Using the gene expression data from 23 breast cancer cell lines, I developed three predictors of dasatinib sensitivity by selecting differentially expressed genes and applying different classification algorithms. The performance of these predictors on independent cell lines with known dasatinib response was tested. The predictor based on weighted voting method has the best overall performance. It correctly predicted dasatinib sensitivity in 11 out of 12 (92%) breast and 17 out of 23 (74%) lung cancer cell lines. These predictors were then applied to the gene expression data from 133 breast cancer patients in an attempt to predict how the patients might respond to dasatinib therapy. Two predictors identified 13 patients in common to be dasatinib sensitive. Sixty two percent of these cases are triple negative (ER-negative, HER2-negative and PR-negative) and 76% are double negative. The result is consistent with the findings from other studies, which identified a target population for dasatinib treatment to be triple negative or basal breast cancer subtype. In conclusion, we think that the cell line-derived dasatinib classifiers can be applied to the human patients. ^

Integrating sequence information in microarray data analysis by free energy modeling

Microarray technology is a high-throughput method for genotyping and gene expression profiling. Limited sensitivity and specificity are one of the essential problems for this technology. Most of existing methods of microarray data analysis have an apparent limitation for they merely deal with the numerical part of microarray data and have made little use of gene sequence information. Because it's the gene sequences that precisely define the physical objects being measured by a microarray, it is natural to make the gene sequences an essential part of the data analysis. This dissertation focused on the development of free energy models to integrate sequence information in microarray data analysis. The models were used to characterize the mechanism of hybridization on microarrays and enhance sensitivity and specificity of microarray measurements. ^ Cross-hybridization is a major obstacle factor for the sensitivity and specificity of microarray measurements. In this dissertation, we evaluated the scope of cross-hybridization problem on short-oligo microarrays. The results showed that cross hybridization on arrays is mostly caused by oligo fragments with a run of 10 to 16 nucleotides complementary to the probes. Furthermore, a free-energy based model was proposed to quantify the amount of cross-hybridization signal on each probe. This model treats cross-hybridization as an integral effect of the interactions between a probe and various off-target oligo fragments. Using public spike-in datasets, the model showed high accuracy in predicting the cross-hybridization signals on those probes whose intended targets are absent in the sample. ^ Several prospective models were proposed to improve Positional Dependent Nearest-Neighbor (PDNN) model for better quantification of gene expression and cross-hybridization. ^ The problem addressed in this dissertation is fundamental to the microarray technology. We expect that this study will help us to understand the detailed mechanism that determines sensitivity and specificity on the microarrays. Consequently, this research will have a wide impact on how microarrays are designed and how the data are interpreted. ^

Expression cloning of a candidate prostate cancer oncogene

Genetic analysis, both karyotyping and comparative genomic hybridization, of prostate cancer cell lines and specimens have revealed multiple areas of concordant increases in DNA content. An increase of DNA in specific regions of the genome in cancer is often associated with the amplification of oncogenes. Based on these observations we have hypothesized that oncogenes are involved in the initiation or progression of prostate cancer. An expression cloning approach was utilized to identify candidate oncogenes in prostate cancer. ^ A full-length, unidirectional cDNA expression library was constructed from DU145 prostate cancer cells. The cDNA library was screened using CP12, a rat prostate epithelial cell line. In soft agarose assays, CP12 (parental or vector transfected) do not form colonies. However, upon the introduction of a number of known oncogenes CP12 becomes anchorage independent in soft agarose. Based on this in-vitro phenotypic shift, a DU145 cDNA library was stably transfected into CP12, and selected for anchorage independence. Two hundred fifty nine anchorage independent clones were isolated. Some colonies contained more than one insert, bringing the candidate oncogene pool to approximately 400. Seven inserts were sequenced at random. Using the sequences obtained, GenBank was screened, and matches were found with p53, PARG1, a mitochondrial ATPase, RNF6, and three unknown genes that mapped to Unigene clusters. As the pool of cDNA inserts appeared promising, overexpressed genes were further selected. From 259 clones, 17 clones were overexpressed more than 6-fold in DU145 compared to Normal Prostate. From the 17 clones, 12 cDNA inserts were strongly expressed in DU145 and were isolated for sequencing. ^ Two of the sequences, 1G6 and 3E9, were identical. Expression of 1G6/2G9/3E9 was tested by RT-PCR. 1G6/2G9/3E9 was not expressed in normal prostate, but was expressed in all prostate cancer cell lines tested as well as six prostate cancer samples. When retransfected into CP12, 1G6/2G9/3E9 induced the formation of foci and anchorage independent colonies. Thus, functional and expression data suggest that 1G6/2G9/3E9 may be a prostate cancer oncogene. ^

Physical and functional mapping of a novel tumor suppressor locus for prostate cancer within chromosome 10p12.31-q11

Prostate cancer remains the second leading cause of male cancer deaths in the United States, yet the molecular mechanisms underlying this disease remain largely unknown. Cytogenetic and molecular analyses of prostate tumors suggest a consistent association with the loss of chromosome 10. Previously, we have defined a novel tumor suppressor locus PAC-1 within chromosome 10pter-q11. Introduction of the short arm of chromosome 10 into a prostatic adenocarcinoma cell line PC-3H resulted in dramatic tumor suppression and restoration of a programmed cell death pathway. Using a combined approach of comparative genomic hybridization and microsatellite analysis of PC-3H, I have identified a region of hemizygosity within 10p12-p15. This region has been shown to be involved in frequent loss of heterozygosity in gliomas and melanoma. To functionally dissect the region within chromosome 10p containing PAC-1, we developed a strategy of serial microcell fusion, a technique that allows the transfer of defined fragments of chromosome 10p into PC-3H. Serial microcell fusion was used to transfer defined 10p fragments into a mouse A9 fibrosarcoma cell line. Once characterized by FISH and microsatellite analyses, the 10p fragments were subsequently transferred into PC-3H to generate a panel of microcell hybrid clones containing overlapping deletions of chromosome 10p. In vivo and microsatellite analyses of these PC hybrids identified a small chromosome 10p fragment (an estimated 31 Mb in size inclusive of the centromere) that when transferred into the PC-3H background, resulted in significant tumor suppression and limited a region of functional tumor suppressor activity to chromosome 10p12.31-q11. This region coincides with a region of LOH demonstrated in prostate cancer. These studies demonstrate the utility of this approach as a powerful tool to limit regions of functional tumor suppressor activity. Furthermore, these data used in conjunction with data generated by the Human Genome Project lent a focused approach to identify candidate tumor suppressor genes involved in prostate cancer. ^

PREDICTION OF DNA METHYLATION BASED ON GENOMIC ARCHITECTURE AND APPLICATIONS OF POSITIONAL WEIGHT MATRICES

Gene silencing due to epigenetic mechanisms shows evidence of significant contributions to cancer development. We hypothesis that the genetic architecture based on retrotransposon elements surrounding the transcription start site, plays an important role in the suppression and promotion of DNA methylation. In our investigation we found a high rate of SINE and LINEs retrotransposon elements near the transcription start site of unmethylated genes when compared to methylated genes. The presence of these elements were positively associated with promoter methylation, contrary to logical expectations, due to the malicious effects of retrotransposon elements which insert themselves randomly into the genome causing possible loss of gene function. In our genome wide analysis of human genes, results suggested that 22% of the genes in cancer were predicted to be methylation-prone; in cancer these genes are generally down-regulated and function in the development process. In summary, our investigation validated our hypothesis and showed that these widespread genomic elements in cancer are highly associated with promoter DNA methylation and may further participate in influencing epigenetic regulation.

From protein microarray to biology: The discovery of epigenetic regulatory molecules

Chromatin, composed of repeating nucleosome units, is the genetic polymer of life. To aid in DNA compaction and organized storage, the double helix wraps around a core complex of histone proteins to form the nucleosome, and is therefore no longer freely accessible to cellular proteins for the processes of transcription, replication and DNA repair. Over the course of evolution, DNA-based applications have developed routes to access DNA bound up in chromatin, and further, have actually utilized the chromatin structure to create another level of complexity and information storage. The histone molecules that DNA surrounds have free-floating tails that extend out of the nucleosome. These tails are post-translationally modified to create docking sites for the proteins involved in transcription, replication and repair, thus providing one prominent way that specific genomic sequences are accessed and manipulated. Adding another degree of information storage, histone tail-modifications paint the genome in precise manners to influence a state of transcriptional activity or repression, to generate euchromatin, containing gene-dense regions, or heterochromatin, containing repeat sequences and low-density gene regions. The work presented here is the study of histone tail modifications, how they are written and how they are read, divided into two projects. Both begin with protein microarray experiments where we discover the protein domains that can bind modified histone tails, and how multiple tail modifications can influence this binding. Project one then looks deeper into the enzymes that lay down the tail modifications. Specifically, we studied histone-tail arginine methylation by PRMT6. We found that methylation of a specific histone residue by PRMT6, arginine 2 of H3, can antagonize the binding of protein domains to the H3 tail and therefore affect transcription of genes regulated by the H3-tail binding proteins. Project two focuses on a protein we identified to bind modified histone tails, PHF20, and was an endeavor to discover the biological role of this protein. Thus, in total, we are looking at a complete process: (1) histone tail modification by an enzyme (here, PRMT6), (2) how this and other modifications are bound by conserved protein domains, and (3) by using PHF20 as an example, the functional outcome of binding through investigating the biological role of a chromatin reader. ^

MOLECULAR AND GENOMIC BASED INSIGHTS INTO THE EVOLUTION OF ENTEROCOCCUS FAECIUM FROM COMMENSAL TO HOSPITAL-ADAPTED PATHOGEN

The basis for the recent transition of Enterococcus faecium from a primarily commensal organism to one of the leading causes of hospital-acquired infections in the United States is not yet understood. To address this, the first part of my project assessed isolates from early outbreaks in the USA and South America using sequence analysis, colony hybridizations, and minimal inhibitory concentrations (MICs) which showed clinical isolates possess virulence and antibiotic resistance determinants that are less abundant or lacking in community isolates. I also revealed that the level of ampicillin resistance increased over time in clinical strains. By sequencing the pbp5 gene, I demonstrated an ~5% difference in the pbp5 gene between strains with MICs <4ug/ml and those with MICs >4µg/ml, but no specific sequence changes correlated with increases in MICs within the latter group. A 3-10% nucleotide difference was also seen in three other genes analyzed, which suggested the existence of two distinct subpopulations of E. faecium. This led to the second part of my project analyzing concatenated core gene sequences, SNPs, the 16S rRNA, and phylogenetics of 21 E. faecium genomes confirming two distinct clades; a community-associated (CA) clade and hospital-associated (HA) clade. Molecular clock calculations indicate that these two clades likely diverged ~ 300,000 to > 1 million years ago, long before the modern antibiotic era. Genomic analysis also showed that, in addition to core genomic differences, HA E. faecium harbor specific accessory genetic elements that may confer selection advantages over CA E. faecium. The third part of my project discovered 6 E. faecium genes with the newly identified “WxL” domain. My analyses, using RT-PCR, western blots, patient sera, whole-cell ELISA, and immunogold electron microscopy, indicated that E. faecium WxL genes exist in operons, encode bacterial cell surface localized proteins, that WxL proteins are antigenic in humans, and are more exposed on the surface of clinical isolates versus community isolates (even though they are ubiquitous in both clades). ELISAs and BIAcore analyses also showed that proteins encoded by these operons bind several different host extracellular matrix proteins, as well as to each other, suggesting a novel cell-surface complex. In summary, my studies provide new insights into the evolution of E. faecium by showing that there are two distantly related clades; one being more successful in the hospital setting. My studies also identified operons encoding WxL proteins whose characteristics could also contribute to colonization and virulence within this species.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

PCR-based assay using occult blood detection cards for detection of diarrheagenic Escherichia coli in specimens from U.S. travelers to Mexico with acute diarrhea.

Large field studies of travelers' diarrhea for multiple destinations are limited by the need to perform stool cultures on site in a timely manner. A method for the collection, transport, and storage of fecal specimens that does not require immediate processing and refrigeration and that is stable for months would be advantageous. This study was designed to determine if enterotoxigenic Escherichia coli (ETEC) and enteroaggregative E. coli (EAEC) DNA could be identified from cards that were processed for the evaluation of fecal occult blood. U.S. students traveling to Mexico during 2005 to 2007 were monitored for the occurrence of diarrheal illness. When ill, students provided a stool specimen for culture and occult blood by the standard methods. Cards then were stored at room temperature prior to DNA extraction. Fecal PCR was performed to identify ETEC and EAEC in DNA extracted from stools and from occult blood cards. Significantly more EAEC cases were identified by PCR that was performed on DNA that was extracted from cards (49%) or from frozen feces (40%) than from culture methods that used HEp-2 adherence assays (13%) (P < 0.001). Similarly, more ETEC cases were detected from card DNA (38%) than from fecal DNA (30%) or by culture that was followed by hybridization (10%) (P < 0.001). The sensitivity and specificity of the card test were 75 and 62%, respectively, compared to those for EAEC by culture and were 50 and 63%, respectively, compared to those for ETEC. DNA extracted from fecal cards that was used for the detection of occult blood is of use in identifying diarrheagenic E. coli.

Aboral ectoderm-specific expression and regulation of the LpS1 genes of {\it Lytechinus pictus\/}

The Spec genes serve as molecular markers for examining the ontogeny of the aboral ectoderm lineage of sea urchin embryos. These genes are activated at late-cleavage stage only in cells contributing to the aboral ectoderm of Strongylocentrotus purpuratus and encode 14,000-17,000 Da calcium-binding proteins. A comparative analysis was undertaken to better understand the mechanisms underlying the activation and function of the Spec genes by investigating Spec homologues from Lytechinus pictus, a distantly related sea urchin. Spec antibodies cross-reacted with 34,000 Da proteins in L. pictus embryos that displayed a similar ontogenetic pattern to that of Spec proteins. One cDNA clone, LpS1, was isolated by hybridization to a synthetic oligonucleotide corresponding to a calcium-binding domain or EF-hand. The LpS1 mRNA has developmental properties similar to those of the Spec mRNAs. LpS1 encodes a 34,000 Da protein containing eight EF-hand domains, which share structural homology with the Spec EF-hands; however, little else in the protein sequence is conserved, implying that calcium-binding is important for Spec protein function. Genomic DNA blot analysis showed two LpS1 genes, LpS1$\alpha$ and LpS1$\beta$, in L. pictus. Partial gene structures for both LpS1$\alpha$ and $\beta$ were constructed based on genomic clones isolated from an L. pictus genomic library. These revealed internal duplications of the LpS1 genes that accounted for the eight EF-hand domains in the LpS1 proteins. Sequencing analysis showed there was little in common among the 5$\sp\prime$-flanking regions of the LpS1 and Spec genes except for the presence of a binding site for the transcription factor USF.^ A sea urchin gene-transfer expression system showed that 762 base pairs (bp) of 5$\sp\prime$-flanking DNA from the LpS1$\beta$ gene were sufficient for correct temporal and spatial expression of reporter genes in sea urchin embryos. Deletions at the 5$\sp\prime$ end to 511, 368, or 108bp resulted in a 3-4 fold decrease in chloramphenicol acetyltransferase (CAT) activity and disrupted the restricted activation of the lac Z gene in aboral ectoderm cells.^ A full-length Spec1 protein and a truncated LpS1 protein were induced and partially purified from an in vitro expression system. (Abstract shortened with permission of author.) ^

The physical map of {\it Escherichia coli\/} K-12 strain MG1655: Construction, comparison with other strains and implications for mechanisms controlling oligonucleotide frequency

Complete NotI, SfiI, XbaI and BlnI cleavage maps of Escherichia coli K-12 strain MG1655 were constructed. Techniques used included: CHEF pulsed field gel electrophoresis; transposon mutagenesis; fragment hybridization to the ordered $\lambda$ library of Kohara et al.; fragment and cosmid hybridization to Southern blots; correlation of fragments and cleavage sites with EcoMap, a sequence-modified version of the genomic restriction map of Kohara et al.; and correlation of cleavage sites with DNA sequence databases. In all, 105 restriction sites were mapped and correlated with the EcoMap coordinate system.^ NotI, SfiI, XbaI and BlnI restriction patterns of five commonly used E. coli K-12 strains were compared to those of MG1655. The variability between strains, some of which are separated by numerous steps of mutagenic treatment, is readily detectable by pulsed-field gel electrophoresis. A model is presented to account for the difference between the strains on the basis of simple insertions, deletions, and in one case an inversion. Insertions and deletions ranged in size from 1 kb to 86 kb. Several of the larger features have previously been characterized and some of the smaller rearrangements can potentially account for previously reported genetic features of these strains.^ Some aspects of the frequency and distribution of NotI, SfiI, XbaI and BlnI cleavage sites were analyzed using a method based on Markov chain theory. Overlaps of Dam and Dcm methylase sites with XbaI and SfiI cleavage sites were examined. The one XbaI-Dam overlap in the database is in accord with the expected frequency of this overlap. The occurrence of certain types of SfiI-Dcm overlaps are overrepresented. Of the four subtypes of SfiI-Dcm overlap, only one has a partial inhibitory effect on the activity of SfiI. Recognition sites for all four enzymes are rarer than expected based on oligonucleotide frequency data, with this effect being much stronger for XbaI and BlnI than for NotI and SfiI. The latter two enzyme sites are rare mainly due to apparent negative selection against GGCC (both) and CGGCCG (NotI). The former two enzyme sites are rare mainly due to effects of the VSP repair system on certain di-tri- and tetranucleotides, most notably CTAG. Models are proposed to explain several of the anomalies of oligonucleotide distribution in E. coli, and the biological significance of the systems that produce these anomalies is discussed. ^

THE TRANSLATION PRODUCTS OF MOLONEY MURINE SARCOMA VIRUS 124 GENOMIC RNA

Murine sarcoma viruses constitute a class of replication-defective retroviruses. Cellular transformation may be induced by these viruses in vitro; whereas, fibrosarcomas may result in animals infected with them in vivo (Tooze, 1973; Bishop, 1978). Hybridization studies suggest that murine sarcoma viruses arose by recombination between nondefective murine leukemia virus sequences and certain cellular sequences present in uninfected mouse cells (Hu et al., 1977). A specific gene product, however, has not been implicated in murine sarcoma virus transformation.^ One line of murine sarcoma virus-producing cells, Mo-MuSV-clone 124, (Ball et al., 1973), was studied biochemically because it mainly produces the sarcoma virus as a pseudotype packaged with helper murine leukemia virus proteins. The sarcoma viral RNA was translated in a sophisticated cell-free protein synthesizing system (Murphy and Arlinghaus, 1978). The translation products were analyzed by a number of techniques, including electrophoresis in denaturing gels of SDS polyacrylamide, immunoprecipitation, and peptide mapping. The major products of the total RNA purified from the virus preparation were shown to have molecular weights of about 63,000 (P63('gag)), 42,000 (P42), 40,000 (P40), 38,000 (P38), and 23,000 (P23). The size class of mRNA coding for each of the cell-free products was estimated using a poly(A) selection technique and sucrose gradient fractionation. These analyses were used to localize the coding information related to each of the in vitro synthesized cell-free products within the sarcoma virus genome.^ The major findings of these studies were: (1) the 5' half of the sarcoma viral RNA codes for the 63,000 dalton polypeptide and 42,000 - 38,000 dalton polypeptides derived from the "gag" gene; and (2) the 3' half of the sarcoma viral RNA codes for a 38,000 dalton polypeptide and possibly derived from the cellular acquired sequences. ^

A PCR based assay for the detection of enteric pathogens from Hemoccult(RTM) cards

Background. Large field studies in travelers' diarrhea (TD) in multiple destinations are limited by the need to perform stool cultures on site in a timely manner. A method for the collection, transport and storage of fecal specimens that does not require immediate processing, refrigeration and is stable for months would be advantageous. ^ Objectives. Determine if enteric pathogen bacterial DNA can be identified in cards routinely used for evaluation of fecal occult blood. ^ Methods. U.S. students traveling to Mexico in 2005-07 were followed for occurrence of diarrheal illness. When ill, students provided a stool specimen for culture and occult blood by the standard method. Cards were then stored at room temperature prior to DNA extraction. A multiplex fecal PCR was performed to identify enterotoxigenic Escherichia coli and enteroaggregative E. coli (EAEC) in DNA extracted from stools and occult blood cards. ^ Results. Significantly more EAEC cases were identified by PCR done in DNA extracted from cards (49%) or from frozen feces (40%) than by culture followed by HEp-2 adherence assays (13%). Similarly more ETEC cases were detected in card DNA (38%) than fecal DNA (30%) or culture followed by hybridization (10%). Sensitivity and specificity of the card test was 75% and 62%, respectively, and 50% and 63%, respectively, when compared to EAEC and ETEC culture, respectively, and 53% and 51%, respectively compared to EAEC multiplex fecal PCR and 56% and 70%, respectively, compared to ETEC multiplex fecal PCR. ^ Conclusions. DNA extracted from fecal cards used for detection of occult blood is of use in detecting enteric pathogens. ^

Identifying and Targeting Molecular Deregulations in Uterine Leiomyosarcoma: A Role for mTOR Inhibition-Based Combination Therapies

Uterine leiomyosarcoma (ULMS) is an aggressive malignancy characterized by marked chemoresistance, frequent relapses, and poor outcome. Despite efforts to improve survival over the past several decades, only minimal advances have been made. Hence, there is an urgent and unmet need for better understanding of the molecular deregulations that underlay ULMS and development of more effective therapeutic strategies. This work identified several common deregulations in a large (n=208) tissue microarray of ULMS compared to GI smooth muscle, myometrium, and leiomyoma controls. Our results suggest that significant loss of smooth muscle and gynecological differentiation markers is common in ULMS, a finding that could help render improved ULMS diagnosis, especially for advanced disease. Similarly to reports in other malignancies, we found that several cancer-related proteins were differentially expressed; these could be useful together as biomarkers for ULMS. Notably, we identified significant upregulation and overexpression of the mTOR pathway in ULMS, examined the possible contribution of tyrosine kinase receptor deregulation promoting mTOR activation, and unraveled a role for pS6RP and p4EBP1 as molecular disease prognosticators. The significance of mTOR activation in ULMS and its potential as a therapeutic target were further investigated. Rapamycin abrogated ULMS cell growth and cell cycle progression in vitro but induced only sight growth delay in vivo. Given that effective mTOR therapies likely require combination mTOR blockade with inhibition of other targets, coupled with recent observations suggesting that Aurora A kinase (Aurk A) deregulations commonly occur in ULMS, the preclinical impact of dually targeting both pathways was evaluated. Combined therapy with rapamycin (an mTORC1 inhibitor) and MLN8237 (an investigational Aurk A inhibitor) profoundly and synergistically abrogated ULMS growth in vitro. Interestingly, the superior effects were noted only when MLN8237 was pre-administered. This novel therapeutic combination and scheduling regimen resulted in marked tumor growth inhibition in vivo. Together, these data support further exploration of dual mTOR and Aurk A blockade for the treatment of human ULMS.