Willem Talloen

FABIA: factor analysis for bicluster acquisition

Motivation: Biclustering of transcriptomic data groups genes and samples simultaneously. It is emerging as a standard tool for extracting knowledge from gene expression measurements. We propose a novel generative approach for biclustering called ‘FABIA: Factor Analysis for Bicluster Acquisition’. FABIA is based on a multiplicative model, which accounts for linear dependencies between gene expression and conditions, and also captures heavy-tailed distributions as observed in real-world transcriptomic data. The generative framework allows to utilize well-founded model selection methods and to apply Bayesian techniques. Results: On 100 simulated datasets with known true, artificially implanted biclusters, FABIA clearly outperformed all 11 competitors. On these datasets, FABIA was able to separate spurious biclusters from true biclusters by ranking biclusters according to their information content. FABIA was tested on three microarray datasets with known subclusters, where it was two times the best and once the second best method among the compared biclustering approaches. Availability: FABIA is available as an R package on Bioconductor (http://www.bioconductor.org). All datasets, results and software are available at http://www.bioinf.jku.at/software/fabia/fabia.html Contact: hochreit@bioinf.jku.at Supplementary information: Supplementary data are available at Bioinformatics online.

Quantitation of pretreatment serum interferon‐γ–inducible protein‐10 improves the predictive value of an IL28B gene polymorphism for hepatitis C treatment response

Polymorphisms of the IL28B gene are highly associated with sustained virological response (SVR) in patients with chronic hepatitis C treated with peginterferon and ribavirin. Quantitation of interferon‐γ–inducible protein‐10 (IP‐10) may also differentiate antiviral response. We evaluated IP‐10 levels in pretreatment serum from 115 nonresponders and 157 sustained responders in the Study of Viral Resistance to Antiviral Therapy of Chronic Hepatitis C cohort, including African American (AA) and Caucasian American (CA) patients. Mean IP‐10 was lower in sustained responders compared with nonresponders (437 ± 31 vs 704 ± 44 pg/mL, P < 0.001), both in AA and CA patients. The positive predictive value of low IP‐10 levels (<600 pg/mL) for SVR was 69%, whereas the negative predictive value of high IP‐10 levels (>600 pg/mL) was 67%. We assessed the combination of pretreatment IP‐10 levels with IL28B genotype as predictors of treatment response. The IL28B polymorphism rs12979860 was tested in 210 participants. The CC, CT, and TT genotypes were found in 30%, 49%, and 21% of patients, respectively, with corresponding SVR rates of 87%, 50%, and 39% (P < 0.0001). Serum IP‐10 levels within the IL28B genotype groups provided additional information regarding the likelihood of SVR (P < 0.0001). CT carriers with low IP‐10 had 64% SVR versus 24% with high IP‐10. Similarly, a higher SVR rate was identified for TT and CC carriers with low versus high IP‐10 (TT, 48% versus 20%; CC, 89% versus 79%). IL28B genotype and baseline IP‐10 levels were additive but independent when predicting SVR in both AA and CA patients. Conclusion: When IL28B genotype is combined with pretreatment serum IP‐10 measurement, the predictive value for discrimination between SVR and nonresponse is significantly improved, especially in non‐CC genotypes. This relationship warrants further investigation to elucidate the mechanisms of antiviral response and prospective validation. (Hepatology 2011;)

Alterations in Mucosal Immunity Identified in the Colon of Patients With Irritable Bowel Syndrome

BACKGROUND & AIMS Irritable bowel syndrome (IBS) has been associated with mucosal dysfunction, mild inflammation, and altered colonic bacteria. We used microarray expression profiling of sigmoid colon mucosa to assess whether there are stably expressed sets of genes that suggest there are objective molecular biomarkers associated with IBS. METHODS Gene expression profiling was performed using Human Genome U133 Plus 2.0 (Affymetrix) GeneChips with RNA from sigmoid colon mucosal biopsy specimens from 36 IBS patients and 25 healthy control subjects. Real-time quantitative polymerase chain reaction was used to confirm the data in 12 genes of interest. Statistical methods for microarray data were applied to search for differentially expressed genes, and to assess the stability of molecular signatures in IBS patients. RESULTS Mucosal gene expression profiles were consistent across different sites within the sigmoid colon and were stable on repeat biopsy over approximately 3 months. Differentially expressed genes suggest functional alterations of several components of the host mucosal immune response to microbial pathogens. The most strikingly increased expression involved a yet uncharacterized gene, DKFZP564O0823. Identified specific genes suggest the hypothesis that molecular signatures may enable distinction of a subset of IBS patients from healthy controls. By using 75% of the biopsy specimens as a validation set to develop a gene profile, the test set (25%) was predicted correctly with approximately 70% accuracy. CONCLUSIONS Mucosal gene expression analysis shows there are relatively stable alterations in colonic mucosal immunity in IBS. These molecular alterations provide the basis to test the hypothesis that objective biomarkers may be identified in IBS and enhance understanding of the disease.

I/NI-calls for the exclusion of non-informative genes

MOTIVATION DNA microarray technology typically generates many measurements of which only a relatively small subset is informative for the interpretation of the experiment. To avoid false positive results, it is therefore critical to select the informative genes from the large noisy data before the actual analysis. Most currently available filtering techniques are supervised and therefore suffer from a potential risk of overfitting. The unsupervised filtering techniques, on the other hand, are either not very efficient or too stringent as they may mix up signal with noise. We propose to use the multiple probes measuring the same target mRNA as repeated measures to quantify the signal-to-noise ratio of that specific probe set. A Bayesian factor analysis with specifically chosen prior settings, which models this probe level information, is providing an objective feature filtering technique, named informative/non-informative calls (I/NI calls). RESULTS Based on 30 real-life data sets (including various human, rat, mice and Arabidopsis studies) and a spiked-in data set, it is shown that I/NI calls is highly effective, with exclusion rates ranging from 70% to 99%. Consequently, it offers a critical solution to the curse of high-dimensionality in the analysis of microarray data. AVAILABILITY This filtering approach is publicly available as a function implemented in the R package FARMS (www.bioinf.jku.at/software/farms/farms.html).

Genome-wide copy number alterations detection in fresh frozen and matched FFPE samples using SNP 6.0 arrays

SNP arrays offer the opportunity to get a genome‐wide view on copy number alterations and are increasingly used in oncology. DNA from formalin‐fixed paraffin‐embedded material (FFPE) is partially degraded which limits the application of those technologies for retrospective studies. We present the use of Affymetrix GeneChip SNP6.0 for identification of copy number alterations in fresh frozen (FF) and matched FFPE samples. Fifteen pairs of adenocarcinomas with both frozen and FFPE embedded material were analyzed. We present an optimization of the sample preparation and show the importance of correcting the measured intensities for fragment length and GC‐content when using FFPE samples. The absence of GC content correction results in a chromosome specific “wave pattern” which may lead to the misclassification of genomic regions as being altered. The highest concordance between FFPE and matched FF were found in samples with the highest call rates. Nineteen of the 23 high level amplifications (83%) seen using FF samples were also detected in the corresponding FFPE material. For limiting the rate of “false positive” alterations, we have chosen a conservative False Discovery Rate (FDR). We observed better results using SNP probes than CNV probes for copy number analysis of FFPE material. This is the first report on the detection of copy number alterations in FFPE samples using Affymetrix GeneChip SNP6.0.

Response prediction to a multitargeted kinase inhibitor in cancer cell lines and xenograft tumors using high-content tyrosine peptide arrays with a kinetic readout

Multitargeted kinase inhibitors have shown clinical efficacy in a range of cancer types. However, two major problems associated with these drugs are the low fraction of patients for which these treatments provide initial clinical benefit and the occurrence of resistance during prolonged therapy. Several types of predictive biomarkers have been suggested, such as expression level and phosphorylation status of the major targeted kinase(s), mutational status of the kinases involved and of key components of the downstream signaling cascades, and gene expression signatures. In this work, we describe the development of a response prediction platform that does not require prior knowledge of the relevant kinases targeted by the inhibitor; instead, a phosphotyrosine peptide profile using peptide arrays with a kinetic readout is derived in lysates in the presence and absence of a kinase inhibitor. We show in a range of cell lines and in xenograft tumors that this approach allows for the stratification of responders and nonresponders to a multitargeted kinase inhibitor. [Mol Cancer Ther 2009;8(7):1846–55]

Using transcriptomics to guide lead optimization in drug discovery projects: Lessons learned from the QSTAR project.

The pharmaceutical industry is faced with steadily declining R&D efficiency which results in fewer drugs reaching the market despite increased investment. A major cause for this low efficiency is the failure of drug candidates in late-stage development owing to safety issues or previously undiscovered side-effects. We analyzed to what extent gene expression data can help to de-risk drug development in early phases by detecting the biological effects of compounds across disease areas, targets and scaffolds. For eight drug discovery projects within a global pharmaceutical company, gene expression data were informative and able to support go/no-go decisions. Our studies show that gene expression profiling can detect adverse effects of compounds, and is a valuable tool in early-stage drug discovery decision making.

cn.FARMS: a latent variable model to detect copy number variations in microarray data with a low false discovery rate

Cost-effective oligonucleotide genotyping arrays like the Affymetrix SNP 6.0 are still the predominant technique to measure DNA copy number variations (CNVs). However, CNV detection methods for microarrays overestimate both the number and the size of CNV regions and, consequently, suffer from a high false discovery rate (FDR). A high FDR means that many CNVs are wrongly detected and therefore not associated with a disease in a clinical study, though correction for multiple testing takes them into account and thereby decreases the studys discovery power. For controlling the FDR, we propose a probabilistic latent variable model, ‘cn.FARMS’, which is optimized by a Bayesian maximum a posteriori approach. cn.FARMS controls the FDR through the information gain of the posterior over the prior. The prior represents the null hypothesis of copy number 2 for all samples from which the posterior can only deviate by strong and consistent signals in the data. On HapMap data, cn.FARMS clearly outperformed the two most prevalent methods with respect to sensitivity and FDR. The software cn.FARMS is publicly available as a R package at http://www.bioinf.jku.at/software/cnfarms/cnfarms.html.

EEG alpha power as an intermediate measure between brain-derived neurotrophic factor Val66Met and depression severity in patients with major depressive disorder.

Summary: Major depressive disorder has a large impact on patients and society and is projected to be the second greatest global burden of disease by 2020. The brain-derived neurotrophic factor (BDNF) gene is considered to be one of the important factors in the etiology of major depressive disorder. In a recent study, alpha power was found to mediate between BDNF Met and subclinical depressed mood. The current study looked at a population of patients with major depressive disorder (N = 107) to examine the association between the BDNF Val66Met polymorphism, resting state EEG alpha power, and depression severity. For this purpose, repeated-measures analysis of variance, partial correlation, and multiple linear models were used. Results indicated a negative association between parietal-occipital alpha power in the eyes open resting state and depression severity. In addition, Met/Met patients showed lower global absolute alpha power in the eyes closed condition compared with Val-carriers. These findings are in accordance with the previously uncovered pathway between BDNF Val66Met, resting state EEG alpha power, and depression severity. Additional research is needed for the clarification of this tentative pathway and its implication in personalized treatment of major depressive disorder.

VirVarSeq: a low-frequency virus variant detection pipeline for Illumina sequencing using adaptive base-calling accuracy filtering

MOTIVATION In virology, massively parallel sequencing (MPS) opens many opportunities for studying viral quasi-species, e.g. in HIV-1- and HCV-infected patients. This is essential for understanding pathways to resistance, which can substantially improve treatment. Although MPS platforms allow in-depth characterization of sequence variation, their measurements still involve substantial technical noise. For Illumina sequencing, single base substitutions are the main error source and impede powerful assessment of low-frequency mutations. Fortunately, base calls are complemented with quality scores (Qs) that are useful for differentiating errors from the real low-frequency mutations. RESULTS A variant calling tool, Q-cpileup, is proposed, which exploits the Qs of nucleotides in a filtering strategy to increase specificity. The tool is imbedded in an open-source pipeline, VirVarSeq, which allows variant calling starting from fastq files. Using both plasmid mixtures and clinical samples, we show that Q-cpileup is able to reduce the number of false-positive findings. The filtering strategy is adaptive and provides an optimized threshold for individual samples in each sequencing run. Additionally, linkage information is kept between single-nucleotide polymorphisms as variants are called at the codon level. This enables virologists to have an immediate biological interpretation of the reported variants with respect to their antiviral drug responses. A comparison with existing SNP caller tools reveals that calling variants at the codon level with Q-cpileup results in an outstanding sensitivity while maintaining a good specificity for variants with frequencies down to 0.5%. AVAILABILITY The VirVarSeq is available, together with a users guide and test data, at sourceforge: http://sourceforge.net/projects/virtools/?source=directory.