Identification and handling of artifactual gene expression profiles emerging in microarray hybridization experiments

Identification of artifactual microarray probe signals constantly present in multiple sample types

The detection, identification, and quantitation of transcripts have evolved from simple Northern analysis, cDNA cloning, and sequencing to RT-PCR, microarrays, and now digital gene expression using ultra-high-throughput RNA sequencing (RNA-Seq). During the course of our studies we observed that some microarray probes show very high signal intensity values yet are discordant when compared with RNA-Seq. A total of 99 probes from approximately 30,000 were identified as consistently discordant in four human tissues or cell lines. Interestingly, this set of discordant probes appears array-dependent. Among the 99 probes identified, 70 constantly exhibited a high signal in all 713 available samples surveyed using the Illumina HumanHT-12v4 platform. Some were discordant with additional probes that annotated the same genes. Absence of a number of these transcripts was confirmed by quantitative RT-PCR (qRT-PCR). Our findings suggest that one must be cautious, as some array probes do not capture the level of the target.

chip artifact CORRECTion (caCORRECT): a bioinformatics system for quality assurance of genomics and proteomics array data

Quality assurance of high throughput "-omics" data is a major concern for biomedical discovery and translational medicine, and is considered a top priority in bioinformatics and systems biology. Here, we report a web-based bioinformatics tool called caCORRECT for chip artifact detection, analysis, and CORRECTion, which removes systematic artifactual noises that are commonly observed in microarray gene expression data. Despite the development of major databases such as GEO arrayExpress, caArray, and the SMD to manage and distribute microarray data to the public, reproducibility has been questioned in many cases, including high-profile papers and datasets. Based on both archived and synthetic data, we have designed the caCORRECT to have several advanced features: (1) to uncover significant, correctable artifacts that affect reproducibility of experiments; (2) to improve the integrity and quality of public archives by removing artifacts; (3) to provide a universal quality score to aid users in their selection of suitable microarray data; and (4) to improve the true-positive rate of biomarker selection verified by test data. These features are expected to improve the reproducibility of Microarray study. caCORRECT is freely available at: http://caCORRECT.bme.gatech.edu.

Changes in gene expression during pegylated interferon and ribavirin therapy of chronic hepatitis C virus distinguish responders from nonresponders to antiviral therapy

Treating chronic hepatitis C virus (HCV) infection using pegylated alpha interferon and ribavirin leads to sustained clearance of virus and clinical improvement in approximately 50% of patients. Response rates are lower among patients with genotype 1 than with genotypes 2 and 3 and among African-American (AA) patients compared to Caucasian (CA) patients. Using DNA microarrays, gene expression was assessed for a group of 33 African-American and 36 Caucasian American patients with chronic HCV genotype 1 infection during the first 28 days of treatment. Results were examined with respect to treatment responses and to race. Patients showed a response to treatment at the gene expression level in RNA isolated from peripheral blood mononuclear cells irrespective of degree of decrease in HCV RNA levels. However, gene expression responses were relatively blunted in patients with poor viral response (<1.5 log(10)-IU/ml decrease at 28 days) compared to those in patients with a marked (>3.5 log(10)-IU/ml decrease) or intermediate (1.5 to 3.5 log(10)-IU/ml decrease) response. The number of genes that were up- or down-regulated by pegylated interferon and ribavirin treatment was fewer in patients with a poor response than in those with an intermediate or marked viral response. However AA patients had a stronger interferon response than CA patients in general. The induced levels of known interferon-stimulated genes such as the 2'5'-oligoadenylate synthetase, MX1, IRF-7, and toll-like receptor TLR-7 genes was lower in poor-response patients than in marked- or intermediate-response patients. Thus, the relative lack of viral response to interferon therapy of hepatitis C virus infection is associated with blunted interferon cell signaling. No specific regulatory gene could be identified as responsible for this global blunting or the racial differences.

Derivation of species-specific hybridization-like knowledge out of cross-species hybridization results

BACKGROUND

One of the approaches for conducting genomics research in organisms without extant microarray platforms is to profile their expression patterns by using Cross-Species Hybridization (CSH). Several different studies using spotted microarray and CSH produced contradicting conclusions in the ability of CSH to reflect biological processes described by species-specific hybridization (SSH).

RESULTS

We used a tomato-spotted cDNA microarray to examine the ability of CSH to reflect SSH data. Potato RNA was hybridized to spotted cDNA tomato and potato microarrays to generate CSH and SSH data, respectively. Difficulties arose in obtaining transcriptomic data from CSH that reflected those obtained from SSH. Nevertheless, once the data was filtered for those corresponding to matching probe sets, by restricting proper cutoffs of probe homology, the CSH transcriptome data showed improved reflection of those of the SSH.

CONCLUSIONS

This study evaluated the relative performance of CSH compared to SSH, and proposes methods to ensure that CSH closely reflects the biological process analyzed by SSH.

SimArray: a user-friendly and user-configurable microarray design tool

BACKGROUND

Microarrays were first developed to assess gene expression but are now also used to map protein-binding sites and to assess allelic variation between individuals. Regardless of the intended application, efficient production and appropriate array design are key determinants of experimental success. Inefficient production can make larger-scale studies prohibitively expensive, whereas poor array design makes normalisation and data analysis problematic.

RESULTS

We have developed a user-friendly tool, SimArray, which generates a randomised spot layout, computes a maximum meta-grid area, and estimates the print time, in response to user-specified design decisions. Selected parameters include: the number of probes to be printed; the microtitre plate format; the printing pin configuration, and the achievable spot density. SimArray is compatible with all current robotic spotters that employ 96-, 384- or 1536-well microtitre plates, and can be configured to reflect most production environments. Print time and maximum meta-grid area estimates facilitate evaluation of each array design for its suitability. Randomisation of the spot layout facilitates correction of systematic biases by normalisation.

CONCLUSION

SimArray is intended to help both established researchers and those new to the microarray field to develop microarray designs with randomised spot layouts that are compatible with their specific production environment. SimArray is an open-source program and is available from http://www.flychip.org.uk/SimArray/.

Evolutionary regulation of the blind subterranean mole rat, Spalax, revealed by genome-wide gene expression

We applied genome-wide gene expression analysis to the evolutionary processes of adaptive speciation of the Israeli blind subterranean mole rats of the Spalax ehrenbergi superspecies. The four Israeli allospecies climatically and adaptively radiated into the cooler, mesic northern domain (N) and warmer, xeric southern domain (S). The kidney and brain mRNAs of two N and two S animals were examined through cross-species hybridizations with two types of Affymetrix arrays (mouse and rat) and muscle mRNA of six N and six S animals with spotted cDNA mouse arrays. The initial microarray analysis was hypothesis-free, i.e., conducted without reference to the origin of animals. Principal component analysis revealed that 20-30% of the expression signal variability could be explained by the differentiation of N-S species. Similar N-S effects were obtained for all tissues and types of arrays: two Affymetrix microarrays using probe oligomer signals and the spotted array. Likewise, ANOVA and t test statistics demonstrated significant N-S ecogeographic divergence and region-tissue specificity in gene expression. Analysis of differential gene expression between species corroborates previous results deduced by allozymes and DNA molecular polymorphisms. Functional categories show significant N-S ecologic putative adaptive divergent up-regulation of genes highlighting a higher metabolism in N, and potential adaptive brain activity and kidney urine cycle pathways in S. The present results confirm ecologic-genomic separation of blind mole rats into N and S. Gene expression regulation appears to be central to the evolution of blind mole rats.