SUBATOMIC: a SUbgraph BAsed mulTi-OMIcs clustering framework to analyze integrated multi-edge networks.

BACKGROUND: Representing the complex interplay between different types of biomolecules across different omics layers in multi-omics networks bears great potential to gain a deep mechanistic understanding of gene regulation and disease. However, multi-omics networks easily grow into giant hairball structures that hamper biological interpretation. Module detection methods can decompose these networks into smaller interpretable modules. However, these methods are not adapted to deal with multi-omics data nor consider topological features. When deriving very large modules or ignoring the broader network context, interpretability remains limited. To address these issues, we developed a SUbgraph BAsed mulTi-OMIcs Clustering framework (SUBATOMIC), which infers small and interpretable modules with a specific topology while keeping track of connections to other modules and regulators. RESULTS: SUBATOMIC groups specific molecular interactions in composite network subgraphs of two and three nodes and clusters them into topological modules. These are functionally annotated, visualized and overlaid with expression profiles to go from static to dynamic modules. To preserve the larger network context, SUBATOMIC investigates statistically the connections in between modules as well as between modules and regulators such as miRNAs and transcription factors. We applied SUBATOMIC to analyze a composite Homo sapiens network containing transcription factor-target gene, miRNA-target gene, protein-protein, homologous and co-functional interactions from different databases. We derived and annotated 5586 modules with diverse topological, functional and regulatory properties. We created novel functional hypotheses for unannotated genes. Furthermore, we integrated modules with condition specific expression data to study the influence of hypoxia in three cancer cell lines. We developed two prioritization strategies to identify the most relevant modules in specific biological contexts: one considering GO term enrichments and one calculating an activity score reflecting the degree of differential expression. Both strategies yielded modules specifically reacting to low oxygen levels. CONCLUSIONS: We developed the SUBATOMIC framework that generates interpretable modules from integrated multi-omics networks and applied it to hypoxia in cancer. SUBATOMIC can infer and contextualize modules, explore condition or disease specific modules, identify regulators and functionally related modules, and derive novel gene functions for uncharacterized genes. The software is available at https://github.com/CBIGR/SUBATOMIC .

Nearby transposable elements impact plant stress gene regulatory networks: a meta-analysis in A. thaliana and S. lycopersicum.

BACKGROUND: Transposable elements (TE) make up a large portion of many plant genomes and are playing innovative roles in genome evolution. Several TEs can contribute to gene regulation by influencing expression of nearby genes as stress-responsive regulatory motifs. To delineate TE-mediated plant stress regulatory networks, we took a 2-step computational approach consisting of identifying TEs in the proximity of stress-responsive genes, followed by searching for cis-regulatory motifs in these TE sequences and linking them to known regulatory factors. Through a systematic meta-analysis of RNA-seq expression profiles and genome annotations, we investigated the relation between the presence of TE superfamilies upstream, downstream or within introns of nearby genes and the differential expression of these genes in various stress conditions in the TE-poor Arabidopsis thaliana and the TE-rich Solanum lycopersicum. RESULTS: We found that stress conditions frequently expressed genes having members of various TE superfamilies in their genomic proximity, such as SINE upon proteotoxic stress and Copia and Gypsy upon heat stress in A. thaliana, and EPRV and hAT upon infection, and Harbinger, LINE and Retrotransposon upon light stress in S. lycopersicum. These stress-specific gene-proximal TEs were mostly located within introns and more detected near upregulated than downregulated genes. Similar stress conditions were often related to the same TE superfamily. Additionally, we detected both novel and known motifs in the sequences of those TEs pointing to regulatory cooption of these TEs upon stress. Next, we constructed the regulatory network of TFs that act through binding these TEs to their target genes upon stress and discovered TE-mediated regulons targeted by TFs such as BRB/BPC, HD, HSF, GATA, NAC, DREB/CBF and MYB factors in Arabidopsis and AP2/ERF/B3, NAC, NF-Y, MYB, CXC and HD factors in tomato. CONCLUSIONS: Overall, we map TE-mediated plant stress regulatory networks using numerous stress expression profile studies for two contrasting plant species to study the regulatory role TEs play in the response to stress. As TE-mediated gene regulation allows plants to adapt more rapidly to new environmental conditions, this study contributes to the future development of climate-resilient plants.

Function, dynamics and evolution of network motif modules in integrated gene regulatory networks of worm and plant.

Gene regulatory networks (GRNs) consist of different molecular interactions that closely work together to establish proper gene expression in time and space. Especially in higher eukaryotes, many questions remain on how these interactions collectively coordinate gene regulation. We study high quality GRNs consisting of undirected protein-protein, genetic and homologous interactions, and directed protein-DNA, regulatory and miRNA-mRNA interactions in the worm Caenorhabditis elegans and the plant Arabidopsis thaliana. Our data-integration framework integrates interactions in composite network motifs, clusters these in biologically relevant, higher-order topological network motif modules, overlays these with gene expression profiles and discovers novel connections between modules and regulators. Similar modules exist in the integrated GRNs of worm and plant. We show how experimental or computational methodologies underlying a certain data type impact network topology. Through phylogenetic decomposition, we found that proteins of worm and plant tend to functionally interact with proteins of a similar age, while at the regulatory level TFs favor same age, but also older target genes. Despite some influence of the duplication mode difference, we also observe at the motif and module level for both species a preference for age homogeneity for undirected and age heterogeneity for directed interactions. This leads to a model where novel genes are added together to the GRNs in a specific biological functional context, regulated by one or more TFs that also target older genes in the GRNs. Overall, we detected topological, functional and evolutionary properties of GRNs that are potentially universal in all species.

Arabidopsis ensemble reverse-engineered gene regulatory network discloses interconnected transcription factors in oxidative stress.

The abiotic stress response in plants is complex and tightly controlled by gene regulation. We present an abiotic stress gene regulatory network of 200,014 interactions for 11,938 target genes by integrating four complementary reverse-engineering solutions through average rank aggregation on an Arabidopsis thaliana microarray expression compendium. This ensemble performed the most robustly in benchmarking and greatly expands upon the availability of interactions currently reported. Besides recovering 1182 known regulatory interactions, cis-regulatory motifs and coherent functionalities of target genes corresponded with the predicted transcription factors. We provide a valuable resource of 572 abiotic stress modules of coregulated genes with functional and regulatory information, from which we deduced functional relationships for 1966 uncharacterized genes and many regulators. Using gain- and loss-of-function mutants of seven transcription factors grown under control and salt stress conditions, we experimentally validated 141 out of 271 predictions (52% precision) for 102 selected genes and mapped 148 additional transcription factor-gene regulatory interactions (49% recall). We identified an intricate core oxidative stress regulatory network where NAC13, NAC053, ERF6, WRKY6, and NAC032 transcription factors interconnect and function in detoxification. Our work shows that ensemble reverse-engineering can generate robust biological hypotheses of gene regulation in a multicellular eukaryote that can be tested by medium-throughput experimental validation.

The membrane-bound NAC transcription factor ANAC013 functions in mitochondrial retrograde regulation of the oxidative stress response in Arabidopsis.

Upon disturbance of their function by stress, mitochondria can signal to the nucleus to steer the expression of responsive genes. This mitochondria-to-nucleus communication is often referred to as mitochondrial retrograde regulation (MRR). Although reactive oxygen species and calcium are likely candidate signaling molecules for MRR, the protein signaling components in plants remain largely unknown. Through meta-analysis of transcriptome data, we detected a set of genes that are common and robust targets of MRR and used them as a bait to identify its transcriptional regulators. In the upstream regions of these mitochondrial dysfunction stimulon (MDS) genes, we found a cis-regulatory element, the mitochondrial dysfunction motif (MDM), which is necessary and sufficient for gene expression under various mitochondrial perturbation conditions. Yeast one-hybrid analysis and electrophoretic mobility shift assays revealed that the transmembrane domain-containing no apical meristem/Arabidopsis transcription activation factor/cup-shaped cotyledon transcription factors (ANAC013, ANAC016, ANAC017, ANAC053, and ANAC078) bound to the MDM cis-regulatory element. We demonstrate that ANAC013 mediates MRR-induced expression of the MDS genes by direct interaction with the MDM cis-regulatory element and triggers increased oxidative stress tolerance. In conclusion, we characterized ANAC013 as a regulator of MRR upon stress in Arabidopsis thaliana.

Reciprocal responses in the interaction between Arabidopsis and the cell-content-feeding chelicerate herbivore spider mite.

Most molecular-genetic studies of plant defense responses to arthropod herbivores have focused on insects. However, plant-feeding mites are also pests of diverse plants, and mites induce different patterns of damage to plant tissues than do well-studied insects (e.g. lepidopteran larvae or aphids). The two-spotted spider mite (Tetranychus urticae) is among the most significant mite pests in agriculture, feeding on a staggering number of plant hosts. To understand the interactions between spider mite and a plant at the molecular level, we examined reciprocal genome-wide responses of mites and its host Arabidopsis (Arabidopsis thaliana). Despite differences in feeding guilds, we found that transcriptional responses of Arabidopsis to mite herbivory resembled those observed for lepidopteran herbivores. Mutant analysis of induced plant defense pathways showed functionally that only a subset of induced programs, including jasmonic acid signaling and biosynthesis of indole glucosinolates, are central to Arabidopsis's defense to mite herbivory. On the herbivore side, indole glucosinolates dramatically increased mite mortality and development times. We identified an indole glucosinolate dose-dependent increase in the number of differentially expressed mite genes belonging to pathways associated with detoxification of xenobiotics. This demonstrates that spider mite is sensitive to Arabidopsis defenses that have also been associated with the deterrence of insect herbivores that are very distantly related to chelicerates. Our findings provide molecular insights into the nature of, and response to, herbivory for a representative of a major class of arthropod herbivores.

Evaluation of single-sample network inference methods for precision oncology.

A major challenge in precision oncology is to detect targetable cancer vulnerabilities in individual patients. Modeling high-throughput omics data in biological networks allows identifying key molecules and processes of tumorigenesis. Traditionally, network inference methods rely on many samples to contain sufficient information for learning, resulting in aggregate networks. However, to implement patient-tailored approaches in precision oncology, we need to interpret omics data at the level of individual patients. Several single-sample network inference methods have been developed that infer biological networks for an individual sample from bulk RNA-seq data. However, only a limited comparison of these methods has been made and many methods rely on 'normal tissue' samples as reference, which are not always available. Here, we conducted an evaluation of the single-sample network inference methods SSN, LIONESS, SWEET, iENA, CSN and SSPGI using transcriptomic profiles of lung and brain cancer cell lines from the CCLE database. The methods constructed functional gene networks with distinct network characteristics. Hub gene analyses revealed different degrees of subtype-specificity across methods. Single-sample networks were able to distinguish between tumor subtypes, as exemplified by node strength clustering, enrichment of known subtype-specific driver genes among hubs and differential node strength. We also showed that single-sample networks correlated better to other omics data from the same cell line as compared to aggregate networks. We conclude that single-sample network inference methods can reflect sample-specific biology when 'normal tissue' samples are absent and we point out peculiarities of each method.

Vagus nerve stimulation enhances remyelination and decreases innate neuroinflammation in lysolecithin-induced demyelination.

BACKGROUND: Current treatments for Multiple Sclerosis (MS) poorly address chronic innate neuroinflammation nor do they offer effective remyelination. The vagus nerve has a strong regulatory role in inflammation and Vagus Nerve Stimulation (VNS) has potential to affect both neuroinflammation and remyelination in MS. OBJECTIVE: This study investigated the effects of VNS on demyelination and innate neuroinflammation in a validated MS rodent model. METHODS: Lysolecithin (LPC) was injected in the corpus callosum (CC) of 46 Lewis rats, inducing a demyelinated lesion. 33/46 rats received continuously-cycled VNS (cVNS) or one-minute per day VNS (1minVNS) or sham VNS from 2 days before LPC-injection until perfusion at 3 days post-injection (dpi) (corresponding with a demyelinated lesion with peak inflammation). 13/46 rats received cVNS or sham from 2 days before LPC-injection until perfusion at 11 dpi (corresponding with a partial remyelinated lesion). Immunohistochemistry and proteomics analyses were performed to investigate the extend of demyelination and inflammation. RESULTS: Immunohistochemistry showed that cVNS significantly reduced microglial and astrocytic activation in the lesion and lesion border, and significantly reduced the Olig2+ cell count at 3 dpi. Furthermore, cVNS significantly improved remyelination with 57.4 % versus sham at 11 dpi. Proteomic gene set enrichment analyses showed increased activation of (glutamatergic) synapse pathways in cVNS versus sham, most pronounced at 3 dpi. CONCLUSION: cVNS improved remyelination of an LPC-induced lesion. Possible mechanisms might include modulation of microglia and astrocyte activity, increased (glutamatergic) synapses and enhanced oligodendrocyte clearance after initial injury.

HTSplotter: An end-to-end data processing, analysis and visualisation tool for chemical and genetic in vitro perturbation screening.

In biomedical research, high-throughput screening is often applied as it comes with automatization, higher-efficiency, and more and faster results. High-throughput screening experiments encompass drug, drug combination, genetic perturbagen or a combination of genetic and chemical perturbagen screens. These experiments are conducted in real-time assays over time or in an endpoint assay. The data analysis consists of data cleaning and structuring, as well as further data processing and visualisation, which, due to the amount of data, can easily become laborious, time-consuming and error-prone. Therefore, several tools have been developed to aid researchers in this process, but these typically focus on specific experimental set-ups and are unable to process data of several time points and genetic-chemical perturbagen screens. To meet these needs, we developed HTSplotter, a web tool and Python module that performs automatic data analysis and visualization of visualization of eitherendpoint or real-time assays from different high-throughput screening experiments: drug, drug combination, genetic perturbagen and genetic-chemical perturbagen screens. HTSplotter implements an algorithm based on conditional statements to identify experiment types and controls. After appropriate data normalization, including growth rate normalization, HTSplotter executes downstream analyses such as dose-response relationship and drug synergism assessment by the Bliss independence (BI), Zero Interaction Potency (ZIP) and Highest Single Agent (HSA) methods. All results are exported as a text file and plots are saved in a PDF file. The main advantage of HTSplotter over other available tools is the automatic analysis of genetic-chemical perturbagen screens and real-time assays where growth rate and perturbagen effect results are plotted over time. In conclusion, HTSplotter allows for the automatic end-to-end data processing, analysis and visualisation of various high-throughput in vitro cell culture screens, offering major improvements in terms of versatility, efficiency and time over existing tools.

Human iPSC modeling recapitulates in vivo sympathoadrenal development and reveals an aberrant developmental subpopulation in familial neuroblastoma.

Studies defining normal and disrupted human neural crest cell development have been challenging given its early timing and intricacy of development. Consequently, insight into the early disruptive events causing a neural crest related disease such as pediatric cancer neuroblastoma is limited. To overcome this problem, we developed an in vitro differentiation model to recapitulate the normal in vivo developmental process of the sympathoadrenal lineage which gives rise to neuroblastoma. We used human in vitro pluripotent stem cells and single-cell RNA sequencing to recapitulate the molecular events during sympathoadrenal development. We provide a detailed map of dynamically regulated transcriptomes during sympathoblast formation and illustrate the power of this model to study early events of the development of human neuroblastoma, identifying a distinct subpopulation of cell marked by SOX2 expression in developing sympathoblast obtained from patient derived iPSC cells harboring a germline activating mutation in the anaplastic lymphoma kinase (ALK) gene.

Opportunities and challenges in interpretable deep learning for drug sensitivity prediction of cancer cells.

In precision oncology, therapy stratification is done based on the patients' tumor molecular profile. Modeling and prediction of the drug response for a given tumor molecular type will further improve therapeutic decision-making for cancer patients. Indeed, deep learning methods hold great potential for drug sensitivity prediction, but a major problem is that these models are black box algorithms and do not clarify the mechanisms of action. This puts a limitation on their clinical implementation. To address this concern, many recent studies attempt to overcome these issues by developing interpretable deep learning methods that facilitate the understanding of the logic behind the drug response prediction. In this review, we discuss strengths and limitations of recent approaches, and suggest future directions that could guide further improvement of interpretable deep learning in drug sensitivity prediction in cancer research.

Whole transcriptome profiling of liquid biopsies from tumour xenografted mouse models enables specific monitoring of tumour-derived extracellular RNA.

While cell-free DNA (cfDNA) is widely being investigated, free circulating RNA (extracellular RNA, exRNA) has the potential to improve cancer therapy response monitoring and detection due to its dynamic nature. However, it remains unclear in which blood subcompartment tumour-derived exRNAs primarily reside. We developed a host-xenograft deconvolution framework, exRNAxeno, with mapping strategies to either a combined human-mouse reference genome or both species genomes in parallel, applicable to exRNA sequencing data from liquid biopsies of human xenograft mouse models. The tool enables to distinguish (human) tumoural RNA from (murine) host RNA, to specifically analyse tumour-derived exRNA. We applied the combined pipeline to total exRNA sequencing data from 95 blood-derived liquid biopsy samples from 30 mice, xenografted with 11 different tumours. Tumoural exRNA concentrations are not determined by plasma platelet levels, while host exRNA concentrations increase with platelet content. Furthermore, a large variability in exRNA abundance and transcript content across individual mice is observed. The tumoural gene detectability in plasma is largely correlated with the RNA expression levels in the tumour tissue or cell line. These findings unravel new aspects of tumour-derived exRNA biology in xenograft models and open new avenues to further investigate the role of exRNA in cancer.

RRM2 enhances MYCN-driven neuroblastoma formation and acts as a synergistic target with CHK1 inhibition.

High-risk neuroblastoma, a pediatric tumor originating from the sympathetic nervous system, has a low mutation load but highly recurrent somatic DNA copy number variants. Previously, segmental gains and/or amplifications allowed identification of drivers for neuroblastoma development. Using this approach, combined with gene dosage impact on expression and survival, we identified ribonucleotide reductase subunit M2 (RRM2) as a candidate dependency factor further supported by growth inhibition upon in vitro knockdown and accelerated tumor formation in a neuroblastoma zebrafish model coexpressing human RRM2 with MYCN. Forced RRM2 induction alleviates excessive replicative stress induced by CHK1 inhibition, while high RRM2 expression in human neuroblastomas correlates with high CHK1 activity. MYCN-driven zebrafish tumors with RRM2 co-overexpression exhibit differentially expressed DNA repair genes in keeping with enhanced ATR-CHK1 signaling activity. In vitro, RRM2 inhibition enhances intrinsic replication stress checkpoint addiction. Last, combinatorial RRM2-CHK1 inhibition acts synergistic in high-risk neuroblastoma cell lines and patient-derived xenograft models, illustrating the therapeutic potential.

From DNA Copy Number Gains and Tumor Dependencies to Novel Therapeutic Targets for High-Risk Neuroblastoma.

Neuroblastoma is a pediatric tumor arising from the sympatho-adrenal lineage and a worldwide leading cause of childhood cancer-related deaths. About half of high-risk patients die from the disease while survivors suffer from multiple therapy-related side-effects. While neuroblastomas present with a low mutational burden, focal and large segmental DNA copy number aberrations are highly recurrent and associated with poor survival. It can be assumed that the affected chromosomal regions contain critical genes implicated in neuroblastoma biology and behavior. More specifically, evidence has emerged that several of these genes are implicated in tumor dependencies thus potentially providing novel therapeutic entry points. In this review, we briefly review the current status of recurrent DNA copy number aberrations in neuroblastoma and provide an overview of the genes affected by these genomic variants for which a direct role in neuroblastoma has been established. Several of these genes are implicated in networks that positively regulate MYCN expression or stability as well as cell cycle control and apoptosis. Finally, we summarize alternative approaches to identify and prioritize candidate copy-number driven dependency genes for neuroblastoma offering novel therapeutic opportunities.

A quantitative in silico analysis calculates the angiotensin I converting enzyme (ACE) inhibitory activity in pea and whey protein digests.

Angiotensin I converting enzyme (ACE) inhibitory peptides can induce antihypertensive effects after oral administration. By means of an ACE inhibitory peptide database, containing about 500 reported sequences and their IC(50) values, the different proteins in pea and whey were quantitatively evaluated as precursors for ACE inhibitory peptides. This analysis was combined with experimental data from the evolution in ACE inhibitory activity and protein degradation during in vitro gastrointestinal digestion. Pea proteins produced similar in silico scores and were degraded early in the in vitro digestion. High ACE inhibitory activity was observed after the simulated stomach phase and augmented slightly in the simulated small intestine phase. The major whey protein beta-lactoglobulin obtained the highest in silico scores, which corresponded with the fact that degradation of this protein in vitro only occurred from the simulated small intestine phase on and resulted in a 10-fold increase in the ACE inhibitory activity. Whey protein obtained total in silico scores of about 124 ml/mg, compared to 46 ml/mg for pea protein, indicating that whey protein would be a richer source of ACE inhibitory peptides than pea protein. Although beta-lactoglobulin is only partially digested, a higher ACE inhibitory activity was indeed found in the whey (IC(50) = 0.048 mg/ml) compared to the pea digest (IC(50) = 0.076 mg/ml). In silico gastrointestinal digestion of the highest scoring proteins in pea and whey, vicilin and albumin PA2, and beta-lactoglobulin, respectively, directly released a number of potent ACE inhibitory peptides. Several other ACE inhibitory sequences resisted in silico digestion by gastrointestinal proteases. Briefly, the quantitative in silico analysis will facilitate the study of precursor proteins on a large scale and the specific release of bioactive peptides.

Influence of the lactokinin Ala-Leu-Pro-Met-His-Ile-Arg (ALPMHIR) on the release of endothelin-1 by endothelial cells.

Milk protein-derived peptides with angiotensin-converting enzyme (ACE) inhibitory activity can reduce blood pressure in hypertensive subjects. The lactokinin Ala-Leu-Pro-Met-His-Ile-Arg (ALPMHIR) is an ACE-inhibitory peptide released by tryptic digestion from the milk protein beta-lactoglobulin. Its ACE-inhibitory activity is 100 times lower than that of captopril. The latter is known to inhibit the release of the vasoconstrictor endothelin-1 (ET-1) by endothelial cells. The effects of ALPMHIR on the endothelium are currently unknown. In this study, the influence of ALPMHIR on release of ET-1 by endothelial cells was investigated. The basal ET-1 release of the cells was reduced by 29% (p<0.01) in the presence of 1 mM ALPMHIR, compared to 42% (p<0.01) for 0.1 mM captopril. Addition of 10 U/ml thrombin to the incubation medium increased the release of ET-1 by 66% (p<0.01). Co-incubation of 10 U/ml thrombin with 1 microM captopril or with 0.1 mM ALPMHIR inhibited the stimulated ET-1 release by 45% (p<0.01) and by 32% (p<0.01), respectively. These data indicate that dietary peptides, such as ALPMHIR, can modulate ET-1 release by endothelial cells. These effects, among other mechanisms, may play a role in the anti-hypertensive effect of milk protein-derived peptides.

Optimisation and validation of an angiotensin-converting enzyme inhibition assay for the screening of bioactive peptides.

Angiotensin-converting enzyme (ACE) plays a major role in the regulation of blood pressure. A diagnostic assay to measure angiotensin-converting enzyme (ACE) activity was transformed into an enzyme inhibition assay and optimised, which led to a more sensitive and less expensive assay. By this spectrophotometric method, ACE inhibition is measured using the substrate furanacryloyl-Phe-Gly-Gly and as ACE source rabbit lung acetone extract. The optimised as well as the original ACE inhibition assay were used to verify the ACE inhibitory activity of captopril. The ACE inhibition assay was further validated by enalapril, its active derivative enalaprilat and the ACE-inhibitory peptide Ala-Leu-Pro-Met-His-Ile-Arg, corresponding to a tryptic fragment of bovine beta-lactoglobulin. Sigmoid curves could be fit adequately to the data points representing ACE inhibition in function of inhibitor concentration. IC(50) values for these compounds corresponded well with literature data. Furthermore, pea and whey protein hydrolysates obtained by digestion with trypsin showed ACE inhibitory activity in the ACE inhibition assay. Hence, this optimised assay is suitable to screen for ACE inhibitory peptides derived from food proteins with a possible antihypertensive effect in vivo.

Release of angiotensin I converting enzyme (ACE) inhibitory activity during in vitro gastrointestinal digestion: from batch experiment to semicontinuous model.

Gastrointestinal digestion is of major importance in the bioavailability of angiotensin I converting enzyme (ACE) inhibitory peptides, bioactive peptides with possible antihypertensive effects. In this study, the conditions of in vitro gastrointestinal digestion leading to the formation and degradation of ACE inhibitory peptides were investigated for pea and whey protein. In batch experiments, the digestion simulating the physiological conditions sufficed to achieve the highest ACE inhibitory activity, with IC(50) values of 0.076 mg/mL for pea and 0.048 mg/mL for whey protein. The degree of proteolysis did not correlate with the ACE inhibitory activity and was always higher for pea than whey. In a semicontinuous model of gastrointestinal digestion, response surface methodology studied the influence of temperature and incubation time in both the stomach and small intestine phases on the ACE inhibitory activity and degree of proteolysis. For pea protein, a linear model for the degree of proteolysis and a quadratic model for the ACE inhibitory activity could be constituted. Within the model, a maximal degree of proteolysis was observed at the highest temperature and the longest incubation time in the small intestine phase, while maximal ACE inhibitory activity was obtained at the longest incubation times in the stomach and small intestine phase. These results show that ACE inhibitory activity of pea and whey hydrolysates can be controlled by the conditions of in vitro gastrointestinal digestion.

Identification of cis-regulatory elements specific for different types of reactive oxygen species in Arabidopsis thaliana.

The type of reactive oxygen species (ROS) is a major factor that determines the specificity of biological responses. These responses may be elicited by activation of transcription factors that recognize ROS-specific cis-regulatory elements in target genes. In search for Arabidopsis promoter motifs specific for particular types of ROS, genome-wide microarray expression profiles for 283 abiotic stress-related conditions were subjected to cluster analysis to identify gene groups induced by singlet oxygen, superoxide radicals, and H(2)O(2). Promoters of these gene groups were analyzed to identify cis-regulatory elements that are associated with specific types of ROS. Eleven ROS-specific de novo identified elements, seven known promoter motifs and several sequences enriched in ROS-responsive clusters but lacking in specificity are reported. The conservation of the identified motifs was determined in orthologous genes in C. papaya, V. vinifera and P. trichocarpa. Finally, biological functions were attributed to the motifs by calculation of GO-term enrichment for genes with conserved ROS-responsive elements.

Transcription regulatory networks in Caenorhabditis elegans inferred through reverse-engineering of gene expression profiles constitute biological hypotheses for metazoan development.

Differential gene expression governs the development, function and pathology of multicellular organisms. Transcription regulatory networks study differential gene expression at a systems level by mapping the interactions between regulatory proteins and target genes. While microarray transcription profiles are the most abundant data for gene expression, it remains challenging to correctly infer the underlying transcription regulatory networks. The reverse-engineering algorithm LeMoNe (learning module networks) uses gene expression profiles to extract ensemble transcription regulatory networks of coexpression modules and their prioritized regulators. Here we apply LeMoNe to a compendium of microarray studies of the worm Caenorhabditis elegans. We obtain 248 modules with a regulation program for 5020 genes and 426 regulators and a total of 24 012 predicted transcription regulatory interactions. Through GO enrichment analysis, comparison with the gene-gene association network WormNet and integration of other biological data, we show that LeMoNe identifies functionally coherent coexpression modules and prioritizes regulators that relate to similar biological processes as the module genes. Furthermore, we can predict new functional relationships for uncharacterized genes and regulators. Based on modules involved in molting, meiosis and oogenesis, ciliated sensory neurons and mitochondrial metabolism, we illustrate the value of LeMoNe as a biological hypothesis generator for differential gene expression in greater detail. In conclusion, through reverse-engineering of C. elegans expression data, we obtained transcription regulatory networks that can provide further insight into metazoan development.