Utility of Extrapolating Human S1500+ Genes to the Whole Transcriptome: Tunicamycin Case Study

Gene-environment interactions within a precision environmental health framework

Understanding the complex interplay of genetic and environmental factors in disease etiology and the role of gene-environment interactions (GEIs) across human development stages is important. We review the state of GEI research, including challenges in measuring environmental factors and advantages of GEI analysis in understanding disease mechanisms. We discuss the evolution of GEI studies from candidate gene-environment studies to genome-wide interaction studies (GWISs) and the role of multi-omics in mediating GEI effects. We review advancements in GEI analysis methods and the importance of large-scale datasets. We also address the translation of GEI findings into precision environmental health (PEH), showcasing real-world applications in healthcare and disease prevention. Additionally, we highlight societal considerations in GEI research, including environmental justice, the return of results to participants, and data privacy. Overall, we underscore the significance of GEI for disease prediction and prevention and advocate for integrating the exposome into PEH omics studies.

Substantiate a read-across hypothesis by using transcriptome data-A case study on volatile diketones

This case study explores the applicability of transcriptome data to characterize a common mechanism of action within groups of short-chain aliphatic α-, β-, and γ-diketones. Human reference in vivo data indicate that the α-diketone diacetyl induces bronchiolitis obliterans in workers involved in the preparation of microwave popcorn. The other three α-diketones induced inflammatory responses in preclinical in vivo animal studies, whereas beta and gamma diketones in addition caused neuronal effects. We investigated early transcriptional responses in primary human bronchiolar (PBEC) cell cultures after 24 h and 72 h of air-liquid exposure. Differentially expressed genes (DEGs) were assessed based on transcriptome data generated with the EUToxRisk gene panel of Temp-O-Seq^®. For each individual substance, genes were identified displaying a consistent differential expression across dose and exposure duration. The log fold change values of the DEG profiles indicate that α- and β-diketones are more active compared to γ-diketones. α-diketones in particular showed a highly concordant expression pattern, which may serve as a first indication of the shared mode of action. In order to gain a better mechanistic understanding, the resultant DEGs were submitted to a pathway analysis using ConsensusPathDB. The four α-diketones showed very similar results with regard to the number of activated and shared pathways. Overall, the number of signaling pathways decreased from α-to β-to γ-diketones. Additionally, we reconstructed networks of genes that interact with one another and are associated with different adverse outcomes such as fibrosis, inflammation or apoptosis using the TRANSPATH-database. Transcription factor enrichment and upstream analyses with the geneXplain platform revealed highly interacting gene products (called master regulators, MRs) per case study compound. The mapping of the resultant MRs on the reconstructed networks, visualized similar gene regulation with regard to fibrosis, inflammation and apoptosis. This analysis showed that transcriptome data can strengthen the similarity assessment of compounds, which is of particular importance, e.g., in read-across approaches. It is one important step towards grouping of compounds based on biological profiles.

The human hepatocyte TXG-MAPr: gene co-expression network modules to support mechanism-based risk assessment

Mechanism-based risk assessment is urged to advance and fully permeate into current safety assessment practices, possibly at early phases of drug safety testing. Toxicogenomics is a promising source of mechanisms-revealing data, but interpretative analysis tools specific for the testing systems (e.g. hepatocytes) are lacking. In this study, we present the TXG-MAPr webtool (available at https://txg-mapr.eu/WGCNA_PHH/TGGATEs_PHH/ ), an R-Shiny-based implementation of weighted gene co-expression network analysis (WGCNA) obtained from the Primary Human Hepatocytes (PHH) TG-GATEs dataset. The 398 gene co-expression networks (modules) were annotated with functional information (pathway enrichment, transcription factor) to reveal their mechanistic interpretation. Several well-known stress response pathways were captured in the modules, were perturbed by specific stressors and showed preservation in rat systems (rat primary hepatocytes and rat in vivo liver), with the exception of DNA damage and oxidative stress responses. A subset of 87 well-annotated and preserved modules was used to evaluate mechanisms of toxicity of endoplasmic reticulum (ER) stress and oxidative stress inducers, including cyclosporine A, tunicamycin and acetaminophen. In addition, module responses can be calculated from external datasets obtained with different hepatocyte cells and platforms, including targeted RNA-seq data, therefore, imputing biological responses from a limited gene set. As another application, donors' sensitivity towards tunicamycin was investigated with the TXG-MAPr, identifying higher basal level of intrinsic immune response in donors with pre-existing liver pathology. In conclusion, we demonstrated that gene co-expression analysis coupled to an interactive visualization environment, the TXG-MAPr, is a promising approach to achieve mechanistic relevant, cross-species and cross-platform evaluation of toxicogenomic data.

Reproducibility and robustness of high-throughput S1500+ transcriptomics on primary rat hepatocytes for chemical-induced hepatotoxicity assessment

•

TempO-Seq assays of rat hepatocytes in collagen sandwich are highly reproducible.

•

Gene expression analysis shows S1500+ is representative of the whole transcriptome.

•

Connectivity mapping shows consistency between TempO-Seq and Affymetrix data.

•

Gene set enrichment shows consistency between S1500+ and the whole transcriptome.

•

Gene set enrichment using hallmark gene sets informs hepatotoxicity.

Cell-based in vitro models coupled with high-throughput transcriptomics (HTTr) are increasingly utilized as alternative methods to animal-based toxicity testing. Here, using a panel of 14 chemicals with different risks of human drug-induced liver injury (DILI) and two dosing concentrations, we evaluated an HTTr platform comprised of collagen sandwich primary rat hepatocyte culture and the TempO-Seq surrogate S1500+ (ST) assay. First, the HTTr platform was found to exhibit high reproducibility between technical and biological replicates (r greater than 0.85). Connectivity mapping analysis further demonstrated a high level of inter-platform reproducibility between TempO-Seq data and Affymetrix GeneChip data from the Open TG-GATES project. Second, the TempO-Seq ST assay was shown to be a robust surrogate to the whole transcriptome (WT) assay in capturing chemical-induced changes in gene expression, as evident from correlation analysis, PCA and unsupervised hierarchical clustering. Gene set enrichment analysis (GSEA) using the Hallmark gene set collection also demonstrated consistency in enrichment scores between ST and WT assays. Lastly, unsupervised hierarchical clustering of hallmark enrichment scores broadly divided the samples into hepatotoxic, intermediate, and non-hepatotoxic groups. Xenobiotic metabolism, bile acid metabolism, apoptosis, p53 pathway, and coagulation were found to be the key hallmarks driving the clustering. Taken together, our results established the reproducibility and performance of collagen sandwich culture in combination with TempO-Seq S1500+ assay, and demonstrated the utility of GSEA using the hallmark gene set collection to identify potential hepatotoxicants for further validation.