Methods to improve gene signal: Application to cDNA microarrays

Microarrays are high throughput biological assays that allow the screening of thousands of genes for their expression. The main idea behind microarrays is to compute for each gene a unique signal that is directly proportional to the quantity of mRNA that was hybridized on the chip. A large number of steps and errors associated with each step make the generated expression signal noisy. As a result, microarray data need to be carefully pre-processed before their analysis can be assumed to lead to reliable and biologically relevant conclusions. This thesis focuses on developing methods for improving gene signal and further utilizing this improved signal for higher level analysis. To achieve this, first, approaches for designing microarray experiments using various optimality criteria, considering both biological and technical replicates, are described. A carefully designed experiment leads to signal with low noise, as the effect of unwanted variations is minimized and the precision of the estimates of the parameters of interest are maximized. Second, a system for improving the gene signal by using three scans at varying scanner sensitivities is developed. A novel Bayesian latent intensity model is then applied on these three sets of expression values, corresponding to the three scans, to estimate the suitably calibrated true signal of genes. Third, a novel image segmentation approach that segregates the fluorescent signal from the undesired noise is developed using an additional dye, SYBR green RNA II. This technique helped in identifying signal only with respect to the hybridized DNA, and signal corresponding to dust, scratch, spilling of dye, and other noises, are avoided. Fourth, an integrated statistical model is developed, where signal correction, systematic array effects, dye effects, and differential expression, are modelled jointly as opposed to a sequential application of several methods of analysis. The methods described in here have been tested only for cDNA microarrays, but can also, with some modifications, be applied to other high-throughput technologies. Keywords: High-throughput technology, microarray, cDNA, multiple scans, Bayesian hierarchical models, image analysis, experimental design, MCMC, WinBUGS.

Identification of novel genes expressed during rhabdomyosarcoma differentiation using cDNA microarrays

Rhabdomyosarcomas (RMS) are highly aggressive tumors that are thought to arise as a consequence of the regulatory disruption of the growth and differentiation of skeletal muscle progenitor cells. Normal myogenesis is characterized by the expression of the myogenic regulatory factor gene family but, despite their expression in RMS, these tumor cells fail to complete the latter stages of myogenesis. The RMS cell line RD-A was treated with 12-O-tetradecanoylphorbol-13-acetate to induce differentiation and cultured for 10 days. RNA was extracted on days 1, 3, 6, 8 and 10. A human skeletal muscle cDNA microarray was developed and used to analyze the global gene expression of RMS tumors over the time-course of differentiation. As a comparison, the genes identified were subsequently examined during the differentiated primary human skeletal muscle cultures. Prothymosin alpha (PTMA), and translocase of inner mitochondrial membrane 10 (Tim10), two genes not previously implicated in RMS, showed reduced expression during differentiation. Marked differences in the expression of PTMA and Tim10 were observed during the differentiation of human primary skeletal muscle cells. These results identify several new genes with potential roles in the myogenic arrest present in rhabdomyosarcoma. PTMA expression in RMS biopsy samples might prove to be an effective diagnostic marker for this disease.

Identification of genes differentially regulated by the P210 BCR/ABL1 fusion oncogene using cDNA microarrays

Objective: The t(9;22) translocation is associated with more than 95% of cases of chronic myeloid leukemia. The resulting fusion of the BCR and ABL1 loci produces the constitutively active BCR/ABL1 tyrosine kinase. A wide range of signal transduction molecules are activated by BCR/ABL1, including MYC, PI-3 kinase, and different STAT molecules. In contrast, relatively few genes are known to be regulated by BCR/ABL1 at the level of transcription.

Materials and Methods: In an effort to better understand the transcriptional program activated by BCR/ABL1, we used cDNA microarrays to evaluate the relative expression of approximately 6450 human genes in U937 myelomonocytic cells expressing P210 BCR/ABL1 via a tetracycline-inducible promoter.

Results: We confirmed the previously reported up-regulation of the PIM1 and JUN oncogenes by BCR/ABL1. In addition, we identified 59 more genes up-regulated by BCR/ABL1. Interestingly, roughly one third of these were genes previously reported to be interferon (IFN)-responsive, including the OAS1, IFIT1, IFI16, ISGF3G, and STAT1 genes. An additional seven BCR/ABL1-regulated genes were found to be IFN-responsive in U937 cells. The expression profile also included genes encoding transcription factors, kinases, and signal transduction molecules, as well as genes regulating cell growth, differentiation, apoptosis, and cell adhesion, features previously suggested to be affected by BCR/ABL1.

Conclusion: These observations shed novel insight into the mechanism of BCR/ABL1 action and provide a range of targets for further investigation.

Disease fingerprinting with cDNA microarrays reveals distinct gene expression profiles in lethal type 1 and type 2 cytokine-mediated inflammatory reactions

Development of polarized immune responses controls resistance and susceptibility to many microorganisms. However, studies of several infectious, allergic, and autoimmune diseases have shown that chronic type-1 and type-2 cytokine responses can also cause significant morbidity and mortality if left unchecked. We used mouse cDNA microarrays to molecularly phenotype the gene expression patterns that characterize two disparate but equally lethal forms of liver pathology that develop in Schistosoma mansoni infected mice polarized for type-1 and type-2 cytokine responses. Hierarchical clustering analysis identified at least three groups of genes associated with a polarized type-2 response and two linked with an extreme type-1 cytokine phenotype. Predictions about liver fibrosis,  apoptosis, and granulocyte recruitment and activation generated by the microarray studies were confirmed later by traditional biological assays. The data show that cDNA microarrays are useful not only for determining  coordinated gene expression profiles but are also highly effective for molecularly “fingerprinting” diseased tissues. Moreover, they illustrate the potential of genome-wide approaches for generating comprehensive views on the molecular and biochemical mechanisms regulating infectious  disease pathogenesis.

Evaluation of reference-based two-color methods for measurement of gene expression ratios using spotted cDNA microarrays

Abstract Background Spotted cDNA microarrays generally employ co-hybridization of fluorescently-labeled RNA targets to produce gene expression ratios for subsequent analysis. Direct comparison of two RNA samples in the same microarray provides the highest level of accuracy; however, due to the number of combinatorial pair-wise comparisons, the direct method is impractical for studies including large number of individual samples (e.g., tumor classification studies). For such studies, indirect comparisons using a common reference standard have been the preferred method. Here we evaluated the precision and accuracy of reconstructed ratios from three indirect methods relative to ratios obtained from direct hybridizations, herein considered as the gold-standard. Results We performed hybridizations using a fixed amount of Cy3-labeled reference oligonucleotide (RefOligo) against distinct Cy5-labeled targets from prostate, breast and kidney tumor samples. Reconstructed ratios between all tissue pairs were derived from ratios between each tissue sample and RefOligo. Reconstructed ratios were compared to (i) ratios obtained in parallel from direct pair-wise hybridizations of tissue samples, and to (ii) reconstructed ratios derived from hybridization of each tissue against a reference RNA pool (RefPool). To evaluate the effect of the external references, reconstructed ratios were also calculated directly from intensity values of single-channel (One-Color) measurements derived from tissue sample data collected in the RefOligo experiments. We show that the average coefficient of variation of ratios between intra- and inter-slide replicates derived from RefOligo, RefPool and One-Color were similar and 2 to 4-fold higher than ratios obtained in direct hybridizations. Correlation coefficients calculated for all three tissue comparisons were also similar. In addition, the performance of all indirect methods in terms of their robustness to identify genes deemed as differentially expressed based on direct hybridizations, as well as false-positive and false-negative rates, were found to be comparable. Conclusion RefOligo produces ratios as precise and accurate as ratios reconstructed from a RNA pool, thus representing a reliable alternative in reference-based hybridization experiments. In addition, One-Color measurements alone can reconstruct expression ratios without loss in precision or accuracy. We conclude that both methods are adequate options in large-scale projects where the amount of a common reference RNA pool is usually restrictive.

cDNA microarrays detect activation of a myogenic transcription program by the PAX3-FKHR fusion oncogene

Alveolar rhabdomyosarcoma is an aggressive pediatric cancer of striated muscle characterized in 60% of cases by a t(2;13)(q35;q14). This results in the fusion of PAX3, a developmental transcription factor required for limb myogenesis, with FKHR, a member of the forkhead family of transcription factors. The resultant PAX3-FKHR gene possesses transforming properties; however, the effects of this chimeric oncogene on gene expression are largely unknown. To investigate the actions of these transcription factors, both Pax3 and PAX3-FKHR were introduced into NIH 3T3 cells, and the resultant gene expression changes were analyzed with a murine cDNA microarray containing 2,225 elements. We found that PAX3-FKHR but not PAX3 activated a myogenic transcription program including the induction of transcription factors MyoD, Myogenin, Six1, and Slug as well as a battery of genes involved in several aspects of muscle function. Notable among this group were the growth factor gene Igf2 and its binding protein Igfbp5. Relevance of this model was suggested by verification that three of these genes (IGFBP5, HSIX1, and Slug) were also expressed in alveolar rhabdomyosarcoma cell lines. This study utilizes cDNA microarrays to elucidate the pattern of gene expression induced by an oncogenic transcription factor and demonstrates the profound myogenic properties of PAX3-FKHR in NIH 3T3 cells.

An evaluation of the performance of cDNA microarrays for detecting changes in global mRNA expression

The cDNA microarray is one technological approach that has the potential to accurately measure changes in global mRNA expression levels. We report an assessment of an optimized cDNA microarray platform to generate accurate, precise and reliable data consistent with the objective of using microarrays as an acquisition platform to populate gene expression databases. The study design consisted of two independent evaluations with 70 arrays from two different manufactured lots and used three human tissue sources as samples: placenta, brain and heart. Overall signal response was linear over three orders of magnitude and the sensitivity for any element was estimated to be 2 pg mRNA. The calculated coefficient of variation for differential expression for all non-differentiated elements was 12–14% across the entire signal range and did not vary with array batch or tissue source. The minimum detectable fold change for differential expression was 1.4. Accuracy, in terms of bias (observed minus expected differential expression ratio), was less than 1 part in 10 000 for all non-differentiated elements. The results presented in this report demonstrate the reproducible performance of the cDNA microarray technology platform and the methods provide a useful framework for evaluating other technologies that monitor changes in global mRNA expression.

Microarrays in molecular profiling of cancer - focus on head and neck squamous cell carcinoma

Microarrays have a wide range of applications in the biomedical field. From the beginning, arrays have mostly been utilized in cancer research, including classification of tumors into different subgroups and identification of clinical associations. In the microarray format, a collection of small features, such as different oligonucleotides, is attached to a solid support. The advantage of microarray technology is the ability to simultaneously measure changes in the levels of multiple biomolecules. Because many diseases, including cancer, are complex, involving an interplay between various genes and environmental factors, the detection of only a single marker molecule is usually insufficient for determining disease status. Thus, a technique that simultaneously collects information on multiple molecules allows better insights into a complex disease. Since microarrays can be custom-manufactured or obtained from a number of commercial providers, understanding data quality and comparability between different platforms is important to enable the use of the technology to areas beyond basic research. When standardized, integrated array data could ultimately help to offer a complete profile of the disease, illuminating mechanisms and genes behind disorders as well as facilitating disease diagnostics. In the first part of this work, we aimed to elucidate the comparability of gene expression measurements from different oligonucleotide and cDNA microarray platforms. We compared three different gene expression microarrays; one was a commercial oligonucleotide microarray and the others commercial and custom-made cDNA microarrays. The filtered gene expression data from the commercial platforms correlated better across experiments (r=0.78-0.86) than the expression data between the custom-made and either of the two commercial platforms (r=0.62-0.76). Although the results from different platforms correlated reasonably well, combining and comparing the measurements were not straightforward. The clone errors on the custom-made array and annotation and technical differences between the platforms introduced variability in the data. In conclusion, the different gene expression microarray platforms provided results sufficiently concordant for the research setting, but the variability represents a challenge for developing diagnostic applications for the microarrays. In the second part of the work, we performed an integrated high-resolution microarray analysis of gene copy number and expression in 38 laryngeal and oral tongue squamous cell carcinoma cell lines and primary tumors. Our aim was to pinpoint genes for which expression was impacted by changes in copy number. The data revealed that especially amplifications had a clear impact on gene expression. Across the genome, 14-32% of genes in the highly amplified regions (copy number ratio >2.5) had associated overexpression. The impact of decreased copy number on gene underexpression was less clear. Using statistical analysis across the samples, we systematically identified hundreds of genes for which an increased copy number was associated with increased expression. For example, our data implied that FADD and PPFIA1 were frequently overexpressed at the 11q13 amplicon in HNSCC. The 11q13 amplicon, including known oncogenes such as CCND1 and CTTN, is well-characterized in different type of cancers, but the roles of FADD and PPFIA1 remain obscure. Taken together, the integrated microarray analysis revealed a number of known as well as novel target genes in altered regions in HNSCC. The identified genes provide a basis for functional validation and may eventually lead to the identification of novel candidates for targeted therapy in HNSCC.

Importance of replication in microarray gene expression studies: Statistical methods and evidence from repetitive cDNA hybridizations

We present statistical methods for analyzing replicated cDNA microarray expression data and report the results of a controlled experiment. The study was conducted to investigate inherent variability in gene expression data and the extent to which replication in an experiment produces more consistent and reliable findings. We introduce a statistical model to describe the probability that mRNA is contained in the target sample tissue, converted to probe, and ultimately detected on the slide. We also introduce a method to analyze the combined data from all replicates. Of the 288 genes considered in this controlled experiment, 32 would be expected to produce strong hybridization signals because of the known presence of repetitive sequences within them. Results based on individual replicates, however, show that there are 55, 36, and 58 highly expressed genes in replicates 1, 2, and 3, respectively. On the other hand, an analysis by using the combined data from all 3 replicates reveals that only 2 of the 288 genes are incorrectly classified as expressed. Our experiment shows that any single microarray output is subject to substantial variability. By pooling data from replicates, we can provide a more reliable analysis of gene expression data. Therefore, we conclude that designing experiments with replications will greatly reduce misclassification rates. We recommend that at least three replicates be used in designing experiments by using cDNA microarrays, particularly when gene expression data from single specimens are being analyzed.

Delineating developmental and metabolic pathways in vivo by expression profiling using the RIKEN set of 18,816 full-length enriched mouse cDNA arrays

We have systematically characterized gene expression patterns in 49 adult and embryonic mouse tissues by using cDNA microarrays with 18,816 mouse cDNAs. Cluster analysis defined sets of genes that were expressed ubiquitously or in similar groups of tissues such as digestive organs and muscle. Clustering of expression profiles was observed in embryonic brain, postnatal cerebellum, and adult olfactory bulb, reflecting similarities in neurogenesis and remodeling. Finally, clustering genes coding for known enzymes into 78 metabolic pathways revealed a surprising coordination of expression within each pathway among different tissues. On the other hand, a more detailed examination of glycolysis revealed tissue-specific differences in profiles of key regulatory enzymes. Thus, by surveying global gene expression by using microarrays with a large number of elements, we provide insights into the commonality and diversity of pathways responsible for the development and maintenance of the mammalian body plan.

Screening insertion libraries for mutations in many genes simultaneously using DNA microarrays

We describe a method to screen pools of DNA from multiple transposon lines for insertions in many genes simultaneously. We use thermal asymmetric interlaced–PCR, a hemispecific PCR amplification protocol that combines nested, insertion-specific primers with degenerate primers, to amplify DNA flanking the transposons. In reconstruction experiments with previously characterized Arabidopsis lines carrying insertions of the maize Dissociation (Ds) transposon, we show that fluorescently labeled, transposon-flanking fragments overlapping ORFs hybridize to cognate expressed sequence tags (ESTs) on a DNA microarray. We further show that insertions can be detected in DNA pools from as many as 100 plants representing different transposon lines and that all of the tested, transposon-disrupted genes whose flanking fragments can be amplified individually also can be detected when amplified from the pool. The ability of a transposon-flanking fragment to hybridize declines rapidly with decreasing homology to the spotted DNA fragment, so that only ESTs with >90% homology to the transposon-disrupted gene exhibit significant cross-hybridization. Because thermal asymmetric interlaced–PCR fragments tend to be short, use of the present method favors recovery of insertions in and near genes. We apply the technique to screening pools of new Ds lines using cDNA microarrays containing ESTs for ≈1,000 stress-induced and -repressed Arabidopsis genes.

Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis

We have constructed cDNA microarrays for soybean (Glycine max L. Merrill), containing approximately 4,100 Unigene ESTs derived from axenic roots, to evaluate their application and utility for functional genomics of organ differentiation in legumes. We assessed microarray technology by conducting studies to evaluate the accuracy of microarray data and have found them to be both reliable and reproducible in repeat hybridisations. Several ESTs showed high levels (>50 fold) of differential expression in either root or shoot tissue of soybean. A small number of physiologically interesting, and differentially expressed sequences found by microarray analysis were verified by both quantitative real-time RT-PCR and Northern blot analysis. There was a linear correlation (r(2) = 0.99, over 5 orders of magnitude) between microarray and quantitative real-time RT-PCR data. Microarray analysis of soybean has enormous potential not only for the discovery of new genes involved in tissue differentiation and function, but also to study the expression of previously characterised genes, gene networks and gene interactions in wild-type, mutant or transgenic; plants.

Systematic expression profiling of the mouse transcriptome using RIKEN cDNA mircoarrays

The number of known mRNA transcripts in the mouse has been greatly expanded by the RIKEN Mouse Gene Encyclopedia project. Validation of their reproducible expression in a tissue is an important contribution to the study of functional genomics. In this report, we determine the expression profile of 57,931 clones on 20 mouse tissues using cDNA microarrays. Of these 57,931 clones, 22,928 clones correspond to the FANTOM2 clone set. The set represents 20,234 transcriptional units (TUs) out of 33,409 TUs in the FANTOM2 set. We identified 7206 separate clones that satisfied stringent criteria for tissue-specific expression. Gene Ontology terms were assigned for these 7206 clones, and the proportion of 'molecular function' ontology for each tissue-specific clone was examined. These data will provide insights into the function of each tissue. Tissue-specific gene expression profiles obtained using our cDNA microarrays were also compared with the data extracted from the GNF Expression Atlas based on Affymetrix microarrays. One major outcome of the RIKEN transcriptome analysis is the identification of numerous nonprotein-coding mRNAs. The expression profile was also used to obtain evidence of expression for putative noncoding RNAs. In addition, 1926 clones (70%) of 2768 clones that were categorized as unknown EST, and 1969 (58%) clones of 3388 clones that were categorized as unclassifiable were also shown to be reproducibly expressed.