Meta-analysis of public RNA sequencing data of abscisic acid-related abiotic stresses in Arabidopsis thaliana.

Abiotic stresses such as drought, salinity, and cold negatively affect plant growth and crop productivity. Understanding the molecular mechanisms underlying plant responses to these stressors is essential for stress tolerance in crops. The plant hormone abscisic acid (ABA) is significantly increased upon abiotic stressors, inducing physiological responses to adapt to stress and regulate gene expression. Although many studies have examined the components of established stress signaling pathways, few have explored other unknown elements. This study aimed to identify novel stress-responsive genes in plants by performing a meta-analysis of public RNA sequencing (RNA-Seq) data in Arabidopsis thaliana, focusing on five ABA-related stress conditions (ABA, Salt, Dehydration, Osmotic, and Cold). The meta-analysis of 216 paired datasets from five stress conditions was conducted, and differentially expressed genes were identified by introducing a new metric, called TN [stress-treated (T) and non-treated (N)] score. We revealed that 14 genes were commonly upregulated and 8 genes were commonly downregulated across all five treatments, including some that were not previously associated with these stress responses. On the other hand, some genes regulated by salt, dehydration, and osmotic treatments were not regulated by exogenous ABA or cold stress, suggesting that they may be involved in the plant response to dehydration independent of ABA. Our meta-analysis revealed a list of candidate genes with unknown molecular mechanisms in ABA-dependent and ABA-independent stress responses. These genes could be valuable resources for selecting genome editing targets and potentially contribute to the discovery of novel stress tolerance mechanisms and pathways in plants.

Transcriptome Analysis Reveals Enhancement of Cardiogenesis-Related Signaling Pathways by S-nitroso-N-pivaloyl-D-penicillamine (SNPiP): Implications for Improved Diastolic Function and Cardiac Performance.

We previously reported a novel compound called S-nitroso-N-pivaloyl-D-penicillamine (SNPiP), which was screened from a group of nitric oxide (NO) donor compounds with a basic chemical structure of S-nitroso-N-acetylpenicillamine (SNAP), to activate the non-neuronal acetylcholine (NNA) system. SNPiP-treated mice exhibited improved cardiac output and enhanced diastolic function, without an increase in heart rate. The NNA-activating effects included increased resilience to ischemia, modulation of energy metabolism preference, and activation of angiogenesis. Here, we performed transcriptome analysis of SNPiP-treated mice ventricles to elucidate how SNPiP exerts beneficial effects on cardiac function. A time-course study (24 and 48 h after SNPiP administration) revealed that SNPiP initially induced Wnt and cGMP-protein kinase G (PKG) signaling pathways, along with upregulation of genes involved in cardiac muscle tissue development and oxytocin signaling pathway. We also observed enrichment of glycolysis-related genes in response to SNPiP treatment, resulting in a metabolic shift from oxidative phosphorylation to glycolysis, which was suggested by reduced cardiac glucose contents while maintaining ATP levels. Additionally, SNPiP significantly upregulated atrial natriuretic peptide (ANP) and sarcolipin (SLN), which play crucial roles in calcium handling and cardiac performance. These findings suggest that SNPiP may have therapeutic potential based on the pleiotropic mechanisms elucidated in this study.

An extract from the frass of swallowtail butterfly (Papilio machaon) larvae inhibits HCT116 colon cancer cell proliferation but not other cancer cell types.

BACKGROUND: The frass of several herbivorous insect species has been utilised as natural medicines in Asia; however, the metabolite makeup and pharmaceutical activities of insect frass have yet to be investigated. Oligophagous Papilionidae insects utilise specific kinds of plants, and it has been suggested that the biochemicals from the plants may be metabolised by cytochrome P450 (CYP) in Papilionidae insects. In this study, we extracted the components of the frass of Papilio machaon larvae reared on Angelica keiskei, Oenanthe javanica or Foeniculum vulgare and examined the biological activity of each component. Then, we explored the expression of CYP genes in the midgut of P. machaon larvae and predicted the characteristics of their metabolic system. RESULTS: The components that were extracted using hexane, chloroform or methanol were biochemically different between larval frass and the host plants on which the larvae had fed. Furthermore, a fraction obtained from the chloroform extract from frass of A. keiskei-fed larvae specifically inhibited the cell proliferation of the human colon cancer cell line HCT116, whereas fractions obtained from the chloroform extracts of O. javanica- or F. vulgare-fed larval frass did not affect HCT116 cell viability. The metabolites from the chloroform extract from frass of A. keiskei-fed larvae prevented cell proliferation and induced apoptosis in HCT116 cells. Next, we explored the metabolic enzyme candidates in A. keiskei-fed larvae by RNA-seq analysis. We found that the A. keiskei-fed larval midgut might have different characteristics from the O. javanica- or F. vulgare-fed larval metabolic systems, and we found that the CYP6B2 transcript was highly expressed in the A. keiskei-fed larval midgut. CONCLUSIONS: These findings indicate that P. machaon metabolites might be useful as pharmaceutical agents against human colon cancer subtypes. Importantly, our findings show that it might be possible to use insect metabolic enzymes for the chemical structural conversion of plant-derived compounds with complex structures.

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 Heat-Stressed Transcriptomes Using the Public Gene Expression Database from Human and Mouse Samples.

Climate change has significantly increased the frequency of our exposure to heat, adversely affecting human health and industries. Heat stress is an environmental stress defined as the exposure of organisms and cells to abnormally high temperatures. To comprehensively explain the mechanisms underlying an organism's response to heat stress, it is essential to investigate and analyze genes that have been under-represented or less well-known in previous studies. In this study, we analyzed heat stress-responsive genes using a meta-analysis of numerous gene expression datasets from the public database. We obtained 322 human and 242 mouse pairs as the heat exposure and control data. The meta-analysis of these data identified 76 upregulated and 37 downregulated genes common to both humans and mice. We performed enrichment, protein-protein interaction network, and transcription factor target gene analyses for these genes. Furthermore, we conducted an integrated analysis of these genes using publicly available chromatin immunoprecipitation sequencing (ChIP-seq) data for HSF1, HSF2, and PPARGC1A (PGC-1α) as well as gene2pubmed data from the existing literature. The results identified previously overlooked genes, such as ABHD3, ZFAND2A, and USPL1, as commonly upregulated genes. Further functional analysis of these genes can contribute to coping with climate change and potentially lead to technological advancements.

Cellular senescence induction leads to progressive cell death via the INK4a-RB pathway in naked mole-rats.

Naked mole-rats (NMRs) have exceptional longevity and are resistant to age-related physiological decline and diseases. Given the role of cellular senescence in aging, we postulated that NMRs possess unidentified species-specific mechanisms to prevent senescent cell accumulation. Here, we show that upon induction of cellular senescence, NMR fibroblasts underwent delayed and progressive cell death that required activation of the INK4a-retinoblastoma protein (RB) pathway (termed "INK4a-RB cell death"), a phenomenon not observed in mouse fibroblasts. Naked mole-rat fibroblasts uniquely accumulated serotonin and were inherently vulnerable to hydrogen peroxide (H2 O2 ). After activation of the INK4a-RB pathway, NMR fibroblasts increased monoamine oxidase levels, leading to serotonin oxidization and H2 O2 production, which resulted in increased intracellular oxidative damage and cell death activation. In the NMR lung, induction of cellular senescence caused delayed, progressive cell death mediated by monoamine oxidase activation, thereby preventing senescent cell accumulation, consistent with in vitro results. The present findings indicate that INK4a-RB cell death likely functions as a natural senolytic mechanism in NMRs, providing an evolutionary rationale for senescent cell removal as a strategy to resist aging.

Meta-Analysis of Public RNA Sequencing Data Revealed Potential Key Genes Associated with Reproductive Division of Labor in Social Hymenoptera and Termites.

Eusociality in insects has evolved independently many times. One of the most notable characteristics of eusociality is the reproductive division of labor. In social insects, the reproductive division of labor is accomplished by queens and workers. Transcriptome analyses of queens and workers have been conducted for various eusocial species. However, the genes that regulate the reproductive division of labor across all or multiple eusocial species have not yet been fully elucidated. Therefore, we conducted a meta-analysis using publicly available RNA-sequencing data from four major groups of social insects. In this meta-analysis, we collected 258 pairs (queen vs. worker) of RNA-sequencing data from 34 eusocial species. The meta-analysis identified a total of 20 genes that were differentially expressed in queens or workers. Out of these, 12 genes have not previously been reported to be involved in the reproductive division of labor. Functional annotation of these 20 genes in other organisms revealed that they could be regulators of behaviors and physiological states related to the reproductive division of labor. These 20 genes, revealed using massive datasets of numerous eusocial insects, may be key regulators of the reproductive division of labor.

Exploratory meta-analysis of hypoxic transcriptomes using a precise transcript reference sequence set.

Gene expression studies are intrinsically biased, with many studies influenced by concomitant information such as gene-disease associations. This limitation can be overcome using a data-driven analysis approach without relying on ancillary information. The FANTOM CAGE-Associated Transcriptome project provides a comprehensive meta-assembly of the human transcriptome using coding and noncoding genes. Hypoxia strongly influences gene expression; in addition, noncoding RNA (ncRNA) metabolism is down-regulated in response to hypoxic stimuli. We evaluated the differential response of various transcripts to hypoxia by determining their hypoxia responsiveness scores. Enrichment analysis revealed that several genes associated with ncRNA metabolism, particularly those involved in ribosomal RNA processing, were down-regulated in response to hypoxia. Previously published information from the FANTOM CAGE-Associated Transcriptome project was suitable for meta-analysis of the transcriptome sequencing data from both coding and ncRNAs and to evaluate the hypoxia responsiveness of target transcripts and relationship between sense-antisense transcripts from the same locus. Our results may facilitate functional annotation of various transcripts including ncRNAs, allowing for both sense and antisense and coding and noncoding evaluations.

A highly contiguous genome assembly of red perilla (Perilla frutescens) domesticated in Japan.

Perilla frutescens (Lamiaceae) is an important herbal plant with hundreds of bioactive chemicals, among which perillaldehyde and rosmarinic acid are the two major bioactive compounds in the plant. The leaves of red perilla are used as traditional Kampo medicine or food ingredients. However, the medicinal and nutritional uses of this plant could be improved by enhancing the production of valuable metabolites through the manipulation of key enzymes or regulatory genes using genome editing technology. Here, we generated a high-quality genome assembly of red perilla domesticated in Japan. A near-complete chromosome-level assembly of P. frutescens was generated contigs with N50 of 41.5 Mb from PacBio HiFi reads. 99.2% of the assembly was anchored into 20 pseudochromosomes, among which seven pseudochromosomes consisted of one contig, while the rest consisted of less than six contigs. Gene annotation and prediction of the sequences successfully predicted 86,258 gene models, including 76,825 protein-coding genes. Further analysis showed that potential targets of genome editing for the engineering of anthocyanin pathways in P. frutescens are located on the late-stage pathways. Overall, our genome assembly could serve as a valuable reference for selecting target genes for genome editing of P. frutescens.

Reference Genome Sequences of the Oriental Armyworm, Mythimna separata (Lepidoptera: Noctuidae).

Lepidopteran insects are an important group of animals, including those used as biochemical and physiological model species in the insect and silk industries as well as others that are major agricultural pests. Therefore, the genome sequences of several lepidopteran insects have been reported. The oriental armyworm, Mythimna separata, is an agricultural pest commonly used to study insect immune reactions and interactions with parasitoid wasps as hosts. To improve our understanding of these research topics, reference genome sequences were constructed in the present study. Using long-read and short-read sequence data, de novo assembly and polishing were performed and haplotigs were purged. Subsequently, gene predictions and functional annotations were performed. To search for orthologs of the Toll and Immune Deficiency (IMD) pathways and for C-type lectins, annotation data analysis, BLASTp, and Hummer scans were performed. The M. separata genome is 682 Mbp; its contig N50 was 2.7 Mbp, with 21,970 genes and 24,452 coding sites predicted. All orthologs of the core components of the Toll and IMD pathways and 105 C-type lectins were identified. These results suggest that the genome data were of sufficient quality for use as reference genome data and could contribute to promoting M. separata and lepidopteran research at the molecular and genome levels.

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.

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.

Revealing Landscapes of Transposable Elements in Apis Species by Meta-Analysis.

Transposable elements (TEs) are grouped into several families with diverse sequences. Owing to their diversity, studies involving the detection, classification, and annotation of TEs are difficult tasks. Moreover, simple comparisons of TEs among different species with different methods can lead to misinterpretations. The genome data of several honey bee (Apis) species are available in public databases. Therefore, we conducted a meta-analysis of TEs, using 11 sets of genome data for Apis species, in order to establish data of "landscape of TEs". Consensus TE sequences were constructed and their distributions in the Apis genomes were determined. Our results showed that TEs belonged to four to seven TE families among 13 and 15 families of TEs detected in classes I and II respectively mainly consisted of Apis TEs and that more DNA/TcMar-Mariner consensus sequences and copies were present in all Apis genomes tested. In addition, more consensus sequences and copy numbers of DNA/TcMar-Mariner were detected in Apis mellifera than in other Apis species. These results suggest that TcMar-Mariner might exert A. mellifera-specific effects on the host A. mellifera species. In conclusion, our unified approach enabled comparison of Apis genome sequences to determine the TE landscape, which provide novel evolutionary insights into Apis species.

Prediction of sex-determination mechanisms in avian primordial germ cells using RNA-seq analysis.

In birds, sex is determined through cell-autonomous mechanisms and various factors, such as the dosage of DMRT1. While the sex-determination mechanism in gonads is well known, the mechanism in germ cells remains unclear. In this study, we explored the gene expression profiles of male and female primordial germ cells (PGCs) during embryogenesis in chickens to predict the mechanism underlying sex determination. Male and female PGCs were isolated from blood and gonads with a purity > 96% using flow cytometry and analyzed using RNA-seq. Prior to settlement in the gonads, female circulating PGCs (cPGCs) obtained from blood displayed sex-biased expression. Gonadal PGCs (gPGCs) also exhibited sex-biased expression, and the number of female-biased genes detected was higher than that of male-biased genes. The female-biased genes in gPGCs were enriched in some metabolic processes. To reveal the mechanisms underlying the transcriptional regulation of female-biased genes in gPGCs, we performed stimulation tests. Retinoic acid stimulation of cultured gPGCs derived from male embryos resulted in the upregulation of several female-biased genes. Overall, our results suggest that sex determination in avian PGCs involves aspects of both cell-autonomous and somatic-cell regulation. Moreover, it appears that sex determination occurs earlier in females than in males.

Genome Sequence of the Edible Green Alga Ulva prolifera, Originating from the Yoshinogawa River in Japan.

We report the genome sequence of Ulva prolifera, which originated from the Yoshinogawa River in Japan, using Oxford Nanopore Technologies MinION and Illumina sequencing reads. The genome assembly size is 103.8 Mbp, consisting of 142 scaffolds with an N50 value of 4.11 Mbp.

Meta-Analysis of RNA Sequencing Data of Arabidopsis and Rice under Hypoxia.

Hypoxia is an abiotic stress in plants. Flooding resulting from climate change is a major crop threat that increases the risk of hypoxic stress. The molecular mechanisms underlying hypoxia in plants were elucidated in recent years, but new genes related to this stress remain to be discovered. Thus, we aimed to perform a meta-analysis of the RNA sequencing (RNA-Seq) data of Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) under hypoxia. We collected 29 (Arabidopsis) and 26 (rice) pairs of RNA-Seq data involving hypoxic (including submergence) and normoxic (control) treatments and extracted the genes that were commonly upregulated or downregulated in the majority of the experiments. The meta-analysis revealed 40 and 19 commonly upregulated and downregulated genes, respectively, in the two species. Several WRKY transcription factors and cinnamate-4-hydroxylase were commonly upregulated, but their involvement in hypoxia remains unclear. Our meta-analysis identified candidate genes for novel molecular mechanisms in plants under hypoxia.

Systematic Functional Annotation Workflow for Insects.

Next-generation sequencing has revolutionized entomological study, rendering it possible to analyze the genomes and transcriptomes of non-model insects. However, use of this technology is often limited to obtaining the nucleotide sequences of target or related genes, with many of the acquired sequences remaining unused because other available sequences are not sufficiently annotated. To address this issue, we have developed a functional annotation workflow for transcriptome-sequenced insects to determine transcript descriptions, which represents a significant improvement over the previous method (functional annotation pipeline for insects). The developed workflow attempts to annotate not only the protein sequences obtained from transcriptome analysis but also the ncRNA sequences obtained simultaneously. In addition, the workflow integrates the expression-level information obtained from transcriptome sequencing for application as functional annotation information. Using the workflow, functional annotation was performed on the sequences obtained from transcriptome sequencing of the stick insect (Entoria okinawaensis) and silkworm (Bombyx mori), yielding richer functional annotation information than that obtained in our previous study. The improved workflow allows the more comprehensive exploitation of transcriptome data and is applicable to other insects because the workflow has been openly developed on GitHub.

Activation of transcription factor HIF inhibits IL-1ß-induced NO production in primary cultured rat hepatocytes.

Roxadustat and other hypoxia-inducible factor-prolyl hydroxylase inhibitors (HIF-PHIs) have recently been approved for the treatment of chronic renal anemia. In macrophages and monocytes, the activation of HIF-1 by pro-inflammatory cytokines induces iNOS expression and activity through the NF-κB pathway to produce nitric oxide (NO), which causes liver injury when excessively produced. Few studies have reported a relationship between HIF activity and iNOS induction in hepatocytes. We investigated the effect of drug- and hypoxia-induced HIF activations on NO production in primary cultured rat hepatocytes. Roxadustat treatment and hypoxic conditions activated HIF. Contrary to expectations, HIF-PHI treatment and hypoxia inhibited IL-1ß-induced NO production. RNA-Seq analysis of mRNA expression in rat hepatocytes showed that roxadustat treatment decreased the expression of genes related to inflammation, and genes in the NF-κB signaling pathway were induced by IL-1ß. Moreover, roxadustat suppressed IL-1ß-activated signaling pathways in an HIF-dependent manner. GalN/LPS-treated rats were used as in vivo models of hepatic injury, and roxadustat treatment showed a tendency to suppress the death of rats. Therefore, exogenous HIF-1 activation, including HIF-PHI and hypoxia exposures, suppressed IL-1ß-induced iNOS mRNA expression and subsequent NO production in hepatocytes, by suppressing the NF-κB signaling pathway. Roxadustat treatment suppresses the expression of pro-inflammatory genes by activating HIF, and thus may exhibit hepatoprotective effects.

Resistance to chemical carcinogenesis induction via a dampened inflammatory response in naked mole-rats.

Naked mole-rats (NMRs) have a very low spontaneous carcinogenesis rate, which has prompted studies on the responsible mechanisms to provide clues for human cancer prevention. However, it remains unknown whether and how NMR tissues respond to experimental carcinogenesis induction. Here, we show that NMRs exhibit extraordinary resistance against potent chemical carcinogenesis induction through a dampened inflammatory response. Although carcinogenic insults damaged skin cells of both NMRs and mice, NMR skin showed markedly lower immune cell infiltration. NMRs harbour loss-of-function mutations in RIPK3 and MLKL genes, which are essential for necroptosis, a type of necrotic cell death that activates strong inflammation. In mice, disruption of Ripk3 reduced immune cell infiltration and delayed carcinogenesis. Therefore, necroptosis deficiency may serve as a cancer resistance mechanism via attenuating the inflammatory response in NMRs. Our study sheds light on the importance of a dampened inflammatory response as a non-cell-autonomous cancer resistance mechanism in NMRs.

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

Analysis of RNA-sequencing (RNA-seq) data is an effective means to analyze the gene expression levels under specific conditions and discover new biological knowledge. More than 74,000 experimental series with RNA-seq have been stored in public databases as of 20 October 2021. Since this huge amount of expression data accumulated from past studies is a promising source of new biological insights, we focused on a meta-analysis of 1783 runs of RNA-seq data under the conditions of two types of stressors: oxidative stress (OS) and hypoxia. The collected RNA-seq data of OS were organized as the OS dataset to retrieve and analyze differentially expressed genes (DEGs). The OS-induced DEGs were compared with the hypoxia-induced DEGs retrieved from a previous study. The results from the meta-analysis of OS transcriptomes revealed two genes, CRIP1 and CRIP3, which were particularly downregulated, suggesting a relationship between OS and zinc homeostasis. The comparison between meta-analysis of OS and hypoxia showed that several genes were differentially expressed under both stress conditions, and it was inferred that the downregulation of cell cycle-related genes is a mutual biological process in both OS and hypoxia.