Modelling the regulatory mechanisms that determine cell fate, response to external perturbation, and disease state depends on measuring many factors, a task made more difficult by the plasticity of the epigenome. Scanning the genome for the sequence patterns defined by Position Weight Matrices (PWM) can be used to estimate transcription factor (TF) binding locations. However, this approach does not incorporate information regarding the epigenetic context necessary for TF binding. CpG methylation is an epigenetic mark influenced by environmental factors that is commonly assayed in human cohort studies. We developed a framework to score inferred TF binding locations using methylation data. We intersected motif locations identified using PWMs with methylation information captured in both whole-genome bisulfite sequencing and Illumina EPIC array data for six cell lines, scored motif locations based on these data, and compared with experimental data characterizing TF binding (ChIP-seq). We found that for most TFs, binding prediction improves using methylation-based scoring compared to standard PWM-scores. We also illustrate that our approach can be generalized to infer TF binding when methylation information is only proximally available, i.e. measured for nearby CpGs that do not directly overlap with a motif location. Overall, our approach provides a framework for inferring context-specific TF binding using methylation data. Importantly, the availability of DNA methylation data in existing patient populations provides an opportunity to use our approach to understand the impact of methylation on gene regulatory processes in the context of human disease.

Single-cell data is increasingly used for gene regulatory network (GRN) inference, and benchmarks for this have been developed based on simulated data. However, existing single-cell simulators cannot model the effects of gene perturbations. A further challenge lies in generating large-scale GRNs that often struggle with computational and stability issues. We present GeneSPIDER2, an update of the GeneSPIDER MATLAB toolbox for GRN benchmarking, inference, and analysis. Several software modules have improved capabilities and performance, and new functionalities have been added. A major improvement is the ability to generate large GRNs with biologically realistic topological properties in terms of scale-free degree distribution and modularity. Another major addition is a simulation of single-cell data, which is becoming increasingly popular as input for GRN inference. Specifically, we introduced the unique feature to generate single-cell data based on genetic perturbations. Finally, the simulated single-cell data was compared to real single-cell Perturb-seq data from two cell lines, showing that the synthetic and real data exhibit similar properties.

Neflamapimod, a selective inhibitor of the alpha isoform of p38 mitogen-activated protein kinase (MAPKα), was investigated for its potential to inhibit lipopolysaccharide (LPS)-induced activation of endothelial cells (ECs), adhesion molecule induction, and subsequent leukocyte attachment to EC monolayers. These events are known to contribute to vascular inflammation and cardiovascular dysfunction. Our results demonstrate that LPS treatment of cultured ECs and rats leads to significant upregulation of adhesion molecules, both in vitro and in vivo, which can be effectively inhibited by neflamapimod treatment. Western blotting data further reveals that neflamapimod inhibits LPS-induced phosphorylation of p38 MAPKα and the activation of NF-κB signaling in ECs. Additionally, leukocyte adhesion assays demonstrate a substantial reduction in leukocyte attachment to cultured ECs and the aorta lumen of rats treated with neflamapimod. Consistent with vascular inflammation, LPS-treated rat arteries exhibit significantly diminished vasodilation response to acetylcholine, however, arteries from rats treated with neflamapimod maintain their vasodilation capacity, demonstrating its ability to limit LPS-induced vascular inflammation. Overall, our data demonstrate that neflamapimod effectively inhibits endothelium activation, adhesion molecule expression, and leukocyte attachment, thereby reducing vascular inflammation.

Inference and analysis of gene regulatory networks (GRNs) require software that integrates multi-omic data from various sources. The Network Zoo (netZoo; netzoo.github.io) is a collection of open-source methods to infer GRNs, conduct differential network analyses, estimate community structure, and explore the transitions between biological states. The netZoo builds on our ongoing development of network methods, harmonizing the implementations in various computing languages and between methods to allow better integration of these tools into analytical pipelines. We demonstrate the utility using multi-omic data from the Cancer Cell Line Encyclopedia. We will continue to expand the netZoo to incorporate additional methods.

Neflamapimod, a selective inhibitor of p38 mitogen activated protein kinase alpha (MAPKα), is under clinical investigation for its efficacy in Alzheimer's disease (AD) and dementia with Lewy Bodies (DLB). Here, we investigated if neflamapimod-mediated acute inhibition of p38 MAPKα could induce vasodilation in resistance-size rat mesenteric arteries. Our pressure myography data demonstrated that neflamapimod produced a dose-dependent vasodilation in mesenteric arteries. Our Western blotting data revealed that acute neflamapimod treatment significantly reduced the phosphorylation of p38 MAPKα and its downstream target heat-shock protein 27 (Hsp27) involved in cytoskeletal reorganization and smooth muscle contraction. Likewise, non-selective inhibition of p38 MAPK by SB203580 attenuated p38 MAPKα and Hsp27 phosphorylation, and induced vasodilation. Endothelium denudation or pharmacological inhibition of endothelium-derived vasodilators such as nitric oxide (NO) and prostacyclin (PGI2) had no effect on such vasodilation. Neflamapimod-evoked vasorelaxation remained unaltered by the inhibition of smooth muscle cell K+ channels. Altogether, our data for the first time demonstrates that in resistance mesenteric arteries, neflamapimod inhibits p38 MAPKα and phosphorylation of its downstream actin-associated protein Hsp27, leading to vasodilation. This novel finding may be clinically significant and is likely to improve systemic blood pressure and cognitive deficits in AD and DLB patients for which neflamapimod is being investigated.

Neflamapimod, a selective inhibitor of p38 MAPKα, is under clinical investigation for its efficacy in Alzheimer’s disease (AD) and dementia with Lewy Bodies (DLB). Here, we investigated if neflamapimod-mediated acute inhibition of p38 MAPKα could induce vasodilation in resistance-size rat mesenteric arteries. Our pressure myography data demonstrated that neflamapimod produced a dose-dependent vasodilation in mesenteric arteries. Our Western blotting data revealed that acute neflamapimod treatment significantly reduced the phosphorylation of p38 MAPKα and its downstream target heat-shock protein 27 (Hsp27) that is involved in cytoskeletal reorganization and smooth muscle contraction. Likewise, non-selective inhibition of p38 MAPK by SB203580 attenuated p38 MAPKα and Hsp27 phosphorylation, and induced vasodilation. Endothelium denudation or pharmacological inhibition of endothelium-derived vasodilators such as nitric oxide (NO) and prostacyclin (PGI2) had no effect on such vasodilation. Neflamapimod-evoked vasorelaxation remained unaltered by the inhibition of smooth muscle cell K+ channels. Altogether, our data for the first time demonstrates that in resistance mesenteric arteries, neflamapimod inhibits p38 MAPKα and phosphorylation of its downstream actin-associated protein Hsp27, leading to vasodilation. This novel finding may be clinically significant and is likely to improve systemic blood pressure and cognitive deficits in AD and DLB patients for which neflamapimod is being investigated.

Human exposure to pentabromodiphenyl ether (PBDE) mixture (DE-71) and its PBDE-47 congener can occur both in utero and during lactation. Here, we tested the hypothesis that PBDE-induced neonatal hepatic transcriptomic alterations in Wistar Han rat pups can inform on potential toxicity and carcinogenicity after longer term PBDE exposures. Wistar Han rat dams were exposed to either DE-71 or PBDE-47 daily from gestation day (GD 6) through postnatal day 4 (PND 4). Total plasma thyroxine (T4) was decreased in PND 4 pups. In liver, transcripts for CYPs and conjugation enzymes, Nrf2, and ABC transporters were upregulated. In general, the hepatic transcriptomic alterations after exposure to DE-71 or PBDE-47 were similar and provided early indicators of oxidative stress and metabolic alterations, key characteristics of toxicity processes. The transcriptional benchmark dose lower confidence limits of the most sensitive biological processes were lower for PBDE-47 than for the PBDE mixture. Neonatal rat liver transcriptomic data provide early indicators on molecular pathway alterations that may lead to toxicity and/or carcinogenicity if the exposures continue for longer durations. These early toxicogenomic indicators may be used to help prioritize chemicals for a more complete toxicity and cancer risk evaluation.

The interactions among the components of a living cell that constitute the gene regulatory network (GRN) can be inferred from perturbation-based gene expression data. Such networks are useful for providing mechanistic insights of a biological system. In order to explore the feasibility and quality of GRN inference at a large scale, we used the L1000 data where ~1000 genes have been perturbed and their expression levels have been quantified in 9 cancer cell lines. We found that these datasets have a very low signal-to-noise ratio (SNR) level causing them to be too uninformative to infer accurate GRNs. We developed a gene reduction pipeline in which we eliminate uninformative genes from the system using a selection criterion based on SNR, until reaching an informative subset. The results show that our pipeline can identify an informative subset in an overall uninformative dataset, allowing inference of accurate subset GRNs. The accurate GRNs were functionally characterized and potential novel cancer-related regulatory interactions were identified.

The gene regulatory network (GRN) of human cells encodes mechanisms to ensure proper functioning. However, if this GRN is dysregulated, the cell may enter into a disease state such as cancer. Understanding the GRN as a system can therefore help identify novel mechanisms underlying disease, which can lead to new therapies. To deduce regulatory interactions relevant to cancer, we applied a recent computational inference framework to data from perturbation experiments in squamous carcinoma cell line A431. GRNs were inferred using several methods, and the false discovery rate was controlled by the NestBoot framework. We developed a novel approach to assess the predictiveness of inferred GRNs against validation data, despite the lack of a gold standard. The best GRN was significantly more predictive than the null model, both in cross-validated benchmarks and for an independent dataset of the same genes under a different perturbation design. The inferred GRN captures many known regulatory interactions central to cancer-relevant processes in addition to predicting many novel interactions, some of which were experimentally validated, thus providing mechanistic insights that are useful for future cancer research.

Motivation: The interactions among the components of a living cell that constitute the gene regulatory network (GRN) can be inferred from perturbation-based gene expression data. Such networks are ...

2,3-Hexanedione is a 6-carbon α-diketone flavoring agent that is a potential replacement for 2,3-butanedione. 2,3-Butanedione inhalation has been associated with obliterative bronchiolitis in workers. The inhalation toxicity data for 2,3-hexanedione are limited to one study conducted in rats. In this study, we provide additional data on the respiratory toxicity of 2,3-hexanedione in male and female B6C3F1 mice. Mice (six per group) were exposed (whole body) to 0, 100, 150, or 200 ppm 2,3-hexanedione 6 hours/day, 5 days/week for 2 weeks plus 2 days (12 exposure days total). Animals were euthanized the morning after the last exposure and respiratory tract tissues collected for histopathological evaluation. All mice survived to the end of the study. Body weights of female mice exposed to 150 and 200 ppm 2,3-hexanedione were significantly less than controls, and relative lung weights were significantly greater following exposure to 100, 150, and 200 ppm. Squamous metaplasia was the most prevalent lesion observed in the nasal cavity, larynx, and trachea of exposed male and female mice. 2,3-Hexanedione caused atypical hyperplasia of the bronchial respiratory epithelium in male and female mice exposed to 200 ppm 2,3-hexanedione. 2,3-Hexanedione exposure did not cause airway fibrosis or obliterative bronchiolitis in mice.

2,3-Hexanedione is a 6-carbon -diketone flavoring agent that is a potential replacement for 2,3-butanedione.2,3-Butanedione inhalation has been associated with obliterative bronchiolitis in workers.The inhalation toxicity data for 2,3-hexanedione are limited to one study conducted in rats.In this study, we provide additional data on the respiratory toxicity of 2,3-hexanedione in male and female B6C3F1 mice.Mice (six per group) were exposed (whole body) to 0, 100, 150, or 200 ppm 2,3-hexanedione 6 hours/day, 5 days/week for 2 weeks plus 2 days (12 exposure days total).Animals were euthanized the morning after the last exposure and respiratory tract tissues collected for histopathological evaluation.All mice survived to the end of the study.Body weights of female mice exposed to 150 and 200 ppm 2,3-hexanedione were significantly less than controls, and relative lung weights were significantly greater following exposure to 100, 150, and 200 ppm.Squamous metaplasia was the most prevalent lesion observed in the nasal cavity, larynx, and trachea of exposed male and female mice.2,3-Hexanedione caused atypical hyperplasia of the bronchial respiratory epithelium in male and female mice exposed to 200 ppm 2,3-hexanedione.2,3-Hexanedione exposure did not cause airway fibrosis or obliterative bronchiolitis in mice.

MOTIVATION:Inference of gene regulatory networks (GRNs) from perturbation data can give detailed mechanistic insights of a biological system. Many inference methods exist, but the resulting GRN is generally sensitive to the choice of method-specific parameters. Even though the inferred GRN is optimal given the parameters, many links may be wrong or missing if the data is not informative. To make GRN inference reliable, a method is needed to estimate the support of each predicted link as the method parameters are varied.RESULTS:To achieve this we have developed a method called nested bootstrapping, which applies a bootstrapping protocol to GRN inference, and by repeated bootstrap runs assesses the stability of the estimated support values. To translate bootstrap support values to false discovery rates we run the same pipeline with shuffled data as input. This provides a general method to control the false discovery rate of GRN inference that can be applied to any setting of inference parameters, noise level, or data properties. We evaluated nested bootstrapping on a simulated dataset spanning a range of such properties, using the LASSO, Least Squares, RNI, GENIE3 and CLR inference methods. An improved inference accuracy was observed in almost all situations. Nested bootstrapping was incorporated into the GeneSPIDER package, which was also used for generating the simulated networks and data, as well as running and analyzing the inferences.AVAILABILITY AND IMPLEMENTATION:https://bitbucket.org/sonnhammergrni/genespider/src/NB/%2B Methods/NestBoot.m.

Motivation: The interactions among the components of a living cell that constitute the gene regulatory network (GRN) can be inferred from perturbation-based gene expression data. Such networks are ...

Nrf2 is essential to antioxidant response element (ARE)-mediated host defense. Sulforaphane (SFN) is a phytochemical antioxidant known to affect multiple cellular targets including Nrf2-ARE pathway in chemoprevention. However, the role of SFN in non-malignant airway disorders remain unclear. To test if pre-activation of Nrf2-ARE signaling protects lungs from oxidant-induced acute injury, wild-type (Nrf2+/+) and Nrf2-deficient (Nrf2−/−) mice were given SFN orally or as standardized broccoli sprout extract diet (SBE) before hyperoxia or air exposure. Hyperoxia-induced pulmonary injury and oxidation indices were significantly reduced by SFN or SBE in Nrf2+/+ mice but not in Nrf2−/− mice. SFN upregulated a large cluster of basal lung genes that are involved in mitochondrial oxidative phosphorylation, energy metabolism, and cardiovascular protection only in Nrf2+/+ mice. Bioinformatic analysis elucidated ARE-like motifs on these genes. Transcript abundance of the mitochondrial machinery genes remained significantly higher after hyperoxia exposure in SFN-treated Nrf2+/+ mice than in SFN-treated Nrf2−/− mice. Nuclear factor-κB was suggested to be a central molecule in transcriptome networks affected by SFN. Minor improvement of hyperoxia-caused lung histopathology and neutrophilia by SFN in Nrf2−/− mice implies Nrf2-independent or alternate effector mechanisms. In conclusion, SFN is suggested to be as a preventive intervention in a preclinical model of acute lung injury by linking mitochondria and Nrf2. Administration of SFN alleviated acute lung injury-like pathogenesis in a Nrf2-dependent manner. Potential AREs in the SFN-inducible transcriptome for mitochondria bioenergetics provided a new insight into the downstream mechanisms of Nrf2-mediated pulmonary protection.

DNA methylation plays a key role in X-chromosome inactivation (XCI), a process that achieves dosage compensation for X-encoded gene products between mammalian female and male cells. However, differential sex chromosome dosage complicates genome-wide epigenomic assessments, and the X chromosome is frequently excluded from female-to-male comparative analyses. Using the X chromosome in the sexually dimorphic mouse liver as a model, we provide a general framework for comparing base-resolution DNA methylation patterns across samples that have different chromosome numbers and ask at a systematic level if predictions by historical analyses of X-linked DNA methylation hold true at a base-resolution chromosome-wide level. We demonstrate that sex-specific methylation patterns on the X chromosome largely reflect the effects of XCI. While our observations concur with longstanding observations of XCI at promoter-proximal CpG islands, we provide evidence that sex-specific DNA methylation differences are not limited to CpG island boundaries. Moreover, these data support a model in which maintenance of CpG islands in the inactive state does not require complete regional methylation. Further, we validate an intragenic non-CpG methylation signature in genes escaping XCI in mouse liver. Our analyses provide insight into underlying methylation patterns that should be considered when assessing sex differences in genome-wide methylation analyses.

Pentabromodiphenyl ethers (PBDE) are found in human tissue, in household dust, and in the environment, and a particular concern is the potential for the induction of cancer pathways from these fat-soluble persistent organic pollutants. Only one PBDE cancer study has been conducted and that was for a PBDE mixture (DE-71). Because it is not feasible to test all PBDE congeners in the environment for cancer potential, it is important to develop a set of biological endpoints that can be used in short-term toxicity studies to predict disease outcome after long-term exposures. In this study, PBDE-47 was selected as the test PBDE congener to evaluate and compare toxicity to that of the carcinogenic PBDE mixture. The toxicities of PBDE-47 and the PBDE mixture were evaluated at PND 22 in Wistar Han rat (Crl: WI (Han)) pups after in utero/postnatal exposure (0, 0.1, 15, or 50 mg/kg; dams, GD6-21; pups, PND 12-PND 21; oral gavage daily dosing). By PND 22, PBDE-47 caused centrilobular hypertrophy and fatty change in liver, and reduced serum thyroxin (T4) levels; similar effects were also observed after PBDE mixture exposure. Transcriptomic changes in the liver included induction of cytochrome p450 transcripts and up-regulation of Nrf2 antioxidant pathway transcripts and ABC membrane transport transcripts. Decreases in other transport transcripts (ABCG5 & 8) provided a plausible mechanism for lipid accumulation, characterized by a treatment-related liver fatty change after PBDE-47 and PBDE mixture exposure. The benchmark dose calculation based on liver transcriptomic data was generally lower for PBDE-47 than for the PBDE mixture. The up-regulation of the Nrf2 antioxidant pathway and changes in metabolic transcripts after PBDE-47 and PBDE mixture exposure suggest that PBDE-47, like the PBDE mixture (NTP 2016, TR 589), could be a liver toxin/carcinogen after long-term exposure.