A versatile workflow to integrate RNA-seq genomic and transcriptomic data into mechanistic models of signaling pathways

Refined variant calling pipeline on RNA-seq data of breast cancer cell lines without matched-normal samples

OBJECTIVE

RNA-seq delivers valuable insights both to transcriptional patterns and mutational landscapes for transcribed genes. However, as tumour cell lines frequently lack their matched-normal counterpart, variant calling without the paired normal sample is still challenging. In order to exclude variants of common genetic variation without a matched-normal control, filtering strategies need to be developed to identify tumour relevant variants in cell lines.

RESULTS

Here, variants of 29 breast cancer cell lines were called on RNA-seq data via HaplotypeCaller. Low read depth sites, RNA-edit sites, and low complexity regions in coding regions were excluded. Common variants were filtered using 1000 genomes, gnomAD, and dbSNP data. Starting from hundred thousands of single nucleotide variants and small insertions and deletions, about thousand variants remained after filtering for each sample. Extracted variants were validated against the Catalogue of Somatic Mutations in Cancer (COSMIC) for 10 cell lines included in both data sets. Approximately half of the COSMIC variants were successfully called. Importantly, missing variants could mainly be attributed to sites with low read depth. Moreover, filtered variants also included all 10 cancer gene census COSMIC variants, a condensed hallmark variant set.

Polystyrene nanoplastics affect transcriptomic and epigenomic signatures of human fibroblasts and derived induced pluripotent stem cells: Implications for human health

Plastic pollution is increasing at an alarming rate yet the impact of this pollution on human health is poorly understood. Because human induced pluripotent stem cells (hiPSC) are frequently derived from dermal fibroblasts, these cells offer a powerful platform for the identification of molecular biomarkers of environmental pollution in human cells. Here, we describe a novel proof-of-concept for deriving hiPSC from human dermal fibroblasts deliberately exposed to polystyrene (PS) nanoplastic particles; unexposed hiPSC served as controls. In parallel, unexposed hiPSC were exposed to low and high concentrations of PS nanoparticles. Transcriptomic and epigenomic signatures of all fibroblasts and hiPSCs were defined using RNA-seq and whole genome methyl-seq, respectively. Both PS-treated fibroblasts and derived hiPSC showed alterations in expression of ESRRB and HNF1A genes and circuits involved in the pluripotency of stem cells, as well as in pathways involved in cancer, inflammatory disorders, gluconeogenesis, carbohydrate metabolism, innate immunity, and dopaminergic synapse. Similarly, the expression levels of identified key transcriptional and DNA methylation changes (DNMT3A, ESSRB, FAM133CP, HNF1A, SEPTIN7P8, and TTC34) were significantly affected in both PS-exposed fibroblasts and hiPSC. This study illustrates the power of human cellular models of environmental pollution to narrow down and prioritize the list of candidate molecular biomarkers of environmental pollution. This knowledge will facilitate the deciphering of the origins of environmental diseases.

A toolkit for enhanced reproducibility of RNASeq analysis for synthetic biologists

Sequencing technologies, in particular RNASeq, have become critical tools in the design, build, test and learn cycle of synthetic biology. They provide a better understanding of synthetic designs, and they help identify ways to improve and select designs. While these data are beneficial to design, their collection and analysis is a complex, multistep process that has implications on both discovery and reproducibility of experiments. Additionally, tool parameters, experimental metadata, normalization of data and standardization of file formats present challenges that are computationally intensive. This calls for high-throughput pipelines expressly designed to handle the combinatorial and longitudinal nature of synthetic biology. In this paper, we present a pipeline to maximize the analytical reproducibility of RNASeq for synthetic biologists. We also explore the impact of reproducibility on the validation of machine learning models. We present the design of a pipeline that combines traditional RNASeq data processing tools with structured metadata tracking to allow for the exploration of the combinatorial design in a high-throughput and reproducible manner. We then demonstrate utility via two different experiments: a control comparison experiment and a machine learning model experiment. The first experiment compares datasets collected from identical biological controls across multiple days for two different organisms. It shows that a reproducible experimental protocol for one organism does not guarantee reproducibility in another. The second experiment quantifies the differences in experimental runs from multiple perspectives. It shows that the lack of reproducibility from these different perspectives can place an upper bound on the validation of machine learning models trained on RNASeq data.

Graphical Abstract

SEAseq: a portable and cloud-based chromatin occupancy analysis suite

Background

Genome-wide protein-DNA binding is popularly assessed using specific antibody pulldown in Chromatin Immunoprecipitation Sequencing (ChIP-Seq) or Cleavage Under Targets and Release Using Nuclease (CUT&RUN) sequencing experiments. These technologies generate high-throughput sequencing data that necessitate the use of multiple sophisticated, computationally intensive genomic tools to make discoveries, but these genomic tools often have a high barrier to use because of computational resource constraints.

Results

We present a comprehensive, infrastructure-independent, computational pipeline called SEAseq, which leverages field-standard, open-source tools for processing and analyzing ChIP-Seq/CUT&RUN data. SEAseq performs extensive analyses from the raw output of the experiment, including alignment, peak calling, motif analysis, promoters and metagene coverage profiling, peak annotation distribution, clustered/stitched peaks (e.g. super-enhancer) identification, and multiple relevant quality assessment metrics, as well as automatic interfacing with data in GEO/SRA. SEAseq enables rapid and cost-effective resource for analysis of both new and publicly available datasets as demonstrated in our comparative case studies.

Conclusions

The easy-to-use and versatile design of SEAseq makes it a reliable and efficient resource for ensuring high quality analysis. Its cloud implementation enables a broad suite of analyses in environments with constrained computational resources. SEAseq is platform-independent and is aimed to be usable by everyone with or without programming skills. It is available on the cloud at https://platform.stjude.cloud/workflows/seaseq and can be locally installed from the repository at https://github.com/stjude/seaseq.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04588-z.

Human Embryos, Induced Pluripotent Stem Cells, and Organoids: Models to Assess the Effects of Environmental Plastic Pollution

For a long time, animal models were used to mimic human biology and diseases. However, animal models are not an ideal solution due to numerous interspecies differences between humans and animals. New technologies, such as human-induced pluripotent stem cells and three-dimensional (3D) cultures such as organoids, represent promising solutions for replacing, refining, and reducing animal models. The capacity of organoids to differentiate, self-organize, and form specific, complex, biologically suitable structures makes them excellent in vitro models of development and disease pathogenesis, as well as drug-screening platforms. Despite significant potential health advantages, further studies and considerable nuances are necessary before their clinical use. This article summarizes the definition of embryoids, gastruloids, and organoids and clarifies their appliance as models for early development, diseases, environmental pollution, drug screening, and bioinformatics.

Genome-scale mechanistic modeling of signaling pathways made easy: A bioconductor/cytoscape/web server framework for the analysis of omic data

Genome-scale mechanistic models of pathways are gaining importance for genomic data interpretation because they provide a natural link between genotype measurements (transcriptomics or genomics data) and the phenotype of the cell (its functional behavior). Moreover, mechanistic models can be used to predict the potential effect of interventions, including drug inhibitions. Here, we present the implementation of a mechanistic model of cell signaling for the interpretation of transcriptomic data as an R/Bioconductor package, a Cytoscape plugin and a web tool with enhanced functionality which includes building interpretable predictors, estimation of the effect of perturbations and assessment of the effect of mutations in complex scenarios.