Quantification of biological network perturbations for mechanistic insight and diagnostics using two-layer causal models

Tobacco Alkaloid Assessment in a DSS-Induced Colitis Mouse Model with a Fully Humanized Immune System

Inflammatory bowel disease (IBD) refers to chronic intestinal immune-mediated diseases including two main disease manifestations: ulcerative colitis (UC) and Crohn’s disease (CD). Epidemiological, clinical, and preclinical evidence has highlighted the potential anti-inflammatory properties of naturally occurring alkaloids. In the present study, we investigated the potential anti-inflammatory activities of the tobacco alkaloids nicotine and anatabine in a dextran sulfate sodium (DSS)-induced UC mouse model with a fully humanized immune system. Our results show that nicotine significantly reduced all acute colitis symptoms and improved colitis-specific endpoints, including histopathologically assessed colon inflammation, tissue damage, and mononuclear cell infiltration. The tobacco alkaloid anatabine showed similar effectiveness trends, although they were generally weaker or not significant. Gene expression analysis in the context of biological network models of IBD further pinpointed a possible mechanism by which nicotine attenuated DSS-induced colitis in humanized mice. The current study enables further investigation of possible molecular mechanisms by which tobacco alkaloids attenuate UC symptoms.

Causal biological network models for reactive astrogliosis: a systems approach to neuroinflammation

Astrocytes play a central role in the neuroimmune response by responding to CNS pathologies with diverse molecular and morphological changes during the process of reactive astrogliosis. Here, we used a computational biological network model and mathematical algorithms that allow the interpretation of high-throughput transcriptomic datasets in the context of known biology to study reactive astrogliosis. We gathered available mechanistic information from the literature into a comprehensive causal biological network (CBN) model of astrocyte reactivity. The CBN model was built in the Biological Expression Language, which is both human-readable and computable. We characterized the CBN with a network analysis of highly connected nodes and demonstrated that the CBN captures relevant astrocyte biology. Subsequently, we used the CBN and transcriptomic data to identify key molecular pathways driving the astrocyte phenotype in four CNS pathologies: samples from mouse models of lipopolysaccharide-induced endotoxemia, Alzheimer’s disease, and amyotrophic lateral sclerosis; and samples from multiple sclerosis patients. The astrocyte CBN provides a new tool to identify causal mechanisms and quantify astrogliosis based on transcriptomic data.

Causal reasoning over knowledge graphs leveraging drug-perturbed and disease-specific transcriptomic signatures for drug discovery

Network-based approaches are becoming increasingly popular for drug discovery as they provide a systems-level overview of the mechanisms underlying disease pathophysiology. They have demonstrated significant early promise over other methods of biological data representation, such as in target discovery, side effect prediction and drug repurposing. In parallel, an explosion of -omics data for the deep characterization of biological systems routinely uncovers molecular signatures of disease for similar applications. Here, we present RPath, a novel algorithm that prioritizes drugs for a given disease by reasoning over causal paths in a knowledge graph (KG), guided by both drug-perturbed as well as disease-specific transcriptomic signatures. First, our approach identifies the causal paths that connect a drug to a particular disease. Next, it reasons over these paths to identify those that correlate with the transcriptional signatures observed in a drug-perturbation experiment, and anti-correlate to signatures observed in the disease of interest. The paths which match this signature profile are then proposed to represent the mechanism of action of the drug. We demonstrate how RPath consistently prioritizes clinically investigated drug-disease pairs on multiple datasets and KGs, achieving better performance over other similar methodologies. Furthermore, we present two case studies showing how one can deconvolute the predictions made by RPath as well as predict novel targets.

Causal Biological Network Model for Inflammasome Signaling Applied for Interpreting Transcriptomic Changes in Various Inflammatory States

Virtually any stressor that alters the cellular homeostatic state may result in an inflammatory response. As a critical component of innate immunity, inflammasomes play a prominent role in the inflammatory response. The information on inflammasome biology is rapidly growing, thus creating the need for structuring it into a model that can help visualize and enhance the understanding of underlying biological processes. Causal biological network (CBN) models provide predictive power for novel disease mechanisms and treatment outcomes. We assembled the available literature information on inflammasome activation into the CBN model and scored it with publicly available transcriptomic datasets that address viral infection of the lungs, osteo- and rheumatoid arthritis, psoriasis, and aging. The scoring inferred pathway activation leading to NLRP3 inflammasome activation in these diverse conditions, demonstrating that the CBN model provides a platform for interpreting transcriptomic data in the context of inflammasome activation.

A Functional Network Model of the Metastasis Suppressor PEBP1/RKIP and Its Regulators in Breast Cancer Cells

Drug screening strategies focus on quantifying the phenotypic effects of different compounds on biological systems. High-throughput technologies have the potential to understand further the mechanisms by which these drugs produce the desired outcome. Reverse causal reasoning integrates existing biological knowledge and measurements of gene and protein abundances to infer their function. This approach can be employed to appraise the existing biological knowledge and data to prioritize targets for cancer therapies. We applied text mining and a manual literature search to extract known interactions between several metastasis suppressors and their regulators. We then identified the relevant interactions in the breast cancer cell line MCF7 using a knockdown dataset. We finally adopted a reverse causal reasoning approach to evaluate and prioritize pathways that are most consistent and responsive to drugs that inhibit cell growth. We evaluated this model in terms of agreement with the observations under treatment of several drugs that produced growth inhibition of cancer cell lines. In particular, we suggested that the metastasis suppressor PEBP1/RKIP is on the receiving end of two significant regulatory mechanisms. One involves RELA (transcription factor p65) and SNAI1, which were previously reported to inhibit PEBP1. The other involves the estrogen receptor (ESR1), which induces PEBP1 through the kinase NME1. Our model was derived in the specific context of breast cancer, but the observed responses to drug treatments were consistent in other cell lines. We further validated some of the predicted regulatory links in the breast cancer cell line MCF7 experimentally and highlighted the points of uncertainty in our model. To summarize, our model was consistent with the observed changes in activity with drug perturbations. In particular, two pathways, including PEBP1, were highly responsive and would be likely targets for intervention.

LINPS: a database for cancer-cell-specific perturbations of biological networks

Screening for potential cancer therapies using existing large datasets of drug perturbations requires expertise and resources not available to all. This is often a barrier for lab scientists to tap into these valuable resources. To address these issues, one can take advantage of prior knowledge especially those coded in standard formats such as causal biological networks (CBN). Large datasets can be converted into appropriate structures, analyzed once and the results made freely available in easy-to-use formats. We used the Library of Integrated Cellular Signatures to model the cell-specific effect of hundreds of drug treatments on gene expression. These signatures were then used to predict the effect of the treatments on several CBN using the network perturbation amplitudes analysis. We packaged the pre-computed scores in a database with an interactive web interface. The intuitive user-friendly interface can be used to query the database for drug perturbations and quantify their effect on multiple key biological functions in cancer cell lines. In addition to describing the process of building the database and the interface, we provide a realistic use case to explain how to use and interpret the results. To sum, we pre-computed cancer-cell-specific perturbation amplitudes of several biological networks and made the output available in a database with an interactive web interface.

Database URL https://mahshaaban.shinyapps.io/LINPSAPP/

Identification of Mucus-Associated Molecular Subtypes of Chronic Obstructive Pulmonary Disease: A Latent Profile Analysis Based on MUC5B-Associated Genes

BACKGROUND Chronic obstructive pulmonary disease (COPD) is a disease with high heterogeneity, which is a major challenge in clinical individualized treatment. A mucus phenotype is one of the main characteristics of COPD. MATERIAL AND METHODS Gene expression profiles of lung tissue samples were from the Lung Genomics Research Consortium. MUC5B-associated gene signatures were obtained based on a nonlinear feature screening algorithm. These signatures were used to fit a latent profile analysis (LPA) model to identify COPD molecular subtypes and build a subtype classifier to verify the subtypes. Then, we explored the characteristics of cilium assembly and beating signatures, transcriptome features, immune infiltration among the 3 subtypes by xCell, single-sample gene set enrichment analysis, network perturbation amplitude, and weighted gene co-expression network analysis algorithms. An external dataset was used to verify the above COPD subtypes. RESULTS Three subtypes associated with mucus were identified by LPA and verified in an external dataset. Subtype 1 displayed higher T helper type 1 (Th1) and basophil infiltration, higher Th17/regulatory T cells (Tregs) ratio, a higher level of cilium assembly and beating, and lower mast cell and Treg infiltration. The subtypes 2 and 3 demonstrated higher macrophage M2 infiltration in lung tissue, while subtype 3 had higher neutrophil and eosinophil infiltration than subtype 2. CONCLUSIONS Overall, this work identified 3 mucus-associated molecular subtypes related to MUC5B expression, which deepens the understanding of airway mucus secretion in COPD and potentially provides valuable information for precision therapy.

Respiratory Effects of Exposure to Aerosol From the Candidate Modified-Risk Tobacco Product THS 2.2 in an 18-Month Systems Toxicology Study With A/J Mice

Smoking cessation is the most effective measure for reducing the risk of smoking-related diseases. However, switching to less harmful products (modified-risk tobacco products [MRTP]) can be an alternative to help reduce the risk for adult smokers who would otherwise continue to smoke. In an 18-month chronic carcinogenicity/toxicity study in A/J mice (OECD Test Guideline 453), we assessed the aerosol of Tobacco Heating System 2.2 (THS 2.2), a candidate MRTP based on the heat-not-burn principle, compared with 3R4F cigarette smoke (CS). To capture toxicity- and disease-relevant mechanisms, we complemented standard toxicology endpoints with in-depth systems toxicology analyses. In this part of our publication series, we report on integrative assessment of the apical and molecular exposure effects on the respiratory tract (nose, larynx, and lungs). Across the respiratory tract, we found changes in inflammatory response following 3R4F CS exposure (eg, antimicrobial peptide response in the nose), with both shared and distinct oxidative and xenobiotic responses. Compared with 3R4F CS, THS 2.2 aerosol exerted far fewer effects on respiratory tract histology, including adaptive tissue changes in nasal and laryngeal epithelium and inflammation and emphysematous changes in the lungs. Integrative analysis of molecular changes confirmed the substantially lower impact of THS 2.2 aerosol than 3R4F CS on toxicologically and disease-relevant molecular processes such as inflammation, oxidative stress responses, and xenobiotic metabolism. In summary, this work exemplifies how apical and molecular endpoints can be combined effectively for toxicology assessment and further supports findings on the reduced respiratory health risks of THS 2.2 aerosol.

Thirteen-week nose-only inhalation exposures of propylene glycol aerosols in Sprague Dawley rats with a lung systems toxicology analysis

The objectives of this study were to increase PG exposure above concentrations tested by Suber et al. and use systems toxicology analysis of lung tissue to understand molecular events. Sprague Dawley rats were exposed to filtered air (sham), propylene glycol/water (PG/W; 90:10) or a propylene glycol/vegetable glycerin/water (PG/VG/W; 50:40:10) reference. The reference group was added at the high dose to observe any changes that might be associated with a carrier more in line with e-vapor products. Macroscopic examinations and terminal organ weights revealed no observations associated with exposure to PG/W or reference. Food consumption and body weights were unaffected by PG/W or reference when compared to sham. No exposure related alterations were observed in serum chemistry, hematology, coagulation, urinalysis or BALF cytology and clinical chemistry. Although clinical observations of dried red material around the nose in the high dose PG/W group were reported, histopathology showed no nasal hemorrhaging which was previously reported by Suber et al. Non-adverse PG/W and reference related findings of minimal mucous cell hyperplasia were noted in nasal cavity section II. No other exposure-related findings were noted in the primary or recovery necropsies. A systems toxicology analysis on lung tissue showed no statistically significant differentially expressed transcripts or proteins compared to the sham group. The endpoints measured from the PG/W high dose group did not differ significantly from those in the more common carrier PG/VG/W. As anticipated, exposure to PG aerosols was slightly irritating but well tolerated. Accordingly, the highest PG exposure (5 mg/L, 6 hrs/day) was regarded as the NOAEC, corresponding to a PG delivered dose of 1,152 mg/kg/day in rats.

Systems Toxicology Approach for Assessing Developmental Neurotoxicity in Larval Zebrafish

Adverse outcomes that result from chemical toxicity are rarely caused by dysregulation of individual proteins; rather, they are often caused by system-level perturbations in networks of molecular events. To fully understand the mechanisms of toxicity, it is necessary to recognize the interactions of molecules, pathways, and biological processes within these networks. The developing brain is a prime example of an extremely complex network, which makes developmental neurotoxicity one of the most challenging areas in toxicology. We have developed a systems toxicology method that uses a computable biological network to represent molecular interactions in the developing brain of zebrafish larvae. The network is curated from scientific literature and describes interactions between biological processes, signaling pathways, and adverse outcomes associated with neurotoxicity. This allows us to identify important signaling hubs, pathway interactions, and emergent adverse outcomes, providing a more complete understanding of neurotoxicity. Here, we describe the construction of a zebrafish developmental neurotoxicity network and its validation by integration with publicly available neurotoxicity-related transcriptomic datasets. Our network analysis identified consistent regulation of tumor suppressors p53 and retinoblastoma 1 (Rb1) as well as the oncogene Krüppel-like factor (Klf8) in response to chemically induced developmental neurotoxicity. The developed network can be used to interpret transcriptomic data in a neurotoxicological context.

Systems biology approach highlights mechanistic differences between Crohn's disease and ulcerative colitis

The molecular mechanisms of IBD have been the subject of intensive exploration. We, therefore, assembled the available information into a suite of causal biological network models, which offer comprehensive visualization of the processes underlying IBD. Scientific text was curated by using Biological Expression Language (BEL) and compiled with OpenBEL 3.0.0. Network properties were analysed by Cytoscape. Network perturbation amplitudes were computed to score the network models with transcriptomic data from public data repositories. The IBD network model suite consists of three independent models that represent signalling pathways that contribute to IBD. In the “intestinal permeability” model, programmed cell death factors were downregulated in CD and upregulated in UC. In the “inflammation” model, PPARG, IL6, and IFN-associated pathways were prominent regulatory factors in both diseases. In the “wound healing” model, factors promoting wound healing were upregulated in CD and downregulated in UC. Scoring of publicly available transcriptomic datasets onto these network models demonstrated that the IBD models capture the perturbation in each dataset accurately. The IBD network model suite can provide better mechanistic insights of the transcriptional changes in IBD and constitutes a valuable tool in personalized medicine to further understand individual drug responses in IBD.

A 6-month inhalation toxicology study in Apoe-/- mice demonstrates substantially lower effects of e-vapor aerosol compared with cigarette smoke in the respiratory tract

Cigarette smoking is the major cause of chronic obstructive pulmonary disease. Considerable attention has been paid to the reduced harm potential of nicotine-containing inhalable products such as electronic cigarettes (e-cigarettes). We investigated the effects of mainstream cigarette smoke (CS) and e-vapor aerosols (containing nicotine and flavor) generated by a capillary aerosol generator on emphysematous changes, lung function, and molecular alterations in the respiratory system of female Apoe-/- mice. Mice were exposed daily (3 h/day, 5 days/week) for 6 months to aerosols from three different e-vapor formulations-(1) carrier (propylene glycol and vegetable glycerol), (2) base (carrier and nicotine), or (3) test (base and flavor)-or to CS from 3R4F reference cigarettes. The CS and base/test aerosol concentrations were matched at 35 µg nicotine/L. CS exposure, but not e-vapor exposure, led to impairment of lung function (pressure-volume loop area, A and K parameters, quasi-static elastance and compliance) and caused marked lung inflammation and emphysematous changes, which were confirmed histopathologically and morphometrically. CS exposure caused lung transcriptome (activation of oxidative stress and inflammatory responses), lipidome, and proteome dysregulation and changes in DNA methylation; in contrast, these effects were substantially reduced in response to the e-vapor aerosol exposure. Compared with sham, aerosol exposure (carrier, base, and test) caused a slight impact on lung inflammation and epithelia irritation. Our results demonstrated that, in comparison with CS, e-vapor aerosols induced substantially lower biological and pathological changes in the respiratory tract associated with chronic inflammation and emphysema.

Impact of whole-body versus nose-only inhalation exposure systems on systemic, respiratory, and cardiovascular endpoints in a 2-month cigarette smoke exposure study in the ApoE-/- mouse model

Cigarette smoking is one major modifiable risk factor in the development and progression of chronic obstructive pulmonary disease and cardiovascular disease. To characterize and compare cigarette smoke (CS)-induced disease endpoints after exposure in either whole-body (WB) or nose-only (NO) exposure systems, we exposed apolipoprotein E-deficient mice to filtered air (Sham) or to the same total particulate matter (TPM) concentration of mainstream smoke from 3R4F reference cigarettes in NO or WB exposure chambers (EC) for 2 months. At matching TPM concentrations, we observed similar concentrations of carbon monoxide, acetaldehyde, and acrolein, but higher concentrations of nicotine and formaldehyde in NOEC than in WBEC. In both exposure systems, CS exposure led to the expected adaptive changes in nasal epithelia, altered lung function, lung inflammation, and pronounced changes in the nasal epithelial transcriptome and lung proteome. Exposure in the NOEC caused generally more severe histopathological changes in the nasal epithelia and a higher stress response as indicated by body weight decrease and lower blood lymphocyte counts compared with WB exposed mice. Erythropoiesis, and increases in total plasma triglyceride levels and atherosclerotic plaque area were observed only in CS-exposed mice in the WBEC group but not in the NOEC group. Although the composition of CS in the breathing zone is not completely comparable in the two exposure systems, the CS-induced respiratory disease endpoints were largely confirmed in both systems, with a higher magnitude of severity after NO exposure. CS-accelerated atherosclerosis and other pro-atherosclerotic factors were only significant in WBEC.

Group Behavior and Emergence of Cancer Drug Resistance

Drug resistance is a major impediment in cancer. Although it is generally thought that acquired drug resistance is due to genetic mutations, emerging evidence indicates that nongenetic mechanisms also play an important role. Resistance emerges through a complex interplay of clonal groups within a heterogeneous tumor and the surrounding microenvironment. Traits such as phenotypic plasticity, intercellular communication, and adaptive stress response, act in concert to ensure survival of intermediate reversible phenotypes, until permanent, resistant clones can emerge. Understanding the role of group behavior, and the underlying nongenetic mechanisms, can lead to more efficacious treatment designs and minimize or delay emergence of resistance.

Approaches to Test the Neurotoxicity of Environmental Contaminants in the Zebrafish Model: From Behavior to Molecular Mechanisms

The occurrence of neuroactive chemicals in the aquatic environment is on the rise and poses a potential threat to aquatic biota of currently unpredictable outcome. In particular, subtle changes caused by these chemicals to an organism's sensation or behavior are difficult to tackle with current test systems that focus on rodents or with in vitro test systems that omit whole-animal responses. In recent years, the zebrafish (Danio rerio) has become a popular model organism for toxicological studies and testing strategies, such as the standardized use of zebrafish early life stages in the Organisation for Economic Co-operation and Development's guideline 236. In terms of neurotoxicity, the zebrafish provides a powerful model to investigate changes to the nervous system from several different angles, offering the ability to tackle the mechanisms of action of chemicals in detail. The mechanistic understanding gained through the analysis of this model species provides a good basic knowledge of how neuroactive chemicals might interact with a teleost nervous system. Such information can help infer potential effects occurring to other species exposed to neuroactive chemicals in their aquatic environment and predicting potential risks of a chemical for the aquatic ecosystem. In the present article, we highlight approaches ranging from behavioral to structural, functional, and molecular analysis of the larval zebrafish nervous system, providing a holistic view of potential neurotoxic outcomes. Environ Toxicol Chem 2021;40:989-1006. © 2020 SETAC.

A 7-month inhalation toxicology study in C57BL/6 mice demonstrates reduced pulmonary inflammation and emphysematous changes following smoking cessation or switching to e-vapor products

Cigarette smoking causes serious diseases, including lung cancer, atherosclerotic coronary artery disease, peripheral vascular disease, chronic bronchitis, and emphysema. While cessation remains the most effective approach to minimize smoking-related disease, alternative non-combustible tobacco-derived nicotine-containing products may reduce disease risks among those unable or unwilling to quit. E-vapor aerosols typically contain significantly lower levels of smoke-related harmful and potentially harmful constituents; however, health risks of long-term inhalation exposures are unknown. We designed a 7-month inhalation study in C57BL/6 mice to evaluate long-term respiratory toxicity of e-vapor aerosols compared to cigarette smoke and to assess the impact of smoking cessation (Cessation group) or switching to an e-vapor product (Switching group) after 3 months of exposure to 3R4F cigarette smoke (CS). There were no significant changes in in-life observations (body weights, clinical signs) in e-vapor groups compared to the Sham Control. The 3R4F CS group showed reduced respiratory function during exposure and had lower body weight and showed transient signs of distress post-exposure. Following 7 months of exposure, e-vapor aerosols resulted in no or minimal increase in pulmonary inflammation, while exposure to 3R4F CS led to impairment of lung function and caused marked lung inflammation and emphysematous changes. Biological changes observed in the Switching group were similar to the Cessation group. 3R4F CS exposure dysregulated the lung and nasal tissue transcriptome, while these molecular effects were substantially lower in the e-vapor group. Results from this study demonstrate that in comparison with 3R4F CS, e-vapor aerosols induce substantially lower biological responses including pulmonary inflammation and emphysematous changes, and that complete switching from CS to e-vapor products significantly reduces biological changes associated with CS in C57BL/6 mice.

Comparison of the biological impact of aerosol of e-vapor device with MESH® technology and cigarette smoke on human bronchial and alveolar cultures

Exposure to aerosol from electronic vapor (e-vapor) products has been suggested to result in less risk of harm to smokers than cigarette smoke (CS) exposure. Although many studies on e-vapor products have tested the effects of liquid formulations on cell cultures, few have evaluated the effects of aerosolized formulations. We examined the effects of acute exposure to the aerosol of an e-vapor device that uses the MESH® technology (IQOS® MESH, Philip Morris International) and to CS from the 3R4F reference cigarette on human organotypic bronchial epithelial culture and alveolar triculture models. In contrast to 3R4F CS exposure, exposure to the IQOS MESH aerosol (Classic Tobacco flavor) did not cause cytotoxicity in bronchial epithelial cultures or alveolar tricultures despite its greater concentrations of deposited nicotine (3- and 4-fold, respectively). CS exposure caused a marked decrease in the frequency and active area of ciliary beating in bronchial cultures, whereas IQOS MESH aerosol exposure did not. Global mRNA expression and secreted protein profiles revealed a significantly lower impact of IQOS MESH aerosol exposure than 3R4F CS exposure. Overall, our whole aerosol exposure study shows a clearly reduced impact of IQOS MESH aerosol relative to CS in bronchial and alveolar cultures, even at greater nicotine doses.

Anatabine ameliorates intestinal inflammation and reduces the production of pro-inflammatory factors in a dextran sulfate sodium mouse model of colitis

Background

Inflammatory bowel disease (IBD) is the collective term for chronic immune-mediated diseases of unknown, multifactorial etiology, arising from the interplay between genetic and environmental factors and including two main disease manifestations: ulcerative colitis (UC) and Crohn’s disease. In the last few decades, naturally occurring alkaloids have gained interest because of their substantial anti-inflammatory effects in several animal models of disease. Studies on mouse models of IBD have demonstrated the anti-inflammatory action of the main tobacco alkaloid, nicotine. In addition, anatabine, a minor tobacco alkaloid also present in peppers, tomato, and eggplant presents anti-inflammatory properties in vivo and in vitro. In this study, we aimed to evaluate the anti-inflammatory properties of nicotine and anatabine in a dextran sulfate sodium (DSS) mouse model of UC.

Results

Oral administration of anatabine, but not nicotine, reduced the clinical symptoms of DSS-induced colitis. The result of gene expression analysis suggested that anatabine had a restorative effect on global DSS-induced gene expression profiles, while nicotine only had limited effects. Accordingly, MAP findings revealed that anatabine reduced the colonic abundance of DSS-associated cytokines and increased IL-10 abundance.

Conclusions

Our results support the amelioration of inflammatory effects by anatabine in the DSS mouse model of UC, and suggest that anatabine constitutes a promising therapeutic agent for IBD treatment.

Systems Toxicology Approach for Testing Chemical Cardiotoxicity in Larval Zebrafish

Transcriptomic approaches can give insight into molecular mechanisms underlying chemical toxicity and are increasingly being used as part of toxicological assessments. To aid the interpretation of transcriptomic data, we have developed a systems toxicology method that relies on a computable biological network model. We created the first network model describing cardiotoxicity in zebrafish larvae-a valuable emerging model species in testing cardiotoxicity associated with drugs and chemicals. The network is based on scientific literature and represents hierarchical molecular pathways that lead from receptor activation to cardiac pathologies. To test the ability of our approach to detect cardiotoxic outcomes from transcriptomic data, we have selected three publicly available data sets that reported chemically induced heart pathologies in zebrafish larvae for five different chemicals. Network-based analysis detected cardiac perturbations for four out of five chemicals tested, for two of them using transcriptomic data collected up to 3 days before the onset of a visible phenotype. Additionally, we identified distinct molecular pathways that were activated by the different chemicals. The results demonstrate that the proposed integrational method can be used for evaluating the effects of chemicals on the zebrafish cardiac function and, together with observed cardiac apical end points, can provide a comprehensive method for connecting molecular events to organ toxicity. The computable network model is freely available and may be used to generate mechanistic hypotheses and quantifiable perturbation values from any zebrafish transcriptomic data.

Multi-omics systems toxicology study of mouse lung assessing the effects of aerosols from two heat-not-burn tobacco products and cigarette smoke

•

Multi-omics systems toxicology study, comprising five omics data modalities.

•

Multi-Omics Factor Analysis and multi-modality functional network interpretation.

•

Cigarettes smoke (CS) induced complex immunoregulatory interactions across molecular layers.

•

Aerosols from two heat-not-burn tobacco products had less impact on lungs than CS.

Cigarette smoke (CS) causes adverse health effects and, for smoker who do not quit, modified risk tobacco products (MRTPs) can be an alternative to reduce the risk of developing smoking-related diseases. Standard toxicological endpoints can lack sensitivity, with systems toxicology approaches yielding broader insights into toxicological mechanisms. In a 6-month systems toxicology study on ApoE^−/− mice, we conducted an integrative multi-omics analysis to assess the effects of aerosols from the Carbon Heated Tobacco Product (CHTP) 1.2 and Tobacco Heating System (THS) 2.2—a potential and a candidate MRTP based on the heat-not-burn (HnB) principle—compared with CS at matched nicotine concentrations. Molecular exposure effects in the lungs were measured by mRNA/microRNA transcriptomics, proteomics, metabolomics, and lipidomics. Integrative data analysis included Multi-Omics Factor Analysis and multi-modality functional network interpretation. Across all five data modalities, CS exposure was associated with an increased inflammatory and oxidative stress response, and lipid/surfactant alterations. Upon HnB aerosol exposure these effects were much more limited or absent, with reversal of CS-induced effects upon cessation and switching to CHTP 1.2. Functional network analysis revealed CS-induced complex immunoregulatory interactions across the investigated molecular layers (e.g., itaconate, quinolinate, and miR-146) and highlighted the engagement of the heme–Hmox–bilirubin oxidative stress axis by CS. This work exemplifies how multi-omics approaches can be leveraged within systems toxicology studies and the generated multi-omics data set can facilitate the development of analysis methods and can yield further insights into the effects of toxicological exposures on the lung of mice.

The extraction of complex relationships and their conversion to biological expression language (BEL) overview of the BioCreative VI (2017) BEL track

Knowledge of the molecular interactions of biological and chemical entities and their involvement in biological processes or clinical phenotypes is important for data interpretation. Unfortunately, this knowledge is mostly embedded in the literature in such a way that it is unavailable for automated data analysis procedures. Biological expression language (BEL) is a syntax representation allowing for the structured representation of a broad range of biological relationships. It is used in various situations to extract such knowledge and transform it into BEL networks. To support the tedious and time-intensive extraction work of curators with automated methods, we developed the BEL track within the framework of BioCreative Challenges. Within the BEL track, we provide training data and an evaluation environment to encourage the text mining community to tackle the automatic extraction of complex BEL relationships. In 2017 BioCreative VI, the 2015 BEL track was repeated with new test data. Although only minor improvements in text snippet retrieval for given statements were achieved during this second BEL task iteration, a significant increase of BEL statement extraction performance from provided sentences could be seen. The best performing system reached a 32% F-score for the extraction of complete BEL statements and with the given named entities this increased to 49%. This time, besides rule-based systems, new methods involving hierarchical sequence labeling and neural networks were applied for BEL statement extraction.

The reduction of DSS-induced colitis severity in mice exposed to cigarette smoke is linked to immune modulation and microbial shifts

Exposure to cigarette smoke (CS) causes detrimental health effects, increasing the risk of cardiovascular, pulmonary diseases and carcinogenesis in exposed individuals. The impact of CS on Inflammatory Bowel Disease (IBD) has been established by a number of epidemiological and clinical studies. In fact, CS is associated with a higher risk of developing Crohn’s disease (CD) while inversely correlates with the development, disease risks, and relapse rate of ulcerative colitis (UC). To investigate the effect of CS exposure on experimental colitis, we performed a comprehensive and integrated comparative analysis of colon transcriptome and microbiome in mice exposed to dextran sodium sulfate (DSS) and CS. Colon transcriptome analysis revealed that CS downregulated specific pathways in a concentration-dependent manner, affecting both the inflammatory state and composition of the gut microbiome. Metagenomics analysis demonstrated that CS can modulate DSS-induced dysbiosis of specific bacterial genera, contributing to resolve the inflammation or accelerate recovery. The risks of smoking far outweigh any possible benefit, thus smoking cessation must always be encouraged because of its significant health benefits. However, the inverse association between active smoking and the development of UC cannot be ignored and the present study lays the foundation for investigating potential molecular mechanisms responsible for the attenuation of colitis by certain compounds of tobacco when decoupled from combustion.

Systems toxicology assessment of a representative e-liquid formulation using human primary bronchial epithelial cells

The development of reduced-risk products aims to provide alternatives to cigarettes that present less risk of harm for adult smokers. Responsible use of flavoring substances in these products may fulfill an important role in product acceptance. While most flavoring substances used in such products are also used by the food industry and are considered safe when ingested, their impact when inhaled may require further assessment. To aid in such an assessment, a three-step approach combining real-time cellular analysis, phenotypic high-content screening assays, and gene expression analysis was developed and tested in normal human bronchial epithelial cells with 28 flavoring substances commonly used in e-liquid formulations, dissolved individually or as a mixture in a base solution composed of propylene glycol, vegetable glycerin, and 0.6% nicotine. By employing this approach, we identified individual flavoring substances that potentially contribute greatly to the overall mixture effect (citronellol and alpha-pinene). By assessing modified mixtures, we showed that, although cytotoxic effects were found when assessed individually, alpha-pinene did not contribute to the overall mixture cytotoxicity. Most of the cytotoxic effect appeared to be attributable to citronellol, with the remaining substances contributing due to synergistic effects. We developed and used different scoring methods (Tox-Score, Phenotypic Score, and Biological Impact Factor/Network Perturbation Amplitude), ultimately enabling a ranking based on cytotoxicity, phenotypic outcome, and molecular network perturbations. This case study highlights the benefits of testing both individual flavoring substances and mixtures for e-liquid flavor assessment and emphasized the importance of data sharing for the benefit of consumer safety.

Application of a multi-layer systems toxicology framework for in vitro assessment of the biological effects of Classic Tobacco e-liquid and its corresponding aerosol using an e-cigarette device with MESH™ technology

We previously proposed a systems toxicology framework for in vitro assessment of e-liquids. The framework starts with the first layer aimed at screening the potential toxicity of e-liquids, followed by the second layer aimed at investigating the toxicity-related mechanism of e-liquids, and finally, the third layer aimed at evaluating the toxicity-related mechanism of the corresponding aerosols. In this work, we applied this framework to assess the impact of the e-liquid MESH Classic Tobacco and its aerosol compared with that of cigarette smoke (CS) from the 3R4F reference cigarette. In the first layer, we evaluated the cytotoxicity profile of the MESH Classic Tobacco e-liquid (containing humectants, nicotine, and flavors) and its Base e-liquid (containing humectant and nicotine only) in comparison with total particulate matter (TPM) of 3R4F CS using primary bronchial epithelial cell cultures. In the second layer, the same culture model was used to explore changes in specific markers using high-content screening assays to identify potential toxicity-related mechanisms induced by the MESH Classic Tobacco and Base e-liquids beyond cell viability in comparison with the 3R4F CS TPM-induced effects. Finally, in the third layer, we compared the impact of exposure to the MESH Classic Tobacco or Base aerosols with 3R4F CS using human organotypic air-liquid interface buccal and small airway epithelial cultures. The results showed that the cytotoxicity of the MESH Classic Tobacco liquid was similar to the Base liquid but lower than 3R4F CS TPM at comparable nicotine concentrations. Relative to 3R4F CS exposure, MESH Classic Tobacco aerosol exposure did not cause tissue damage and elicited lower changes in the mRNA, microRNA, and protein markers. In the context of tobacco harm reduction strategy, the framework is suitable to assess the potential-reduced impact of electronic cigarette aerosol relative to CS.

NPA: an R package for computing network perturbation amplitudes using gene expression data and two-layer networks

Background

High-throughput gene expression technologies provide complex datasets reflecting mechanisms perturbed in an experiment, typically in a treatment versus control design. Analysis of these information-rich data can be guided based on a priori knowledge, such as networks of related proteins or genes. Assessing the response of a specific mechanism and investigating its biological basis is extremely important in systems toxicology; as compounds or treatment need to be assessed with respect to a predefined set of key mechanisms that could lead to toxicity. Two-layer networks are suitable for this task, and a robust computational methodology specifically addressing those needs was previously published.

The NPA package (https://github.com/philipmorrisintl/NPA) implements the algorithm, and a data package of eight two-layer networks representing key mechanisms, such as xenobiotic metabolism, apoptosis, or epithelial immune innate activation, is provided.

Results

Gene expression data from an animal study are analyzed using the package and its network models. The functionalities are implemented using R6 classes, making the use of the package seamless and intuitive. The various network responses are analyzed using the leading node analysis, and an overall perturbation, called the Biological Impact Factor, is computed.

Conclusions

The NPA package implements the published network perturbation amplitude methodology and provides a set of two-layer networks encoded in the Biological Expression Language.

A lower impact of an acute exposure to electronic cigarette aerosols than to cigarette smoke in human organotypic buccal and small airway cultures was demonstrated using systems toxicology assessment

In the context of tobacco harm-reduction strategy, the potential reduced impact of electronic cigarette (EC) exposure should be evaluated relative to the impact of cigarette smoke exposure. We conducted a series of in vitro studies to compare the biological impact of an acute exposure to aerosols of "test mix" (flavors, nicotine, and humectants), "base" (nicotine and humectants), and "carrier" (humectants) formulations using MarkTen® EC devices with the impact of exposure to smoke of 3R4F reference cigarettes, at a matching puff number, using human organotypic air-liquid interface buccal and small airway cultures. We measured the concentrations of nicotine and carbonyls deposited in the exposure chamber after each exposure experiment. The deposited carbonyl concentrations were used as representative measures to assess the reduced exposure to potentially toxic volatile substances. We followed a systems toxicology approach whereby functional biological endpoints, such as histopathology and ciliary beating frequency, were complemented by multiplex and omics assays to measure secreted inflammatory proteins and whole-genome transcriptomes, respectively. Among the endpoints analyzed, the only parameters that showed a significant response to EC exposure were secretion of proteins and whole-genome transcriptomes. Based on the multiplex and omics analyzes, the cellular responses to EC aerosol exposure were tissue type-specific; however, those alterations were much smaller than those following cigarette smoke exposure, even when the EC aerosol exposure under the testing conditions resulted in a deposited nicotine concentration approximately 200 times that in saliva of EC users.

A six-month systems toxicology inhalation/cessation study in ApoE-/- mice to investigate cardiovascular and respiratory exposure effects of modified risk tobacco products, CHTP 1.2 and THS 2.2, compared with conventional cigarettes

Smoking is one of the major modifiable risk factors in the development and progression of chronic obstructive pulmonary disease (COPD) and cardiovascular disease (CVD). Modified-risk tobacco products (MRTP) are being developed to provide substitute products for smokers who are unable or unwilling to quit, to lessen the smoking-related health risks. In this study, the ApoE-/- mouse model was used to investigate the impact of cigarette smoke (CS) from the reference cigarette 3R4F, or aerosol from two potential MRTPs based on the heat-not-burn principle, carbon heated tobacco product 1.2 (CHTP1.2) and tobacco heating system 2.2 (THS 2.2), on the cardiorespiratory system over a 6-month period. In addition, cessation or switching to CHTP1.2 after 3 months of CS exposure was assessed. A systems toxicology approach combining physiology, histology and molecular measurements was used to evaluate the impact of MRTP aerosols in comparison to CS. CHTP1.2 and THS2.2 aerosols, compared with CS, demonstrated lower impact on the cardiorespiratory system, including low to absent lung inflammation and emphysematous changes, and reduced atherosclerotic plaque formation. Molecular analyses confirmed the lower engagement of pathological mechanisms by MRTP aerosols than CS. Both cessation and switching to CHTP1.2 reduced the observed CS effects to almost sham exposure levels.

Construction of a Suite of Computable Biological Network Models Focused on Mucociliary Clearance in the Respiratory Tract

Mucociliary clearance (MCC), considered as a collaboration of mucus secreted from goblet cells, the airway surface liquid layer, and the beating of cilia of ciliated cells, is the airways’ defense system against airborne contaminants. Because the process is well described at the molecular level, we gathered the available information into a suite of comprehensive causal biological network (CBN) models. The suite consists of three independent models that represent (1) cilium assembly, (2) ciliary beating, and (3) goblet cell hyperplasia/metaplasia and that were built in the Biological Expression Language, which is both human-readable and computable. The network analysis of highly connected nodes and pathways demonstrated that the relevant biology was captured in the MCC models. We also show the scoring of transcriptomic data onto these network models and demonstrate that the models capture the perturbation in each dataset accurately. This work is a continuation of our approach to use computational biological network models and mathematical algorithms that allow for the interpretation of high-throughput molecular datasets in the context of known biology. The MCC network model suite can be a valuable tool in personalized medicine to further understand heterogeneity and individual drug responses in complex respiratory diseases.

In vitro systems toxicology-based assessment of the potential modified risk tobacco product CHTP 1.2 for vascular inflammation- and cytotoxicity-associated mechanisms promoting adhesion of monocytic cells to human coronary arterial endothelial cells

Cigarette smoking causes cardiovascular diseases. Heating tobacco instead of burning it reduces the amount of toxic compounds in the aerosol and may exert a reduced impact on health compared with cigarette smoke. Aqueous extract from the aerosol of a potential modified risk tobacco product, the Carbon Heated Tobacco Product (CHTP) 1.2, was compared in vitro with aqueous extract from the smoke of a 3R4F reference cigarette for its impact on the adhesion of monocytic cells to artery endothelial cells. Human coronary artery endothelial cells (HCAEC) were treated for 4 h with conditioned media from human monocytic Mono Mac 6 (MM6) cells exposed to CHTP1.2 or 3R4F extracts for 2 h or directly with those extracts freshly generated. In vitro monocyte-endothelial cell adhesion was measured concomitantly with inflammatory, oxidative stress, cytotoxicity, and death markers. Furthermore, transcriptomics analyses enabled to quantify the level of perturbation in HCAECs, and provide biological interpretation for the underlying molecular changes following exposure to 3R4F or CHTP1.2 extract. Our systems toxicology study demonstrated that approximately 10-15-fold higher concentrations of the CHTP 1.2 aerosol extract were needed to elicit similar effects as the 3R4F smoke extract on cardiovascular disease-relevant inflammation and cytotoxicity-related mechanisms and markers investigated in vitro.

Comparative biological impacts of an aerosol from carbon-heated tobacco and smoke from cigarettes on human respiratory epithelial cultures: A systems toxicology assessment

The biological impact of an aerosol of a potential modified-risk tobacco product, carbon heated tobacco product 1.2 (CHTP1.2), was comprehensively assessed for the first time in vitro using human small airway and nasal epithelial models following a systems toxicology approach. The potentially reduced effects of CHTP1.2 aerosol exposure were benchmarked against those of 3R4F cigarette smoke at similar nicotine concentrations. Experimental repetitions were conducted for which new batches of small airway and nasal cultures were exposed to CHTP1.2 aerosol or 3R4F smoke for 28 minutes. The biological impacts were determined based on a collection of endpoints including morphology, cytotoxicity, proinflammatory mediator profiles, cytochrome P450 1A1/1B1 activity, global mRNA and microRNA changes and proteome profiles. Alterations in mRNA expression were detected in cultures exposed to CHTP1.2 aerosol, without noticeable morphological changes and cytotoxicity, and minimal impact on proinflammatory mediator and proteome profiles. The changes linked to CHTP1.2 aerosol exposure, when observed, were transient. However, the impact of 3R4F smoke exposure persisted long post-exposure and greater than CHTP1.2 aerosol. Morphological changes were observed only in cultures exposed to 3R4F smoke. The lower biological effects of CHTP1.2 aerosol than 3R4F smoke exposure were observed similarly in both small airway and nasal epithelial cultures.

Assessment of the impact of aerosol from a potential modified risk tobacco product compared with cigarette smoke on human organotypic oral epithelial cultures under different exposure regimens

Cigarette smoke (CS) is affecting considerably the oral mucosa. Heating, instead of burning, tobacco reduces consistently the amount of toxic compounds and may exert a lower impact on oral health than combusted cigarettes. The carbon-heated tobacco product 1.2 (CHTP1.2) is a potential modified risk tobacco product (MRTP) based on heat-not-burn technology. Using a systems toxicology assessment framework, we compared the effects of exposure to CHTP1.2 aerosol with those of CS from a reference cigarette (3R4F). Human organotypic cultures derived from buccal and gingival epithelia were exposed acutely (28-min) or repeatedly (28 min/day for 3 days), respectively, to two matching concentrations of CHTP1.2 aerosol or 3R4F CS, and a non-diluted (100%) CHTP1.2 aerosol. The results showed an absence of cytotoxicity, reduction in pathophysiological alterations, toxicological marker proteins, and inflammatory mediators following exposure to CHTP1.2 aerosol compared with 3R4F CS. Changes in mRNA and miRNA expression were linked by an integrative analysis approach, suggesting a regulatory role of miRNAs in several smoke/disease-relevant biological processes induced by 3R4F CS. The identification of mechanisms by which potential MRTPs can reduce the impact of tobacco use on biological systems is of great importance in understanding the molecular basis of the smoking harm reduction paradigm.

A 90-day OECD TG 413 rat inhalation study with systems toxicology endpoints demonstrates reduced exposure effects of the aerosol from the carbon heated tobacco product version 1.2 (CHTP1.2) compared with cigarette smoke. II. Systems toxicology assessment

Modified risk tobacco products (MRTPs) have the potential to reduce smoking-related health risks. The Carbon Heated Tobacco Product 1.2 (CHTP1.2) is a potential MRTP that uses a pressed carbon heat source to generate an aerosol by heating tobacco. Here, we report the results from the systems toxicology arm of a 90-day rat inhalation study (OECD test guideline 413) to assess the effects of CHTP1.2 aerosol compared with cigarette smoke (CS). Transcriptomics, proteomics, and lipidomics analyses complemented the standard endpoints. In the respiratory nasal epithelium, CS induced an adaptive tissue and inflammatory response, which was much weaker after CHTP1.2 aerosol exposure, mostly limited to the highest CHTP1.2 concentration (at twice the 3R4F CS concentration: 50 vs. 23 μg nicotine/L), in female rats. In the lungs, the effects of CS exposure included inflammatory and cellular stress responses, which were absent or much lower after CHTP1.2 aerosol exposure. Outside of the respiratory tract, CS and CHTP1.2 aerosol induced effects that were previously associated with exposure to any nicotine-containing aerosol, e.g., lower lipid concentrations in serum. Overall, this systems toxicology analysis complements and confirms the results from classical toxicological endpoints and further suggests potentially reduced respiratory health risks of CHTP1.2.

Network-based technologies for early drug discovery

Although the traditional drug discovery approach has led to the development of many successful drugs, the attrition rates remain high. Recent advances in systems-oriented approaches (systems-biology and/or pharmacology) and 'omics technologies has led to a plethora of new computational tools that promise to enable a more-informed and successful implementation of the reductionist, one drug for one target for one disease, approach. These tools, based on biomolecular pathways and interaction networks, offer a systematic approach to unravel the mechanism(s) of a disease and link them to the chemical space and network footprint of a drug. Drug discovery can draw upon this holistic approach to identify the most-promising targets and compounds during the early phases of development.

Comparative effects of a candidate modified-risk tobacco product Aerosol and cigarette smoke on human organotypic small airway cultures: a systems toxicology approach

Using an in vitro human small airway epithelium model, we assessed the biological impact of an aerosol from a candidate modified-risk tobacco product, the tobacco heating system (THS) 2.2, to investigate the potential reduced risk of THS2.2 aerosol exposure compared with cigarette smoke. Following the recommendations of the Institute of Medicine and the Tobacco Product Assessment Consortium, in which modified-risk tobacco products assessment should be performed in comparison with standard conventional products, the effects of the THS2.2 aerosol exposure on the small airway cultures were compared with those of 3R4F cigarette smoke. We used a systems toxicology approach whereby elucidation of toxic effects is derived not only from functional assay readouts but also from omics technologies. Cytotoxicity, ciliary beating function, secretion of pro-inflammatory mediators and histological assessment represented functional assays. The omics data included transcriptomic and miRNA profiles. Exposure-induced perturbations of causal biological networks were computed from the transcriptomic data. The results showed that THS2.2 aerosol exposure at the tested doses elicited lower cytotoxicity levels and lower changes in the secreted pro-inflammatory mediators than 3R4F smoke. Although THS2.2 exposure elicited alterations in the gene expression, a higher transcriptome-induced biological impact was observed following 3R4F smoke: The effects of THS2.2 aerosol exposure, if observed, were mostly transient and diminished more rapidly after exposure than those of 3R4F smoke. The study demonstrated that the systems toxicology approach can reveal changes at the cellular level that would be otherwise not detected from functional assays, thus increasing the sensitivity to detect potential toxicity of a treatment/exposure.

Systems toxicology meta-analysis of in vitro assessment studies: biological impact of a candidate modified-risk tobacco product aerosol compared with cigarette smoke on human organotypic cultures of the aerodigestive tract

Systems biology combines comprehensive molecular analyses with quantitative modeling to understand the characteristics of a biological system as a whole. Leveraging a similar approach, systems toxicology aims to decipher complex biological responses following exposures. This work reports a systems toxicology meta-analysis in the context of in vitro assessment of a candidate modified-risk tobacco product (MRTP) using three human organotypic cultures of the aerodigestive tract (buccal, bronchial, and nasal epithelia). Complementing a series of functional measures, a causal network enrichment analysis of transcriptomic data was used to compare quantitatively the biological impact of aerosol from the Tobacco Heating System (THS) 2.2, a candidate MRTP, with 3R4F cigarette smoke (CS) at similar nicotine concentrations. Lower toxicity was observed in all cultures following exposure to THS2.2 aerosol compared with 3R4F CS. Because of their morphological differences, a smaller exposure impact was observed in the buccal (stratified epithelium) compared with the bronchial and nasal (pseudostratified epithelium). However, the causal network enrichment approach supported a similar mechanistic impact of CS across the three cultures, including the impact on xenobiotic, oxidative stress, and inflammatory responses. At comparable nicotine concentrations, THS2.2 aerosol elicited reduced and more transient effects on these processes. To demonstrate the benefits of additional data modalities, we employed a newly established targeted mass-spectrometry marker panel to further confirm the reduced cellular stress responses elicited by THS2.2 aerosol compared with 3R4F CS in the nasal culture. Overall, this work demonstrates the applicability and robustness of the systems toxicology approach for in vitro inhalation toxicity assessment.

Logic programming reveals alteration of key transcription factors in multiple myeloma

Innovative approaches combining regulatory networks (RN) and genomic data are needed to extract biological information for a better understanding of diseases, such as cancer, by improving the identification of entities and thereby leading to potential new therapeutic avenues. In this study, we confronted an automatically generated RN with gene expression profiles (GEP) from a cohort of multiple myeloma (MM) patients and normal individuals using global reasoning on the RN causality to identify key-nodes. We modeled each patient by his or her GEP, the RN and the possible automatically detected repairs needed to establish a coherent flow of the information that explains the logic of the GEP. These repairs could represent cancer mutations leading to GEP variability. With this reasoning, unmeasured protein states can be inferred, and we can simulate the impact of a protein perturbation on the RN behavior to identify therapeutic targets. We showed that JUN/FOS and FOXM1 activities are altered in almost all MM patients and identified two survival markers for MM patients. Our results suggest that JUN/FOS-activation has a strong impact on the RN in view of the whole GEP, whereas FOXM1-activation could be an interesting way to perturb an MM subgroup identified by our method.

Novel approaches to develop community-built biological network models for potential drug discovery

INTRODUCTION

Hundreds of thousands of data points are now routinely generated in clinical trials by molecular profiling and NGS technologies. A true translation of this data into knowledge is not possible without analysis and interpretation in a well-defined biology context. Currently, there are many public and commercial pathway tools and network models that can facilitate such analysis. At the same time, insights and knowledge that can be gained is highly dependent on the underlying biological content of these resources. Crowdsourcing can be employed to guarantee the accuracy and transparency of the biological content underlining the tools used to interpret rich molecular data. Areas covered: In this review, the authors describe crowdsourcing in drug discovery. The focal point is the efforts that have successfully used the crowdsourcing approach to verify and augment pathway tools and biological network models. Technologies that enable the building of biological networks with the community are also described. Expert opinion: A crowd of experts can be leveraged for the entire development process of biological network models, from ontologies to the evaluation of their mechanistic completeness. The ultimate goal is to facilitate biomarker discovery and personalized medicine by mechanistically explaining patients' differences with respect to disease prevention, diagnosis, and therapy outcome.

Mechanistic Evaluation of the Impact of Smoking and Chronic Obstructive Pulmonary Disease on the Nasal Epithelium

Chronic obstructive pulmonary disease (COPD) is one of the major causes of chronic morbidity and mortality worldwide. The development of markers of COPD onset is hampered by the lack of accessibility to the primary target tissue, and there is a need to consider other sample sources as surrogates for biomarker research. Airborne toxicants pass through the nasal epithelium before reaching the lower airways, and the similarity with bronchial histology makes it an attractive surrogate for lower airways. In this work, we describe the transcriptomics findings from the nasal epithelia of subjects enrolled in a clinical study focusing on the identification of COPD biomarkers. Transcriptomic data were analyzed using the biological network approach that enabled us to pinpoint the biological processes affected in the upper respiratory tract in response to smoking and mild-to-moderate COPD. Our results indicated that nasal and lower airway immune responses were considerably different in COPD subjects and caution should be exercised when using upper airway samples as a surrogate for the lower airway. Nevertheless, the network approach described here could present a sensitive means of identifying smokers at risk of developing COPD.

A systems toxicology approach for comparative assessment: Biological impact of an aerosol from a candidate modified-risk tobacco product and cigarette smoke on human organotypic bronchial epithelial cultures

This study reports a comparative assessment of the biological impact of a heated tobacco aerosol from the tobacco heating system (THS) 2.2 and smoke from a combustible 3R4F cigarette. Human organotypic bronchial epithelial cultures were exposed to an aerosol from THS2.2 (a candidate modified-risk tobacco product) or 3R4F smoke at similar nicotine concentrations. A systems toxicology approach was applied to enable a comprehensive exposure impact assessment. Culture histology, cytotoxicity, secreted pro-inflammatory mediators, ciliary beating, and genome-wide mRNA/miRNA profiles were assessed at various time points post-exposure. Series of experimental repetitions were conducted to increase the robustness of the assessment. At similar nicotine concentrations, THS2.2 aerosol elicited lower cytotoxicity compared with 3R4F smoke. No morphological change was observed following exposure to THS2.2 aerosol, even at nicotine concentration three times that of 3R4F smoke. Lower levels of secreted mediators and fewer miRNA alterations were observed following exposure to THS2.2 aerosol than following 3R4F smoke. Based on the computational analysis of the gene expression changes, 3R4F (0.13 mg nicotine/L) elicited the highest biological impact (100%) in the context of Cell Fate, Cell Proliferation, Cell Stress, and Inflammatory Network Models at 4 h post-exposure. Whereas, the corresponding impact of THS2.2 (0.14 mg nicotine/L) was 7.6%.

Comparative systems toxicology analysis of cigarette smoke and aerosol from a candidate modified risk tobacco product in organotypic human gingival epithelial cultures: A 3-day repeated exposure study

Smoking is one of the major lifestyle-related risk factors for periodontal diseases. Modified risk tobacco products (MRTP) offer a promising alternative in the harm reduction strategy for adult smokers unable to quit. Using a systems toxicology approach, we investigated and compared the exposure effects of a reference cigarette (3R4F) and a heat-not-burn technology-based candidate MRTP, the Tobacco Heating System (THS) 2.2. Human gingival epithelial organotypic cultures were repeatedly exposed (3 days) for 28 min at two matching concentrations of cigarette smoke (CS) or THS2.2 aerosol. Results showed only minor histopathological alterations and minimal cytotoxicity upon THS2.2 aerosol exposure compared to CS (1% for THS2.2 aerosol vs. 30% for CS, at the high concentration). Among the 14 proinflammatory mediators analyzed, only 5 exhibited significant alterations with THS2.2 exposure compared with 11 upon CS exposure. Transcriptomic and metabolomic analysis indicated a general reduction of the impact in THS2.2 aerosol-exposed samples with respect to CS (∼79% lower biological impact for the high THS2.2 aerosol concentration compared to CS, and 13 metabolites significantly perturbed for THS2.2 vs. 181 for CS). This study indicates that exposure to THS2.2 aerosol had a lower impact on the pathophysiology of human gingival organotypic cultures than CS.

Evaluation of the Tobacco Heating System 2.2 (THS2.2). Part 5: microRNA expression from a 90-day rat inhalation study indicates that exposure to THS2.2 aerosol causes reduced effects on lung tissue compared with cigarette smoke

Modified-risk tobacco products (MRTP) are designed to reduce the individual risk of tobacco-related disease as well as population harm compared to smoking cigarettes. Experimental proof of their benefit needs to be provided at multiple levels in research fields. Here, we examined microRNA (miRNA) levels in the lungs of rats exposed to a candidate modified-risk tobacco product, the Tobacco Heating System 2.2 (THS2.2) in a 90-day OECD TG-413 inhalation study. Our aim was to assess the miRNA response to THS2.2 aerosol compared with the response to combustible cigarettes (CC) smoke from the reference cigarette 3R4F. CC smoke exposure, but not THS2.2 aerosol exposure, caused global miRNA downregulation, which may be explained by the interference of CC smoke constituents with the miRNA processing machinery. Upregulation of specific miRNA species, such as miR-146a/b and miR-182, indicated that they are causal elements in the inflammatory response in CC-exposed lungs, but they were reduced after THS2.2 aerosol exposure. Transforming transcriptomic data into protein activity based on corresponding downstream gene expression, we identified potential mechanisms for miR-146a/b and miR-182 that were activated by CC smoke but not by THS2.2 aerosol and possibly involved in the regulation of those miRNAs. The inclusion of miRNA profiling in systems toxicology approaches increases the mechanistic understanding of the complex exposure responses.