PubMed Health⌕ Search

Biomedical subjects

Mario Medvedovic

Publications and source records attributed to Mario Medvedovic.

At least 19 recordsLinked to original sources

Conditional activation of RET/PTC3 and BRAFV600E in thyroid cells is associated with gene expression profiles that predict a preferential role of BRAF in extracellular matrix remodeling.

Papillary thyroid cancers (PTC) are associated with nonoverlapping mutations of genes coding for mitogen-activated protein kinase signaling effectors (i.e., the TK receptors RET or NTRK and the signaling proteins RAS and BRAF). We examined the pattern of gene expression after activation of these oncoproteins in thyroid PCCL3 cells, with the goal of identifying pathways or gene subsets that may account for the phenotypic differences observed in human cancers. We hybridized cDNA from cells treated with or without doxycycline to induce expression of BRAF(V600E), RET/PTC3, or RET/PTC3 with small interfering RNA-mediated knockdown of BRAF, respectively, to slides arrayed with a rat 70-mer oligonucleotide library consisting of 27,342 oligos. Among the RET/PTC3-induced genes, 2,552 did not require BRAF as they were similarly regulated by RET/PTC3 with or without BRAF knockdown and not by expression of BRAF(V600E). Immune response and IFN-related genes were highly represented in this group. About 24% of RET/PTC3-regulated genes were BRAF dependent, as they were similarly modified by RET/PTC3 and BRAF(V600E) but not in cells expressing RET/PTC3 with knockdown of BRAF. A gene cluster coding for components of the mitochondrial electron transport chain pathway was down-regulated in this group, potentially altering regulation of cell viability. Metalloproteinases were also preferentially induced by BRAF, particularly matrix metalloproteinase 3 (MMP3), MMP9, and MMP13. Accordingly, conditional expression of BRAF was associated with markedly increased invasion into Matrigel compared with cells expressing RET/PTC3. The preferential induction of MMPs by BRAF could explain in part the more invasive behavior of thyroid cancers with BRAF mutations.

Animals↗

Context-specific infinite mixtures for clustering gene expression profiles across diverse microarray dataset.

MOTIVATION: Identifying groups of co-regulated genes by monitoring their expression over various experimental conditions is complicated by the fact that such co-regulation is condition-specific. Ignoring the context-specific nature of co-regulation significantly reduces the ability of clustering procedures to detect co-expressed genes due to additional 'noise' introduced by non-informative measurements. RESULTS: We have developed a novel Bayesian hierarchical model and corresponding computational algorithms for clustering gene expression profiles across diverse experimental conditions and studies that accounts for context-specificity of gene expression patterns. The model is based on the Bayesian infinite mixtures framework and does not require a priori specification of the number of clusters. We demonstrate that explicit modeling of context-specificity results in increased accuracy of the cluster analysis by examining the specificity and sensitivity of clusters in microarray data. We also demonstrate that probabilities of co-expression derived from the posterior distribution of clusterings are valid estimates of statistical significance of created clusters. AVAILABILITY: The open-source package gimm is available at http://eh3.uc.edu/gimm.

Algorithms↗

Gene expression and discovery during lens regeneration in mouse: regulation of epithelial to mesenchymal transition and lens differentiation.

PURPOSE: It has been shown that after extracapsular lens removal by anterior capsulotomy in the mouse, the lens can be regenerated. However, as the capsular bag is filled with fibers, epithelial to mesenchymal transition (EMT), an event which is common after cataract surgery as well, takes place during early stages. This study, using a unique mouse model, was undertaken to identify novel regulators and networks in order to more clearly understand secondary cataracts at the molecular level. METHODS: We examined global gene expression via microarray analysis of mouse lens regeneration after extracapsular surgery. Gene expression at different times after surgery was correlated with the processes of EMT, which is seen in the initial stages of regeneration, and lens fiber differentiation, which occurs later. RESULTS: Several notable patterns were observed from the gene clustering data. It was obvious from the analysis that initially there is a response to injury, extensive matrix remodeling, and severe downregulation of genes encoding lens structural proteins. The patterns returned gradually to normal three weeks after surgery. New genes were identified from the clustering results that might be potential regulators of EMT and lens differentiation. CONCLUSIONS: With this approach, we demonstrated the utility of a mouse model to study secondary cataracts at the molecular level. Extension of these studies in mice with known mutations affecting EMT or lens differentiation should allow the identification of the crucial molecular players that could lead to better treatments of secondary cataracts.

Animals↗

A new method to remove hybridization bias for interspecies comparison of global gene expression profiles uncovers an association between mRNA sequence divergence and differential gene expression in Xenopus.

The recent sequencing of a large number of Xenopus tropicalis expressed sequences has allowed development of a high-throughput approach to study Xenopus global RNA gene expression. We examined the global gene expression similarities and differences between the historically significant Xenopus laevis model system and the increasingly used X.tropicalis model system and assessed whether an X.tropicalis microarray platform can be used for X.laevis. These closely related species were also used to investigate a more general question: is there an association between mRNA sequence divergence and differences in gene expression levels? We carried out a comprehensive comparison of global gene expression profiles using microarrays of different tissues and developmental stages of X.laevis and X.tropicalis. We (i) show that the X.tropicalis probes provide an efficacious microarray platform for X.laevis, (ii) describe methods to compare interspecies mRNA profiles that correct differences in hybridization efficiency and (iii) show independently of hybridization bias that as mRNA sequence divergence increases between X.laevis and X.tropicalis differences in mRNA expression levels also increase.

Animals↗

Developmental characteristics of adapting mouse small intestine crypt cells.

BACKGROUND & AIMS: Following massive small bowel resection (SBR), the remnant intestine undergoes an adaptive process characterized by increases in a number of physiologic and morphologic parameters. These changes are the result of a stimulus that increases crypt cell mitosis and augments cellular progression along the villus axis. To better define this process, we identified patterns of gene expression specifically within adapting intestinal crypt cells following SBR. METHODS: Laser capture microdissection was used to isolate mouse intestinal crypt cells following SBR or sham operation. Multiple biological and technical complementary DNA microarray replicates allowed rigorous statistical analyses for identification of important expression profiles. Major groups of genes were classified as to site of action, functional pathway, and possible regulatory groups. RESULTS: A total of 300 genes differentially expressed at significant levels within adapting crypt enterocytes were analyzed. Comparison of this list of differentially expressed adapting crypt cell genes with a generalized mouse gene expression database (from 82 developing and adult mouse tissues) showed the greatest overlap with developing and immature intestinal tissues. We identified prominent groups of genes involved with cell growth, signal transduction, and nucleic acid binding. Genes not previously shown to be involved with adaptation or development and maturation were identified. CONCLUSIONS: Identification of similar genes coordinately regulated during both adaptation and development, processes that share key morphologic features, provides a basis for new mechanistic insights into these shared characteristics.

Adaptation, Physiological↗

The role of metallothionein in the pathogenesis of acute lung injury.

Often fatal, acute lung injury has a complicated etiology. Previous studies from our laboratory in mice have demonstrated that survival during acute lung injury is a complex trait governed by multiple loci. We also found that the increase in metallothionein (MT) is one of the greatest noted in transcriptome-wide analyses of gene expression. To assess the role of MT in nickel-induced acute lung injury, the survival of Mt-transgenic, Mt1/2(+/+), and Mt1/2(-/-) mice was compared. Pulmonary inflammation and global gene expression were compared in Mt1/2(+/+) and Mt1/2(-/-) mice. Gene-targeted Mt1/2(-/-) mice were more susceptible than Mt1/2(+/+) mice to nickel-induced inflammation, surfactant-associated protein B transcript loss, and lethality. Similarly, Mt-transgenic mice exhibited increased survival. MAPPFinder analyses also noted significant decreases in genes involved in protein processing (e.g., ubiquitination, folding), which were greater in Mt1/2(-/-) mice as compared with Mt1/2(+/+) mice early in the progression of acute lung injury, possibly due to a zinc-mediated transcript destabilization. In contrast, transcript levels of genes associated with the inflammatory response, extracellular matrix regulation, and coagulation/fibrinolysis were increased more in Mt1/2(-/-) mice as compared with Mt1/2(+/+) mice late in the development of acute lung injury. Thus, MT ultimately improves survival in the progression of acute lung injury in mice. Transcriptome-wide analysis suggests that this survival may be mediated through changes in the destabilization of transcripts associated with protein processing, the subsequent augmentation of transcripts controlling inflammation, extracellular matrix regulation, coagulation/fibrinolysis, and disruption of surfactant homeostasis.

Animals↗

Gene expression changes during the development of acute lung injury: role of transforming growth factor beta.

RATIONALE: Acute lung injury can occur from multiple causes, resulting in high mortality. The pathophysiology of nickel-induced acute lung injury in mice is remarkably complex, and the molecular mechanisms are uncertain. OBJECTIVES: To integrate molecular pathways and investigate the role of transforming growth factor beta (TGF-beta) in acute lung injury in mice. METHODS: cDNA microarray analyses were used to identify lung gene expression changes after nickel exposure. MAPPFinder analysis of the microarray data was used to determine significantly altered molecular pathways. TGF-beta1 protein in bronchoalveolar lavage fluid, as well as the effect of inhibition of TGF-beta, was assessed in nickel-exposed mice. The effect of TGF-beta on surfactant-associated protein B (Sftpb) promoter activity was measured in mouse lung epithelial cells. MEASUREMENTS AND MAIN RESULTS: Genes that decreased the most after nickel exposure play important roles in lung fluid absorption or surfactant and phospholipid synthesis, and genes that increased the most were involved in TGF-beta signaling. MAPPFinder analysis further established TGF-beta signaling to be significantly altered. TGF-beta-inducible genes involved in the regulation of extracellular matrix function and fibrinolysis were significantly increased after nickel exposure, and TGF-beta1 protein was also increased in the lavage fluid. Pharmacologic inhibition of TGF-beta attenuated nickel-induced protein in bronchoalveolar lavage. In addition, treatment with TGF-beta1 dose-dependently repressed Sftpb promoter activity in vitro, and a novel TGF-beta-responsive region in the Sftpb promoter was identified. CONCLUSIONS: These data suggest that TGF-beta acts as a central mediator of acute lung injury through the alteration of several different molecular pathways.

Animals↗

Competitive coexistence of two Pneumocystis species.

Pneumocystis are fungal pathogens of mammalian lungs that can cause lethal pneumonia in immunocompromised hosts. In some mammals, coinfections of genetically distinct Pneumocystis populations have been identified, but the nature of their interaction and its significance are unknown. Two species that infect rats, Pneumocystis carinii and Pneumocystis wakefieldiae, were studied over a 6-year period, representing approximately 700 generations of Pneumocystis. Population densities of each species were analyzed within the framework of the Lotka-Volterra competition model, which revealed the two species were in competition and predicted competitive exclusion of one species. However, stable coexistence was observed in 460 replicate populations. Selected extrinsic factors that might mitigate the extinction were evaluated. Logistic-regression analyses showed that higher relative humidity and higher organism lung burdens were associated with infections comprised of P. carinii alone, while lower temperatures and an increased rat census were associated with the presence of P. wakefieldiae. PCR and immunofluorescent analysis of rat lung tissue showed that both species were present within the same alveoli, excluding habitat heterogeneity as a mechanism of coexistence. These data suggest that P. carinii and P. wakefieldiae were in competitive coexistence, which was influenced in part by extrinsic factors. To our knowledge, this is the first report to evaluate interactions of pathogenic fungal species within a mammalian host using ecological models.

Acquired Immunodeficiency Syndrome↗

Characterization of novel non-clonal intrachromosomal rearrangements between the H4 and PTEN genes (H4/PTEN) in human thyroid cell lines and papillary thyroid cancer specimens.

The two main forms of RET rearrangement in papillary thyroid carcinomas (PTC) arise from intrachromosomal inversions fusing the tyrosine kinase domain of RET with either the H4 (RET/PTC1) or the ELE1/RFG genes (RET/PTC3). PTEN codes for a dual-specificity phosphatase and maps to chromosome 10q22-23. Germline mutations confer susceptibility to Cowden syndrome whereas somatic mutations or deletions are common in several sporadic human tumors. Decreased PTEN expression has been implicated in thyroid cancer development. We report the characterization of a new chromosome 10 rearrangement involving H4 and PTEN. The initial H4/PTEN rearrangement was discovered as a non-specific product of RT-PCR for RET/PTC1 in irradiated thyroid cell lines. Sequencing revealed a transcript consisting of exon 1 and 2 of H4 fused with exons 3-6 of PTEN. Nested RT-PCR with specific primers bracketing the breakpoints confirmed the H4/PTEN rearrangements in irradiated KAT-1 and KAT-50 cells. Additional H4/PTEN variants, generated by recombination of either exon 1 or exon 2 of H4 with exon 6 of PTEN, were found in non-irradiated KAK-1, KAT-50, ARO and NPA cells. Their origin through chromosomal recombination was confirmed by detection of the reciprocal PTEN/H4 product. H4/PTEN recombination was not a clonal event in any of the cell lines, as Southern blots with appropriate probes failed to demonstrate aberrant bands, and multicolor FISH of KAK1 cells with BAC probes for H4 and PTEN did not show a signal overlap in all cells. Based on PCR of serially diluted samples, the minimal frequency of spontaneous recombination between these loci was estimated to be approximately 1/10(6) cells. H4/PTEN products were found by nested RT-PCR in 4/14 normal thyroid tissues (28%) and 14/18 PTC (78%) (P<0.01). H4/PTEN is another example of recombination involving the H4 locus, and points to the high susceptibility of thyroid cells to intrachromosomal gene rearrangements. As this also represents a plausible mechanism for loss-of-function of PTEN, other thyroid neoplastic phenotypes and eventually other cancer types need to be screened for clonal H4/PTEN rearrangements.

Chromosome Mapping↗

Modeling variation in tumors in vivo.

Transgenic mice that allow mutant cells to be visualized in situ were used to study variation in tumors. These mice carry the G11 placental alkaline phosphatase (PLAP) transgene, a mutant allele rendered incapable of producing its enzyme product by a frameshift caused by insertion of a tract of G:C base pairs in a coding region. Spontaneous deletion of one G:C base pair from this tract restores gene function, and cells with PLAP activity can be detected histochemically. To study tumors, the G11 PLAP transgene was introduced into the polyoma virus middle T antigen mammary tumor model. Tumors in these mice exhibited up to 300 times more PLAP+ cells than normal tissues. PLAP+ cells were located throughout each tumor. Many of the PLAP+ cells were singlets, but clusters also were common, with one cluster containing >30,000 cells. Comparison of these data to simulations produced by computer models suggested that multiple factors were involved in generating mutant cells in tumors. Although genetic instability appeared to have occurred in most tumors, large clusters were much more common than expected based on instability alone.

Alkaline Phosphatase↗

Neural system-enriched gene expression: relationship to biological pathways and neurological diseases.

To understand the commitment of the genome to nervous system differentiation and function, we sought to compare nervous system gene expression to that of a wide variety of other tissues by gene expression database construction and mining. Gene expression profiles of 10 different adult nervous tissues were compared with that of 72 other tissues. Using ANOVA, we identified 1,361 genes whose expression was higher in the nervous system than other organs and, separately, 600 genes whose expression was at least threefold higher in one or more regions of the nervous system compared with their median expression across all organs. Of the 600 genes, 381 overlapped with the 1,361-gene list. Limited in situ gene expression analysis confirmed that identified genes did represent nervous system-enriched gene expression, and we therefore sought to evaluate the validity and significance of these top-ranked nervous system genes using known gene literature and gene ontology categorization criteria. Diverse functional categories were present in the 381 genes, including genes involved in intracellular signaling, cytoskeleton structure and function, enzymes, RNA metabolism and transcription, membrane proteins, as well as cell differentiation, death, proliferation, and division. We searched existing public sites and identified 110 known genes related to mental retardation, neurological disease, and neurodegeneration. Twenty-one of the 381 genes were within the 110-gene list, compared with a random expectation of 5. This suggests that the 381 genes provide a candidate set for further analyses in neurological and psychiatric disease studies and that as a field, we are as yet, far from a large-scale understanding of the genes that are critical for nervous system structure and function. Together, our data indicate the power of profiling an individual biologic system in a multisystem context to gain insight into the genomic basis of its structure and function.

Animals↗

From co-expression to co-regulation: how many microarray experiments do we need?

BACKGROUND: Cluster analysis is often used to infer regulatory modules or biological function by associating unknown genes with other genes that have similar expression patterns and known regulatory elements or functions. However, clustering results may not have any biological relevance. RESULTS: We applied various clustering algorithms to microarray datasets with different sizes, and we evaluated the clustering results by determining the fraction of gene pairs from the same clusters that share at least one known common transcription factor. We used both yeast transcription factor databases (SCPD, YPD) and chromatin immunoprecipitation (ChIP) data to evaluate our clustering results. We showed that the ability to identify co-regulated genes from clustering results is strongly dependent on the number of microarray experiments used in cluster analysis and the accuracy of these associations plateaus at between 50 and 100 experiments on yeast data. Moreover, the model-based clustering algorithm MCLUST consistently outperforms more traditional methods in accurately assigning co-regulated genes to the same clusters on standardized data. CONCLUSIONS: Our results are consistent with respect to independent evaluation criteria that strengthen our confidence in our results. However, when one compares ChIP data to YPD, the false-negative rate is approximately 80% using the recommended p-value of 0.001. In addition, we showed that even with large numbers of experiments, the false-positive rate may exceed the true-positive rate. In particular, even when all experiments are included, the best results produce clusters with only a 28% true-positive rate using known gene transcription factor interactions.

Algorithms↗

Mutational biases associated with potential iron-binding DNA motifs in rodent lacI and human p53 mutational databases.

The role of Fenton oxidants in DNA damage, aging, and cancer is appreciated, but not well understood. Six potential iron-binding (PIB) DNA motifs were previously identified as sites of preferential strand cleavage. Since DNA-metal binding domains are a known determinant of oxidative DNA damage, and the location of strand breaks explains where oxidant attack occurs, we sought to determine whether the likelihood of base change mutations is a function of neighboring PIB motifs. We developed a sliding window function that computes the density of PIB motifs on both strands, within 4-12bp, for each location along a target gene. This range of window sizes reflects known diffusion distances of Fenton reaction products. Using mutational databases, odds of mutation at each base were calculated relative to PIB motif density, for all PIB motif types in aggregate, or for individual PIB motifs. Using mutational data from lacI transgenic animals, we observed a non-random distribution of PIB motifs, associated with increased odds of mutation, showing a strand bias. Sensitivity analysis confirmed that the optimum association between PIB motif density and mutations occurs when a 7bp radius is used for the window size. Randomly simulated mutations showed no association with PIB motif density. When the method was applied to human TP53 mutation data, we saw similar results, but no strand bias. As PIB motif density rises, linear trends are observed for increasing odds of mutation. Sensitivity analysis revealed associations between PIB motifs and GC --> AT transitions and GC --> TA transversions-the most commonly observed types of mutations arising from oxidative DNA damage. DNA-metal binding motifs are found in a wide variety of biological contexts, including many where conformational sensitivity to redox state is important. These techniques can help elucidate how DNA-iron-binding may affect lesions and subsequent mutations from multiple agents.

Amino Acid Motifs↗

Expression of genes in the TGF-beta signaling pathway is significantly deregulated in smooth muscle cells from aorta of aryl hydrocarbon receptor knockout mice.

The molecular basis for the adverse biological effects of dioxin (2,3,7,8-tetrachlorodibenzo-p-dioxin; TCDD), a pervasive environmental toxin, is largely unknown. TCDD is a ligand for the cytosolic aromatic hydrocarbon receptor (AHR) which mediates the transcriptional induction of the xenobiotic metabolizing genes in the CYP1 family of cytochromes P450. Previous studies have suggested that the AHR may carry out important functions in the cell in addition to metabolizing toxins. We present gene expression profiles of smooth muscle cells from wild type and Ahr(-/-) mice that show significant changes in the RNA levels of the transforming growth factor-beta3 (Tgfb3) gene and genes involved in the modulation and processing of TGF-beta. The RNA expression profiles support a hypothesis that in the wild type, the AHR represses Tgfb gene expression and affects the gene expression of several TGF-beta-modulating and processing genes. We also observed that RNA levels increased for TGF-beta2, CYP1b1, and TGF-beta-related genes in Ahr(-/-) smooth muscle cells exposed to TCDD. These data are consistent with a hypothesis that TCDD stimulates the TGF-beta2 signaling pathway in the absence of the AHR to activate the Cyp1b1 gene. The above results provide a possible explanation for some of the multiple biological effects of TCDD and the physiological role played by the AHR in the absence of environmental agents.

Animals↗

Critical regulation of genes for tumor cell migration by AP-1.

The AP-1 transcription factor plays a critical role in regulating tumor cell proliferation and has been implicated in controlling a program of gene expression that mediates cell motility and invasion in vitro. We have utilized two dominant negative AP-1 constructs, TAM67 and aFos, each fused to GFP, to investigate the role of AP-1 complexes in an invasive, clinically derived human tumor cell line, HT-1080. As expected, high levels of both GFP-TAM67 and GFP-aFos arrested HT-1080 cells in the G1 phase of the cell cycle. Strikingly, at low levels GFP-aFos, but not GFP-TAM67, caused a change in colony morphology, impairment of directional motility in a monolayer wound healing assay, as well as inhibition of chemotaxis and haptotaxis. Microarray analysis identified a novel set of AP-1 target genes, including the tumor suppressor TSCL-1 and regulators of actin cytoskeletal dynamics, including the gelsolin-like actin capping protein CapG. The demonstration that AP-1 regulates the expression of genes involved in tumor cell motility and cytoskeletal dynamics in a clinically derived human tumor cell line identifies new pathways of control for tumor cell motility.

Cell Adhesion Molecule-1↗

Different global gene expression profiles in benzo[a]pyrene- and dioxin-treated vascular smooth muscle cells of AHR-knockout and wild-type mice.

Benzo[a]pyrene (B[a]P) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) are potent ligands for the aryl hydrocarbon receptor (AHR). High-density oligonucleotide microarrays were used to generate global gene expression profiles of wild-type and Ahr(-/-) vascular smooth muscle cells (SMCs) from mouse aorta. To determine whether there are signaling pathways other than the AHR involved in B[a]P metabolism, wild-type and AHR knockout (Ahr(-/-) SMCs were exposed to B[a]P. Two signaling pathways, represented by TGF-beta2 and IGF-1, were identified as potential candidates of an AHR alternate pathway for cells to respond to B[a]P. The wild-type SMCs responded similarly to B[a]P and TCDD in the regulation of a small set of common genes known to respond to the activated AHR (e.g., glutamine S-transferase). However, wild-type SMCs responded in a way that involves many additional genes, suggesting that a very divergent cellular response may be involved when SMCs are exposed to the two classic inducers of the AHR. In contrast, many more genes in the Ahr(-/-) cells responded similarly to B[a]P and TCDD, including Cyp1b1, than responded differently, which indicates that eliminating the AHR is effective for investigating potential alternate cellular mechanisms that respond to B[a]P and TCDD.

Animals↗

Gene expression profiles of mouse aorta and cultured vascular smooth muscle cells differ widely, yet show common responses to dioxin exposure.

Exposure to environmental toxicants may play a role in the onset and progression of cardiovascular disease. Many environmental agents, such as dioxin, are risk factors for atherosclerosis because they may exacerbate an underlying disease by altering gene expression patterns. Expression profiling of vascular tissues allows the simultaneous analysis of thousands of genes and may provide predictive information particularly useful in early disease stages. Often, however, in vivo experiments are unfeasible for material or ethical reasons, and data from cultured cells must be used instead, even though it may not be known whether cultured cells and live tissues share common global responses to the same toxicant. In a search for genes responsive to dioxin exposure, we used oligonucleotide microarrays with DNA sequences from 13,433 genes to compare global gene expression profiles of C57BL/6 mice aortas with cultured vascular smooth muscle cells (vSMCs) of the same mice. Aorta segments and vSMCs differed in the expression of more than 4500 genes, many showing expression differences greater than 1000-fold. Integration of microarray data into Gene Ontology Project annotations showed that many of the genes differentially expressed belonged to the same biological process or metabolic pathway. Notwithstanding these results, a subset of 35 genes responded in the same fashion to dioxin exposure in both systems. Genes in this subset encoded phase I and phase II detoxification enzymes, signal transduction kinases and phosphatases, and proteins involved in DNA repair and the cell cycle. We conclude that vSMCS may be useful aorta surrogates to study early gene expression responses to dioxin exposure, provided that analyses focus on this subset of genes.

Animals↗

Microarray results improve significantly as hybridization approaches equilibrium.

Dual-channel long oligonucleotide microarrays are in widespread use. Although much attention has been given to proper experimental design and analysis regarding long oligonucleotide microarrays, relatively little information is available concerning the optimization of protocols. We carried out a series of microarray experiments designed to investigate the effects of different levels of target concentration and hybridization times using a long oligonucleotide library. Based on principles developed from nucleic acid renaturation kinetics studies, we show that increasing the time of hybridization from 18 h to 42 h and 66 h, especially when lower than optimal concentrations of target were used, significantly improved the quality of the microarray results. Longer hybridization times significantly increased the number of spots detected, signal-to-noise ratios, and the number of differentially expressed genes and correlations among replicate arrays. We conclude that at 18 h of incubation, target-to-probe hybridization has not reached equilibrium and that a relatively high proportion of nonspecific hybridization occurs. This result is striking, given that most, if not all, published microarray protocols stipulate 8-24 h for hybridization. Using shorter than optimal hybridization times (i.e., not allowing hybridization to reach equilibrium) has the consequence of underestimating the fold change of differentially expressed genes and of missing less represented sequences.

Algorithms↗