Principal component analysis of binary genomics data.

Motivation: Genome-wide measurements of genetic and epigenetic alterations are generating more and more high-dimensional binary data. The special mathematical characteristics of binary data make the direct use of the classical principal component analysis (PCA) model to explore low-dimensional structures less obvious. Although there are several PCA alternatives for binary data in the psychometric, data analysis and machine learning literature, they are not well known to the bioinformatics community. Results: In this article, we introduce the motivation and rationale of some parametric and nonparametric versions of PCA specifically geared for binary data. Using both realistic simulations of binary data as well as mutation, CNA and methylation data of the Genomic Determinants of Sensitivity in Cancer 1000 (GDSC1000), the methods were explored for their performance with respect to finding the correct number of components, overfit, finding back the correct low-dimensional structure, variable importance, etc. The results show that if a low-dimensional structure exists in the data, that most of the methods can find it. When assuming a probabilistic generating process is underlying the data, we recommend to use the parametric logistic PCA model, while when such an assumption is not valid and the data are considered as given, the nonparametric Gifi model is recommended. Availability: The codes to reproduce the results in this article are available at the homepage of the Biosystems Data Analysis group (www.bdagroup.nl).

iTOP: inferring the topology of omics data.

Motivation: In biology, we are often faced with multiple datasets recorded on the same set of objects, such as multi-omics and phenotypic data of the same tumors. These datasets are typically not independent from each other. For example, methylation may influence gene expression, which may, in turn, influence drug response. Such relationships can strongly affect analyses performed on the data, as we have previously shown for the identification of biomarkers of drug response. Therefore, it is important to be able to chart the relationships between datasets. Results: We present iTOP, a methodology to infer a topology of relationships between datasets. We base this methodology on the RV coefficient, a measure of matrix correlation, which can be used to determine how much information is shared between two datasets. We extended the RV coefficient for partial matrix correlations, which allows the use of graph reconstruction algorithms, such as the PC algorithm, to infer the topologies. In addition, since multi-omics data often contain binary data (e.g. mutations), we also extended the RV coefficient for binary data. Applying iTOP to pharmacogenomics data, we found that gene expression acts as a mediator between most other datasets and drug response: only proteomics clearly shares information with drug response that is not present in gene expression. Based on this result, we used TANDEM, a method for drug response prediction, to identify which variables predictive of drug response were distinct to either gene expression or proteomics. Availability and implementation: An implementation of our methodology is available in the R package iTOP on CRAN. Additionally, an R Markdown document with code to reproduce all figures is provided as Supplementary Material. Supplementary information: Supplementary data are available at Bioinformatics online.

The ISCB Student Council Internship Program: Expanding computational biology capacity worldwide.

Education and training are two essential ingredients for a successful career. On one hand, universities provide students a curriculum for specializing in one's field of study, and on the other, internships complement coursework and provide invaluable training experience for a fruitful career. Consequently, undergraduates and graduates are encouraged to undertake an internship during the course of their degree. The opportunity to explore one's research interests in the early stages of their education is important for students because it improves their skill set and gives their career a boost. In the long term, this helps to close the gap between skills and employability among students across the globe and balance the research capacity in the field of computational biology. However, training opportunities are often scarce for computational biology students, particularly for those who reside in less-privileged regions. Aimed at helping students develop research and academic skills in computational biology and alleviating the divide across countries, the Student Council of the International Society for Computational Biology introduced its Internship Program in 2009. The Internship Program is committed to providing access to computational biology training, especially for students from developing regions, and improving competencies in the field. Here, we present how the Internship Program works and the impact of the internship opportunities so far, along with the challenges associated with this program.

The BRCA1ness signature is associated significantly with response to PARP inhibitor treatment versus control in the I-SPY 2 randomized neoadjuvant setting.

BACKGROUND: Patients with BRCA1-like tumors correlate with improved response to DNA double-strand break-inducing therapy. A gene expression-based classifier was developed to distinguish between BRCA1-like and non-BRCA1-like tumors. We hypothesized that these tumors may also be more sensitive to PARP inhibitors than standard treatments. METHODS: A diagnostic gene expression signature (BRCA1ness) was developed using a centroid model with 128 triple-negative breast cancer samples from the EU FP7 RATHER project. This BRCA1ness signature was then tested in HER2-negative patients (n = 116) from the I-SPY 2 TRIAL who received an oral PARP inhibitor veliparib in combination with carboplatin (V-C), or standard chemotherapy alone. We assessed the association between BRCA1ness and pathologic complete response in the V-C and control arms alone using Fisher's exact test, and the relative performance between arms (biomarker × treatment interaction, likelihood ratio p < 0.05) using a logistic model and adjusting for hormone receptor status (HR). RESULTS: We developed a gene expression signature to identify BRCA1-like status. In the I-SPY 2 neoadjuvant setting the BRCA1ness signature associated significantly with response to V-C (p = 0.03), but not in the control arm (p = 0.45). We identified a significant interaction between BRCA1ness and V-C (p = 0.023) after correcting for HR. CONCLUSIONS: A genomic-based BRCA1-like signature was successfully translated to an expression-based signature (BRC1Aness). In the I-SPY 2 neoadjuvant setting, we determined that the BRCA1ness signature is capable of predicting benefit of V-C added to standard chemotherapy compared to standard chemotherapy alone. TRIAL REGISTRATION: I-SPY 2 TRIAL beginning December 31, 2009: Neoadjuvant and Personalized Adaptive Novel Agents to Treat Breast Cancer (I-SPY 2), NCT01042379 .

TANDEM: a two-stage approach to maximize interpretability of drug response models based on multiple molecular data types.

MOTIVATION: Clinical response to anti-cancer drugs varies between patients. A large portion of this variation can be explained by differences in molecular features, such as mutation status, copy number alterations, methylation and gene expression profiles. We show that the classic approach for combining these molecular features (Elastic Net regression on all molecular features simultaneously) results in models that are almost exclusively based on gene expression. The gene expression features selected by the classic approach are difficult to interpret as they often represent poorly studied combinations of genes, activated by aberrations in upstream signaling pathways. RESULTS: To utilize all data types in a more balanced way, we developed TANDEM, a two-stage approach in which the first stage explains response using upstream features (mutations, copy number, methylation and cancer type) and the second stage explains the remainder using downstream features (gene expression). Applying TANDEM to 934 cell lines profiled across 265 drugs (GDSC1000), we show that the resulting models are more interpretable, while retaining the same predictive performance as the classic approach. Using the more balanced contributions per data type as determined with TANDEM, we find that response to MAPK pathway inhibitors is largely predicted by mutation data, while predicting response to DNA damaging agents requires gene expression data, in particular SLFN11 expression. AVAILABILITY AND IMPLEMENTATION: TANDEM is available as an R package on CRAN (for more information, see http://ccb.nki.nl/software/tandem). CONTACT: m.michaut@nki.nl or l.wessels@nki.nl SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A quick guide for building a successful bioinformatics community.

"Scientific community" refers to a group of people collaborating together on scientific-research-related activities who also share common goals, interests, and values. Such communities play a key role in many bioinformatics activities. Communities may be linked to a specific location or institute, or involve people working at many different institutions and locations. Education and training is typically an important component of these communities, providing a valuable context in which to develop skills and expertise, while also strengthening links and relationships within the community. Scientific communities facilitate: (i) the exchange and development of ideas and expertise; (ii) career development; (iii) coordinated funding activities; (iv) interactions and engagement with professionals from other fields; and (v) other activities beneficial to individual participants, communities, and the scientific field as a whole. It is thus beneficial at many different levels to understand the general features of successful, high-impact bioinformatics communities; how individual participants can contribute to the success of these communities; and the role of education and training within these communities. We present here a quick guide to building and maintaining a successful, high-impact bioinformatics community, along with an overview of the general benefits of participating in such communities. This article grew out of contributions made by organizers, presenters, panelists, and other participants of the ISMB/ECCB 2013 workshop "The 'How To Guide' for Establishing a Successful Bioinformatics Network" at the 21st Annual International Conference on Intelligent Systems for Molecular Biology (ISMB) and the 12th European Conference on Computational Biology (ECCB).

Next generation sequencing of triple negative breast cancer to find predictors for chemotherapy response.

INTRODUCTION: In triple negative breast cancers (TNBC) the initial response to chemotherapy is often favorable, but relapse and chemotherapy resistance frequently occur in advanced disease. Hence there is an urgent need for targeted treatments in this breast cancer subtype. In the current study we deep sequenced DNA of tumors prior to chemotherapy to search for predictors of response or resistance. METHODS: Next generation sequencing (NGS) was performed for 1,977 genes involved in tumorigenesis. DNA from 56 pre-treatment TNBC-biopsies was sequenced, as well as matched normal DNA. Following their tumor biopsy, patients started neoadjuvant chemotherapy with doxorubicin and cyclophosphamide. We studied associations between genetic alterations and three clinical variables: chemotherapy response, relapse-free survival and BRCA proficiency. RESULTS: The mutations observed were diverse and few recurrent mutations were detected. Most mutations were in TP53, TTN, and PIK3CA (55 %, 14 %, and 9 %, respectively). The mutation rates were similar between responders and non-responders (average mutation rate 9 vs 8 mutations). No recurrent mutations were associated with chemotherapy response or relapse. Interestingly, PIK3CA mutations were exclusively observed in patients proficient for BRCA1. Samples with a relapse had a higher copy number alteration rate, and amplifications of TTK and TP53BP2 were associated with a poor chemotherapy response. CONCLUSIONS: In this homogenous cohort of TNBCs few recurrent mutations were found. However, PIK3CA mutations were associated with BRCA proficiency, which can have clinical consequences in the near future.

Learn from the best.

What is more inspiring than a discussion with the leading scientists in your field? As a student or a young researcher, you have likely been influenced by mentors guiding you in your career and leading you to your current position. Any discussion with or advice from an expert is certainly very helpful for young people. But how often do we have the opportunity to meet experts? Do we make the most out of these situations? Meetings organized for young scientists are a great opportunity not only for the attendees: they are an opportunity for experts to meet bright students and learn from them in return. In this article, we introduce several successful events organized by Regional Student Groups all around the world, bridging the gap between experts and young scientists. We highlight how rewarding it is for all participants: young researchers, experts, and organizers. We then discuss the various benefits and emphasize the importance of organizing and attending such meetings. As a young researcher, seeking mentorship and additional skills training is a crucial step in career development. Keep in mind that one day, you may be an inspiring mentor, too.

Building the future of bioinformatics through student-facilitated conferencing.

Sharing results, techniques, and challenges is paramount to advance our understanding of any field of science. In the scientific community this exchange of ideas is mainly made possible through national and international conferences. Scientists have the opportunity to showcase their work, receive feedback, and improve their presentation skills. However, conferences can be large and intimidating for young researchers. In addition, for many of the more prestigious conferences, the very high number of submissions and low selection rate are major limitations to aspiring young researchers aiming to present their work to the scientific community. To improve student participation and proliferation of information, regional student groups have successfully organized conferences and symposia specifically aimed at students. This gives more students the opportunity to present their work and receive valuable experience and insight from peers and leaders in the field. At the same time, it is an ideal way for students to gain familiarity with the conference experience. In this paper, we highlight some of the benefits of participating in such student conferences, and we review the challenges we have encountered when organizing them. Both topics are illustrated in detail with examples from different ISCB Student Council Regional Student Groups.

Crossing borders for science.

Exchanging ideas with like-minded, enthusiastic people interested in the same topic is crucial for the advancement of a scientist's career. Several Regional Student Groups (RSGs) of the International Society for Computational Biology (ISCB) Student Council have cooperated in the last six years to organize scientific workshops and conferences. With motivated students, it is possible to create a memorable event for fellow scientists; in doing so, the organizers gain valuable experiences. While collaborating across borders and time zones can be difficult, feedback from event organizers was always positive. When limited resources are juxtaposed with great ideas and a network of contacts, the outcome is always an amazing experience, despite organizers being separated geographically across different countries.

Mapping the Hsp90 genetic interaction network in Candida albicans reveals environmental contingency and rewired circuitry.

The molecular chaperone Hsp90 regulates the folding of diverse signal transducers in all eukaryotes, profoundly affecting cellular circuitry. In fungi, Hsp90 influences development, drug resistance, and evolution. Hsp90 interacts with -10% of the proteome in the model yeast Saccharomyces cerevisiae, while only two interactions have been identified in Candida albicans, the leading fungal pathogen of humans. Utilizing a chemical genomic approach, we mapped the C. albicans Hsp90 interaction network under diverse stress conditions. The chaperone network is environmentally contingent, and most of the 226 genetic interactors are important for growth only under specific conditions, suggesting that they operate downstream of Hsp90, as with the MAPK Hog1. Few interactors are important for growth in many environments, and these are poised to operate upstream of Hsp90, as with the protein kinase CK2 and the transcription factor Ahr1. We establish environmental contingency in the first chaperone network of a fungal pathogen, novel effectors upstream and downstream of Hsp90, and network rewiring over evolutionary time.

Putting genetic interactions in context through a global modular decomposition.

Genetic interactions provide a powerful perspective into gene function, but our knowledge of the specific mechanisms that give rise to these interactions is still relatively limited. The availability of a global genetic interaction map in Saccharomyces cerevisiae, covering ∼30% of all possible double mutant combinations, provides an unprecedented opportunity for an unbiased assessment of the native structure within genetic interaction networks and how it relates to gene function and modular organization. Toward this end, we developed a data mining approach to exhaustively discover all block structures within this network, which allowed for its complete modular decomposition. The resulting modular structures revealed the importance of the context of individual genetic interactions in their interpretation and revealed distinct trends among genetic interaction hubs as well as insights into the evolution of duplicate genes. Block membership also revealed a surprising degree of multifunctionality across the yeast genome and enabled a novel association of VIP1 and IPK1 with DNA replication and repair, which is supported by experimental evidence. Our modular decomposition also provided a basis for testing the between-pathway model of negative genetic interactions and within-pathway model of positive genetic interactions. While we find that most modular structures involving negative genetic interactions fit the between-pathway model, we found that current models for positive genetic interactions fail to explain 80% of the modular structures detected. We also find differences between the modular structures of essential and nonessential genes.

Explain bioinformatics to your grandmother!

What are you working on? You have certainly been asked that question many times, whether it be at a Saturday night party, during a discussion with your neighbors, or at a family gathering. Communicating with a lay audience about scientific subjects and making them attractive is a difficult task. But difficult or not, you will have to do it for many years, not only with your family and friends, but also with your colleagues and collaborators. So, better learn now! Although not usually taught, the ability to explain your work to others is an essential skill in science, where communication plays a key role. Using some examples of the French Regional Student Group activities, we discuss here (i) why it is important to have such communication skills, (ii) how you can get involved in these activities by using existing resources or working with people who have previous experience, and (iii) what you get out of this amazing experience. We aim to motivate you and provide you with tips and ideas to get involved in promoting scientific activities while getting all the benefits.

The regional student group program of the ISCB student council: stories from the road.

The International Society for Computational Biology (ISCB) Student Council was launched in 2004 to facilitate interaction between young scientists in the fields of bioinformatics and computational biology. Since then, the Student Council has successfully run events and programs to promote the development of the next generation of computational biologists. However, in its early years, the Student Council faced a major challenge, in that students from different geographical regions had different needs; no single activity or event could address the needs of all students. To overcome this challenge, the Student Council created the Regional Student Group (RSG) program. The program consists of locally organised and run student groups that address the specific needs of students in their region. These groups usually encompass a given country, and, via affiliation with the international Student Council, are provided with financial support, organisational support, and the ability to share information with other RSGs. In the last five years, RSGs have been created all over the world and organised activities that have helped develop dynamic bioinformatics student communities. In this article series, we present common themes emerging from RSG initiatives, explain their goals, and highlight the challenges and rewards through specific examples. This article, the first in the series, introduces the Student Council and provides a high-level overview of RSG activities. Our hope is that the article series will be a valuable source of information and inspiration for initiating similar activities in other regions and scientific communities.

Distinct types of disorder in the human proteome: functional implications for alternative splicing.

Intrinsically disordered regions have been associated with various cellular processes and are implicated in several human diseases, but their exact roles remain unclear. We previously defined two classes of conserved disordered regions in budding yeast, referred to as "flexible" and "constrained" conserved disorder. In flexible disorder, the property of disorder has been positionally conserved during evolution, whereas in constrained disorder, both the amino acid sequence and the property of disorder have been conserved. Here, we show that flexible and constrained disorder are widespread in the human proteome, and are particularly common in proteins with regulatory functions. Both classes of disordered sequences are highly enriched in regions of proteins that undergo tissue-specific (TS) alternative splicing (AS), but not in regions of proteins that undergo general (i.e., not tissue-regulated) AS. Flexible disorder is more highly enriched in TS alternative exons, whereas constrained disorder is more highly enriched in exons that flank TS alternative exons. These latter regions are also significantly more enriched in potential phosphosites and other short linear motifs associated with cell signaling. We further show that cancer driver mutations are significantly enriched in regions of proteins associated with TS and general AS. Collectively, our results point to distinct roles for TS alternative exons and flanking exons in the dynamic regulation of protein interaction networks in response to signaling activity, and they further suggest that alternatively spliced regions of proteins are often functionally altered by mutations responsible for cancer.

A systems biology approach reveals the role of a novel methyltransferase in response to chemical stress and lipid homeostasis.

Using small molecule probes to understand gene function is an attractive approach that allows functional characterization of genes that are dispensable in standard laboratory conditions and provides insight into the mode of action of these compounds. Using chemogenomic assays we previously identified yeast Crg1, an uncharacterized SAM-dependent methyltransferase, as a novel interactor of the protein phosphatase inhibitor cantharidin. In this study we used a combinatorial approach that exploits contemporary high-throughput techniques available in Saccharomyces cerevisiae combined with rigorous biological follow-up to characterize the interaction of Crg1 with cantharidin. Biochemical analysis of this enzyme followed by a systematic analysis of the interactome and lipidome of CRG1 mutants revealed that Crg1, a stress-responsive SAM-dependent methyltransferase, methylates cantharidin in vitro. Chemogenomic assays uncovered that lipid-related processes are essential for cantharidin resistance in cells sensitized by deletion of the CRG1 gene. Lipidome-wide analysis of mutants further showed that cantharidin induces alterations in glycerophospholipid and sphingolipid abundance in a Crg1-dependent manner. We propose that Crg1 is a small molecule methyltransferase important for maintaining lipid homeostasis in response to drug perturbation. This approach demonstrates the value of combining chemical genomics with other systems-based methods for characterizing proteins and elucidating previously unknown mechanisms of action of small molecule inhibitors.

Genetic interaction maps in Escherichia coli reveal functional crosstalk among cell envelope biogenesis pathways.

As the interface between a microbe and its environment, the bacterial cell envelope has broad biological and clinical significance. While numerous biosynthesis genes and pathways have been identified and studied in isolation, how these intersect functionally to ensure envelope integrity during adaptive responses to environmental challenge remains unclear. To this end, we performed high-density synthetic genetic screens to generate quantitative functional association maps encompassing virtually the entire cell envelope biosynthetic machinery of Escherichia coli under both auxotrophic (rich medium) and prototrophic (minimal medium) culture conditions. The differential patterns of genetic interactions detected among > 235,000 digenic mutant combinations tested reveal unexpected condition-specific functional crosstalk and genetic backup mechanisms that ensure stress-resistant envelope assembly and maintenance. These networks also provide insights into the global systems connectivity and dynamic functional reorganization of a universal bacterial structure that is both broadly conserved among eubacteria (including pathogens) and an important target.

Multiple genetic interaction experiments provide complementary information useful for gene function prediction.

Genetic interactions help map biological processes and their functional relationships. A genetic interaction is defined as a deviation from the expected phenotype when combining multiple genetic mutations. In Saccharomyces cerevisiae, most genetic interactions are measured under a single phenotype - growth rate in standard laboratory conditions. Recently genetic interactions have been collected under different phenotypic readouts and experimental conditions. How different are these networks and what can we learn from their differences? We conducted a systematic analysis of quantitative genetic interaction networks in yeast performed under different experimental conditions. We find that networks obtained using different phenotypic readouts, in different conditions and from different laboratories overlap less than expected and provide significant unique information. To exploit this information, we develop a novel method to combine individual genetic interaction data sets and show that the resulting network improves gene function prediction performance, demonstrating that individual networks provide complementary information. Our results support the notion that using diverse phenotypic readouts and experimental conditions will substantially increase the amount of gene function information produced by genetic interaction screens.

A robust genomic signature for the detection of colorectal cancer patients with microsatellite instability phenotype and high mutation frequency.

Microsatellite instability (MSI) occurs in 10-20% of colorectal tumours and is associated with good prognosis. Here we describe the development and validation of a genomic signature that identifies colorectal cancer patients with MSI caused by DNA mismatch repair deficiency with high accuracy. Microsatellite status for 276 stage II and III colorectal tumours has been determined. Full-genome expression data was used to identify genes that correlate with MSI status. A subset of these samples (n = 73) had sequencing data for 615 genes available. An MSI gene signature of 64 genes was developed and validated in two independent validation sets: the first consisting of frozen samples from 132 stage II patients; and the second consisting of FFPE samples from the PETACC-3 trial (n = 625). The 64-gene MSI signature identified MSI patients in the first validation set with a sensitivity of 90.3% and an overall accuracy of 84.8%, with an AUC of 0.942 (95% CI, 0.888-0.975). In the second validation, the signature also showed excellent performance, with a sensitivity 94.3% and an overall accuracy of 90.6%, with an AUC of 0.965 (95% CI, 0.943-0.988). Besides correct identification of MSI patients, the gene signature identified a group of MSI-like patients that were MSS by standard assessment but MSI by signature assessment. The MSI-signature could be linked to a deficient MMR phenotype, as both MSI and MSI-like patients showed a high mutation frequency (8.2% and 6.4% of 615 genes assayed, respectively) as compared to patients classified as MSS (1.6% mutation frequency). The MSI signature showed prognostic power in stage II patients (n = 215) with a hazard ratio of 0.252 (p = 0.0145). Patients with an MSI-like phenotype had also an improved survival when compared to MSS patients. The MSI signature was translated to a diagnostic microarray and technically and clinically validated in FFPE and frozen samples.

Protein complexes are central in the yeast genetic landscape.

If perturbing two genes together has a stronger or weaker effect than expected, they are said to genetically interact. Genetic interactions are important because they help map gene function, and functionally related genes have similar genetic interaction patterns. Mapping quantitative (positive and negative) genetic interactions on a global scale has recently become possible. This data clearly shows groups of genes connected by predominantly positive or negative interactions, termed monochromatic groups. These groups often correspond to functional modules, like biological processes or complexes, or connections between modules. However it is not yet known how these patterns globally relate to known functional modules. Here we systematically study the monochromatic nature of known biological processes using the largest quantitative genetic interaction data set available, which includes fitness measurements for ∼5.4 million gene pairs in the yeast Saccharomyces cerevisiae. We find that only 10% of biological processes, as defined by Gene Ontology annotations, and less than 1% of inter-process connections are monochromatic. Further, we show that protein complexes are responsible for a surprisingly large fraction of these patterns. This suggests that complexes play a central role in shaping the monochromatic landscape of biological processes. Altogether this work shows that both positive and negative monochromatic patterns are found in known biological processes and in their connections and that protein complexes play an important role in these patterns. The monochromatic processes, complexes and connections we find chart a hierarchical and modular map of sensitive and redundant biological systems in the yeast cell that will be useful for gene function prediction and comparison across phenotypes and organisms. Furthermore the analysis methods we develop are applicable to other species for which genetic interactions will progressively become more available.