Metabolomics and integrated network analysis reveal roles of endocannabinoids and large neutral amino acid balance in the ayahuasca experience.

There has been a renewed interest in the potential use of psychedelics for the treatment of psychiatric conditions. Nevertheless, little is known about the mechanism of action and molecular pathways influenced by ayahuasca use in humans. Therefore, for the first time, our study aims to investigate the human metabolomics signature after consumption of a psychedelic, ayahuasca, and its connection with both the psychedelic-induced subjective effects and the plasma concentrations of ayahuasca alkaloids. Plasma samples of 23 individuals were collected both before and after ayahuasca consumption. Samples were analysed through targeted metabolomics and further integrated with subjective ratings of the ayahuasca experience (i.e., using the 5-Dimension Altered States of Consciousness Rating Scale [ASC]), and plasma ayahuasca-alkaloids using integrated network analysis. Metabolic pathways enrichment analysis using diffusion algorithms for specific KEGG modules was performed on the metabolic output. Compared to baseline, the consumption of ayahuasca increased N-acyl-ethanolamine endocannabinoids, decreased 2-acyl-glycerol endocannabinoids, and altered several large-neutral amino acids (LNAAs). Integrated network results indicated that most of the LNAAs were inversely associated with 9 out of the 11 subscales of the ASC, except for tryptophan which was positively associated. Several endocannabinoids and hexosylceramides were directly associated with the ayahuasca alkaloids. Enrichment analysis confirmed dysregulation in several pathways involved in neurotransmission such as serotonin and dopamine synthesis. In conclusion, a crosstalk between the circulating LNAAs and the subjective effects is suggested, which is independent of the alkaloid concentrations and provides insights into the specific metabolic fingerprint and mechanism of action underlying ayahuasca experiences.

Identification of consistent post-translational regulatory triplets related to oncogenic and tumour suppressive modulators in childhood acute lymphoblastic leukemia.

BACKGROUND: Acute lymphoblastic leukemia (ALL) is the most common type of childhood cancer. It can be caused by mutations that turn on oncogenes or turn off tumour suppressor genes. For instance, changes in certain genes including Rb and p53 are common in ALL cells. Oncogenes and TSGs may serve as a modulator gene to regulate the gene expression level via their respective target genes. To investigate the regulatory relationship between oncogenes, tumour suppressor genes and transcription factors at the post translational level in childhood ALL, we performed an integrative network analysis on the gene regulation in the post-translational level for childhood ALL based on many publicly available cancer gene expression data including TARGET and GEO database. METHODS: We collected 259 childhood ALL-related genes from the latest online leukemia database, Leukemia Gene Literature Database. These 259 genes were selected from a comprehensive systematic literature with experimental evidences. The identified and curated genes were also associated with patient survival cases and we incorporated this pediatric ALL-related gene list into our analysis. We extracted the known human TFs from the TRRUST database. Among 259 childhood ALL-related genes, 101 unique regulators were mapped to the list of oncogene and tumour suppressor genes (TSGs) from the ONGene and the TSGene databases, and these included 74 TSGs, 62 oncogenes and 46 TF genes. RESULTS: The resulted regulation was presented as a hierarchical regulatory network with transcription factors (TFs) as intermediate regulators connecting the top modulators (oncogene and TSGs) to the common target genes. Cross-validation was applied to the results from the TARGET dataset by identifying the consistent regulatory motifs based on three independent ALL expression datasets. A three-layer regulatory network of consistent positive modulators in childhood ALL was constructed in which 74 modulators (40 oncogenes, 34 TSGs) are considered as the most important regulators. The middle layer and the bottom layer contain 34 TFs and 176 target genes, respectively. Oncogenes mostly participated in positive regulation of gene expression and the transcription process of RNA II polymerase, while TSGs were mainly involved in the negative regulation of gene expression. In addition, the oncogene-specific targets were enriched with regulators of the MAPK cascade while tumour suppressor-specific targets were associated with cell death. CONCLUSION: The results revealed that oncogenes and TSGs possess a different functional regulatory pattern with regard to not only their biological functions but also their specific target genes in childhood ALL cancer progression. Taken together, our findings could contribute to a better understanding of the important regulatory mechanisms and this method could be used to analyse the targeted genes at the post-translational level in childhood ALL through integrative network analysis.

Comprehensive Network Analysis Reveals Alternative Splicing-Related lncRNAs in Hepatocellular Carcinoma.

It is increasingly appreciated that long non-coding RNAs (lncRNAs) associated with alternative splicing (AS) could be involved in aggressive hepatocellular carcinoma. Although many recent studies show the alteration of RNA alternative splicing by deregulated lncRNAs in cancer, the extent to which and how lncRNAs impact alternative splicing at the genome scale remains largely elusive. We analyzed RNA-seq data obtained from 369 hepatocellular carcinomas (HCCs) and 160 normal liver tissues, quantified 198,619 isoform transcripts, and identified a total of 1,375 significant AS events in liver cancer. In order to predict novel AS-associated lncRNAs, we performed an integration of co-expression, protein-protein interaction (PPI) and epigenetic interaction networks that links lncRNA modulators (such as splicing factors, transcript factors, and miRNAs) along with their targeted AS genes in HCC. We developed a random walk-based multi-graphic (RWMG) model algorithm that prioritizes functional lncRNAs with their associated AS targets to computationally model the heterogeneous networks in HCC. RWMG shows a good performance evaluated by the ROC curve based on cross-validation and bootstrapping strategies. As a conclusion, our robust network-based framework has derived 31 AS-related lncRNAs that not only validates known cancer-associated cases MALAT1 and HOXA11-AS, but also reveals new players such as DNM1P35 and DLX6-AS1with potential functional implications. Survival analysis further provides insights into the clinical significance of identified lncRNAs.

3D Cellular Architecture Affects MicroRNA and Protein Cargo of Extracellular Vesicles.

The success of malignant tumors is conditioned by the intercellular communication between tumor cells and their microenvironment, with extracellular vesicles (EVs) acting as main mediators. While the value of 3D conditions to study tumor cells is well established, the impact of cellular architecture on EV content and function is not investigated yet. Here, a recently developed 3D cell culture microwell array is adapted for EV production and a comprehensive comparative analysis of biochemical features, RNA and proteomic profiles of EVs secreted by 2D vs 3D cultures of gastric cancer cells, is performed. 3D cultures are significantly more efficient in producing EVs than 2D cultures. Global upregulation of microRNAs and downregulation of proteins in 3D are observed, indicating their dynamic coregulation in response to cellular architecture, with the ADP-ribosylation factor 6 signaling pathway significantly downregulated in 3D EVs. The data strengthen the biological relevance of cellular architecture for production and cargo of EVs.

The Immune Landscape of Cancer.

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-ß dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.