Comparison of multiple comparison methods for identifying differential gene expression in simulated and real papillary thyroid cancer microarray data

The difficulty of detecting differential gene expression in microarray data has existed for many years. Several correction procedures try to avoid the family-wise error rate in multiple comparison process, including the Bonferroni and Sidak single-step p-value adjustments, Holm's step-down correction method, and Benjamini and Hochberg's false discovery rate (FDR) correction procedure. Each multiple comparison technique has its advantages and weaknesses. We studied each multiple comparison method through numerical studies (simulations) and applied the methods to the real exploratory DNA microarray data, which detect of molecular signatures in papillary thyroid cancer (PTC) patients. According to our results of simulation studies, Benjamini and Hochberg step-up FDR controlling procedure is the best process among these multiple comparison methods and we discovered 1277 potential biomarkers among 54675 probe sets after applying the Benjamini and Hochberg's method to PTC microarray data.^

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

Repression of tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) but not its receptors during oral cancer progression.

BACKGROUND: TRAIL plays an important role in host immunosurveillance against tumor progression, as it induces apoptosis of tumor cells but not normal cells, and thus has great therapeutic potential for cancer treatment. TRAIL binds to two cell-death-inducing (DR4 and DR5) and two decoy (DcR1, and DcR2) receptors. Here, we compare the expression levels of TRAIL and its receptors in normal oral mucosa (NOM), oral premalignancies (OPM), and primary and metastatic oral squamous cell carcinomas (OSCC) in order to characterize the changes in their expression patterns during OSCC initiation and progression. METHODS: DNA microarray, immunoblotting and immunohistochemical analyses were used to examine the expression levels of TRAIL and its receptors in oral epithelial cell lines and in archival tissues of NOM, OPM, primary and metastatic OSCC. Apoptotic rates of tumor cells and tumor-infiltrating lymphocytes (TIL) in OSCC specimens were determined by cleaved caspase 3 immunohistochemistry. RESULTS: Normal oral epithelia constitutively expressed TRAIL, but expression was progressively lost in OPM and OSCC. Reduction in DcR2 expression levels was noted frequently in OPM and OSCC compared to respective patient-matched uninvolved oral mucosa. OSCC frequently expressed DR4, DR5 and DcR1 but less frequently DcR2. Expression levels of DR4, DR5 and DcR1 receptors were not significantly altered in OPM, primary OSCC and metastatic OSCC compared to patient-matched normal oral mucosa. Expression of proapoptotic TRAIL-receptors DR4 and DR5 in OSCC seemed to depend, at least in part, on whether or not these receptors were expressed in their parental oral epithelia. High DR5 expression in primary OSCC correlated significantly with larger tumor size. There was no significant association between TRAIL-R expression and OSSC histology grade, nodal status or apoptosis rates of tumor cells and TIL. CONCLUSION: Loss of TRAIL expression is an early event during oral carcinogenesis and may be involved in dysregulation of apoptosis and contribute to the molecular carcinogenesis of OSCC. Differential expressions of TRAIL receptors in OSCC do not appear to play a crucial role in their apoptotic rate or metastatic progression.

New aspects of estrogen action in the endometrium: Reciprocal interplay with retinoic acid pathway and a novel estrogen-regulated gene

The uterine endometrium is a major target for the estrogen. However, the molecular basis of estrogen action in the endometrium is largely unknown. I have used two approaches to study the effects of estrogen on the endometrium. One approach involved the study of the interaction between estrogen and retinoic acid (RA) pathways in the endometrium. I have demonstrated that estrogen administration to rodents and estrogen replacement therapy (ERT) in postmenopausal women selectively induced the endometrial expression of retinaldehyde dehydrogenase II (RALDH2), a critical enzyme of RA biosynthesis. RALDH2 was expressed exclusively in the stromal cells, especially in the stroma adjacent to the luminal and glandular epithelia. The induction of RALDH2 by estrogen required estrogen receptor and occurred via a direct increase in RALDH2 transcription. Among the three RA receptors, estrogen selectively induced the expression of RARα. In parallel, estrogen also increased the utilization of all-trans retinol (the substrate for RA biosynthesis) and the expression of two RA-regulated marker genes, cellular retinoic acid binding protein II (CRABP2) and tissue transglutaminase (tTG) in the endometrium. Thus estrogen coordinately upregulated both the production and signaling of RA in both the rodent and human endometrium. This coordinate upregulation of RA system appeared to play a role in counterbalancing the stimulatory effects of estrogen on the endometrium, since the depletion of endogenous RA in mice led to an increase in estrogen-stimulated stromal proliferation and endometrial Akt phosphorylation. In addition, I have also used a systematic approach (DNA microarray) to categorize genes and pathways affected by the ERT in the endometrium of postmenopausal women and identified a novel estrogen-regulated gene EIG121. EIG121 was exclusively expressed in the glandular epithelial cells of the endometrium and induced by estrogen in vivo and in cultured cell lines. Compared with the normal endometrium, EIG121 was highly overexpressed in type 1 endometrial cancer, but profoundly suppressed in type 2 endometrial tumors. Taken together, these studies suggested that estrogen regulates the expression of many genes of both the pro-proliferative and anti-proliferative pathways and the abnormality of these pathways may increase the risks for estrogen-dependent endometrial hyperplasia and endometrial cancer. ^

Regulation and transcript analysis of the ceroperon responsible for quorum sensing in Rhodobacter sphaeroides 2.4.1

Rhodobacter sphaeroides 2.4.1 is a Gram negative facultative photoheterotrophic bacterium that has been shown to have an N-acyl homoserine lactone-based quorum sensing system called cer for c&barbelow;ommunity e&barbelow;scape r&barbelow;esponse. The cer ORFs are cerR, the transcriptional regulator, cerI, the autoinducer synthase and cerA , whose function is unknown. The autoinducer molecule, 7,8- cis-N-(tetradecenoyl) homoserine lactone, has been characterized. The objective of this study was to identify an environmental stimulus that influences the regulation of cerRAI and, to characterize transcription of the cer operon. ^ A cerR::lacZ transcriptional fusion was made and β-Galactosidase assays were performed in R. sphaeroides 2.4.1 strains, wild type, AP3 (CerI−) and AP4 (CerR−). The cerR::lacZ β-Galactosidase assays were used as an initial survey of the mode of regulation of the Cer system. A cerA::lacZ translational fusion was created and was used to show that cerA can be translated. The presence of 7,8-cis-N-(tetradecenoyl) homoserine lactone was detected from R. sphaeroides strains wild type and AP4 (CerR−) using a lasR::lacZ translational fusion autoinducer bioassay. The cerR::lacZ transcriptional fusion in R. sphaeroides 2.4.1 wild type was tested under different environmental stimuli, such as various carbon sources, oxygen tensions, light intensities and culture media to determine if they influence transcription of the cer ORFs. Although lacZ assay data implicated high light intensity at 100 W/m2 to stimulate cer transcription, quantitative Northern RNA data of the cerR transcript showed that low light intensity at 3 W/m2 is at least one environmental stimulus that induces cer transcription. This finding was supported by DNA microarray analysis. Northern analysis of the cerRAI transcript provided evidence that the cer ORFs are co-transcribed, and that the cer operon contains two additional genes. Bioinformatics was used to identify genes that may be regulated by the Cer system by identifying putative lux box homologue sequences in the presumed promoter region of these genes. Genes that were identified were fliQ, celB and calsymin, all implicated in interacting with plants. Primer extension was used to help localize cis-elements in the promoter region. The cerR::lacZ transcriptional fusion was monitored in a subset of different global DNA binding transcriptional regulator mutant strains of R. sphaeroides 2.4.1. Those regulators involved in maintaining an anaerobic photosynthetic lifestyle appeared to have an effect. Collectively, the data imply that R. sphaeroides 2.4.1 activates the Cer system when grown anaerobic photosynthetically at low light intensity, 3 W/m2, and it may be involved in an interaction with plants. ^

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. ^

DNA methylation inhibits p53 mediated survivin repression in cancer cells

Survivin (BIRC5) is a member of the Inhibitor of Apoptosis (IAP) gene family and functions as a chromosomal passenger protein as well as a mediator of cell survival. Survivin is widely expressed during embryonic development then becomes transcriptionally silent in most highly differentiated adult tissues. It is also overexpressed in virtually every type of tumor. The survivin promoter contains a canonical CpG island that has been described as epigenetically regulated by DNA methylation. We observed that survivin is overexpressed in high grade, poorly differentiated endometrial tumors, and we hypothesized that DNA hypomethylation could explain this expression pattern. Surprisingly, methylation specific PCR and bisulfite pyrosequencing analysis showed that survivin was hypermethylated in endometrial tumors and that this hypermethylation correlated with increased survivin expression. We proposed that methylation could activate survivin expression by inhibit the binding of a transcriptional repressor. ^ The tumor suppressor protein p53 is a well documented transcriptional repressor of survivin and examination of the survivin promoter showed that the p53 binding site contains 3 CpG sites which often become methylated in endometrial tumors. To determine if methylation regulates survivin expression, we treated HCT116 cells with decitabine, a demethylation agent, and observed that survivin transcript and protein levels were significantly repressed following demethylation in a p53 dependent manner. Subsequent binding studies confirmed that DNA methylation inhibited the binding of p53 protein to its binding site in the survivin promoter. ^ We are the first to report this novel mechanism of epigenetic regulation of survivin. We also conducted microarray analysis which showed that many other cancer relevant genes may also be regulated in this manner. While demethylation agents are traditionally thought to inhibit cancer cell growth by reactivating tumor suppressors, our results indicate that an additional important mechanism is to decrease the expression of oncogenes. ^

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. ^