Curation of over 10 000 transcriptomic studies to enable data reuse

Resource: A curated database of brain-related functional gene sets (Brain.GMT)

Transcriptional profiling has become a common tool for investigating the nervous system. During analysis, differential expression results are often compared to functional ontology databases, which contain curated gene sets representing well-studied pathways. This dependence can cause neuroscience studies to be interpreted in terms of functional pathways documented in better studied tissues (e.g., liver) and topics (e.g., cancer), and systematically emphasizes well-studied genes, leaving other findings in the obscurity of the brain "ignorome". To address this issue, we compiled a curated database of 918 gene sets related to nervous system function, tissue, and cell types ("Brain.GMT") that can be used within common analysis pipelines (GSEA, limma, edgeR) to interpret results from three species (rat, mouse, human). Brain.GMT includes brain-related gene sets curated from the Molecular Signatures Database (MSigDB) and extracted from public databases (GeneWeaver, Gemma, DropViz, BrainInABlender, HippoSeq) and published studies containing differential expression results. Although Brain.GMT is still undergoing development and currently only represents a fraction of available brain gene sets, "brain ignorome" genes are already better represented than in traditional Gene Ontology databases. Moreover, Brain.GMT substantially improves the quantity and quality of gene sets identified as enriched with differential expression in neuroscience studies, enhancing interpretation. •We compiled a curated database of 918 gene sets related to nervous system function, tissue, and cell types ("Brain.GMT").•Brain.GMT can be used within common analysis pipelines (GSEA, limma, edgeR) to interpret neuroscience transcriptional profiling results from three species (rat, mouse, human).•Although Brain.GMT is still undergoing development, it substantially improved the interpretation of differential expression results within our initial use cases.

In silico approaches for drug repurposing in oncology: a scoping review

Introduction: Cancer refers to a group of diseases characterized by the uncontrolled growth and spread of abnormal cells in the body. Due to its complexity, it has been hard to find an ideal medicine to treat all cancer types, although there is an urgent need for it. However, the cost of developing a new drug is high and time-consuming. In this sense, drug repurposing (DR) can hasten drug discovery by giving existing drugs new disease indications. Many computational methods have been applied to achieve DR, but just a few have succeeded. Therefore, this review aims to show in silico DR approaches and the gap between these strategies and their ultimate application in oncology. Methods: The scoping review was conducted according to the Arksey and O'Malley framework and the Joanna Briggs Institute recommendations. Relevant studies were identified through electronic searching of PubMed/MEDLINE, Embase, Scopus, and Web of Science databases, as well as the grey literature. We included peer-reviewed research articles involving in silico strategies applied to drug repurposing in oncology, published between 1 January 2003, and 31 December 2021. Results: We identified 238 studies for inclusion in the review. Most studies revealed that the United States, India, China, South Korea, and Italy are top publishers. Regarding cancer types, breast cancer, lymphomas and leukemias, lung, colorectal, and prostate cancer are the top investigated. Additionally, most studies solely used computational methods, and just a few assessed more complex scientific models. Lastly, molecular modeling, which includes molecular docking and molecular dynamics simulations, was the most frequently used method, followed by signature-, Machine Learning-, and network-based strategies. Discussion: DR is a trending opportunity but still demands extensive testing to ensure its safety and efficacy for the new indications. Finally, implementing DR can be challenging due to various factors, including lack of quality data, patient populations, cost, intellectual property issues, market considerations, and regulatory requirements. Despite all the hurdles, DR remains an exciting strategy for identifying new treatments for numerous diseases, including cancer types, and giving patients faster access to new medications.

OMD Curation Toolkit: a workflow for in-house curation of public omics datasets

BACKGROUND

Major advances in sequencing technologies and the sharing of data and metadata in science have resulted in a wealth of publicly available datasets. However, working with and especially curating public omics datasets remains challenging despite these efforts. While a growing number of initiatives aim to re-use previous results, these present limitations that often lead to the need for further in-house curation and processing.

RESULTS

Here, we present the Omics Dataset Curation Toolkit (OMD Curation Toolkit), a python3 package designed to accompany and guide the researcher during the curation process of metadata and fastq files of public omics datasets. This workflow provides a standardized framework with multiple capabilities (collection, control check, treatment and integration) to facilitate the arduous task of curating public sequencing data projects. While centered on the European Nucleotide Archive (ENA), the majority of the provided tools are generic and can be used to curate datasets from different sources.

CONCLUSIONS

Thus, it offers valuable tools for the in-house curation previously needed to re-use public omics data. Due to its workflow structure and capabilities, it can be easily used and benefit investigators in developing novel omics meta-analyses based on sequencing data.

Resource: A Curated Database of Brain-Related Functional Gene Sets (Brain.GMT)

Transcriptional profiling has become a common tool for investigating the nervous system. During analysis, differential expression results are often compared to functional ontology databases, which contain curated gene sets representing well-studied pathways. This dependence can cause neuroscience studies to be interpreted in terms of functional pathways documented in better studied tissues (e.g., liver) and topics (e.g., cancer), and systematically emphasizes well-studied genes, leaving other findings in the obscurity of the brain "ignorome". To address this issue, we compiled a curated database of 918 gene sets related to nervous system function, tissue, and cell types ("Brain.GMT") that can be used within common analysis pipelines (GSEA, limma, edgeR) to interpret results from three species (rat, mouse, human). Brain.GMT includes brain-related gene sets curated from the Molecular Signatures Database (MSigDB) and extracted from public databases (GeneWeaver, Gemma, DropViz, BrainInABlender, HippoSeq) and published studies containing differential expression results. Although Brain.GMT is still undergoing development and currently only represents a fraction of available brain gene sets, "brain ignorome" genes are already better represented than in traditional Gene Ontology databases. Moreover, Brain.GMT substantially improves the quantity and quality of gene sets identified as enriched with differential expression in neuroscience studies, enhancing interpretation.

TidyGEO: preparing analysis-ready datasets from Gene Expression Omnibus

TidyGEO is a Web-based tool for downloading, tidying, and reformatting data series from Gene Expression Omnibus (GEO). As a freely accessible repository with data from over 6 million biological samples across more than 4000 organisms, GEO provides diverse opportunities for secondary research. Although scientists may find assay data relevant to a given research question, most analyses require sample-level annotations. In GEO, such annotations are stored alongside assay data in delimited, text-based files. However, the structure and semantics of the annotations vary widely from one series to another, and many annotations are not useful for analysis purposes. Thus, every GEO series must be tidied before it is analyzed. Manual approaches may be used, but these are error prone and take time away from other research tasks. Custom computer scripts can be written, but many scientists lack the computational expertise to create such scripts. To address these challenges, we created TidyGEO, which supports essential data-cleaning tasks for sample-level annotations, such as selecting informative columns, renaming columns, splitting or merging columns, standardizing data values, and filtering samples. Additionally, users can integrate annotations with assay data, restructure assay data, and generate code that enables others to reproduce these steps.

Benchmarking enrichment analysis methods with the disease pathway network

Enrichment analysis (EA) is a common approach to gain functional insights from genome-scale experiments. As a consequence, a large number of EA methods have been developed, yet it is unclear from previous studies which method is the best for a given dataset. The main issues with previous benchmarks include the complexity of correctly assigning true pathways to a test dataset, and lack of generality of the evaluation metrics, for which the rank of a single target pathway is commonly used. We here provide a generalized EA benchmark and apply it to the most widely used EA methods, representing all four categories of current approaches. The benchmark employs a new set of 82 curated gene expression datasets from DNA microarray and RNA-Seq experiments for 26 diseases, of which only 13 are cancers. In order to address the shortcomings of the single target pathway approach and to enhance the sensitivity evaluation, we present the Disease Pathway Network, in which related Kyoto Encyclopedia of Genes and Genomes pathways are linked. We introduce a novel approach to evaluate pathway EA by combining sensitivity and specificity to provide a balanced evaluation of EA methods. This approach identifies Network Enrichment Analysis methods as the overall top performers compared with overlap-based methods. By using randomized gene expression datasets, we explore the null hypothesis bias of each method, revealing that most of them produce skewed P-values.

Speos: an ensemble graph representation learning framework to predict core gene candidates for complex diseases

Understanding phenotype-to-genotype relationships is a grand challenge of 21st century biology with translational implications. The recently proposed "omnigenic" model postulates that effects of genetic variation on traits are mediated by core-genes and -proteins whose activities mechanistically influence the phenotype, whereas peripheral genes encode a regulatory network that indirectly affects phenotypes via core gene products. Here, we develop a positive-unlabeled graph representation-learning ensemble-approach based on a nested cross-validation to predict core-like genes for diverse diseases using Mendelian disorder genes for training. Employing mouse knockout phenotypes for external validations, we demonstrate that core-like genes display several key properties of core genes: Mouse knockouts of genes corresponding to our most confident predictions give rise to relevant mouse phenotypes at rates on par with the Mendelian disorder genes, and all candidates exhibit core gene properties like transcriptional deregulation in disease and loss-of-function intolerance. Moreover, as predicted for core genes, our candidates are enriched for drug targets and druggable proteins. In contrast to Mendelian disorder genes the new core-like genes are enriched for druggable yet untargeted gene products, which are therefore attractive targets for drug development. Interpretation of the underlying deep learning model suggests plausible explanations for our core gene predictions in form of molecular mechanisms and physical interactions. Our results demonstrate the potential of graph representation learning for the interpretation of biological complexity and pave the way for studying core gene properties and future drug development.

Characterizing the targets of transcription regulators by aggregating ChIP-seq and perturbation expression data sets

Mapping the gene targets of chromatin-associated transcription regulators (TRs) is a major goal of genomics research. ChIP-seq of TRs and experiments that perturb a TR and measure the differential abundance of gene transcripts are a primary means by which direct relationships are tested on a genomic scale. It has been reported that there is a poor overlap in the evidence across gene regulation strategies, emphasizing the need for integrating results from multiple experiments. Although research consortia interested in gene regulation have produced a valuable trove of high-quality data, there is an even greater volume of TR-specific data throughout the literature. In this study, we show a workflow for the identification, uniform processing, and aggregation of ChIP-seq and TR perturbation experiments for the ultimate purpose of ranking human and mouse TR-target interactions. Focusing on an initial set of eight regulators (ASCL1, HES1, MECP2, MEF2C, NEUROD1, PAX6, RUNX1, and TCF4), we identified 497 experiments suitable for analysis. We used this corpus to examine data concordance, to identify systematic patterns of the two data types, and to identify putative orthologous interactions between human and mouse. We build upon commonly used strategies to forward a procedure for aggregating and combining these two genomic methodologies, assessing these rankings against independent literature-curated evidence. Beyond a framework extensible to other TRs, our work also provides empirically ranked TR-target listings, as well as transparent experiment-level gene summaries for community use.

Overnight Corticosterone and Gene Expression in Mouse Hippocampus: Time Course during Resting Period

The aim of the experiment was to test the effect of an elevated level of glucocorticoids on the mouse hippocampal transcriptome after 12 h of treatment with corticosterone that was administered during an active phase of the circadian cycle. Additionally, we also tested the circadian changes in gene expression and the decay time of transcriptomic response to corticosterone. Gene expression was analyzed using microarrays. Obtained results show that transcriptomic responses to glucocorticoids are heterogeneous in terms of the decay time with some genes displaying persistent changes in expression during 9 h of rest. We have also found a considerable overlap between genes regulated by corticosterone and genes implicated previously in stress response. The examples of such genes are Acer2, Agt, Apod, Aqp4, Etnppl, Fabp7, Fam107a, Fjx1, Fmo2, Galnt15, Gjc2, Heph, Hes5, Htra1, Jdp2, Kif5a, Lfng, Lrg1, Mgp, Mt1, Pglyrp1, Pla2g3, Plin4, Pllp, Ptgds, Ptn, Slc2a1, Slco1c1, Sult1a1, Thbd and Txnip. This indicates that the applied model is a useful tool for the investigation of mechanisms underlying the stress response.

The Neuroscience Multi-Omic Archive: a BRAIN Initiative resource for single-cell transcriptomic and epigenomic data from the mammalian brain

Scalable technologies to sequence the transcriptomes and epigenomes of single cells are transforming our understanding of cell types and cell states. The Brain Research through Advancing Innovative Neurotechnologies (BRAIN) Initiative Cell Census Network (BICCN) is applying these technologies at unprecedented scale to map the cell types in the mammalian brain. In an effort to increase data FAIRness (Findable, Accessible, Interoperable, Reusable), the NIH has established repositories to make data generated by the BICCN and related BRAIN Initiative projects accessible to the broader research community. Here, we describe the Neuroscience Multi-Omic Archive (NeMO Archive; nemoarchive.org), which serves as the primary repository for genomics data from the BRAIN Initiative. Working closely with other BRAIN Initiative researchers, we have organized these data into a continually expanding, curated repository, which contains transcriptomic and epigenomic data from over 50 million brain cells, including single-cell genomic data from all of the major regions of the adult and prenatal human and mouse brains, as well as substantial single-cell genomic data from non-human primates. We make available several tools for accessing these data, including a searchable web portal, a cloud-computing interface for large-scale data processing (implemented on Terra, terra.bio), and a visualization and analysis platform, NeMO Analytics (nemoanalytics.org).

Stress and the brain transcriptome: Identifying commonalities and clusters in standardized data from published experiments

Interpretation of transcriptomic experiments is hindered by many problems including false positives/negatives inherent to big-data methods and changes in gene nomenclature. To find the most consistent effect of stress on brain transcriptome, we retrieved data from 79 studies applying animal models and 3 human studies investigating post-traumatic stress disorder (PTSD). The analyzed data were obtained either with microarrays or RNA sequencing applied to samples collected from more than 1887 laboratory animals and from 121 human subjects. Based on the initial database containing a quarter million differential expression effect sizes representing transcripts in three species, we identified the most frequently reported genes in 223 stress-control comparisons. Additionally, the analysis considers sex, individual vulnerability and contribution of glucocorticoids. We also found an overlap between gene expression in PTSD patients and animals which indicates relevance of laboratory models for human stress response. Our analysis points to genes that, as far as we know, were not specifically tested for their role in stress response (Pllp, Arrdc2, Midn, Mfsd2a, Ccn1, Htra1, Csrnp1, Tenm4, Tnfrsf25, Sema3b, Fmo2, Adamts4, Gjb1, Errfi1, Fgf18, Galnt6, Slc25a42, Ifi30, Slc4a1, Cemip, Klf10, Tom1, Dcdc2c, Fancd2, Luzp2, Trpm1, Abcc12, Osbpl1a, Ptp4a2). Provided transcriptomic resource will be useful for guiding the new research.

Elucidating gene expression patterns across multiple biological contexts through a large-scale investigation of transcriptomic datasets

Distinct gene expression patterns within cells are foundational for the diversity of functions and unique characteristics observed in specific contexts, such as human tissues and cell types. Though some biological processes commonly occur across contexts, by harnessing the vast amounts of available gene expression data, we can decipher the processes that are unique to a specific context. Therefore, with the goal of developing a portrait of context-specific patterns to better elucidate how they govern distinct biological processes, this work presents a large-scale exploration of transcriptomic signatures across three different contexts (i.e., tissues, cell types, and cell lines) by leveraging over 600 gene expression datasets categorized into 98 subcontexts. The strongest pairwise correlations between genes from these subcontexts are used for the construction of co-expression networks. Using a network-based approach, we then pinpoint patterns that are unique and common across these subcontexts. First, we focused on patterns at the level of individual nodes and evaluated their functional roles using a human protein-protein interactome as a referential network. Next, within each context, we systematically overlaid the co-expression networks to identify specific and shared correlations as well as relations already described in scientific literature. Additionally, in a pathway-level analysis, we overlaid node and edge sets from co-expression networks against pathway knowledge to identify biological processes that are related to specific subcontexts or groups of them. Finally, we have released our data and scripts at https://zenodo.org/record/5831786 and https://github.com/ContNeXt/ , respectively and developed ContNeXt ( https://contnext.scai.fraunhofer.de/ ), a web application to explore the networks generated in this work.

Atlas of clinically distinct cell states and ecosystems across human solid tumors

Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.

Construction of a Potentially Functional circRNA-miRNA-mRNA Network in Intervertebral Disc Degeneration by Bioinformatics Analysis

Background

The competing endogenous RNA- (ceRNA-) mediated regulatory mechanisms are known to play a pivotal role in intervertebral disc degeneration (IDD). Our research intended to establish a ceRNA regulatory network related to IDD through bioinformatics analyses.

Methods

The expression profiles of circRNA, miRNA, and mRNA were obtained from the public Gene Expression Omnibus (GEO) datasets. Then, we use sequence-based bioinformatics methods to select differentially expressed mRNAs (DEmRNAs), microRNAs (DEmiRNAs), or circRNAs (DEcircRNAs) related to IDD. We used ChEA3 to verify the targets of transcription factors (TFs). Then, we used DAVID to annotate the DEmRNAs. Finally, we constructed a potentially circRNA-miRNA-mRNA network related to IDD by predicting in the database (ENCORI, TargetScan, miRecords, miRmap, and circBank).

Results

We identified 31 common DEmRNAs by Venn analysis, of which MMP2 was regarded as the key hub genes. Simultaneously, miR-423-5p and miR-185-5p were predicted as the upstream molecules of MMP2. Furthermore, a total of six DEcircRNAs were predicted as the upstream circRNAs of miR-423-5p and miR-185-5p. Then, a potential circRNA-miRNA-mRNA network related to IDD was constructed by bioinformatics analysis.

Conclusion

A comprehensive ceRNA regulatory network was constructed, which was found to be significant in IDD progression.

A pH-eQTL Interaction at the RIT2-SYT4 Parkinson's Disease Risk Locus in the Substantia Nigra

Parkinson's disease causes severe motor and cognitive disabilities that result from the progressive loss of dopamine neurons in the substantia nigra. The rs12456492 variant in the RIT2 gene has been repeatedly associated with increased risk for Parkinson's disease. From a transcriptomic perspective, a meta-analysis found that RIT2 gene expression is correlated with pH in the human brain. To assess these pH associations in relation to Parkinson's disease risk, we examined the two datasets that assayed rs12456492, gene expression, and pH in the postmortem human brain. Using the BrainEAC dataset, we replicate the positive correlation between RIT2 gene expression and pH in the human brain (n = 100). Furthermore, we found that the relationship between expression and pH is influenced by rs12456492. When tested across ten brain regions, this interaction is specifically found in the substantia nigra. A similar association was found for the co-localized SYT4 gene. In addition, SYT4 associations are stronger in a combined model with both genes, and the SYT4 interaction appears to be specific to males. In the Genotype-Tissue Expression (GTEx) dataset, the pH associations involving rs12456492 and expression of either SYT4 and RIT2 were not seen. This null finding may be due to the short postmortem intervals of the GTEx tissue samples. In the BrainEAC data, we tested the effect of postmortem interval and only observed the interactions in samples with the longer intervals. These previously unknown associations suggest novel roles for rs12456492, RIT2, and SYT4 in the regulation and response to pH in the substantia nigra.

Characterization of Human Genes Modulated by Porphyromonas gingivalis Highlights the Ribosome, Hypothalamus, and Cholinergic Neurons

Porphyromonas gingivalis, a bacterium associated with periodontal disease, is a suspected cause of Alzheimer's disease. This bacterium is reliant on gingipain proteases, which cleave host proteins after arginine and lysine residues. To characterize gingipain susceptibility, we performed enrichment analyses of arginine and lysine proportion proteome-wide. Genes differentially expressed in brain samples with detected P. gingivalis reads were also examined. Genes from these analyses were tested for functional enrichment and specific neuroanatomical expression patterns. Proteins in the SRP-dependent cotranslational protein targeting to membrane pathway were enriched for these residues and previously associated with periodontal and Alzheimer's disease. These ribosomal genes are up-regulated in prefrontal cortex samples with detected P. gingivalis sequences. Other differentially expressed genes have been previously associated with dementia (ITM2B, MAPT, ZNF267, and DHX37). For an anatomical perspective, we characterized the expression of the P. gingivalis associated genes in the mouse and human brain. This analysis highlighted the hypothalamus, cholinergic neurons, and the basal forebrain. Our results suggest markers of neural P. gingivalis infection and link the cholinergic and gingipain hypotheses of Alzheimer's disease.

Estimating and Correcting for Off-Target Cellular Contamination in Brain Cell Type Specific RNA-Seq Data

Transcriptionally profiling minor cellular populations remains an ongoing challenge in molecular genomics. Single-cell RNA sequencing has provided valuable insights into a number of hypotheses, but practical and analytical challenges have limited its widespread adoption. A similar approach, which we term single-cell type RNA sequencing (sctRNA-seq), involves the enrichment and sequencing of a pool of cells, yielding cell type-level resolution transcriptomes. While this approach offers benefits in terms of mRNA sampling from targeted cell types, it is potentially affected by off-target contamination from surrounding cell types. Here, we leveraged single-cell sequencing datasets to apply a computational approach for estimating and controlling the amount of off-target cell type contamination in sctRNA-seq datasets. In datasets obtained using a number of technologies for cell purification, we found that most sctRNA-seq datasets tended to show some amount of off-target mRNA contamination from surrounding cells. However, using covariates for cellular contamination in downstream differential expression analyses increased the quality of our models for differential expression analysis in case/control comparisons and typically resulted in the discovery of more differentially expressed genes. In general, our method provides a flexible approach for detecting and controlling off-target cell type contamination in sctRNA-seq datasets.