MOTIVATION: A fundamental problem for disease treatment is that while antibiotics are a powerful counter to bacteria, they are ineffective against viruses. Often, bacterial and viral infections are confused due to their similar symptoms and lack of rapid diagnostics. With many clinicians relying primarily on symptoms for diagnosis, overuse and misuse of modern antibiotics are rife, contributing to the growing pool of antibiotic resistance. To ensure an individual receives optimal treatment given their disease state and to reduce over-prescription of antibiotics, the host response can in theory be measured quickly to distinguish between the two states. To establish a predictive biomarker panel of disease state (viral/bacterial/no-infection), we conducted a meta-analysis of human blood infection studies using machine learning. RESULTS: We focused on publicly available gene expression data from two widely used platforms, Affymetrix and Illumina microarrays as they represented a significant proportion of the available data. We were able to develop multi-class models with high accuracies with our best model predicting 93% of bacterial and 89% viral samples correctly. To compare the selected features in each of the different technologies, we reverse-engineered the underlying molecular regulatory network and explored the neighbourhood of the selected features. The networks highlighted that although on the gene-level the models differed, they contained genes from the same areas of the network. Specifically, this convergence was to pathways including the Type I interferon Signalling Pathway, Chemotaxis, Apoptotic Processes and Inflammatory/Innate Response. AVAILABILITY: Data and code are available on the Gene Expression Omnibus and github. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MOTIVATION: The molecular processes regulating molluscan shell production remain relatively uncharacterized, despite the clear evolutionary and societal importance of biomineralization. RESULTS: Here we built the first computationally predicted gene regulatory network (GRN) for molluscan biomineralization using Antarctic clam (Laternula elliptica) mantle gene expression data produced over an age-categorized shell damage-repair time-course. We used previously published in vivo in situ hybridization expression data to ground truth gene interactions predicted by the GRN and show that candidate biomineralization genes from different shell layers, and hence microstructures, were connected in unique modules. We characterized two biomineralization modules of the GRN and hypothesize that one module is responsible for translating the extracellular proteins required for growing, repairing or remodelling the nacreous shell layer, whereas the second module orchestrates the transport of both ions and proteins to the shell secretion site, which are required during normal shell growth, and repair. Our findings demonstrate that unbiased computational methods are particularly valuable for studying fundamental biological processes and gene interactions in non-model species where rich sources of gene expression data exist, but annotation rates are poor and the ability to carry out true functional tests are still lacking. AVAILABILITY AND IMPLEMENTATION: The raw RNA-Seq data is freely available for download from NCBI SRA (Accession: PRJNA398984), the assembled and annotated transcriptome can be viewed and downloaded from molluscDB (ensembl.molluscdb.org) and in addition, the assembled transcripts, reconstructed GRN, modules and detailed annotations are all available as Supplementary Files. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MOTIVATION: Childhood B-precursor lymphoblastic leukaemia (ALL) is the most common paediatric malignancy. Despite the fact that 80% of ALL patients respond to anti-cancer drugs, the patho-physiology of this disease is still not fully understood. mRNA expression-profiling studies that have been performed have not yet provided novel insights into the mechanisms behind cellular response to DNA damage. More powerful data analysis techniques may be required for identifying novel functional pathways involved in the cellular responses to DNA damage. RESULTS: In order to explore the possibility that unforeseen biological processes may be involved in the response to DNA damage, we have developed and applied a novel procedure for the identification of functional modules in ALL cells. We have discovered that the overall activity of functional modules integrating protein degradation and mRNA processing is predictive of response to DNA damage. AVAILABILITY: Supplementary material including R code, additional results, experimental datasets, as well as a detailed description of the methodology are available at http://www.bip.bham.ac.uk/vivo/fumo.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.