Characteristics of oligonucleotide tiling arrays measured by hybridizing full-length cDNA clones: causes of signal variation and false positive signals

An assessment of the hybridization characteristics of oligonucleotide tiling arrays was carried out using 162 full-length sequenced cDNA clones in spike-in experiments. The properties of array probes that influence signal intensity were investigated, and their capability in the detection of the cDNA exons was evaluated. The signal intensities detected in exonic and nonexonic genomic regions were examined by focusing on the features of probe sequences that raise or lower the level of intensity and on the causes of false positive signals found in nonexonic regions. The effectiveness of measures used in published protocols to improve the separation between signal and background intensity distributions, including the use of replicates and threshold parameterization of signal intensity, was assessed. Sensitivity and specificity in the detection of exons were measured using various sets of threshold parameters, and the effects of each parameter on the detection efficiency and the rate of false positives were evaluated. It was also demonstrated that hybridization of full-length cDNA clones is an excellent method to investigate the characteristics of oligonucleotide tiling arrays.     

Comparison of T Cell Receptor-Induced Proximal Signaling and Downstream Functions in Immortalized and Primary T Cells

Background

Human T cells play an important role in pathogen clearance, but their aberrant activation is also linked to numerous diseases. T cells are activated by the concurrent induction of the T cell receptor (TCR) and one or more costimulatory receptors. The characterization of signaling pathways induced by TCR and/or costimulatory receptor activation is critical, since these pathways are excellent targets for novel therapies for human disease. Although studies using human T cell lines have provided substantial insight into these signaling pathways, no comprehensive, direct comparison of these cell lines to activated peripheral blood T cells (APBTs) has been performed to validate their usefulness as a model of primary T cells.

Methodology/Principal Findings

We used quantitative biochemical techniques to compare the activation of two widely used human T cell lines, Jurkat E6.1 and HuT78 T cells, to APBTs. We found that HuT78 cells were similar to APBTs in proximal TCR-mediated signaling events. In contrast, Jurkat E6.1 cells had significantly increased site-specific phosphorylation of Pyk2, PLCγ1, Vav1, and Erk1/Erk2 and substantially more Ca²⁺ flux compared to HuT78 cells and APBTs. In part, these effects appear to be due to an overexpression of Itk in Jurkat E6.1 cells compared to HuT78 cells and APBTs. Both cell lines differ from APBTs in the expression and function of costimulatory receptors and in the range of cytokines and chemokines released upon TCR and costimulatory receptor activation.

Conclusions/Significance

Both Jurkat E6.1 and HuT78 T cells had distinct similarities and differences compared to APBTs. Both cell lines have advantages and disadvantages, which must be taken into account when choosing them as a model T cell line.     

Rosetta error model for gene expression analysis

MOTIVATION: In microarray gene expression studies, the number of replicated microarrays is usually small because of cost and sample availability, resulting in unreliable variance estimation and thus unreliable statistical hypothesis tests. The unreliable variance estimation is further complicated by the fact that the technology-specific variance is intrinsically intensity-dependent. RESULTS: The Rosetta error model captures the variance-intensity relationship for various types of microarray technologies, such as single-color arrays and two-color arrays. This error model conservatively estimates intensity error and uses this value to stabilize the variance estimation. We present two commonly used error models: the intensity error-model for single-color microarrays and the ratio error model for two-color microarrays or ratios built from two single-color arrays. We present examples to demonstrate the strength of our error models in improving statistical power of microarray data analysis, particularly, in increasing expression detection sensitivity and specificity when the number of replicates is limited.     

Indirect measurements of differential gene expression with cDNA microarrays

The use of universal RNA reference sets is an increasingly common approach to molecular classification studies with cDNA microarrays. Here we evaluated the reliability of indirect measurements of fluorescence ratios with a common RNA reference as a means of identifying differentially expressed genes. Comparisons of direct and indirect measures of differential gene expression showed a strong overall correlation in fluorescence ratio measurements but also a high degree of false positives in our indirect measurements. These results indicated that the application of more stringent ratio filters may be required when assessing differential gene expression utilizing a common RNA reference in classification studies.     

Microarrays for bacterial detection and microbial community analysis

Several types of microarrays have recently been developed and evaluated for bacterial detection and microbial community analysis. These studies demonstrated that specific, sensitive and quantitative detection could be obtained with both functional gene arrays and community genome arrays. Although single-base mismatch can be differentiated with phylogenetic oligonucleotide arrays, reliable specific detection at the single-base level is still problematic. Microarray-based hybridization approaches are also useful for defining genome diversity and bacterial relatedness. However, more rigorous and systematic assessment and development are needed to realize the full potential of microarrays for microbial detection and community analysis.     

Exogenous expression of SAMHD1 inhibits proliferation and induces apoptosis in cutaneous T-cell lymphoma-derived HuT78 cells

Sterile α motif and HD domain-containing protein 1 (SAMHD1) is a mammalian dNTP hydrolase (dNTPase) that regulates intracellular dNTP balance. We have previously reported that SAMHD1 mRNA and protein levels are significantly downregulated in CD4⁺ T-cells of patients with cutaneous T-cell lymphoma (CTCL), a disease characterized by infiltration of neoplastic CD4⁺ T-lymphocytes into the skin. However, functional significance of SAMHD1 in CTCL development and progression remains unknown. Here we investigate the mechanism by which SAMHD1 induces apoptosis in CTCL-derived CD4⁺ T-cells. We stably expressed exogenous SAMHD1 in the CTCL-derived HuT78 T-cell line containing a very low level of endogenous SAMHD1 protein. We found that low-level exogenous expression of SAMHD1 led to a significant reduction in HuT78 cell growth, proliferation, and colony formation. Exogenous SAMHD1 expression in HuT78 cells also resulted in increased spontaneous and Fas ligand (Fas-L)-induced apoptosis levels via activation of the extrinsic pathway, including caspase-8, ?3 and ?7. Additionally, increased SAMHD1 significantly reduced the protein and mRNA expression of the short isoform of cFLIP (cFLIP_S), an important negative regulator of Fas-L-mediated apoptotic signaling. Our results indicate that exogenous SAMHD1 expression inhibits HuT78 cell growth and proliferation in part by increasing apoptosis. These findings implicate that SAMHD1 acts as an inhibitor in CTCL cell growth, suggesting that downregulation of SAMHD1 expression in neoplastic T-cells can facilitate uncontrolled cell proliferation.     

A quality-controlled microarray method for gene expression profiling

Gene expression profiling on microarrays is widely used to measure the expression of large numbers of genes in a single experiment. Because of the high cost of this method, feasible numbers of replicates are limited, thus impairing the power of statistical analysis. As a step toward reducing technically induced variation, we developed a procedure of sample preparation and analysis that minimizes the number of sample manipulation steps, introduces quality control before array hybridization, and allows recovery of the prepared mRNA for independent validation of results. Sample preparation is based on mRNA separation using oligo(dT) magnetic beads, which are subsequently used for first-strand cDNA synthesis on the beads. cDNA covalently bound to the magnetic beads is used as template for second-strand cDNA synthesis, leaving the intact mRNA in solution for further analysis. The quality of the synthesized cDNA can be assessed by quantitative polymerase chain reaction using 3'- and 5'-specific primer pairs for housekeeping genes such as glyceraldehyde-3-phosphate dehydrogenase. Second-strand cDNA is chemically labeled with fluorescent dyes to avoid dye bias in enzymatic labeling reactions. After hybridization of two differently labeled samples to microarray slides, arrays are scanned and images analyzed automatically with high reproducibility. Quantile-normalized data from five biological replica display a coefficient of variation 45% for 90% of profiled genes, allowing detection of twofold changes with false positive and false negative rates of 10% each. We demonstrate successful application of the procedure for expression profiling in plant leaf tissue. However, the method could be easily adapted for samples from animal including human or from microbial origin.     

Postinternalization inhibition of adenovirus gene expression and infectious virus production in human T-cell lines

Detection of adenovirus DNA in human tonsillar T cells in the absence of active virus replication suggests that T cells may be a site of latency or of attenuated virus replication in persistently infected individuals. The lytic replication cycle of Ad5 in permissive epithelial cells (A549) was compared to the behavior of Ad5 in four human T-cell lines, Jurkat, HuT78, CEM, and KE37. All four T-cell lines expressed the integrin coreceptors for Ad2 and Ad5, but only Jurkat and HuT78 express detectable surface levels of the coxsackie adenovirus receptor (CAR). Jurkat and HuT78 cells supported full lytic replication of Ad5, albeit at a level approximately 10% of that of A549, while CAR-transduced CEM and KE37 cells (CEM-CARhi and KE37-CARhi, respectively) produced no detectable virus following infection. All four T-cell lines bind and internalize fluorescently labeled virus. In A549, Jurkat, and HuT78 cells, viral proteins were detected in 95% of cells. In contrast, only a small subpopulation of CEM-CARhi and KE37-CARhi cells contained detectable viral proteins. Interestingly, Jurkat and HuT78 cells synthesize four to six times more copies of viral DNA per cell than did A549 cells, indicating that these cells produce infectious virions with much lower efficiency than A549. Similarly, CEM-CARhi and KE37-CARhi cells, which produce no detectable infectious virus, synthesize three times more viral genomes per cell than A549. The observed blocks to adenovirus gene expression and replication in all four human T-cell lines may contribute to the maintenance of naturally occurring persistent adenovirus infections in human T cells.     

Reduction of autofluorescence on DNA microarrays and slide surfaces by treatment with sodium borohydride

Microarray technology is currently being used extensively in functional genomics research and modern drug discovery and development. Henceforward, tremendous application potential for this technology exists in the fields of clinical diagnostics and prognostics, pathology, and toxicology for high-throughput analysis of "disease" gene expression. However, the major hurdle now in this technology is not the performance of the arrays but rather the efficient reproducibility of the hybridization signal intensity in a fluorescence-based analysis. The sensitivity of fluorescence detection on an array is to a large extent limited by the amount of background signal arising due to nonspecifically bound probes and fluorescence that is intrinsically associated with the chip substrate and/or the attached target DNA, the so-called autofluorescence. Here, we describe a simple and efficient method to reduce autofluorescence from undetermined sources on coated glass slides with and without DNA arrays. This sodium borohydride-mediated reduction process resulted in significantly lower and more even background fluorescence. This in turn extended the dynamic range of detection and reduced the average coefficient of variation of fluorescent signal ratios on DNA microarrays in addition to improving the detection of genes that are expressed at a low level.     

Kinetic studies and unstructured models of lymphocyte metabolism in fed-batch culture

The growth of two lymphocyte cell lines, a hybridoma cell line and a human cutaneous T cell lymphoma (HuT78), was studied in fed-batch culture, and unstructured models of growth developed. A criteria was established to insure that the growth rate varied by less than a specified tolerance throughout the culture period. Glutamine and serum were growth-limiting nutrients for both cell lines with half-maximal growth rates at 0. 53 mM glutamine and 0. 55%(v/v) serum for the hybridoma cells and 0. 21 mM glutamine and 1. 5% serum for the HuT-78 cells. Over the range of glucose concentrations from 5. 5 mM to 28 mM, the specific growth rate of hybridoma cells was independent of glucose concentration, whereas glucose concentrations above 5. 5 mM inhibited HuT-78 growth. For both cell lines, the growth rate was significantly inhibited by the addition of ammonium, although the hybridoma cell line was more affected by ammonia than was the HuT-78 cell line. Growth of HuT-78 cells increased in the presence of interleukin-2. Unstructured models for the hybridoma cells were similar to other models presented in the literature. Applications of these models to adoptive immunotherapy are discussed.     

Effects of filtering by Present call on analysis of microarray experiments

Background

Affymetrix GeneChips^® are widely used for expression profiling of tens of thousands of genes. The large number of comparisons can lead to false positives. Various methods have been used to reduce false positives, but they have rarely been compared or quantitatively evaluated. Here we describe and evaluate a simple method that uses the detection (Present/Absent) call generated by the Affymetrix microarray suite version 5 software (MAS5) to remove data that is not reliably detected before further analysis, and compare this with filtering by expression level. We explore the effects of various thresholds for removing data in experiments of different size (from 3 to 10 arrays per treatment), as well as their relative power to detect significant differences in expression.

Results

Our approach sets a threshold for the fraction of arrays called Present in at least one treatment group. This method removes a large percentage of probe sets called Absent before carrying out the comparisons, while retaining most of the probe sets called Present. It preferentially retains the more significant probe sets (p ≤ 0.001) and those probe sets that are turned on or off, and improves the false discovery rate. Permutations to estimate false positives indicate that probe sets removed by the filter contribute a disproportionate number of false positives. Filtering by fraction Present is effective when applied to data generated either by the MAS5 algorithm or by other probe-level algorithms, for example RMA (robust multichip average). Experiment size greatly affects the ability to reproducibly detect significant differences, and also impacts the effect of filtering; smaller experiments (3–5 samples per treatment group) benefit from more restrictive filtering (≥50% Present).

Conclusion

Use of a threshold fraction of Present detection calls (derived by MAS5) provided a simple method that effectively eliminated from analysis probe sets that are unlikely to be reliable while preserving the most significant probe sets and those turned on or off; it thereby increased the ratio of true positives to false positives.     

Development and characterization of a protein kinase C beta-isozyme-deficient T-cell line.

In the human T-cell lymphoma line, HuT 78, proliferation and phorbol ester-induced growth arrest and differentiation were inhibited by the protein kinase C (PKC) inhibitor, staurosporine. By contrast, an alternative PKC inhibitor, H-7, inhibited proliferation but not phorbol ester-induced growth arrest. The cell line was found to contain both alpha and beta isoforms of PKC by Western blot techniques. A cell line, K-4, was cloned from HuT 78 in the presence of H-7 and this clone was found to be positive for PKC-alpha only. PKC-beta did not return on cultivation in the absence of H-7. Proliferation of K-4 was insensitive to inhibition with both H-7 and staurosporine. However, phorbol ester-induced growth arrest remained staurosporine sensitive. Phorbol-stimulated IL-2 secretion was minimal in the PKC-beta-deficient cell line. These data suggest that PKC-beta may be a regulatory enzyme for proliferation and stimulated interleukin-2 secretion in HuT 78 cells. Heterogeneity of responses to PKC activation may reflect the use of different isozymes in different intracellular pathways. The K-4 cell line should provide a useful tool in the dissection of involvement of PKC isozymes in cellular function.     

Interleukin-2 induces the activities of DNA topoisomerase I and DNA topoisomerase II in HuT 78 cells

The induction by interleukin-2 of DNA topoisomerase I and DNA topoisomerase II activities in the human T cell line HuT 78 was investigated. HuT 78 cells were treated with 1000 U of interleukin-2/ml, and extracts of the HuT 78 nuclei were prepared over a 24 h period. The extracts were assayed quantitatively for the activities of DNA topoisomerase I and DNA topoisomerase II. Three concomitant, transient increases of 3- to 11-fold in the specific activities of both DNA topoisomerase I and DNA topoisomerase II were observed following treatment with IL-2 at 0.5, 4, and 10 h after treatment with interleukin-2. The specific activities of both enzymes returned to base-line values after each of these transient increases. These results reveal that the activities of DNA topoisomerase I and DNA topoisomerase II are highly regulated in HuT 78 cells upon treatment with IL-2.     

Data rotation improves genomotyping efficiency

Unsequenced bacterial strains can be characterized by comparing their genomic DNA to a sequenced reference genome of the same species. This comparative genomic approach, also called genomotyping, is leading to an increased understanding of bacterial evolution and pathogenesis. It is efficiently accomplished by comparative genomic hybridization on custom-designed cDNA microarrays. The microarray experiment results in fluorescence intensities for reference and sample genome for each gene. The log-ratio of these intensities is usually compared to a cut-off, classifying each gene of the sample genome as a candidate for an absent or present gene with respect to the reference genome. Reducing the usually high rate of false positives in the list of candidates for absent genes is decisive for both time and costs of the experiment. We propose a novel method to improve efficiency of genomotyping experiments in this sense, by rotating the normalized intensity data before setting up the list of candidate genes. We analyze simulated genomotyping data and also re-analyze an experimental data set for comparison and illustration. We approximately halve the proportion of false positives in the list of candidate absent genes for the example comparative genomic hybridization experiment as well as for the simulation experiments.     

Dissecting systems-wide data using mixture models: application to identify affected cellular processes

Background  

Functional analysis of data from genome-scale experiments, such as microarrays, requires an extensive selection of differentially expressed genes. Under many conditions, the proportion of differentially expressed genes is considerable, making the selection criteria a balance between the inclusion of false positives and the exclusion of false negatives.