Search

INTRODUCTION: Preeclampsia is a life-threatening disorder of pregnancy unique to humans. Interleukin (IL)11 is elevated in serum from pregnancies that subsequently develop early-onset preeclampsia and pharmacological elevation of IL11 in pregnant mice causes the development of early-onset preeclampsia-like features (hypertension, proteinuria, and fetal growth restriction). However, the mechanism by which IL11 drives preeclampsia is unknown. METHOD: Pregnant mice were administered PEGylated (PEG)IL11 or control (PEG) from embryonic day (E)10-16 and the effect on inflammasome activation, systolic blood pressure (during gestation and at 50/90 days post-natal), placental development, and fetal/post-natal pup growth measured. RNAseq analysis was performed on E13 placenta. Human 1st trimester placental villi were treated with IL11 and the effect on inflammasome activation and pyroptosis identified by immunohistochemistry and ELISA. RESULT: PEGIL11 activated the placental inflammasome causing inflammation, fibrosis, and acute and chronic hypertension in wild-type mice. Global and placental-specific loss of the inflammasome adaptor protein Asc and global loss of the Nlrp3 sensor protein prevented PEGIL11-induced fibrosis and hypertension in mice but did not prevent PEGIL11-induced fetal growth restriction or stillbirths. RNA-sequencing and histology identified that PEGIL11 inhibited trophoblast differentiation towards spongiotrophoblast and syncytiotrophoblast lineages in mice and extravillous trophoblast lineages in human placental villi. DISCUSSION: Inhibition of ASC/NLRP3 inflammasome activity could prevent IL11-induced inflammation and fibrosis in various disease states including preeclampsia.

The microbiome is a complex and dynamic community of microorganisms that co-exist interdependently within an ecosystem, and interact with its host or environment. Longitudinal studies can capture temporal variation within the microbiome to gain mechanistic insights into microbial systems; however, current statistical methods are limited due to the complex and inherent features of the data. We have identified three analytical objectives in longitudinal microbial studies: (1) differential abundance over time and between sample groups, demographic factors or clinical variables of interest; (2) clustering of microorganisms evolving concomitantly across time and (3) network modelling to identify temporal relationships between microorganisms. This review explores the strengths and limitations of current methods to fulfill these objectives, compares different methods in simulation and case studies for objectives (1) and (2), and highlights opportunities for further methodological developments. R tutorials are provided to reproduce the analyses conducted in this review.

Evolving in sync with the computation revolution over the past 30 years, computational biology has emerged as a mature scientific field. While the field has made major contributions toward improving scientific knowledge and human health, individual computational biology practitioners at various institutions often languish in career development. As optimistic biologists passionate about the future of our field, we propose solutions for both eager and reluctant individual scientists, institutions, publishers, funding agencies, and educators to fully embrace computational biology. We believe that in order to pave the way for the next generation of discoveries, we need to improve recognition for computational biologists and better align pathways of career success with pathways of scientific progress. With 10 outlined steps, we call on all adjacent fields to move away from the traditional individual, single-discipline investigator research model and embrace multidisciplinary, data-driven, team science.

BACKGROUND: Vaccination remains one of the most effective means of reducing the burden of infectious diseases globally. Improving our understanding of the molecular basis for effective vaccine response is of paramount importance if we are to ensure the success of future vaccine development efforts. METHODS: We applied cutting edge multi-omics approaches to extensively characterize temporal molecular responses following vaccination with hepatitis B virus (HBV) vaccine. Data were integrated across cellular, epigenomic, transcriptomic, proteomic, and fecal microbiome profiles, and correlated to final HBV antibody titres. RESULTS: Using both an unsupervised molecular-interaction network integration method (NetworkAnalyst) and a data-driven integration approach (DIABLO), we uncovered baseline molecular patterns and pathways associated with more effective vaccine responses to HBV. Biological associations were unravelled, with signalling pathways such as JAK-STAT and interleukin signalling, Toll-like receptor cascades, interferon signalling, and Th17 cell differentiation emerging as important pre-vaccination modulators of response. CONCLUSION: This study provides further evidence that baseline cellular and molecular characteristics of an individual's immune system influence vaccine responses, and highlights the utility of integrating information across many parallel molecular datasets.

Ultraviolet radiation-induced DNA mutations are a primary environmental driver of melanoma. The reason for this very high level of unrepaired DNA lesions leading to these mutations is still poorly understood. The primary DNA repair mechanism for UV-induced lesions, that is, the nucleotide excision repair pathway, appears intact in most melanomas. We have previously reported a postreplication repair mechanism that is commonly defective in melanoma cell lines. Here we have used a genome-wide approach to identify the components of this postreplication repair mechanism. We have used differential transcript polysome loading to identify transcripts that are associated with UV response, and then functionally assessed these to identify novel components of this repair and cell cycle checkpoint network. We have identified multiple interaction nodes, including global genomic nucleotide excision repair and homologous recombination repair, and previously unexpected MASTL pathway, as components of the response. Finally, we have used bioinformatics to assess the contribution of dysregulated expression of these pathways to the UV signature mutation load of a large melanoma cohort. We show that dysregulation of the pathway, especially the DNA damage repair components, are significant contributors to UV mutation load, and that dysregulation of the MASTL pathway appears to be a significant contributor to high UV signature mutation load.

Gene expression atlases have transformed our understanding of the development, composition and function of human tissues. New technologies promise improved cellular or molecular resolution, and have led to the identification of new cell types, or better defined cell states. But as new technologies emerge, information derived on old platforms becomes obsolete. We demonstrate that it is possible to combine a large number of different profiling experiments summarised from dozens of laboratories and representing hundreds of donors, to create an integrated molecular map of human tissue. As an example, we combine 850 samples from 38 platforms to build an integrated atlas of human blood cells. We achieve robust and unbiased cell type clustering using a variance partitioning method, selecting genes with low platform bias relative to biological variation. Other than an initial rescaling, no other transformation to the primary data is applied through batch correction or renormalisation. Additional data, including single-cell datasets, can be projected for comparison, classification and annotation. The resulting atlas provides a multi-scaled approach to visualise and analyse the relationships between sets of genes and blood cell lineages, including the maturation and activation of leukocytes in vivo and in vitro. In allowing for data integration across hundreds of studies, we address a key reproduciblity challenge which is faced by any new technology. This allows us to draw on the deep phenotypes and functional annotations that accompany traditional profiling methods, and provide important context to the high cellular resolution of single cell profiling. Here, we have implemented the blood atlas in the open access Stemformatics.org platform, drawing on its extensive collection of curated transcriptome data. The method is simple, scalable and amenable for rapid deployment in other biological systems or computational workflows.

BACKGROUND: Stimulated C-peptide measurement after a mixed meal tolerance test (MMTT) is the accepted gold standard for assessing residual beta-cell function in type 1 diabetes (T1D); however, this approach is impractical outside of clinical trials. OBJECTIVE: To develop an improved estimate of residual beta-cell function in children with T1D using commonly measured clinical variables. SUBJECTS/METHODS: A clinical model to predict 90-minute MMTT stimulated C-peptide in children with recent-onset T1D was developed from the combined AbATE, START, and TIDAL placebo subjects (n = 46) 6 months post-recruitment using multiple linear regression. This model was then validated in a clinical cohort (Hvidoere study group, n = 262). RESULTS: A model of estimated C-peptide at 6 months post-diagnosis, which included age, gender, body mass index (BMI), hemoglobin A1c (HbA1c), and insulin dose predicted 90-minute stimulated C-peptide measurements (adjusted R2 = 0.63, P < 0.0001). The predictive value of insulin dose and HbA1c alone (IDAA1c) for 90-minute stimulated C-peptide was significantly lower (R2 = 0.37, P < 0.0001). The slopes of linear regression lines of the estimated and stimulated 90-minute C-peptide levels obtained at 6 and 12 months post diagnosis in the Hvidoere clinical cohort were R2 = 0.36, P < 0.0001 at 6 months and R2 = 0.37, P < 0.0001 at 12 months. CONCLUSIONS: A clinical model including age, gender, BMI, HbA1c, and insulin dose predicts stimulated C-peptide levels in children with recent-onset T1D. Estimated C-peptide is an improved surrogate to monitor residual beta-cell function outside clinical trial settings.

Early life nutrition is a vital determinant of an individual's life-long health and also directly influences the ecological and functional development of the gut microbiota. However, there are limited longitudinal studies examining the effect of diet on the gut microbiota development in early childhood. Here, up to seven stool samples were collected from each of 48 healthy children during their second year of life, and microbiota dynamics were assessed using 16S rRNA gene amplicon sequencing. Children's dietary information was also collected during the same period using a validated food frequency questionnaire designed for this age group, over five time points. We observed significant changes in gut microbiota community, concordant with changes in the children's dietary pattern over the 12-month period. In particular, we found differential effects on specific Firmicutes-affiliated lineages in response to frequent intake of either processed or unprocessed foods. Additionally, the consumption of fortified milk supplemented with a Bifidobacterium probiotic and prebiotics (synbiotics) further increased the presence of Bifidobacterium spp., highlighting the potential use of synbiotics to prolong and sustain changes in these lineages and shaping the gut microbiota community in young children.

Human microbiomes are predicted to assemble in a reproducible and ordered manner yet there is limited knowledge on the development of the complex bacterial communities that constitute the oral microbiome. The oral microbiome plays major roles in many oral diseases including early childhood caries (ECC), which afflicts up to 70% of children in some countries. Saliva contains oral bacteria that are indicative of the whole oral microbiome and may have the ability to reflect the dysbiosis in supragingival plaque communities that initiates the clinical manifestations of ECC. The aim of this study was to determine the assembly of the oral microbiome during the first four years of life and compare it with the clinical development of ECC. The oral microbiomes of 134 children enrolled in a birth cohort study were determined at six ages between two months and four years-of-age and their mother's oral microbiome was determined at a single time point. We identified and quantified 356 operational taxonomic units (OTUs) of bacteria in saliva by sequencing the V4 region of the bacterial 16S RNA genes. Bacterial alpha diversity increased from a mean of 31 OTUs in the saliva of infants at 1.9 months-of-age to 84 OTUs at 39 months-of-age. The oral microbiome showed a distinct shift in composition as the children matured. The microbiome data were compared with the clinical development of ECC in the cohort at 39, 48, and 60 months-of-age as determined by ICDAS-II assessment. Streptococcus mutans was the most discriminatory oral bacterial species between health and current disease, with an increased abundance in disease. Overall our study demonstrates an ordered temporal development of the oral microbiome, describes a limited core oral microbiome and indicates that saliva testing of infants may help predict ECC risk.

Cutaneous squamous cell carcinoma (SCC) is the second-most-common cancer in Australia. The majority of SCCs progress from premalignant actinic keratosis (AK) lesions that form on chronically sun-exposed skin. The role of skin microbiota in this progression is not well understood; therefore, we performed a longitudinal microbiome analysis of AKs and SCCs using a cohort of 13 SCC-prone immunocompetent men. The majority of variability in microbial profiles was attributable to subject, followed by time and lesion type. Propionibacterium and Malassezia organisms were relatively more abundant in nonlesional photodamaged skin than in AKs and SCCs. Staphylococcus was most commonly associated with lesional skin, in particular, sequences most closely related to Staphylococcus aureus Of 11 S. aureus-like operational taxonomic units (OTUs), six were significantly associated with SCC lesions across seven subjects, suggesting their specific involvement with AK-to-SCC progression. If a causative link exists between certain S. aureus-like OTUs and SCC etiology, therapeutic approaches specifically targeting these bacteria could be used to reduce SCC.IMPORTANCE Actinic keratosis (AK) and cutaneous squamous cell carcinoma (SCC) are two of the most common dermatologic conditions in Western countries and cause substantial morbidity worldwide. The role of human papillomaviruses under these conditions has been well studied yet remains inconclusive. One PCR-based study has investigated bacteria in the etiology of these conditions; however, no study has investigated the microbiomes of AK and SCC more broadly. We longitudinally profiled the microbiomes of 112 AK lesions, profiled cross sections of 32 spontaneously arising SCC lesions, and compared these to matching nonlesional photodamaged control skin sites. We identified commonly occurring strains of Propionibacterium and Malassezia at higher relative abundances on nonlesional skin than in AK and SCC lesions, and strains of Staphylococcus

The advent of high throughput technologies has led to a wealth of publicly available 'omics data coming from different sources, such as transcriptomics, proteomics, metabolomics. Combining such large-scale biological data sets can lead to the discovery of important biological insights, provided that relevant information can be extracted in a holistic manner. Current statistical approaches have been focusing on identifying small subsets of molecules (a 'molecular signature') to explain or predict biological conditions, but mainly for a single type of 'omics. In addition, commonly used methods are univariate and consider each biological feature independently. We introduce mixOmics, an R package dedicated to the multivariate analysis of biological data sets with a specific focus on data exploration, dimension reduction and visualisation. By adopting a systems biology approach, the toolkit provides a wide range of methods that statistically integrate several data sets at once to probe relationships between heterogeneous 'omics data sets. Our recent methods extend Projection to Latent Structure (PLS) models for discriminant analysis, for data integration across multiple 'omics data or across independent studies, and for the identification of molecular signatures. We illustrate our latest mixOmics integrative frameworks for the multivariate analyses of 'omics data available from the package.

BACKGROUND: Molecular signatures identified from high-throughput transcriptomic studies often have poor reliability and fail to reproduce across studies. One solution is to combine independent studies into a single integrative analysis, additionally increasing sample size. However, the different protocols and technological platforms across transcriptomic studies produce unwanted systematic variation that strongly confounds the integrative analysis results. When studies aim to discriminate an outcome of interest, the common approach is a sequential two-step procedure; unwanted systematic variation removal techniques are applied prior to classification methods. RESULTS: To limit the risk of overfitting and over-optimistic results of a two-step procedure, we developed a novel multivariate integration method, MINT, that simultaneously accounts for unwanted systematic variation and identifies predictive gene signatures with greater reproducibility and accuracy. In two biological examples on the classification of three human cell types and four subtypes of breast cancer, we combined high-dimensional microarray and RNA-seq data sets and MINT identified highly reproducible and relevant gene signatures predictive of a given phenotype. MINT led to superior classification and prediction accuracy compared to the existing sequential two-step procedures. CONCLUSIONS: MINT is a powerful approach and the first of its kind to solve the integrative classification framework in a single step by combining multiple independent studies. MINT is computationally fast as part of the mixOmics R CRAN package, available at http://www.mixOmics.org/mixMINT/ and http://cran.r-project.org/web/packages/mixOmics/ .

Dynamic changes in biological systems can be captured by measuring molecular expression from different levels (e.g., genes and proteins) across time. Integration of such data aims to identify molecules that show similar expression changes over time; such molecules may be co-regulated and thus involved in similar biological processes. Combining data sources presents a systematic approach to study molecular behaviour. It can compensate for missing data in one source, and can reduce false positives when multiple sources highlight the same pathways. However, integrative approaches must accommodate the challenges inherent in 'omics' data, including high-dimensionality, noise, and timing differences in expression. As current methods for identification of co-expression cannot cope with this level of complexity, we developed a novel algorithm called DynOmics. DynOmics is based on the fast Fourier transform, from which the difference in expression initiation between trajectories can be estimated. This delay can then be used to realign the trajectories and identify those which show a high degree of correlation. Through extensive simulations, we demonstrate that DynOmics is efficient and accurate compared to existing approaches. We consider two case studies highlighting its application, identifying regulatory relationships across 'omics' data within an organism and for comparative gene expression analysis across organisms.

Mesenchymal stromal cells (MSC) are widely used for the study of mesenchymal tissue repair, and increasingly adopted for cell therapy, despite the lack of consensus on the identity of these cells. In part this is due to the lack of specificity of MSC markers. Distinguishing MSC from other stromal cells such as fibroblasts is particularly difficult using standard analysis of surface proteins, and there is an urgent need for improved classification approaches. Transcriptome profiling is commonly used to describe and compare different cell types; however, efforts to identify specific markers of rare cellular subsets may be confounded by the small sample sizes of most studies. Consequently, it is difficult to derive reproducible, and therefore useful markers. We addressed the question of MSC classification with a large integrative analysis of many public MSC datasets. We derived a sparse classifier (The Rohart MSC test) that accurately distinguished MSC from non-MSC samples with >97% accuracy on an internal training set of 635 samples from 41 studies derived on 10 different microarray platforms. The classifier was validated on an external test set of 1,291 samples from 65 studies derived on 15 different platforms, with >95% accuracy. The genes that contribute to the MSC classifier formed a protein-interaction network that included known MSC markers. Further evidence of the relevance of this new MSC panel came from the high number of Mendelian disorders associated with mutations in more than 65% of the network. These result in mesenchymal defects, particularly impacting on skeletal growth and function. The Rohart MSC test is a simple in silico test that accurately discriminates MSC from fibroblasts, other adult stem/progenitor cell types or differentiated stromal cells. It has been implemented in the www.stemformatics.org resource, to assist researchers wishing to benchmark their own MSC datasets or data from the public domain. The code is available from the CRAN reposit

DNA repair genes and pathways that are transcriptionally dysregulated in cancer provide the first line of evidence for the altered DNA repair status in tumours, and hence have been explored intensively as a source for biomarker discovery. The molecular mechanisms underlying DNA repair dysregulation, however, have not been systematically investigated in any cancer type. In this study, we performed a statistical analysis to dissect the roles of DNA copy number alteration (CNA), DNA methylation (DM) at gene promoter regions and the expression changes of transcription factors (TFs) in the differential expression of individual DNA repair genes in normal versus tumour breast samples. These gene-level results were summarised at pathway level to assess whether different DNA repair pathways are affected in distinct manners. Our results suggest that CNA and expression changes of TFs are major causes of DNA repair dysregulation in breast cancer, and that a subset of the identified TFs may exert global impacts on the dysregulation of multiple repair pathways. Our work hence provides novel insights into DNA repair dysregulation in breast cancer. These insights improve our understanding of the molecular basis of the DNA repair biomarkers identified thus far, and have potential to inform future biomarker discovery.

Time course 'omics' experiments are becoming increasingly important to study system-wide dynamic regulation. Despite their high information content, analysis remains challenging. 'Omics' technologies capture quantitative measurements on tens of thousands of molecules. Therefore, in a time course 'omics' experiment molecules are measured for multiple subjects over multiple time points. This results in a large, high-dimensional dataset, which requires computationally efficient approaches for statistical analysis. Moreover, methods need to be able to handle missing values and various levels of noise. We present a novel, robust and powerful framework to analyze time course 'omics' data that consists of three stages: quality assessment and filtering, profile modelling, and analysis. The first step consists of removing molecules for which expression or abundance is highly variable over time. The second step models each molecular expression profile in a linear mixed model framework which takes into account subject-specific variability. The best model is selected through a serial model selection approach and results in dimension reduction of the time course data. The final step includes two types of analysis of the modelled trajectories, namely, clustering analysis to identify groups of correlated profiles over time, and differential expression analysis to identify profiles which differ over time and/or between treatment groups. Through simulation studies we demonstrate the high sensitivity and specificity of our approach for differential expression analysis. We then illustrate how our framework can bring novel insights on two time course 'omics' studies in breast cancer and kidney rejection. The methods are publicly available, implemented in the R CRAN package lmms.

BACKGROUND: A key question when analyzing high throughput data is whether the information provided by the measured biological entities (gene, metabolite expression for example) is related to the experimental conditions, or, rather, to some interfering signals, such as experimental bias or artefacts. Visualization tools are therefore useful to better understand the underlying structure of the data in a 'blind' (unsupervised) way. A well-established technique to do so is Principal Component Analysis (PCA). PCA is particularly powerful if the biological question is related to the highest variance. Independent Component Analysis (ICA) has been proposed as an alternative to PCA as it optimizes an independence condition to give more meaningful components. However, neither PCA nor ICA can overcome both the high dimensionality and noisy characteristics of biological data. RESULTS: We propose Independent Principal Component Analysis (IPCA) that combines the advantages of both PCA and ICA. It uses ICA as a denoising process of the loading vectors produced by PCA to better highlight the important biological entities and reveal insightful patterns in the data. The result is a better clustering of the biological samples on graphical representations. In addition, a sparse version is proposed that performs an internal variable selection to identify biologically relevant features (sIPCA). CONCLUSIONS: On simulation studies and real data sets, we showed that IPCA offers a better visualization of the data than ICA and with a smaller number of components than PCA. Furthermore, a preliminary investigation of the list of genes selected with sIPCA demonstrate that the approach is well able to highlight relevant genes in the data with respect to the biological experiment.IPCA and sIPCA are both implemented in the R package mixomics dedicated to the analysis and exploration of high dimensional biological data sets, and on mixomics' web-interface.

Somatic cell nuclear transfer (SCNT) is the most efficient cell reprogramming technique available, especially when working with bovine species. Although SCNT blastocysts performed equally well or better than controls in the weeks following embryo transfer at Day 7, elongation and gastrulation defects were observed prior to implantation. To understand the developmental implications of embryonic/extra-embryonic interactions, the morphological and molecular features of elongating and gastrulating tissues were analysed. At Day 18, 30 SCNT conceptuses were compared to 20 controls (AI and IVP: 10 conceptuses each); one-half of the SCNT conceptuses appeared normal while the other half showed signs of atypical elongation and gastrulation. SCNT was also associated with a high incidence of discordance in embryonic and extra-embryonic patterns, as evidenced by morphological and molecular "uncoupling". Elongation appeared to be secondarily affected; only 3 of 30 conceptuses had abnormally elongated shapes and there were very few differences in gene expression when they were compared to the controls. However, some of these differences could be linked to defects in microvilli formation or extracellular matrix composition and could thus impact extra-embryonic functions. In contrast to elongation, gastrulation stages included embryonic defects that likely affected the hypoblast, the epiblast, or the early stages of their differentiation. When taking into account SCNT conceptus somatic origin, i.e. the reprogramming efficiency of each bovine ear fibroblast (Low: 0029, Med: 7711, High: 5538), we found that embryonic abnormalities or severe embryonic/extra-embryonic uncoupling were more tightly correlated to embryo loss at implantation than were elongation defects. Alternatively, extra-embryonic differences between SCNT and control conceptuses at Day 18 were related to molecular plasticity (high efficiency/high plasticity) and subsequent pregnancy loss. Finally, because it alters re-d

BACKGROUND: High throughput 'omics' experiments are usually designed to compare changes observed between different conditions (or interventions) and to identify biomarkers capable of characterizing each condition. We consider the complex structure of repeated measurements from different assays where different conditions are applied on the same subjects. RESULTS: We propose a two-step analysis combining a multilevel approach and a multivariate approach to reveal separately the effects of conditions within subjects from the biological variation between subjects. The approach is extended to two-factor designs and to the integration of two matched data sets. It allows internal variable selection to highlight genes able to discriminate the net condition effect within subjects. A simulation study was performed to demonstrate the good performance of the multilevel multivariate approach compared to a classical multivariate method. The multilevel multivariate approach outperformed the classical multivariate approach with respect to the classification error rate and the selection of relevant genes. The approach was applied to an HIV-vaccine trial evaluating the response with gene expression and cytokine secretion. The discriminant multilevel analysis selected a relevant subset of genes while the integrative multilevel analysis highlighted clusters of genes and cytokines that were highly correlated across the samples. CONCLUSIONS: Our combined multilevel multivariate approach may help in finding signatures of vaccine effect and allows for a better understanding of immunological mechanisms activated by the intervention. The integrative analysis revealed clusters of genes, that were associated with cytokine secretion. These clusters can be seen as gene signatures to predict future cytokine response. The approach is implemented in the R package mixOmics (http://cran.r-project.org/) with associated tutorials to perform the analysis(a).

BACKGROUND: Each omics platform is now able to generate a large amount of data. Genomics, proteomics, metabolomics, interactomics are compiled at an ever increasing pace and now form a core part of the fundamental systems biology framework. Recently, several integrative approaches have been proposed to extract meaningful information. However, these approaches lack of visualisation outputs to fully unravel the complex associations between different biological entities. RESULTS: The multivariate statistical approaches 'regularized Canonical Correlation Analysis' and 'sparse Partial Least Squares regression' were recently developed to integrate two types of highly dimensional 'omics' data and to select relevant information. Using the results of these methods, we propose to revisit few graphical outputs to better understand the relationships between two 'omics' data and to better visualise the correlation structure between the different biological entities. These graphical outputs include Correlation Circle plots, Relevance Networks and Clustered Image Maps. We demonstrate the usefulness of such graphical outputs on several biological data sets and further assess their biological relevance using gene ontology analysis. CONCLUSIONS: Such graphical outputs are undoubtedly useful to aid the interpretation of these promising integrative analysis tools and will certainly help in addressing fundamental biological questions and understanding systems as a whole. AVAILABILITY: The graphical tools described in this paper are implemented in the freely available R package mixOmics and in its associated web application.

OBJECTIVE: Neonatal adiposity is a well-recognized complication of gestational diabetes mellitus (GDM). This study aimed to identify factors influencing adiposity in male and female infants of women treated for GDM. RESEARCH DESIGN AND METHODS: This was a prospective study of 84 women with GDM. Daily blood glucose levels (BGLs) were retrieved from glucose meters, and overall mean fasting and mean 2-h postprandial BGLs were calculated for each woman. Infant body composition was measured at birth, and regression analysis was used to identify significant predictors of infant body fat separately in male and female infants. RESULTS: Maternal fasting BGL was the major predictor of adiposity in male infants but had little relationship to adiposity in female infants. In male infants, percent fat was increased by 0.44% for each 0.1 mmol/L increase in mean maternal fasting BGL. Maternal BMI was the primary predictor in female infants but had little effect in males. In female infants, percent fat was increased by 0.11% for each 1 kg/m(2) increase in maternal prepregnancy BMI. CONCLUSIONS: Fetal sex may influence the impact that treatment strategies for GDM have on infant adiposity.

MOTIVATION: Microarrays are being increasingly used in cancer research to better characterize and classify tumors by selecting marker genes. However, as very few of these genes have been validated as predictive biomarkers so far, it is mostly conventional clinical and pathological factors that are being used as prognostic indicators of clinical course. Combining clinical data with gene expression data may add valuable information, but it is a challenging task due to their categorical versus continuous characteristics. We have further developed the mixture of experts (ME) methodology, a promising approach to tackle complex non-linear problems. Several variants are proposed in integrative ME as well as the inclusion of various gene selection methods to select a hybrid signature. RESULTS: We show on three cancer studies that prediction accuracy can be improved when combining both types of variables. Furthermore, the selected genes were found to be of high relevance and can be considered as potential biomarkers for the prognostic selection of cancer therapy. AVAILABILITY: Integrative ME is implemented in the R package integrativeME (http://cran.r-project.org/).

BACKGROUND: In the context of systems biology, few sparse approaches have been proposed so far to integrate several data sets. It is however an important and fundamental issue that will be widely encountered in post genomic studies, when simultaneously analyzing transcriptomics, proteomics and metabolomics data using different platforms, so as to understand the mutual interactions between the different data sets. In this high dimensional setting, variable selection is crucial to give interpretable results. We focus on a sparse Partial Least Squares approach (sPLS) to handle two-block data sets, where the relationship between the two types of variables is known to be symmetric. Sparse PLS has been developed either for a regression or a canonical correlation framework and includes a built-in procedure to select variables while integrating data. To illustrate the canonical mode approach, we analyzed the NCI60 data sets, where two different platforms (cDNA and Affymetrix chips) were used to study the transcriptome of sixty cancer cell lines. RESULTS: We compare the results obtained with two other sparse or related canonical correlation approaches: CCA with Elastic Net penalization (CCA-EN) and Co-Inertia Analysis (CIA). The latter does not include a built-in procedure for variable selection and requires a two-step analysis. We stress the lack of statistical criteria to evaluate canonical correlation methods, which makes biological interpretation absolutely necessary to compare the different gene selections. We also propose comprehensive graphical representations of both samples and variables to facilitate the interpretation of the results. CONCLUSION: sPLS and CCA-EN selected highly relevant genes and complementary findings from the two data sets, which enabled a detailed understanding of the molecular characteristics of several groups of cell lines. These two approaches were found to bring similar results, although they highlighted the same phenomenons with a different

MOTIVATION: With the availability of many 'omics' data, such as transcriptomics, proteomics or metabolomics, the integrative or joint analysis of multiple datasets from different technology platforms is becoming crucial to unravel the relationships between different biological functional levels. However, the development of such an analysis is a major computational and technical challenge as most approaches suffer from high data dimensionality. New methodologies need to be developed and validated. RESULTS: integrOmics efficiently performs integrative analyses of two types of 'omics' variables that are measured on the same samples. It includes a regularized version of canonical correlation analysis to enlighten correlations between two datasets, and a sparse version of partial least squares (PLS) regression that includes simultaneous variable selection in both datasets. The usefulness of both approaches has been demonstrated previously and successfully applied in various integrative studies. AVAILABILITY: integrOmics is freely available from http://CRAN.R-project.org/ or from the web site companion (http://math.univ-toulouse.fr/biostat) that provides full documentation and tutorials. CONTACT: k.lecao@uq.edu.au SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The aim of this paper was to describe, and when possible compare, the multivariate methods used by the participants in the EADGENE WP1.4 workshop. The first approach was for class discovery and class prediction using evidence from the data at hand. Several teams used hierarchical clustering (HC) or principal component analysis (PCA) to identify groups of differentially expressed genes with a similar expression pattern over time points and infective agent (E. coli or S. aureus). The main result from these analyses was that HC and PCA were able to separate tissue samples taken at 24 h following E. coli infection from the other samples. The second approach identified groups of differentially co-expressed genes, by identifying clusters of genes highly correlated when animals were infected with E. coli but not correlated more than expected by chance when the infective pathogen was S. aureus. The third approach looked at differential expression of predefined gene sets. Gene sets were defined based on information retrieved from biological databases such as Gene Ontology. Based on these annotation sources the teams used either the GlobalTest or the Fisher exact test to identify differentially expressed gene sets. The main result from these analyses was that gene sets involved in immune defence responses were differentially expressed.

Microarray analyses have become an important tool in animal genomics. While their use is becoming widespread, there is still a lot of ongoing research regarding the analysis of microarray data. In the context of a European Network of Excellence, 31 researchers representing 14 research groups from 10 countries performed and discussed the statistical analyses of real and simulated 2-colour microarray data that were distributed among participants. The real data consisted of 48 microarrays from a disease challenge experiment in dairy cattle, while the simulated data consisted of 10 microarrays from a direct comparison of two treatments (dye-balanced). While there was broader agreement with regards to methods of microarray normalisation and significance testing, there were major differences with regards to quality control. The quality control approaches varied from none, through using statistical weights, to omitting a large number of spots or omitting entire slides. Surprisingly, these very different approaches gave quite similar results when applied to the simulated data, although not all participating groups analysed both real and simulated data. The workshop was very successful in facilitating interaction between scientists with a diverse background but a common interest in microarray analyses.

A large variety of methods has been proposed in the literature for microarray data analysis. The aim of this paper was to present techniques used by the EADGENE (European Animal Disease Genomics Network of Excellence) WP1.4 participants for data quality control, normalisation and statistical methods for the detection of differentially expressed genes in order to provide some more general data analysis guidelines. All the workshop participants were given a real data set obtained in an EADGENE funded microarray study looking at the gene expression changes following artificial infection with two different mastitis causing bacteria: Escherichia coli and Staphylococcus aureus. It was reassuring to see that most of the teams found the same main biological results. In fact, most of the differentially expressed genes were found for infection by E. coli between uninfected and 24 h challenged udder quarters. Very little transcriptional variation was observed for the bacteria S. aureus. Lists of differentially expressed genes found by the different research teams were, however, quite dependent on the method used, especially concerning the data quality control step. These analyses also emphasised a biological problem of cross-talk between infected and uninfected quarters which will have to be dealt with for further microarray studies.