Reactome: a database of reactions, pathways and biological processes

A Computational Framework to Identify Transcriptional and Network Differences between Hepatocellular Carcinoma and Normal Liver Tissue and Their Applications in Repositioning Drugs

Hepatocellular carcinoma (HCC) is one of the most common and lethal malignancies worldwide. Although there have been extensive studies on the molecular mechanisms of its carcinogenesis, FDA-approved drugs for HCC are rare. Side effects, development time, and cost of these drugs are the major bottlenecks, which can be partially overcome by drug repositioning. In this study, we developed a computational framework to study the mechanisms of HCC carcinogenesis, in which drug perturbation-induced gene expression signatures were utilized for repositioning of potential drugs. Specifically, we first performed differential expression analysis and coexpression network module analysis on the HCC dataset from The Cancer Genome Atlas database. Differential gene expression analysis identified 1,337 differentially expressed genes between HCC and adjacent normal tissues, which were significantly enriched in functions related to various pathways, including α-adrenergic receptor activity pathway and epinephrine binding pathway. Weighted gene correlation network analysis (WGCNA) suggested that the number of coexpression modules was higher in HCC tissues than in normal tissues. Finally, by correlating differentially expressed genes with drug perturbation-related signatures, we prioritized a few potential drugs, including nutlin and eribulin, for the treatment of hepatocellular carcinoma. The drugs have been reported by a few experimental studies to be effective in killing cancer cells.

Mapping the Secretome of Dental Pulp Stem Cells Under Variable Microenvironmental Conditions

There is substantial evidence supporting the anti-inflammatory and regenerative potential of dental pulp stem cells (DPSCs) through direct cell transplantation or paracrine action. However, DPSC secretome profile remains inadequately studied. This study provides proteomic profiling of the human DPSC secretome by comparatively analysising cell lysates and respective culture supernatants (i.e. conditioned media-CM) under variable oxygen tension conditions (normoxia-20% O₂/CM_Norm vs. hypoxia 2% O₂/CM_Hyp) and/or stimulation with Tumor Necrosis Factor alpha (TNF-α). DPSC-CM samples and respective crude lysates (DPSC-CL) were collected and subjected to SDS-PAGE, followed by LC-MS/MS analysis. The identified proteins were analyzed by Gene Ontology, Reactome, and String databases. The anti-inflammatory properties of DPSC-CMs were validated via an in vitro RAW_246.7 murine macrophages model through evaluation of the expression of pro-and anti-inflammatory markers by real-time PCR. Results showed a total of 2413 proteins identified in CM_Norm, 2479 in CM_Norm+TNF-α, 1642 in CM_Hyp, and 2002 in CM_Hyp + TNF-α samples. CM_Norm contained 122 proteins statistically significantly upregulated compared to the CM_Hyp and involved in pathways related to "ECM organization", "cellular response to hypoxia", and "IL signaling". Functional network analysis showed that TGFβ1, TIMP1 and TIMP2 were key nodes among proteins significantly upregulated in the CM_Norm compared to the CM_Hyp, interacting with more than 10 proteins, each. DPSC-CM application in the in vitro RAW_246.7 model decreased the expression of pro-inflammatory markers (MMP-3, MMP-9, MMP-13, MCP-1), while increasing anti-inflammatory markers (IL-10). Overall, DPSC-CM collected under normoxic conditions is enriched with anti-inflammatory, tissue repair and regenerative factors, which prompts further investigation on its therapeutic applications.

Alternative Splicing Regulation of Low-Frequency Genetic Variants in Exon 2 of TREM2 in Alzheimer's Disease by Splicing-Based Aggregation

TREM2 is among the most well-known Alzheimer’s disease (AD) risk genes; however, the functional roles of its AD-associated variants remain to be elucidated, and most known risk alleles are low-frequency variants whose investigation is challenging. Here, we utilized a splicing-guided aggregation method in which multiple low-frequency TREM2 variants were bundled together to investigate the functional impact of those variants on alternative splicing in AD. We analyzed whole genome sequencing (WGS) and RNA-seq data generated from cognitively normal elderly controls (CN) and AD patients in two independent cohorts, representing three regions in the frontal lobe of the human brain: the dorsolateral prefrontal cortex (CN = 213 and AD = 376), frontal pole (CN = 72 and AD = 175), and inferior frontal (CN = 63 and AD = 157). We observed an exon skipping event in the second exon of TREM2, with that exon tending to be more frequently skipped (p = 0.0012) in individuals having at least one low-frequency variant that caused loss-of-function for a splicing regulatory element. In addition, genes differentially expressed between AD patients with high vs. low skipping of the second exon (i.e., loss of a TREM2 functional domain) were significantly enriched in immune-related pathways. Our splicing-guided aggregation method thus provides new insight into the regulation of alternative splicing of the second exon of TREM2 by low-frequency variants and could be a useful tool for further exploring the potential molecular mechanisms of multiple, disease-associated, low-frequency variants.

An integrative study of genetic variants with brain tissue expression identifies viral etiology and potential drug targets of multiple sclerosis

Multiple sclerosis (MS) is a neuroinflammatory disorder leading to chronic disability. Brain lesions in MS commonly arise in normal-appearing white matter (NAWM). Genome-wide association studies (GWAS) have identified genetic variants associated with MS. Transcriptome alterations have been observed in case-control studies of NAWM. We developed a Cross-Dataset Evaluation (CDE) function for our network-based tool, Edge-Weighted Dense Module Search of GWAS (EW_dmGWAS). We applied CDE to integrate publicly available MS GWAS summary statistics of 41,505 cases and controls with collectively 38 NAWM expression samples, using the human protein interactome as the reference network, to investigate biological underpinnings of MS etiology. We validated the resulting modules with colocalization of GWAS and expression quantitative trait loci (eQTL) signals, using GTEx Consortium expression data for MS-relevant tissues: 14 brain tissues and 4 immune-related tissues. Other network assessments included a drug target query and functional gene set enrichment analysis. CDE prioritized a MS NAWM network containing 55 unique genes. The gene list was enriched (p-value = 2.34 × 10-7) with GWAS-eQTL colocalized genes: CDK4, IFITM3, MAPK1, MAPK3, METTL12B and PIK3R2. The resultant network also included drug signatures of FDA-approved medications. Gene set enrichment analysis revealed the top functional term "intracellular transport of virus", among other viral pathways. We prioritize critical genes from the resultant network: CDK4, IFITM3, MAPK1, MAPK3, METTL12B and PIK3R2. Enriched drug signatures suggest potential drug targets and drug repositioning strategies for MS. Finally, we propose mechanisms of potential MS viral onset, based on prioritized gene set and functional enrichment analysis.

Identification and Characterization of Robust Hepatocellular Carcinoma Prognostic Subtypes Based on an Integrative Metabolite-Protein Interaction Network

Metabolite-protein interactions (MPIs) play key roles in cancer metabolism. However, our current knowledge about MPIs in cancers remains limited due to the complexity of cancer cells. Herein, the authors construct an integrative MPI network and propose a MPI network based hepatocellular carcinoma (HCC) subtyping and mechanism exploration workflow. Based on the expressions of hub proteins on the MPI network, two prognosis-distinctive HCC subtypes are identified. Meanwhile, multiple interdependent features of the poor prognostic subtype are observed, including hypoxia, DNA hypermethylation of metabolic pathways, fatty acid accumulation, immune pathway up-regulation, and exhausted T-cell infiltration. Notably, the immune pathway up-regulation is probably induced by accumulated unsaturated fatty acids which are predicted to interact with multiple immune regulators like SRC and TGFB1. Moreover, based on tumor microenvironment compositions, the poor prognostic subtype is further divided into two sub-populations showing remarkable differences in metabolism. The subtyping shows a strong consistency across multiple HCC cohorts including early-stage HCC. Overall, the authors redefine robust HCC prognosis subtypes and identify potential MPIs linking metabolism to immune regulations, thus promoting understanding and clinical applications about HCC metabolism heterogeneity.

A Novel Pathway Network Analytics Method Based on Graph Theory

A biological pathway is an ordered set of interactions between intracellular molecules having collective activity that impacts cellular function, for example, by controlling metabolite synthesis or by regulating the expression of sets of genes. They play a key role in advanced studies of genomics. However, existing pathway analytics methods are inadequate to extract meaningful biological structure underneath the network of pathways. They also lack automation. Given these circumstances, we have come up with a novel graph theoretic method to analyze disease-related genes through weighted network of biological pathways. The method automatically extracts biological structures, such as clusters of pathways and their relevance, significance of each pathway and gene, and so forth hidden in the complex network. We have demonstrated the effectiveness of the proposed method on a set of genes associated with coronavirus disease 2019.

Identifying causal models between genetically regulated methylation patterns and gene expression in healthy colon tissue

BACKGROUND

DNA methylation is involved in the regulation of gene expression and phenotypic variation, but the inter-relationship between genetic variation, DNA methylation and gene expression remains poorly understood. Here we combine the analysis of genetic variants related to methylation markers (methylation quantitative trait loci: mQTLs) and gene expression (expression quantitative trait loci: eQTLs) with methylation markers related to gene expression (expression quantitative trait methylation: eQTMs), to provide novel insights into the genetic/epigenetic architecture of colocalizing molecular markers.

RESULTS

Normal mucosa from 100 patients with colon cancer and 50 healthy donors included in the Colonomics project have been analyzed. Linear models have been used to find mQTLs and eQTMs within 1 Mb of the target gene. From 32,446 eQTLs previously detected, we found a total of 6850 SNPs, 114 CpGs and 52 genes interrelated, generating 13,987 significant combinations of co-occurring associations (meQTLs) after Bonferromi correction. Non-redundant meQTLs were 54, enriched in genes involved in metabolism of glucose and xenobiotics and immune system. SNPs in meQTLs were enriched in regulatory elements (enhancers and promoters) compared to random SNPs within 1 Mb of genes. Three colorectal cancer GWAS SNPs were related to methylation changes, and four SNPs were related to chemerin levels. Bayesian networks have been used to identify putative causal relationships among associated SNPs, CpG and gene expression triads. We identified that most of these combinations showed the canonical pathway of methylation markers causes gene expression variation (60.1%) or non-causal relationship between methylation and gene expression (33.9%); however, in up to 6% of these combinations, gene expression was causing variation in methylation markers.

CONCLUSIONS

In this study we provided a characterization of the regulation between genetic variants and inter-dependent methylation markers and gene expression in a set of 150 healthy colon tissue samples. This is an important finding for the understanding of molecular susceptibility on colon-related complex diseases.

Spatiotemporal Immune Landscape of Colorectal Cancer Liver Metastasis at Single-Cell Level

Liver metastasis, the leading cause of colorectal cancer mortality, exhibits a highly heterogeneous and suppressive immune microenvironment. Here, we sequenced 97 matched samples by using single-cell RNA-seq and Spatial Transcriptomics. Strikingly, metastatic microenvironment underwent remarkable spatial reprogramming of immunosuppressive cells such as MRC1+ CCL18+ M2-like macrophages. We further developed scMetabolism, a computational pipeline for quantifying single-cell metabolism, and observed that those macrophages harbored enhanced metabolic activity. Interestingly, neoadjuvant chemotherapy could block this status and restore the antitumor immune balance in responsive patients, while the non-responsive patients deteriorated into a more suppressive one. Our work described the immune evolution of metastasis and uncovered the black box of how tumors respond to neoadjuvant chemotherapy.

A Novel Approach for the Discovery of Biomarkers of Radiotherapy Response in Breast Cancer

Radiotherapy (RT) is an important treatment modality for the local control of breast cancer (BC). Unfortunately, not all patients that receive RT will obtain a therapeutic benefit, as cancer cells that either possess intrinsic radioresistance or develop resistance during treatment can reduce its efficacy. For RT treatment regimens to become personalised, there is a need to identify biomarkers that can predict and/or monitor a tumour's response to radiation. Here we describe a novel method to identify such biomarkers. Liquid chromatography-mass spectrometry (LC-MS) was used on conditioned media (CM) samples from a radiosensitive oestrogen receptor positive (ER⁺) BC cell line (MCF-7) to identify cancer-secreted biomarkers which reflected a response to radiation. A total of 33 radiation-induced secreted proteins that had higher (up to 12-fold) secretion levels at 24 h post-2 Gy radiation were identified. Secretomic results were combined with whole-transcriptome gene expression experiments, using both radiosensitive and radioresistant cells, to identify a signature related to intrinsic radiosensitivity. Gene expression analysis assessing the levels of the 33 proteins showed that 5 (YBX3, EIF4EBP2, DKK1, GNPNAT1 and TK1) had higher expression levels in the radiosensitive cells compared to their radioresistant derivatives; 3 of these proteins (DKK1, GNPNAT1 and TK1) underwent in-lab and initial clinical validation. Western blot analysis using CM samples from cell lines confirmed a significant increase in the release of each candidate biomarker from radiosensitive cells 24 h after treatment with a 2 Gy dose of radiation; no significant increase in secretion was observed in the radioresistant cells after radiation. Immunohistochemistry showed that higher intracellular protein levels of the biomarkers were associated with greater radiosensitivity. Intracellular levels were further assessed in pre-treatment biopsy tissues from patients diagnosed with ER⁺ BC that were subsequently treated with breast-conserving surgery and RT. High DKK1 and GNPNAT1 intracellular levels were associated with significantly increased recurrence-free survival times, indicating that these two candidate biomarkers have the potential to predict sensitivity to RT. We suggest that the methods highlighted in this study could be utilised for the identification of biomarkers that may have a potential clinical role in personalising and optimising RT dosing regimens, whilst limiting the administration of RT to patients who are unlikely to benefit.

Mapping gene and gene pathways associated with coronary artery disease: a CARDIoGRAM exome and multi-ancestry UK biobank analysis

Coronary artery disease (CAD) genome-wide association studies typically focus on single nucleotide variants (SNVs), and many potentially associated SNVs fail to reach the GWAS significance threshold. We performed gene and pathway-based association (GBA) tests on publicly available Coronary ARtery DIsease Genome wide Replication and Meta-analysis consortium Exome (n = 120,575) and multi ancestry pan UK Biobank study (n = 442,574) summary data using versatile gene-based association study (VEGAS2) and Multi-marker analysis of genomic annotation (MAGMA) to identify novel genes and pathways associated with CAD. We included only exonic SNVs and excluded regulatory regions. VEGAS2 and MAGMA ranked genes and pathways based on aggregated SNV test statistics. We used Bonferroni corrected gene and pathway significance threshold at 3.0 × 10^-6 and 1.0 × 10^-5, respectively. We also report the top one percent of ranked genes and pathways. We identified 17 top enriched genes with four genes (PCSK9, FAM177, LPL, ARGEF26), reaching statistical significance (p ≤ 3.0 × 10^-6) using both GBA tests in two GWAS studies. In addition, our analyses identified ten genes (DUSP13, KCNJ11, CD300LF/RAB37, SLCO1B1, LRRFIP1, QSER1, UBR2, MOB3C, MST1R, and ABCC8) with previously unreported associations with CAD, although none of the single SNV associations within the genes were genome-wide significant. Among the top 1% non-lipid pathways, we detected pathways regulating coagulation, inflammation, neuronal aging, and wound healing.

Studying the effect of alpha-synuclein and Parkinson's disease linked mutants on inter pathway connectivities

Parkinson's disease is a common neurodegenerative disease. The differential expression of alpha-synuclein within Lewy Bodies leads to this disease. Some missense mutations of alpha-synuclein may resultant in functional aberrations. In this study, our objective is to verify the functional adaptation due to early and late-onset mutation which can trigger or control the rate of alpha-synuclein aggregation. In this regard, we have proposed a computational model to study the difference and similarities among the Wild type alpha-synuclein and mutants i.e., A30P, A53T, G51D, E46K, and H50Q. Evolutionary sequence space analysis is also performed in this experiment. Subsequently, a comparative study has been performed between structural information and sequence space outcomes. The study shows the structural variability among the selected subtypes. This information assists inter pathway modeling due to mutational aberrations. Based on the structural variability, we have identified the protein-protein interaction partners for each protein that helps to increase the robustness of the inter-pathway connectivity. Finally, few pathways have been identified from 12 semantic networks based on their association with mitochondrial dysfunction and dopaminergic pathways.

NICEdrug.ch, a workflow for rational drug design and systems-level analysis of drug metabolism

The discovery of a drug requires over a decade of intensive research and financial investments - and still has a high risk of failure. To reduce this burden, we developed the NICEdrug.ch resource, which incorporates 250,000 bioactive molecules, and studied their enzymatic metabolic targets, fate, and toxicity. NICEdrug.ch includes a unique fingerprint that identifies reactive similarities between drug-drug and drug-metabolite pairs. We validated the application, scope, and performance of NICEdrug.ch over similar methods in the field on golden standard datasets describing drugs and metabolites sharing reactivity, drug toxicities, and drug targets. We use NICEdrug.ch to evaluate inhibition and toxicity by the anticancer drug 5-fluorouracil, and suggest avenues to alleviate its side effects. We propose shikimate 3-phosphate for targeting liver-stage malaria with minimal impact on the human host cell. Finally, NICEdrug.ch suggests over 1300 candidate drugs and food molecules to target COVID-19 and explains their inhibitory mechanism for further experimental screening. The NICEdrug.ch database is accessible online to systematically identify the reactivity of small molecules and druggable enzymes with practical applications in lead discovery and drug repurposing.

Regulatory Network Analysis in Estradiol-Treated Human Endothelial Cells

Background/Aims: Estrogen has been reported to have beneficial effects on vascular biology through direct actions on endothelium. Together with transcription factors, miRNAs are the major drivers of gene expression and signaling networks. The objective of this study was to identify a comprehensive regulatory network (miRNA–transcription factor–downstream genes) that controls the transcriptomic changes observed in endothelial cells exposed to estradiol. Methods: miRNA/mRNA interactions were assembled using our previous microarray data of human umbilical vein endothelial cells (HUVEC) treated with 17β-estradiol (E2) (1 nmol/L, 24 h). miRNA–mRNA pairings and their associated canonical pathways were determined using Ingenuity Pathway Analysis software. Transcription factors were identified among the miRNA-regulated genes. Transcription factor downstream target genes were predicted by consensus transcription factor binding sites in the promoter region of E2-regulated genes by using JASPAR and TRANSFAC tools in Enrichr software. Results: miRNA–target pairings were filtered by using differentially expressed miRNAs and mRNAs characterized by a regulatory relationship according to miRNA target prediction databases. The analysis identified 588 miRNA–target interactions between 102 miRNAs and 588 targets. Specifically, 63 upregulated miRNAs interacted with 295 downregulated targets, while 39 downregulated miRNAs were paired with 293 upregulated mRNA targets. Functional characterization of miRNA/mRNA association analysis highlighted hypoxia signaling, integrin, ephrin receptor signaling and regulation of actin-based motility by Rho among the canonical pathways regulated by E2 in HUVEC. Transcription factors and downstream genes analysis revealed eight networks, including those mediated by JUN and REPIN1, which are associated with cadherin binding and cell adhesion molecule binding pathways. Conclusion: This study identifies regulatory networks obtained by integrative microarray analysis and provides additional insights into the way estradiol could regulate endothelial function in human endothelial cells.

Differential Protein Interactome in Esophageal Squamous Cell Carcinoma Offers Novel Systems Biomarker Candidates with High Diagnostic and Prognostic Performance

Esophageal squamous cell carcinoma (ESCC) is among the most dangerous cancers with high mortality and lack of robust diagnostics and personalized/precision therapeutics. To achieve a systems-level understanding of tumorigenesis, unraveling of variations in the protein interactome and determination of key proteins exhibiting significant alterations in their interaction patterns during tumorigenesis are crucial. To this end, we have described differential protein-protein interactions and differentially interacting proteins (DIPs) in ESCC by utilizing the human protein interactome and transcriptome. Furthermore, DIP-centered modules were analyzed according to their potential in elucidation of disease mechanisms and improvement of efficient diagnostic, prognostic, and treatment strategies. Seven modules were presented as potential diagnostic, and 16 modules were presented as potential prognostic biomarker candidates. Importantly, our findings also suggest that 30 out of the 53 repurposed drugs were noncancer drugs, which could be used in the treatment of ESCC. Interestingly, 25 of these, proposed as novel drug candidates here, have not been previously associated in a context of esophageal cancer. In this context, risperidone and clozapine were validated for their growth inhibitory potential in three ESCC lines. Our findings offer a high potential for the development of innovative diagnostic, prognostic, and therapeutic strategies for further experimental studies in line with predictive diagnostics, targeted prevention, and personalization of medical services in ESCC specifically, and personalized cancer care broadly.

Genome-Wide Association Study Identifies Risk Loci for Cluster Headache

OBJECTIVE

This study was undertaken to identify susceptibility loci for cluster headache and obtain insights into relevant disease pathways.

METHODS

We carried out a genome-wide association study, where 852 UK and 591 Swedish cluster headache cases were compared with 5,614 and 1,134 controls, respectively. Following quality control and imputation, single variant association testing was conducted using a logistic mixed model for each cohort. The 2 cohorts were subsequently combined in a merged analysis. Downstream analyses, such as gene-set enrichment, functional variant annotation, prediction and pathway analyses, were performed.

RESULTS

Initial independent analysis identified 2 replicable cluster headache susceptibility loci on chromosome 2. A merged analysis identified an additional locus on chromosome 1 and confirmed a locus significant in the UK analysis on chromosome 6, which overlaps with a previously known migraine locus. The lead single nucleotide polymorphisms were rs113658130 (p = 1.92 × 10^-17 , odds ratio [OR] = 1.51, 95% confidence interval [CI] = 1.37-1.66) and rs4519530 (p = 6.98 × 10^-17 , OR = 1.47, 95% CI = 1.34-1.61) on chromosome 2, rs12121134 on chromosome 1 (p = 1.66 × 10^-8 , OR = 1.36, 95% CI = 1.22-1.52), and rs11153082 (p = 1.85 × 10^-8 , OR = 1.30, 95% CI = 1.19-1.42) on chromosome 6. Downstream analyses implicated immunological processes in the pathogenesis of cluster headache.

INTERPRETATION

We identified and replicated several genome-wide significant associations supporting a genetic predisposition in cluster headache in a genome-wide association study involving 1,443 cases. Replication in larger independent cohorts combined with comprehensive phenotyping, in relation to, for example, treatment response and cluster headache subtypes, could provide unprecedented insights into genotype-phenotype correlations and the pathophysiological pathways underlying cluster headache. ANN NEUROL 2021;90:193-202.

Control of tissue development and cell diversity by cell cycle-dependent transcriptional filtering

Cell cycle duration changes dramatically during development, starting out fast to generate cells quickly and slowing down over time as the organism matures. The cell cycle can also act as a transcriptional filter to control the expression of long gene transcripts, which are partially transcribed in short cycles. Using mathematical simulations of cell proliferation, we identify an emergent property that this filter can act as a tuning knob to control gene transcript expression, cell diversity, and the number and proportion of different cell types in a tissue. Our predictions are supported by comparison to single-cell RNA-seq data captured over embryonic development. Additionally, evolutionary genome analysis shows that fast-developing organisms have a narrow genomic distribution of gene lengths while slower developers have an expanded number of long genes. Our results support the idea that cell cycle dynamics may be important across multicellular animals for controlling gene transcript expression and cell fate.

Bioinformatic Analysis of Temporal and Spatial Proteome Alternations During Infections

Microbial pathogens have evolved numerous mechanisms to hijack host's systems, thus causing disease. This is mediated by alterations in the combined host-pathogen proteome in time and space. Mass spectrometry-based proteomics approaches have been developed and tailored to map disease progression. The result is complex multidimensional data that pose numerous analytic challenges for downstream interpretation. However, a systematic review of approaches for the downstream analysis of such data has been lacking in the field. In this review, we detail the steps of a typical temporal and spatial analysis, including data pre-processing steps (i.e., quality control, data normalization, the imputation of missing values, and dimensionality reduction), different statistical and machine learning approaches, validation, interpretation, and the extraction of biological information from mass spectrometry data. We also discuss current best practices for these steps based on a collection of independent studies to guide users in selecting the most suitable strategies for their dataset and analysis objectives. Moreover, we also compiled the list of commonly used R software packages for each step of the analysis. These could be easily integrated into one's analysis pipeline. Furthermore, we guide readers through various analysis steps by applying these workflows to mock and host-pathogen interaction data from public datasets. The workflows presented in this review will serve as an introduction for data analysis novices, while also helping established users update their data analysis pipelines. We conclude the review by discussing future directions and developments in temporal and spatial proteomics and data analysis approaches. Data analysis codes, prepared for this review are available from https://github.com/BabuLab-UofR/TempSpac, where guidelines and sample datasets are also offered for testing purposes.

Genetic architectures of proximal and distal colorectal cancer are partly distinct

OBJECTIVE

An understanding of the etiologic heterogeneity of colorectal cancer (CRC) is critical for improving precision prevention, including individualized screening recommendations and the discovery of novel drug targets and repurposable drug candidates for chemoprevention. Known differences in molecular characteristics and environmental risk factors among tumors arising in different locations of the colorectum suggest partly distinct mechanisms of carcinogenesis. The extent to which the contribution of inherited genetic risk factors for CRC differs by anatomical subsite of the primary tumor has not been examined.

DESIGN

To identify new anatomical subsite-specific risk loci, we performed genome-wide association study (GWAS) meta-analyses including data of 48 214 CRC cases and 64 159 controls of European ancestry. We characterised effect heterogeneity at CRC risk loci using multinomial modelling.

RESULTS

We identified 13 loci that reached genome-wide significance (p<5×10-8) and that were not reported by previous GWASs for overall CRC risk. Multiple lines of evidence support candidate genes at several of these loci. We detected substantial heterogeneity between anatomical subsites. Just over half (61) of 109 known and new risk variants showed no evidence for heterogeneity. In contrast, 22 variants showed association with distal CRC (including rectal cancer), but no evidence for association or an attenuated association with proximal CRC. For two loci, there was strong evidence for effects confined to proximal colon cancer.

CONCLUSION

Genetic architectures of proximal and distal CRC are partly distinct. Studies of risk factors and mechanisms of carcinogenesis, and precision prevention strategies should take into consideration the anatomical subsite of the tumour.

Newt: a comprehensive web-based tool for viewing, constructing and analyzing biological maps

MOTIVATION

Visualization of cellular processes and pathways is a key recurring requirement for effective biological data analysis. There is a considerable need for sophisticated web-based pathway viewers and editors operating with widely accepted standard formats, using the latest visualization techniques and libraries.

RESULTS

We developed a web-based tool named Newt for viewing, constructing and analyzing biological maps in standard formats such as SBGN, SBML and SIF.

AVAILABILITY AND IMPLEMENTATION

Newt's source code is publicly available on GitHub and freely distributed under the GNU LGPL. Ample documentation on Newt can be found on http://newteditor.org and on YouTube.

StrongestPath: a Cytoscape application for protein-protein interaction analysis

BACKGROUND

StrongestPath is a Cytoscape 3 application that enables the analysis of interactions between two proteins or groups of proteins in a collection of protein-protein interaction (PPI) network or signaling network databases. When there are different levels of confidence over the interactions, the application is able to process them and identify the cascade of interactions with the highest total confidence score. Given a set of proteins, StrongestPath can extract a set of possible interactions between the input proteins, and expand the network by adding new proteins that have the most interactions with highest total confidence to the current network of proteins. The application can also identify any activating or inhibitory regulatory paths between two distinct sets of transcription factors and target genes. This application can be used on the built-in human and mouse PPI or signaling databases, or any user-provided database for some organism.

RESULTS

Our results on 12 signaling pathways from the NetPath database demonstrate that the application can be used for indicating proteins which may play significant roles in a pathway by finding the strongest path(s) in the PPI or signaling network.

CONCLUSION

Easy access to multiple public large databases, generating output in a short time, addressing some key challenges in one platform, and providing a user-friendly graphical interface make StrongestPath an extremely useful application.

Systems Biomedicine of Primary and Metastatic Colorectal Cancer Reveals Potential Therapeutic Targets

Colorectal cancer (CRC) is one of the major causes of cancer deaths across the world. Patients' survival at time of diagnosis depends mainly on stage of the tumor. Therefore, understanding the molecular mechanisms from low-grade to high-grade stages of cancer that lead to cellular migration from one tissue/organ to another tissue/organ is essential for implementing therapeutic approaches. To this end, we performed a unique meta-analysis flowchart by identifying differentially expressed genes (DEGs) between normal, primary (primary sites), and metastatic samples (Colorectal metastatic lesions in liver and lung) in some Test datasets. DEGs were employed to construct a protein-protein interaction (PPI) network. A smaller network containing 39 DEGs was then extracted from the PPI network whose nodes expression induction or suppression alone or in combination with each other would inhibit tumor progression or metastasis. These DEGs were then verified by gene expression profiling, survival analysis, and multiple Validation datasets. We suggested for the first time that downregulation of mitochondrial genes, including ETHE1, SQOR, TST, and GPX3, would help colorectal cancer cells to produce more energy under hypoxic conditions through mechanisms that are different from "Warburg Effect". Augmentation of given antioxidants and repression of P4HA1 and COL1A2 genes could be a choice of CRC treatment. Moreover, promoting active GSK-3β together with expression control of EIF2B would prevent EMT. We also proposed that OAS1 expression enhancement can induce the anti-cancer effects of interferon-gamma, while suppression of CTSH hinders formation of focal adhesions. ATF5 expression suppression sensitizes cancer cells to anchorage-dependent death signals, while LGALS4 induction recovers cell-cell junctions. These inhibitions and inductions would be another combinatory mechanism that inhibits EMT and cell migration. Furthermore, expression inhibition of TMPO, TOP2A, RFC3, GINS1, and CKS2 genes could prevent tumor growth. Besides, TRIB3 suppression would be a promising target for anti-angiogenic therapy. SORD is a poorly studied enzyme in cancer, found to be upregulated in CRC. Finally, TMEM131 and DARS genes were identified in this study whose roles have never been interrogated in any kind of cancer, neither as a biomarker nor curative target. All the mentioned mechanisms must be further validated by experimental wet-lab techniques.

Evolving scenario of big data and Artificial Intelligence (AI) in drug discovery

The accumulation of massive data in the plethora of Cheminformatics databases has made the role of big data and artificial intelligence (AI) indispensable in drug design. This has necessitated the development of newer algorithms and architectures to mine these databases and fulfil the specific needs of various drug discovery processes such as virtual drug screening, de novo molecule design and discovery in this big data era. The development of deep learning neural networks and their variants with the corresponding increase in chemical data has resulted in a paradigm shift in information mining pertaining to the chemical space. The present review summarizes the role of big data and AI techniques currently being implemented to satisfy the ever-increasing research demands in drug discovery pipelines.

Genetic and behavioral adaptation of Candida parapsilosis to the microbiome of hospitalized infants revealed by in situ genomics, transcriptomics, and proteomics

BACKGROUND

Candida parapsilosis is a common cause of invasive candidiasis, especially in newborn infants, and infections have been increasing over the past two decades. C. parapsilosis has been primarily studied in pure culture, leaving gaps in understanding of its function in a microbiome context.

RESULTS

Here, we compare five unique C. parapsilosis genomes assembled from premature infant fecal samples, three of which are newly reconstructed, and analyze their genome structure, population diversity, and in situ activity relative to reference strains in pure culture. All five genomes contain hotspots of single nucleotide variants, some of which are shared by strains from multiple hospitals. A subset of environmental and hospital-derived genomes share variants within these hotspots suggesting derivation of that region from a common ancestor. Four of the newly reconstructed C. parapsilosis genomes have 4 to 16 copies of the gene RTA3, which encodes a lipid translocase and is implicated in antifungal resistance, potentially indicating adaptation to hospital antifungal use. Time course metatranscriptomics and metaproteomics on fecal samples from a premature infant with a C. parapsilosis blood infection revealed highly variable in situ expression patterns that are distinct from those of similar strains in pure cultures. For example, biofilm formation genes were relatively less expressed in situ, whereas genes linked to oxygen utilization were more highly expressed, indicative of growth in a relatively aerobic environment. In gut microbiome samples, C. parapsilosis co-existed with Enterococcus faecalis that shifted in relative abundance over time, accompanied by changes in bacterial and fungal gene expression and proteome composition.

CONCLUSIONS

The results reveal potentially medically relevant differences in Candida function in gut vs. laboratory environments, and constrain evolutionary processes that could contribute to hospital strain persistence and transfer into premature infant microbiomes. Video abstract.

Identification of Diagnostic Markers for Major Depressive Disorder Using Machine Learning Methods

Background

Major depressive disorder (MDD) is a global health challenge that impacts the quality of patients’ lives severely. The disorder can manifest in many forms with different combinations of symptoms, which makes its clinical diagnosis difficult. Robust biomarkers are greatly needed to improve diagnosis and to understand the etiology of the disease. The main purpose of this study was to create a predictive model for MDD diagnosis based on peripheral blood transcriptomes.

Materials and Methods

We collected nine RNA expression datasets for MDD patients and healthy samples from the Gene Expression Omnibus database. After a series of quality control and heterogeneity tests, 302 samples from six studies were deemed suitable for the study. R package “MetaOmics” was applied for systematic meta-analysis of genome-wide expression data. Receiver operating characteristic (ROC) curve analysis was used to evaluate the diagnostic effectiveness of individual genes. To obtain a better diagnostic model, we also adopted the support vector machine (SVM), random forest (RF), k-nearest neighbors (kNN), and naive Bayesian (NB) tools for modeling, with the RF method being used for feature selection.

Results

Our analysis revealed six differentially expressed genes (AKR1C3, ARG1, KLRB1, MAFG, TPST1, and WWC3) with a false discovery rate (FDR) < 0.05 between MDD patients and control subjects. We then evaluated the diagnostic ability of these genes individually. With single gene prediction, we achieved a corresponding area under the curve (AUC) value of 0.63 ± 0.04, 0.67 ± 0.07, 0.70 ± 0.11, 0.64 ± 0.08, 0.68 ± 0.07, and 0.62 ± 0.09, respectively, for these genes. Next, we constructed the classifiers of SVM, RF, kNN, and NB with an AUC of 0.84 ± 0.09, 0.81 ± 0.10, 0.73 ± 0.11, and 0.83 ± 0.09, respectively, in validation datasets, suggesting that the SVM classifier might be superior for constructing an MDD diagnostic model. The final SVM classifier including 70 feature genes was capable of distinguishing MDD samples from healthy controls and yielded an AUC of 0.78 in an independent dataset.

Conclusion

This study provides new insights into potential biomarkers through meta-analysis of GEO data. Constructing different machine learning models based on these biomarkers could be a valuable approach for diagnosing MDD in clinical practice.

Identification of Dysregulated Complement Activation Pathways Driven by N-Glycosylation Alterations in T2D Patients

Diabetes has become a major public health concern worldwide, most of which are type 2 diabetes (T2D). The diagnosis of T2D is commonly based on plasma glucose levels, and there are no reliable clinical biomarkers available for early detection. Recent advances in proteome technologies offer new opportunity for the understanding of T2D; however, the underlying proteomic characteristics of T2D have not been thoroughly investigated yet. Here, using proteomic and glycoproteomic profiling, we provided a comprehensive landscape of molecular alterations in the fasting plasma of the 24 Chinese participants, including eight T2D patients, eight prediabetic (PDB) subjects, and eight healthy control (HC) individuals. Our analyses identified a diverse set of potential biomarkers that might enhance the efficiency and accuracy based on current existing biological indicators of (pre)diabetes. Through integrative omics analysis, we showed the capability of glycoproteomics as a complement to proteomics or metabolomics, to provide additional insights into the pathogenesis of (pre)diabetes. We have newly identified systemic site-specific N-glycosylation alterations underlying T2D patients in the complement activation pathways, including decreased levels of N-glycopeptides from C1s, MASP1, and CFP proteins, and increased levels of N-glycopeptides from C2, C4, C4BPA, C4BPB, and CFH. These alterations were not observed at proteomic levels, suggesting new opportunities for the diagnosis and treatment of this disease. Our results demonstrate a great potential role of glycoproteomics in understanding (pre)diabetes and present a new direction for diabetes research which deserves more attention.

Clinicopathological-Associated Regulatory Network of Deregulated circRNAs in Hepatocellular Carcinoma

Simple Summary

Here, we present a novel strategy to identify key signatures of clinically-relevant co-expressed circRNA-mRNA networks in pertinent cancer-pathways that modulate the prognosis of HCC patients, by integrating clinicopathological features, circRNA and mRNA expression profiles. Five master circRNAs were identified and experimentally demonstrated to upregulate proliferate and promote transformation. Through further integration with miRNA-expression profiles, clinically-relevant competing-endogenous-RNA (ceRNA) networks of circRNA-miRNA-mRNAs were constructed. The most up-regulated nodal-circRNA, circGPC3 was experimentally demonstrated to up-regulate cell-cycle, migration and invasion. circGPC3 was found to act as a sponge of miR-378a-3p to regulate ASPM expression and modulate cell transformation. These 5 nodal circRNAs has potential to be good prognostic biomarkers with good prognostic performance. circGPC3 has great potential to be a promising non-invasive prognostic biomarker for early HCC. We have thus demonstrated the robustness of bioinformatically-predicted master circRNAs in clinically-relevant, circRNA-mRNA networks, underscoring the important roles that these identified deregulated key/master circRNAs play in HCC.

Abstract

Hepatocellular carcinoma (HCC) is one of the most common and lethal cancers worldwide. Here, we present a novel strategy to identify key circRNA signatures of clinically relevant co-expressed circRNA-mRNA networks in pertinent cancer-pathways that modulate prognosis of HCC patients, by integrating clinic-pathological features, circRNA and mRNA expression profiles. Through further integration with miRNA expression profiles, clinically relevant competing-endogenous-RNA (ceRNA) networks of circRNA-miRNA-mRNAs were constructed. At least five clinically relevant nodal-circRNAs, co-expressed with numerous genes, were identified from the circRNA-mRNA networks. These nodal circRNAs upregulated proliferation (except circRaly) and transformation in cells. The most upregulated nodal-circRNA, circGPC3, associated with higher-grade tumors and co-expressed with 33 genes, competes with 11 mRNAs for two shared miRNAs. circGPC3 was experimentally demonstrated to upregulate cell-cycle and migration/invasion in both transformed and non-transformed liver cell-lines. circGPC3 was further shown to act as a sponge of miR-378a-3p to regulate APSM (Abnormal spindle-like microcephaly associated) expression and modulate cell transformation. This study identifies 5 key nodal master circRNAs in a clinically relevant circRNA-centric network that are significantly associated with poorer prognosis of HCC patients and promotes tumorigenesis in cell-lines. The identification and characterization of these key circRNAs in clinically relevant circRNA-mRNA and ceRNA networks may facilitate the design of novel strategies targeting these important regulators for better HCC prognosis.

CladeOScope: functional interactions through the prism of clade-wise co-evolution

Mapping co-evolved genes via phylogenetic profiling (PP) is a powerful approach to uncover functional interactions between genes and to associate them with pathways. Despite many successful endeavors, the understanding of co-evolutionary signals in eukaryotes remains partial. Our hypothesis is that 'Clades', branches of the tree of life (e.g. primates and mammals), encompass signals that cannot be detected by PP using all eukaryotes. As such, integrating information from different clades should reveal local co-evolution signals and improve function prediction. Accordingly, we analyzed 1028 genomes in 66 clades and demonstrated that the co-evolutionary signal was scattered across clades. We showed that functionally related genes are frequently co-evolved in only parts of the eukaryotic tree and that clades are complementary in detecting functional interactions within pathways. We examined the non-homologous end joining pathway and the UFM1 ubiquitin-like protein pathway and showed that both demonstrated distinguished co-evolution patterns in specific clades. Our research offers a different way to look at co-evolution across eukaryotes and points to the importance of modular co-evolution analysis. We developed the 'CladeOScope' PP method to integrate information from 16 clades across over 1000 eukaryotic genomes and is accessible via an easy to use web server at http://cladeoscope.cs.huji.ac.il.

Whole-Exome Sequencing Analysis of Oral Squamous Cell Carcinoma Delineated by Tobacco Usage Habits

Oral squamous cell carcinoma (OSCC) is a common cancer of the oral cavity in India. Cigarette smoking and chewing tobacco are known risk factors associated with OSCC. However, genomic alterations in OSCC with varied tobacco consumption history are not well-characterized. In this study, we carried out whole-exome sequencing to characterize the mutational landscape of OSCC tumors from subjects with different tobacco consumption habits. We identified several frequently mutated genes, including TP53, NOTCH1, CASP8, RYR2, LRP2, CDKN2A, and ATM. TP53 and HRAS exhibited mutually exclusive mutation patterns. We identified recurrent amplifications in the 1q31, 7q35, 14q11, 22q11, and 22q13 regions and observed amplification of EGFR in 25% of samples with tobacco consumption history. We observed genomic alterations in several genes associated with PTK6 signaling. We observed alterations in clinically actionable targets including ERBB4, HRAS, EGFR, NOTCH1, NOTCH4, and NOTCH3. We observed enrichment of signature 29 in 40% of OSCC samples from tobacco chewers. Signature 15 associated with defective DNA mismatch repair was enriched in 80% of OSCC samples. NOTCH1 was mutated in 36% of samples and harbored truncating as well as missense variants. We observed copy number alterations in 67% of OSCC samples. Several genes associated with non-receptor tyrosine kinase signaling were affected in OSCC. These molecules can serve as potential candidates for therapeutic targeting in OSCC.

Viral status, immune microenvironment and immunological response to checkpoint inhibitors in hepatocellular carcinoma

Background and aims

Immune checkpoint inhibitors (ICIs) targeting the programmed cell death protein 1 (PD-1)/programmed death-ligand 1 (PD-L1) pathway have clinical activity in hepatocellular carcinoma (HCC), but only a subset of patients respond to these therapies, highlighting a need for novel biomarkers to improve clinical benefit. HCC usually occurs in the setting of liver cirrhosis from chronic hepatitis B or C viral infection, but the effects of viral status on the tumor immune microenvironment and clinical responses to ICIs in HCC remains unclear.

Methods

We conducted a meta-analysis to estimate the objective response rates for PD-1/PD-L1 inhibitors in virally-infected and uninfected patients, and examined the effects of viral etiology on the tumor microenvironment using data from The Cancer Genome Atlas, as well as peripheral blood responses using an independent cohort of patients studied by mass cytometry (cytometry by time-of-flight (CyTOF)).

Results

Meta-analysis comparing objective response rates (ORR) between virally-infected and uninfected patients showed no clinically meaningful difference (absolute difference of ORR in virally-infected vs uninfected=−1.4%, 95% CI: −13.5% to 10.6%). There was no relationship between viral etiology on features of the tumor immune microenvironment that are known to modulate responses to PD-1/PD-L1 inhibitors, and the tumor mutational burden was similar between virally-infected and uninfected HCC. RNA sequencing of tissue-resident T cell and B cell repertoires similarly showed no effect of viral status on their diversity. CyTOF analysis of peripheral blood specimens further demonstrated similar expression of immune-related markers in response to PD-1 inhibitor therapy in virally-infected and uninfected HCC.

Conclusion

There is no significant effect of viral etiology on the tumor immune microenvironment in HCC, and viral status should not be used as a criterion to select patients for PD-1/PD-L1 therapy.

Proteomic Analyses of Vitreous in Proliferative Diabetic Retinopathy: Prior Studies and Future Outlook

Vitreous fluid is becoming an increasingly popular medium for the study of retinal disease. Numerous studies have demonstrated that proteomic analysis of the vitreous from patients with proliferative diabetic retinopathy yields valuable molecular information regarding known and novel proteins and pathways involved in this disease. However, there is no standardized methodology for vitreous proteomic studies. Here, we share a suggested protocol for such studies and outline the various experimental and analytic methods that are currently available. We also review prior mass spectrometry-based proteomic studies of the vitreous from patients with proliferative diabetic retinopathy, discuss common pitfalls of these studies, and propose next steps for moving the field forward.

Understanding the Molecular Drivers of Disease Heterogeneity in Crohn's Disease Using Multi-omic Data Integration and Network Analysis

Crohn's disease (CD), a form of inflammatory bowel disease (IBD), is characterized by heterogeneity along multiple clinical axes, which in turn impacts disease progression and treatment modalities. Using advanced data integration approaches and systems biology tools, we studied the contribution of CD susceptibility variants and gene expression in distinct peripheral immune cell subsets (CD14+ monocytes and CD4+ T cells) to relevant clinical traits. Our analyses revealed that most clinical traits capturing CD heterogeneity could be associated with CD14+ and CD4+ gene expression rather than disease susceptibility variants. By disentangling the sources of variation, we identified molecular features that could potentially be driving the heterogeneity of various clinical traits of CD patients. Further downstream analyses identified contextual hub proteins such as genes encoding barrier functions, antimicrobial peptides, chemokines, and their receptors, which are either targeted by drugs used in CD or other inflammatory diseases or are relevant to the biological functions implicated in disease pathology. These hubs could be used as cell type-specific targets to treat specific subtypes of CD patients in a more individualized approach based on the underlying biology driving their disease subtypes. Our study highlights the importance of data integration and systems approaches to investigate complex and heterogeneous diseases such as IBD.

Drug-Drug Interaction between Metformin and Sorafenib Alters Antitumor Effect in Hepatocellular Carcinoma Cells

Hepatocellular carcinoma (HCC) is the most common primary liver malignancy and is one of the leading causes of cancer-related deaths worldwide. The multitarget inhibitor sorafenib is a first-line treatment of patients with advanced unresectable HCC. Recent clinical studies have evidenced that patients treated with sorafenib together with the antidiabetic drug metformin have a survival disadvantage compared with patients receiving sorafenib only. Here, we examined whether a clinically relevant dose of metformin (50 mg/kg per day) could influence the antitumoral effects of sorafenib (15 mg/kg per day) in a subcutaneous xenograft model of human HCC growth using two different sequences of administration, i.e., concomitant versus sequential dosing regimens. We observed that the administration of metformin 6 hours prior to sorafenib was significantly less effective in inhibiting tumor growth (15.4% tumor growth inhibition) than concomitant administration of the two drugs (59.5% tumor growth inhibition). In vitro experiments confirmed that pretreatment of different human HCC cell lines with metformin reduced the effects of sorafenib on cell viability, proliferation, and signaling. Transcriptomic analysis confirmed significant differences between xenografted tumors obtained under the concomitant and the sequential dosing regimens. Taken together, these observations call into question the benefit of parallel use of metformin and sorafenib in patients with advanced HCC and diabetes, as the interaction between the two drugs could ultimately compromise patient survival. SIGNIFICANCE STATEMENT: When drugs are administered sequentially, metformin alters the antitumor effect of sorafenib, the reference treatment for advanced hepatocellular carcinoma, in a preclinical murine xenograft model of liver cancer progression as well as in hepatic cancer cell lines. Defective activation of the AMP-activated protein kinase pathway as well as major transcriptomic changes are associated with the loss of the antitumor effect. These results echo recent clinical work reporting a poorer prognosis for patients with liver cancer who were cotreated with metformin and sorafenib.

Towards defining the Mechanisms of Alzheimer's disease based on a contextual analysis of molecular pathways

Alzheimer's disease (AD) is posing an increasingly profound problem to society. Our genuine understanding of the pathogenesis of AD is inadequate and as a consequence, diagnostic and therapeutic strategies are currently insufficient. The understandable focus of many studies is the identification of molecules with high diagnostic utility however the opportunity to obtain a further understanding of the mechanistic origins of the disease from such putative biomarkers is often overlooked. This study examines the involvement of biomarkers in AD to shed light on potential mechanisms and pathways through which they are implicated in the pathology of this devastating neurodegenerative disorder. The computational tools required to analyse ever-growing datasets in the context of AD are also discussed.

An Immunometabolic Shift Modulates Cytotoxic Lymphocyte Activation During Melanoma Progression in TRPA1 Channel Null Mice

Melanoma skin cancer is extremely aggressive with increasing incidence and mortality. Among the emerging therapeutic targets in the treatment of cancer, the family of transient receptor potential channels (TRPs) has been reported as a possible pharmacological target. Specifically, the ankyrin subfamily, representing TRPA1 channels, can act as a pro-inflammatory hub. These channels have already been implicated in the control of intracellular metabolism in several cell models, but little is known about their role in immune cells, and how it could affect tumor progression in a process known as immune surveillance. Here, we investigated the participation of the TRPA1 channel in the immune response against melanoma tumor progression in a mouse model. Using Trpa1 ^+/+ and Trpa1 ^-/- animals, we evaluated tumor progression using murine B16-F10 cells and assessed isolated CD8+ T cells for respiratory and cytotoxic functions. Tumor growth was significantly reduced in Trpa1 ^-/- animals. We observed an increase in the frequency of circulating lymphocytes. Using a dataset of CD8+ T cells isolated from metastatic melanoma patients, we found that TRPA1 reduction correlates with several immunological pathways. Naïve CD8+ T cells from Trpa1 ^+/+ and Trpa1 ^-/- animals showed different mitochondrial respiration and glycolysis profiles. However, under CD3/CD28 costimulatory conditions, the absence of TRPA1 led to an even more extensive metabolic shift, probably linked to a greater in vitro killling ability of Trpa1 ^-/- CD8+ T cells. Therefore, these data demonstrate an unprecedented role of TRPA1 channel in the metabolism control of the immune system cells during carcinogenesis.

Performance assessment of sample-specific network control methods for bulk and single-cell biological data analysis

In the past few years, a wealth of sample-specific network construction methods and structural network control methods has been proposed to identify sample-specific driver nodes for supporting the Sample-Specific network Control (SSC) analysis of biological networked systems. However, there is no comprehensive evaluation for these state-of-the-art methods. Here, we conducted a performance assessment for 16 SSC analysis workflows by using the combination of 4 sample-specific network reconstruction methods and 4 representative structural control methods. This study includes simulation evaluation of representative biological networks, personalized driver genes prioritization on multiple cancer bulk expression datasets with matched patient samples from TCGA, and cell marker genes and key time point identification related to cell differentiation on single-cell RNA-seq datasets. By widely comparing analysis of existing SSC analysis workflows, we provided the following recommendations and banchmarking workflows. (i) The performance of a network control method is strongly dependent on the up-stream sample-specific network method, and Cell-Specific Network construction (CSN) method and Single-Sample Network (SSN) method are the preferred sample-specific network construction methods. (ii) After constructing the sample-specific networks, the undirected network-based control methods are more effective than the directed network-based control methods. In addition, these data and evaluation pipeline are freely available on https://github.com/WilfongGuo/Benchmark_control.

Targeting the tumor stroma: integrative analysis reveal GATA2 and TORYAIP1 as novel prognostic targets in breast and ovarian cancer

Tumor stroma interaction is known to take a crucial role in cancer growth and progression. In the present study, it was performed gene expression analysis of stroma samples with ovarian and breast cancer through an integrative analysis framework to identify common critical biomolecules at multiomics levels. Gene expression datasets were statistically analyzed to identify common differentially expressed genes (DEGs) by comparing tumor stroma and normal stroma samples. The integrative analyses of DEGs indicated that there were 59 common core genes, which might be feasible to be potential marks for cancer stroma targeted strategies. Reporter molecules (i.e. receptor, transcription factors and miRNAs) were determined through a statistical test employing the hypergeometric probability density function. Afterward, the tumor microenvironment protein-protein interaction and the generic network were reconstructed by using identified reporter molecules and common core DEGs. Through a systems medicine approach, it was determined that hub biomolecules, AR, GATA2, miR-124, TOR1AIP1, ESR1, EGFR, STAT1, miR-192, GATA3, COL1A1, in tumor microenvironment generic network. These molecules were also identified as prognostic signatures in breast and ovarian tumor samples via survival analysis. According to literature searching, GATA2 and TORYAIP1 might represent potential biomarkers and candidate drug targets for the stroma targeted cancer therapy applications.

A Population-Specific Major Allele Reference Genome From The United Arab Emirates Population

The ethnic composition of the population of a country contributes to the uniqueness of each national DNA sequencing project and, ideally, individual reference genomes are required to reduce the confounding nature of ethnic bias. This work represents a representative Whole Genome Sequencing effort of an understudied population. Specifically, high coverage consensus sequences from 120 whole genomes and 33 whole exomes were used to construct the first ever population specific major allele reference genome for the United Arab Emirates (UAE). When this was applied and compared to the archetype hg19 reference, assembly of local Emirati genomes was reduced by ∼19% (i.e., some 1 million fewer calls). In compiling the United Arab Emirates Reference Genome (UAERG), sets of annotated 23,038,090 short (novel: 1,790,171) and 137,713 structural (novel: 8,462) variants; their allele frequencies (AFs) and distribution across the genome were identified. Population-specific genetic characteristics including loss-of-function variants, admixture, and ancestral haplogroup distribution were identified and reported here. We also detect a strong correlation between F ST and admixture components in the UAE. This baseline study was conceived to establish a high-quality reference genome and a genetic variations resource to enable the development of regional population specific initiatives and thus inform the application of population studies and precision medicine in the UAE.

RNA-Seq Reveals Function of Bta-miR-149-5p in the Regulation of Bovine Adipocyte Differentiation

Intramuscular fat is a real challenge for the experts of animal science to improve meat quality traits. Research on the mechanism of adipogenesis provides invaluable information for the improvement of meat quality traits. This study investigated the effect of bta-miR-149-5p and its underlying mechanism on lipid metabolism in bovine adipocytes. Bovine adipocytes were differentiated and transfected with bta-miR-149-5p mimics or its negative control (NC). A total of 115 DEGs including 72 upregulated and 43 downregulated genes were identified in bovine adipocytes. The unigenes and GO term biological processes were the most annotated unigene contributor parts at 80.08%, followed by cellular component at 13.4% and molecular function at 6.7%. The KEGG pathways regulated by the DEGs were PI3K-Akt signaling pathway, calcium signaling pathway, pathways in cancer, MAPK signaling pathway, lipid metabolism/metabolic pathway, PPAR signaling pathway, AMPK signaling pathway, TGF-beta signaling pathway, cAMP signaling pathway, cholesterol metabolism, Wnt signaling pathway, and FoxO signaling pathway. In addition to this, the most important reactome enrichment pathways were R-BTA-373813 receptor CXCR2 binding ligands CXCL1 to 7, R-BTA-373791 receptor CXCR1 binding CXCL6 and CXCL8 ligands, R-BTA-210991 basigin interactions, R-BTA-380108 chemokine receptors binding chemokines, R-BTA-445704 calcium binding caldesmon, and R-BTA-5669034 TNFs binding their physiological receptors. Furthermore, the expression trend of the DEGs in these pathways were also exploited. Moreover, the bta-miR-149-5p significantly (p < 0.01) downregulated the mRNA levels of adipogenic marker genes such as CCND2, KLF6, ACSL1, Cdk2, SCD, SIK2, and ZEB1 in bovine adipocytes. In conclusion, our results suggest that bta-miR-149-5p regulates lipid metabolism in bovine adipocytes. The results of this study provide a basis for studying the function and molecular mechanism of the bta-miR-149-5p in regulating bovine adipogenesis.

Network controllability-based algorithm to target personalized driver genes for discovering combinatorial drugs of individual patients

Multiple driver genes in individual patient samples may cause resistance to individual drugs in precision medicine. However, current computational methods have not studied how to fill the gap between personalized driver gene identification and combinatorial drug discovery for individual patients. Here, we developed a novel structural network controllability-based personalized driver genes and combinatorial drug identification algorithm (CPGD), aiming to identify combinatorial drugs for an individual patient by targeting personalized driver genes from network controllability perspective. On two benchmark disease datasets (i.e. breast cancer and lung cancer datasets), performance of CPGD is superior to that of other state-of-the-art driver gene-focus methods in terms of discovery rate among prior-known clinical efficacious combinatorial drugs. Especially on breast cancer dataset, CPGD evaluated synergistic effect of pairwise drug combinations by measuring synergistic effect of their corresponding personalized driver gene modules, which are affected by a given targeting personalized driver gene set of drugs. The results showed that CPGD performs better than existing synergistic combinatorial strategies in identifying clinical efficacious paired combinatorial drugs. Furthermore, CPGD enhanced cancer subtyping by computationally providing personalized side effect signatures for individual patients. In addition, CPGD identified 90 drug combinations candidates from SARS-COV2 dataset as potential drug repurposing candidates for recently spreading COVID-19.

Epigenetic alteration contributes to the transcriptional reprogramming in T-cell prolymphocytic leukemia

T cell prolymphocytic leukemia (T-PLL) is a rare disease with aggressive clinical course. Cytogenetic analysis, whole-exome and whole-genome sequencing have identified primary structural alterations in T-PLL, including inversion, translocation and copy number variation. Recurrent somatic mutations were also identified in genes encoding chromatin regulators and those in the JAK-STAT signaling pathway. Epigenetic alterations are the hallmark of many cancers. However, genome-wide epigenomic profiles have not been reported in T-PLL, limiting the mechanistic study of its carcinogenesis. We hypothesize epigenetic mechanisms also play a key role in T-PLL pathogenesis. To systematically test this hypothesis, we generated genome-wide maps of regulatory regions using H3K4me3 and H3K27ac ChIP-seq, as well as RNA-seq data in both T-PLL patients and healthy individuals. We found that genes down-regulated in T-PLL are mainly associated with defense response, immune system or adaptive immune response, while up-regulated genes are enriched in developmental process, as well as WNT signaling pathway with crucial roles in cell fate decision. In particular, our analysis revealed a global alteration of regulatory landscape in T-PLL, with differential peaks highly enriched for binding motifs of immune related transcription factors, supporting the epigenetic regulation of oncogenes and genes involved in DNA damage response and T-cell activation. Together, our work reveals a causal role of epigenetic dysregulation in T-PLL.

Microbial Metabolomics: From Methods to Translational Applications

Most microbe-associated infectious diseases severely affect human health. However, clinical diagnosis of pathogenic diseases remains challenging due to the lack of specific and highly reliable methods. To better understand the diagnosis, pathogenesis, and treatment of these diseases, systems biology-driven metabolomics goes beyond the annotated phenotype and better targets the functions than conventional approaches. As a novel strategy for analysis of metabolomes in microbes, microbial metabolomics has been recently used to study many diseases, such as obesity, urinary tract infection (UTI), and hepatitis C. In this chapter, we attempt to introduce various microbial metabolomics methods to better interpret the microbial metabolism underlying a diversity of infectious diseases and inspire scientists to pay more attention to microbial metabolomics, enabling broadly and efficiently its translational applications to infectious diseases, from molecular diagnosis to therapeutic discovery.

Comprehensive functional annotation of susceptibility variants identifies genetic heterogeneity between lung adenocarcinoma and squamous cell carcinoma

Although genome-wide association studies have identified more than eighty genetic variants associated with non-small cell lung cancer (NSCLC) risk, biological mechanisms of these variants remain largely unknown. By integrating a large-scale genotype data of 15 581 lung adenocarcinoma (AD) cases, 8350 squamous cell carcinoma (SqCC) cases, and 27 355 controls, as well as multiple transcriptome and epigenomic databases, we conducted histology-specific meta-analyses and functional annotations of both reported and novel susceptibility variants. We identified 3064 credible risk variants for NSCLC, which were overrepresented in enhancer-like and promoter-like histone modification peaks as well as DNase I hypersensitive sites. Transcription factor enrichment analysis revealed that USF1 was AD-specific while CREB1 was SqCC-specific. Functional annotation and gene-based analysis implicated 894 target genes, including 274 specifics for AD and 123 for SqCC, which were overrepresented in somatic driver genes (ER = 1.95, P = 0.005). Pathway enrichment analysis and Gene-Set Enrichment Analysis revealed that AD genes were primarily involved in immune-related pathways, while SqCC genes were homologous recombination deficiency related. Our results illustrate the molecular basis of both well-studied and new susceptibility loci of NSCLC, providing not only novel insights into the genetic heterogeneity between AD and SqCC but also a set of plausible gene targets for post-GWAS functional experiments.

A genome-wide association study of severe asthma exacerbations in Latino children and adolescents

Severe asthma exacerbations are a major cause of school absences and healthcare costs in children, particularly those in high-risk racial/ethnic groups.To identify susceptibility genes for severe asthma exacerbations in Latino children and adolescents, we conducted a meta-analysis of genome-wide association studies (GWAS) in 4010 Latino youth with asthma in four independent cohorts, including 1693 Puerto Ricans, 1019 Costa Ricans, 640 Mexicans, 256 Brazilians and 402 members of other Latino subgroups. We then conducted methylation quantitative trait locus, expression quantitative trait locus and expression quantitative trait methylation analyses to assess whether the top single nucleotide polymorphism (SNP) in the meta-analysis is linked to DNA methylation and gene expression in nasal (airway) epithelium in separate cohorts of Puerto Rican and Dutch children and adolescents.In the meta-analysis of GWAS, an SNP in FLJ22447 (rs2253681) was significantly associated with 1.55 increased odds of severe asthma exacerbation (95% CI 1.34-1.79, p=6.3×10-9). This SNP was significantly associated with DNA methylation of a CpG site (cg25024579) at the FLJ22447 locus, which was in turn associated with increased expression of KCNJ2-AS1 in nasal airway epithelium from Puerto Rican children and adolescents (β=0.10, p=2.18×10-7).SNP rs2253681 was significantly associated with both DNA methylation of a cis-CpG in FLJ22447 and severe asthma exacerbations in Latino youth. This may be partly explained by changes in airway epithelial expression of a gene recently implicated in atopic asthma in Puerto Rican children and adolescents (KCNJ2-AS1).

Highly deregulated lncRNA LOC is associated with overall worse prognosis in Hepatocellular Carcinoma patients

Although numerous long non-coding RNAs (lncRNAs) were reported to be deregulated in Hepatocellular Carcinoma (HCC), experimentally characterized, and/or associated with patient's clinical characteristics, there is, thus far, minimal concerted research strategy to identify deregulated lncRNAs that modulate prognosis of HCC patients. Here, we present a novel strategy where we identify lncRNAs, which are not only de-regulated in HCC patients, but are also associated with pertinent clinical characteristics, potentially contributing to the prognosis of HCC patients. LOC101926913 (LOC) was further characterized because it is the most highly differentially expressed amongst those that are associated with the most number of clinical features (tumor-stage, vascular and tumor invasion and poorer overall survival). Experimental gain- and loss-of-function manipulation of LOC in liver cell-lines highlight LOC as a potential onco-lncRNA promoting cell proliferation, anchorage independent growth and invasion. LOC expression in cells up-regulated genes involved in GTPase-activities and downregulated genes associated with cellular detoxification, oxygen- and drug-transport. Hence, LOC may represent a novel therapeutic target, modulating prognosis of HCC patients through up-regulating GTPase-activities and down-regulating detoxification, oxygen- and drug-transport. This strategy may thus be useful for the identification of clinically relevant lncRNAs as potential biomarkers/targets that modulate prognosis in other cancers as well.

Destination Amyotrophic Lateral Sclerosis

Amyotrophic Lateral Sclerosis (ALS) is a prototypical neurodegenerative disease characterized by progressive degeneration of motor neurons both in the brain and spinal cord. The constantly evolving nature of ALS represents a fundamental dimension of individual differences that underlie this disorder, yet it involves multiple levels of functional entities that alternate in different directions and finally converge functionally to define ALS disease progression. ALS may start from a single entity and gradually becomes multifactorial. However, the functional convergence of these diverse entities in eventually defining ALS progression is poorly understood. Various hypotheses have been proposed without any consensus between the for-and-against schools of thought. The present review aims to capture explanatory hierarchy both in terms of hypotheses and mechanisms to provide better insights on how they functionally connect. We can then integrate them within a common functional frame of reference for a better understanding of ALS and defining future treatments and possible therapeutic strategies. Here, we provide a philosophical understanding of how early leads are crucial to understanding the endpoints in ALS, because invariably, all early symptomatic leads are underpinned by neurodegeneration at the cellular, molecular and genomic levels. Consolidation of these ideas could be applied to other neurodegenerative diseases (NDs) and guide further critical thinking to unveil their roadmap of destination ALS.

Machine Learning Reduced Gene/Non-Coding RNA Features That Classify Schizophrenia Patients Accurately and Highlight Insightful Gene Clusters

RNA-seq has been a powerful method to detect the differentially expressed genes/long non-coding RNAs (lncRNAs) in schizophrenia (SCZ) patients; however, due to overfitting problems differentially expressed targets (DETs) cannot be used properly as biomarkers. This study used machine learning to reduce gene/non-coding RNA features. Dorsolateral prefrontal cortex (dlpfc) RNA-seq data from 254 individuals was obtained from the CommonMind consortium. The average predictive accuracy for SCZ patients was 67% based on coding genes, and 96% based on long non-coding RNAs (lncRNAs). Machine learning is a powerful algorithm to reduce functional biomarkers in SCZ patients. The lncRNAs capture the characteristics of SCZ tissue more accurately than mRNA as the former regulate every level of gene expression, not limited to mRNA levels.

Identification of IFN-Induced Transmembrane Protein 1 With Prognostic Value in Pancreatic Cancer Using Network Module-Based Analysis

Despite improvements reported in diagnosis and treatments in recent decades, pancreatic cancer is still characterized by poor prognosis and low survival rate among solid tumors. Intensive interests have grown in exploring novel predictive biomarkers, aiming to enhance the efficiency in early detection and treatment prognosis. In this study, we identified the differentially expressed genes (DEGs) in pancreatic cancer by analyzing five gene expression profiles and established the functional modules according to the functional interaction (FI) network between the DEGs. A significant upregulation of the selected DEG, interferon (IFN)-induced transmembrane protein 1 (IFITM1), was evaluated in several bioinformatics online tools and verified with immunohistochemistry staining from samples of 90 patients with pancreatic cancer. Prognostic data showed that high expression of IFITM1 associated with poor survival, and multivariate Cox regression analysis showed IFITM1 was one of the independent prognostic factors for overall survival. Meanwhile, significant correlations of the expression of IFITM1 and the infiltration of immune cells were found by TIMER. Furthermore, a higher level of IFITM1 was assessed in pancreatic cancer cell lines compared to normal human pancreatic duct epithelial cells, and silencing IFITM1 in tumor cells remarkedly inhibited cancer tumorigenicity. Collectively, our findings suggested that IFITM1 might have promising utility for pancreatic cancer.

Unveiling COVID-19-associated organ-specific cell types and cell-specific pathway cascade

The novel coronavirus or COVID-19 has first been found in Wuhan, China, and became pandemic. Angiotensin-converting enzyme 2 (ACE2) plays a key role in the host cells as a receptor of Spike-I Glycoprotein of COVID-19 which causes final infection. ACE2 is highly expressed in the bladder, ileum, kidney and liver, comparing with ACE2 expression in the lung-specific pulmonary alveolar type II cells. In this study, the single-cell RNAseq data of the five tissues from different humans are curated and cell types with high expressions of ACE2 are identified. Subsequently, the protein-protein interaction networks have been established. From the network, potential biomarkers which can form functional hubs, are selected based on k-means network clustering. It is observed that angiotensin PPAR family proteins show important roles in the functional hubs. To understand the functions of the potential markers, corresponding pathways have been researched thoroughly through the pathway semantic networks. Subsequently, the pathways have been ranked according to their influence and dependency in the network using PageRank algorithm. The outcomes show some important facts in terms of infection. Firstly, renin-angiotensin system and PPAR signaling pathway can play a vital role for enhancing the infection after its intrusion through ACE2. Next, pathway networks consist of few basic metabolic and influential pathways, e.g. insulin resistance. This information corroborate the fact that diabetic patients are more vulnerable to COVID-19 infection. Interestingly, the key regulators of the aforementioned pathways are angiontensin and PPAR family proteins. Hence, angiotensin and PPAR family proteins can be considered as possible therapeutic targets. Contact: sagnik.sen2008@gmail.com, umaulik@cse.jdvu.ac.in Supplementary information: Supplementary data are available online.