permGPU: Using graphics processing units in RNA microarray association studies.

BACKGROUND: Many analyses of microarray association studies involve permutation, bootstrap resampling and cross-validation, that are ideally formulated as embarrassingly parallel computing problems. Given that these analyses are computationally intensive, scalable approaches that can take advantage of multi-core processor systems need to be developed. RESULTS: We have developed a CUDA based implementation, permGPU, that employs graphics processing units in microarray association studies. We illustrate the performance and applicability of permGPU within the context of permutation resampling for a number of test statistics. An extensive simulation study demonstrates a dramatic increase in performance when using permGPU on an NVIDIA GTX 280 card compared to an optimized C/C++ solution running on a conventional Linux server. CONCLUSIONS: permGPU is available as an open-source stand-alone application and as an extension package for the R statistical environment. It provides a dramatic increase in performance for permutation resampling analysis in the context of microarray association studies. The current version offers six test statistics for carrying out permutation resampling analyses for binary, quantitative and censored time-to-event traits.

Robust test method for time-course microarray experiments.

BACKGROUND: In a time-course microarray experiment, the expression level for each gene is observed across a number of time-points in order to characterize the temporal trajectories of the gene-expression profiles. For many of these experiments, the scientific aim is the identification of genes for which the trajectories depend on an experimental or phenotypic factor. There is an extensive recent body of literature on statistical methodology for addressing this analytical problem. Most of the existing methods are based on estimating the time-course trajectories using parametric or non-parametric mean regression methods. The sensitivity of these regression methods to outliers, an issue that is well documented in the statistical literature, should be of concern when analyzing microarray data. RESULTS: In this paper, we propose a robust testing method for identifying genes whose expression time profiles depend on a factor. Furthermore, we propose a multiple testing procedure to adjust for multiplicity. CONCLUSIONS: Through an extensive simulation study, we will illustrate the performance of our method. Finally, we will report the results from applying our method to a case study and discussing potential extensions.

Comparison of pattern detection methods in microarray time series of the segmentation clock.

While genome-wide gene expression data are generated at an increasing rate, the repertoire of approaches for pattern discovery in these data is still limited. Identifying subtle patterns of interest in large amounts of data (tens of thousands of profiles) associated with a certain level of noise remains a challenge. A microarray time series was recently generated to study the transcriptional program of the mouse segmentation clock, a biological oscillator associated with the periodic formation of the segments of the body axis. A method related to Fourier analysis, the Lomb-Scargle periodogram, was used to detect periodic profiles in the dataset, leading to the identification of a novel set of cyclic genes associated with the segmentation clock. Here, we applied to the same microarray time series dataset four distinct mathematical methods to identify significant patterns in gene expression profiles. These methods are called: Phase consistency, Address reduction, Cyclohedron test and Stable persistence, and are based on different conceptual frameworks that are either hypothesis- or data-driven. Some of the methods, unlike Fourier transforms, are not dependent on the assumption of periodicity of the pattern of interest. Remarkably, these methods identified blindly the expression profiles of known cyclic genes as the most significant patterns in the dataset. Many candidate genes predicted by more than one approach appeared to be true positive cyclic genes and will be of particular interest for future research. In addition, these methods predicted novel candidate cyclic genes that were consistent with previous biological knowledge and experimental validation in mouse embryos. Our results demonstrate the utility of these novel pattern detection strategies, notably for detection of periodic profiles, and suggest that combining several distinct mathematical approaches to analyze microarray datasets is a valuable strategy for identifying genes that exhibit novel, interesting transcriptional patterns.

Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

Validation of the detection of Alexandrium species using specific RNA probes tested in a microarray format: Calibration of signal using variability of RNA content with environmental conditions.

The dinoflagellate genus Alexandrium contains several toxin producing species and strains, which can cause major economic losses to the shell fish industry. It is therefore important to be able to detect these toxin producers and also distinguish toxic strains from some of the morphologically identical non-toxic strains. To facilitate this DNA probes to be used in a microarray format were designed in silico or developed from existing published probes. These probes targeted either the 18S or 28S ribosomal ribonucleic acid (rRNA) gene in Alexandrium tamarense Group I, Group III and Group IV, Alexandrium ostenfeldii and Alexandrium minutum. Three strains of A. tamarense Group I, A. tamarense Group III, A. minutum and two strains of A. ostenfeldii were grown at optimal conditions and transferred into new environmental conditions changing either the light intensity, salinity, temperature or nutrient concentrations, to check if any of these environmental conditions induced changes in the cellular ribonucleic acid (RNA) concentration or growth rate. The aim of this experiment was the calibration of several species-specific probes for the quantification of the toxic Alexandrium strains. Growth rates were highly variable but only elevated or lowered salinity significantly lowered growth rate for A. tamarense Group I and Group III; differences in RNA content were not significant for the majority of the treatments. Only light intensity seemed to affect significantly the RNA content in A. tamarense Group I and Group III, but this was still within the same range as for the other treatments meaning that a back calibration from RNA to cell numbers was possible. The designed probes allow the production of quantitative information for Alexandrium species for the microarray chip.

Validation of the detection of Pseudo-nitzschia spp. using specific RNA probes tested in a microarray format: Calibration of signal based on variability of RNA content with environmental conditions.

Harmful algal blooms (HAB) occur worldwide and cause health problems and economic damage to fisheries and tourism. Monitoring for toxic algae is therefore essential but is based primarily on light microscopy, which is time consuming and can be limited by insufficient morphological characters such that more time is needed to examine critical features with electron microscopy. Monitoring with molecular tools is done in only a few places world-wide. EU FP7 MIDTAL (Microarray Detection of Toxic Algae) used SSU and LSU rRNA genes as targets on microarrays to identify toxic species. In order to comply with current monitoring requirements to report cell numbers as the relevant threshold measurement to trigger closure of fisheries, it was necessary to calibrate our microarray to convert the hybridisation signal obtained to cell numbers. Calibration curves for two species of Pseudo-nitzschia for use with the MIDTAL microarray are presented to obtain cell numbers following hybridisation. It complements work presented by Barra et al. (2012b. Environ. Sci. Pollut. Res. doi: 10.1007/s11356-012-1330-1v) for two other Pseudo-nitzschia spp., Dittami and Edvardsen (2012a. J. Phycol. 48, 1050) for Pseudochatonella, Blanco et al. (2013. Harmful Algae 24, 80) for Heterosigma, McCoy et al. (2013. FEMS. doi: 10.1111/1574-6941.12277) for Prymnesium spp., Karlodinium veneficum, and cf. Chatonella spp. and Taylor et al. (2014. Harmful Algae, in press) for Alexandrium.

Identification of 5-Fluorouracil-inducible Target Genes Using cDNA Microarray Profiling

The fluoropyrimidine 5-Fluorouracil (5-FU) is widely used in the treatment of cancer. To identify novel downstream mediators of tumor cell response to 5-FU, we used DNA microarray technology to identify genes that are transcriptionally activated by 5-FU treatment in the MCF-7 breast cancer cell line. Of 2400 genes analyzed, 619 were up-regulated by >3-fold. Highly up-regulated genes (>6-fold) with signal intensities of >3000 were analyzed by Northern blot. Genes that were consistently found to be up-regulated were spermine/spermidine acetyl transferase (SSAT), annexin II, thymosin-beta-10, chaperonin-10, and MAT-8. Treatment of MCF-7 cells with the antifolate tomudex and DNA-damaging agent oxaliplatin also resulted in up-regulation of each of these targets. The 5-FU-induced activation of MAT-8, thymosin-beta-10, and chaperonin-10 was abrogated by inactivation of p53 in MCF-7 cells, whereas induction of SSAT and annexin II was significantly reduced in the absence of p53. Moreover, each of these genes contained more than one potential p53-binding site, suggesting that p53 may play an important regulatory role in 5-FU-induced expression of these genes. In addition, we found that basal expression levels of SSAT, annexin II, thymosin beta-10, and chaperonin-10 were increased (by approximately 2-3-fold), and MAT-8 expression dramatically increased (by approximately 10-fold) in a 5-FU-resistant colorectal cancer cell line (H630-R10) compared with the parental H630 cell line, suggesting these genes may be useful biomarkers of resistance. These results demonstrate the potential of DNA microarrays to identify novel genes involved in mediating the response of tumor cells to chemotherapy.

Generation of a non-small cell lung cancer transcriptome microarray

Background: Non-small cell lung cancer (NSCLC) is the leading cause of cancer mortality worldwide. At present no reliable biomarkers are available to guide the management of this condition. Microarray technology may allow appropriate biomarkers to be identified but present platforms are lacking disease focus and are thus likely to miss potentially vital information contained in patient tissue samples.