Ultrafast functional profiling of RNA-seq data for nonmodel organisms

MangroveDB: A Comprehensive Online Database for Mangroves Based on Multi-Omics Data

Mangroves are dominant flora of intertidal zones along tropical and subtropical coastline around the world that offer important ecological and economic value. Recently, the genomes of mangroves have been decoded, and massive omics data were generated and deposited in the public databases. Reanalysis of multi-omics data can provide new biological insights excluded in the original studies. However, the requirements for computational resource and lack of bioinformatics skill for experimental researchers limit the effective use of the original data. To fill this gap, we uniformly processed 942 transcriptome data, 386 whole-genome sequencing data, and provided 13 reference genomes and 40 reference transcriptomes for 53 mangroves. Finally, we built an interactive web-based database platform MangroveDB (https://github.com/Jasonxu0109/MangroveDB), which was designed to provide comprehensive gene expression datasets to facilitate their exploration and equipped with several online analysis tools, including principal components analysis, differential gene expression analysis, tissue-specific gene expression analysis, GO and KEGG enrichment analysis. MangroveDB not only provides query functions about genes annotation, but also supports some useful visualization functions for analysis results, such as volcano plot, heatmap, dotplot, PCA plot, bubble plot, population structure, and so on. In conclusion, MangroveDB is a valuable resource for the mangroves research community to efficiently use the massive public omics datasets.

DNA-protein quasi-mapping for rapid differential gene expression analysis in non-model organisms

BACKGROUND

Conventional differential gene expression analysis pipelines for non-model organisms require computationally expensive transcriptome assembly. We recently proposed an alternative strategy of directly aligning RNA-seq reads to a protein database, and demonstrated drastic improvements in speed, memory usage, and accuracy in identifying differentially expressed genes.

RESULT

Here we report a further speed-up by replacing DNA-protein alignment by quasi-mapping, making our pipeline > 1000× faster than assembly-based approach, and still more accurate. We also compare quasi-mapping to other mapping techniques, and show that it is faster but at the cost of sensitivity.

CONCLUSION

We provide a quick-and-dirty differential gene expression analysis pipeline for non-model organisms without a reference transcriptome, which directly quasi-maps RNA-seq reads to a reference protein database, avoiding computationally expensive transcriptome assembly.

Comparative Transcriptome Analysis of Bovine, Porcine, and Sheep Muscle Using Interpretable Machine Learning Models

The growth and development of muscle tissue play a pivotal role in the economic value and quality of meat in agricultural animals, garnering close attention from breeders and researchers. The quality and palatability of muscle tissue directly determine the market competitiveness of meat products and the satisfaction of consumers. Therefore, a profound understanding and management of muscle growth is essential for enhancing the overall economic efficiency and product quality of the meat industry. Despite this, systematic research on muscle development-related genes across different species still needs to be improved. This study addresses this gap through extensive cross-species muscle transcriptome analysis, combined with interpretable machine learning models. Utilizing a comprehensive dataset of 275 publicly available transcriptomes derived from porcine, bovine, and ovine muscle tissues, encompassing samples from ten distinct muscle types such as the semimembranosus and longissimus dorsi, this study analyzes 113 porcine (n = 113), 94 bovine (n = 94), and 68 ovine (n = 68) specimens. We employed nine machine learning models, such as Support Vector Classifier (SVC) and Support Vector Machine (SVM). Applying the SHapley Additive exPlanations (SHAP) method, we analyzed the muscle transcriptome data of cattle, pigs, and sheep. The optimal model, adaptive boosting (AdaBoost), identified key genes potentially influencing muscle growth and development across the three species, termed SHAP genes. Among these, 41 genes (including NANOG, ADAMTS8, LHX3, and TLR9) were consistently expressed in all three species, designated as homologous genes. Specific candidate genes for cattle included SLC47A1, IGSF1, IRF4, EIF3F, CGAS, ZSWIM9, RROB1, and ABHD18; for pigs, DRP2 and COL12A1; and for sheep, only COL10A1. Through the analysis of SHAP genes utilizing Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, relevant pathways such as ether lipid metabolism, cortisol synthesis and secretion, and calcium signaling pathways have been identified, revealing their pivotal roles in muscle growth and development.

Propofol Protects the Blood-Brain Barrier After Traumatic Brain Injury by Stabilizing the Extracellular Matrix via Prrx1: From Neuroglioma to Neurotrauma

The purpose of this study is to explore the shared molecular pathogenesis of traumatic brain injury (TBI) and high-grade glioma and investigate the mechanism of propofol (PF) as a potential protective agent. By analyzing the Chinese glioma genome atlas (CGGA) and The Cancer Genome Atlas (TCGA) databases, we compared the transcriptomic data of high-grade glioma and TBI patients to identify common pathological mechanisms. Through bioinformatics analysis, in vitro experiments and in vivo TBI model, we investigated the regulatory effect of PF on extracellular matrix (ECM)-related genes through Prrx1 under oxidative stress. The impact of PF on BBB integrity under oxidative stress was investigated using a dual-layer BBB model, and we explored the protective effect of PF on tight junction proteins and ECM-related genes in mice after TBI. The study found that high-grade glioma and TBI share ECM instability as an important molecular pathological mechanism. PF stabilizes the ECM and protects the BBB by directly binding to Prrx1 or indirectly regulating Prrx1 through miRNAs. In addition, PF reduces intracellular calcium ions and ROS levels under oxidative stress, thereby preserving BBB integrity. In a TBI mouse model, PF protected BBB integrity through up-regulated tight junction proteins and stabilized the expression of ECM-related genes. Our study reveals the shared molecular pathogenesis between TBI and glioblastoma and demonstrate the potential of PF as a protective agent of BBB. This provides new targets and approaches for the development of novel neurotrauma therapeutic drugs.

Simple protocol for combined extraction of exocrine secretions and RNA in small arthropods

The integration of data from multiple sources and analytical techniques to obtain novel insights and answer challenging questions is a hallmark of modern science. In arthropods, exocrine secretions may act as pheromones, defensive substances, antibiotics, as well as surface protectants, and as such they play a crucial role in ecology and evolution. Exocrine chemical compounds are frequently characterized by gas chromatography-mass spectrometry. Technological advances of recent years now allow us to routinely characterize the total gene complement transcribed in a particular biological tissue, often in the context of experimental treatment, via RNAseq. We here introduce a novel methodological approach to successfully characterize exocrine secretions and full transcriptomes of one and the same individual of oribatid mites. We found that chemical extraction prior to RNA extraction had only minor effects on the total RNA integrity. De novo transcriptomes obtained from such combined extractions were of comparable quality to those assembled for samples that were subject to RNA extraction only, indicating that combined chemical/RNA extraction is perfectly suitable for phylotranscriptomic studies. However, in-depth analysis of RNA expression analysis indicates that chemical extraction prior to RNAseq may affect transcript degradation rates, similar to the effects reported in previous studies comparing RNA extraction protocols. With this pilot study, we demonstrate that profiling chemical secretions and RNA expression levels from the same individual is methodologically feasible, paving the way for future research to understand the genes and pathways underlying the syntheses of biogenic chemical compounds. Our approach should be applicable broadly to most arachnids, insects, and other arthropods.

Wildlife ecological risk assessment in the 21st century: Promising technologies to assess toxicological effects

Despite advances in toxicity testing and the development of new approach methodologies (NAMs) for hazard assessment, the ecological risk assessment (ERA) framework for terrestrial wildlife (i.e., air-breathing amphibians, reptiles, birds, and mammals) has remained unchanged for decades. While survival, growth, and reproductive endpoints derived from whole-animal toxicity tests are central to hazard assessment, nonstandard measures of biological effects at multiple levels of biological organization (e.g., molecular, cellular, tissue, organ, organism, population, community, ecosystem) have the potential to enhance the relevance of prospective and retrospective wildlife ERAs. Other factors (e.g., indirect effects of contaminants on food supplies and infectious disease processes) are influenced by toxicants at individual, population, and community levels, and need to be factored into chemically based risk assessments to enhance the "eco" component of ERAs. Regulatory and logistical challenges often relegate such nonstandard endpoints and indirect effects to postregistration evaluations of pesticides and industrial chemicals and contaminated site evaluations. While NAMs are being developed, to date, their applications in ERAs focused on wildlife have been limited. No single magic tool or model will address all uncertainties in hazard assessment. Modernizing wildlife ERAs will likely entail combinations of laboratory- and field-derived data at multiple levels of biological organization, knowledge collection solutions (e.g., systematic review, adverse outcome pathway frameworks), and inferential methods that facilitate integrations and risk estimations focused on species, populations, interspecific extrapolations, and ecosystem services modeling, with less dependence on whole-animal data and simple hazard ratios. Integr Environ Assess Manag 2024;20:725-748. © 2023 His Majesty the King in Right of Canada and The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). Reproduced with the permission of the Minister of Environment and Climate Change Canada. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Web-based multi-omics integration using the Analyst software suite

The growing number of multi-omics studies demands clear conceptual workflows coupled with easy-to-use software tools to facilitate data analysis and interpretation. This protocol covers three key components involved in multi-omics analysis, including single-omics data analysis, knowledge-driven integration using biological networks and data-driven integration through joint dimensionality reduction. Using the dataset from a recent multi-omics study of human pancreatic islet tissue and plasma samples, the first section introduces how to perform transcriptomics/proteomics data analysis using ExpressAnalyst and lipidomics data analysis using MetaboAnalyst. On the basis of significant features detected in these workflows, the second section demonstrates how to perform knowledge-driven integration using OmicsNet. The last section illustrates how to perform data-driven integration from the normalized omics data and metadata using OmicsAnalyst. The complete protocol can be executed in ~2 h. Compared with other available options for multi-omics integration, the Analyst software suite described in this protocol enables researchers to perform a wide range of omics data analysis tasks via a user-friendly web interface.

Adverse Environmental Perturbations May Threaten Kelp Farming Sustainability by Exacerbating Enterobacterales Diseases

Globally kelp farming is gaining attention to mitigate land-use pressures and achieve carbon neutrality. However, the influence of environmental perturbations on kelp farming remains largely unknown. Recently, a severe disease outbreak caused extensive kelp mortality in Sanggou Bay, China, one of the world's largest high-density kelp farming areas. Here, through in situ investigations and simulation experiments, we find indications that an anomalously dramatic increase in elevated coastal seawater light penetration may have contributed to dysbiosis in the kelp Saccharina japonica's microbiome. This dysbiosis promoted the proliferation of opportunistic pathogenic Enterobacterales, mainly including the genera Colwellia and Pseudoalteromonas. Using transcriptomic analyses, we revealed that high-light conditions likely induced oxidative stress in kelp, potentially facilitating opportunistic bacterial Enterobacterales attack that activates a terrestrial plant-like pattern recognition receptor system in kelp. Furthermore, we uncover crucial genotypic determinants of Enterobacterales dominance and pathogenicity within kelp tissue, including pathogen-associated molecular patterns, potential membrane-damaging toxins, and alginate and mannitol lysis capability. Finally, through analysis of kelp-associated microbiome data sets under the influence of ocean warming and acidification, we conclude that such Enterobacterales favoring microbiome shifts are likely to become more prevalent in future environmental conditions. Our study highlights the need for understanding complex environmental influences on kelp health and associated microbiomes for the sustainable development of seaweed farming.

Hepatic Transcriptomic Responses to Ethinylestradiol in Embryonic Japanese Quail and Double-Crested Cormorant

Understanding species differences in sensitivity to toxicants is a critical issue in ecotoxicology. We recently established that double-crested cormorant (DCCO) embryos are more sensitive than Japanese quail (JQ) to the developmental effects of ethinylestradiol (EE2). We explored how this difference in sensitivity between species is reflected at a transcriptomic level. The EE2 was dissolved in dimethyl sulfoxide and injected into the air cell of eggs prior to incubation at nominal concentrations of 0, 3.33, and 33.3 µg/g egg weight. At midincubation (JQ 9 days; DCCO 16 days), livers were collected from five embryos/treatment group for RNA sequencing. Data were processed and analyzed using EcoOmicsAnalyst and ExpressAnalyst. The EE2 exposure dysregulated 238 and 1,987 genes in JQ and DCCO, respectively, with 78 genes in common between the two species. These included classic biomarkers of estrogen exposure such as vitellogenin and apovitellenin. We also report DCCO-specific dysregulation of Phase I/II enzyme-coding genes and species-specific transcriptional ontogeny of vitellogenin-2. Twelve Kyoto Encyclopedia of Genes and Genomes pathways and two EcoToxModules were dysregulated in common in both species including the peroxisome proliferator-activated receptor (PPAR) signaling pathway and fatty acid metabolism. Similar to previously reported differences at the organismal level, DCCO were more responsive to EE2 exposure than JQ at the gene expression level. Our description of differences in transcriptional responses to EE2 in early life stage birds may contribute to a better understanding of the molecular basis for species differences. Environ Toxicol Chem 2024;00:1-12. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Comparing Transcriptomic Responses to Chemicals Across Six Species Using the EcoToxChip RNASeq Database

The EcoToxChip project includes RNA-sequencing data from experiments involving model (Japanese quail, fathead minnow, African clawed frog) and ecological (double-crested cormorant, rainbow trout, northern leopard frog) species at multiple life stages (whole embryo and adult) exposed to eight chemicals of environmental concern known to perturb a wide range of biological systems (ethinyl estradiol, hexabromocyclododecane, lead, selenomethionine, 17β trenbolone, chlorpyrifos, fluoxetine, and benzo[a]pyrene). The objectives of this short communication were to (1) present and make available this RNA-sequencing database (i.e., 724 samples from 49 experiments) under the FAIR principles (FAIR data are data which meet principles of findability, accessibility, interoperability, and reusability), while also summarizing key meta-data attributes and (2) use ExpressAnalyst (including the Seq2Fun algorithm and EcoOmicsDB) to perform a comparative transcriptomics analysis of this database focusing on baseline and differential transcriptomic changes across species-life stage-chemical combinations. The database is available in NCBI GEO under accession number GSE239776. Across all species, the number of raw reads per sample ranged between 13 and 58 million, with 30% to 79% of clean reads mapped to the "vertebrate" subgroup database in EcoOmicsDB. Principal component analyses of the reads illustrated separation across the three taxonomic groups as well as some between tissue types. The most common differentially expressed gene was CYP1A1 followed by CTSE, FAM20CL, MYC, ST1S3, RIPK4, VTG1, and VIT2. The most common enriched pathways were metabolic pathways, biosynthesis of cofactors and biosynthesis of secondary metabolites, and chemical carcinogenesis, drug metabolism, and metabolism of xenobiotics by cytochrome P450. The RNA-sequencing database in the present study may be used by the research community for multiple purposes, including, for example, cross-species investigations, in-depth analyses of a particular test compound, and transcriptomic meta-analyses. Environ Toxicol Chem 2024;00:1-6. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Using ExpressAnalyst for Comprehensive Gene Expression Analysis in Model and Non-Model Organisms

ExpressAnalyst is a web-based platform that enables intuitive, end-to-end transcriptomics and proteomics data analysis. Users can start from FASTQ files, gene/protein abundance tables, or gene/protein lists. ExpressAnalyst will perform read quantification, gene expression table processing and normalization, differential expression analysis, or meta-analysis with complex study designs. The results are presented via various interactive visualizations such as volcano plots, heatmaps, networks, and ridgeline charts, with built-in functional enrichment analysis to allow flexible data exploration and understanding. ExpressAnalyst currently contains built-in support for 29 common organisms. For non-model organisms without good reference genomes, it can perform comprehensive transcriptome profiling directly from RNA-seq reads. These common tasks are covered in 11 Basic Protocols. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: RNA-seq count table uploading, processing, and normalization Basic Protocol 2: Differential expression analysis with linear models Basic Protocol 3: Functional analysis with volcano plot, enrichment network, and ridgeline visualization Basic Protocol 4: Hierarchical clustering analysis of transcriptomics data using interactive heatmaps Basic Protocol 5: Cross-species gene expression analysis based on ortholog mapping results Basic Protocol 6: Proteomics and microarray data processing and normalization Basic Protocol 7: Preparing multiple gene expression tables for meta-analysis Basic Protocol 8: Statistical and functional meta-analysis of gene expression data Basic Protocol 9: Functional analysis of transcriptomics signatures Basic Protocol 10: Dose-response and time-series data analysis Basic Protocol 11: RNA-seq reads processing and quantification with and without reference transcriptomes.

ExpressAnalyst: A unified platform for RNA-sequencing analysis in non-model species

The increasing application of RNA sequencing to study non-model species demands easy-to-use and efficient bioinformatics tools to help researchers quickly uncover biological and functional insights. We developed ExpressAnalyst ( www.expressanalyst.ca ), a web-based platform for processing, analyzing, and interpreting RNA-sequencing data from any eukaryotic species. ExpressAnalyst contains a series of modules that cover from processing and annotation of FASTQ files to statistical and functional analysis of count tables or gene lists. All modules are integrated with EcoOmicsDB, an ortholog database that enables comprehensive analysis for species without a reference transcriptome. By coupling ultra-fast read mapping algorithms with high-resolution ortholog databases through a user-friendly web interface, ExpressAnalyst allows researchers to obtain global expression profiles and gene-level insights from raw RNA-sequencing reads within 24 h. Here, we present ExpressAnalyst and demonstrate its utility with a case study of RNA-sequencing data from multiple non-model salamander species, including two that do not have a reference transcriptome.

Transcriptomic Points of Departure Calculated from Rainbow Trout Gill, Liver, and Gut Cell Lines Exposed to Methylmercury and Fluoxetine

Ethical and resource limitation concerns are pushing chemicals management to develop alternatives to animal testing strategies. The objective of our study was to determine whether transcriptomic point of departure (tPOD) values could be derived from studies that followed Organisation for Economic Co-operation and Development (OECD) Test No. 249 (rainbow trout gill cell line), as well as from studies on trout liver and gut cells. Gill, liver, and gut cell lines were exposed to methylmercury and fluoxetine. Concentrations causing 50% cytotoxicity (LC50) were derived, the whole transcriptome was sequenced, and gene tPOD and pathway benchmark dose (BMD) values were derived from transcriptomic dose-response analysis. Differences in LC50 and transcriptomic responses across the cell lines were noted. For methylmercury, the tPOD_mode values were 14.5, 20.5, and 17.8 ppb for the gill, liver, and gut cells, respectively. The most sensitive pathway (pathway BMDs in parentheses) was ferroptosis in the gill (3.1 ppb) and liver (3.5 ppb), and glutathione metabolism in the gut (6.6 ppb). For fluoxetine, the tPOD_mode values were 109.4, 108.4, and 97.4 ppb for the gill, liver, and gut cells, respectively. The most sensitive pathway was neurotrophin signaling in the gill (147 ppb) and dopaminergic signaling in the gut (86.3 ppb). For both chemicals, the gene tPOD and pathway BMD values were lower than cytotoxic concentrations in vitro, and within 10-fold below the in vivo LC50s. By bringing together transcriptomics and dose-response analysis with an OECD test method in three cell lines, the results help to establish an in vitro method yielding tPOD values that are hypothesized to be protective of in vivo concentrations associated with adverse outcomes, and also give insights into mechanisms of action. Environ Toxicol Chem 2022;41:1982-1992. © 2022 SETAC.

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.

De novo transcriptome assembly and annotation of parthenogenetic lizard Darevskia unisexualis and its parental ancestors Darevskia valentini and Darevskia raddei nairensis

Darevskia rock lizards include 29 sexual and seven parthenogenetic species of hybrid origin distributed in the Caucasus. All seven parthenogenetic species of the genus Darevskia were formed as a result of interspecific hybridization of only four sexual species. It remains unknown what are the main advantages of interspecific hybridization along with switching on parthenogenetic reproduction in evolution of reptiles. Data on whole transcriptome sequencing of parthenogens and their parental ancestors can provide value impact in solving this problem. Here we have sequenced ovary tissue transcriptomes from unisexual parthenogenetic lizard D. unisexualis and its parental bisexual ancestors to facilitate the subsequent annotation and to obtain the collinear characteristics for comparison with other lizard species. Here we report generated RNAseq data from total mRNA of ovary tissues of D. unisexualis, D. valentini and D. raddei with 58932755, 51634041 and 62788216 reads. Obtained RNA reads were assembled by Trinity assembler and 95141, 62123, 61836 contigs were identified with N50 values of 2409, 2801 and 2827 respectively. For further analysis top Gene Ontology terms were annotated for all species and transcript number was calculated. The raw data were deposited in the NCBI SRA database (BioProject PRJNA773939). The assemblies are available in Mendeley Data and can be accessed via doi:10.17632/rtd8cx7zc3.1.

Using Transcriptomics and Metabolomics to Understand Species Differences in Sensitivity to Chlorpyrifos in Japanese Quail and Double-Crested Cormorant Embryos

Modern 21st-century toxicity testing makes use of omics technologies to address critical questions in toxicology and chemical management. Of interest are questions relating to chemical mechanisms of toxicity, differences in species sensitivity, and translation of molecular effects to observable apical endpoints. Our study addressed these questions by comparing apical outcomes and multiple omics responses in early-life stage exposure studies with Japanese quail (Coturnix japonica) and double-crested cormorant (Phalacrocorax auritus), representing a model and ecological species, respectively. Specifically, we investigated the dose-dependent response of apical outcomes as well as transcriptomics and metabolomics in the liver of each species exposed to chlorpyrifos, a widely used organophosphate pesticide. Our results revealed a clear pattern of dose-dependent disruption of gene expression and metabolic profiles in Japanese quail but not double-crested cormorant at similar chlorpyrifos exposure concentrations. The difference in sensitivity between species was likely due to higher metabolic transformation of chlorpyrifos in Japanese quail compared to double-crested cormorant. The most impacted biological pathways after chlorpyrifos exposure in Japanese quail included hepatic metabolism, oxidative stress, endocrine disruption (steroid and nonsteroid hormones), and metabolic disease (lipid and fatty acid metabolism). Importantly, we show consistent responses across biological scales, suggesting that significant disruption at the level of gene expression and metabolite profiles leads to observable apical responses at the organism level. Our study demonstrates the utility of evaluating effects at multiple biological levels of organization to understand how modern toxicity testing relates to outcomes of regulatory relevance, while also highlighting important, yet poorly understood, species differences in sensitivity to chemical exposure. Environ Toxicol Chem 2021;40:3019-3033. © 2021 SETAC.