The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins

modEnrichr: a suite of gene set enrichment analysis tools for model organisms

High-throughput experiments produce increasingly large datasets that are difficult to analyze and integrate. While most data integration approaches focus on aligning metadata, data integration can be achieved by abstracting experimental results into gene sets. Such gene sets can be made available for reuse through gene set enrichment analysis tools such as Enrichr. Enrichr currently only supports gene sets compiled from human and mouse, limiting accessibility for investigators that study other model organisms. modEnrichr is an expansion of Enrichr for four model organisms: fish, fly, worm and yeast. The gene set libraries within FishEnrichr, FlyEnrichr, WormEnrichr and YeastEnrichr are created from the Gene Ontology, mRNA expression profiles, GeneRIF, pathway databases, protein domain databases and other organism-specific resources. Additionally, libraries were created by predicting gene function from RNA-seq co-expression data processed uniformly from the gene expression omnibus for each organism. The modEnrichr suite of tools provides the ability to convert gene lists across species using an ortholog conversion tool that automatically detects the species. For complex analyses, modEnrichr provides API access that enables submitting batch queries. In summary, modEnrichr leverages existing model organism databases and other resources to facilitate comprehensive hypothesis generation through data integration.

GraphGR: A graph neural network to predict the effect of pharmacotherapy on the cancer cell growth.. [DOI: 10.1101/2020.05.20.107458]

Genomic profiles of cancer cells provide valuable information on genetic alterations in cancer. Several recent studies employed these data to predict the response of cancer cell lines to treatment with drugs. Nonetheless, due to the multifactorial phenotypes and intricate mechanisms of cancer, the accurate prediction of the effect of pharmacotherapy on a specific cell line based on the genetic information alone is problematic. High prediction accuracies reported in the literature likely result from significant overlaps among training, validation, and testing sets, making many predictors inapplicable to new data. To address these issues, we developed GraphGR, a graph neural network with sophisticated attention propagation mechanisms to predict the therapeutic effects of kinase inhibitors across various tumors. Emphasizing on the system-level complexity of cancer, GraphGR integrates multiple heterogeneous data, such as biological networks, genomics, inhibitor profiling, and genedisease associations, into a unified graph structure. In order to construct diverse and information-rich cancer-specific networks, we devised a novel graph reduction protocol based on not only the topological information, but also the biological knowledge. The performance of GraphGR, properly cross-validated at the tissue level, is 0.83 in terms of the area under the receiver operating characteristics, which is notably higher than those measured for other approaches on the same data. Finally, several new predictions are validated against the biomedical literature demonstrating that GraphGR generalizes well to unseen data, i.e. it can predict therapeutic effects across a variety of cancer cell lines and inhibitors. GraphGR is freely available to the academic community at https://github.com/pulimeng/GraphGR.

Ethanol Extract of the Infructescence of Platycarya strobilacea Sieb. et Zucc. Induces Methuosis of Human Nasopharyngeal Carcinoma Cells

The infructescence of Platycarya strobilacea Sieb. et Zucc. (PS) has been used in the treatment of rhinitis and sinusitis in clinical practice. Our preliminary study showed that an ethanol extract of the infructescence of PS (EPS) had significant antinasopharyngeal carcinoma (NPC) effects in vitro. However, the mechanism underlying the NPS cell death induced by EPS remains unclear. The aim of the present study was to investigate the inhibitory effects of EPS on NPC cells and to elucidate the underlying mechanism. The effects of EPS on NPC cells were investigated in CNE1 and CNE2 cells in vitro. In EPS-treated cells, the cell morphological changes were evaluated through fluorescence microscope, transmission electron microscopy, and flow cytometry. The underlying mechanism was analyzed via network pharmacology and further verified by western blot analysis. The anticancer effects of EPS were associated with the generation of CNE1 and CNE2 cell fusion and vacuoles, the perturbation of lysosomal vesicle transportation, and the induction of methuosis. The network pharmacology and western blot results indicated that the effect of EPS in NPC cells might be achieved via regulation of the Ras proto-oncogene (RAS)/mitogen-activated protein kinase (MAPK) signaling pathway and the transcription factor c-Fos proto-oncogene (c-FOS) and its downstream genes. EPS induces NPC cell death through methuosis. The mechanism might be related to regulation of the transcription factor c-FOS and its downstream genes.

Geneshot: search engine for ranking genes from arbitrary text queries

The frequency by which genes are studied correlates with the prior knowledge accumulated about them. This leads to an imbalance in research attention where some genes are highly investigated while others are ignored. Geneshot is a search engine developed to illuminate this gap and to promote attention to the under-studied genome. Through a simple web interface, Geneshot enables researchers to enter arbitrary search terms, to receive ranked lists of genes relevant to the search terms. Returned ranked gene lists contain genes that were previously published in association with the search terms, as well as genes predicted to be associated with the terms based on data integration from multiple sources. The search results are presented with interactive visualizations. To predict gene function, Geneshot utilizes gene–gene similarity matrices from processed RNA-seq data, or from gene–gene co-occurrence data obtained from multiple sources. In addition, Geneshot can be used to analyze the novelty of gene sets and augment gene sets with additional relevant genes. The Geneshot web-server and API are freely and openly available from https://amp.pharm.mssm.edu/geneshot.

BAHD1 haploinsufficiency results in anxiety-like phenotypes in male mice

BAHD1 is a heterochomatinization factor recently described as a component of a multiprotein complex associated with histone deacetylases HDAC1/2. The physiological and patho-physiological functions of BAHD1 are not yet well characterized. Here, we examined the consequences of BAHD1 deficiency in the brains of male mice. While Bahd1 knockout mice had no detectable defects in brain anatomy, RNA sequencing profiling revealed about 2500 deregulated genes in Bahd1-/- brains compared to Bahd1+/+ brains. A majority of these genes were involved in nervous system development and function, behavior, metabolism and immunity. Exploration of the Allen Brain Atlas and Dropviz databases, assessing gene expression in the brain, revealed that expression of the Bahd1 gene was limited to a few territories and cell subtypes, particularly in the hippocampal formation, the isocortex and the olfactory regions. The effect of partial BAHD1 deficiency on behavior was then evaluated on Bahd1 heterozygous male mice, which have no lethal or metabolic phenotypes. Bahd1+/- mice showed anxiety-like behavior and reduced prepulse inhibition (PPI) of the startle response. Altogether, these results suggest that BAHD1 plays a role in chromatin-dependent gene regulation in a subset of brain cells and support recent evidence linking genetic alteration of BAHD1 to psychiatric disorders in a human patient.

The COVID-19 Gene and Drug Set Library

The coronavirus (CoV) severe acute respiratory syndrome (SARS)-CoV-2 (COVID-19) pandemic has received rapid response by the research community to offer suggestions for repurposing of approved drugs as well as to improve our understanding of the COVID-19 viral life cycle molecular mechanisms. In a short period, tens of thousands of research preprints and other publications have emerged including those that report lists of experimentally validated drugs and compounds as potential COVID-19 therapies. In addition, gene sets from interacting COVID-19 virus-host proteins and differentially expressed genes when comparing infected to uninfected cells are being published at a fast rate. To organize this rapidly accumulating knowledge, we developed the COVID-19 Gene and Drug Set Library (https://amp.pharm.mssm.edu/covid19/), a collection of gene and drug sets related to COVID-19 research from multiple sources. The COVID-19 Gene and Drug Set Library is delivered as a web-based interface that enables users to view, download, analyze, visualize, and contribute gene and drug sets related to COVID-19 research. To evaluate the content of the library, we performed several analyses including comparing the results from 6 in-vitro drug screens for COVID-19 repurposing candidates. Surprisingly, we observe little overlap across these initial screens. The most common and unique hit across these screen is mefloquine, a malaria drug that should receive more attention as a potential therapeutic for COVID-19. Overall, the library of gene and drug sets can be used to identify community consensus, make researchers and clinicians aware of the development of new potential therapies, as well as allow the research community to work together towards a cure for COVID-19.

Naked mole-rat very-high-molecular-mass hyaluronan exhibits superior cytoprotective properties

Naked mole-rat (NMR), the longest-living rodent, produces very-high-molecular-mass hyaluronan (vHMM-HA), compared to other mammalian species. However, it is unclear if exceptional polymer length of vHMM-HA is important for longevity. Here, we show that vHMM-HA (>6.1 MDa) has superior cytoprotective properties compared to the shorter HMM-HA. It protects not only NMR cells, but also mouse and human cells from stress-induced cell-cycle arrest and cell death in a polymer length-dependent manner. The cytoprotective effect is dependent on the major HA-receptor, CD44. We find that vHMM-HA suppresses CD44 protein-protein interactions, whereas HMM-HA promotes them. As a result, vHMM-HA and HMM-HA induce opposing effects on the expression of CD44-dependent genes, which are associated with the p53 pathway. Concomitantly, vHMM-HA partially attenuates p53 and protects cells from stress in a p53-dependent manner. Our results implicate vHMM-HA in anti-aging mechanisms and suggest the potential applications of vHMM-HA for enhancing cellular stress resistance.

A Deregulated HOX Gene Axis Confers an Epigenetic Vulnerability in KRAS-Mutant Lung Cancers

While KRAS mutations are common in non-small cell lung cancer (NSCLC), effective treatments are lacking. Here, we report that half of KRAS-mutant NSCLCs aberrantly express the homeobox protein HOXC10, largely due to unappreciated defects in PRC2, which confers sensitivity to combined BET/MEK inhibitors in xenograft and PDX models. Efficacy of the combination is dependent on suppression of HOXC10 by BET inhibitors. We further show that HOXC10 regulates the expression of pre-replication complex (pre-RC) proteins in sensitive tumors. Accordingly, BET/MEK inhibitors suppress pre-RC proteins in cycling cells, triggering stalled replication, DNA damage, and death. These studies reveal a promising therapeutic strategy for KRAS-mutant NSCLCs, identify a predictive biomarker of response, and define a subset of NSCLCs with a targetable epigenetic vulnerability.

Linker histone H1.2 and H1.4 affect the neutrophil lineage determination

Neutrophils are important innate immune cells that tackle invading pathogens with different effector mechanisms. They acquire this antimicrobial potential during their maturation in the bone marrow, where they differentiate from hematopoietic stem cells in a process called granulopoiesis. Mature neutrophils are terminally differentiated and short-lived with a high turnover rate. Here, we show a critical role for linker histone H1 on the differentiation and function of neutrophils using a genome-wide CRISPR/Cas9 screen in the human cell line PLB-985. We systematically disrupted expression of somatic H1 subtypes to show that individual H1 subtypes affect PLB-985 maturation in opposite ways. Loss of H1.2 and H1.4 induced an eosinophil-like transcriptional program, thereby negatively regulating the differentiation into the neutrophil lineage. Importantly, H1 subtypes also affect neutrophil differentiation and the eosinophil-directed bias of murine bone marrow stem cells, demonstrating an unexpected subtype-specific role for H1 in granulopoiesis.

Epigenome environment interactions accelerate epigenomic aging and unlock metabolically restricted epigenetic reprogramming in adulthood

Our early-life environment has a profound influence on developing organs that impacts metabolic function and determines disease susceptibility across the life-course. Using a rat model for exposure to an endocrine disrupting chemical (EDC), we show that early-life chemical exposure causes metabolic dysfunction in adulthood and reprograms histone marks in the developing liver to accelerate acquisition of an adult epigenomic signature. This epigenomic reprogramming persists long after the initial exposure, but many reprogrammed genes remain transcriptionally silent with their impact on metabolism not revealed until a later life exposure to a Western-style diet. Diet-dependent metabolic disruption was largely driven by reprogramming of the Early Growth Response 1 (EGR1) transcriptome and production of metabolites in pathways linked to cholesterol, lipid and one-carbon metabolism. These findings demonstrate the importance of epigenome:environment interactions, which early in life accelerate epigenomic aging, and later in adulthood unlock metabolically restricted epigenetic reprogramming to drive metabolic dysfunction.

Early life exposure to environmental stressors, including endocrine disrupting chemicals (EDCs), can impact health later in life. Here, the authors show that neonatal EDC exposure in rats causes epigenetic reprogramming in the liver, which is transcriptionally silent until animals are placed on a Western-style diet.

Helicobacter pylori infection impairs chaperone-assisted maturation of Na-K-ATPase in gastric epithelium

Helicobacter pylori infection always induces gastritis, which may progress to ulcer disease or cancer. The mechanisms underlying mucosal injury by the bacteria are incompletely understood. Here, we identify a novel pathway for H. pylori-induced gastric injury, the impairment of maturation of the essential transport enzyme and cell adhesion molecule, Na-K-ATPase. Na-K-ATPase comprises α- and β-subunits that assemble in the endoplasmic reticulum (ER) before trafficking to the plasma membrane. Attachment of H. pylori to gastric epithelial cells increased Na-K-ATPase ubiquitylation, decreased its surface and total levels, and impaired ion balance. H. pylori did not alter degradation of plasmalemma-resident Na-K-ATPase subunits or their mRNA levels. Infection decreased association of α- and β-subunits with ER chaperone BiP and impaired assembly of α/β-heterodimers, as was revealed by quantitative mass spectrometry and immunoblotting of immunoprecipitated complexes. The total level of BiP was not altered, and the decrease in interaction with BiP was not observed for other BiP client proteins. The H. pylori-induced decrease in Na-K-ATPase was prevented by BiP overexpression, stopping protein synthesis, or inhibiting proteasomal, but not lysosomal, protein degradation. The results indicate that H. pylori impairs chaperone-assisted maturation of newly made Na-K-ATPase subunits in the ER independently of a generalized ER stress and induces their ubiquitylation and proteasomal degradation. The decrease in Na-K-ATPase levels is also seen in vivo in the stomachs of gerbils and chronically infected children. Further understanding of H. pylori-induced Na-K-ATPase degradation will provide insights for protection against advanced disease.NEW & NOTEWORTHY This work provides evidence that Helicobacter pylori decreases levels of Na-K-ATPase, a vital transport enzyme, in gastric epithelia, both in acutely infected cultured cells and in chronically infected patients and animals. The bacteria interfere with BiP-assisted folding of newly-made Na-K-ATPase subunits in the endoplasmic reticulum, accelerating their ubiquitylation and proteasomal degradation and decreasing efficiency of the assembly of native enzyme. Decreased Na-K-ATPase expression contributes to H. pylori-induced gastric injury.

Imaging-AMARETTO: An Imaging Genomics Software Tool to Interrogate Multiomics Networks for Relevance to Radiography and Histopathology Imaging Biomarkers of Clinical Outcomes

PURPOSE

The availability of increasing volumes of multiomics, imaging, and clinical data in complex diseases such as cancer opens opportunities for the formulation and development of computational imaging genomics methods that can link multiomics, imaging, and clinical data.

METHODS

Here, we present the Imaging-AMARETTO algorithms and software tools to systematically interrogate regulatory networks derived from multiomics data within and across related patient studies for their relevance to radiography and histopathology imaging features predicting clinical outcomes.

RESULTS

To demonstrate its utility, we applied Imaging-AMARETTO to integrate three patient studies of brain tumors, specifically, multiomics with radiography imaging data from The Cancer Genome Atlas (TCGA) glioblastoma multiforme (GBM) and low-grade glioma (LGG) cohorts and transcriptomics with histopathology imaging data from the Ivy Glioblastoma Atlas Project (IvyGAP) GBM cohort. Our results show that Imaging-AMARETTO recapitulates known key drivers of tumor-associated microglia and macrophage mechanisms, mediated by STAT3, AHR, and CCR2, and neurodevelopmental and stemness mechanisms, mediated by OLIG2. Imaging-AMARETTO provides interpretation of their underlying molecular mechanisms in light of imaging biomarkers of clinical outcomes and uncovers novel master drivers, THBS1 and MAP2, that establish relationships across these distinct mechanisms.

CONCLUSION

Our network-based imaging genomics tools serve as hypothesis generators that facilitate the interrogation of known and uncovering of novel hypotheses for follow-up with experimental validation studies. We anticipate that our Imaging-AMARETTO imaging genomics tools will be useful to the community of biomedical researchers for applications to similar studies of cancer and other complex diseases with available multiomics, imaging, and clinical data.

Downregulated AKT-mTOR signaling pathway proteins in dorsolateral prefrontal cortex in Schizophrenia

Abnormal neurotransmission is central to schizophrenia (SZ). Alterations across multiple neurotransmitter systems in SZ suggest that this illness may be associated with dysregulation of core intracellular processes such as signaling pathways that underlie the regulation and integration of these systems. The AKT-mTOR signaling cascade has been implicated in SZ by gene association, postmortem brain and animal studies. AKT and mTOR are serine/threonine kinases which play important roles in cell growth, proliferation, survival, and differentiation. Both AKT and mTOR require phosphorylation at specific sites for their complete activation. mTOR forms two functionally distinct multiprotein complexes, mTOR Complex 1 (mTORC1) and Complex 2 (mTORC2). mTORC1 mediates ribosome biogenesis, protein translation, and autophagy, whereas mTORC2 contributes to actin dynamics. Altered protein synthesis and actin dynamics can lead to an abnormal neuronal morphology resulting in deficits in learning and memory. Currently, there is a lack of direct evidence to support the hypothesis of disrupted mTOR signaling in SZ, and we have addressed this by characterizing this signaling pathway in SZ brain. We found a reduction in AKT and mTOR protein expression and/or phosphorylation state in dorsolateral prefrontal cortex (DLPFC) from 22 pairs of SZ and matched comparison subjects. We also found reduced protein expression of GβL, a subunit protein common to both mTOR complexes. We further investigated mTOR complex-specific subunit composition and phosphorylation state, and found abnormal mTOR expression in both complexes in SZ DLPFC. These findings provide evidence that proteins associated with the AKT-mTOR signaling cascade are downregulated in SZ DLPFC.

Effect for Human Genomic Variation During the BMP4-Induced Conversion From Pluripotent Stem Cells to Trophoblast

The role of genomic variation in differentiation is currently not well understood. Here, the genomic variations were determined with the whole-genome sequencing for three pairs of pluripotent stem cell lines and their corresponding BMP4-induced trophoblast cell lines. We identified ∼3,500 single nucleotide variations and ∼4,500 indels by comparing the genome sequenced data between the stem cell lines and the matched BMP4-induced trophoblast cell lines and annotated them by integrating the epigenomic and transcriptomic datasets. Relatively, introns enrich more variations. We found ∼45% (42 genes) of the differentially expressed genes in trophoblasts that associate genomic variations. Six variations, located at transcription factor binding sites where H3K4me3 and H3K27ac are enriched in both H1 and H1_BMP4, were identified. The epigenetic status around the genomic variations in H1 was similar to that in H1_BMP4. This means that the variation-associated gene’s expression change can not be attributed to epigenetic alteration. The genes associated with the six variations were upregulated in differentiation. We inferred that during the differentiation, an increased in the expression level of the MEF2C gene is due to a genomic variation in chromosomes 5: 88179358 A > G, which is at a binding site of TFs KLF16, NR2C2, and ZNF740 to MEF2C. Allele G shows a higher affinity to the TFs in the induced cells. The increased expression of MEF2C leads to an increased expression of TF MEF2C’s target genes, subsequently affecting the differentiation. Although genomic variation should not be a dominant factor in differentiation, we believe that genomic variation could indeed play a role in the differentiation from stem cells into trophoblast.

Maternal blood EBF1-based microRNA transcripts as biomarkers for detecting risk of spontaneous preterm birth: a nested case-control study

Objective: Both genetic variants and maternal blood mRNA levels of EBF1 gene have been linked to sPTB. Animal and human studies suggest that specific EBF1-based miRNAs are involved in various physiological and pathophysiological processes. However, to date, we did not find any reports of EBF1-based miRNAs or miRNA transcripts in relation to sPTB. We therefore aimed to examine whether maternal blood early B cell factor 1 (EBF1) gene-based microRNA (miRNA) transcripts can be used for detecting risk of spontaneous preterm birth (sPTB).Methods: We conducted a nested case-control study within a Canadian cohort consisting of 1878 singleton pregnancies enrolled from May 2008 to December 2010 in Calgary, Alberta, Canada. We used a public gene expression dataset (GSE59491) derived from maternal blood in trimesters 2-3 that included women with sPTB (n = 51) and term births (n = 106) matched for maternal age, race/ethnicity, pre-pregnancy body mass index, smoking during pregnancy, and parity within the Canadian cohort. Two bioinformatics tools, miRWalk and STarMirDB, with different algorithms were applied to retrieve miRNA transcripts that putatively target the EBF1 gene (i.e. EBF1-based). Limma moderated t-tests were used to examine differentially expressed (DE) miRNA transcripts (sPTB vs term) within trimesters. Logistic regression models with miRNA transcript tertiles were applied to assess threshold associations between candidate miRNA transcripts' levels and sPTB. Receiver operating characteristic (ROC) analyses were used to identify the maximum Youden Index and its corresponding optimal sensitivity/specificity cut-point of EBF1-based miRNA transcripts for classifying sPTB, and to compare the classification performance of a linear combination (score) of miRNA transcripts with that of individual miRNA transcripts. A five-fold cross-validation was applied to examine the possible overfitting problem of the final ROC model.Results: Four maternal blood EBF1-based miRNA transcripts (MIR4266, MIR1251, MIR601, MIR3612) in the 3rd trimester were significantly associated with sPTB. The odds ratios (95%CIs) for highest versus lowest tertile of the four miRNA transcripts were 3.01-5.25(1.21-13.14, p ≤ .018). The combined 4-miRNA transcripts' score significantly improved the classification of sPTB compared to individual miRNA transcripts (AUC increased from 0.65-0.69 to 0.82, p ≤ .0034) and showed a sensitivity for sPTB of 0.81 and a specificity of 0.72. The final ROC model of the EBF1-based 4 miRNA transcripts' score in cases and controls had no significant overfitting issue.Conclusions: Maternal blood EBF1-based miRNA transcripts may, along with other biomarkers, be useful in screening for sPTB risk in 3rd trimester. Our results also provide clues for further study of potential molecular mechanisms underlying the relationship between EBF1 gene and sPTB, e.g. connecting genetic variants, mRNA expression, and miRNA regulation.

Overexpression of transcription factor Foxa1 and target genes remediate therapeutic protein production bottlenecks in Chinese hamster ovary cells

Despite extensive research conducted to increase protein production from Chinese hamster ovary (CHO) cells, cellular bottlenecks often remain, hindering high yields. In this study, a transcriptomic analysis led to the identification of 32 genes that are consistently upregulated in high producer clones and thus might mediate high productivity. Candidate genes were associated with functions such as signaling, protein folding, cytoskeleton organization, and cell survival. We focused on two engineering targets, Erp27, which binds unfolded proteins and the Erp57 disulfide isomerase in the endoplasmic reticulum, and Foxa1, a pioneering transcription factor involved in organ development. Erp27 moderate overexpression increased production of an easy-to-express antibody, whereas Erp27 and Erp57 co-overexpression increased cell density, viability, and the yield of difficult-to-express proteins. Foxa1 overexpression increased cell density, cell viability, and easy- and difficult-to-express protein yields, whereas it decreased reactive oxygen species late in fed-batch cultures. Foxa1 overexpression upregulated two other candidate genes that increased the production of difficult- and/or easy-to-express proteins, namely Ca3, involved in protecting cells from oxidative stress, and Tagap, involved in signaling and cytoskeleton remodeling. Overall, several genes allowing to overcome CHO cell bottlenecks were identified, including Foxa1, which mediated multiple favorable metabolic changes that improve therapeutic protein yields.

Exercise rejuvenates quiescent skeletal muscle stem cells in old mice through restoration of Cyclin D1

Aging impairs tissue repair. This is pronounced in skeletal muscle, whose regeneration by muscle stem cells (MuSCs) is robust in young adult animals but inefficient in older organisms. Despite this functional decline, old MuSCs are amenable to rejuvenation through strategies that improve the systemic milieu, such as heterochronic parabiosis. One such strategy, exercise, has long been appreciated for its benefits on healthspan, but its effects on aged stem cell function in the context of tissue regeneration are incompletely understood. Here we show that exercise in the form of voluntary wheel running accelerates muscle repair in old animals and improves old MuSC function. Through transcriptional profiling and genetic studies, we discovered that the restoration of old MuSC activation ability hinges on restoration of Cyclin D1, whose expression declines with age in MuSCs. Pharmacologic studies revealed that Cyclin D1 maintains MuSC activation capacity by repressing TGFβ signaling. Taken together, these studies demonstrate that voluntary exercise is a practicable intervention for old MuSC rejuvenation. Furthermore, this work highlights the distinct role of Cyclin D1 in stem cell quiescence.

Ligand-Mediated Targeting of Cytokine Interleukin-27 Enhances Its Bioactivity In Vivo

We have examined the role of a novel targeted cytokine, interleukin-27 (IL-27), modified at the C terminus with a dual targeting and therapeutic heptapeptide, in treating prostate cancer. IL-27 has shown promise in halting tumor growth and mediating tumor regression in several cancer models, including prostate cancer. We describe our findings on the effects of targeted IL-27 gene delivery on prostate cancer cells in vitro and in vivo and how the targeting enhances bioactivity of the IL-27 cytokine. We applied the IL-27 gene delivery protocol utilizing sonoporation (sonodelivery) with the goal of reducing prostate tumor growth in an immunocompetent TC2R C57/BL6 model. The reduction in tumor growth and effector cellular profiles implicate targeted IL-27 as more effective than an untargeted version of IL-27 in promoting bioactivity, as assessed by STAT1 and IFN-γ reporter genes. Moreover, enhanced antitumor effects and significantly higher accumulation of natural killer T (NKT) and CD8 effector cells in the tumors were observed. These results support a novel IL-27-based targeting strategy that is promising since it shows improved therapeutic efficacy while utilizing simple and effective sonodelivery methods.

Minority Gene Expression Profiling: Probing the Genetic Signatures of Pathogenesis Using Ribosome Profiling

Minority Gene Expression Profiling (MGEP) refers to a scenario where the expression profiles of specific genes of interest are concentrated in a small cellular pool that is embedded within a larger, non-expressive pool. An example of this is the analysis of disease-related genes within sub-populations of blood or biopsied tissues. These systems are characterized by low signal-to-noise ratios that make it difficult, if not impossible, to uncover the desired signatures of pathogenesis in the absence of lengthy, and often problematic, technical manipulations. We have adapted ribosome profiling (RP) workflows from the Illumina to the Ion Proton platform and used them to analyze signatures of pathogenesis in an MGEP model system consisting of human cells eliciting <3% productive dengue infection. We find that RP is powerful enough to identify relevant responses of differentially expressed genes, even in the presence of significant noise. We discuss how to deal with sources of unwanted variation, and propose ways to further improve this powerful approach to the study of pathogenic signatures within MGEP systems.

mitoXplorer, a visual data mining platform to systematically analyze and visualize mitochondrial expression dynamics and mutations

Mitochondria participate in metabolism and signaling. They adapt to the requirements of various cell types. Publicly available expression data permit to study expression dynamics of genes with mitochondrial function (mito-genes) in various cell types, conditions and organisms. Yet, we lack an easy way of extracting these data for mito-genes. Here, we introduce the visual data mining platform mitoXplorer, which integrates expression and mutation data of mito-genes with a manually curated mitochondrial interactome containing ∼1200 genes grouped in 38 mitochondrial processes. User-friendly analysis and visualization tools allow to mine mitochondrial expression dynamics and mutations across various datasets from four model species including human. To test the predictive power of mitoXplorer, we quantify mito-gene expression dynamics in trisomy 21 cells, as mitochondrial defects are frequent in trisomy 21. We uncover remarkable differences in the regulation of the mitochondrial transcriptome and proteome in one of the trisomy 21 cell lines, caused by dysregulation of the mitochondrial ribosome and resulting in severe defects in oxidative phosphorylation. With the newly developed Fiji plugin mitoMorph, we identify mild changes in mitochondrial morphology in trisomy 21. Taken together, mitoXplorer (http://mitoxplorer.ibdm.univ-mrs.fr) is a user-friendly, web-based and freely accessible software, aiding experimental scientists to quantify mitochondrial expression dynamics.

Melanogenic Difference Consideration in Ethnic Skin Type: A Balance Approach Between Skin Brightening Applications and Beneficial Sun Exposure

Human skin demonstrates a striking variation in tone and color that is evident among multiple demographic populations. Such characteristics are determined predominantly by the expression of the genes controlling the quantity and quality of melanin, which can alter significantly due to the presence of small nucleotide polymorphism affecting various steps of the melanogenesis process and generally linked to the lighter skin phenotypes. Genetically determined, constitutive skin color is additionally complemented by the facultative melanogenesis and tanning responses; with high levels of melanin and melanogenic factors broadly recognized to have a protective effect against the UVR-induced molecular damage in darker skin. Long-term sun exposure, together with a genetic makeup responsible for the ability to tan or the activity of constitutive melanogenic factors, triggers defects in pigmentation across all ethnic skin types. However, sun exposure also has well documented beneficial effects that manifest at both skin homeostasis and the systemic level, such as synthesis of vitamin D, which is thought to be less efficient in the presence of high levels of melanin or potentially linked to the polymorphism in the genes responsible for skin darkening triggered by UVR. In this review, we discuss melanogenesis in a context of constitutive pigmentation, defined by gene polymorphism in ethnic skin types, and facultative pigmentation that is not only associated with the capacity to protect the skin against photo-damage but could also have an impact on vitamin D synthesis through gene polymorphism. Modulating the activities of melanogenic genes, with the focus on the markers specifically altered by polymorphism combined with differential requirements of sun exposure in ethnic skin types, could enhance the applications of already existing skin brightening factors and provide a novel approach toward improved skin tone and health in personalized skincare.

Temporal changes in the gene expression heterogeneity during brain development and aging

Cells in largely non-mitotic tissues such as the brain are prone to stochastic (epi-)genetic alterations that may cause increased variability between cells and individuals over time. Although increased inter-individual heterogeneity in gene expression was previously reported, whether this process starts during development or if it is restricted to the aging period has not yet been studied. The regulatory dynamics and functional significance of putative aging-related heterogeneity are also unknown. Here we address these by a meta-analysis of 19 transcriptome datasets from three independent studies, covering diverse human brain regions. We observed a significant increase in inter-individual heterogeneity during aging (20 + years) compared to postnatal development (0 to 20 years). Increased heterogeneity during aging was consistent among different brain regions at the gene level and associated with lifespan regulation and neuronal functions. Overall, our results show that increased expression heterogeneity is a characteristic of aging human brain, and may influence aging-related changes in brain functions.

Expression of DOCK10.1 protein revealed with a specific antiserum: insights into regulation of first exon isoforms of DOCK10

DOCK10, a guanine-nucleotide exchange factor (GEF) for Rho GTPases, represents the example of a gene that gives rise to alternative first exon mRNA isoforms, named DOCK10.1 and DOCK10.2. Expression of human DOCK10.2 protein in cell lines, and its induction by interleukin-4 (IL-4) in normal B lymphocytes and chronic lymphocytic leukemia (CLL) cells, were previously demonstrated using an antiserum raised against a peptide encoded by sequences from exon 1.2. Here, expression of human DOCK10.1 protein was demonstrated using an antiserum raised against a peptide encoded by sequences from exon 1.1. Specificity of the DOCK10.1 and DOCK10.2 antisera for their respective isoforms was demonstrated using transfected human 293 T cells. Their specificity for endogenous DOCK10 was strongly suggested by the high significance of the correlations between the levels of their expected signals at the molecular size of 250 kDa and the levels of DOCK10.1 and DOCK10.2 mRNAs, respectively, in human hematopoietic cell lines. Specificity of the DOCK10.1 antiserum for DOCK10 was also demostrated in mouse using the DOCK10 knockout model. The DOCK10.1 protein was induced by IL-4 in CLL cells, which demonstrates that the mechanism by which IL-4 regulates DOCK10 is not isoform-specific. Last, to get insights into differential regulation of the DOCK10 isoforms, their protein levels in cell lines were compared with their gene expression profiles retrieved from the Cancer Cell Line Encyclopedia (CCLE), leading to the identification of BCL3 and KLF12 as potential transcriptional regulators of DOCK10.1 and DOCK10.2, respectively.

DNA Methylation Changes are Associated with the Programming of White Adipose Tissue Browning Features by Resveratrol and Nicotinamide Riboside Neonatal Supplementations in Mice

Neonatal supplementation with resveratrol (RSV) or nicotinamide riboside (NR) programs in male mice brown adipocyte-like features in white adipose tissue (WAT browning) together with improved metabolism in adulthood. We tested the involvement in this programming of long-term epigenetic changes in two browning-related genes that are overexpressed in WAT of supplemented mice, Slc27a1 and Prdm16. Suckling mice received orally the vehicle, RSV or NR from postnatal days 2-to-20. After weaning (d21) onto a chow diet, male mice were habituated to a normal-fat diet (NFD) starting d75, and split on d90 into continuation on the NFD or switching to a high-fat diet (HFD) until euthanization on d164. CpG methylation by bisulfite-sequencing was analyzed on inguinal WAT. Both treatments modified methylation marks in Slc27a1 and Prdm16 and the HFD-dependent dynamics of these marks in the adult WAT, with distinct and common effects. The treatments also affected gene expression of de novo DNA methylases in WAT of young animals (euthanized at d35 in independent experiments). Studies in 3T3-L1 adipocytes indicated the direct effects of RSV and NR on the DNA methylation machinery and favoring browning features. The results support epigenetic effects being involved in WAT programming by neonatal RSV or NR supplementation in male mice.

Cell-autonomous light sensitivity via Opsin3 regulates fuel utilization in brown adipocytes

Opsin3 (Opn3) is a transmembrane heptahelical G protein–coupled receptor (GPCR) with the potential to produce a nonvisual photoreceptive effect. Interestingly, anatomical profiling of GPCRs reveals that Opn3 mRNA is highly expressed in adipose tissue. The photosensitive functions of Opn3 in mammals are poorly understood, and whether Opn3 has a role in fat is entirely unknown. In this study, we found that Opn3-knockout (Opn3-KO) mice were prone to diet-induced obesity and insulin resistance. At the cellular level, Opn3-KO brown adipocytes cultured in darkness had decreased glucose uptake and lower nutrient-induced mitochondrial respiration than wild-type (WT) cells. Light exposure promoted mitochondrial activity and glucose uptake in WT adipocytes but not in Opn3-KO cells. Brown adipocytes carrying a defective mutation in Opn3’s putative G protein–binding domain also exhibited a reduction in glucose uptake and mitochondrial respiration in darkness. Using RNA-sequencing, we identified several novel light-sensitive and Opn3-dependent molecular signatures in brown adipocytes. Importantly, direct exposure of brown adipose tissue (BAT) to light in living mice significantly enhanced thermogenic capacity of BAT, and this effect was diminished in Opn3-KO animals. These results uncover a previously unrecognized cell-autonomous, light-sensing mechanism in brown adipocytes via Opn3-GPCR signaling that can regulate fuel metabolism and mitochondrial respiration. Our work also provides a molecular basis for developing light-based treatments for obesity and its related metabolic disorders.

Brown adipose tissue plays a pivotal role in energy homeostasis and serves as a metabolic sink for glucose and fatty acid. This study demonstrates a novel light-sensing mechanism in mice via the photoreceptor Opsin3 that regulates fuel utilization and mitochondrial activity of brown adipocytes.

ER-resident sensor PERK is essential for mitochondrial thermogenesis in brown adipose tissue

Mitochondria play a central role in the function of brown adipocytes (BAs). Although mitochondrial biogenesis, which is indispensable for thermogenesis, is regulated by coordination between nuclear DNA transcription and mitochondrial DNA transcription, the molecular mechanisms of mitochondrial development during BA differentiation are largely unknown. Here, we show the importance of the ER-resident sensor PKR-like ER kinase (PERK) in the mitochondrial thermogenesis of brown adipose tissue. During BA differentiation, PERK is physiologically phosphorylated independently of the ER stress. This PERK phosphorylation induces transcriptional activation by GA-binding protein transcription factor α subunit (GABPα), which is required for mitochondrial inner membrane protein biogenesis, and this novel role of PERK is involved in maintaining the body temperatures of mice during cold exposure. Our findings demonstrate that mitochondrial development regulated by the PERK-GABPα axis is indispensable for thermogenesis in brown adipose tissue.

Weighted gene co-expression network analysis identifies specific modules and hub genes related to subsyndromal symptomatic depression

Objectives: The identification of the potential molecule targets for subsyndromal symptomatic depression (SSD) is critical for improving the effective clinical treatment on the mental illness. In the current study, we mined the genome-wide expression profiling and investigated the novel biological pathways associated with SSD.Methods: Expression of differentially expressed genes (DEGs) were analysed with microarrays of blood tissue cohort of eight SSD patients and eight healthy subjects. The gene co-expression is calculated by WGCNA, an R package software. The function of the genes was annotated by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis.Results: We identified 11 modules from the 9,427 DEGs. Three co-expression modules (blue, cyan and red) showed striking correlation with the phenotypic trait between SSD and healthy controls. Gene ontology and KEGG pathway analysis demonstrated that the function of these three modules was enriched with the pathway of inflammatory response and type II diabetes mellitus. Finally, three hub genes, NT5DC1, SGSM2 and MYCBP, were identified from the blue module as significant genes.Conclusions: This first blood gene expression study in SSD observed distinct patterns between cases and controls which may provide novel insight into understanding the molecular mechanisms of SSD.

Whole exome sequencing in familial isolated primary hyperparathyroidism

PURPOSE

Familial isolated hyperparathyroidism (FIHP) is a rare inherited disease accounting for 1% of all cases of primary hyperparathyroidism (PHPT). It is genetically heterogeneous being associated with mutations in different genes, including MEN1, CDC73, CASR, and recently GCM2. The aim of the study was to further investigate the molecular pathogenesis in Italian FIHP kindreds.

METHODS

We used whole exome sequencing (WES) in the probands of seven unrelated FIHP kindreds. We carried out a separate family-based exome analysis in a large family characterized by the co-occurrence of PHPT with multiple tumors apparently unrelated to the disease. Selected variants were also screened in 18 additional FIHP kindreds. The clinical, biochemical, and pathological characteristics of the families were also investigated.

RESULTS

Three different variants in GCM2 gene were found in two families, but only one (p.Tyr394Ser), already been shown to be pathogenic in vitro, segregated with the disease. Six probands carried seven heterozygous missense mutations segregating with the disease in the FAT3, PARK2, HDAC4, ITPR2 and TBCE genes. A genetic variant in the APC gene co-segregating with PHPT (p.Val530Ala) was detected in a family whose affected relatives had additional tumors, including colonic polyposis.

CONCLUSION

We confirm the role of GCM2 germline mutations in the pathogenesis of FIHP, although at a lower rate than in the previous WES study. Further studies are needed to establish the prevalence and the role in the predisposition to FIHP of the novel variants in additional genes.

The Role of Dynamic miRISC During Neuronal Development

Activity-dependent protein synthesis plays an important role during neuronal development by fine-tuning the formation and function of neuronal circuits. Recent studies have shown that miRNAs are integral to this regulation because of their ability to control protein synthesis in a rapid, specific and potentially reversible manner. miRNA mediated regulation is a multistep process that involves inhibition of translation before degradation of targeted mRNA, which provides the possibility to store and reverse the inhibition at multiple stages. This flexibility is primarily thought to be derived from the composition of miRNA induced silencing complex (miRISC). AGO2 is likely the only obligatory component of miRISC, while multiple RBPs are shown to be associated with this core miRISC to form diverse miRISC complexes. The formation of these heterogeneous miRISC complexes is intricately regulated by various extracellular signals and cell-specific contexts. In this review, we discuss the composition of miRISC and its functions during neuronal development. Neurodevelopment is guided by both internal programs and external cues. Neuronal activity and external signals play an important role in the formation and refining of the neuronal network. miRISC composition and diversity have a critical role at distinct stages of neurodevelopment. Even though there is a good amount of literature available on the role of miRNAs mediated regulation of neuronal development, surprisingly the role of miRISC composition and its functional dynamics in neuronal development is not much discussed. In this article, we review the available literature on the heterogeneity of the neuronal miRISC composition and how this may influence translation regulation in the context of neuronal development.

Statistical analysis of spatial expression patterns for spatially resolved transcriptomic studies

Identifying genes that display spatial expression patterns in spatially resolved transcriptomic studies is an important first step toward characterizing the spatial transcriptomic landscape of complex tissues. Here we present a statistical method, SPARK, for identifying spatial expression patterns of genes in data generated from various spatially resolved transcriptomic techniques. SPARK directly models spatial count data through generalized linear spatial models. It relies on recently developed statistical formulas for hypothesis testing, providing effective control of type I errors and yielding high statistical power. With a computationally efficient algorithm, which is based on penalized quasi-likelihood, SPARK is also scalable to datasets with tens of thousands of genes measured on tens of thousands of samples. Analyzing four published spatially resolved transcriptomic datasets using SPARK, we show it can be up to ten times more powerful than existing methods and disclose biological discoveries that otherwise cannot be revealed by existing approaches.

TrkB Signaling Influences Gene Expression in Cortistatin-Expressing Interneurons

Brain-derived neurotrophic factor (BDNF) signals through its cognate receptor tropomyosin receptor kinase B (TrkB) to promote the function of several classes of inhibitory interneurons. We previously reported that loss of BDNF-TrkB signaling in cortistatin (Cort)-expressing interneurons leads to behavioral hyperactivity and spontaneous seizures in mice. We performed bulk RNA sequencing (RNA-seq) from the cortex of mice with disruption of BDNF-TrkB signaling in cortistatin interneurons, and identified differential expression of genes important for excitatory neuron function. Using translating ribosome affinity purification and RNA-seq, we define a molecular profile for Cort-expressing inhibitory neurons and subsequently compare the translatome of normal and TrkB-depleted Cort neurons, revealing alterations in calcium signaling and axon development. Several of the genes enriched in Cort neurons and differentially expressed in TrkB-depleted neurons are also implicated in autism and epilepsy. Our findings highlight TrkB-dependent molecular pathways as critical for the maturation of inhibitory interneurons and support the hypothesis that loss of BDNF signaling in Cort interneurons leads to altered excitatory/inhibitory balance.

Machine learning approaches and databases for prediction of drug-target interaction: a survey paper

The task of predicting the interactions between drugs and targets plays a key role in the process of drug discovery. There is a need to develop novel and efficient prediction approaches in order to avoid costly and laborious yet not-always-deterministic experiments to determine drug–target interactions (DTIs) by experiments alone. These approaches should be capable of identifying the potential DTIs in a timely manner. In this article, we describe the data required for the task of DTI prediction followed by a comprehensive catalog consisting of machine learning methods and databases, which have been proposed and utilized to predict DTIs. The advantages and disadvantages of each set of methods are also briefly discussed. Lastly, the challenges one may face in prediction of DTI using machine learning approaches are highlighted and we conclude by shedding some lights on important future research directions.

Protein expression of prenyltransferase subunits in postmortem schizophrenia dorsolateral prefrontal cortex

The pathophysiology of schizophrenia includes altered neurotransmission, dysregulated intracellular signaling pathway activity, and abnormal dendritic morphology that contribute to deficits of synaptic plasticity in the disorder. These processes all require dynamic protein-protein interactions at cell membranes. Lipid modifications target proteins to membranes by increasing substrate hydrophobicity by the addition of a fatty acid or isoprenyl moiety, and recent evidence suggests that dysregulated posttranslational lipid modifications may play a role in multiple neuropsychiatric disorders, including schizophrenia. Consistent with these emerging findings, we have recently reported decreased protein S-palmitoylation in schizophrenia. Protein prenylation is a lipid modification that occurs upstream of S-palmitoylation on many protein substrates, facilitating membrane localization and activity of key intracellular signaling proteins. Accordingly, we hypothesized that, in addition to palmitoylation, protein prenylation may be abnormal in schizophrenia. To test this, we assayed protein expression of the five prenyltransferase subunits (FNTA, FNTB, PGGT1B, RABGGTA, and RABGGTB) in postmortem dorsolateral prefrontal cortex from patients with schizophrenia and paired comparison subjects (n = 13 pairs). We found decreased levels of FNTA (14%), PGGT1B (13%), and RABGGTB (8%) in schizophrenia. To determine whether upstream or downstream factors may be driving these changes, we also assayed protein expression of the isoprenoid synthases FDPS and GGPS1 and prenylation-dependent processing enzymes RCE and ICMT. We found these upstream and downstream enzymes to have normal protein expression. To rule out effects from chronic antipsychotic treatment, we assayed FNTA, PGGT1B, and RABGGTB in the cortex from rats treated long-term with haloperidol decanoate and found no change in the expression of these proteins. Given the role prenylation plays in localization of key signaling proteins found at the synapse, these data offer a potential mechanism underlying abnormal protein-protein interactions and protein localization in schizophrenia.

Pathway Commons 2019 Update: integration, analysis and exploration of pathway data

Pathway Commons (https://www.pathwaycommons.org) is an integrated resource of publicly available information about biological pathways including biochemical reactions, assembly of biomolecular complexes, transport and catalysis events and physical interactions involving proteins, DNA, RNA, and small molecules (e.g. metabolites and drug compounds). Data is collected from multiple providers in standard formats, including the Biological Pathway Exchange (BioPAX) language and the Proteomics Standards Initiative Molecular Interactions format, and then integrated. Pathway Commons provides biologists with (i) tools to search this comprehensive resource, (ii) a download site offering integrated bulk sets of pathway data (e.g. tables of interactions and gene sets), (iii) reusable software libraries for working with pathway information in several programming languages (Java, R, Python and Javascript) and (iv) a web service for programmatically querying the entire dataset. Visualization of pathways is supported using the Systems Biological Graphical Notation (SBGN). Pathway Commons currently contains data from 22 databases with 4794 detailed human biochemical processes (i.e. pathways) and ∼2.3 million interactions. To enhance the usability of this large resource for end-users, we develop and maintain interactive web applications and training materials that enable pathway exploration and advanced analysis.

RNA-binding protein Musashi2 stabilizing androgen receptor drives prostate cancer progression

The androgen receptor (AR) pathway is critical for prostate cancer carcinogenesis and development; however, after 18‐24 months of AR blocking therapy, patients invariably progress to castration‐resistant prostate cancer (CRPC), which remains an urgent problem to be solved. Therefore, finding key molecules that interact with AR as novel strategies to treat prostate cancer and even CRPC is desperately needed. In the current study, we focused on the regulation of RNA‐binding proteins (RBPs) associated with AR and determined that the mRNA and protein levels of AR were highly correlated with Musashi2 (MSI2) levels. MSI2 was upregulated in prostate cancer specimens and significantly correlated with advanced tumor grades. Downregulation of MSI2 in both androgen sensitive and insensitive prostate cancer cells inhibited tumor formation in vivo and decreased cell growth in vitro, which could be reversed by AR overexpression. Mechanistically, MSI2 directly bound to the 3′‐untranslated region (UTR) of AR mRNA to increase its stability and, thus, enhanced its transcriptional activity. Our findings illustrate a previously unknown regulatory mechanism in prostate cancer cell proliferation regulated by the MSI2‐AR axis and provide novel evidence towards a strategy against prostate cancer.

Transcriptional characterisation of human lung cells identifies novel mesenchymal lineage markers

RATIONALE

The lung mesenchyme gives rise to multiple distinct lineages of cells in the mature respiratory system, including smooth muscle cells of the airway and vasculature. However, a thorough understanding of the specification and mesenchymal cell diversity in the human lung is lacking.

METHODS

We completed single-cell RNA sequencing analysis of fetal human lung tissues. Canonical correlation analysis, clustering, cluster marker gene identification and t-distributed stochastic neighbour embedding representation was performed in Seurat. Cell populations were annotated using ToppFun. Immunohistochemistry and in situ hybridisation were used to validate spatiotemporal gene expression patterns for key marker genes.

RESULTS

We identified molecularly distinct populations representing "committed" fetal human lung endothelial cells, pericytes and smooth muscle cells. Early endothelial lineages expressed "classic" endothelial cell markers (platelet endothelial cell adhesion molecule/CD31 and claudin 5), while pericytes expressed platelet-derived growth factor receptor-β, Thy-1 membrane glycoprotein and basement membrane molecules (collagen IV, laminin and proteoglycans). We observed a large population of "nonspecific" human lung mesenchymal progenitor cells characterised by expression of collagen I and multiple elastin fibre genes (ELN, MFAP2 and FBN1). We closely characterised the diversity of mesenchymal lineages defined by α2-smooth muscle actin (ACTA2) expression. Two cell populations, with the highest levels of ACTA2 transcriptional activity, expressed unique sets of markers associated with airway or vascular smooth muscle cells. Spatiotemporal analysis of these marker genes confirmed early and persistent spatial specification of airway (HHIP, MYLK and IGF1) and vascular (NTRK3 and MEF2C) smooth muscle cells in the developing human lung.

CONCLUSION

Our data suggest that specification of distinct airway and vascular smooth muscle cell phenotypes is established early in development and can be identified using the markers we provide.

TRPV6 as A Target for Cancer Therapy

Two decades ago a class of ion channels, hitherto unsuspected, was discovered. In mammals these Transient Receptor Potential channels (TRPs) have not only expanded in number (to 26 functional channels) but also expanded the view of our interface with the physical and chemical environment. Some are heat and cold sensors while others monitor endogenous and/or exogenous chemical signals. Some TRP channels monitor osmotic potential, and others measure cell movement, stretching, and fluid flow. Many TRP channels are major players in nociception and integration of pain signals. One member of the vanilloid sub-family of channels is TRPV6. This channel is highly selective for divalent cations, particularly calcium, and plays a part in general whole-body calcium homeostasis, capturing calcium in the gut from the diet. TRPV6 can be greatly elevated in a number of cancers deriving from epithelia and considerable study has been made of its role in the cancer phenotype where calcium control is dysfunctional. This review compiles and updates recent published work on TRPV6 as a promising drug target in a number of cancers including those afflicting breast, ovarian, prostate and pancreatic tissues.

Human chromatin remodeler cofactor, RNA interactor, eraser and writer sperm RNAs responding to obesity

In the United States almost 33% of adults are considered obese (BMI > 30 kg/m2). Both animal models and to a lesser extent human studies, have associated BMI, a measure of obesity, with alterations in sperm DNA methylation and RNAs. Sperm RNAs from the Assessment of Multiple Gestations from Ovarian Stimulation trial, were isolated and sequenced. A Generalized Linear Model identified 487 BMI associated human sperm RNA elements (short exon-sized sequences). They partitioned into four patterns; a continual increase with BMI, increase once obese (BMI>30 kg/m2); a steady decrease with BMI; and decrease once overweight (BMI 25 - 30 kg/m2). Gene Ontology revealed a unique relationship between BMI and transcripts associated with chromosome organization, adipogenesis, cellular stress and obesity-related inflammation. Coregulatory networks linked by Chromatin remodeler cofactors, RNA interactors, Erasers and Writers (CREWs) were uncovered to reveal a hierarchical epigenetic response pathway.

New candidate genes for ST-elevation myocardial infarction

BACKGROUND

Despite extensive research in atherosclerosis, the mechanisms of coronary atherothrombosis in ST-elevation myocardial infarction (STEMI) patients are undetermined.

OBJECTIVES

Our aim was to find candidate genes involved in STEMI by analysing leucocyte gene expression in STEMI patients, without the influence of secondary inflammation from innate immunity, which was assumed to be a consequence rather than the cause of coronary atherothrombosis.

METHODS

Fifty-one patients were included at coronary angiography because of STEMI. Arterial blood was sampled in the acute phase (P1), at 24-48 h (P2) and at 3 months (P3). Leucocyte RNA was isolated and gene expression analysis was performed by Affymetrix Human Transcriptome Array 2.0. By omission of up- or downregulated genes at P2, secondary changes from innate immunity were excluded. Genes differentially expressed in P1 when compared to the convalescent sample in P3 were determined as genes involved in STEMI.

RESULTS

Three genes were upregulated at P1 compared to P3; ABCG1 (P = 5.81 × 10^-5 ), RAB20 (P = 3.69 × 10^-5 ) and TMEM2 (P = 7.75 × 10^-6 ) whilst four were downregulated; ACVR1 (P = 9.01 × 10^-5 ), NFATC2IP (P = 8.86 × 10^-5 ), SUN1 (P = 3.87 × 10^-5 ) and TTC9C (P = 7.18 × 10^-6 ). These genes were also highly expressed in carotid atherosclerotic plaques.

CONCLUSIONS

We found seven genes involved in STEMI. The study is unique regarding the blood sampling in the acute phase and omission of secondary expressed genes from innate immunity. However, the results need to be replicated by future studies.

Identification and comprehensive characterization of lncRNAs with copy number variations and their driving transcriptional perturbed subpathways reveal functional significance for cancer

Numerous studies have shown that copy number variation (CNV) in lncRNA regions play critical roles in the initiation and progression of cancer. However, our knowledge about their functionalities is still limited. Here, we firstly provided a computational method to identify lncRNAs with copy number variation (lncRNAs-CNV) and their driving transcriptional perturbed subpathways by integrating multidimensional omics data of cancer. The high reliability and accuracy of our method have been demonstrated. Then, the method was applied to 14 cancer types, and a comprehensive characterization and analysis was performed. LncRNAs-CNV had high specificity in cancers, and those with high CNV level may perturb broad biological functions. Some core subpathways and cancer hallmarks widely perturbed by lncRNAs-CNV were revealed. Moreover, subpathways highlighted the functional diversity of lncRNAs-CNV in various cancers. Survival analysis indicated that functional lncRNAs-CNV could be candidate prognostic biomarkers for clinical applications, such as ST7-AS1, CDKN2B-AS1 and EGFR-AS1. In addition, cascade responses and a functional crosstalk model among lncRNAs-CNV, impacted genes, driving subpathways and cancer hallmarks were proposed for understanding the driving mechanism of lncRNAs-CNV. Finally, we developed a user-friendly web interface-LncCASE (http://bio-bigdata.hrbmu.edu.cn/LncCASE/) for exploring lncRNAs-CNV and their driving subpathways in various cancer types. Our study identified and systematically characterized lncRNAs-CNV and their driving subpathways and presented valuable resources for investigating the functionalities of non-coding variations and the mechanisms of tumorigenesis.

Glucose-dependent control of leucine metabolism by leucyl-tRNA synthetase 1

Despite the importance of glucose and amino acids for energy metabolism, interactions between the two nutrients are not well understood. We provide evidence for a role of leucyl-tRNA synthetase 1 (LARS1) in glucose-dependent control of leucine usage. Upon glucose starvation, LARS1 was phosphorylated by Unc-51 like autophagy activating kinase 1 (ULK1) at the residues crucial for leucine binding. The phosphorylated LARS1 showed decreased leucine binding, which may inhibit protein synthesis and help save energy. Leucine that is not used for anabolic processes may be available for catabolic pathway energy generation. The LARS1-mediated changes in leucine utilization might help support cell survival under glucose deprivation. Thus, depending on glucose availability, LARS1 may help regulate whether leucine is used for protein synthesis or energy production.

TGFB1-Mediated Gliosis in Multiple Sclerosis Spinal Cords Is Favored by the Regionalized Expression of HOXA5 and the Age-Dependent Decline in Androgen Receptor Ligands

In multiple sclerosis (MS) patients with a progressive form of the disease, spinal cord (SC) functions slowly deteriorate beyond age 40. We previously showed that in the SC of these patients, large areas of incomplete demyelination extend distance away from plaque borders and are characterized by a unique progliotic TGFB1 (Transforming Growth Factor Beta 1) genomic signature. Here, we attempted to determine whether region- and age-specific physiological parameters could promote the progression of SC periplaques in MS patients beyond age 40. An analysis of transcriptomics databases showed that, under physiological conditions, a set of 10 homeobox (HOX) genes are highly significantly overexpressed in the human SC as compared to distinct brain regions. Among these HOX genes, a survey of the human proteome showed that only HOXA5 encodes a protein which interacts with a member of the TGF-beta signaling pathway, namely SMAD1 (SMAD family member 1). Moreover, HOXA5 was previously found to promote the TGF-beta pathway. Interestingly, SMAD1 is also a protein partner of the androgen receptor (AR) and an unsupervised analysis of gene ontology terms indicates that the AR pathway antagonizes the TGF-beta/SMAD pathway. Retrieval of promoter analysis data further confirmed that AR negatively regulates the transcription of several members of the TGF-beta/SMAD pathway. On this basis, we propose that in progressive MS patients, the physiological SC overexpression of HOXA5 combined with the age-dependent decline in AR ligands may favor the slow progression of TGFB1-mediated gliosis. Potential therapeutic implications are discussed.

Prechondrogenic ATDC5 Cell Attachment and Differentiation on Graphene Foam; Modulation by Surface Functionalization with Fibronectin

Graphene foam holds promise for tissue engineering applications. In this study, graphene foam was used as a three-dimension scaffold to evaluate cell attachment, cell morphology, and molecular markers of early differentiation. The aim of this study was to determine if cell attachment and elaboration of an extracellular matrix would be modulated by functionalization of graphene foam with fibronectin, an extracellular matrix protein that cells adhere well to, prior to the establishment of three-dimensional cell culture. The molecular dynamic simulation demonstrated that the fibronectin-graphene interaction was stabilized predominantly through interaction between the graphene and arginine side chains of the protein. Quasi-static and dynamic mechanical testing indicated that fibronectin functionalization of graphene altered the mechanical properties of graphene foam. The elastic strength of the scaffold increased due to fibronectin, but the viscoelastic mechanical behavior remained unchanged. An additive effect was observed in the mechanical stiffness when the graphene foam was both coated with fibronectin and cultured with cells for 28 days. Cytoskeletal organization assessed by fluorescence microscopy demonstrated a fibronectin-dependent reorganization of the actin cytoskeleton and an increase in actin stress fibers. Gene expression assessed by quantitative real-time polymerase chain reaction of 9 genes encoding cell attachment proteins (Cd44, Ctnna1, Ctnnb1, Itga3, Itga5, Itgav, Itgb1, Ncam1, Sgce), 16 genes encoding extracellular matrix proteins (Col1a1, Col2a1, Col3a1, Col5a1, Col6a1, Ecm1, Emilin1, Fn1, Hapln1, Lamb3, Postn, Sparc, Spp1, Thbs1, Thbs2, Tnc), and 9 genes encoding modulators of remodeling (Adamts1, Adamts2, Ctgf, Mmp14, Mmp2, Tgfbi, Timp1, Timp2, Timp3) indicated that graphene foam provided a microenvironment conducive to expression of genes that are important in early chondrogenesis. Functionalization of graphene foam with fibronectin modified the cellular response to graphene foam, demonstrated by decreases in relative gene expression levels. These findings illustrate the combinatorial factors of microscale materials properties and nanoscale molecular features to consider in the design of three-dimensional graphene scaffolds for tissue engineering applications.

ATRX In-Frame Fusion Neuroblastoma Is Sensitive to EZH2 Inhibition via Modulation of Neuronal Gene Signatures

ATRX alterations occur at high frequency in neuroblastoma of adolescents and young adults. Particularly intriguing are the large N-terminal deletions of ATRX (Alpha Thalassemia/Mental Retardation, X-linked) that generate in-frame fusion (IFF) proteins devoid of key chromatin interaction domains, while retaining the SWI/SNF-like helicase region. We demonstrate that ATRX IFF proteins are redistributed from H3K9me3-enriched chromatin to promoters of active genes and identify REST as an ATRX IFF target whose activation promotes silencing of neuronal differentiation genes. We further show that ATRX IFF cells display sensitivity to EZH2 inhibitors, due to derepression of neurogenesis genes, including a subset of REST targets. Taken together, we demonstrate that ATRX structural alterations are not loss-of-function and put forward EZH2 inhibitors as a potential therapy for ATRX IFF neuroblastoma.

A Parametric Targetability Evaluation Approach for Vitiligo Proteome Extracted through Integration of Gene Ontologies and Protein Interaction Topologies

Vitiligo is a well-known skin disorder with complex etiology. Vitiligo pathogenesis is multifaceted with many ramifications. A computational systemic path was designed to first propose candidate disease proteins by merging properties from protein interaction networks and gene ontology terms. All in all, 109 proteins were identified and suggested to be involved in the onset of disease or its progression. Later, a composite approach was employed to prioritize vitiligo disease proteins by comparing and benchmarking the properties against standard target identification criteria. This includes sequence-based, structural, functional, essentiality, protein-protein interaction, vulnerability, secretability, assayability, and druggability information. The existing information was seamlessly integrated into efficient pipelines to propose a novel protocol for assessment of targetability of disease proteins. Using the online data resources and the scripting, an illustrative list of 68 potential drug targets was generated for vitiligo. While this list is broadly consistent with the research community's current interest in certain specific proteins, and suggests novel target candidates that may merit further study, it can still be modified to correspond to a user-specific environment, either by adjusting the weights for chosen criteria (i.e., a quantitative approach) or by changing the considered criteria (i.e., a qualitative approach).

Cutaneous Epithelial to Mesenchymal Transition Activator ZEB1 Regulates Wound Angiogenesis and Closure in a Glycemic Status-Dependent Manner

Epithelial to mesenchymal transition (EMT) and wound vascularization are two critical interrelated processes that enable cutaneous wound healing. Zinc finger E-box binding homeobox 1 (ZEB1), primarily studied in the context of tumor biology, is a potent EMT activator. ZEB1 is also known to contribute to endothelial cell survival as well as stimulate tumor angiogenesis. The role of ZEB1 in cutaneous wounds was assessed using Zeb1^+/- mice, as Zeb1^-/- mice are not viable. Quantitative stable isotope labeling by amino acids in cell culture (SILAC) proteomics was used to elucidate the effect of elevated ZEB1, as noted during hyperglycemia. Under different glycemic conditions, ZEB1 binding to E-cadherin promoter was investigated using chromatin immunoprecipitation. Cutaneous wounding resulted in loss of epithelial marker E-cadherin with concomitant gain of ZEB1. The dominant proteins downregulated after ZEB1 overexpression functionally represented adherens junction pathway. Zeb1^+/- mice exhibited compromised wound closure complicated by defective EMT and poor wound angiogenesis. Under hyperglycemic conditions, ZEB1 lost its ability to bind E-cadherin promoter. Keratinocyte E-cadherin, thus upregulated, resisted EMT required for wound healing. Diabetic wound healing was improved in ZEB^+/- as well as in db/db mice subjected to ZEB1 knockdown. This work recognizes ZEB1 as a key regulator of cutaneous wound healing that is of particular relevance to diabetic wound complication.

Building gene regulatory networks from scATAC-seq and scRNA-seq using Linked Self Organizing Maps

Rapid advances in single-cell assays have outpaced methods for analysis of those data types. Different single-cell assays show extensive variation in sensitivity and signal to noise levels. In particular, scATAC-seq generates extremely sparse and noisy datasets. Existing methods developed to analyze this data require cells amenable to pseudo-time analysis or require datasets with drastically different cell-types. We describe a novel approach using self-organizing maps (SOM) to link scATAC-seq regions with scRNA-seq genes that overcomes these challenges and can generate draft regulatory networks. Our SOMatic package generates chromatin and gene expression SOMs separately and combines them using a linking function. We applied SOMatic on a mouse pre-B cell differentiation time-course using controlled Ikaros over-expression to recover gene ontology enrichments, identify motifs in genomic regions showing similar single-cell profiles, and generate a gene regulatory network that both recovers known interactions and predicts new Ikaros targets during the differentiation process. The ability of linked SOMs to detect emergent properties from multiple types of highly-dimensional genomic data with very different signal properties opens new avenues for integrative analysis of heterogeneous data.

Gene regulatory network inference resources: A practical overview

Transcriptional regulation is a fundamental molecular mechanism involved in almost every aspect of life, from homeostasis to development, from metabolism to behavior, from reaction to stimuli to disease progression. In recent years, the concept of Gene Regulatory Networks (GRNs) has grown popular as an effective applied biology approach for describing the complex and highly dynamic set of transcriptional interactions, due to its easy-to-interpret features. Since cataloguing, predicting and understanding every GRN connection in all species and cellular contexts remains a great challenge for biology, researchers have developed numerous tools and methods to infer regulatory processes. In this review, we catalogue these methods in six major areas, based on the dominant underlying information leveraged to infer GRNs: Coexpression, Sequence Motifs, Chromatin Immunoprecipitation (ChIP), Orthology, Literature and Protein-Protein Interaction (PPI) specifically focused on transcriptional complexes. The methods described here cover a wide range of user-friendliness: from web tools that require no prior computational expertise to command line programs and algorithms for large scale GRN inferences. Each method for GRN inference described herein effectively illustrates a type of transcriptional relationship, with many methods being complementary to others. While a truly holistic approach for inferring and displaying GRNs remains one of the greatest challenges in the field of systems biology, we believe that the integration of multiple methods described herein provides an effective means with which experimental and computational biologists alike may obtain the most complete pictures of transcriptional relationships. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Single-Cell RNA Sequencing of Plant-Associated Bacterial Communities

Plants in soil are not solitary, hence continually interact with and obtain benefits from a community of microbes ("microbiome"). The meta-functional output from the microbiome results from complex interactions among the different community members with distinct taxonomic identities and metabolic capacities. Particularly, the bacterial communities of the root surface are spatially organized structures composed of root-attached biofilms and planktonic cells arranged in complex layers. With the distinct but coordinated roles among the different member cells, bacterial communities resemble properties of a multicellular organism. High throughput sequencing technologies have allowed rapid and large-scale analysis of taxonomic composition and metabolic capacities of bacterial communities. However, these methods are generally unable to reconstruct the assembly of these communities, or how the gene expression patterns in individual cells/species are coordinated within these communities. Single-cell transcriptomes of community members can identify how gene expression patterns vary among members of the community, including differences among different cells of the same species. This information can be used to classify cells based on functional gene expression patterns, and predict the spatial organization of the community. Here we discuss strategies for the isolation of single bacterial cells, mRNA enrichment, library construction, and analysis and interpretation of the resulting single-cell RNA-Seq datasets. Unraveling regulatory and metabolic processes at the single cell level is expected to yield an unprecedented discovery of mechanisms involved in bacterial recruitment, attachment, assembly, organization of the community, or in the specific interactions among the different members of these communities.