Comparison of Oxidative and Hypoxic Stress Responsive Genes from Meta-Analysis of Public Transcriptomes

A systematic exploration of unexploited genes for oxidative stress in Parkinson's disease

Human disease-associated gene data are accessible through databases, including the Open Targets Platform, DisGeNET, miRTex, RNADisease, and PubChem. However, missing data entries in such databases are anticipated because of curational errors, biases, and text-mining failures. Additionally, the extensive research on human diseases has led to challenges in registering comprehensive data. The lack of essential data in databases hinders knowledge sharing and should be addressed. Therefore, we propose an analysis pipeline to explore missing entries of unexploited genes in the human disease-associated gene databases. Using this pipeline for genes in Parkinson's disease with oxidative stress revealed two unexploited genes: nuclear protein 1 (NUPR1) and ubiquitin-like with PHD and ring finger domains 2 (UHRF2). This methodology enhances the identification of underrepresented disease-associated genes, facilitating easier access to potential human disease-related functional genes. This study aims to identify unexploited genes for further research and does not include independent experimental validation.

Identification of drug responsible glycogene signature in liver carcinoma from meta-analysis using RNA-seq data

Glycans have attracted much attention in cancer therapeutic strategies, and cell surface proteins and lipids with glycans are known to be altered during the carcinogenic process. However, our understanding of how the glycogenes profile responds to drug stimulation remains incomplete. In this study, we search public databases for Sequence Read Archive data on drug-treated liver cancer cells, with the aim to comprehensively analyze the drug responses of glycogenes via bioinformatic meta-analysis. The study comprised 86 datasets, encompassing eight distinct liver cancer cell lines and 13 different drugs. Differentially expressed genes were quantified, and 399 glycogenes were identified. The glycogenes signature was then analyzed using bioinformatics methodologies. In the Protein-protein interaction network analysis, we identified drug-responsive glycogenes such as Beta-1,4-Galactosyltransferase 1, GDP-Mannose 4,6-Dehydratase, UDP-Glucose Ceramide Glucosyltransferase, and Solute Carrier Family 2 Member 4 as key glycan biomarkers. In the enrichment analysis using the pathway list of glycogenes, the results also demonstrated that drug stimulation resulted in alterations to glycopathway-related genes involved in several processes, namely O-Mannosylation, POMGNT2 Type, Capping, Heparan Sulfate Sulfation, and Glucuronidation pathways. These genes and pathways commonly exhibit variable expression across multiple liver cancer cells in response to the same drug, making them potential targets for new cancer therapies. In addition to their primary roles, drugs may also participate in the regulation of glycans. The insights from this study could pave the way for the development of liver cancer therapies that target the regulation of gene profiles involved in the biosynthesis of glycans.

DANGER analysis: risk-averse on/off-target assessment for CRISPR editing without a reference genome

Motivation

The CRISPR-Cas9 system has successfully achieved site-specific gene editing in organisms ranging from humans to bacteria. The technology efficiently generates mutants, allowing for phenotypic analysis of the on-target gene. However, some conventional studies did not investigate whether deleterious off-target effects partially affect the phenotype.

Results

Herein, we present a novel phenotypic assessment of CRISPR-mediated gene editing: Deleterious and ANticipatable Guides Evaluated by RNA-sequencing (DANGER) analysis. Using RNA-seq data, this bioinformatics pipeline can elucidate genomic on/off-target sites on mRNA-transcribed regions related to expression changes and then quantify phenotypic risk at the gene ontology term level. We demonstrated the risk-averse on/off-target assessment in RNA-seq data from gene-edited samples of human cells and zebrafish brains. Our DANGER analysis successfully detected off-target sites, and it quantitatively evaluated the potential contribution of deleterious off-targets to the transcriptome phenotypes of the edited mutants. Notably, DANGER analysis harnessed de novo transcriptome assembly to perform risk-averse on/off-target assessments without a reference genome. Thus, our resources would help assess genome editing in non-model organisms, individual human genomes, and atypical genomes from diseases and viruses. In conclusion, DANGER analysis facilitates the safer design of genome editing in all organisms with a transcriptome.

Availability and implementation

The Script for the DANGER analysis pipeline is available at https://github.com/KazukiNakamae/DANGER_analysis. In addition, the software provides a tutorial on reproducing the results presented in this article on the Readme page. The Docker image of DANGER_analysis is also available at https://hub.docker.com/repository/docker/kazukinakamae/dangeranalysis/general.

Meta-Analysis of the Public RNA-Seq Data of the Western Honeybee Apis mellifera to Construct Reference Transcriptome Data

The Western honeybee (Apis mellifera) is valuable in biological research and agriculture. Its genome sequence was published before those for other insect species. RNA-Seq data for A. mellifera have been applied in several recently published studies. Nevertheless, these data have not been prepared for use in subsequent meta-analyses. To promote A. mellifera transcriptome analysis, we constructed reference transcriptome data using the reference genome sequence and RNA-Seq data curated from about 1,000 runs of public databases. The new reference transcriptome data construct comprised 149,685 transcripts, and 194,174 protein sequences were predicted. Approximately 50-60% of the predicted protein sequences were functionally annotated using the protein sequence data for several model and insect species. Novel candidate immune-related transcripts were searched by meta-analysis using immune-response-related RNA-Seq and reference transcriptome data. Three to twenty candidate transcripts including autophagy-related protein 3 were upregulated or downregulated in response to both viral and bacterial infections. The constructed reference transcriptome data may facilitate future transcriptome analyses of A. mellifera.

Meta-Analysis of Transcriptomes in Insects Showing Density-Dependent Polyphenism

With increasing public data, a statistical analysis approach called meta-analysis, which combines transcriptome results obtained from multiple studies, has succeeded in providing novel insights into targeted biological processes. Locusts and aphids are representative of insect groups that exhibit density-dependent plasticity. Although the physiological mechanisms underlying density-dependent polyphenism have been identified in aphids and locusts, the underlying molecular mechanisms remain largely unknown. In this study, we performed a meta-analysis of public transcriptomes to gain additional insights into the molecular underpinning of density-dependent plasticity. We collected RNA sequencing data of aphids and locusts from public databases and detected differentially expressed genes (DEGs) between crowded and isolated conditions. Gene set enrichment analysis was performed to reveal the characteristics of the DEGs. DNA replication (GO:0006260), DNA metabolic processes (GO:0006259), and mitotic cell cycle (GO:0000278) were enriched in response to crowded conditions. To date, these processes have scarcely been the focus of research. The importance of the oxidative stress response and neurological system modifications under isolated conditions has been highlighted. These biological processes, clarified by meta-analysis, are thought to play key roles in the regulation of density-dependent plasticity.