A coalescent analysis for modeling the mutation process in colorectal cancer

Colorectal cancer is the forth most common diagnosed cancer in the United States. Every year about a hundred forty-seven thousand people will be diagnosed with colorectal cancer and fifty-six thousand people lose their lives due to this disease. Most of the hereditary nonpolyposis colorectal cancer (HNPCC) and 12% of the sporadic colorectal cancer show microsatellite instability. Colorectal cancer is a multistep progressive disease. It starts from a mutation in a normal colorectal cell and grows into a clone of cells that further accumulates mutations and finally develops into a malignant tumor. In terms of molecular evolution, the process of colorectal tumor progression represents the acquisition of sequential mutations. ^ Clinical studies use biomarkers such as microsatellite or single nucleotide polymorphisms (SNPs) to study mutation frequencies in colorectal cancer. Microsatellite data obtained from single genome equivalent PCR or small pool PCR can be used to infer tumor progression. Since tumor progression is similar to population evolution, we used an approach known as coalescent, which is well established in population genetics, to analyze this type of data. Coalescent theory has been known to infer the sample's evolutionary path through the analysis of microsatellite data. ^ The simulation results indicate that the constant population size pattern and the rapid tumor growth pattern have different genetic polymorphic patterns. The simulation results were compared with experimental data collected from HNPCC patients. The preliminary result shows the mutation rate in 6 HNPCC patients range from 0.001 to 0.01. The patients' polymorphic patterns are similar to the constant population size pattern which implies the tumor progression is through multilineage persistence instead of clonal sequential evolution. The results should be further verified using a larger dataset. ^

Studies of pseudoachondroplasia chondrocytes and knockdown of the cartilage oligomeric matrix protein

Cartilage oligomeric matrix protein (COMP) is a large, homopentameric, extracellular matrix glycoprotein. Mutations in COMP cause two skeletal dysplasias: pseudoachondroplasia (PSACH) and multiple epiphyseal dysplasia (EMD1). These dwarfing conditions are caused by retention of misfolded mutant COMP with type IX collagen and matrilin-3 (MATN3) in the rough endoplasmic reticulum (rER) of the chondrocyte. These proteins form a matrix in the rER that continues to expand until it fills the entire cell, eventually causing cell death. Interestingly, loss of COMP in COMP null mice does not affect normal bone development or growth, suggesting that elimination of COMP (wildtype and mutant) expression may prevent PSACH. The hypothesis of these studies was that a hammerhead ribozyme could eliminate or knockdown COMP mRNA expression in PSACH chondrocytes . To test this hypothesis, a human chondrocyte model system that recapitulates the PSACH chondrocyte phenotype was developed by over-expressing mutant (mt-) COMP in normal chondrocytes using a recombinant adenovirus. Chondrocytes over-expressing mt-COMP developed giant rER cisternae containing COMP, type IX collagen and MATN3. Deconvolution microscopy and computer modeling showed that these proteins formed an ordered matrix surrounding a type II pro-collagen core. Additionally, the results show that a hammerhead ribozyme, ribozyme 56 (Ribo56) reduced over-expressed mt-COMP in COS cells and endogenous COMP in normal chondrocytes and mt-COMP in three PSACH chondrocytes cell line (with different mutations) by 40-70%. Altogether, these studies show that the PSACH cellular phenotype can be created in vitro and that the mt-COMP protein burden can be reduced by the presence of a COMP-specific ribozyme. Future studies will focus on designing ribozymes or short interfering RNA (siRNA) technologies that will result in better knockdown of COMP expression as well as the temporal constraints imposed by the PSACH phenotype. ^

Impact of COP9 signalosome subunit 6 (CSN6) on MDM2-p53 axis in DNA damage-mediated apoptosis and tumorigenesis

Mammalian COP9 signalosome, which connects signaling with the ubiquitin-mediated proteasome degradation pathway, is implicated in cell cycle regulation and DNA damage response. However, whether COP9 is dysregulated in cancers has not been well established. Here, we showed that COP9 subunit 6 (CSN6) was upregulated in malignant breast and thyroid tumors and positively correlated with MDM2 expression. Investigation of the underlying mechanism suggested that CSN6 stabilized MDM2, thereby accelerating the degradation of p53. We generated mice carrying a targeted disruption of the Csn6 gene, and found that the mice with both alleles disrupted (Csn6-/- ) died in early embryogenesis (E7.5). Csn6+/- mice were sensitized to undergo γ-radiation-induced p53-dependent apoptosis in both thymus and developing central nervous system. Consequently. Csn6 +/- mice were more susceptible to the lethal effects of high-dose γ-radiation than wild-type mice. Notably, Csn6+/- mice were less susceptible to γ-radiation-induced tumorigenesis and had better long-term survival after low-dose γ-radiation exposure compared with wild-type animals, indicating that loss of CSN6 enhanced p53-mediated tumor suppression in vivo. In summary, the regulation of MDM2-p53 signaling by CSN6 plays a significant role in DNA damage-mediated apoptosis and tumorigenesis, which suggests that CSN6 may potentially be a valuable diagnostic marker for cancers with a dysregulated MDM2-p53 axis. ^

EGFR and Notch signaling pathways are regulated by distinct isoforms of Drosophila cbl

Activation of cell surface receptors transduces extracellular signals into cellular responses such as proliferation, differentiation and survival. However, the appropriate spatial and temporal down-regulation of signaling receptors is essential for normal development and homeostasis. The Cbl family of E3-ubiquitin ligases plays a major role for the ligand-dependent inactivation of growth factor receptors through ubiquitin-mediated endocytosis and lysosomal degradation. Here, we report the D-cbl mutant phenotypes in the Drosophila eye. D-cbl mutants display overgrowth, inhibition of apoptosis, differentiation defects and increased ommatidial spacing. Many of these phenotypes are caused by lack of down-regulation of the Drosophila EGFR signaling. However, not all D-cbl phenotypes can be explained by inappropriate EGFR activity. We found that D-Cbl also negatively regulates Notch activity during eye and wing development. D-cbl produces two isoforms by alternative splicing. Strikingly, the long isoform, D-CblL, preferentially regulates the EGFR, whereas the short isoform, D-CblS, preferentially regulates Notch. Taken together, these data suggest that D-Cbl controls at least two signaling pathways, EGFR and Notch, through production of two alternatively spliced isoforms during development in Drosophila.^

Identification of germ cell tumor modifier genes from the 129.MOLF-Chr19 consomic mouse strain

Naturally occurring genetic variants confer susceptibility to disease in the human population, including in testicular germ cell tumor development. Disease susceptibility loci for testicular germ cell tumors have been identified by genetic mapping in humans and mice. However, the identity of many of the susceptibility genes remains unclear. My study utilized a chromosome substitution strain, the 129.MOLF-Chr 19 (or M19 strain), to identify candidate testicular germ cell tumor susceptibility genes. Males of this strain have a high incidence of germ cell tumors in the testes. By forward genetic approaches, five susceptibility loci were fine-mapped and the genetic interactions were dissected. In addition, I identified three protein-coding genes and one micro-RNA as testicular tumor susceptibility genes by genomic screening. Using reverse genetic approaches, I verified one of the candidates, Splicing factor 1, as a modifier of testicular tumor. Deficiency of SF1 significantly reduces the incidence of testicular tumors in mice. This study highlights the advantage of the 129.MOLF-Chr 19 consomic strain in disease gene identification and validation. It also sets the stage to elucidate the molecular mechanisms of tumorigenesis in the testis. ^

The INO80 chromatin remodeling complex is involved in the nucleotide excision repair pathway

The genomic DNA of eukaryotic cells is well organized into chromatin structures. However, these repressed structures present barriers that block the access of regulatory factors to the genome during various nuclear events. To overcome the obstacle, two major cellular processes, post-modification of histone tails and ATP-dependent chromatin remodeling, are involved in reconfiguring chromatin structure and creating accessible DNA. Despite the current research progress, much remains to be explored concerning the relationship between chromatin remodeling and DNA repair. Recently, one member of the ATP-dependent chromatin remodeling complexes, INO80, has been found to play a crucial role in DNA damage repair. However, the functions of this complex in higher eukaryotes have yet to be determined. The goal of my study is to generate a human somatic INO80 conditional knockout model and investigate the functions of Ino80 in damage repair.^ By homologous targeting of the INO80 locus in human HCT116 colon epithelial cells, I established a human somatic INO80 conditional knockout model. I have demonstrated that the conditional INO80 cells exhibited a sufficiently viable period when the INO80 protein is removed. Moreover, I found that loss of INO80 resulted in deficient UV lesion repair in response to UV while the protein levels of the NER factors such as XPC, XPA, XPD were not affected. And in vitro repair synthesis assay showed that the NER incision and repair synthesis activities were intact in the absence of INO80. Examination on the damage recognition factor XPC showed its recruitment to damage sites was impaired in the INO80 mutant cells. Loss of INO80 also led to reduced enrichment of XPA at the site of UV lesions. Despite the reduced recruitment of XPC and XPA observed in INO80 mutants, no direct interaction was detected. Meanwhile, direct interaction between INO80 and DDB1, the initial UV lesion detector, was detected by coimmunoprecipitation. UV-induced chromosome relaxation was reduced in cells devoid of INO80. These results demonstrate the INO80 complex may participates in the NER by interacting with DDB1 and having a critical role of in creating DNA accessibility for the nucleotide excision pathway. ^

The critical role and regulation of transcription factor FoxM1 in gastric cancer development and progression

The mammalian Forkhead Box (Fox) transcription factor (FoxM1) is implicated in tumorgenesis. However, the role and regulation of FoxM1 in gastric cancer remain unknown.^ I examined FoxM1 expression in 86 cases of primary gastric cancer and 57 normal gastric tissue specimens. I found weak expression of FoxM1 protein in normal gastric mucosa, whereas I observed strong staining for FoxM1 in tumor-cell nuclei in various gastric tumors and lymph node metastases. The aberrant FoxM1 expression is associated with VEGF expression and increased angiogenesis in human gastric cancer. A Cox proportional hazards model revealed that FoxM1 expression was an independent prognostic factor in multivariate analysis. Furthermore, overexpression of FoxM1 by gene transfer significantly promoted the growth and metastasis of gastric cancer cells in orthotopic mouse models, whereas knockdown of FoxM1 expression by small interfering RNA did the opposite. Next, I observed that alteration of tumor growth and metastasis by elevated FoxM1 expression was directly correlated with alteration of VEGF expression and angiogenesis. In addition, promotion of gastric tumorigenesis by FoxM1 directly and significantly correlated with transactivation of vascular endothelial growth factor (VEGF) expression and elevation of angiogenesis. ^ To further investigate the underlying mechanisms that result in FoxM1 overexpression in gastric cancer, I investigated FoxM1 and Krüppel-like factor 4 (KLF4) expressions in primary gastric cancer and normal gastric tissue specimens. Concomitance of increased expression of FoxM1 protein and decreased expression of KLF4 protein was evident in human gastric cancer. Enforced KLF4 expression suppressed FoxM1 protein expression. Moreover, a region within the proximal FoxM1 promoter was identified to have KLF4-binding sites. Finally, I found an increased FoxM1 expression in gastric mucosa of villin-Cre -directed tissue specific Klf4-null mice.^ In summary, I offered both clinical and mechanistic evidence that dysregulated expression of FoxM1 play an important role in gastric cancer development and progression, while KLF4 mediates negative regulation of FoxM1 expression and its loss significantly contributes to FoxM1 dysregulation. ^

Identification and characterization of virulence determinants in the Lyme disease spirochete Borrelia burgdorferi

Borrelia burgdorferi is the etiological agent of Lyme disease, the most common tick-borne disease in the United States. Although the most frequently reported symptom is arthritis, patients can also experience severe cardiac, neurologic, and dermatologic abnormalities. The identification of virulence determinants in infectious B. burgdorferi strains has been limited by their slow growth rate, poor transformability, and general lack of genetic tools. The present study demonstrates the use of transposon mutagenesis for the identification of infectivity-related factors in infectious B. burgdorferi, examines the potential role for chemotaxis in mammalian infection, and describes the development of a novel method for the analysis of recombination events at the Ids antigenic variation locus. A pool of Himar1 mutants was isolated using an infectious B. burgdorferi clone and the transposon vector pMarGent. Clones exhibiting reduced infectivity in mice possessed insertions in virulence determinants putatively involved in host survival and dissemination. These results demonstrated the feasibility of extensive transposon mutagenesis studies for the identification of additional infectivity-related factors. mcp-5 mutants were chosen for further study to determine the role of chemotaxis during infection. Animal studies indicated that mcp-5 mutants exhibited a reduced infectivity potential, and suggested a role for mcp-5 during the early stages of infection. An in vitro phenotype for an mcp-5 mutant was not detected. Genetic complementation of an mcp-5 mutant resulted in restoration of Mcp-5 expression in the complemented clone, as demonstrated by western blotting, but the organisms were not infectious in mice. We believe this result is a consequence of differences in expression between genes located on the linear chromosome and genes present on the circular plasmid used for trans-complementation. Overall, this work implicates mcp-5 as an important determinant of mammalian infectivity. Finally, the development of a computer-assisted method for the analysis of recombination events occurring at the B. burgdorferi vls antigenic variation locus has proven highly valuable for the detailed examination of vls gene conversion. The studies described here provide evidence for the importance of chemotaxis during infection in mice and demonstrate advances in both genetic and computational approaches for the further characterization of the Lyme disease spirochete. ^

Transcriptional regulation of lung tumor suppressor GPRC5A in malignant and normal cells

Chronic exposure of the airways to cigarette smoke induces inflammatory response and genomic instability that play important roles in lung cancer development. Nuclear factor kappa B (NF-κB), the major intracellular mediator of inflammatory signals, is frequently activated in preneoplastic and malignant lung lesions. ^ Previously, we had shown that a lung tumor suppressor GPRC5A is frequently repressed in human non-small cell lung cancers (NSCLC) cells and lung tumor specimens. Recently, other groups have shown that human GPRC5A transcript levels are higher in bronchial samples of former than of current smokers. These results suggested that smoking represses GPRC5A expression and thus promotes the occurrence of lung cancer. We hypothesized that cigarette smoking or associated inflammatory response repressed GPRC5A expression through NF-κB signaling. ^ To determine the effect of inflammation, we examined GPRC5A protein expression in several lung cell lines following by TNF-α treatment. TNF-α significantly suppressed GPRC5A expression in normal small airway epithelial cells (SAEC) as well as in Calu-1 cells. Real-time PCR analysis indicated that TNF-α inhibits GPRC5A expression at the transcriptional level. NF-κB, the major downstream effectors of TNF-α signaling, mediates TNF-α-induced repression of GPRC5A because over-expression of NF-κB suppressed GPRC5A. To determine the region in the GPRC5A promoter through which NF-κB acts, we examined the ability of TNF-α to inhibit a series of reporter constructs with different deletions of GPRC5A promoter. The luciferase assay showed that the potential NF-κB binding sites containing region are irresponsible for TNF-α-induced suppression. Further analysis using constructs with different deletions in p65 revealed that NF-κB-mediated repression of GPRC5A is transcription-independent. Co-immunoprecipitation assays revealed that NF-κB could form a complex with RAR/RXR heterodimer. Moreover, the inhibitory effect of NF-κB has been found to be proportional to NF-κB/RAR ratio in luciferase assay. Finally, Chromatin IP demonstrated that NF-κB/p65 bound to GPRC5A promoter as well as RAR/RXR and suppressed transcription. Taken together, we propose that inflammation-induced NF-κB activation disrupts the RA signaling and suppresses GPRC5A expression and thus contributes to the oncogenesis of lung cancer. Our studies shed new light on the pathogenesis of lung cancer and potentially provide novel interventions for preventing and treating this disease. ^

Role of ovarian cancer antigen CA125/MUC16 in development and cancer

Cancer antigen 125 (CA125) is a tumor antigen that is routinely used to monitor the disease progress and the outcome of treatment in ovarian cancer patients. Elevated serum levels of CA125 are detected in over 80% of epithelial ovarian cancer patients. CA125 is a high molecular weight (>1M Dalton) mucin-type glycoprotein encoded by the MUC16 gene on human chromosome 19. Although MUC16 has served as the best serum marker for monitoring growth of ovarian cancer, roles for MUC16 in normal physiology and ovarian cancer are largely unknown. To understand the biological functions of MUC16, I characterized a mouse Muc16 homolog on chromosome 9 by means of expression pattern profiling, phenotype analysis of Muc16 knockout mice, and in vitro and in vivo studies of Muc16 null transformed ovarian surface epithelial (OSE) cells. ^ The mouse Muc16 homolog shares a conserved genomic structure with human MUC16. In addition to being expressed in mouse ovarian cancer, mouse Muc16 mRNA and protein were expressed in the mesothelia covering the heart, lung, ovary, oviduct, spleen, testis, and uterus. The conserved genomic structure and expression pattern of mouse Muc16 to human MUC16 suggests that mouse Muc16 is the ortholog of human MUC16. To understand the biological functions of Muc16, I generated Muc16 knockout mice. Muc16 knockout mice were viable, fertile and normal by one year of age. However, between 18 and 24 months of age, Muc16 knockout mice developed various tissue abnormalities such as ovarian cysts and tumors of the liver and other peritoneal organs. To determine the role of MUC16 in ovarian cancer progression, I established Muc16 null transformed ovarian surface epithelial (OSE) cell lines, following the same method to develop mouse model of epithelial ovarian cancer (Orsulic et al., 2002). Loss of Muc16 did not affect cell morphology, cell proliferation rate, or tumorigenic potential. However, Muc16-null OSE cells showed decreased attachment to extracellular matrix proteins as well as to primary mouse peritoneal mesothelial cells. Peritoneal mesothelia are the most frequent implantation sites of ovarian cancer. Furthermore, a pilot transplantation assay suggests that Muc16 null transformed OSE cells formed less disseminated tumors in the peritoneal cavity compared to wild-type OSE cells. ^ In conclusion, these results demonstrate that MUC16 is not required for normal mouse development or reproduction, but plays important roles in tissue homeostasis, ovarian cancer cell adhesion and dissemination. This study provides the first in vivo evidence of the roles of MUC16 in development, as well as ovarian cancer progression and dissemination. These studies offer valuable insights into possible mechanisms of ovarian cancer development and potential molecular targets for ovarian cancer treatment. ^

A systematic analysis of INO80 in DNA replication

ATP-dependent chromatin remodeling has been shown to be critical for transcription and DNA repair. However, the involvement of ATP-dependent chromatin remodeling in DNA replication remains poorly defined. Interestingly, we found that the INO80 chromatin-remodeling complex is directly involved in the DNA damage tolerance pathways activated during DNA replication. DNA damage tolerance is important for genomic stability and is controlled by formation of either mono-ubiquitinated or multi-ubiquitinated PCNA, which respectively induce error prone or error-free replication bypass of the lesions. In addition, homologous recombination (HR) mediated by the Rad51 pathway is also involved in the DNA damage tolerance pathways. ^ We found that INO80 is specifically recruited to replication origins during S phase in a genome-wide fashion. In addition, DNA combing analysis shows INO80 is required for the resumption of replication at stalled forks induced by methyl methane-sulfonate (MMS). Mechanistically, we find that INO80 is required for PCNA ubiquitination as well as for Rad51 mediated processing of replication forks after MMS treatment. Furthermore, chromatin immunoprecipitation at specific ARSs indicates INO80 is necessary for Rad18 and Rad51 recruitment to replication forks after MMS treatment. Moreover, 2D gel analysis shows INO80 is necessary to process Rad51 mediated intermediates at impeded replication forks. ^ In conclusion, our findings establish a novel role of a chromatin-remodeling complex in DNA damage tolerance pathways and suggest that chromatin remodeling is fundamentally important to ensure faithful replication of DNA and genome stability in eukaryotes. ^

Genetic variations in the PI3K-AKT-mTOR pathway and bladder cancer susceptibility and clinical outcomes

Bladder cancer is the fourth most common cancer in men in the United States. There is compelling evidence supporting that genetic variations contribute to the risk and outcomes of bladder cancer. The PI3K-AKT-mTOR pathway is a major cellular pathway involved in proliferation, invasion, inflammation, tumorigenesis, and drug response. Somatic aberrations of PI3K-AKT-mTOR pathway are frequent events in several cancers including bladder cancer; however, no studies have investigated the role of germline genetic variations in this pathway in bladder cancer. In this project, we used a large case control study to evaluate the associations of a comprehensive catalogue of SNPs in this pathway with bladder cancer risk and outcomes. Three SNPs in RAPTOR were significantly associated with susceptibility: rs11653499 (OR: 1.79, 95%CI: 1.24–2.60), rs7211818 (OR: 2.13, 95%CI: 1.35–3.36), and rs7212142 (OR: 1.57, 95%CI: 1.19–2.07). Two haplotypes constructed from these 3 SNPs were also associated with bladder cancer risk. In combined analysis, a significant trend was observed for increased risk with an increase in the number of unfavorable genotypes (P for trend<0.001). Classification and regression tree analysis identified potential gene-environment interactions between RPS6KA5 rs11653499 and smoking. In superficial bladder cancer, we found that PTEN rs1234219 and rs11202600, TSC1 rs7040593, RAPTOR rs901065, and PIK3R1 rs251404 were significantly associated with recurrence in patients receiving BCG. In muscle invasive and metastatic bladder cancer, AKT2 rs3730050, PIK3R1 rs10515074, and RAPTOR rs9906827 were associated with survival. Survival tree analysis revealed potential gene-gene interactions: patients carrying the unfavorable genotypes of PTEN rs1234219 and TSC1 rs704059 exhibited a 5.24-fold (95% CI: 2.44–11.24) increased risk of recurrence. In combined analysis, with the increasing number of unfavorable genotypes, there was a significant trend of higher risk of recurrence and death (P for trend<0.001) in Cox proportional hazard regression analysis, and shorter event (recurrence and death) free survival in Kaplan-Meier estimates (P log rank<0.001). This study strongly suggests that genetic variations in PI3K-AKT-mTOR pathway play an important role in bladder cancer development. The identified SNPs, if validated in further studies, may become valuable biomarkers in assessing an individual's cancer risk, predicting prognosis and treatment response, and facilitating physicians to make individualized treatment decisions. ^

Association of folate metabolic pathway genes with human spina bifida meningomyelocele

Neural tube defects including spina bifida meningomyelocele (SBMM) are common malformations of the brain and spinal cord, and include all abnormalities resulting from lack of closure of the developing neural tube during embryological development.^ The specific aims of this study were to determine if single nucleotide polymorphic variants (SNPs) in the folate/homocysteine metabolic pathway genes confer a risk for NTD susceptibility within this SBMM population.^ In completion of the first specific aim, two novel SNPs were identified in the FOLR1 gene in Chromosome 11of patients including one in non-coding exon 1 with a C → T transition at nucleotide position 71578317 and another in non-coding exon 3 with a T → G transversion at nucleotide position 71579123. It will be important to determine if these variants are present in the respective parents of these individuals. If they are in fact de novo variants, then these SNPs may be more likely to contribute to the birth defect.^ The second project aim was to analyze genotypes associated with SBMM risk by transmission disequilibrium tests (TDT) and association was detected on several SNPs across the folate metabolic pathway genes in this population. SNPs with significant RC-TDT values were found within the DHFR gene (rs1650723), the MTRR gene (rs327592), the FOLR2 gene (rs13908), four tightly linked variants in the FOLR3 gene (rs7925545, rs7926875, rs7926987, rs7926360) and a variant in the SLC19A1 gene (rs1888530). The product of each of these genes performs a vital function in the folate metabolic pathway. It is conceivable, therefore, that if the individual SNP or SNPs can be proven to perturb the function in some way that they may be involved in the disruption of folate metabolism and in the resulting birth defect. Validating the results of this study in other independent populations will further strengthen the evidence that dysfunction of folate enzymes and receptors may confer SBMM risk in humans. ^

A computerized statistical framework for coalescent analysis

Coalescent theory represents the most significant progress in theoretical population genetics in the past three decades. The coalescent theory states that all genes or alleles in a given population are ultimately inherited from a single ancestor shared by all members of the population, known as the most recent common ancestor. It is now widely recognized as a cornerstone for rigorous statistical analyses of molecular data from population [1]. The scientists have developed a large number of coalescent models and methods[2,3,4,5,6], which are not only applied in coalescent analysis and process, but also in today’s population genetics and genome studies, even public health. The thesis aims at completing a statistical framework based on computers for coalescent analysis. This framework provides a large number of coalescent models and statistic methods to assist students and researchers in coalescent analysis, whose results are presented in various formats as texts, graphics and printed pages. In particular, it also supports to create new coalescent models and statistical methods. ^

Statistical Methods for Differential Expressions of Genes Detected in Multiple-Condition Experiment of Microarray

Most studies of differential gene-expressions have been conducted between two given conditions. The two-condition experimental (TCE) approach is simple in that all genes detected display a common differential expression pattern responsive to a common two-condition difference. Therefore, the genes that are differentially expressed under the other conditions other than the given two conditions are undetectable with the TCE approach. In order to address the problem, we propose a new approach called multiple-condition experiment (MCE) without replication and develop corresponding statistical methods including inference of pairs of conditions for genes, new t-statistics, and a generalized multiple-testing method for any multiple-testing procedure via a control parameter C. We applied these statistical methods to analyze our real MCE data from breast cancer cell lines and found that 85 percent of gene-expression variations were caused by genotypic effects and genotype-ANAX1 overexpression interactions, which agrees well with our expected results. We also applied our methods to the adenoma dataset of Notterman et al. and identified 93 differentially expressed genes that could not be found in TCE. The MCE approach is a conceptual breakthrough in many aspects: (a) many conditions of interests can be conducted simultaneously; (b) study of association between differential expressions of genes and conditions becomes easy; (c) it can provide more precise information for molecular classification and diagnosis of tumors; (d) it can save lot of experimental resources and time for investigators.^