Functional analysis of OMICs data and small molecule compounds in an integrated "knowledge-based" platform

Bayesian-frequentist hybrid inference framework for single cell RNA-seq analyses

BACKGROUND

Single cell RNA sequencing technology (scRNA-seq) has been proven useful in understanding cell-specific disease mechanisms. However, identifying genes of interest remains a key challenge. Pseudo-bulk methods that pool scRNA-seq counts in the same biological replicates have been commonly used to identify differentially expressed genes. However, such methods may lack power due to the limited sample size of scRNA-seq datasets, which can be prohibitively expensive.

RESULTS

Motivated by this, we proposed to use the Bayesian-frequentist hybrid (BFH) framework to increase the power and we showed in simulated scenario, the proposed BFH would be an optimal method when compared with other popular single cell differential expression methods if both FDR and power were considered. As an example, the method was applied to an idiopathic pulmonary fibrosis (IPF) case study.

CONCLUSION

In our IPF example, we demonstrated that with a proper informative prior, the BFH approach identified more genes of interest. Furthermore, these genes were reasonable based on the current knowledge of IPF. Thus, the BFH offers a unique and flexible framework for future scRNA-seq analyses.

Bruton's Tyrosine Kinase (BTK) Inhibitors and Autoimmune Diseases: Making Sense of BTK Inhibitor Specificity Profiles and Recent Clinical Trial Successes and Failures

Clinical development of BTK kinase inhibitors for treating autoimmune diseases has lagged behind development of these drugs for treating cancers, due in part from concerns over the lack of selectivity and associated toxicity profiles of first generation drug candidates when used in the long term treatment of immune mediated diseases. Second generation BTK inhibitors have made great strides in limiting off-target activities for distantly related kinases, though they have had variable success at limiting cross-reactivity within the more closely related TEC family of kinases. We investigated the BTK specificity and toxicity profiles, drug properties, disease associated signaling pathways, clinical indications, and trial successes and failures for the 13 BTK inhibitor drug candidates tested in phase 2 or higher clinical trials representing 7 autoimmune and 2 inflammatory immune-mediated diseases. We focused on rheumatoid arthritis (RA), multiple sclerosis (MS), and systemic lupus erythematosus (SLE) where the majority of BTK nonclinical and clinical studies have been reported, with additional information for pemphigus vulgaris (PV), Sjogren’s disease (SJ), chronic spontaneous urticaria (CSU), graft versus host disease (GVHD), and asthma included where available. While improved BTK selectivity versus kinases outside the TEC family improved clinical toxicity profiles, less profile distinction was evident within the TEC family. Analysis of genetic associations of RA, MS, and SLE biomarkers with TEC family members revealed that BTK and TEC family members may not be drivers of disease. They are, however, mediators of signaling pathways associated with the pathophysiology of autoimmune diseases. BTK in particular may be associated with B cell and myeloid differentiation as well as autoantibody development implicated in immune mediated diseases. Successes in the clinic for treating RA, MS, PV, ITP, and GVHD, but not for SLE and SJ support the concept that BTK plays an important role in mediating pathogenic processes amenable to therapeutic intervention, depending on the disease. Based on the data collected in this study, we propose that current compound characteristics of BTK inhibitor drug candidates for the treatment of autoimmune diseases have achieved the selectivity, safety, and coverage requirements necessary to deliver therapeutic benefit.

Towards realizing the vision of precision medicine: AI based prediction of clinical drug response

Accurate and individualized prediction of response to therapies is central to precision medicine. However, because of the generally complex and multifaceted nature of clinical drug response, realizing this vision is highly challenging, requiring integrating different data types from the same individual into one prediction model. We used the anti-epileptic drug brivaracetam as a case study and combine a hybrid data/knowledge-driven feature extraction with machine learning to systematically integrate clinical and genetic data from a clinical discovery dataset (n = 235 patients). We constructed a model that successfully predicts clinical drug response [area under the curve (AUC) = 0.76] and show that even with limited sample size, integrating high-dimensional genetics data with clinical data can inform drug response prediction. After further validation on data collected from an independently conducted clinical study (AUC = 0.75), we extensively explore our model to gain insights into the determinants of drug response, and identify various clinical and genetic characteristics predisposing to poor response. Finally, we assess the potential impact of our model on clinical trial design and demonstrate that, by enriching for probable responders, significant reductions in clinical study sizes may be achieved. To our knowledge, our model represents the first retrospectively validated machine learning model linking drug mechanism of action and the genetic, clinical and demographic background in epilepsy patients to clinical drug response. Hence, it provides a blueprint for how machine learning-based multimodal data integration can act as a driver in achieving the goals of precision medicine in fields such as neurology.

Linking Coregulated Gene Modules with Polycyclic Aromatic Hydrocarbon-Related Cancer Risk in the 3D Human Bronchial Epithelium

Exposure to polycyclic aromatic hydrocarbons (PAHs) often occurs as complex chemical mixtures, which are linked to numerous adverse health outcomes in humans, with cancer as the greatest concern. The cancer risk associated with PAH exposures is commonly evaluated using the relative potency factor (RPF) approach, which estimates PAH mixture carcinogenic potential based on the sum of relative potency estimates of individual PAHs, compared to benzo[a]pyrene (BAP), a reference carcinogen. The present study evaluates molecular mechanisms related to PAH cancer risk through integration of transcriptomic and bioinformatic approaches in a 3D human bronchial epithelial cell model. Genes with significant differential expression from human bronchial epithelium exposed to PAHs were analyzed using a weighted gene coexpression network analysis (WGCNA) two-tiered approach: first to identify gene sets comodulated to RPF and second to link genes to a more comprehensive list of regulatory values, including inhalation-specific risk values. Over 3000 genes associated with processes of cell cycle regulation, inflammation, DNA damage, and cell adhesion processes were found to be comodulated with increasing RPF with pathways for cell cycle S phase and cytoskeleton actin identified as the most significantly enriched biological networks correlated to RPF. In addition, comodulated genes were linked to additional cancer-relevant risk values, including inhalation unit risks, oral cancer slope factors, and cancer hazard classifications from the World Health Organization's International Agency for Research on Cancer (IARC). These gene sets represent potential biomarkers that could be used to evaluate cancer risk associated with PAH mixtures. Among the values tested, RPF values and IARC categorizations shared the most similar responses in positively and negatively correlated gene modules. Together, we demonstrated a novel manner of integrating gene sets with chemical toxicity equivalence estimates through WGCNA to understand potential mechanisms.

E7386, a Selective Inhibitor of the Interaction between β-Catenin and CBP, Exerts Antitumor Activity in Tumor Models with Activated Canonical Wnt Signaling

The Wnt/β-catenin signaling pathway plays crucial roles in embryonic development and the development of multiple types of cancer, and its aberrant activation provides cancer cells with escape mechanisms from immune checkpoint inhibitors. E7386, an orally active selective inhibitor of the interaction between β-catenin and CREB binding protein, which is part of the Wnt/β-catenin signaling pathway, disrupts the Wnt/β-catenin signaling pathway in HEK293 and adenomatous polyposis coli (APC)-mutated human gastric cancer ECC10 cells. It also inhibited tumor growth in an ECC10 xenograft model and suppressed polyp formation in the intestinal tract of Apc^{Min /+} mice, in which mutation of Apc activates the Wnt/β-catenin signaling pathway. E7386 demonstrated antitumor activity against mouse mammary tumors developed in mouse mammary tumor virus (MMTV)-Wnt1 transgenic mice. Gene expression profiling using RNA sequencing data of MMTV-Wnt1 tumor tissue from mice treated with E7386 showed that E7386 downregulated genes in the hypoxia signaling pathway and immune responses related to the CCL2, and IHC analysis showed that E7386 induced infiltration of CD8⁺ cells into tumor tissues. Furthermore, E7386 showed synergistic antitumor activity against MMTV-Wnt1 tumor in combination with anti-PD-1 antibody. In conclusion, E7386 demonstrates clear antitumor activity via modulation of the Wnt/β-catenin signaling pathway and alteration of the tumor and immune microenvironments, and its antitumor activity can be enhanced in combination with anti-PD-1 antibody. SIGNIFICANCE: These findings demonstrate that the novel anticancer agent, E7386, modulates Wnt/β-catenin signaling, altering the tumor immune microenvironment and exhibiting synergistic antitumor activity in combination with anti-PD-1 antibody.

Classifying polycyclic aromatic hydrocarbons by carcinogenic potency using in vitro biosignatures

One of the most difficult challenges for risk assessment is evaluation of chemicals that predominately co-occur in mixtures like polycyclic aromatic hydrocarbons (PAHs). We previously developed a classification model in which systems biology data collected from mice short-term after chemical exposure accurately predict tumor outcome. The present study demonstrates translation of this approach into a human in vitro model in which chemical-specific bioactivity profiles from 3D human bronchial epithelial cells (HBEC) classify PAHs by carcinogenic potency. Gene expression profiles were analyzed from HBEC exposed to carcinogenic and non-carcinogenic PAHs and classification accuracies were identified for individual pathway-based gene sets. Posterior probabilities of best performing gene sets were combined via Bayesian integration resulting in a classifier with four gene sets, including aryl hydrocarbon receptor signaling, regulation of epithelial mesenchymal transition, regulation of angiogenesis, and cell cycle G2-M. In addition, transcriptional benchmark dose modeling of benzo[a]pyrene (BAP) showed that the most sensitive gene sets to BAP regulation were largely dissimilar from those that best classified PAH carcinogenicity challenging current assumptions that BAP carcinogenicity (and subsequent mode of action) is reflective of overall PAH carcinogenicity. These results illustrate utility of using systems toxicology approaches to analyze global gene expression towards carcinogenic hazard assessment.

Chemocentric Informatics Analysis: Dexamethasone Versus Combination Therapy for COVID-19

COVID-19 is a biphasic infectious disease with no approved vaccine or pharmacotherapy. The first drug that has shown promise in reducing COVID-19 mortality in severely-ill patients is dexamethasone, a cheap, well-known anti-inflammatory glucocorticoid, approved for the treatment of inflammatory conditions including respiratory diseases such as asthma and tuberculosis. However, about 80% of COVID-19 patients requiring oxygenation, and about 67% of patients on ventilators, are not responsive to dexamethasone therapy mainly. Additionally, using higher doses of dexamethasone for prolonged periods of time can lead to severe side effects and some patients may develop corticosteroid resistance leading to treatment failure. In order to increase the therapeutic efficacy of dexamethasone in COVID-19 patients, while minimizing dexamethasone-related complications that could result from using higher doses of the drug, we applied a chemocentric informatics approach to identify combination therapies. Our results indicated that combining dexamethasone with fast long-acting beta-2 adrenergic agonists (LABAs), such as formoterol and salmeterol, can ease respiratory symptoms hastily, until dexamethasone's anti-inflammatory and immunosuppressant effects kick in. Our studies demonstrated that LABAs and dexamethasone (or other glucocorticoids) exert synergistic effects that will augment both anti-inflammatory and fibronectin-mediated anticoagulant effects. We also propose other alternatives to LABAs that are supported by sound systems biology evidence, such as nitric oxide. Other drugs such as sevoflurane and treprostinil interact with the SARS-CoV-2 interactome and deserve further exploration. Moreover, our chemocentric informatics approach provides systems biology evidence that combination therapies for COVID-19 will have higher chances of perturbing the SARS-CoV-2 human interactome, which may negatively impact COVID-19 disease pathways.

A Systems Biology Workflow for Drug and Vaccine Repurposing: Identifying Small-Molecule BCG Mimics to Reduce or Prevent COVID-19 Mortality

PURPOSE

Coronavirus disease 2019 (COVID-19) is expected to continue to cause worldwide fatalities until the World population develops 'herd immunity', or until a vaccine is developed and used as a prevention. Meanwhile, there is an urgent need to identify alternative means of antiviral defense. Bacillus Calmette-Guérin (BCG) vaccine that has been recognized for its off-target beneficial effects on the immune system can be exploited to boast immunity and protect from emerging novel viruses.

METHODS

We developed and employed a systems biology workflow capable of identifying small-molecule antiviral drugs and vaccines that can boast immunity and affect a wide variety of viral disease pathways to protect from the fatal consequences of emerging viruses.

RESULTS

Our analysis demonstrates that BCG vaccine affects the production and maturation of naïve T cells resulting in enhanced, long-lasting trained innate immune responses that can provide protection against novel viruses. We have identified small-molecule BCG mimics, including antiviral drugs such as raltegravir and lopinavir as high confidence hits. Strikingly, our top hits emetine and lopinavir were independently validated by recent experimental findings that these compounds inhibit the growth of SARS-CoV-2 in vitro.

CONCLUSIONS

Our results provide systems biology support for using BCG and small-molecule BCG mimics as putative vaccine and drug candidates against emergent viruses including SARS-CoV-2.

TRPS1 Is a Lineage-Specific Transcriptional Dependency in Breast Cancer

Perturbed epigenomic programs play key roles in tumorigenesis, and chromatin modulators are candidate therapeutic targets in various human cancer types. To define singular and shared dependencies on DNA and histone modifiers and transcription factors in poorly differentiated adult and pediatric cancers, we conducted a targeted shRNA screen across 59 cell lines of 6 cancer types. Here, we describe the TRPS1 transcription factor as a strong breast cancer-specific hit, owing largely to lineage-restricted expression. Knockdown of TRPS1 resulted in perturbed mitosis, apoptosis, and reduced tumor growth. Integrated analysis of TRPS1 transcriptional targets, chromatin binding, and protein interactions revealed that TRPS1 is associated with the NuRD repressor complex. These findings uncover a transcriptional network that is essential for breast cancer cell survival and propagation.

Witwicki et al. use a targeted shRNA screening strategy to identify transcriptional and epigenomic dependencies in poorly differentiated human cancers. TRPS1 is a lineage-specific transcription factor that is required for mitosis in breast cancer cells. TRPS1 is associated with the NuRD complex, and it regulates cell adhesion, cytoskeleton, and G2-M phase-related genes.

Comparative mechanisms of PAH toxicity by benzo[a]pyrene and dibenzo[def,p]chrysene in primary human bronchial epithelial cells cultured at air-liquid interface

Current assumption for assessing carcinogenic risk of polycyclic aromatic hydrocarbons (PAHs) is that they function through a common mechanism of action; however, recent studies demonstrate that PAHs can act through unique mechanisms potentially contributing to cancer outcomes in a non-additive manner. Using a primary human 3D bronchial epithelial culture (HBEC) model, we assessed potential differences in mechanism of toxicity for two PAHs, benzo[a]pyrene (BAP) and dibenzo[def,p]chrysene (DBC), compared to a complex PAH mixture based on short-term biosignatures identified from transcriptional profiling. Differentiated bronchial epithelial cells were treated with BAP (100-500 μg/ml), DBC (10 μg/ml), and coal tar extract (CTE 500-1500 μg/ml, SRM1597a) for 48 h and gene expression was measured by RNA sequencing or quantitative PCR. Comparison of BAP and DBC gene signatures showed that the majority of genes (~60%) were uniquely regulated by treatment, including signaling pathways for inflammation and DNA damage by DBC and processes for cell cycle, hypoxia and oxidative stress by BAP. Specifically, BAP upregulated targets of AhR, NRF2, and KLF4, while DBC downregulated these same targets, suggesting a chemical-specific pattern in transcriptional regulation involved in antioxidant response, potentially contributing to differences in PAH potency. Other processes were regulated in common by all PAH treatments, BAP, DBC and CTE, including downregulation of genes involved in cell adhesion and reduced functional measurements of barrier integrity. This work supports prior in vivo studies and demonstrates the utility of profiling short-term biosignatures in an organotypic 3D model to identify mechanisms linked to carcinogenic risk of PAHs in humans.

Pathway analysis of GWAS loci identifies novel drug targets and repurposing opportunities

Genome-wide association studies (GWAS) have made considerable progress and there is emerging evidence that genetics-based targets can lead to 28% more launched drugs. We analyzed 1589 GWAS across 1456 pathways to translate these often imprecise genetic loci into therapeutic hypotheses for 182 diseases. These pathway-based genetic targets were validated by testing whether current drug targets were enriched in the pathway space for the same indication. Remarkably, 30% of diseases had significantly more targets in these pathways than expected by chance; the comparable number for GWAS alone (without pathway analysis) was zero. This study shows that a systematic global pathway analysis can translate genetic findings into therapeutic hypotheses for both new drug discovery and repositioning opportunities for current drugs.

Phenotypic and Functional Characteristics of Human Schwann Cells as Revealed by Cell-Based Assays and RNA-SEQ

This study comprehensively addresses the phenotype, function, and whole transcriptome of primary human and rodent Schwann cells (SCs) and highlights key species-specific features beyond the expected donor variability that account for the differential ability of human SCs to proliferate, differentiate, and interact with axons in vitro. Contrary to rat SCs, human SCs were insensitive to mitogenic factors other than neuregulin and presented phenotypic variants at various stages of differentiation, along with a mixture of proliferating and senescent cells, under optimal growth-promoting conditions. The responses of human SCs to cAMP-induced differentiation featured morphological changes and cell cycle exit without a concomitant increase in myelin-related proteins and lipids. Human SCs efficiently extended processes along those of other SCs (human or rat) but failed to do so when placed in co-culture with sensory neurons under conditions supportive of myelination. Indeed, axon contact-dependent human SC alignment, proliferation, and differentiation were not observed and could not be overcome by growth factor supplementation. Strikingly, RNA-seq data revealed that ~ 44 of the transcriptome contained differentially expressed genes in human and rat SCs. A bioinformatics approach further highlighted that representative SC-specific transcripts encoding myelin-related and axon growth-promoting proteins were significantly affected and that a deficient expression of key transducers of cAMP and adhesion signaling explained the fairly limited potential of human SCs to differentiate and respond to axonal cues. These results confirmed the significance of combining traditional bioassays and high-resolution genomics methods to characterize human SCs and identify genes predictive of cell function and therapeutic value.

Histamine Regulates the Inflammatory Profile of SOD1-G93A Microglia and the Histaminergic System Is Dysregulated in Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a late-onset motor neuron disease where activated glia release pro-inflammatory cytokines that trigger a vicious cycle of neurodegeneration in the absence of resolution of inflammation. Given the well-established role of histamine as a neuron-to-glia alarm signal implicated in brain disorders, the aim of this study was to investigate the expression and regulation of the histaminergic pathway in microglial activation in ALS mouse model and in humans. By examining the contribution of the histaminergic system to ALS, we found that particularly via H1 and H4 receptors, histamine promoted an anti-inflammatory profile in microglia from SOD1-G93A mice by modulating their activation state. A decrease in NF-κB and NADPH oxidase 2 with an increase in arginase 1 and P2Y12 receptor was induced by histamine only in the ALS inflammatory environment, but not in the healthy microglia, together with an increase in IL-6, IL-10, CD163, and CD206 phenotypic markers in SOD1-G93A cells. Moreover, histaminergic H1, H2, H3, and H4 receptors, and histamine metabolizing enzymes histidine decarboxylase, histamine N-methyltransferase, and diamine oxidase were found deregulated in spinal cord, cortex, and hypothalamus of SOD1-G93A mice during disease progression. Finally, by performing a meta-analysis study, we found a modulated expression of histamine-related genes in cortex and spinal cord from sporadic ALS patients. Our findings disclose that histamine acts as anti-inflammatory agent in ALS microglia and suggest a dysregulation of the histaminergic signaling in ALS.

Agonist immunotherapy restores T cell function following MEK inhibition improving efficacy in breast cancer

The presence of tumor-infiltrating lymphocytes in triple-negative breast cancers is correlated with improved outcomes. Ras/MAPK pathway activation is associated with significantly lower levels of tumor-infiltrating lymphocytes in triple-negative breast cancers and while MEK inhibition can promote recruitment of tumor-infiltrating lymphocytes to the tumor, here we show that MEK inhibition adversely affects early onset T-cell effector function. We show that α-4-1BB and α-OX-40 T-cell agonist antibodies can rescue the adverse effects of MEK inhibition on T cells in both mouse and human T cells, which results in augmented anti-tumor effects in vivo. This effect is dependent upon increased downstream p38/JNK pathway activation. Taken together, our data suggest that although Ras/MAPK pathway inhibition can increase tumor immunogenicity, the negative impact on T-cell activity is functionally important. This undesirable impact is effectively prevented by combination with T-cell immune agonist immunotherapies resulting in superior therapeutic efficacy.MEK inhibition in breast cancer is associated with increased tumour infiltrating lymphocytes (TILs), however, MAPK activity is required for T cells function. Here the authors show that TILs activity following MEK inhibition can be enhanced by agonist immunotherapy resulting in synergic therapeutic effects.

Selection and Prioritization of Candidate Drug Targets for Amyotrophic Lateral Sclerosis Through a Meta-Analysis Approach

Amyotrophic lateral sclerosis (ALS) is a progressive and incurable neurodegenerative disease. Although several compounds have shown promising results in preclinical studies, their translation into clinical trials has failed. This clinical failure is likely due to the inadequacy of the animal models that do not sufficiently reflect the human disease. Therefore, it is important to optimize drug target selection by identifying those that overlap in human and mouse pathology. We have recently characterized the transcriptional profiles of motor cortex samples from sporadic ALS (SALS) patients and differentiated these into two subgroups based on differentially expressed genes, which encode 70 potential therapeutic targets. To prioritize drug target selection, we investigated their degree of conservation in superoxide dismutase 1 (SOD1) G93A transgenic mice, the most widely used ALS animal model. Interspecies comparison of our human expression data with those of eight different SOD1^G93A datasets present in public repositories revealed the presence of commonly deregulated targets and related biological processes. Moreover, deregulated expression of the majority of our candidate targets occurred at the onset of the disease, offering the possibility to use them for an early and more effective diagnosis and therapy. In addition to highlighting the existence of common key drivers in human and mouse pathology, our study represents the basis for a rational preclinical drug development.

Identification of Differential Gene Expression Patterns after Acute Exposure to High and Low Doses of Low-LET Ionizing Radiation in a Reconstituted Human Skin Tissue

In this study we utilized a systems biology approach to identify dose- (0.1, 2.0 and 10 Gy) and time- (3 and 8 h) dependent responses to acute ionizing radiation exposure in a complex tissue, reconstituted human skin. The low dose used here (0.1 Gy) falls within the range of certain medical diagnostic procedures. Of the two higher doses used, 2.0 Gy is typically administered for radiotherapy, while 10 Gy is lethal. Because exposure to any of these doses is possible after an intentional or accidental radiation events, biomarkers are needed to rapidly and accurately triage potentially exposed individuals. Here, tissue samples were acutely exposed to X-ray-generated low-linear-energy transfer (LET) ionizing radiation, and direct RNA sequencing (RNA-seq) was used to quantify altered transcripts. The time points used for this study aid in assessing early responses to exposure, when key signaling pathways and biomarkers can be identified, which precede and regulate later phenotypic alterations that occur at high doses, including cell death. We determined that a total of 1,701 genes expressed were significantly affected by high-dose radiation, with the majority of genes affected at 10 Gy. Expression levels of a group of 29 genes, including GDF15, BBC3, PPM1D, FDXR, GADD45A, MDM2, CDKN1A, TP53INP1, CYCSP27, SESN1, SESN2, PCNA and AEN, were similarly altered at both 2 and 10 Gy, but not 0.1 Gy, at both time points. A much larger group of upregulated genes, including those involved in inflammatory responses, was significantly altered only after 10 Gy irradiation. At high doses, downregulated genes were associated with cell cycle regulation and exhibited an apparent linear response between 2 and 10 Gy. While only a few genes were significantly affected by 0.1 Gy irradiation, using stringent statistical filters, groups of related genes regulating cell cycle progression and inflammatory responses consistently exhibited opposite trends in their regulation compared to high-dose irradiated groups. Differential regulation of PLK1 signaling at low- and high-dose irradiation was confirmed using qRT-PCR. These results indicate that some alterations in gene expression are qualitatively different at low and high doses of ionizing radiation in this model system. They also highlight potential biomarkers for radiation exposure that may precede the development of overt physiological symptoms in exposed individuals.

Spatial Proximity to Fibroblasts Impacts Molecular Features and Therapeutic Sensitivity of Breast Cancer Cells Influencing Clinical Outcomes

Using a three-dimensional coculture model, we identified significant subtype-specific changes in gene expression, metabolic, and therapeutic sensitivity profiles of breast cancer cells in contact with cancer-associated fibroblasts (CAF). CAF-induced gene expression signatures predicted clinical outcome and immune-related differences in the microenvironment. We found that fibroblasts strongly protect carcinoma cells from lapatinib, attributable to its reduced accumulation in carcinoma cells and an elevated apoptotic threshold. Fibroblasts from normal breast tissues and stromal cultures of brain metastases of breast cancer had similar effects as CAFs. Using synthetic lethality approaches, we identified molecular pathways whose inhibition sensitizes HER2⁺ breast cancer cells to lapatinib both in vitro and in vivo, including JAK2/STAT3 and hyaluronic acid. Neoadjuvant lapatinib therapy in HER2⁺ breast tumors lead to a significant increase of phospho-STAT3⁺ cancer cells and a decrease in the spatial proximity of proliferating (Ki67⁺) cells to CAFs impacting therapeutic responses. Our studies identify CAF-induced physiologically and clinically relevant changes in cancer cells and offer novel approaches for overcoming microenvironment-mediated therapeutic resistance. Cancer Res; 76(22); 6495-506. ©2016 AACR.

Phenotypically anchored transcriptome profiling of developmental exposure to the antimicrobial agent, triclosan, reveals hepatotoxicity in embryonic zebrafish

Triclosan (TCS) is an antimicrobial agent commonly found in a variety of personal care products and cosmetics. TCS readily enters the environment through wastewater and is detected in human plasma, urine, and breast milk due to its widespread use. Studies have implicated TCS as a disruptor of thyroid and estrogen signaling; therefore, research examining the developmental effects of TCS is warranted. In this study, we used embryonic zebrafish to investigate the developmental toxicity and potential mechanism of action of TCS. Embryos were exposed to graded concentrations of TCS from 6 to 120hours post-fertilization (hpf) and the concentration where 80% of the animals had mortality or morbidity at 120hpf (EC80) was calculated. Transcriptomic profiling was conducted on embryos exposed to the EC80 (7.37μM). We identified a total of 922 significant differentially expressed transcripts (FDR adjusted P-value≤0.05; fold change ≥2). Pathway and gene ontology enrichment analyses identified biological networks and transcriptional hubs involving normal liver functioning, suggesting TCS may be hepatotoxic in zebrafish. Tissue-specific gene enrichment analysis further supported the role of the liver as a target organ for TCS toxicity. We also examined the in vitro bioactivity profile of TCS reported by the ToxCast screening program. TCS had a diverse bioactivity profile and was a hit in 217 of the 385 assay endpoints we identified. We observed similarities in gene expression and hepatic steatosis assays; however, hit data for TCS were more concordant with the hypothesized CAR/PXR activity of TCS from rodent and human in vitro studies.

Ligand-Specific Transcriptional Mechanisms Underlie Aryl Hydrocarbon Receptor-Mediated Developmental Toxicity of Oxygenated PAHs

Polycyclic aromatic hydrocarbons (PAHs) are priority environmental contaminants that exhibit mutagenic, carcinogenic, proinflammatory, and teratogenic properties. Oxygen-substituted PAHs (OPAHs) are formed during combustion processes and via phototoxidation and biological degradation of parent (unsubstituted) PAHs. Despite their prevalence both in contaminated industrial sites and in urban air, OPAH mechanisms of action in biological systems are relatively understudied. Like parent PAHs, OPAHs exert structure-dependent mutagenic activities and activation of the aryl hydrocarbon receptor (AHR) and cytochrome p450 metabolic pathway. Four-ring OPAHs 1,9-benz-10-anthrone (BEZO) and benz(a)anthracene-7,12-dione (7,12-B[a]AQ) cause morphological aberrations and induce markers of oxidative stress in developing zebrafish with similar potency, but only 7,12-B[a]AQ induces robust Cyp1a protein expression. We investigated the role of the AHR in mediating the toxicity of BEZO and 7,12-B[a]AQ, and found that knockdown of AHR2 rescued developmental effects caused by both compounds. Using RNA-seq and molecular docking, we identified transcriptional responses that precede developmental toxicity induced via differential interaction with AHR2. Redox-homeostasis genes were affected similarly by these OPAHs, while 7,12-B[a]AQ preferentially activated phase 1 metabolism and BEZO uniquely decreased visual system genes. Analysis of biological functions and upstream regulators suggests that BEZO is a weak AHR agonist, but interacts with other transcriptional regulators to cause developmental toxicity in an AHR-dependent manner. Identifying ligand-dependent AHR interactions and signaling pathways is essential for understanding toxicity of this class of environmentally relevant compounds.

Mechanism-Based Classification of PAH Mixtures to Predict Carcinogenic Potential

We have previously shown that relative potency factors and DNA adduct measurements are inadequate for predicting carcinogenicity of certain polycyclic aromatic hydrocarbons (PAHs) and PAH mixtures, particularly those that function through alternate pathways or exhibit greater promotional activity compared to benzo[a]pyrene (BaP). Therefore, we developed a pathway-based approach for classification of tumor outcome after dermal exposure to PAH/mixtures. FVB/N mice were exposed to dibenzo[def,p]chrysene (DBC), BaP, or environmental PAH mixtures (Mix 1-3) following a 2-stage initiation/promotion skin tumor protocol. Resulting tumor incidence could be categorized by carcinogenic potency as DBC >> BaP = Mix2 = Mix3 > Mix1 = Control, based on statistical significance. Gene expression profiles measured in skin of mice collected 12 h post-initiation were compared with tumor outcome for identification of short-term bioactivity profiles. A Bayesian integration model was utilized to identify biological pathways predictive of PAH carcinogenic potential during initiation. Integration of probability matrices from four enriched pathways (P < .05) for DNA damage, apoptosis, response to chemical stimulus, and interferon gamma signaling resulted in the highest classification accuracy with leave-one-out cross validation. This pathway-driven approach was successfully utilized to distinguish early regulatory events during initiation prognostic for tumor outcome and provides proof-of-concept for using short-term initiation studies to classify carcinogenic potential of environmental PAH mixtures. These data further provide a 'source-to-outcome' model that could be used to predict PAH interactions during tumorigenesis and provide an example of how mode-of-action-based risk assessment could be employed for environmental PAH mixtures.

Development of a Drug-Response Modeling Framework to Identify Cell Line Derived Translational Biomarkers That Can Predict Treatment Outcome to Erlotinib or Sorafenib

Development of drug responsive biomarkers from pre-clinical data is a critical step in drug discovery, as it enables patient stratification in clinical trial design. Such translational biomarkers can be validated in early clinical trial phases and utilized as a patient inclusion parameter in later stage trials. Here we present a study on building accurate and selective drug sensitivity models for Erlotinib or Sorafenib from pre-clinical in vitro data, followed by validation of individual models on corresponding treatment arms from patient data generated in the BATTLE clinical trial. A Partial Least Squares Regression (PLSR) based modeling framework was designed and implemented, using a special splitting strategy and canonical pathways to capture robust information for model building. Erlotinib and Sorafenib predictive models could be used to identify a sub-group of patients that respond better to the corresponding treatment, and these models are specific to the corresponding drugs. The model derived signature genes reflect each drug’s known mechanism of action. Also, the models predict each drug’s potential cancer indications consistent with clinical trial results from a selection of globally normalized GEO expression datasets.

Association between perinatal methylation of the neuronal differentiation regulator HES1 and later childhood neurocognitive function and behaviour

Background Early life environments induce long-term changes in neurocognitive development and behaviour. In animal models, early environmental cues affect neuropsychological phenotypes via epigenetic processes but, as yet, there is little direct evidence for such mechanisms in humans.

Method We examined the relation between DNA methylation at birth and child neuropsychological outcomes in two culturally diverse populations using a genome-wide methylation analysis and validation by pyrosequencing.

Results Within the UK Southampton Women’s Survey (SWS) we first identified 41 differentially methylated regions of interest (DMROI) at birth associated with child’s full-scale IQ at age 4 years. Associations between HES1 DMROI methylation and later cognitive function were confirmed by pyrosequencing in 175 SWS children. Consistent with these findings, higher HES1 methylation was associated with higher executive memory function in a second independent group of 200 SWS 7-year-olds. Finally, we examined a pathway for this relationship within a Singaporean cohort (n = 108). Here, HES1 DMROI methylation predicted differences in early infant behaviour, known to be associated with academic success. In vitro, methylation of HES1 inhibited ETS transcription factor binding, suggesting a functional role of this site.

Conclusions Thus, our findings suggest that perinatal epigenetic processes mark later neurocognitive function and behaviour, providing support for a role of epigenetic processes in mediating the long-term consequences of early life environment on cognitive development.

Reconstitution of the ERG Gene Expression Network Reveals New Biomarkers and Therapeutic Targets in ERG Positive Prostate Tumors

BACKGROUND

Despite a growing number of studies evaluating cancer of prostate (CaP) specific gene alterations, oncogenic activation of the ETS Related Gene (ERG) by gene fusions remains the most validated cancer gene alteration in CaP. Prevalent gene fusions have been described between the ERG gene and promoter upstream sequences of androgen-inducible genes, predominantly TMPRSS2 (transmembrane protease serine 2). Despite the extensive evaluations of ERG genomic rearrangements, fusion transcripts and the ERG oncoprotein, the prognostic value of ERG remains to be better understood. Using gene expression dataset from matched prostate tumor and normal epithelial cells from an 80 GeneChip experiment examining 40 tumors and their matching normal pairs in 40 patients with known ERG status, we conducted a cancer signaling-focused functional analysis of prostatic carcinoma representing moderate and aggressive cancers stratified by ERG expression.

RESULTS

In the present study of matched pairs of laser capture microdissected normal epithelial cells and well-to-moderately differentiated tumor epithelial cells with known ERG gene expression status from 20 patients with localized prostate cancer, we have discovered novel ERG associated biochemical networks.

CONCLUSIONS

Using causal network reconstruction methods, we have identified three major signaling pathways related to MAPK/PI3K cascade that may indeed contribute synergistically to the ERG dependent tumor development. Moreover, the key components of these pathways have potential as biomarkers and therapeutic target for ERG positive prostate tumors.

Data integration reveals key homeostatic mechanisms following low dose radiation exposure

The goal of this study was to define pathways regulated by low dose radiation to understand how biological systems respond to subtle perturbations in their environment and prioritize pathways for human health assessment. Using an in vitro 3-D human full thickness skin model, we have examined the temporal response of dermal and epidermal layers to 10 cGy X-ray using transcriptomic, proteomic, phosphoproteomic and metabolomic platforms. Bioinformatics analysis of each dataset independently revealed potential signaling mechanisms affected by low dose radiation, and integrating data shed additional insight into the mechanisms regulating low dose responses in human tissue. We examined direct interactions among datasets (top down approach) and defined several hubs as significant regulators, including transcription factors (YY1, MYC and CREB1), kinases (CDK2, PLK1) and a protease (MMP2). These data indicate a shift in response across time - with an increase in DNA repair, tissue remodeling and repression of cell proliferation acutely (24-72h). Pathway-based integration (bottom up approach) identified common molecular and pathway responses to low dose radiation, including oxidative stress, nitric oxide signaling and transcriptional regulation through the SP1 factor that would not have been identified by the individual data sets. Significant regulation of key downstream metabolites of nitrative stress was measured within these pathways. Among the features identified in our study, the regulation of MMP2 and SP1 was experimentally validated. Our results demonstrate the advantage of data integration to broadly define the pathways and networks that represent the mechanisms by which complex biological systems respond to perturbation.

Mice lacking the circadian modulators SHARP1 and SHARP2 display altered sleep and mixed state endophenotypes of psychiatric disorders

Increasing evidence suggests that clock genes may be implicated in a spectrum of psychiatric diseases, including sleep and mood related disorders as well as schizophrenia. The bHLH transcription factors SHARP1/DEC2/BHLHE41 and SHARP2/DEC1/BHLHE40 are modulators of the circadian system and SHARP1/DEC2/BHLHE40 has been shown to regulate homeostatic sleep drive in humans. In this study, we characterized Sharp1 and Sharp2 double mutant mice (S1/2-/-) using online EEG recordings in living animals, behavioral assays and global gene expression profiling. EEG recordings revealed attenuated sleep/wake amplitudes and alterations of theta oscillations. Increased sleep in the dark phase is paralleled by reduced voluntary activity and cortical gene expression signatures reveal associations with psychiatric diseases. S1/2-/- mice display alterations in novelty induced activity, anxiety and curiosity. Moreover, mutant mice exhibit impaired working memory and deficits in prepulse inhibition resembling symptoms of psychiatric diseases. Network modeling indicates a connection between neural plasticity and clock genes, particularly for SHARP1 and PER1. Our findings support the hypothesis that abnormal sleep and certain (endo)phenotypes of psychiatric diseases may be caused by common mechanisms involving components of the molecular clock including SHARP1 and SHARP2.

Quantitative iTRAQ LC-MS/MS Proteomics Reveals Transcription Factor Crosstalk and Regulatory Networks in Hypopharyngeal Squamous Cell Carcinoma

To date, no effective therapeutic treatments have been developed for hypopharyngeal squamous cell carcinoma (HPSCC), a disease that has a five-year survival rate of approximately 31% because of its late diagnosis and aggressive nature. Despite recent improvements in diagnostic methods, there are no effective measures to prevent or detect HPSCC in an early stage. The goal of the current study was to identify molecular biomarkers and networks that can facilitate the speedy identification of HPSCC patients who could benefit from individualized treatment. Isobaric tags for relative and absolute quantification (iTRAQ) labeling was employed with two-dimensional liquid chromatography-tandem mass spectrometry to identify quantitatively the differentially expressed proteins among three types of HPSCC disease stages. The iTRAQ results were evaluated by literature searches and western blot analysis. For example, FUBP1, one of 412 proteins with significantly altered expression profiles, was confirmed to have elevated expression in fresh HPSCC tissues. Integrin-mediated cell matrix adhesion and actin filament-inducing cytoskeleton remodeling were the cellular events that were the most relevant to HPSCC tumorigenesis and the metastatic process. The construction of transcriptional regulation networks led to the identification of key transcriptional regulators of tumor development and lymph node metastasis of HPSCC, including Sp1, c-Myc and p53. Additionally, our study indicated that the interactions among Sp1, c-Myc and p53 may play vital roles in the carcinogenesis and metastasis of HPSCC.

Integrative transcriptomic and proteomic analysis of osteocytic cells exposed to fluid flow reveals novel mechano-sensitive signaling pathways

Osteocytes, positioned within bone׳s porous structure, are subject to interstitial fluid flow upon whole bone loading. Such fluid flow is widely theorized to be a mechanical signal transduced by osteocytes, initiating a poorly understood cascade of signaling events mediating bone adaptation to mechanical load. The objective of this study was to examine the time course of flow-induced changes in osteocyte gene transcript and protein levels using high-throughput approaches. Osteocyte-like MLO-Y4 cells were subjected to 2h of oscillating fluid flow (1Pa peak shear stress) and analyzed following 0, 2, 8, and 24h post-flow incubation. Transcriptomic microarray analysis, followed by gene ontology pathway analysis, demonstrated fluid flow regulation of genes consistent with both known and unknown metabolic and inflammatory responses in bone. Additionally, two of the more highly up-regulated gene products - chemokines Cxcl1 and Cxcl2, supported by qPCR - have not previously been reported as responsive to fluid flow. Proteomic analysis demonstrated greatest up-regulation of the ATP-producing enzyme NDK, calcium-binding Calcyclin, and G protein-coupled receptor kinase 6. Finally, an integrative pathway analysis merging fold changes in transcript and protein levels predicted signaling nodes not directly detected at the sampled time points, including transcription factors c-Myc, c-Jun, and RelA/NF-κB. These results extend our knowledge of the osteocytic response to fluid flow, most notably up-regulation of Cxcl1 and Cxcl2 as possible paracrine agents for osteoblastic and osteoclastic recruitment. Moreover, these results demonstrate the utility of integrative, high-throughput approaches in place of a traditional candidate approach for identifying novel mechano-sensitive signaling molecules.

Genetic factors involved in risk for methamphetamine intake and sensitization

Lines of mice were created by selective breeding for the purpose of identifying genetic mechanisms that influence the magnitude of the selected trait and to explore genetic correlations for additional traits thought to be influenced by shared mechanisms. DNA samples from high and low methamphetamine-drinking (MADR) and high and low methamphetamine-sensitization lines were used for quantitative trait locus (QTL) mapping. Significant additive genetic correlations between the two traits indicated a common genetic influence, and a QTL on chromosome X was detected for both traits, suggesting one source of this commonality. For MADR mice, a QTL on chromosome 10 accounted for more than 50 % of the genetic variance in that trait. Microarray gene expression analyses were performed for three brain regions for methamphetamine-naïve MADR line mice: nucleus accumbens, prefrontal cortex, and ventral midbrain. Many of the genes that were differentially expressed between the high and low MADR lines were shared in common across the three brain regions. A gene network highly enriched in transcription factor genes was identified as being relevant to genetically determined differences in methamphetamine intake. When the mu opioid receptor gene (Oprm1), located on chromosome 10 in the QTL region, was added to this top-ranked transcription factor network, it became a hub in the network. These data are consistent with previously published findings of opioid response and intake differences between the MADR lines and suggest that Oprm1, or a gene that impacts activity of the opioid system, plays a role in genetically determined differences in methamphetamine intake.

Molecular analysis of the inhibitory effect of N-acetyl-L-cysteine on the proliferation and invasiveness of pancreatic cancer cells

Preliminary studies have suggested that the reactive oxygen species (ROS) scavenger N-acetyl-L-cysteine (NAC) may be effective in inhibiting the growth of pancreatic cancer cells. In-depth cellular and molecular analyses were carried out to determine NAC's mode of action in inhibiting the growth of a well-characterized pancreatic cancer cell line (AsPC-1). Standardized assays were used to monitor cellular growth, apoptosis, levels of ROS, cellular senescence, migration, and invasiveness. Cell stiffness was measured using atomic force microscopy. Gene expression was monitored by quantitative PCR. NAC significantly inhibits the growth and metastatic potential of AsPC-1 cells by inducing cell-cycle arrest in G1 and subsequent cellular senescence and decreased invasiveness. These anticancer properties are associated with an unexpected increase in the intracellular concentrations of ROS. NAC does not decrease the susceptibility of AsPC-1 cells to the anticancer drugs gemcitabine, mitomycin C, and doxorubicin. NAC-induced changes in gene expression are consistent with the onset of mesenchymal-to-epithelial transition. In conclusion, our findings indicate that NAC induces an integrated series of responses in AsPC-1 cells that make it a highly promising candidate for development as a pancreatic cancer therapeutic.

Molecular profiling of human mammary gland links breast cancer risk to a p27(+) cell population with progenitor characteristics

Early full-term pregnancy is one of the most effective natural protections against breast cancer. To investigate this effect, we have characterized the global gene expression and epigenetic profiles of multiple cell types from normal breast tissue of nulliparous and parous women and carriers of BRCA1 or BRCA2 mutations. We found significant differences in CD44(+) progenitor cells, where the levels of many stem cell-related genes and pathways, including the cell-cycle regulator p27, are lower in parous women without BRCA1/BRCA2 mutations. We also noted a significant reduction in the frequency of CD44(+)p27(+) cells in parous women and showed, using explant cultures, that parity-related signaling pathways play a role in regulating the number of p27(+) cells and their proliferation. Our results suggest that pathways controlling p27(+) mammary epithelial cells and the numbers of these cells relate to breast cancer risk and can be explored for cancer risk assessment and prevention.

Structurally distinct polycyclic aromatic hydrocarbons induce differential transcriptional responses in developing zebrafish

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous in the environment as components of fossil fuels and by-products of combustion. These multi-ring chemicals differentially activate the aryl hydrocarbon receptor (AHR) in a structurally dependent manner, and induce toxicity via both AHR-dependent and -independent mechanisms. PAH exposure is known to induce developmental malformations in zebrafish embryos, and recent studies have shown cardiac toxicity induced by compounds with low AHR affinity. Unraveling the potentially diverse molecular mechanisms of PAH toxicity is essential for understanding the hazard posed by complex PAH mixtures present in the environment. We analyzed transcriptional responses to PAH exposure in zebrafish embryos exposed to benz(a)anthracene (BAA), dibenzothiophene (DBT) and pyrene (PYR) at concentrations that induced developmental malformations by 120 h post-fertilization (hpf). Whole genome microarray analysis of mRNA expression at 24 and 48 hpf identified genes that were differentially regulated over time and in response to the three PAH structures. PAH body burdens were analyzed at both time points using GC-MS, and demonstrated differences in PAH uptake into the embryos. This was important for discerning dose-related differences from those that represented unique molecular mechanisms. While BAA misregulated the least number of transcripts, it caused strong induction of cyp1a and other genes known to be downstream of the AHR, which were not induced by the other two PAHs. Analysis of functional roles of misregulated genes and their predicted regulatory transcription factors also distinguished the BAA response from regulatory networks disrupted by DBT and PYR exposure. These results indicate that systems approaches can be used to classify the toxicity of PAHs based on the networks perturbed following exposure, and may provide a path for unraveling the toxicity of complex PAH mixtures.

Release of severe acute respiratory syndrome coronavirus nuclear import block enhances host transcription in human lung cells

The severe acute respiratory syndrome coronavirus accessory protein ORF6 antagonizes interferon signaling by blocking karyopherin-mediated nuclear import processes. Viral nuclear import antagonists, expressed by several highly pathogenic RNA viruses, likely mediate pleiotropic effects on host gene expression, presumably interfering with transcription factors, cytokines, hormones, and/or signaling cascades that occur in response to infection. By bioinformatic and systems biology approaches, we evaluated the impact of nuclear import antagonism on host expression networks by using human lung epithelial cells infected with either wild-type virus or a mutant that does not express ORF6 protein. Microarray analysis revealed significant changes in differential gene expression, with approximately twice as many upregulated genes in the mutant virus samples by 48 h postinfection, despite identical viral titers. Our data demonstrated that ORF6 protein expression attenuates the activity of numerous karyopherin-dependent host transcription factors (VDR, CREB1, SMAD4, p53, EpasI, and Oct3/4) that are critical for establishing antiviral responses and regulating key host responses during virus infection. Results were confirmed by proteomic and chromatin immunoprecipitation assay analyses and in parallel microarray studies using infected primary human airway epithelial cell cultures. The data strongly support the hypothesis that viral antagonists of nuclear import actively manipulate host responses in specific hierarchical patterns, contributing to the viral pathogenic potential in vivo. Importantly, these studies and modeling approaches not only provide templates for evaluating virus antagonism of nuclear import processes but also can reveal candidate cellular genes and pathways that may significantly influence disease outcomes following severe acute respiratory syndrome coronavirus infection in vivo.

Proteomic profile identifies dysregulated pathways in Cornelia de Lange syndrome cells with distinct mutations in SMC1A and SMC3 genes

Mutations in cohesin genes have been identified in Cornelia de Lange syndrome (CdLS), but its etiopathogenetic mechanisms are still poorly understood. To define biochemical pathways that are affected in CdLS, we analyzed the proteomic profile of CdLS cell lines carrying mutations in the core cohesin genes, SMC1A and SMC3. Dysregulated protein expression was found in CdLS probands compared to controls. The proteomics analysis was able to discriminate between probands harboring mutations in the different domains of the SMC proteins. In particular, proteins involved in the response to oxidative stress were specifically down-regulated in hinge mutated probands. In addition, the finding that CdLS cell lines show an increase in global oxidative stress argues that it could contribute to some CdLS phenotypic features such as premature physiological aging and genome instability. Finally, the c-MYC gene represents a convergent hub lying at the center of dysregulated pathways, and is down-regulated in CdLS. This study allowed us to highlight, for the first time, specific biochemical pathways that are affected in CdLS, providing plausible causal evidence for some of the phenotypic features seen in CdLS.

Cross-species genome wide expression analysis during pluripotent cell determination in mouse and rat preimplantation embryos

The transition between morula and blastocyst stage during preimplantation development represents the first differentiation event of embryogenesis. Morula cells undergo the first cellular specialization and produce two well-defined populations of cells, the trophoblast and the inner cell mass (ICM). Embryonic stem cells (ESCs) with unlimited self-renewal capacity are believed to represent the in vitro counterpart of the ICM. Both mouse and rat ESCs can be derived from the ICM cells, but their in vitro stability differs. In this study we performed a microarray analysis in which we compared the transcriptome of mouse and rat morula, blastocyst, and ICM. This cross-species comparison represents a good model for understanding the differences in derivation and cultivation of ESCs observed in the two species. In order to identify alternative regulation of important molecular mechanisms the investigation of differential gene expression between the two species was extended at the level of signaling pathways, gene families, and single selected genes of interest. Some of the genes differentially expressed between the two species are already known to be important factors in the maintenance of pluripotency in ESCs, like for example Sox2 or Stat3, or play a role in reprogramming somatic cells to pluripotency like c-Myc, Klf4 and p53 and therefore represent interesting candidates to further analyze in vitro in the rat ESCs. This is the first study investigating the gene expression changes during the transition from morula to blastocyst in the rat preimplantation development. Our data show that in the pluripotent pool of cells of the rat and mouse preimplantation embryo substantial differential regulation of genes is present, which might explain the difficulties observed for the derivation and culture of rat ESCs using mouse conditions.

Ouabain Mimics Low Temperature Rescue of F508del-CFTR in Cystic Fibrosis Epithelial Cells

Most cases of cystic fibrosis (CF) are caused by the deletion of a single phenylalanine residue at position 508 of the cystic fibrosis transmembrane conductance regulator (CFTR). The mutant F508del-CFTR is retained in the endoplasmic reticulum and degraded, but can be induced by low temperature incubation (29°C) to traffic to the plasma membrane where it functions as a chloride channel. Here we show that, cardiac glycosides, at nanomolar concentrations, can partially correct the trafficking of F508del-CFTR in human CF bronchial epithelial cells (CFBE41o-) and in an F508del-CFTR mouse model. Comparison of the transcriptional profiles obtained with polarized CFBE41o-cells after treatment with ouabain and by low temperature has revealed a striking similarity between the two corrector treatments that is not shared with other correctors. In summary, our study shows a novel function of ouabain and its analogs in the regulation of F508del-CFTR trafficking and suggests that compounds that mimic this low temperature correction of trafficking will provide new avenues for the development of therapeutics for CF.

Transcriptome changes affecting Hedgehog and cytokine signalling in the umbilical cord: implications for disease risk

BACKGROUND

Babies born at lower gestational ages or smaller birthweights have a greater risk of poorer health in later life. Both the causes of these sub-optimal birth outcomes and the mechanism by which the effects are transmitted over decades are the subject of extensive study. We investigated whether a transcriptomic signature of either birthweight or gestational age could be detected in umbilical cord RNA.

METHODS

The gene expression patterns of 32 umbilical cords from Singaporean babies of Chinese ethnicity across a range of birthweights (1698-4151 g) and gestational ages (35-41 weeks) were determined. We confirmed the differential expression pattern by gestational age for 12 genes in a series of 127 umbilical cords of Chinese, Malay and Indian ethnicity.

RESULTS

We found that the transcriptome is substantially influenced by gestational age; but less so by birthweight. We show that some of the expression changes dependent on gestational age are enriched in signal transduction pathways, such as Hedgehog and in genes with roles in cytokine signalling and angiogenesis. We show that some of the gene expression changes we report are reflected in the epigenome.

CONCLUSIONS

We studied the umbilical cord which is peripheral to disease susceptible tissues. The results suggest that soma-wide transcriptome changes, preserved at the epigenetic level, may be a mechanism whereby birth outcomes are linked to the risk of adult metabolic and arthritic disease and suggest that greater attention be given to the association between premature birth and later disease risk.

Genomic characterization of explant tumorgraft models derived from fresh patient tumor tissue

Background

There is resurgence within drug and biomarker development communities for the use of primary tumorgraft models as improved predictors of patient tumor response to novel therapeutic strategies. Despite perceived advantages over cell line derived xenograft models, there is limited data comparing the genotype and phenotype of tumorgrafts to the donor patient tumor, limiting the determination of molecular relevance of the tumorgraft model. This report directly compares the genomic characteristics of patient tumors and the derived tumorgraft models, including gene expression, and oncogenic mutation status.

Methods

Fresh tumor tissues from 182 cancer patients were implanted subcutaneously into immune-compromised mice for the development of primary patient tumorgraft models. Histological assessment was performed on both patient tumors and the resulting tumorgraft models. Somatic mutations in key oncogenes and gene expression levels of resulting tumorgrafts were compared to the matched patient tumors using the OncoCarta (Sequenom, San Diego, CA) and human gene microarray (Affymetrix, Santa Clara, CA) platforms respectively. The genomic stability of the established tumorgrafts was assessed across serial in vivo generations in a representative subset of models. The genomes of patient tumors that formed tumorgrafts were compared to those that did not to identify the possible molecular basis to successful engraftment or rejection.

Results

Fresh tumor tissues from 182 cancer patients were implanted into immune-compromised mice with forty-nine tumorgraft models that have been successfully established, exhibiting strong histological and genomic fidelity to the originating patient tumors. Comparison of the transcriptomes and oncogenic mutations between the tumorgrafts and the matched patient tumors were found to be stable across four tumorgraft generations. Not only did the various tumors retain the differentiation pattern, but supporting stromal elements were preserved. Those genes down-regulated specifically in tumorgrafts were enriched in biological pathways involved in host immune response, consistent with the immune deficiency status of the host. Patient tumors that successfully formed tumorgrafts were enriched for cell signaling, cell cycle, and cytoskeleton pathways and exhibited evidence of reduced immunogenicity.

Conclusions

The preservation of the patient’s tumor genomic profile and tumor microenvironment supports the view that primary patient tumorgrafts provide a relevant model to support the translation of new therapeutic strategies and personalized medicine approaches in oncology.

Identification of common biological pathways and drug targets across multiple respiratory viruses based on human host gene expression analysis

BACKGROUND

Pandemic and seasonal respiratory viruses are a major global health concern. Given the genetic diversity of respiratory viruses and the emergence of drug resistant strains, the targeted disruption of human host-virus interactions is a potential therapeutic strategy for treating multi-viral infections. The availability of large-scale genomic datasets focused on host-pathogen interactions can be used to discover novel drug targets as well as potential opportunities for drug repositioning.

METHODS/RESULTS

In this study, we performed a large-scale analysis of microarray datasets involving host response to infections by influenza A virus, respiratory syncytial virus, rhinovirus, SARS-coronavirus, metapneumonia virus, coxsackievirus and cytomegalovirus. Common genes and pathways were found through a rigorous, iterative analysis pipeline where relevant host mRNA expression datasets were identified, analyzed for quality and gene differential expression, then mapped to pathways for enrichment analysis. Possible repurposed drugs targets were found through database and literature searches. A total of 67 common biological pathways were identified among the seven different respiratory viruses analyzed, representing fifteen laboratories, nine different cell types, and seven different array platforms. A large overlap in the general immune response was observed among the top twenty of these 67 pathways, adding validation to our analysis strategy. Of the top five pathways, we found 53 differentially expressed genes affected by at least five of the seven viruses. We suggest five new therapeutic indications for existing small molecules or biological agents targeting proteins encoded by the genes F3, IL1B, TNF, CASP1 and MMP9. Pathway enrichment analysis also identified a potential novel host response, the Parkin-Ubiquitin Proteasomal System (Parkin-UPS) pathway, which is known to be involved in the progression of neurodegenerative Parkinson's disease.

CONCLUSIONS

Our study suggests that multiple and diverse respiratory viruses invoke several common host response pathways. Further analysis of these pathways suggests potential opportunities for therapeutic intervention.

Phosphoproteomic analysis: an emerging role in deciphering cellular signaling in human embryonic stem cells and their differentiated derivatives

Cellular signaling is largely controlled by protein phosphorylation. This post-translational modification (PTM) has been extensively analyzed when examining one or a few protein phosphorylation events that effect cell signaling. However, protein kinase-driven signaling networks, comprising total (phospho)proteomes, largely control cell fate. Therefore, large-scale analysis of differentially regulated protein phosphorylation is central to elucidating complex cellular events, including maintenance of pluripotency and differentiation of embryonic stem cells (ESCs). The current technology of choice for total phosphoproteome and combined total proteome plus total phosphoproteome (termed (phospho)proteome) analyses is multidimensional liquid chromatography-(MDLC) tandem mass spectrometry (MS/MS). Advances in the use of MDLC for separation of peptides comprising total (phospho)proteomes, phosphopeptide enrichment, separation of enriched fractions, and quantitative peptide identification by MS/MS have been rapid in recent years, as have improvements in the sensitivity, speed, and accuracy of mass spectrometers. Increasingly deep coverage of (phospho)proteomes is allowing an improved understanding of changes in protein phosphorylation networks as cells respond to stimuli and progress from one undifferentiated or differentiated state to another. Although MDLC-MS/MS studies are powerful, understanding the interpretation of the data is important, and targeted experimental pursuit of biological predictions provided by total (phospho)proteome analyses is needed. (Phospho)proteomic analyses of pluripotent stem cells are in their infancy at this time. However, such studies have already begun to contribute to an improved and accelerated understanding of basic pluripotent stem cell signaling and fate control, especially at the systems-biology level.

Deregulated hepatic metabolism exacerbates impaired testosterone production in Mrp4-deficient mice

The physiological role of multidrug resistance protein 4 (Mrp4, Abcc4) in the testes is unknown. We found that Mrp4 is expressed primarily in mouse and human Leydig cells; however, there is no current evidence that Mrp4 regulates testosterone production. We investigated its role in Leydig cells, where testosterone production is regulated by cAMP, an intracellular messenger formed when the luteinizing hormone (LH) receptor is activated. Because Mrp4 regulates cAMP, we compared testosterone levels in Mrp4(-/-) and Mrp4(+/+) mice. Young Mrp4(-/-) mice had significantly impaired gametogenesis, reduced testicular testosterone, and disruption of Leydig cell cAMP homeostasis. Both young and adult mice had impaired testosterone production. In Mrp4(-/-) primary Leydig cells treated with LH, intracellular cAMP production was impaired and cAMP-response element-binding protein (CREB) phosphorylation was strongly attenuated. Notably, expression of CREB target genes that regulate testosterone biosynthesis was reduced in Mrp4(-/-) Leydig cells in vivo. Therefore, Mrp4 is required for normal Leydig cell testosterone production. However, adult Mrp4(-/-) mice are fertile, with a normal circulating testosterone concentration. The difference is that in 3-week-old Mrp4(-/-) mice, disruption of gonadal testosterone production up-regulates hepatic Cyp2b10, a known testosterone-metabolizing enzyme. Therefore, defective testicular testosterone production de-regulates hepatic Cyp-mediated testosterone metabolism to disrupt gametogenesis. These findings have important implications for understanding the side effects of therapeutics that disrupt Mrp4 function and are reported to alter androgen production.

Putative multifunctional signature of lung metastases in dedifferentiated chondrosarcoma

Chondrosarcomas are among the most malignant skeletal tumors. Dedifferentiated chondrosarcoma is a highly aggressive subtype of chondrosarcoma, with lung metastases developing within a few months of diagnosis in 90% of patients. In this paper we performed comparative analyses of the transcriptomes of five individual metastatic lung lesions that were surgically resected from a patient with dedifferentiated chondrosarcoma. We document for the first time a high heterogeneity of gene expression profiles among the individual lung metastases. Moreover, we reveal a signature of “multifunctional” genes that are expressed in all metastatic lung lesions. Also, for the first time, we document the occurrence of massive macrophage infiltration in dedifferentiated chondrosarcoma lung metastases.

Analysis of transcriptional factors and regulation networks in laryngeal squamous cell carcinoma patients with lymph node metastasis

The present study was to identify and quantitate differentially expressed proteins in laryngeal squamous cell carcinoma (LSCC) tissues with or without lymph node metastasis and to explore transcriptional factors and regulation networks associated with the process. Tissue specimens were taken from 20 patients with LSCC, including 10 cases of LSCC without metastasis LSCC (N0) and 10 cases of LSCC with metastasis LSCC (Nx). Among the 643 unique proteins identified by using iTRAQ labeling and quantitative proteomic technology, 389 proteins showed an abundance change in LSCC (Nx) as compared to LSCC (N0). Cytoskeleton remodeling, cell adhesion, and immune response activation were found to be the main processes in LSCC metastasis. The construction of transcription regulation networks identified key transcription regulators for lymph node metastasis of LSCC, including Sp1, c-myc, and p53, which may affect LSCC metastasis through the epithelial-mesenchymal transition. Furthermore, our results suggest that ubiquitination may be a critical factor in the networks. The present study provides insights into transcriptional factors and regulation networks involved in LSCC metastasis, which may lead to new strategies for treatment of LSCC metastasis.

Optimal deconvolution of transcriptional profiling data using quadratic programming with application to complex clinical blood samples

Large-scale molecular profiling technologies have assisted the identification of disease biomarkers and facilitated the basic understanding of cellular processes. However, samples collected from human subjects in clinical trials possess a level of complexity, arising from multiple cell types, that can obfuscate the analysis of data derived from them. Failure to identify, quantify, and incorporate sources of heterogeneity into an analysis can have widespread and detrimental effects on subsequent statistical studies.We describe an approach that builds upon a linear latent variable model, in which expression levels from mixed cell populations are modeled as the weighted average of expression from different cell types. We solve these equations using quadratic programming, which efficiently identifies the globally optimal solution while preserving non-negativity of the fraction of the cells. We applied our method to various existing platforms to estimate proportions of different pure cell or tissue types and gene expression profilings of distinct phenotypes, with a focus on complex samples collected in clinical trials. We tested our methods on several well controlled benchmark data sets with known mixing fractions of pure cell or tissue types and mRNA expression profiling data from samples collected in a clinical trial. Accurate agreement between predicted and actual mixing fractions was observed. In addition, our method was able to predict mixing fractions for more than ten species of circulating cells and to provide accurate estimates for relatively rare cell types (<10% total population). Furthermore, accurate changes in leukocyte trafficking associated with Fingolomid (FTY720) treatment were identified that were consistent with previous results generated by both cell counts and flow cytometry. These data suggest that our method can solve one of the open questions regarding the analysis of complex transcriptional data: namely, how to identify the optimal mixing fractions in a given experiment.

A model of early human embryonic stem cell differentiation reveals inter- and intracellular changes on transition to squamous epithelium

The molecular events leading to human embryonic stem cell (hESC) differentiation are the subject of considerable scrutiny. Here, we characterize an in vitro model that permits analysis of the earliest steps in the transition of hESC colonies to squamous epithelium on basic fibroblast growth factor withdrawal. A set of markers (GSC, CK18, Gata4, Eomes, and Sox17) point to a mesendodermal nature of the epithelial cells with subsequent commitment to definitive endoderm (Sox17, Cdx2, nestin, and Islet1). We assayed alterations in the transcriptome in parallel with the distribution of immunohistochemical markers. Our results indicate that the alterations of tight junctions in pluripotent culture precede the beginning of differentiation. We defined this cell population as "specified," as it is committed toward differentiation. The transitional zone between "specified" pluripotent and differentiated cells displays significant up-regulation of keratin-18 (CK18) along with a decrease in the functional activity of gap junctions and the down-regulation of 2 gap junction proteins, connexin 43 (Cx43) and connexin 45 (Cx45), which is coincidental with substantial elevation of intracellular Ca2+ levels. These findings reveal a set of cellular changes that may represent the earliest markers of in vitro hESC transition to an epithelial phenotype, before the induction of gene expression networks that guide hESC differentiation. Moreover, we hypothesize that these events may be common during the primary steps of hESC commitment to functionally varied epithelial tissue derivatives of different embryological origins.

The JAK2/STAT3 signaling pathway is required for growth of CD44⁺CD24⁻ stem cell-like breast cancer cells in human tumors

Intratumor heterogeneity is a major clinical problem because tumor cell subtypes display variable sensitivity to therapeutics and may play different roles in progression. We previously characterized 2 cell populations in human breast tumors with distinct properties: CD44+CD24- cells that have stem cell-like characteristics, and CD44-CD24+ cells that resemble more differentiated breast cancer cells. Here we identified 15 genes required for cell growth or proliferation in CD44+CD24- human breast cancer cells in a large-scale loss-of-function screen and found that inhibition of several of these (IL6, PTGIS, HAS1, CXCL3, and PFKFB3) reduced Stat3 activation. We found that the IL-6/JAK2/Stat3 pathway was preferentially active in CD44+CD24- breast cancer cells compared with other tumor cell types, and inhibition of JAK2 decreased their number and blocked growth of xenografts. Our results highlight the differences between distinct breast cancer cell types and identify targets such as JAK2 and Stat3 that may lead to more specific and effective breast cancer therapies.

A systems biology strategy for predicting similarities and differences of drug effects: evidence for drug-specific modulation of inflammation in atherosclerosis

Background

Successful drug development has been hampered by a limited understanding of how to translate laboratory-based biological discoveries into safe and effective medicines. We have developed a generic method for predicting the effects of drugs on biological processes. Information derived from the chemical structure and experimental omics data from short-term efficacy studies are combined to predict the possible protein targets and cellular pathways affected by drugs.

Results

Validation of the method with anti-atherosclerotic compounds (fenofibrate, rosuvastatin, LXR activator T0901317) demonstrated a great conformity between the computationally predicted effects and the wet-lab biochemical effects. Comparative genome-wide pathway mapping revealed that the biological drug effects were realized largely via different pathways and mechanisms. In line with the predictions, the drugs showed differential effects on inflammatory pathways (downstream of PDGF, VEGF, IFNγ, TGFβ, IL1β, TNFα, LPS), transcriptional regulators (NFκB, C/EBP, STAT3, AP-1) and enzymes (PKCδ, AKT, PLA2), and they quenched different aspects of the inflammatory signaling cascade. Fenofibrate, the compound predicted to be most efficacious in inhibiting early processes of atherosclerosis, had the strongest effect on early lesion development.

Conclusion

Our approach provides mechanistic rationales for the differential and common effects of drugs and may help to better understand the origins of drug actions and the design of combination therapies.

Epigenetic regulation of cell type-specific expression patterns in the human mammary epithelium

Differentiation is an epigenetic program that involves the gradual loss of pluripotency and acquisition of cell type–specific features. Understanding these processes requires genome-wide analysis of epigenetic and gene expression profiles, which have been challenging in primary tissue samples due to limited numbers of cells available. Here we describe the application of high-throughput sequencing technology for profiling histone and DNA methylation, as well as gene expression patterns of normal human mammary progenitor-enriched and luminal lineage-committed cells. We observed significant differences in histone H3 lysine 27 tri-methylation (H3K27me3) enrichment and DNA methylation of genes expressed in a cell type–specific manner, suggesting their regulation by epigenetic mechanisms and a dynamic interplay between the two processes that together define developmental potential. The technologies we developed and the epigenetically regulated genes we identified will accelerate the characterization of primary cell epigenomes and the dissection of human mammary epithelial lineage-commitment and luminal differentiation.

Cellular differentiation is a precisely controlled and largely irreversible process orchestrated by cell type–specific epigenetic programs. Abnormalities in these programs lead to developmental disorders and play a key role in tumorigenesis. To better understand the regulation of human mammary epithelial cell type specification, we analyzed the gene expression, DNA methylation, and histone H3 K4 and K27 trimethylation profiles of progenitor-enriched and more differentiated luminal epithelial cell populations from multiple individuals. Network analysis of these profiles and their comparison to that of human embryonic stem cells identified key regulators of mammary epithelial and luminal lineage commitment. The list of genes epigenetically regulated in a cell type–specific manner provides a rich resource for the further analysis of human breast development and the role of epigenetic mechanisms in breast tumorigenesis.

The role of computational methods in the identification of bioactive compounds

Computational methods play an ever increasing role in lead finding. A vast repertoire of molecular design and virtual screening methods emerged in the past two decades and are today routinely used. There is increasing awareness that there is no single best computational protocol and correspondingly there is a shift recommending the combination of complementary methods. A promising trend for the application of computational methods in lead finding is to take advantage of the vast amounts of HTS (High Throughput Screening) data to allow lead assessment by detailed systems-based data analysis, especially for phenotypic screens where the identification of compound-target pairs is the primary goal. Herein, we review trends and provide examples of successful applications of computational methods in lead finding.