Application of stellar photometry to the analysis of microarray images

Improvements in the analysis of microarray images are critical for accurately quantifying gene expression levels. The acquisition of accurate spot intensities directly influences the results and interpretation of statistical analyses. This dissertation discusses the implementation of a novel approach to the analysis of cDNA microarray images. We use a stellar photometric model, the Moffat function, to quantify microarray spots from nylon microarray images. The inherent flexibility of the Moffat shape model makes it ideal for quantifying microarray spots. We apply our novel approach to a Wilms' tumor microarray study and compare our results with a fixed-circle segmentation approach for spot quantification. Our results suggest that different spot feature extraction methods can have an impact on the ability of statistical methods to identify differentially expressed genes. We also used the Moffat function to simulate a series of microarray images under various experimental conditions. These simulations were used to validate the performance of various statistical methods for identifying differentially expressed genes. Our simulation results indicate that tests taking into account the dependency between mean spot intensity and variance estimation, such as the smoothened t-test, can better identify differentially expressed genes, especially when the number of replicates and mean fold change are low. The analysis of the simulations also showed that overall, a rank sum test (Mann-Whitney) performed well at identifying differentially expressed genes. Previous work has suggested the strengths of nonparametric approaches for identifying differentially expressed genes. We also show that multivariate approaches, such as hierarchical and k-means cluster analysis along with principal components analysis, are only effective at classifying samples when replicate numbers and mean fold change are high. Finally, we show how our stellar shape model approach can be extended to the analysis of 2D-gel images by adapting the Moffat function to take into account the elliptical nature of spots in such images. Our results indicate that stellar shape models offer a previously unexplored approach for the quantification of 2D-gel spots. ^

MCAM/MUC18 REGULATES MELANOMA PROGRESSION BY MODULATING THE EXPRESSION OF ID-1

The acquisition of the metastatic melanoma phenotype is associated with increased expression of the melanoma cell adhesion molecule MCAM/MUC18 (CD146). However, the mechanism by which MUC18 contributes to melanoma metastasis remains unclear. Herein, we stably silenced MUC18 expression utilizing lentivirus-incorporated small hairpin RNA, in two metastatic melanoma cell lines, A375SM and C8161, and conducted cDNA microarray analysis. We identified and validated that the transcriptional regulator, Inhibitor of DNA Binding-1 (Id-1), previously shown to function as an oncogene in several malignancies, was downregulated by 5.6-fold following MUC18 silencing. Additionally, we found that MUC18 regulated Id-1 expression at the transcriptional level via ATF-3. Interestingly, ATF-3 was upregulated by 6.9 fold in our cDNA microarray analysis following MUC18 silencing. ChIP analysis showed increased binding of ATF-3 to the Id-1 promoter after MUC18 silencing, while mutation of the ATF-3 binding site on the Id-1 promoter increased Id-1 promoter activity in MUC18-silenced cells. These Data suggest that MUC18 silencing promotes inhibition of Id-1 expression by increasing ATF-3 expression and binding to the Id-1 promoter. Rescue of MUC18 reverted the expression of Id-1 and ATF-3, thus validating that they are not off-target effects of MUC18. To further assess the role of Id-1 in melanoma invasion and metastasis, we overexpressed Id-1 in MUC18-silenced cells. Overexpression of Id-1 in MUC18-silenced cells resulted in increased cell invasion, as well as increased expression and activity of MMP-2. Our data further reveal that Id-1 regulates MMP-2 at the transcriptional level through Sp1 and Ets-1. This is the first report to demonstrate that MUC18 does not act exclusively in cell adherence, but is also involved in cell signaling that regulates the expression of genes, such as Id-1 and ATF-3, thus contributing to the metastatic melanoma phenotype.

Collagen I modulates growth and gene expression in prostate cancer cells

Prostate cancer is the second leading cause of male cancer-related deaths in the United States. Interestingly, prostate cancer preferentially metastasizes to skeletal tissue. Once in the bone microenvironment, advanced prostate cancer becomes highly resistant to therapeutic modalities. Several factors, such as extracellular matrix (ECM) components, have been implicated in the spread and propagation of prostatic carcinoma. In these studies, we have utilized the PC3 cell line, derived from a human bone metastasis, to investigate the influence of the predominant bone ECM protein, type I collagen, on prostate cancer cell proliferation and gene expression. We have also initiated the design and production of ribozymes to specific gene targets that may influence prostate cancer bone metastasis. ^ Our results demonstrate that PC3 cells rapidly adhere and spread on collagen I to a greater degree than on fibronectin (FN) or poly-L-lysine (PLL). Flow cytometry analysis reveals the presence of the α1, α2 and α3 collagen binding integrin subunits. The use of antibody function blocking studies reveals that PC3 cells can utilize α2β 1 and α3β1 integrins to adhere to collagen I. Once plated on collagen I, the cells exhibit increased rates of proliferation compared with cells plated on FN or tissue culture plastic. Additionally, cells plated on collagen I show increased expression of proteins associated with progression through G1 phase of the cell cycle. Inhibitor studies point to a role for phosphatidylinositol 3-kinase (PI3K), MAP kinase (MAPK), and p70 S6 kinase in collagen I-mediated PC3 cell proliferation and cyclin D1 expression. To further characterize the effect of type I collagen on prostate cancer bone metastasis, we utilized a cDNA microarray strategy to monitor type I collagen-mediated changes in gene expression. Results of this analysis revealed a gene expression profile reflecting the increased proliferation occurring on type I collagen. Microarray analysis also revealed differences in the expression of specific gene targets that may impact on prostate cancer metastasis to bone. ^ As a result of our studies on the interaction of prostate cancer cells and the skeletal ECM, we sought to develop novel molecular tools for future gene therapy of functional knockdown experiments. To this end, we developed a series of ribozymes directed against the α2 integrin and at osteopontin, a protein implicated in the metastasis of various cancers, including prostate. These ribozymes should facilitate the future study of the mechanism of prostate cancer cell proliferation, and disease progression occurring at sites of skeletal metastasis where a type I collagen-based environment predominates. ^ Together these studies demonstrate the involvement of bone ECM proteins on prostate cancer cell proliferation and suggest that they may play a significant role on the growth of prostate metastases once in the bone microenvironment. ^

The role of coactivator associated arginine methyltransferase 1 (CARM1) in nuclear receptor mediated transcription

Arginine methylation has been implicated in the regulation of gene expression. The coactivator-associated arginine methyltransferase 1 (CARMI/PRMT4) binds the p160 family of steroid receptor coactivators (SRCs). This association enhances transcriptional activation by nuclear receptors. Here, we generated and characterized CARM1 knockout mice. Embryos with a targeted disruption of CARM1 are 35% smaller in size than the wild-type littermates and die perinatally. We also generated Carm1-/- and Carm1+/+ mouse embryonic fibroblasts and tested gene expression in response to estrogen. Estrogenresponsive gene expression was aberrant in Carm1-/- fibroblasts and embryos, thus emphasizing the role of arginine methylation as a transcription activation tag. We subsequently studied the role of CARM1 in estrogen signaling in viva in the mammary gland. Conditional knockout of CARM1 in mammary gland and Carml-1-embryonic mammary anlagen transplant experiments did not show any defects in growth and development of the glands. To further dissect the role of CARM1 in estrogen receptor mediated transactivation, we performed cDNA microarray and serial analysis of gene expression on Carm1-/- and Carm1+/+ embryos treated with the estrogen analog, DES. Our results indicate global changes in estrogen regulated genes as well as genes involved in lipid homeostasis. Marker genes for Peroxisome Proliferator Activated Receptor γ (PPARγ) activity, adipsin and aP2, are downregulated in the Carm1-/- embryos. Furthermore, OCT frozen sections of 18.5dpc embryos, processed simultaneously for oil red O staining to look for neutral fat, reveals greatly reduced brown fat accumulation in the Carm1-/- embryos in contrast to wild-type and gain-of-function Carm1 transgenic (ubiquitous) embryo. We used a well-established 3T3-L1 preadipocyte cell line to knockdown CARM1 by short hairpin RNA. 3T3-L1 cells with CARM1 knockdown showed greatly reduced potential to differentiate into mature lipid accumulating adipocytes upon administration of adipogenic stimuli. Ligand-dependent activation of reporter genes by the PPARγ receptor showed that PPRE-luciferase reporter activity was enhanced in the presence of CARM1, additionally, luciferase activity was reduced to background levels when enzyme dead CARM1 (CARM1-VLD) was used. Thus, in this study, we have identified novel pathways that use CARM1 as coactivator and showed that CARM1 functions as a key component of PPARγ receptor mediated gene expression. ^

Protease-activated receptor-1 signaling plays a major role in melanoma growth and metastasis

Overexpression of the thrombin receptor (Protease-Activated-Receptor-1), PAR-1, in cell lines and tissue specimens correlates with the metastatic potential of human melanoma. Utilizing lentiviral shRNA to stably silence PAR-1 in metastatic melanoma cell lines results in decreased tumor growth and lung metastasis in vivo. Since the use of viral technology is not ideal for clinical therapies, neutral liposomes (DOPC) were utilized as a delivery vehicle for PAR-1 siRNA. Our data suggest that PAR-1 siRNA-DOPC treatment by systemic delivery significantly decreases tumor growth and lung metastasis in nude mice. Concomitant decreases in angiogenic and invasive factors (IL-8, VEGF, MMP-2) were observed in PAR-1 siRNA-DOPC-treated mice. Utilizing a cDNA microarray platform, several novel PAR-1 downstream target genes were identified, including Connexin 43 (Cx-43) and Maspin. Cx-43, known to be involved in tumor cell diapedesis and attachment to endothelial cells, is decreased after PAR-1 silencing. Furthermore, the Cx-43 promoter activity was significantly inhibited in PAR-1-silenced cells suggesting transcriptional regulation of Cx-43 by PAR-1. ChIP analysis revealed a reduction in SP-1 and AP-1 binding to the Cx-43 promoter. Moreover, melanoma cell attachment to HUVEC was significantly decreased in PAR-1-silenced cells as well as in Cx-43 shRNA transduced cells. As both SP-1 and AP-1 transcription factors act as positive regulators of Cx-43, our data provide a novel mechanism for the regulation of Cx-43 expression by PAR-1. Maspin, a serine protease inhibitor with tumor-suppressor function, was found to be upregulated after PAR-1 silencing. Our results indicate that PAR-1 transcriptionally regulates Maspin, as the promoter activity was significantly increased after PAR-1 silencing. ChIP analysis revealed that silencing PAR-1 increased binding of Ets and c-Jun to the Maspin promoter. As Maspin was recently found to be a tumor-suppressor in melanoma by reducing the invasive capacity of melanoma cells, invasion assays revealed a decrease in invasion after PAR-1 silencing and in cells transduced with a Maspin expression vector. We propose that PAR-1 is key to the progression and metastasis of melanoma in part by regulating the expression of Cx-43 and Maspin. Taken together, we propose that PAR-1 is an attractive target for the treatment of melanoma.^

The molecular cloning and characterization of human heparanase cDNA and the immunochemical localization of heparanase in metastatic melanomas

Heparanase, an endo-$\beta$-D-glucuronidase, has been associated with melanoma metastasis. Polyclonal antibodies directed against the murine N-terminal heparanase peptide detected a M$\sb{\rm r}\sim 97,000$ protein upon SDS-polyacrylamide gel electrophoresis of mouse melanoma and human melanoma cell lysates. In an indirect immunocytochemical study, metastatic human A375-SM and mouse B16-BL6 melanoma cells were stained with the anti-heparanase antibodies. Heparanase antigen was localized in the cytoplasm of permeabilized melanoma cells as well as at the cell surface of unpermeabilized cells. Immunohistochemical staining of frozen sections from syngeneic mouse organs containing micrometastases of B16-BL6 melanoma demonstrated heparanase localized in metastatic melanoma cells, but not in adjacent normal tissues. Similar studies using frozen sections of malignant melanomas resected from patients indicated that heparanase is localized in invading melanoma cells, but not in adjacent connective tissues.^ Monoclonal antibodies directed against murine heparanase were developed and characterized. Monoclonal antibody 10E5, an IgM, precipitated and inhibitated the enzymatic activity of heparanase. A 2.6 kb cDNA was isolated from a human melanoma $\lambda$gt11 cDNA library using the monoclonal antibody 10E5. Heparan sulfate cleavage activity was detected in the lysogen lysates from E. Coli Y1089 infected with the $\lambda$gt11 cDNA and this activity was inhibited in the presence of 10-fold excess of heparin, a potent inhibitor of heparanase. The nucleotide sequence of the cDNA was determined and insignificant homology was found with the gene sequences currently known. The cDNA hybridized to a 3.2-3.4 kb mRNA in human A375 melanoma, WI-38 fibroblast, and THP-1 leukemia cells using Northern blots.^ Heparanase expression was examined using Western and Northern blots. In comparison to human A375-P melanoma cells, the quantity of 97,000 protein recognized by the polyclonal anti-heparanase antibodies doubled in the metastatic variant A375-SM cells and the quantity of 3.2-3.4 kb mRNA doubled in A375MetMix, a metastatic variant similar to A375-SM cells. In B16 murine melanoma cell, the intensity of the 97,000 protein increased more than 2 times comparing with B16-F1 cells. The extent in the increase of the protein and the mRNA levels is comparable to the change of heparanase activity observed in those cells.^ In summary, the studies suggest that (a) the N-terminus of the heparanase molecule in mouse and human is antigenically related; (b) heparanase antigens are localized at the cell surface and in the cytoplasm of metastatic human and mouse melanoma cells; (c) heparanase antigens are localized in invasive and metastatic murine and human melanomas in vivo, but not in adjacent normal tissues; (d) heparanase molecule appeared to be differentially expressed at the transcriptional as well as at the translational level; and (e) the size of human heparanase mRNA is 3.2-3.4 kilobase. ^

cDNA cloning and expression of a novel cell surface-expressed \b heparan sulfate/heparin \b interacting \b protein (HIP) from human cell lines and functional studies of a synthetic HIP peptide

Heparan sulfate proteoglycans and their corresponding binding sites have been suggested to play an important role during the initial attachment of blastocysts to uterine epithelium and human trophoblastic cell lines to uterine epithelial cell lines. Previous studies on RL95 cells, a human uterine epithelial cell line, characterized a single class of cell surface heparin/heparan sulfate (HP/HS)-binding sites. Three major HP/HS-binding peptide fragments were isolated from RL95 cell surfaces by tryptic digestion and partial amino-terminal amino acid sequence from each peptide fragment was obtained. In the current study, using the approaches of reverse transcription-polymerase chain reaction and cDNA library screening, a novel cell surface $\rm\underline{H}$P/HS $\rm\underline{i}$nteracting $\rm\underline{p}$rotein (HIP) has been isolated from RL95 cells. The full-length cDNA of HIP encodes a protein of 259 amino acids with a calculated molecular weight of 17,754 Da and pI of 11.75. Transfection of HIP cDNA into NIH-3T3 cells demonstrated cell surface expression and a size similar to that of HIP expressed by human cells. Predicted amino acid sequence indicates that HIP lacks a membrane spanning region and has no consensus sites for glycosylation. Northern blot analysis detected a single transcript of 1.3 kb in both total RNA and poly(A$\sp+$) RNA. Examination of human cell lines and normal tissues using both Northern blot and Western blot analysis revealed that HIP is differentially expressed in a variety of human cell lines and normal tissues, but absent in some cell lines examined. HIP has about 80% homology, at the level of both mRNA and protein, to a rodent protein, designated as ribosomal protein L29. Thus, members of the L29 family may be displayed on cell surfaces where they participate in HP/HS binding events. Studies on a synthetic peptide derived from HIP demonstrate that HIP peptide binds HS/HP with high selectivity and has high affinity (Kd = 10 nM) for a subset of polysaccharides found in commercial HIP preparations. Moreover, HIP peptide also binds certain forms of cell surface, but not secreted or intracellular. HS expressed by RL95 and JAR cells. This peptide supports the attachment of several human trophoblastic cell lines and a variety of mammalian adherent cell lines in a HS-dependent fashion. Furthermore, studies on the subset of HP specifically recognized by HIP peptide indicate that this high-affinity HP (HA-HP) has a larger median MW and a greater negative charge density than bulk HP. The minimum size of oligosaccharide required to bind to HIP peptide with high affinity is a septa- or octasaccharide. HA-HP also quantitatively binds to antithrombin-III (AT-III) with high affinity, indicating that HIP peptide and AT-III may recognize the same or similar oligosaccharide structure(s). Furthermore, HIP peptide antagonizes HP action and promotes blood coagulation in both factor Xa- and thrombin-dependent assays. Finally, HA-HP recognized by HP peptide is highly enriched with anticoagulant activity relative to bulk HP. Collectively, these results demonstrate that HIP may play a role in the HP/HS-involved cell-cell and cell-matrix interactions and recognizes a motif in HP similar or identical to that recognized by AT-III and therefore, may modulate blood coagulation. ^

A novel form of CYP3A from rat brain: cDNA cloning, protein characterization and gene regulation of cytochrome P450 3A9

One full length cDNA clone, designated 3aH15, was isolated from a rat brain cDNA library using a fragment of CYP3A2 cDNA as a probe. 3aH15 encoded a protein composed of 503 amino acid residues. The deduced amino acid sequence of 3aH15 was 92% identical to mouse Cyp3a-13 and had a 68.4% to 76.5% homology with the other reported rat CYP3A sequences. Clone 3aH15 was thus named CYP3A9 by Cytochrome P450 Nomenclature Committee. CYP3A9 seems to the major CYP3A isozyme expressed in rat brain. Sexual dimorphism of the expression of CYP3A9 was shown for the first time in rat brain as well as in rat liver. CYP3A9 appears to be female specific in rat liver based on the standards proposed by Kato and Yamazoe who defined sex specific expression of P450s as being a 10-fold or higher expression level in one sex compared with the other. CYP3A9 gene expression was inducible by estrogen treatment both in male and in female rats. Male rats treated with estrogen had a similar expression level of CYP3A9 mRNA both in the liver and brain. Ovariectomy of adult female rats drastically reduced the mRNA level of CYP3A9 which could be fully restored by estrogen replacement. On the other hand, only a two-fold induction of CYP3A9 expression by dexamethasone was observed in male liver and no significant induction of CYP3A9 mRNA was observed in female liver or in the brains. These results suggest that estrogen may play an important role in the female specific expression of the CYP3A9 gene and that CYP3A9 gene expression is regulated differently from other CYP3A isozymes. ^ P450 3A9 recombinant protein was expressed in E. coli using the pCWOri+ expression vector and the MALLLAVF amino terminal sequence modification. This construct gave a high level of expression (130 nmol P450 3A9/liter culture) and the recombinant protein of the modified P450 3A9 was purified to electrophoretic homogeneity (10.1 nmol P450/mg protein) from solubilized fractions using two chromatographic steps. The purified P450 3A9 protein was active towards the metabolism of many clinically important drugs such as imipramine, erythromycin, benzphetamine, ethylmorphine, chlorzoxazone, cyclosporine, rapamycin, etc. in a reconstituted system containing lipid and rat NADPH-P450 reductase. Although P450 3A9 was active towards the catabolism of testosterone, androstenedione, dehydroepiandrosterone (DHEA) and 17β-estradiol, P450 3A9 preferentially catalyzes the metabolism of progesterone to form four different hydroxylated products. Optimal reconstitution conditions for P450 3A9 activities required a lipid mixture and GSH. The possible mechanisms of the stimulatory effects of GSH on P450 3A9 activities are discussed. Sexually dimorphic expression of P450 3A9 in the brain and its involvement in many neuroactive drugs as well as neurosteroids suggest the possible role of P450 3A9 in some mental disorders and brain functions. ^

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

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

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

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

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

Cloning and characterization of the full length cDNA and promoter of mouse tissue transglutaminase

Transglutaminases are a family of enzymes that catalyze the covalent cross-linking of proteins through the formation of $\varepsilon$-($\gamma$-glutaminyl)-lysyl isopeptide bonds. Tissue transglutaminase (Tgase) is an intracellular enzyme which is expressed in terminally differentiated and senescent cells and also in cells undergoing apoptotic cell death. To characterize this enzyme and examine its relationship with other members of the transglutaminase family, cDNAs, the first two exons of the gene and 2 kb of the 5$\sp\prime$ flanking region, including the promoter, were isolated. The full length Tgase transcript consists of 66 bp of 5$\sp\prime$-UTR (untranslated) sequence, an open reading frame which encodes 686 amino acids and 1400 bp of 3$\sp\prime$-UTR sequence. Alignment of the deduced Tgase protein sequence with that of other transglutaminases revealed regions of strong homology, particularly in the active site region.^ The Tgase cDNA was used to isolate and characterize a genomic clone encompassing the 5$\sp\prime$ end of the mouse Tgase gene. The transcription start site was defined using genomic and cDNA clones coupled with S1 protection analysis and anchored PCR. This clone includes 2.3 kb upstream of the transcription start site and two exons that contain the first 256 nucleotides of the mouse Tgase cDNA sequence. The exon intron boundaries have been mapped and compared with the exon intron boundaries of three members of the transglutaminase family: human factor XIIIa, the human keratinocyte transglutaminase and human erythrocyte band 4.1. Tissue Tgase exon II is similar to comparable exons of these genes. However, exon I bears no resemblance with any of the other transglutaminase amino terminus exons.^ Previous work in our laboratory has shown that the transcription of the Tgase gene is directly controlled by retinoic acid and retinoic acid receptors. To identify the region of the Tgase gene responsible for regulating its expression, fragments of the Tgase promoter and 5$\sp\prime$-flanking region were cloned into the chloramphenicol actetyl transferase (CAT) reporter constructs. Transient transfection experiments with these constructs demonstrated that the upstream region of Tgase is a functional promoter which contains a retinoid response element within a 1573 nucleotide region spanning nucleotides $-$252 to $-$1825. ^

The cDNA cloning of murine HIP /L29 and the functional analysis of HIP/L29 domains

Human heparin/heparan sulfate interacting protein/L29 (HIP/L29) is a heparin/heparan sulfate (Hp/HS) binding protein found in many adult human tissues. Potential functions of this protein are promotion of embryo adhesion, modulation of blood coagulation, and control of cell growth. While these activities are diverse, the ability of human HIP/L29 to interact with Hp/HS at the cell surface may be a unifying mechanism of action since Hp/HS influences all of these processes. A murine ortholog has been identified that has 78.8% homology over the entire sequence and identity over the N-terminal 64 amino acids when compared to human HIP/L29. Northern, Western, and immunohistochemical analysis shows that murine HIP/L29 mRNA and protein are expressed in a tissue specific manner. Murine HIP/L29 is enriched in the membrane fraction of NmuMG cells where it is eluted with high salt, suggesting that it is a peripheral membrane protein. The ability of murine HIP/L29 to bind Hp is verified by studies using native and recombinant forms of murine HIP/L29. A synthetic peptide (HIP peptide-2) derived from the identical N-terminal region of HIP/L29 proteins was tested for the ability to bind Hp and support cell adhesion. This peptide was chosen because it conforms to a proposed consensus sequence for Hp/HS binding peptides. HIP peptide-2 binds Hp in a dose-dependent, saturable, and selective manner and supports Hp-dependent cell adhesion. However, a scrambled form of this peptide displayed similar activities indicating a lack of peptide sequence specificity required for activity. Lastly, an unbiased approach was used to identify sequences within human and mouse HIP/L29 proteins necessary for Hp/HS binding. A panel of recombinant proteins was made that collectively are deficient in every human HIP/L29 domain. The activities of these deletion mutants and recombinant murine HIP/L29 were compared to the activity of recombinant human HIP/L29 in a number of assays designed to look at differences in the ability to bind Hp/HS. These studies suggest that each domain within human HIP/L29 is important for binding to Hp/HS and divergences in the C-terminus of human and mouse HIP/L29 account for a decrease in murine HIP/L29 affinity for Hp/HS. It is apparent that multiple domains within human and mouse HIP/L29 contribute to the function of Hp/HS binding. The interaction of multiple HIP/L29 domains with Hp/HS will influence the biological activity of HIP/L29 proteins. ^