Establishing correlates of maternal-fetal cytomegalovirus transmission-one step closer through predictive modeling.

BACKGROUND: Determinants of maternal-fetal cytomegalovirus (CMV) transmission and factors influencing the severity of congenital CMV (cCMV) infection are not well understood. METHODS: We conducted a descriptive, multi-center study in pregnant women ≥18 years old with primary CMV infection and their newborns (NCT01251744) to explore maternal immune responses to CMV and determine potential immunologic/virologic correlates of cCMV following primary infection during pregnancy. We developed alternative approaches looking into univariate/multivariate factors associated with cCMV, including a participant clustering/stratification approach and an interpretable predictive model-based approach using trained decision trees for risk prediction (post-hoc analyses). RESULTS: Pregnant women were grouped in three distinct clusters with similar baseline characteristics, particularly gestational age at diagnosis. We observed a trend for higher viral loads in urine and saliva samples from mothers of infants with cCMV versus without cCMV. When using a trained predictive-model approach that accounts for interaction effects between variables, anti-pentamer IgG antibody concentration and viral load in saliva were identified as biomarkers jointly associated with the risk of maternal-fetal CMV transmission. CONCLUSION: We identified biomarkers of CMV maternal-fetal transmission. After validation in larger studies, our findings will guide the management of primary infection during pregnancy and the development of vaccines against cCMV.

The human cytomegalovirus (CMV) is common and usually causes no symptoms in healthy individuals. However, CMV infections can be life-threatening in individuals with improperly functioning or immature immune systems, such as fetuses. Women can become infected with CMV for the first time (primary infection) during pregnancy. If CMV is transmitted from mother to fetus before the second trimester, the infant can suffer from severe disorders such as hearing loss and delayed development. We aimed to identify characteristics of pregnant women with a primary CMV infection that may increase the likelihood of transmitting CMV to the fetus. We considered demographical, clinical, and behavioral characteristics, as well as immune responses and the quantity of virus detected in the women's blood, urine, saliva, and vaginal mucus. Because we could not identify one single characteristic that could predict a high risk of CMV transmission, we developed new data analysis models to study how they can be combined. We found that antibodies targeting a pentameric antigen of the virus envelope and the presence of virus in saliva can together predict the risk of CMV transmission from mother to fetus. Our results can help improve the care of CMV-infected pregnant women and the design of CMV vaccines.

Subsequent AS01-adjuvanted vaccinations induce similar transcriptional responses in populations with different disease statuses.

Transcriptional responses to adjuvanted vaccines can vary substantially among populations. Interindividual diversity in levels of pathogen exposure, and thus of cell-mediated immunological memory at baseline, may be an important determinant of population differences in vaccine responses. Adjuvant System AS01 is used in licensed or candidate vaccines for several diseases and populations, yet the impact of pre-existing immunity on its adjuvanticity remains to be elucidated. In this exploratory post-hoc analysis of clinical trial samples (clinicalTrials.gov: NCT01424501), we compared gene expression patterns elicited by two immunizations with the candidate tuberculosis (TB) vaccine M72/AS01, between three groups of individuals with different levels of memory responses to TB antigens before vaccination. Analyzed were one group of TB-disease-treated individuals, and two groups of TB-disease-naïve individuals who were (based on purified protein derivative [PPD] skin-test results) stratified into PPD-positive and PPD-negative groups. Although TB-disease-treated individuals displayed slightly stronger transcriptional responses after each vaccine dose, functional gene signatures were overall not distinctly different between groups. Considering the similarities with the signatures found previously for other AS01-adjuvanted vaccines, many features of the response appeared to be adjuvant-driven. Across groups, cell proliferation-related signals at 7 days post-dose 1 were associated with increased anti-M72 antibody response magnitudes. These early signals were stronger in the TB-disease-treated group as compared to both TB-disease-naïve groups. Interindividual homogeneity in gene expression levels was also higher for TB-disease-treated individuals post-dose 1, but increased in all groups post-dose 2 to attain similar levels between the three groups. Altogether, strong cell-mediated memory responses at baseline accelerated and amplified transcriptional responses to a single dose of this AS01-adjuvanted vaccine, resulting in more homogenous gene expression levels among the highly-primed individuals as compared to the disease-naïve individuals. However, after a second vaccination, response heterogeneity decreased and was similar across groups, irrespective of the degree of immune memory acquired at baseline. This information can support the design and analysis of future clinical trials evaluating AS01-adjuvanted vaccines.

Innate Immune Responses to Chimpanzee Adenovirus Vector 155 Vaccination in Mice and Monkeys.

Replication-deficient chimpanzee adenovirus (ChAd) vectors represent an attractive vaccine platform and are thus employed as vaccine candidates against several infectious diseases. Since inducing effective immunity depends on the interplay between innate and adaptive immunity, a deeper understanding of innate immune responses elicited by intramuscularly injected ChAd vectors in tissues can advance the platform's development. Using different candidate vaccines based on the Group C ChAd type 155 (ChAd155) vector, we characterized early immune responses in injected muscles and draining lymph nodes (dLNs) from mice, and complemented these analyses by evaluating cytokine responses and gene expression patterns in peripheral blood from ChAd155-injected macaques. In mice, vector DNA levels gradually decreased post-immunization, but local transgene mRNA expression exhibited two transient peaks [at 6 h and Day (D)5], which were most obvious in dLNs. This dynamic pattern was mirrored by the innate responses in tissues, which developed as early as 1-3 h (cytokines/chemokines) or D1 (immune cells) post-vaccination. They were characterized by a CCL2- and CXCL9/10-dominated chemokine profile, peaking at 6 h (with CXCL10/CCL2 signals also detectable in serum) and D7, and clear immune-cell infiltration peaks at D1/D2 and D6/D7. Experiments with a green fluorescent protein-expressing ChAd155 vector revealed infiltrating hematopoietic cell subsets at the injection site. Cell infiltrates comprised mostly monocytes in muscles, and NK cells, T cells, dendritic cells, monocytes, and B cells in dLNs. Similar bimodal dynamics were observed in whole-blood gene signatures in macaques: most of the 17 enriched immune/innate signaling pathways were significantly upregulated at D1 and D7 and downregulated at D3, and clustering analysis revealed stronger similarities between D1 and D7 signatures versus the D3 signature. Serum cytokine responses (CXCL10, IL1Ra, and low-level IFN-α) in macaques were predominantly observed at D1. Altogether, the early immune responses exhibited bimodal kinetics with transient peaks at D1/D2 and D6/D7, mostly with an IFN-associated signature, and these features were remarkably consistent across most analyzed parameters in murine tissues and macaque blood. These compelling observations reveal a novel aspect of the dynamics of innate immunity induced by ChAd155-vectored vaccines, and contribute to ongoing research to better understand how adenovectors can promote vaccine-induced immunity.

Transcriptional profiles of adjuvanted hepatitis B vaccines display variable interindividual homogeneity but a shared core signature.

The current routine use of adjuvants in human vaccines provides a strong incentive to increase our understanding of how adjuvants differ in their ability to stimulate innate immunity and consequently enhance vaccine immunogenicity. Here, we evaluated gene expression profiles in cells from whole blood elicited in naive subjects receiving the hepatitis B surface antigen formulated with different adjuvants. We identified a core innate gene signature emerging 1 day after the second vaccination and that was shared by the recipients of vaccines formulated with adjuvant systems AS01B, AS01E, or AS03. This core signature associated with the magnitude of the hepatitis B surface-specific antibody response and was characterized by positive regulation of genes associated with interferon-related responses or the innate cell compartment and by negative regulation of natural killer cell-associated genes. Analysis at the individual subject level revealed that the higher immunogenicity of AS01B-adjuvanted vaccine was linked to its ability to induce this signature in most vaccinees even after the first vaccination. Therefore, our data suggest that adjuvanticity is not strictly defined by the nature of the receptors or signaling pathways it activates but by the ability of the adjuvant to consistently induce a core inflammatory signature across individuals.

Modeling adaptive response profiles in a vaccine clinical trial.

BACKGROUND: Vaccine clinical studies typically provide time-resolved data on adaptive response read-outs in response to the administration of that particular vaccine to a cohort of individuals. However, modeling such data is challenged by the properties of these time-resolved profiles such as non-linearity, scarcity of measurement points, scheduling of the vaccine at multiple time points. Linear Mixed Models (LMM) are often used for the analysis of longitudinal data but their use in these time-resolved immunological data is not common yet. Apart from the modeling challenges mentioned earlier, selection of the optimal model by using information-criterion-based measures is far from being straight-forward. The aim of this study is to provide guidelines for the application and selection of LMMs that deal with the challenging characteristics of the typical data sets in the field of vaccine clinical studies. METHODS: We used antibody measurements in response to Hepatitis-B vaccine with five different adjuvant formulations for demonstration purposes. We built piecewise-linear, piecewise-quadratic and cubic models with transformations of the axes with pre-selected or optimized knot locations where time is a numerical variable. We also investigated models where time is categorical and random effects are shared intercepts between different measurement points. We compared all models by using Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), Deviance Information Criterion (DIC), variations of conditional AIC and by visual inspection of the model fit in the light of prior biological information. RESULTS: There are various ways of dealing with the challenges of the data which have their own advantages and disadvantages. We explain these in detail here. Traditional information-criteria-based measures work well for the coarse selection of the model structure and complexity, however are not efficient at fine tuning of the complexity level of the random effects. CONCLUSIONS: We show that common statistical measures for optimal model complexity are not sufficient. Rather, explicitly accounting for model purpose and biological interpretation is needed to arrive at relevant models. TRIAL REGISTRATION: Clinical trial registration number for this study: NCT00805389, date of registration: December 9, 2008 (pro-active registration).

The Ratiometric Transcript Signature MX2/GPR183 Is Consistently Associated With RTS,S-Mediated Protection Against Controlled Human Malaria Infection.

The RTS,S/AS01 vaccine provides partial protection against Plasmodium falciparum infection but determinants of protection and/or disease are unclear. Previously, anti-circumsporozoite protein (CSP) antibody titers and blood RNA signatures were associated with RTS,S/AS01 efficacy against controlled human malaria infection (CHMI). By analyzing host blood transcriptomes from five RTS,S vaccination CHMI studies, we demonstrate that the transcript ratio MX2/GPR183, measured 1 day after third immunization, discriminates protected from non-protected individuals. This ratiometric signature provides information that is complementary to anti-CSP titer levels for identifying RTS,S/AS01 immunized people who developed protective immunity and suggests a role for interferon and oxysterol signaling in the RTS,S mode of action.

Identification of potential biomarkers of vaccine inflammation in mice.

Systems vaccinology approaches have been used successfully to define early signatures of the vaccine-induced immune response. However, the possibility that transcriptomics can also identify a correlate or surrogate for vaccine inflammation has not been fully explored. We have compared four licensed vaccines with known safety profiles, as well as three agonists of Toll-like receptors (TLRs) with known inflammatory potential, to elucidate the transcriptomic profile of an acceptable response to vaccination versus that of an inflammatory reaction. In mice, we looked at the transcriptomic changes in muscle at the injection site, the lymph node that drained the muscle, and the peripheral blood mononuclear cells (PBMCs)isolated from the circulating blood from 4 hr after injection and over the next week. A detailed examination and comparative analysis of these transcriptomes revealed a set of novel biomarkers that are reflective of inflammation after vaccination. These biomarkers are readily measurable in the peripheral blood, providing useful surrogates of inflammation, and provide a way to select candidates with acceptable safety profiles.

Adjuvant-Associated Peripheral Blood mRNA Profiles and Kinetics Induced by the Adjuvanted Recombinant Protein Candidate Tuberculosis Vaccine M72/AS01 in Bacillus Calmette-Guérin-Vaccinated Adults.

Systems biology has the potential to identify gene signatures associated with vaccine immunogenicity and protective efficacy. The main objective of this study was to identify optimal postvaccination time points for evaluating peripheral blood RNA expression profiles in relation to vaccine immunogenicity and potential efficacy in recipients of the candidate tuberculosis vaccine M72/AS01. In this phase II open-label study (NCT01669096; https://clinicaltrials.gov/), healthy Bacillus Calmette-Guérin-primed, HIV-negative adults were administered two doses (30 days apart) of M72/AS01. Twenty subjects completed the study and 18 subjects received two doses. Blood samples were collected pre-dose 1, pre-dose 2, and 1, 7, 10, 14, 17, and 30 days post-dose 2. RNA expression in whole blood (WB) and peripheral blood mononuclear cells (PBMCs) was quantified using microarray technology. Serum interferon-gamma responses and M72-specific CD4+ T cell responses to vaccination, and the observed safety profile were similar to previous trials. Two different approaches were utilized to analyze the RNA expression data. First, a kinetic analysis of RNA expression changes using blood transcription modules revealed early (1 day post-dose 2) activation of several pathways related to innate immune activation, both in WB and PBMC. Second, using a previously identified gene signature as a classifier, optimal postvaccination time points were identified. Since M72/AS01 efficacy remains to be established, a PBMC-derived gene signature associated with the protective efficacy of a similarly adjuvanted candidate malaria vaccine was used as a proxy for this purpose. This approach was based on the assumption that the AS01 adjuvant used in both studies could induce shared innate immune pathways. Subjects were classified as gene signature positive (GS+) or gene signature negative (GS-). Assignments of subjects to GS+ or GS- groups were confirmed by significant differences in RNA expression of the gene signature genes in PBMCs at 14 days post-dose 2 relative to prevaccination and in WB samples at 7, 10, 14, and 17 days post-dose 2 relative to prevaccination. Hence, in comparison with a prevaccination, 7, 10, 14, and 17 days postvaccination appeared to be suitable time points for identifying potentially clinically relevant transcriptome responses to M72/AS01 in WB samples.

Erratum: Author Correction: Cellular and molecular synergy in AS01-adjuvanted vaccines results in an early IFNγ response promoting vaccine immunogenicity.

[This corrects the article DOI: 10.1038/s41541-017-0027-3.].

Different Adjuvants Induce Common Innate Pathways That Are Associated with Enhanced Adaptive Responses against a Model Antigen in Humans.

To elucidate the role of innate responses in vaccine immunogenicity, we compared early responses to hepatitis B virus (HBV) surface antigen (HBsAg) combined with different Adjuvant Systems (AS) in healthy HBV-naïve adults, and included these parameters in multi-parametric models of adaptive responses. A total of 291 participants aged 18-45 years were randomized 1:1:1:1:1 to receive HBsAg with AS01B, AS01E, AS03, AS04, or Alum/Al(OH)3 at days 0 and 30 (ClinicalTrials.gov: NCT00805389). Blood protein, cellular, and mRNA innate responses were assessed at early time-points and up to 7 days after vaccination, and used with reactogenicity symptoms in linear regression analyses evaluating their correlation with HBs-specific CD4+ T-cell and antibody responses at day 44. All AS induced transient innate responses, including interleukin (IL)-6 and C-reactive protein (CRP), mostly peaking at 24 h post-vaccination and subsiding to baseline within 1-3 days. After the second but not the first injection, median interferon (IFN)-γ levels were increased in the AS01B group, and IFN-γ-inducible protein-10 levels and IFN-inducible genes upregulated in the AS01 and AS03 groups. No distinct marker or signature was specific to one particular AS. Innate profiles were comparable between AS01B, AS01E, and AS03 groups, and between AS04 and Alum groups. AS group rankings within adaptive and innate response levels and reactogenicity prevalence were similar (AS01B ≥ AS01E > AS03 > AS04 > Alum), suggesting an association between magnitudes of inflammatory and vaccine responses. Modeling revealed associations between adaptive responses and specific traits of the innate response post-dose 2 (activation of the IFN-signaling pathway, CRP and IL-6 responses). In conclusion, the ability of AS01 and AS03 to enhance adaptive responses to co-administered HBsAg is likely linked to their capacity to activate innate immunity, particularly the IFN-signaling pathway.

Predicting RTS,S Vaccine-Mediated Protection from Transcriptomes in a Malaria-Challenge Clinical Trial.

The RTS,S candidate malaria vaccine can protect against controlled human malaria infection (CHMI), but how protection is achieved remains unclear. Here, we have analyzed longitudinal peripheral blood transcriptome and immunogenicity data from a clinical efficacy trial in which healthy adults received three RTS,S doses 4 weeks apart followed by CHMI 2 weeks later. Multiway partial least squares discriminant analysis (N-PLS-DA) of transcriptome data identified 110 genes that could be used in predictive models of protection. Among the 110 genes, 42 had known immune-related functions, including 29 that were related to the NF-κB-signaling pathway and 14 to the IFN-γ-signaling pathway. Post-dose 3 serum IFN-γ concentrations were also correlated with protection; and N-PLS-DA of IFN-γ-signaling pathway transcriptome data selected almost all (44/45) of the representative genes for predictive models of protection. Hence, the identification of the NF-κB and IFN-γ pathways provides further insight into how vaccine-mediated protection may be achieved.

Systems analysis of protective immune responses to RTS,S malaria vaccination in humans.

RTS,S is an advanced malaria vaccine candidate and confers significant protection against Plasmodium falciparum infection in humans. Little is known about the molecular mechanisms driving vaccine immunity. Here, we applied a systems biology approach to study immune responses in subjects receiving three consecutive immunizations with RTS,S (RRR), or in those receiving two immunizations of RTS,S/AS01 following a primary immunization with adenovirus 35 (Ad35) (ARR) vector expressing circumsporozoite protein. Subsequent controlled human malaria challenge (CHMI) of the vaccinees with Plasmodium-infected mosquitoes, 3 wk after the final immunization, resulted in ∼50% protection in both groups of vaccinees. Circumsporozoite protein (CSP)-specific antibody titers, prechallenge, were associated with protection in the RRR group. In contrast, ARR-induced lower antibody responses, and protection was associated with polyfunctional CD4+ T-cell responses 2 wk after priming with Ad35. Molecular signatures of B and plasma cells detected in PBMCs were highly correlated with antibody titers prechallenge and protection in the RRR cohort. In contrast, early signatures of innate immunity and dendritic cell activation were highly associated with protection in the ARR cohort. For both vaccine regimens, natural killer (NK) cell signatures negatively correlated with and predicted protection. These results suggest that protective immunity against P. falciparum can be achieved via multiple mechanisms and highlight the utility of systems approaches in defining molecular correlates of protection to vaccination.

Not Just a Sum? Identifying Different Types of Interplay between Constituents in Combined Interventions.

MOTIVATION: Experiments in which the effect of combined manipulations is compared with the effects of their pure constituents have received a great deal of attention. Examples include the study of combination therapies and the comparison of double and single knockout model organisms. Often the effect of the combined manipulation is not a mere addition of the effects of its constituents, with quite different forms of interplay between the constituents being possible. Yet, a well-formalized taxonomy of possible forms of interplay is lacking, let alone a statistical methodology to test for their presence in empirical data. RESULTS: Starting from a taxonomy of a broad range of forms of interplay between constituents of a combined manipulation, we propose a sound statistical hypothesis testing framework to test for the presence of each particular form of interplay. We illustrate the framework with analyses of public gene expression data on the combined treatment of dendritic cells with curdlan and GM-CSF and show that these lead to valuable insights into the mode of action of the constituent treatments and their combination. AVAILABILITY AND IMPLEMENTATION: R code implementing the statistical testing procedure for microarray gene expression data is available as supplementary material. The data are available from the Gene Expression Omnibus with accession number GSE32986.

Characterization of Clostridium thermocellum strains with disrupted fermentation end-product pathways.

Clostridium thermocellum is a thermophilic, cellulolytic anaerobe that is a candidate microorganism for industrial biofuels production. Strains with mutations in genes associated with production of L-lactate (Δldh) and/or acetate (Δpta) were characterized to gain insight into the intracellular processes that convert cellobiose to ethanol and other fermentation end-products. Cellobiose-grown cultures of the Δldh strain had identical biomass accumulation, fermentation end-products, transcription profile, and intracellular metabolite concentrations compared to its parent strain (DSM1313 Δhpt Δspo0A). The Δpta-deficient strain grew slower and had 30 % lower final biomass concentration compared to the parent strain, yet produced 75 % more ethanol. A Δldh Δpta double-mutant strain evolved for faster growth had a growth rate and ethanol yield comparable to the parent strain, whereas its biomass accumulation was comparable to Δpta. Free amino acids were secreted by all examined strains, with both Δpta strains secreting higher amounts of alanine, valine, isoleucine, proline, glutamine, and threonine. Valine concentration for Δldh Δpta reached 5 mM by the end of growth, or 2.7 % of the substrate carbon utilized. These secreted amino acid concentrations correlate with increased intracellular pyruvate concentrations, up to sixfold in the Δpta and 16-fold in the Δldh Δpta strain. We hypothesize that the deletions in fermentation end-product pathways result in an intracellular redox imbalance, which the organism attempts to relieve, in part by recycling NADP⁺ through increased production of amino acids.

A flexible framework for sparse simultaneous component based data integration.

BACKGROUND: High throughput data are complex and methods that reveal structure underlying the data are most useful. Principal component analysis, frequently implemented as a singular value decomposition, is a popular technique in this respect. Nowadays often the challenge is to reveal structure in several sources of information (e.g., transcriptomics, proteomics) that are available for the same biological entities under study. Simultaneous component methods are most promising in this respect. However, the interpretation of the principal and simultaneous components is often daunting because contributions of each of the biomolecules (transcripts, proteins) have to be taken into account. RESULTS: We propose a sparse simultaneous component method that makes many of the parameters redundant by shrinking them to zero. It includes principal component analysis, sparse principal component analysis, and ordinary simultaneous component analysis as special cases. Several penalties can be tuned that account in different ways for the block structure present in the integrated data. This yields known sparse approaches as the lasso, the ridge penalty, the elastic net, the group lasso, sparse group lasso, and elitist lasso. In addition, the algorithmic results can be easily transposed to the context of regression. Metabolomics data obtained with two measurement platforms for the same set of Escherichia coli samples are used to illustrate the proposed methodology and the properties of different penalties with respect to sparseness across and within data blocks. CONCLUSION: Sparse simultaneous component analysis is a useful method for data integration: First, simultaneous analyses of multiple blocks offer advantages over sequential and separate analyses and second, interpretation of the results is highly facilitated by their sparseness. The approach offered is flexible and allows to take the block structure in different ways into account. As such, structures can be found that are exclusively tied to one data platform (group lasso approach) as well as structures that involve all data platforms (Elitist lasso approach). AVAILABILITY: The additional file contains a MATLAB implementation of the sparse simultaneous component method.

Simultaneous analysis of coupled data matrices subject to different amounts of noise.

In many areas of science, research questions imply the analysis of a set of coupled data blocks, with, for instance, each block being an experimental unit by variable matrix, and the variables being the same in all matrices. To obtain an overall picture of the mechanisms that play a role in the different data matrices, the information in these matrices needs to be integrated. This may be achieved by applying a data-analytic strategy in which a global model is fitted to all data matrices simultaneously, as in some forms of simultaneous component analysis (SCA). Since such a strategy implies that all data entries, regardless the matrix they belong to, contribute equally to the analysis, it may obfuscate the overall picture of the mechanisms underlying the data when the different data matrices are subject to different amounts of noise. One way out is to downweight entries from noisy data matrices in favour of entries from less noisy matrices. Information regarding the amount of noise that is present in each matrix, however, is, in most cases, not available. To deal with these problems, in this paper a novel maximum-likelihood-based simultaneous component analysis method, referred to as MxLSCA, is proposed. Being a stochastic extension of SCA, in MxLSCA the amount of noise in each data matrix is estimated and entries from noisy data matrices are downweighted. Both in an extensive simulation study and in an application to data stemming from cross-cultural emotion psychology, it is shown that the novel MxLSCA strategy outperforms the SCA strategy with respect to disclosing the mechanisms underlying the coupled data.

Illness and symptom perception: a theoretical approach towards an integrative measurement model.

Several models have been proposed to conceptualize psychological representations of health, illness, and bodily sensations. These models differ as to the cognitive and affective components they include, whether they study the interaction of these components, and whether associations between psychological representations of bodily states and affective and behavioral reactions to these representations are considered conditional. These different conceptualizations and corresponding measurement approaches exist in parallel without resulting in synergistic effects or theoretical advancements within the field. In this paper, we review theoretical models on perception and attitudes and construct an integrative theoretical framework on psychological representation of bodily symptoms as well as more abstract representations of health and disease. The aim of this combination of approaches is to unify the strengths of different research domains in the conceptualization and measurement of mental representations of bodily states. Furthermore, the aim is to specify new, testable predictions and implications about the (conditional) relationship of these mental representations and affective and behavioral consequences. A core element in this integrative model is comparison. We review how comparison processes can change the cognitive and affective reference frame for illness and symptom perception and in turn affective and behavioral reactions. We discuss implications for measurement of illness and symptom representations as well as implications for clinical practice. Finally, we make suggestions for a research agenda to validate the proposed model as well as to address new questions derived from it.

Identification of modules in Aspergillus niger by gene co-expression network analysis.

The fungus Aspergillus niger has been studied in considerable detail with respect to various industrial applications. Although its central metabolic pathways are established relatively well, the mechanisms that control the adaptation of its metabolism are understood rather poorly. In this study, clustering of co-expressed genes has been performed on the basis of DNA microarray data sets from two experimental approaches. In one approach, low amounts of inducer caused a relatively mild perturbation, while in the other approach the imposed environmental conditions including carbon source starvation caused severe perturbed stress. A set of conserved genes was used to construct gene co-expression networks for both the individual and combined data sets. Comparative analysis revealed the existence of modules, some of which are present in all three networks. In addition, experimental condition-specific modules were identified. Module-derived consensus expression profiles enabled the integration of all protein-coding A. niger genes to the co-expression analysis, including hypothetical and poorly conserved genes. Conserved sequence motifs were detected in the upstream region of genes that cluster in some modules, e.g., the binding site for the amino acid metabolism-related transcription factor CpcA as well as for the fatty acid metabolism-related transcription factors, FarA and FarB. Moreover, not previously described putative transcription factor binding sites were discovered for two modules: the motif 5'-CGACAA is overrepresented in the module containing genes encoding cytosolic ribosomal proteins, while the motif 5'-GGCCGCG is overrepresented in genes related to 'gene expression', such as RNA helicases and translation initiation factors.

Integrating functional genomics data using maximum likelihood based simultaneous component analysis.

BACKGROUND: In contemporary biology, complex biological processes are increasingly studied by collecting and analyzing measurements of the same entities that are collected with different analytical platforms. Such data comprise a number of data blocks that are coupled via a common mode. The goal of collecting this type of data is to discover biological mechanisms that underlie the behavior of the variables in the different data blocks. The simultaneous component analysis (SCA) family of data analysis methods is suited for this task. However, a SCA may be hampered by the data blocks being subjected to different amounts of measurement error, or noise. To unveil the true mechanisms underlying the data, it could be fruitful to take noise heterogeneity into consideration in the data analysis. Maximum likelihood based SCA (MxLSCA-P) was developed for this purpose. In a previous simulation study it outperformed normal SCA-P. This previous study, however, did not mimic in many respects typical functional genomics data sets, such as, data blocks coupled via the experimental mode, more variables than experimental units, and medium to high correlations between variables. Here, we present a new simulation study in which the usefulness of MxLSCA-P compared to ordinary SCA-P is evaluated within a typical functional genomics setting. Subsequently, the performance of the two methods is evaluated by analysis of a real life Escherichia coli metabolomics data set. RESULTS: In the simulation study, MxLSCA-P outperforms SCA-P in terms of recovery of the true underlying scores of the common mode and of the true values underlying the data entries. MxLSCA-P further performed especially better when the simulated data blocks were subject to different noise levels. In the analysis of an E. coli metabolomics data set, MxLSCA-P provided a slightly better and more consistent interpretation. CONCLUSION: MxLSCA-P is a promising addition to the SCA family. The analysis of coupled functional genomics data blocks could benefit from its ability to take different noise levels per data block into consideration and improve the recovery of the true patterns underlying the data. Moreover, the maximum likelihood based approach underlying MxLSCA-P could be extended to custom-made solutions to specific problems encountered.

Metabolomics data exploration guided by prior knowledge.

In metabolomics research, it is often important to focus the data analysis to specific areas of interest within the metabolome. In this paper, we describe the application of consensus principal component analysis (CPCA) and canonical correlation analysis (CCA) as a means to explore the relation between metabolome data and (i) biochemically related metabolites and (ii) an amino acid biosynthesis pathway. CPCA searches for major trends in the behavior of metabolite concentrations that are in common for the metabolites of interest and the remainder of the metabolome. CCA identifies the strongest correlations between the metabolites of interest and the remainder of the metabolome. CPCA and CCA were applied to two different microbial metabolomics data sets. The first data set, derived from Pseudomonas putida S12, was relatively simple as it contained metabolomes obtained under four environmental conditions only. The second data set, obtained from Escherichia coli, was much more complex as it consisted of metabolomes obtained under 28 different environmental conditions. In case of the simple and coherent P. putida S12 data set, CCA and CPCA gave similar results as the variation in the subset of the selected metabolites and the remainder of the metabolome was similar. In contrast, CCA and CPCA yielded different results in case of the E. coli data set. With CPCA the trends in the selected subset--the phenylalanine biosynthesis pathway--dominated the results. The main trends were related to high and low phenylalanine productivity, and the metabolites showing a similar behavior in concentration were metabolites regulating the phenylalanine biosynthesis route in the subset and metabolites related to general amino acid metabolism in the remainder of the metabolome. With CCA, neither subset truly dominated the data analysis. CCA described the differences between the wild type and the overproducing strain and the differences between the succinate and glucose grown cells. For the difference between the wild type and the overproducing strain, metabolites from the beginning and the end of aromatic amino acid pathways like erythrose-4-phosphate, tryptophan, and phenylalanine were important for the selected metabolites. CCA and CPCA proved to be complementary data analysis tools that enable the focusing of the data analysis on groups of metabolites that are of specific interest in relation to the remainder of the metabolome. Compared to an ordinary PCA, focusing the data analysis on biologically relevant metabolites lead especially for the complex E. coli data to a better biological interpretation of the data.